Matlab读取mhd文件并转为dcm,做此记录,后续再进行更改。
clc; clear; close all
id = fopen('CT_FCN/CT_predicted_3.raw');
imgs = fread(id,'int16');%以'short'数据类型打开,因为本人存储的raw数据是以short保存的
imgsize = size(imgs);%读出的图像的size为n*1大小
rows = 320;%根据mhd文件的图像大小设定
clos =192;
nums = imgsize(1)/rows/clos; %计算得到图像张数
imgs = reshape(imgs,[rows,clos,nums]);
% img = int16(imgs(:,:,1)');
% imshow(imgs(:,:,1)',[])
for i = 1:26
img = int16(imgs(:,:,5*(i-1)+1)');
infopath = 'sun210816/sun210816_CT1gt_image00001.dcm';
imginfo = dicominfo(infopath);
imginfo.InstanceNumber=imginfo.InstanceNumber+i;%保存时保持层数与序列一致
imginfo.ImagePositionPatient=[imginfo.ImagePositionPatient(1), imginfo.ImagePositionPatient(2),imginfo.ImagePositionPatient(3)+5.0*(i-1)];
imginfo.SliceLocation=imginfo.SliceLocation+5.0*(i-1);
dicomwrite(img, ['fcn_dcm/', sprintf('%02d.dcm', 26-i)], imginfo);
end