我有3D CT图像的二维切片。它们是DICOM格式的,其中有250个。我想用MATLAB重建三维图像。我怎么才能在循环中做到这一点呢?
/home/amadeus/Desktop/images发布于 2012-10-09 13:55:51
显然,有一个函数只用于读取DICOM文件:双读,我建议使用它来加载图像,然后将它们存储在3D矩阵中。斯普林特可以用来构造图像的文件名(使用%04d生成带前导零的四位数字)。
假设所有图像都对齐,大小相同:
N = 250;
img_dir = '/home/amadeus/Desktop/images'
% read the first image separately just to get the size
strfile = 'IM-0001-0001.dcm';
img = dicomread(fullfile(img_dir, strfile));
siz_img = size(img);
% create result matrix:
ct3d = NaN([siz_img N]);
ct3d(:,:,1) = img;
% load all the remaining images and put them in the matrix
for ii=2:N
strfile = sprintf('IM-0001-%04d.dcm',ii);
ct3d(:,:,ii)= dicomread(fullfile(img_dir, strfile));
end编辑:这假设图像是灰度(2d)。如果它们是完全彩色的(宽度x高度x3),您应该在赋值给ct3d中添加另一个冒号运算符。
https://stackoverflow.com/questions/12800333
复制相似问题