醫(yī)學(xué)圖像處理實(shí)驗(yàn)

1、實(shí)驗(yàn)一yq1I=imread('cameraman.tif');%讀黑白圖像subplot(2,2,1);imshow(I) %顯示圖像subplot(2,2,2);imhist(I) %顯示直方圖J=imadjust(I,0.02 0.7,0 1);%對(duì)比度增強(qiáng)subplot(2,2,3);imshow(J)subplot(2,2,4);imhist(J)I1=imresize(I,0.5);imview(I1)%縮小I2=imresize(I,1.5);imview(I2)%放大I3=imrotate(I,45,'bilinear','crop

2、9;);imview(I3)%旋轉(zhuǎn)45°%原圖、直方圖 對(duì)比度增強(qiáng)、直方圖%縮小%放大%旋轉(zhuǎn)45°yq2I=imread('pears.png');imshow(I);I1=rgb2gray(I);%把彩色圖像轉(zhuǎn)換成灰度圖像figure,imshow(I1);info= imfinfo('pears.png')%查詢(xún)文件信息imwrite(I1,'D:yq小小.png'); %寫(xiě)圖像info = Filename:'C:MATLAB7toolboximagesimdemospears.png'FileModDa

3、te: '03-May-2003 13:53:58'FileSize: 554554Format: 'png'FormatVersion: Width: 732Height: 486BitDepth: 24ColorType: 'truecolor'FormatSignature: 137 80 78 71 13 10 26 10Colormap: Histogram: InterlaceType: 'none'Transparency: 'none'SimpleTransparencyData: Backgrou

4、ndColor: RenderingIntent: Chromaticities: Gamma: XResolution: YResolution: ResolutionUnit: XOffset: YOffset: OffsetUnit: SignificantBits: ImageModTime:'20 Feb 2003 20:53:33 +0000'Title: Author: Description: Copyright: 'Copyright Corel'CreationTime: Software: Disclaimer: Warning: Sour

5、ce: Comment: OtherText: yq3I,map=imread('trees.tif');imshow(I,map)I1=ind2gray(I,map);%把索引色轉(zhuǎn)換成灰度圖像figure,imshow(I1,map)%索引色圖像%灰度圖像實(shí)驗(yàn)二yq4%對(duì)比度調(diào)整I=imread('pout.tif');subplot(2,2,1);imshow(I);subplot(2,2,2);imhist(I);J=imadjust(I,0.3 0.7,0 1);subplot(2,2,3);imshow(J);subplot(2,2,4);imhist

6、(J);%原圖、直方圖 對(duì)比度增強(qiáng)、直方圖yq5%直方圖均衡化：I=imread('pout.tif');imshow(I);figure,imhist(I);J,T=histeq(I,64);figure,imshow(J);figure,imhist(J); yq6%線性濾波的MATLAB實(shí)現(xiàn)I=imread('kids.tif');I1=imnoise(I,'salt & pepper',0.02);%加椒鹽噪聲K1=filter2(fspecial('average',3),I1)/255;K2=filter2(f

7、special('average',5),I1)/255;K3=filter2(fspecial('average',7),I1)/255;subplot(2,2,1);imshow(I1); title('噪聲圖像');subplot(2,2,2);imshow(K1);title('3×3');subplot(2,2,3);imshow(K2); title('5×5');subplot(2,2,4);imshow(K2); title('7×7');yq7%中值濾

8、波MATLAB實(shí)現(xiàn)I=imread('eight.tif');imshow(I);J=imnoise(I,'gaussian',0,0.02);figure;imshow(J);K1=medfilt2(J,3,3);K2=medfilt2(J,5,5);K3=medfilt2(J,7,7);figure,imshow(K1);figure,imshow(K2);figure,imshow(K2);%原圖%加高斯噪聲實(shí)驗(yàn)4實(shí)驗(yàn)四：圖像的配準(zhǔn)與融合一、 實(shí)驗(yàn)?zāi)康?1)熟悉MATLAB圖像處理工具箱的使用方法(2)掌握基于特征點(diǎn)的圖像配準(zhǔn)的過(guò)程。(3)掌握常用的圖像融

9、合方. 二、實(shí)驗(yàn)的主要儀器設(shè)備(1)微型計(jì)算機(jī)(2) MATLAB軟件（安裝圖像處理工具箱）(3) 參考MRI圖像與待配準(zhǔn)的CT圖像三、實(shí)驗(yàn)原理 圖像配準(zhǔn)指的是將同一場(chǎng)景的兩幅或多幅圖像進(jìn)行對(duì)準(zhǔn)。一個(gè)典型的應(yīng)用是，將一幅圖像(稱(chēng)為基準(zhǔn)圖像)作為其他圖像(稱(chēng)為輸入圖像)的參照進(jìn)行比較。圖像配準(zhǔn)的目的是，通過(guò)對(duì)輸入圖像進(jìn)行空間變換，使輸入圖像與基準(zhǔn)圖像對(duì)準(zhǔn)。使用點(diǎn)映射的圖像配準(zhǔn)包括以下步驟： 將圖像讀入到MATLAB空間； 指定圖像中的成對(duì)控制點(diǎn)； 保存控制點(diǎn)對(duì)； 指定要使用的變換類(lèi)型，并根據(jù)控制點(diǎn)對(duì)推測(cè)參數(shù)。 對(duì)沒(méi)有配準(zhǔn)的圖像進(jìn)行變換，使之對(duì)準(zhǔn)。四、試驗(yàn)內(nèi)容1、讀入圖像M= imread(

10、9;C:MATLAB7toolboximagesimdemosct.jpg');figure,imshow(M);N= imread('C:MATLAB7toolboximagesimdemosmri.jpg'); figure,imshow(N);2、選擇控制點(diǎn)cpselect(M,N);3、將控制點(diǎn)保存到MATLAB工作區(qū)間input_points =305.8058 208.6547 236.7410 256.5396 246.8705 315.4748 262.5252 395.58994、調(diào)整控制點(diǎn)的位置和指定變換類(lèi)型計(jì)算參數(shù)input_points_corr

11、=cpcorr(input_points,base_points,M,N);mytform=cp2tform(input_points_corr,base_points,'linear conformal');info=imfinfo('C:MATLAB7toolboximagesimdemosmri.jpg');5、配準(zhǔn)圖像，并顯示registered=imtransform(unregistered,mytform,'XData',1 info.Width,'YData',1 info.Height);figure,imsho

12、w(registered),title('配準(zhǔn)后圖像');6、融合M1=imread('C:MATLAB7toolboximagesimdemosmri.jpg'); %讀入圖像M2=imread('C:MATLAB7toolboximagesimdemoshh.jpg');m1 n1=size(M1); %判斷圖像大小是否一致m2 n2=size(M2);if (m1 = m2) | (n1 = n2)error('Input images are not of same size');end;M1=im2double(M1); %數(shù)據(jù)類(lèi)型轉(zhuǎn)換M2=im2double(M2);M3=0.5*M1+0.5*M2; %圖像加權(quán)融合M3=im2uint8(M3);M4=0.3*M1+0.7*M2;M4=im2uint8(M4);subplot(2,2,1);imshow(M1),title('