Check the valid wavefront array index
%% check valid sub-aperture index
% rchen@niaot.ac.cn
clear all;
clc;
number_of_subaperture_in_a_row = 7;
k=1;
m=1;
for i = 0: 6
for j = 0:6
temp = (1- 1/number_of_subaperture_in_a_row - ...
(2/number_of_subaperture_in_a_row * i))^2 + ...
((2/number_of_subaperture_in_a_row * j) + ...
(1/number_of_subaperture_in_a_row -1))^2;
if (temp > 1)
invalid(k) = j +7*i;
k = k+1;
else
valid(m) = j +7*i;
m = m+1;
end
end
end