CT dicom 5mm层厚插值成1mm-20201110

最新推荐文章于 2025-10-09 13:38:28 发布

原创

最新推荐文章于 2025-10-09 13:38:28 发布 · 2k 阅读

14 ·

CC 4.0 BY-SA版权

放射科的CT大多数是5 mm，有时需要插值成1 mm，这里是医工所的曹老师写的代码，在此做个记录，学习一下。

clear;clc;
%程序功能：读取dicom文件，并将其插值为1mm
folderPath = 'Y180521\';
outFolderPath = 'out\';

folderInfo = dir(folderPath);
fileNum = size(folderInfo,1)-2;
outFileNum = 0;
fileIndex = [];
instanceNumber = [];
for i = 1:fileNum
    %读取相邻两张图片
    fileName = folderInfo(i+2).name;
    filePath = strcat(folderPath, fileName);
%     ImgB = dicomread(filePath);  % 这里似乎不需要，只需要头文件信息
    ImgInfoB = dicominfo(filePath);%存储信息
%     disp(ImgInfoB.InstanceNumber);
    instanceNumber(end+1,1) = ImgInfoB.InstanceNumber;  % end+1后的,1可省略
end
startInstanceNumber = min(instanceNumber);  % 细心，注意设置起始number
for i = 0:(fileNum-2)
    %读取相邻两张图片
    fileIndexB = find(instanceNumber == (startInstanceNumber + i));    
    fileNameB = folderInfo(fileIndexB+2).name;
    filePathB = strcat(folderPath, fileNameB);
    ImgB = dicomread(filePathB);
    ImgInfoB = dicominfo(filePathB);%存储信息
    
    fileIndexT = find(instanceNumber == (startInstanceNumbe