0001 function [param_list, warped_points, warped_image, score] = optimise_groupwise(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, model_evaluation_method, n_modes, weighting_normalisation_method, shape_weight, variation_kept, pdf_evalution_method, kernel, 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
0049 [grid, steps] = setup_grid(2, unwarped_image, ref_points, placement_method, n_points);
0050
0051 if (do_plot),
0052 figure(subplot_fig);
0053 subplot(number_of_images, column, pane, 'replace');
0054 plot(unwarped_image);
0055 axis([0 image_width 0 1]);
0056 drawnow;
0057 end
0058
0059 warp_params.green = 'biharmCPS';
0060 param_list = [];
0061 start_image = unwarped_image;
0062 start_points = unwarped_points;
0063
0064 if(strcmp(spline_type,'single_point'))
0065 last_score = eval_groupwise_single_point(zeros(size(grid)),grid, steps, unwarped_image, unwarped_points, ref_image, ref_points, n_modes, start_points, model_evaluation_method, weighting_normalisation_method, shape_weight, variation_kept, pdf_evalution_method, kernel);
0066 [params, score] = fminsearch(@eval_groupwise_single_point,zeros(size(grid)),optimset('Display',verbose,'TolX', precision_required),grid,steps,start_image, start_points, ref_image, ref_points, n_modes, start_points, model_evaluation_method, weighting_normalisation_method, shape_weight, variation_kept, pdf_evalution_method, kernel);
0067 param_list(size(param_list,1)+1).d = params;
0068 [dummy,start_points] = eval_groupwise_single_point(params,grid,steps, start_image, start_points, ref_image, ref_points, n_modes, start_points, model_evaluation_method, weighting_normalisation_method, shape_weight, variation_kept, pdf_evalution_method, kernel);
0069 start_image = interp1(unwarped_points,unwarped_image,start_points, 'linear',0);
0070 elseif(strcmp(spline_type,'multi_point'))
0071 warped_grid = grid;
0072 warp_params.green = 'biharmCPS';
0073 last_score = eval_groupwise_multi_point(zeros(size(grid)),grid, unwarped_image, unwarped_points, ref_image, ref_points, n_modes, start_points, warp_params, model_evaluation_method, weighting_normalisation_method, shape_weight, variation_kept, pdf_evalution_method, kernel);
0074 [params, score] = fminsearch(@eval_groupwise_multi_point,zeros(size(grid)),optimset('Display',verbose,'TolX', precision_required),grid, unwarped_image, unwarped_points, ref_image, ref_points, n_modes, start_points, warp_params, model_evaluation_method, weighting_normalisation_method, shape_weight, variation_kept, pdf_evalution_method, kernel);
0075 param_list(size(param_list,1)+1,:).d = params;
0076 warped_grid = grid + params;
0077 start_points = nrr_trans_1d(start_points, grid, warped_grid, warp_params,[]);
0078 start_image = interp1(unwarped_points,unwarped_image,start_points, 'linear',0);
0079 else
0080 error('Unknown spline type');
0081 end
0082
0083
0084 if (verbose_score == 1),
0085 disp(['Initial score: ',num2str(last_score),' Final score: ', num2str(score), ' Difference: ', num2str(last_score - score)]);
0086 end
0087
0088 warped_points = start_points;
0089 warped_image = start_image;
0090
0091
0092 if (do_plot),
0093 figure(subplot_fig);
0094 subplot(number_of_images, column ,pane, 'replace');
0095 plot(warped_image);
0096 axis([0 image_width 0 1]);
0097 title(num2str(score));
0098 drawnow;
0099 end