showCurrentShotForNewFigure.m

% this script is specifically made for showing the current shot loaded by
% ImgTestGui into a new figure. This is called in ImgTestGui_ShowImage.m

% This script cannot run on its own. Based on image_processing_common_matter

camera = getappdata(0, 'camera');
dataset_str = getappdata(0, 'dataset_str');
bend_struc = getappdata(0, 'bend_struc');
i = getappdata(0, 'i');
j = getappdata(0, 'j');

if camera.energy_camera 
    if (isfield(bend_struc, 'variable_bend') && bend_struc.variable_bend)
        camera.dipole_bend = bend_struc.dipole_multiplier_values(i);
        camera.zero_Gev_px = bend_struc.zero_Gev_px_vector(i);                        
    end
    if camera.dipole_bend ~= 1
        stack_text{i} = ['Dipole=' num2str(camera.dipole_bend,2) ', ' ...
            stack_text{i}];
        setappdata(0, 'stack_text', stack_text);
    end
end
% sets current image to be used by load_noiseless_images_15_edited

% must make figure and axis current and selected when loading image        
f = figure('Name', ['Shot ' num2str(j) ' of ' 'Stack ' num2str(i) ' ' ...
    dataset_str], 'Position',[900 450 500 500]);        
figureImgTestGui = findobj('Tag', 'figureImgTestGui');
pushPreviousShot = findobj('Tag', 'pushPreviousShot');
axesImage = findobj('Tag', 'axesImage');
set(groot, 'CurrentFigure', f);

% load necessary variables not previously loaded
curr_lim = getappdata(0, 'curr_lim');
sebas_colors = sebasColorTable();
linearize = getappdata(0, 'linearize');
remove_xrays = getappdata(0, 'remove_xrays');

% this checkpoint is for right after setting the ROI for linearizing 
if linearize
    check_linearize_ROI = getappdata(0, 'check_linearize_ROI');
    load_linearize_ROI = getappdata(0, 'load_linearize_ROI');
    if check_linearize_ROI
        draw_ROI = getappdata(0, 'draw_ROI'); % CMOSFAR 2014 15 up to 
                                              % 35 GeV 13128 close up
        setappdata(0, 'load_linearize_ROI', 1);
        load_linearize_ROI = getappdata(0, 'load_linearize_ROI');
    end
end

%%%%%%%%%%%%%%%%%
%laser select code
% if (laser_select)
% laser_list = cell(3,1);
% laser_list{1,1} = 'laser on';
% laser_list{2,1} = 'laser off';
% laser_list{3,1} = 'all shots';
% setappdata(0, 'laser_list', laser_list);
% [s, ok] = listdlg('PromptString','Please select image type', 'SelectionMode', 'single', 'ListString',laser_list);
% setappdata(0, 's' , s);
% end
% s  = getappdata(0, 's');
% getappdata(0, 'laser_list');
% switch laser_list{s}
%     case 'laser on'
%         setappdata(0, 'j', 1);
%     case 'laser off'
%         setappdata(0, 'j', 2);
%     case 'all shots'
%         
% end
%         j = getappdata(0, 'j');
num_images = getappdata(0, 'num_images');
curr_image = (i - 1) * num_images + j; 
setappdata(0, 'curr_image', curr_image);
curr_image = getappdata(0, 'curr_image');

load_noiseless_images_15_edited;

this_image = getappdata(0, 'this_image'); % from load_noiseless_images_15_edited

%       %%%%%%%%%%%%%%%%%
        %show the image 
        % if show_image && mod(j,2)==1 % fix For odd(1) or even(0)
        if show_image % fix For odd(1) or even(0)
            image_axis = gca;
           
            if curr_image == 1
%                 colorbar
                if exist('draw_ROI', 'var')
%                     axis(gca, draw_ROI);
%                     display_width = draw_ROI(3) - draw_ROI(1);
%                     display_height = draw_ROI(4) - draw_ROI(2);
                    display_width = draw_ROI(3);
                    display_height = draw_ROI(4);
                    setappdata(0, 'display_width', display_width);
                    setappdata(0, 'display_height', display_height);
%                     width_start = draw_ROI(1);
%                     height_start = draw_ROI(3);
                    width_start = 1;
                    height_start = 1;
                    setappdata(0, 'height_start', height_start);
                    setappdata(0, 'width_start', width_start);
    %                 lineout(i, j, :) = sum(image(:, box_region), 2);
                    if exist('axis_limits_lineout', 'var')
                        axis_limits_lineout(3) = draw_ROI(3);
                        axis_limits_lineout(4) = draw_ROI(4);
                    end
                else
                    horizsize = getappdata(0, 'horizsize');
                    vertsize = getappdata(0, 'vertsize');
                    display_width = horizsize;
                    display_height = vertsize;
                    setappdata(0, 'display_width', display_width);
                    setappdata(0, 'display_height', display_height);
                    width_start = 1;
                    height_start = 1;
                    setappdata(0, 'height_start', height_start);
                    setappdata(0, 'width_start', width_start);
                end
                %TODO:This will need to be set by me eventually!
                contained_p = get(f, 'position');
                sub_p = get(gca, 'position');
                sub_width = sub_p(3) * contained_p(3);
                sub_height = sub_p(4) * contained_p(4);
                r1 = sub_width / sub_height;
                r2 = display_width / display_height;
                
                setappdata(0, 'contained_p', contained_p);
                setappdata(0, 'r1', r1);
                setappdata(0, 'r2', r2);
                
                set(gca, 'fontsize', 15); 
                % if the resizing pushes the window out of the border, this
                % command brings it back
                movegui(f);
            end        
            
            % make the original image if diagnostic is on, otherwise work
            % only with linearized (or otherwise processed?) image
            if log_color
                imagesc(log10(abs(this_image)))
                hC = colorbar;
                set(hC, 'Ytick', log_range);
                set(hC, 'YTicklabel', 10 .^ log_range);
            else
                imagesc(this_image)
                colorbar
            end
            if camera.energy_camera
                energy_pixel = getappdata(0, 'energy_pixel');
                energy_ticks = getappdata(0, 'energy_ticks');
                set(gca, 'YTick', energy_pixel);
                set(gca, 'YTickLabel', energy_ticks);
            end
            colormap(sebas_colors)
            ylabel(camera.y_label_text, 'fontsize', 15); 
            xlabel('Pixels', 'fontsize', 15);
            
            if log_color
                set(gca, 'CLim', [log_range(1) log_range(end)]);
            elseif ~strcmp(curr_lim, 'auto')
                set(gca, 'CLim', curr_lim);
            end
            if exist('axis_limits_image', 'var')
                axis(draw_ROI);
            end
%             grid

            % writing and drawing on the image
            height_start = getappdata(0, 'height_start');
            display_height = getappdata(0, 'display_height');
            width_start = getappdata(0, 'width_start');
            display_width = getappdata(0, 'display_width');
            % draws shot number textbox and stack specific setting if
            % multiple stacks            
            text(.02 * display_width + width_start, .03 * display_height + ...
                height_start, ['Shot # ' int2str(j)], 'color', 'w', ...
                'backgroundcolor', 'b');
            if exist('stack_text', 'var')
                text(.02 * display_width + width_start, .1 * display_height ...
                    + height_start, stack_text{i}, 'color', 'w', ...
                    'backgroundcolor', 'blue');
            end
        end
     
% linearizes by allowing user to select region of interest
if (linearize)
    if (camera.energy_camera && ~check_linearize_ROI)    
        msgbox('Please select region of interest to adjust energy axis');
        set(groot, 'CurrentFigure', f);
        h = imrect;
        fcn = makeConstrainToRectFcn('imrect', get(gca, 'XLim'), ...
            get(gca, 'YLim'));
        setPositionConstraintFcn(h, fcn);
        draw_ROI = wait(h);
        setappdata(0, 'draw_ROI', draw_ROI);
        setappdata(0, 'check_linearize_ROI', 1);
        pushLoadImages_Callback(hObject, eventdata, handles);   
    end    
end