DICOM图像横断面投影.docxVIP

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Share Information %%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%% Tag: Matlab, Medical Imaging, DICOM %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 标题：DICOM图像横断面投影%% 内容：对DICOM图像进行横断面投影%% 作者：Destiny%% 日期：2012-10-05%% 是否转载：否%For initialP = zeros(512,512); %To store the final result %For read all images automaticly filename = dir(*.dcm);[r c]=size(filename);%Projection for X-Rayfor a = 1:(r-1) %For each image, read it by using dicominfo() an dicomread() str = sprintf(filename(a).name); info=dicominfo(str); Y0 = dicomread(info); str = sprintf(filename(a+1).name); info=dicominfo(str); Y1 = dicomread(info); %Get the image size [m,n] = size(Y); %Maximum Intensity Projection(MIP) on transverse plane(heng duan mian) for i=1:m for j=1:n if Y1(i,j)Y0(i,j) P(i,j)=Y1(i,j); else P(i,j)=Y0(i,j); end end endendfigure,imshow(P,[])%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%% Contact: 1139904786@ %%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%% Share and Access - Destiny %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%