Home > rest_20090422 > rest_spm5_files > nic_spm_get_data.m

nic_spm_get_data

PURPOSE ^

gets data from image files at specified locations

SYNOPSIS ^

function [Y] = nic_spm_get_data(V,XYZ)

DESCRIPTION ^

 gets data from image files at specified locations
 FORMAT [Y] = spm_get_data(V,XYZ);

 V    - [1 x n] struct array of file handles (or filename matrix)
 XYZ  - [4 x m] or [3 x m]location matrix (voxel)

 Y    - (n x m) double values

 see spm_sample_vol
__________________________________________________________________________
 Copyright (C) 2005 Wellcome Department of Imaging Neuroscience

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

DOWNLOAD ^

nic_spm_get_data.m

SOURCE CODE ^

0001 function [Y] = nic_spm_get_data(V,XYZ)
0002 % gets data from image files at specified locations
0003 % FORMAT [Y] = spm_get_data(V,XYZ);
0004 %
0005 % V    - [1 x n] struct array of file handles (or filename matrix)
0006 % XYZ  - [4 x m] or [3 x m]location matrix (voxel)
0007 %
0008 % Y    - (n x m) double values
0009 %
0010 % see spm_sample_vol
0011 %__________________________________________________________________________
0012 % Copyright (C) 2005 Wellcome Department of Imaging Neuroscience
0013 
0014 % Karl Friston
0015 % $Id: spm_get_data.m 905 2007-09-04 14:38:28Z guillaume $
0016 
0017 
0018 % ensure V is an array of handle structures
0019 %--------------------------------------------------------------------------
0020 if ~isstruct(V)
0021     V = spm_vol(V);
0022     try
0023         V = cat(2,V{:});
0024     end
0025 end
0026 
0027 % get data
0028 %--------------------------------------------------------------------------
0029 Y     = zeros(length(V),size(XYZ,2));
0030 for i = 1:length(V)
0031     
0032         % check files exists, if not try pwd
0033     %----------------------------------------------------------------------
0034     if exist(V(i).fname,'file') ~=2
0035         [p,n,e]    = fileparts(V(i).fname);
0036         V(i).fname = [n e];
0037     end
0038 
0039     %-Load mask image within current mask & update mask
0040     %----------------------------------------------------------------------
0041     Y(i,:) = spm_sample_vol(V(i),XYZ(1,:),XYZ(2,:),XYZ(3,:),0);
0042 end
0043

Generated on Wed 29-Apr-2009 01:06:38 by m2html © 2005