From e44a0ee9e32819900d1aebaedd938655811b5143 Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Thu, 7 Jul 2011 15:58:24 -0400 Subject: make image matrix with beam traces --- make_beam_trace.m | 39 +++++++++++++++++++++++++++++++++++++++ 1 file changed, 39 insertions(+) create mode 100644 make_beam_trace.m diff --git a/make_beam_trace.m b/make_beam_trace.m new file mode 100644 index 0000000..c2257a6 --- /dev/null +++ b/make_beam_trace.m @@ -0,0 +1,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 + + -- cgit v1.2.3