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