图像处理实验指导书

1、实验图像显示方RGB= 实验图像显示方RGB= loadclown bitmap = imread(clown.bmp,bmp); imwrite(bitmap,clown.png,png);,notruesize%创建仅包含 0/1 的双精度图%逻辑标志置为 figure,imshow(BW,100;00J=filter2(12;-1-load%装载一幅图RGB 图像的显示 %保存为J=filter2(12;-1-load%装载一幅图RGB 图像的显示 %保存为bmp 文% 327mri.tiffor frame=1:27 mri(:,:,:,frame),map=imread(mri.ti

2、f,frame);%36RGB 图像的某一个颜色分量？实验图像运1. %转换为双精实验图像运1. %转换为双精度类%转换为8J=I*1.5-%转换为8backgroundimopen(I,strel(disk,15);估计背景图像 figure, imshow(background);figure,%backgroundimopen(I,strel(disk,15);估计背景图像 figure, imshow(background);figure,%从原始图像中减去背景图读入图像ic.tif1.5 0.5 e(I,-I2=imcrop(I,30601204. e(I,-I2=imcrop(I,

3、30601204. *读入图像tire.tifnlfilterblkproc对图像进行滑动邻域操作和分I2=nlfilter(I,3 3,f); %构造复合函%滑动邻域操8(round(mean2(x)*ones(size(x); %构造复合函I2=blkproc(I,8%滑动邻域操实验图像变换与滤波器设Radon实验图像变换与滤波器设Radonfigure,imshow(F2,-15,notruesize);colormap(jet); t2(f,256,256); %256256 矩阵figure,imshow(log(abs(F),-1 %将图像频谱中心由矩阵原点移至矩阵中figure,

4、imshow(log(abs(F2),-1 figure,imshow(Cthresh-离散余弦变换dct2对图像autumn.tifDCTI=rgb2gray(RGB转换为灰度图像 DCT 10 idct2 figure,imshow(K,0I=rgb2gray(RGB转换为灰度图像 DCT 10 idct2 figure,imshow(K,0255) %转换为灰度图%舍弃系 DCT1110000110000010000000000000000000000000000000000000000B2=blkproc(B,8 I2=blkproc(B2,8 3. Radon RadonP=phan

5、tom(256); %创建256 灰度级大脑图 %18 个P=phantom(256); %创建256 灰度级大脑图 %18 个投%36 个投%90 个投使用逆变换重构图像 %R1 重构图%R2 重构图%R3 重构图4. FIRFIR滤波器(p96)b=remez(10,00.40.61,110plot(w/pi-figure,freqz2(h,32 %阶次,频率向量,%x 0 处对应曲线中实验四 形态学操作与空间变实验四 形态学操作与空间变%正方形结构元改变上述结构元素类型(如：line,diamond,disk 等)%线型结构元se=strel(rectangle,4030); %选择适当

6、大小的矩形结构元%先腐蚀,删除较细的直%再膨胀,恢复矩形的大se=strel(rectangle,40 %开启操 置结构元素大小为4 3,同时观察形态开启(imopen)与闭合(imclose)se=strel(rectangle,43); %开启操%闭合操pearlite.tif（se=strel(disk,5);pearlite.tif（se=strel(disk,5);AimregionalmaxA=1010 10101010 1010 101013 13 1310 10 1110 11 1013 13 1310 10 1011 10 1013 13 1310 1011 10 11101

7、0 10 1010 10 1010 10 1011 10 1010 18 1818 10 1010 10 1110 18 1818 10 1010 11 1010 18 1818 10 1011 10 1110 10 1010 10 1010 10 1010 10 1110 10 利用函数imextendedmax26. *创建一幅图像,求其距离矩阵。 lims=floor(center1-1.2*radius)ceil(center2+1.2*radius); figure,imshow(D,) figure,imshow(D1,lims=floor(center1-1.2*radius)c

8、eil(center2+1.2*radius); figure,imshow(D,) figure,imshow(D1,0000000111000011100000001110000111010000001100000%读入待匹配图像和标准图unregistered=imread(westconcordaerial.png); orthophoto=imread(westconcordorthophoto.png); %读入待匹配图像和标准图unregistered=imread(westconcordaerial.png); orthophoto=imread(westconcordorth

9、ophoto.png); s%使用相关性进一步协调控制点%变换未匹配的图实验图像增强与复实验图像增强与复debye1.tif flowers.tifrice.tifimadjustJ=imadjust(I,0.150.9,0J=imadjust(I,00.2,0.51); J=imadjust(I,00.2,0.51); 利用函数 imnoise 给图像eight.tifI2=imnoise(I,salt & pepper);(gaussian)噪声和椒盐(salt J=imnoise(I,salt&K1=medfilt2(J,33);3*3, *对图像saturn.tif采用K1=medfi

10、lt2(J,33);3*3, *对图像saturn.tif采用对图像flowers.tifPSFPSF%剪切图%运动位%运动角%运动%均值滤波J=imnoise(motionblur,salt&pepper); deconvwnrlen=30;theta=10;%真实J=imnoise(motionblur,salt&pepper); deconvwnrlen=30;theta=10;%真实%PSF %PSF *deconvwnr *设置信噪比和相关函数的维纳滤波复原(p179 9.12)%计算信噪%设置信噪比的复原结%噪声能量谱密%噪声自相关函%图像能量谱密%图像自相关函%设置%设置自相关函数的复原结实验图像分析与理解及图像压DCT实验图像分析与理解及图像压DCTedge函数对图像rice.tifedgecanny方法。优点roicolor, roifilt2 roifilld c=222272300270c=222272300270221121*p对图像alumgrns.tif subplot(231),imshow(afm),title(Original % subplot(232),imshow(Itop, ),title(Tophat % subplot(233),imshow(Ibot, ),title(Bottomhat Image); ;% % % %