0001 function [param_list, warped_points, warped_image, score] = optimise_all_warps_mi(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, number_of_bins, precision_required)
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018
0019
0020
0021
0022
0023
0024
0025
0026
0027
0028
0029
0030
0031
0032
0033
0034
0035
0036
0037
0038
0039
0040
0041
0042
0043
0044
0045
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 = mi(ref_image, unwarped_image, number_of_bins);
0063 if(strcmp(spline_type,'single_point'))
0064 [params, score] = fminsearch(@eval_mi_cps_warp,zeros(size(grid)),optimset('Display',verbose,'TolX', precision_required),grid,steps,start_points, start_image, ref_image, number_of_bins);
0065 param_list(size(param_list,1)+1).d = params;
0066 [dummy,start_points] = eval_mi_cps_warp(params,grid,steps, start_points, start_image, ref_image, number_of_bins);
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_mi_multi_warp(zeros(size(grid)),grid, unwarped_image, unwarped_points, ref_image, ref_points, 20, warp_params, number_of_bins);
0072 [params, score] = fminsearch(@eval_mi_multi_warp,zeros(size(grid)),optimset('Display',verbose,'TolX', precision_required),grid, unwarped_image, unwarped_points, ref_image, ref_points, 20, warp_params, number_of_bins);
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 if (do_plot),
0089 figure(subplot_fig);
0090 subplot(number_of_images, column, pane, 'replace');
0091 plot(warped_image);
0092 axis([0 image_width 0 1]);
0093 title(num2str(score));
0094 drawnow;
0095 end