dicom图像处理：使用Matlab裁剪并输出选定区域的CT图像

1、主要函数简介：

• dicomread(filepath)：matlab自带的读取dicom文件的库函数，filepath是待读取文件的路径；

• dicomwrite(image,filepath)：matlab自带的dicom写入函数，image是存储图像像素值的二维矩阵，filepath是待写入文件的路径；

• dicominfo(filepath)：matlab自带的读取指定路径dicom文件头信息的函数；

• imrect：创建可变形矩形，在已显示的图像中选中矩形区域，双击左键确认；

• wait(rect)：返回矩形的对角线的顶点坐标

• imcrop(image，h)：裁剪指定区域的图像，image为待裁剪图像，h是存有矩形区域对角线顶点的向量

2、代码:

%读取dicom图像并显示
I = dicomread('D:\...\IMG-0002-00260.dcm');%读取图像
metadata = dicominfo('D:\...\IMG-0002-00260.dcm');%存储信息
imagesc(I);

%在显示的图像中选定矩形区域
b = imrect;
h = wait(b);%得到矩形的对角线的顶点坐标

%剪切选定区域
o = imcrop(I,h);
imagesc(o);

%写入新的dcm文件
dicomwrite(o,'D:\...\01.dcm');