mhd格式三维图像显示_[转载]Matlab读取mha格式医学图像

mha格式图像，读取与显示，matlab语言

function info =mha_read_header(filename)

% Function for reading the header of a Insight Meta-Image

(.mha,.mhd) file

%

% info  = mha_read_header(filename);

%

% examples:

% 1,  info=mha_read_header()

% 2,  info=mha_read_header('volume.mha');

if(exist('filename','var')==0)

[filename, pathname] =

uigetfile('*.mha', 'Read mha-file');

filename = [pathname

filename];

end

fid=fopen(filename,'rb');

if(fid<0)

fprintf('could not open file

%sn',filename);

return

end

info.Filename=filename;

info.Format='MHA';

info.CompressedData='false';

readelementdatafile=false;

while(~readelementdatafile)

str=fgetl(fid);

s=find(str=='=',1,'first');

if(~isempty(s))

type=str(1:s-1);

data=str(s+1:end);

while(type(end)==' '); type=type(1:end-1);

end

while(data(1)==' '); data=data(2:end); end

else

type=''; data=str;

end

switch(lower(type))

case 'ndims'

info.NumberOfDimensions=sscanf(data, '%d')';

case 'dimsize'

info.Dimensions=sscanf(data, '%d')';

case 'elementspacing'

info.PixelDimensions=sscanf(data, '%lf')';

case 'elementsize'

info.ElementSize=sscanf(data, '%lf')';

if(~isfield(info,'PixelDimensions'))

info.PixelDimensions=info.ElementSize;

end

case 'elementbyteordermsb'

info.ByteOrder=lower(data);

case 'anatomicalorientation'

info.AnatomicalOrientation=data;

case 'centerofrotation'

info.CenterOfRotation=sscanf(data, '%lf')';

case 'offset'

info.Offset=sscanf(data, '%lf')';

case 'binarydata'

info.BinaryData=lower(data);

case 'compresseddatasize'

info.CompressedDataSize=sscanf(data, '%d')';

case 'objecttype',

info.ObjectType=lower(data);

case 'transformmatrix'

info.TransformMatrix=sscanf(data, '%lf')';

case 'compresseddata';

info.CompressedData=lower(data);

case 'binarydatabyteordermsb'

info.ByteOrder=lower(data);

case 'elementdatafile'

info.DataFile=data;

readelementdatafile=true;

case 'elementtype'

info.DataType=lower(data(5:end));

case 'headersize'

val=sscanf(data, '%d')';

if(val(1)>0), info.HeaderSize=val(1); end

otherwise

info.(type)=data;

end

end

switch(info.DataType)

case 'char',

info.BitDepth=8;

case 'uchar',

info.BitDepth=8;

case 'short',

info.BitDepth=16;

case 'ushort',

info.BitDepth=16;

case 'int',

info.BitDepth=32;

case 'uint',

info.BitDepth=32;

case 'float',

info.BitDepth=32;

case 'double',

info.BitDepth=64;

otherwise,

info.BitDepth=0;

end

if(~isfield(info,'HeaderSize'))

info.HeaderSize=ftell(fid);

end

fclose(fid);

获取mha格式图像文件头，得到像素类型，图像尺度。

function V = mha_read_volume(info)

% Function for reading the volume of a Insight Meta-Image (.mha,

.mhd) file

%

% volume = tk_read_volume(file-header)

%

% examples:

% 1: info = mha_read_header()

%  V =

mha_read_volume(info);

%  imshow(squeeze(V(:,:,round(end/2))),[]);

%

% 2: V = mha_read_volume('test.mha');

if(~isstruct(info)), info=mha_read_header(info); end

switch(lower(info.DataFile))

case 'local'

otherwise

% Seperate file

info.Filename=fullfile(fileparts(info.Filename),info.DataFile);

end

% Open file

switch(info.ByteOrder(1))

case ('true')

fid=fopen(info.Filename','rb','ieee-be');

otherwise

fid=fopen(info.Filename','rb','ieee-le');

end

switch(lower(info.DataFile))

case 'local'

% Skip header

fseek(fid,info.HeaderSize,'bof');

otherwise

fseek(fid,0,'bof');

end

datasize=prod(info.Dimensions)*info.BitDepth/8;

switch(info.CompressedData(1))

case 'f'

% Read the Data

switch(info.DataType)

case

'char'

V =

int8(fread(fid,datasize,'char'));

case

'uchar'

V =

uint8(fread(fid,datasize,'uchar'));

case

'short'

V =

int16(fread(fid,datasize,'short'));

case

'ushort'

V =

uint16(fread(fid,datasize,'ushort'));

case

'int'

V =

int32(fread(fid,datasize,'int'));

case

'uint'

V =

uint32(fread(fid,datasize,'uint'));

case

'float'

V =

single(fread(fid,datasize,'float'));

case

'double'

V =

double(fread(fid,datasize,'double'));

end

case 't'

switch(info.DataType)

case

'char', DataType='int8';

case

'uchar', DataType='uint8';

case

'short', DataType='int16';

case

'ushort', DataType='uint16';

case

'int', DataType='int32';

case

'uint', DataType='uint32';

case

'float', DataType='single';

case

'double', DataType='double';

end

Z  =

fread(fid,inf,'uchar=>uint8');

V = zlib_decompress(Z,DataType);

end

fclose(fid);

V = reshape(V,info.Dimensions);

function M = zlib_decompress(Z,DataType)

import com.mathworks.mlwidgets.io.InterruptibleStreamCopier

a=java.io.ByteArrayInputStream(Z);

b=java.util.zip.InflaterInputStream(a);

isc = InterruptibleStreamCopier.getInterruptibleStreamCopier;

c = java.io.ByteArrayOutputStream;

isc.copyStream(b,c);

M=typecast(c.toByteArray,DataType);