Changeset 77 for an2_roi_stats.m


Ignore:
Timestamp:
Feb 5, 2009, 10:16:41 AM (10 years ago)
Author:
tjniskan
Message:
  • Fixed a strange bug related to resizing Aedes window in Matlab

R2008a (and possibly R2008b and R2009a).

  • Fixed bugs of Voxel TimeSeries? and Overlay Control figure

handles not being properly destroyed on exit.

  • Changed default values for reading VNMR files (now Fast Read is on

by default, output precision is single and data is oriented using
PROCPAR.orient)

  • Fixed a bug in the calculate SNR -plugin
  • Made a temporary workaround for controlling the output in the "View

ROI statistics" -table. A proper GUI implementation will (hopefully) follow
shortly. The output can be controlled by using Aedes preferences. For
example setpref('Aedes','StatPrintDirs?','TXYZV') will print in the
output Total, X-dir, Y-dir, Z-dir and V-dir statistics.

M an2_revision.m
M an2_data_read.m
M an2_resviewer.m
M aedes.m
M an2_res2table.m
M plugins/calculate_snr.m
M an2_readfid.m
M an2_readfidprefs.m
M an2_roi_stats.m

File:
1 edited

Legend:

Unmodified
Added
Removed
  • an2_roi_stats.m

    r74 r77  
    115115      Res.Stat(ii).Max = max(data);
    116116      Res.Stat(ii).Min = min(data);
    117       Res.Stat(ii).PixelCount = length(data);     
    118 if isempty(data)
    119 Res.Stat(ii).Mean = NaN;
    120 Res.Stat(ii).Std = NaN;
    121 Res.Stat(ii).Sum = NaN;
    122 Res.Stat(ii).Max = NaN;
    123 Res.Stat(ii).Min = NaN;
    124 Res.Stat(ii).PixelCount = 0;
    125 end
    126 
     117      Res.Stat(ii).PixelCount = length(data);
     118      if isempty(data)
     119        Res.Stat(ii).Mean = NaN;
     120        Res.Stat(ii).Std = NaN;
     121        Res.Stat(ii).Sum = NaN;
     122        Res.Stat(ii).Max = NaN;
     123        Res.Stat(ii).Min = NaN;
     124        Res.Stat(ii).PixelCount = 0;
     125      end
     126     
    127127      %% Calculate results in X direction
    128128      Res.Stat(ii).XD.Mean = [];
     
    214214        end
    215215      end
     216     
     217      %% Calculate results in V direction
     218      Res.Stat(ii).VD.Mean = [];
     219      Res.Stat(ii).VD.Std = [];
     220      Res.Stat(ii).VD.Sum = [];
     221      Res.Stat(ii).VD.Max = [];
     222      Res.Stat(ii).VD.Min = [];
     223      Res.Stat(ii).VD.PixelCount = [];
     224     
     225      for jj=1:size(ROI(ii).voxels{kk},4)
     226        roiv=ROI(ii).voxels{kk}(:,:,:,jj);
     227        datav = DATA{kk}.FTDATA(:,:,:,jj);
     228        tmpv=datav(roiv);
     229        tmpv=double(tmpv);
     230        if isempty(tmpv)
     231          Res.Stat(ii).VD.Mean(end+1) = NaN;
     232          Res.Stat(ii).VD.Std(end+1) = NaN;
     233          Res.Stat(ii).VD.Sum(end+1) = NaN;
     234          Res.Stat(ii).VD.Max(end+1) = NaN;
     235          Res.Stat(ii).VD.Min(end+1) = NaN;
     236          Res.Stat(ii).VD.PixelCount(end+1) = 0;
     237        else
     238          Res.Stat(ii).VD.Mean(end+1) = mean(tmpv);
     239          Res.Stat(ii).VD.Std(end+1) = std(tmpv);
     240          Res.Stat(ii).VD.Sum(end+1) = sum(tmpv);
     241          Res.Stat(ii).VD.Max(end+1) = max(tmpv);
     242          Res.Stat(ii).VD.Min(end+1) = min(tmpv);
     243          Res.Stat(ii).VD.PixelCount(end+1) = length(tmpv);
     244        end
     245      end
    216246    end
    217247  end
Note: See TracChangeset for help on using the changeset viewer.

Powered by Trac 1.0.9.Copyright © Juha-Pekka Niskanen 2008