blob: c2257a694b890a8825c9060483ca66c5c9f1799b (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
|
function img=make_beam_trace(beam, stop_point, border_limits, img)
% img so far beam traced part
% border_limits has coordinates of left bottom and right top coners
xlb=border_limits(1);
ylb=border_limits(2);
xrt=border_limits(3);
yrt=border_limits(4);
% beam start stop coordinates with respect to lower left border point
xb1=beam.origin(1)-xlb;
yb1=beam.origin(2)-ylb;
xb2=stop_point(1)-xlb;
yb2=stop_point(2)-ylb;
% beam path coordinates
Nc=100; % number of coordinate points
xb=linspace(xb1,xb2, Nc);
yb=linspace(yb1,yb2, Nc);
% beam coordinate to image pixel coordinate
[Ny,Nx]=size(img);
px=round(xb/(xrt-xlb)*(Nx-1) + 1);
py=round(yb/(yrt-ylb)*(Ny-1) + 1);
% truncate to borders
py=py( (1<=px) & (px<=Nx) );
px=px( (1<=px) & (px<=Nx) );
px=px( (1<=py) & (py<=Ny) );
py=py( (1<=py) & (py<=Ny) );
for i=1:length(px)
img(py(i),px(i))=beam.intensity;
end
end
|
