资源描述：
《彩色图片matlab傅里叶小波变换》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、彩色图片傅里叶变换、小波变换傅里叶变换x=imread('1.jpg','jpg');%读取1.jpg彩色图片，并赋值给xr(:,:)=x(:,:,1);%将三维矩阵中代表R的二维矩阵赋值给rg(:,:)=x(:,:,2);%将三维矩阵中代表G的二维矩阵赋值给gb(:,:)=x(:,:,3);%将三维矩阵中代表B的二维矩阵赋值给b对r、g、b二维矩阵进行傅里叶变换fftr=fft2(r);%二维离散傅立叶变换fftg=fft2(g);%二维离散傅立叶变换fftb=fft2(b);%二维离散傅立叶变换sfftr=fftshift(fftr);%直流分量移到频谱
2、中心sfftg=fftshift(fftg);%直流分量移到频谱中心sfftb=fftshift(fftb);%直流分量移到频谱中心%计算频谱幅值A=sqrt(real(sfftr).^2+imag(sfftr).^2);B=sqrt(real(sfftg).^2+imag(sfftg).^2);C=sqrt(real(sfftb).^2+imag(sfftb).^2);%归一化，其中乘以225或255是为了在图灰度范围内显示频谱A=(A-min(min(A)))/(max(max(A))-min(min(A)))*225;B=(B-min(min(B)))
3、/(max(max(B))-min(min(B)))*225;C=(C-min(min(C)))/(max(max(C))-min(min(C)))*225;figure(1);%设定窗口1imshow(A);%显示图像的频谱title('Rfftgraphic')%添加标题figure(2);%设定窗口2imshow(B);%显示图像的频谱title('Bfftgraphic');%添加标题figure(3);%设定窗口3imshow(C);%显示图像的频谱title('Gfftgraphic');%添加标题jpg转为灰度图像的傅里叶变换dd=rgb2gr
4、ay(x);%RGB转为灰度图像figure(4);sfftdd=fftshift(fft2(dd));DD=sqrt(real(sfftdd).^2+imag(sfftdd).^2);DD=(DD-min(min(DD)))/(max(max(DD))-min(min(DD)))*225;mshow(DD);title('jpg2grayfftgraphic');小波变换，注释略x=imread('1.jpg','jpg');%读取1.jpg彩色图片，并赋值给xr(:,:)=x(:,:,1);%将三维矩阵中代表R的二维矩阵赋值给rg(:,:)=x(:,:,
5、2);%将三维矩阵中代表G的二维矩阵赋值给gb(:,:)=x(:,:,3);%将三维矩阵中代表B的二维矩阵赋值给b%对R进行小波变换重构[c,s]=wavedec2(r,2,'sym4');sizec=size(c);fori=1:sizec(2)if(c(i)>350)c(i)=2*c(i);elsec(i)=0.5*c(i);endendxx=waverec2(c,s,'sym4');%对G进行小波变换重构[c,s]=wavedec2(g,2,'sym4');sizec=size(c);fori=1:sizec(2)if(c(i)>350)c(i)=2*
6、c(i);elsec(i)=0.5*c(i);endendyy=waverec2(c,s,'sym4');%对B进行小波变换重构[c,s]=wavedec2(b,2,'sym4');sizec=size(c);fori=1:sizec(2)if(c(i)>350)c(i)=2*c(i);elsec(i)=0.5*c(i);endendzz=waverec2(c,s,'sym4');figure(1);image(xx);title('enhancedR');%设置窗口，绘制小波处理后的R图像figure(2);image(r);title('origiona
7、lR');%设置窗口，绘制小波处理前的R图像figure(1);image(yy);title('enhancedG');%设置窗口，绘制小波处理后的G图像figure(2);image(g);title('origionalG');%设置窗口，绘制小波处理前的G图像figure(1);image(zz);title('enhancedB');%设置窗口，绘制小波处理后的B图像figure(2);image(b);title('origionalB');%设置窗口，绘制小波处理前的B图像%将小波处理后的二维矩阵重新赋值给三维矩阵hh(1:size(xx,1)
8、,1:size(xx,2),1)=xx;hh(1:s
