eval_msd_multi_warp

 ==============================================================
 EVAL_MSD_MULTI_WARP: Applies a multi warp and evaluates the
                        warped image against the reference image.

 Code written by Katherine Smith, 2003

    GENERAL

      score = eval_msd_multi_warp(knot_displacements, init_knots,
              unwarped_image, unwarped_points, ref_image, ref_points,
              n_modes, warp_params)

    INPUT/S

      -knot_displacements:
            Some warp specific parameter.

      -init_knots:
            Some warp specific parameter.

      -unwarped_image:
            Images to be warped.

      -unwarped_points:
            Points to be warped (points define the images in a way).

      -ref_image:
            The reference image.

      -ref_points:
            The reference image points.

      -n_modes:
           The number of modes

      -warp_params:
            The parameters of the warp
           
    OUTPUT/S

      -score:
           Mean-squared difference score.

    PENDING WORK

      -

    KNOWN BUG/S

      -None.

    COMMENT/S

      -

    RELATED FUNCTION/S

      MSD, EVAL_MSD, EVAL_MSD_CPS_WARP

    ABOUT

      -Created:     November 23rd, 2003
      -Last update: Novermber 26th, 2003
      -Revision:    0.0.2
      -Author:      R. S. Schestowitz, University of Manchester
 ==============================================================

0001 function score = eval_msd_multi_warp(knot_displacements, init_knots, unwarped_image, unwarped_points, ref_image, ref_points, n_modes, warp_params)
0002 % ==============================================================
0003 % EVAL_MSD_MULTI_WARP: Applies a multi warp and evaluates the
0004 %                        warped image against the reference image.
0005 %
0006 % Code written by Katherine Smith, 2003
0007 %
0008 %    GENERAL
0009 %
0010 %      score = eval_msd_multi_warp(knot_displacements, init_knots,
0011 %              unwarped_image, unwarped_points, ref_image, ref_points,
0012 %              n_modes, warp_params)
0013 %
0014 %    INPUT/S
0015 %
0016 %      -knot_displacements:
0017 %            Some warp specific parameter.
0018 %
0019 %      -init_knots:
0020 %            Some warp specific parameter.
0021 %
0022 %      -unwarped_image:
0023 %            Images to be warped.
0024 %
0025 %      -unwarped_points:
0026 %            Points to be warped (points define the images in a way).
0027 %
0028 %      -ref_image:
0029 %            The reference image.
0030 %
0031 %      -ref_points:
0032 %            The reference image points.
0033 %
0034 %      -n_modes:
0035 %           The number of modes
0036 %
0037 %      -warp_params:
0038 %            The parameters of the warp
0039 %
0040 %    OUTPUT/S
0041 %
0042 %      -score:
0043 %           Mean-squared difference score.
0044 %
0045 %    PENDING WORK
0046 %
0047 %      -
0048 %
0049 %    KNOWN BUG/S
0050 %
0051 %      -None.
0052 %
0053 %    COMMENT/S
0054 %
0055 %      -
0056 %
0057 %    RELATED FUNCTION/S
0058 %
0059 %      MSD, EVAL_MSD, EVAL_MSD_CPS_WARP
0060 %
0061 %    ABOUT
0062 %
0063 %      -Created:     November 23rd, 2003
0064 %      -Last update: Novermber 26th, 2003
0065 %      -Revision:    0.0.2
0066 %      -Author:      R. S. Schestowitz, University of Manchester
0067 % ==============================================================
0068 
0069 % Smith:scale it all down to try and get it to move more
0070 % RSS: (when warping?)... might need some alteration
0071 
0072 % init_knots = init_knots*0.01;
0073 % unwarped_points = unwarped_points*0.01;
0074 
0075 warped_points = nrr_trans_1d(unwarped_points, init_knots, init_knots + knot_displacements, warp_params, []);
0076                 % warp points using <multi warp>
0077 warped_image = interp1(unwarped_points,unwarped_image,warped_points, 'linear',0);
0078                 % interpolate image to get cleaner warped image
0079 score = msd(ref_image, warped_image);
0080                 % evaluate warped image according to reference image