实验一图像处理的基本操作一,实验目的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;imwriteB=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(:,:,[23])=0;g=D,g(:,:,[13])=0;b=D,b(:,:,[12])=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);fori=1:32:256a(:,i:i+32)=i;endfigure(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);%计算均方误差fori=1:mforj=1:ntemp(i,j)=(If(i,j)-I(i,j))^2;endendwucha1=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);%计算均方误差fori=1:mforj=1:ntemp(i,j)=(If(i,j)-I(i,j))^2;endendwucha2=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);%计算均方误差fori=1:mforj=1:ntemp(i,j)=(Xr(i,j)-I(i,j))^2;endendwucha3=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...