MatLab读写dicom 像素操作

　　使用Matlab读写Dicom比较简单。

一、读DICOM图像

　　主要涉及到两个函数：dicominfo和dicomread

（1）metadata = dicominfo(filename)

　　　　主要用于读取dicom的头文件

　　获取某个tag直接使用metadata.方法即可，例如metadata.SeriesInstanceUID。

　　需要修改tag也简单，例如：metadata.SeriesInstanceUID = '1.3.12.2.1107.5.1.4.73336.3000002109301112409880041433001';

（2）img= dicomread(filename)

　　　　用于读取dicom像素信息，img直接为二维矩阵，使用imshow可以直接显示图像

二、写DICOM图像

主要涉及到的函数为：dicomwrite

dicomwrite(img,newfilename,metadata);

如果需要修改dicom图像的像素信息，只需要修改img中的值即可，img(i,j) = 0;

filename=fullfile('E:\\***\\IMG-0001-00484.dcm');
info = dicominfo(filename);
img = dicomread(filename);
info.SeriesInstanceUID = '1.3.12.2.1107.5.1.4.73336.3000002109301112409880041433001';
for i = 1:512
    for j = 1:512
        if(i>240 && i < 320)
            z = 1;
        else
            img(i,j) = 0;                 
        end
    end
end    
newfilename = 'E:\\***\\IMG-0001-00161.dcm'
dicomwrite(img,newfilename,info);