eval_wavelet_multi_point

 OPTIMISE_WAVELET_MULTI_POINT: 

    GENERAL


 
    INPUT/S

      -X:
          X
           
    OUTPUT/S

      -X:
           X

    PENDING WORK
      
      

    KNOWN BUG/S

      

    COMMENT/S

      

    RELATED FUNCTION/S

      

    ABOUT

      -Created:     February 9th, 2004
      -Last update: February 2004
      -Revision:    0.1.0
      -Author:      R. S. Schestowitz, University of Manchester
 ==============================================================

0001 function [score, warped_points] = eval_wavelet_multi_point(knot_displacements, init_knots, unwarped_image_vec, unwarped_points_vec, image_set, points_set, n_modes, start_points_vec, warp_params, model_evaluation_method, weighting_normalisation_method, shape_weight, variation_kept, wavelet_evaluation_method, filter, filter_level)
0002 % OPTIMISE_WAVELET_MULTI_POINT:
0003 %
0004 %    GENERAL
0005 %
0006 %
0007 %
0008 %    INPUT/S
0009 %
0010 %      -X:
0011 %          X
0012 %
0013 %    OUTPUT/S
0014 %
0015 %      -X:
0016 %           X
0017 %
0018 %    PENDING WORK
0019 %
0020 %
0021 %
0022 %    KNOWN BUG/S
0023 %
0024 %
0025 %
0026 %    COMMENT/S
0027 %
0028 %
0029 %
0030 %    RELATED FUNCTION/S
0031 %
0032 %
0033 %
0034 %    ABOUT
0035 %
0036 %      -Created:     February 9th, 2004
0037 %      -Last update: February 2004
0038 %      -Revision:    0.1.0
0039 %      -Author:      R. S. Schestowitz, University of Manchester
0040 % ==============================================================
0041 
0042 % scale it all down to try and get it to move more!
0043 % init_knots = init_knots*0.01;
0044 % start_points_vec = start_points_vec*0.01;
0045 
0046 warped_points = nrr_trans_1d(start_points_vec, init_knots, init_knots + knot_displacements, warp_params, []);
0047 warped_image = interp1(unwarped_points_vec,unwarped_image_vec,warped_points, 'linear',0);
0048 image_set(:,size(image_set,2)+1) = warped_image;
0049 points_set(:,size(points_set,2)+1) = warped_points;
0050 
0051 wavelet_evaluation_method = 'Default';
0052 
0053 compressed_data = compute_wavelet(image_set, points_set, filter, filter_level);
0054 score = measure_wavelet(compressed_data, wavelet_evaluation_method, warped_image);
0055 
0056 % score = 0; % to be modified w.r.t. eval_groupwise_single_point