matlab 实现马赫带效应,图像上机实验.doc

图像上机实验.doc

实验一 图像处理的基本操作 一，实验目的 1．熟悉有关数字图像处理的MATLAB基本操作和编写； 2．掌握数字图像的基本读写、显示与保存； 3．能够利用MATLAB工具箱完成图像处理的基本运算； 二，实验内容 1．(1)读取图像并显示真彩色图像house.tiff； (2)读取图像并显示索引色图像canoe.tif，并转换为灰度图像； (3)读取图像并显示灰度图像rice.png，并转换为二进制图像； (4)检测内存中的图像和保存图像； 2．查看不同分辨率下图像lena.tiff的显示结果； 3．对彩色图像Peppers.png观看GRB通道各单通道图像； 4．实现马赫带效应并画出灰度直方图。 三，代码 I=imread( house.tiff ); figure(1) imshow(I) [A,map]=imread( canoe.tif ); figure(2) subplot(121) imshow(A,map) D=ind2gray(A,map); subplot(122) imshow(D) G=imread( rice.png ); figure(3) subplot(121) imshow(G) th=graythresh(G);H=im2bw(G,th); subplot(122) imshow(H) whos;imwrite B=imread( lena.tiff ); figure(4), subplot(121) imshow(B) C = imresize(B, 0.2) subplot(122) imshow(C)D=imread( Peppers.png ); r=D,r(:,:,[2 3])=0; g=D,g(:,:,[1 3])=0; b=D,b(:,:,[1 2])=0; figure(5), subplot(1,3,1),imshow(r); subplot(1,3,2),imshow(g); subplot(1,3,3),imshow(b); a=zeros(256,256); for i=1:32:256a(:,i:i+32)=i; end figure(6); subplot(121) imshow(uint8(a)); subplot(122) imhist(uint8(a)); 四，实验结果五，结果分析从运行结果可以看出，实验还是比较成功的。但是在做索引色那个图像的 时候，如果索引色和灰度图像放在同一个figure里，索引色图像显示的也是灰 色，因此我只能把索引色图像和灰度图像分开放到两个figure里。在查看不同 分辨率下图像lena.tiff的显示结果时，必须要对两个图像进行放大查看才能 看出不同。 实验二 图像变换1 一，实验目的 1．熟悉有关MATLAB工具箱中提供的图像变换函数； 2．掌握傅立叶变换和小波变换等常用的正交变换函数； 二，实验内容 1．对图像lena.png进行傅立叶变换和离散余弦变换，并分别求出其逆变换后 重构图像的均方误差； 2．实现图像lena.png的两层小波分解，观察分解系数并重构，求重构图像误 差值； 3．观察图像IM1.BMP经小波分解系数的方向性，实现三层小波分解，分别把 HL,LH,HH子带置零和重构，观察重构图像跟原始图像的差异。 三，代码 %1．对图像lena.png进行傅立叶变换和离散余弦变换，% 并分别求出其逆变换后重构图像的均方误差； clear; clc; I=imread( lena.png ); %傅里叶变换 figure(1); I1=fft2(I); I1=fftshift(I1); If=ifftshift(I1); If=uint8(real(ifft2(If))); subplot(231); imshow(I); subplot(232); imshow(log(abs(I1)),[]); subplot(233); imshow(If); [m,n]=size(I); %计算均方误差 for i=1:mfor j=1:ntemp(i,j)=(If(i,j)-I(i,j))^2;end end wucha1=sum(sum(temp)) %离散余弦变换 I2=dct2(I); I3=fftshift(I2); Id=idct2(I2); Id=uint8(real(Id)); subplot(234); imshow(I); subplot(235); imshow(log(abs(I3)),[]); subplot(236); imshow(Id); [m,n]=size(I); %计算均方误差 for i=1:mfor j=1:ntemp(i,j)=(If(i,j)-I(i,j))^2;end end wucha2=sum(sum(temp))%% % 2．实现图像lena.png的两层小波分解，观察分解系数并重构，求重构图像误差值； X=imread( lena.png ); [c,s]=wavedec2(X,2, bior3.7 ); a1=wrcoef2( a ,c,s, bior3.7 ,1); h1=wrcoef2( h ,c,s, bior3.7 ,1); v1=wrcoef2( v ,c,s, bior3.7 ,1); d1=wrcoef2( d ,c,s, bior3.7 ,1); c1=[a1,h1;v1,d1]; c1=uint8(c1); Xr=waverec2(c,s, bior3.7 ); figure(2); subplot(131); imshow(X); subplot(132); imshow(c1); subplot(133); imshow(uint8(Xr)); [m,n]=size(X); %计算均方误差 for i=1:mfor j=1:ntemp(i,j)=(Xr(i,j)-I(i,j))^2;end end wucha3=sum(sum(temp))%% % 3．观察图像IM1.BMP经小波分解系数的方向性，实现三层小波分解， % 分别把HL,LH,HH子带置零和重构，观察重构图像跟原始图像的差异。 P=imread( C:\shiyan2\IM1.BMP ); [c,s]=wavedec2(P,3, bior3.7 ); a1=wrcoef2( a ,c,s, bior3.7 ,1); h1=wrcoef2( h ,c,s, bior3.7 ,1); v1=wrcoef2( v ,c,s, bior3.7 ,1); d1=wrcoef2( d ,c,s, bior3.7 ,1); c1=[a1,h1;v1,d1]; c1=uint8(c1); figure(3); subplot(131); imshow(P); subplot(132); imshow(c1); subplot(133); imshow(uint8(a1));