Checking if 3D straight line intersects a rectangular grid cell.
显示 更早的评论
I have a 3D logical matrix (each element means: "filled" or "empty") representing a 3D volume (spatial step in each dimension is the same). I need to check, which grid cells is crossing a beam (sraight line) going from a particular point in particular direction (vertical angle and azimuth). Are there any standard scripts for it?
Thanks
采纳的回答
Shadaab Siddiqie
2020-10-7
From my understanding you want to find all the cells in a 3D matrix which is intersected by an 3D line (given vertical angle and azimuth). You could convert spatial coordinates to cartesian coordinate to map it with the 3D matrix and then make small increments to find all the cells in the matrix that are in your path (with some conditions). Here is a simple script that might be helpfull.
A = zeros(3, 3, 3);
teta = pi/10; % vertical angle
phi = pi/6; % azimuth
smallstep = 1; % you can make it smaller for more precision eg:0.5
for r = 0:1:3
x = r*sin(teta)*cos(phi);
y = r*sin(teta)*sin(phi);
z = r*cos(teta);
if A(round(x),round(y),round(z)) == 0: % any condition that you need
disp([round(x) round(y) round(z)]);
end
end
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Cartesian Coordinate System Conversion 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
