optimise_all_warps_msd

 OPTIMISE_ALL_WARPS_MSD: 

 Code written by Katherine Smith, 2003

    GENERAL

      [param_list, warped_points, warped_image, score] = 
          optimise_all_warps_msd(ref_image, ref_points, unwarped_image,
          unwarped_points, spline_type, placement_method, n_points, verbose,
          do_plot, subplot_fig, pane,
          image_width, number_of_images, verbose_score)         

    INPUT/S

      -X:
          X
           
    OUTPUT/S

      -X:
           X

    PENDING WORK
      
      

    KNOWN BUG/S

      

    COMMENT/S

      

    RELATED FUNCTION/S

      

    ABOUT

      -Created:     November 23rd, 2003
      -Last update: March 8th, 2004
      -Revision:    0.1.6
      -Author:      R. S. Schestowitz, University of Manchester
 ==============================================================

0001 function [param_list, warped_points, warped_image, score] = optimise_all_warps_msd(ref_image, ref_points, unwarped_image, unwarped_points, spline_type, placement_method, n_points, verbose, do_plot, subplot_fig, pane, column, image_width, number_of_images, verbose_score, precision_required)
0002 % OPTIMISE_ALL_WARPS_MSD:
0003 %
0004 % Code written by Katherine Smith, 2003
0005 %
0006 %    GENERAL
0007 %
0008 %      [param_list, warped_points, warped_image, score] =
0009 %          optimise_all_warps_msd(ref_image, ref_points, unwarped_image,
0010 %          unwarped_points, spline_type, placement_method, n_points, verbose,
0011 %          do_plot, subplot_fig, pane,
0012 %          image_width, number_of_images, verbose_score)
0013 %
0014 %    INPUT/S
0015 %
0016 %      -X:
0017 %          X
0018 %
0019 %    OUTPUT/S
0020 %
0021 %      -X:
0022 %           X
0023 %
0024 %    PENDING WORK
0025 %
0026 %
0027 %
0028 %    KNOWN BUG/S
0029 %
0030 %
0031 %
0032 %    COMMENT/S
0033 %
0034 %
0035 %
0036 %    RELATED FUNCTION/S
0037 %
0038 %
0039 %
0040 %    ABOUT
0041 %
0042 %      -Created:     November 23rd, 2003
0043 %      -Last update: March 8th, 2004
0044 %      -Revision:    0.1.6
0045 %      -Author:      R. S. Schestowitz, University of Manchester
0046 % ==============================================================
0047 
0048 [grid, steps] = setup_grid(2, unwarped_image, ref_points, placement_method, n_points);
0049 
0050 if (do_plot),
0051   figure(subplot_fig);
0052   subplot(number_of_images, column, pane, 'replace');
0053   plot(unwarped_image);
0054   axis([0 image_width 0 1]);
0055   drawnow;
0056 end
0057 
0058 warp_params.green = 'biharmCPS';
0059 param_list = [];
0060 start_image = unwarped_image;
0061 start_points = unwarped_points;
0062 last_score = msd(ref_image, unwarped_image);
0063 if(strcmp(spline_type,'single_point')),
0064   [params, score] = fminsearch(@eval_msd_cps_warp,zeros(size(grid)),optimset('Display',verbose,'TolX', precision_required),grid,steps,start_points, start_image, ref_image);
0065   param_list(size(param_list,1)+1).d = params;
0066   [dummy,start_points] = eval_msd_cps_warp(params,grid,steps, start_points, start_image, ref_image);
0067   start_image = interp1(ref_points,unwarped_image,start_points, 'linear',0);  
0068 elseif(strcmp(spline_type,'multi_point')),
0069   warped_grid = grid;
0070   warp_params.green = 'biharmCPS';
0071   last_score = eval_msd_multi_warp(zeros(size(grid)),grid, unwarped_image, unwarped_points, ref_image, ref_points, 20, warp_params);
0072   [params, score] = fminsearch(@eval_msd_multi_warp,zeros(size(grid)),optimset('Display',verbose,'TolX', precision_required),grid, unwarped_image, unwarped_points, ref_image, ref_points, 20, warp_params);
0073   param_list(size(param_list,1)+1,:).d = params;
0074   warped_grid = grid + params;
0075   start_points = nrr_trans_1d(start_points, grid, warped_grid, warp_params,[]);
0076   start_image = interp1(unwarped_points,unwarped_image,start_points, 'linear',0);
0077 else
0078   error('Unknown spline type');      
0079 end
0080 
0081 if (verbose_score == 1),
0082       disp(['Initial score: ',num2str(last_score),' Final score: ', num2str(score), ' Difference: ', num2str(last_score - score)]);
0083 end
0084   
0085 warped_points = start_points;
0086 warped_image = start_image;
0087 
0088 
0089 if (do_plot),
0090   figure(subplot_fig);
0091   subplot(number_of_images, column, pane, 'replace');
0092   plot(warped_image);
0093   axis([0 image_width 0 1]);
0094   title(num2str(score));
0095   drawnow;
0096 end