% ss_write_edge_images( filestem )
% -----------------------------------------------------------------------------------
%
% Component of Scale Space Edges, version 1.2b.
%
% Write image files of the original and various summaries of edge information.                         
%
% Mark Dow, last modified May 2003
% University of Oregon
% Brain Development Lab
% dow@braindev.uoregon.edu
%


function ss_write_edge_images( filestem )


if ~(nargin == 1 )
    
    msg = ['\n\nThe ss_write_edge_images function needs 2 arguements, the input and output  file stems.  You have given ' num2str(nargin) ' arguments.\n\n'];
    fprintf( msg );
    return
end

load( [ filestem '_edge_info' ] )

sz_Image = size( flim_gray_cropped );


% Write grayscale image.
imwrite( flim_gray_cropped, [ filestem '_gray.tif' ] );


% Write gradient of edges image.
flim_gradient = flim_gradient*10;
flim_gradient(flim_gradient>1) = 1;
imwrite( flim_gradient, [ filestem '_edges.tif' ] );

% Generate average gradient of segment image.
flim_ave_gradient = zeros( sz_Image );
for i = 1:seg_n    
    for n_ep = seg(i,1) : seg(i,2)
        
        if seg(i,5) > flim_ave_gradient( ep(n_ep, 1), ep(n_ep, 2) )
            flim_ave_gradient( ep(n_ep, 1), ep(n_ep, 2) ) = seg(i,5);
        end
    end
end
% Write average gradient of segment image.
flim_ave_gradient = flim_ave_gradient*10;
flim_ave_gradient( flim_ave_gradient > 1 ) = 1;
imwrite( flim_ave_gradient, [ filestem '_edges_ave.tif' ] );