I'm attempting to take values in space and if they fall within a predetermined circle to change the value of that point in a blank matrix to one with zeros elsewhere.
After the first x,y plane is checked for points falling within the circle, the for loop index's to the next 'z' point and reexamines the x,y plane at this new value of z to see what points of the xy plane fall within the circle (that has expanded due to it being a cross section of a cone).
When I attempt to run it, I get the error that B(1,15,1) cannot be accessed as the index needs to be real or logical yet 1,15 and 1 are all logical and the for loop has managed to index up to 15 in the y dimension. Any suggestions?
Here is my FOR loop;
B=zeros(50,50,length(Z));
for z=1:length(Z)
c=(z*tan(15))^2;
for x=1:0.1:5
for y=1:0.1:5
if c>(x-x_0)^2+(y-y_0)^2
i=(10*x)-9
j=(10*y)-9
B(i,j,z)=1;
else
i=(10*x)-9
j=(10*y)-9
B(i,j,z)=0;
end
end
end
I use 'i' and 'j' to remove the fact x and y are indexing in 0.1 increments, x_0,y_0 and z_0 are the location of the apex of the cone and falls at (25,25,0). It is probably not the most efficient way of doing this.
