数字图像处理实验指导书.pptVIP

%Matlab中灰度变换法伪彩色图象的实现clear;I=imread(高对比度：精美高清风光壁纸5.jpg);%I=rgb2gray(I);subplot(1,2,1);imshow(I);title(原始灰度图像);I=double(I);[MN]=size(I);L=256;fori=1:Nforj=1:MifI(j,i)L/4R(j,i)=0;G(j,i)=4*I(j,i);B(j,i)=L;elseifI(j,i)=L/2R(j,i)=0;G(j,i)=L;B(j,i)=-4*I(j,i)+2*L;elseifI(j,i)=3*L/4R(j,i)=4*I(j,i)-2*L;G(j,i)=L;B(j,i)=-0;elseR(j,i)=L;G(j,i)=-4*I(j,i)+4*L;B(j,i)=-0;endendendOUT=zeros(M,N,3);OUT(:,:,1)=R;OUT(:,:,2)=G;OUT(:,:,3)=B;OUT=uint8(OUT);subplot(1,2,2);imshow(OUT);title(灰度分级伪彩色图像);图3-6“频域伪彩色处理法”：读入一幅灰度图像，进行FFT变换，分别进行低通滤波、带通滤波、高通滤波，得到三个独立的频率分量，人为地将它们分别归属为R、G、B分量，再分别将R、G、B分量进行快速傅立叶逆变换，得到R、G、B三个位平面，最后将三个位平面合成为一幅伪彩色图像。%Matlab频域伪彩色图象的实现clear;I=imread(风景04.jpg);I=rgb2gray(I);I=imadjust(I);subplot(1,2,1);imshow(I);title(原始灰度图像);I=double(I);[MN]=size(I);F=fft2(I);fftshift(F);Red_cut=5;Green_cut=205;Blue_center=100;Blue_width=200;Blue_u0=10.5;Blue_v0=10.3;foru=1:Mforv=1:ND(u,v)=sqrt(u^2+v^2);Red_L(u,v)=1/(1+(sqrt(2)-1)*(D(u,v)/Red_cut)^2);Green_H(u,v)=1/(1+(sqrt(2)-1)*(Green_cut/D(u,v))^2);Blue_D(u,v)=sqrt((u-Blue_u0)^2+(v-Blue_v0)^2);Blue_H(u,v)=1-1/(1+Blue_D(u,v)*Blue_width/((Blue_D(u,v))^2-(Blue_center)^2)^2);endEndRed=Red_L.*F;Red_color=ifft2(Red);Green=Green_H.*F;Green_color=ifft2(Green);Blue=Blue_H.*F;Blue_color=ifft2(Blue);Red_color=real(Red_color)/256;Red_color=imadjust(Red_color);Green_color=real(Green_color)/256;Green_color=imadjust(Green_color);Blue_color=real(Blue_color)/256;Blue_color=imadjust(Blue_color);OUT=zeros(M,N,3);OUT(:,:,1)=Red_color;OUT(:,:,2)=Green_color;OUT(:,:,3)=Blue_color;OUT=abs(OUT);subplot(1,2,2);imshow(OUT);title(频率域伪彩色图像);图3-7读入一幅图像，对其作DWT分解。Matlab下对