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

放射科的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 == (startInstanceNumber + i + 1));fileNameT = folderInfo(fileIndexT+2).name;filePathT = strcat(folderPath, fileNameT);ImgT = dicomread(filePathT);ImgInfoT = dicominfo(filePathT);%存储信息%线性插值ImgD = double(ImgT - ImgB) * 0.2;Img1 = int16(round(ImgD * 1 + double(ImgB)));Img2 = int16(round(ImgD * 2 + double(ImgB)));Img3 = int16(round(ImgD * 3 + double(ImgB)));Img4 = int16(round(ImgD * 4 + double(ImgB)));%保存结果filePath = strcat( outFolderPath,num2str(i*5),'.DCM');ImgInfoB.InstanceNumber = i * 5 + 1;ImgInfoB.SliceThickness = 1;ImgInfoB.Filename = filePath;dicomwrite(ImgB,filePath,ImgInfoB);%写入Dicom图像filePath = strcat( outFolderPath,num2str(i*5+1),'.DCM');ImgInfoB.InstanceNumber = i * 5 + 2;ImgInfoB.SliceThickness = 1;ImgInfoB.Filename = filePath;ImgInfoB.SliceLocation = ImgInfoB.SliceLocation - 1;ImgInfoB.ImagePositionPatient(3,1) = ImgInfoB.ImagePositionPatient(3,1)-1;dicomwrite(Img1,filePath,ImgInfoB);%写入Dicom图像filePath = strcat( outFolderPath,num2str(i*5+2),'.DCM');ImgInfoB.InstanceNumber = i * 5 + 3;ImgInfoB.SliceThickness = 1;ImgInfoB.Filename = filePath;ImgInfoB.SliceLocation = ImgInfoB.SliceLocation - 1;ImgInfoB.ImagePositionPatient(3,1) = ImgInfoB.ImagePositionPatient(3,1)-1;dicomwrite(Img2,filePath,ImgInfoB);%写入Dicom图像filePath = strcat( outFolderPath,num2str(i*5+3),'.DCM');ImgInfoB.InstanceNumber = i * 5 + 4;ImgInfoB.SliceThickness = 1;ImgInfoB.Filename = filePath;ImgInfoB.SliceLocation = ImgInfoB.SliceLocation - 1;ImgInfoB.ImagePositionPatient(3,1) = ImgInfoB.ImagePositionPatient(3,1)-1;dicomwrite(Img3,filePath,ImgInfoB);%写入Dicom图像filePath = strcat( outFolderPath,num2str(i*5+4),'.DCM');ImgInfoB.InstanceNumber = i * 5 + 5;ImgInfoB.SliceThickness = 1;ImgInfoB.Filename = filePath;ImgInfoB.SliceLocation = ImgInfoB.SliceLocation - 1;ImgInfoB.ImagePositionPatient(3,1) = ImgInfoB.ImagePositionPatient(3,1)-1;dicomwrite(Img4,filePath,ImgInfoB);%写入Dicom图像filePath = strcat( outFolderPath,num2str((i+1)*5+1),'.DCM');ImgInfoT.InstanceNumber = (startInstanceNumber + i) * 5 + 1;ImgInfoT.SliceThickness = 1;ImgInfoT.Filename = filePath;dicomwrite(ImgT,filePath,ImgInfoT);%写入Dicom图像
end

CT dicom 5mm层厚插值成1mm-20201110相关推荐

层间距离对ct图像的影响_多层螺旋CT螺距、层厚和重建间隔关系及对图像质量影响...
多层螺旋 CT 螺距.层厚和重建间隔关系及对图像质量影响周泽俊 ; 胡永胜 ; 高斌 ; 巢惠民 [期刊名称] <中国医疗设备> [年 ( 卷 ), 期] 2008(023)008 [摘 ...
【转】CT层厚、层间距、层间隔的概念是什么，MRI的层厚、层间距、曾间隔是什么
CT层厚指扫描层的厚度.CT层间距指两个扫描层面中心之间的距离.CT层间隔指两层之间的距离. 做CT打比方就是把胡萝卜切成片,观察每片的结构. 那胡萝卜需要切多厚呢?这个葫芦卜片的厚度就是这个扫描层厚 ...
【转】CT中的“层“与“排“的区别
转自:CT中的"层"与"排"的区别自从离开医院,专业培训磁共振,有差不多5年时间没有接触CT了.影像技术及影像医学诊断说到底是一门经验技术累积学科,这种学科的 ...
CT中的层与排的区别
转载自https://www.sohu.com/a/253957045_100012192 今天在整理一些CT的内容,自己回顾一遍,再查一些资料,重新写一遍,也是对知识的温故知新. CT术语中的&qu ...
VTK根据几个点插值成样条曲线，然后以样条曲线为中心线，生成圆柱体管道
VTK根据几个点插值成样条曲线,然后以样条曲线为中心线,生成圆柱体管道 #include "vtkActor.h" #include "vtkCamera.h" ...
SimpleITK使用深度学习识别肺癌CT DICOM数据集
肺癌数据集DICOM :https://wiki.cancerimagingarchive.net/display/Public/LIDC-IDRI 首先用SimpleITK把mhd图片读入,对每个切 ...
ITK读取CT DICOM,如何设置窗位和窗宽?
关键字: dicom itk CT 不清晰窗位窗宽碎碎念: 最近忽然要用到ITK读取DICOM,在听到任务之前,我连ITK和DICOM是个啥东西都不知道.当然,刚听到任务的时候,我脑子里也只有& ...
dicom文件_DICOM数据转成NIfTI数据
使用MRIcron工具中的dcm2nii命令,还可以可视化NIfTI数据. $ sudo apt install mricron 使用命令行操作: ## 查看帮助文档使用dcm2nii工具转换成NI ...
链路层基本问题 : 封装成帧、差错检测、流量控制
一.封装成帧 1.MAC帧类型字段 (2个字节): 用来标志上一层使用的是什么协议,以便把收到的MAC帧的数据上交给上一层的这个协议. 数据字段 (46-1500): 正式名称是MAC客户数据字段最 ...

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

CT dicom 5mm层厚插值成1mm-20201110相关推荐

最新文章

热门文章