reslice_tms.m

mri_toolbox是一个工具用来MRI. 来自于SourceForge, 我上传这个软件,希望能结识对医疗软件感兴趣的兄弟.

function reslice_tms(fname)
%
% function reslice_tms(fname)
%
% Takes an analyze file acquired sagitally on the GE scanner,
% converted to SPM format using GE2SPM (version 3.2 or later)
% and rewrites it axially in radiological coordinates
% 
% Use for the highres and for also for the blobs resliced to
% space of the highres.
% makes tms_fname.img and tms_fname.hdr which are standard analyze
% format in radiological convention.
%
% Souheil Inati
% Dartmouth College
% Jan. 2002
%

% Get the spm_vol structure
V = spm_vol(fname);

% Get the info in the header
[DIM VOX SCALE TYPE OFFSET ORIGIN DESCRIP] = spm_hread(fname);

% Make a new header with things swapped around
[dirname,imgname] = fileparts(V.fname);
ax_fname = fullfile(dirname,['tms_' imgname '.img']);
ax_DIM = [DIM(3) DIM(1) DIM(2)];
ax_VOX = [VOX(3) VOX(1) VOX(2)];
ax_SCALE = SCALE;
ax_TYPE = TYPE;
ax_OFFSET = OFFSET;
ax_DESCRIP = ['axialized for TMS ' DESCRIP];
ax_ORIGIN = [DIM(3)-ORIGIN(3) DIM(1)-ORIGIN(1) DIM(2)-ORIGIN(2)];
spm_hwrite(ax_fname,ax_DIM,ax_VOX,ax_SCALE,ax_TYPE,ax_OFFSET,ax_ORIGIN,ax_DESCRIP);

V2 = spm_vol(ax_fname);

% Read in the data
[Y,XYZ] = spm_read_vols(V);

% Reorient
data = permute(Y,[3 1 2]);
% Deal with GE image flip
data = flipdim(data,2);  % flip y
data = flipdim(data,3);  % flip z
% Deal with radiological flip
data = flipdim(data,1);  % flip x  (make it radiological)

% Write out the reoriented data
spm_write_vol(V2,data);

return