show_intensity_model

 ==============================================================
 SHOW_INTENSITY_MODEL: Shows the intensity model in the form of
                       a figure.

 Code written by Katherine Smith, 2003

    GENERAL

      show_intensity_model(intensity_model, n_sds, max_modes, fig_title)

    INPUT/S

      -intensity_model:
           The intensity model to be shown.

      -n_sds:
           Number of standard deviations.

      -max_modes: 
           Maximum modes, meaning what??

      -fig_title:
           Title of the figure to be displayed?

           
    OUTPUT/S

      -None.

    PENDING WORK

      -

    KNOWN BUG/S

      -None.

    COMMENT/S

      -

    RELATED FUNCTION/S

      SHOW_SHAPE_MODEL, SHOW_COMBINED_MODEL

    ABOUT

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

0001 function figure_handle = show_intensity_model(intensity_model, n_sds, max_modes, fig_title, image_width)
0002 % ==============================================================
0003 % SHOW_INTENSITY_MODEL: Shows the intensity model in the form of
0004 %                       a figure.
0005 %
0006 % Code written by Katherine Smith, 2003
0007 %
0008 %    GENERAL
0009 %
0010 %      show_intensity_model(intensity_model, n_sds, max_modes, fig_title)
0011 %
0012 %    INPUT/S
0013 %
0014 %      -intensity_model:
0015 %           The intensity model to be shown.
0016 %
0017 %      -n_sds:
0018 %           Number of standard deviations.
0019 %
0020 %      -max_modes:
0021 %           Maximum modes, meaning what??
0022 %
0023 %      -fig_title:
0024 %           Title of the figure to be displayed?
0025 %
0026 %
0027 %    OUTPUT/S
0028 %
0029 %      -None.
0030 %
0031 %    PENDING WORK
0032 %
0033 %      -
0034 %
0035 %    KNOWN BUG/S
0036 %
0037 %      -None.
0038 %
0039 %    COMMENT/S
0040 %
0041 %      -
0042 %
0043 %    RELATED FUNCTION/S
0044 %
0045 %      SHOW_SHAPE_MODEL, SHOW_COMBINED_MODEL
0046 %
0047 %    ABOUT
0048 %
0049 %      -Created:     November 23rd, 2003
0050 %      -Last update: Novermber 26th, 2003
0051 %      -Revision:    0.0.2
0052 %      -Author:      R. S. Schestowitz, University of Manchester
0053 % ==============================================================
0054 
0055 %n_modes = size(intensity_model.pcs,2);
0056 
0057 n_modes = min(size(intensity_model.pcs,2),max_modes); 
0058                % get number of modes. The one above (in comments)
0059                % is probably the remains of an older version.
0060 
0061 fig = figure('Name',([fig_title,'Intensity modes'])); 
0062                     % prepare figure and assign name input
0063 figure_handle = figure(fig);
0064 
0065 
0066 for m=1:n_modes % for all modes as resolved above
0067   for i=-n_sds:n_sds  % for range of negative deviations to possitive ones
0068                       % with (currently) one unit steps (too large??)
0069     % Smith: what are the parameters corresponding to this sd
0070     % create example with these parameters
0071     example_params = repmat(i*intensity_model.sd(1,m),size(intensity_model.pcs,1),1);
0072                        % get the parameters of the models
0073     example = example_params.*intensity_model.pcs(:,m) + intensity_model.mean';
0074                        % apply the changes according to the mode parameters. Use
0075                        % them as an offset to the mean intensity model
0076     figure(fig);
0077 
0078     subplot(n_modes, n_sds * 2 + 1, (m - 1) * (n_sds * 2 + 1) + i + n_sds + 1);
0079     plot(example);
0080                      % plot the different modes of intensity variation
0081     axis([0 image_width -1 1.4]);
0082     if (m == 1 & i ==1),
0083       text(1,1,'Intensity Model');
0084     end
0085   end
0086 end