eval_groupwise_single_point
PURPOSE
OPTIMISE_GROUPWISE_SINGLE_POINT: Uses PDF's to minimise group-wise complexity
SYNOPSIS
function [score, warped_points] = eval_groupwise_single_point (params, warp_centre, r, unwarped_image_vec, unwarped_points_vec, image_set, points_set, n_modes, start_points_vec, model_evaluation_method, weighting_normalisation_method, shape_weight, variation_kept, pdf_evaluation_method, kernel)
DESCRIPTION
CROSS-REFERENCE INFORMATION
This function calls:
This function is called by:
SOURCE CODE
0001 function [score, warped_points] = eval_groupwise_single_point (params, warp_centre, r, unwarped_image_vec, unwarped_points_vec, image_set, points_set, n_modes, start_points_vec, model_evaluation_method, weighting_normalisation_method, shape_weight, variation_kept, pdf_evaluation_method, kernel)
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 start_points = unwarped_points_vec;
0048 max_d = 0.67;
0049
0050 if (any((abs(params) > max_d)) | any(warp_centre-r <-1) | any(warp_centre+r > 1))
0051 score = realmax;
0052 warped_points = unwarped_points_vec;
0053 else
0054 for i=1:size(params,2)
0055 warped_points = cps_warp_1d(start_points_vec,warp_centre(i),r(i),params(i));
0056 start_points = warped_points;
0057 end
0058 warped_image = interp1(unwarped_points_vec,unwarped_image_vec,warped_points, 'linear',0);
0059 image_set(:,size(image_set,2)+1) = warped_image;
0060 points_set(:,size(points_set,2)+1) = warped_points;
0061
0062
0063
0064
0065 groupwise_pdf = compute_pdf(image_set, points_set, kernel, shape_weight);
0066 score = measure_pdf(groupwise_pdf, pdf_evaluation_method, warped_image);
0067
0068
0069 end
Generated on Fri 14-May-2004 10:05:30 by m2html © 2003