diff options
Diffstat (limited to 'beam_trace.m')
| -rw-r--r-- | beam_trace.m | 10 |
1 files changed, 7 insertions, 3 deletions
diff --git a/beam_trace.m b/beam_trace.m index 0896975..90359d4 100644 --- a/beam_trace.m +++ b/beam_trace.m @@ -1,10 +1,11 @@ -function img = beam_trace(beams, faces, borders, border_limits, img) +function [processed_beams] = beam_trace(beams, faces, borders, border_limits) %% trace beam and all it reflections and refractions on faces % for beams and faces structure description see face_beam_interaction.m % border similar to faces cell array which enclose the beam % i.e. it does not leave the polygon consisting of border faces intensity_threshold = 1e-3; % intensity of the weakest beam which we still trace + processed_beams = {}; Nbeams=size(beams)(2); if ( Nbeams == 0 ) @@ -17,6 +18,7 @@ function img = beam_trace(beams, faces, borders, border_limits, img) beam=beams{i}; [is_face_hit, hit_position, hit_distance, new_beams] = face_beam_interaction(beam, faces); if (is_face_hit) + beams{i}.hit_position=hit_position; %% remove beams which are to weak to trace N_new_beams=size(new_beams)(2); intense_enough_beams={}; @@ -27,7 +29,8 @@ function img = beam_trace(beams, faces, borders, border_limits, img) intense_enough_beams{intense_beams_counter}=new_beams{k}; end end - img = beam_trace(intense_enough_beams, faces, borders, border_limits, img); + [new_processed_beams] = beam_trace(intense_enough_beams, faces, borders, border_limits); + processed_beams=cat(2,processed_beams, new_processed_beams); else % beam does not hit face but it should stop and borders beam.face=NA; @@ -35,9 +38,10 @@ function img = beam_trace(beams, faces, borders, border_limits, img) if (!is_face_hit) error('borders are badly defined, the beam misses them'); end + beams{i}.hit_position=hit_position; end - img=make_beam_trace(beam, hit_position, border_limits, img); end + processed_beams=cat(2,beams, processed_beams); end |