图像数字处理部分程序

.-.例3.1-----------------------------------------------------------------I=imread('');%模拟高斯噪声J=imnoise(I,'gaussian',0,0.02);将多个图画到一个平面上subplot(1,2,1),imshow(I);subplot(1,2,2),imshow(J);K=zeros(242,308);fori=1:100J=imnoise(I,'gaussian',0,0.02);J1=im2double(J);K=K+J1;endK=K/100;figure;imshow(K);-----------------------------------------------------------------例3.2I=imread();J=imread();K=imsubtract(I,J);%实现两幅图片相减K1=255-K;figure;imshow(I);title('有噪声的图');figure;imshow(J);title('原图');figure;imshow(K1);title('提取的噪声');-----------------------------------------------------------------例3.3I=imread(); %读取一张图片J=immultiply(I,1.2); %将此图片乘以1.2K=immultiply(I,2); %将此图片乘以2subplot(1,3,1),imshow(I);subplot(1,3,2),imshow(J);subplot(1,3,3),imshow(K);-----------------------------------------------------------------例3.4moon=imread('');I=double(moon);J=I*0.43+90;K=I*0.1+90;L=I*0.01+90;moon2=uint8(J);moon3=uint8(K);moon4=uint8(L);J=imdivide(moon,moon2);K=imdivide(moon,moon3);L=imdivide(moon,moon4);subplot(2,2,1),imshow(moon);subplot(2,2,2),inshow(J,[]);subplot(2,2,3),imshow(K,[]);subplot(2,2,4),imshow(L,[]);-----------------------------------------------------------------例3.5A=zeros(128);A(40:67,60:100)=1;figure(1);imshow(A);B=zeros(128);B(50:80,40:70)=1;figure(2);imshow(B);C=and(A,B);figure(3);imshow(C);D=or(A,B);figure(4);imshow(D);E=not(A);figure(5);imshow(E);-----------------------------------------------------------------例3.6I=imread('');J=imrotate(i,60,'bilinear');K=imroatte(i,60,'bilinear','crop');subplot(1,3,1),imshow(I);subplot(1,3,2),imshow(J);subplot(1,3,3),imshow(K);-----------------------------------------------------------------例3.7I=imread('');J1=imresize(I,10,'nearest');J2=imresize(I,10,'bilinear');J3=imresize(I,10,'bicubic');figureimshow(I);title('原图');figureimshow(J1);title('最近邻法');figureimshow(J2);title('双线性插值法');figureimshow(J3);title('三次内插法');-----------------------------------------------------------------P54例4.2-----------------------------------------------------------------%读入图像I=imread('');imshow(I)%求离散傅里叶频谱J=fftshift(fft2(I));%对原始图像进行二维傅里叶变换，并将其坐标原点移到频谱图中央位置figure(2)imshow(log(abs(J)),[8,10])p58例4.4I=zeros(256,256);I(28:228,108:148)=1imshow(I)%求原始图像的傅里叶频谱J=fft2(I);F=abs(J);J1=fftshift(F);figureimshow(J1,[550])%构造原始图像I=zero(256,256)I(28:228,108:148)=1;%对原始图像进行旋转J=imrotate(I,315,'bilinear','crop');figureimshow(J)%求旋转后图像的傅里叶频谱J1=fft2(J);F=abs(J1);J2=fftshift(F);figureimshow(J2,[550]);P59例4.4.6clc;clear;I=zeros(256,256);I(8:248,110:136)=5figure(1)imshow(I)a=0.1;%原始图像的傅里叶频谱J3=fft2(I);F2=abs(J3);J4=fftshift(F2);figure(2)imshow(J4,[530])%乘以比例尺度fori=1:256forj=1:256I(i,j)=I(i,j)*a;endend%比例尺度展宽后的傅里叶频谱J2=fft2(I);F1=abs(J2);J3=fftshift(F1);figure(3)imshow(J3,[530]);p62例4.7I=imread('');figure(1);imshow(I);P=I*exp(-1);figure(2);imshow(P)P1=fftshift(fft2(P));figure(3);imshow(log(abs(P1)),[8,10])p63例4.8I=imread('');I=imnoise(I,'gaussian',0,0.01);figiure;imshow(I);P=fftshift(fft2(I));figure;imshow(log(abs(P)),[8,10])p64例4.9I=imread('');figure(1)imshow(real(I));I=I(:,:,3);fftI=fft2(I);sfftI=fftshift(fftI);RRfdp1=real(sfftI);IIfdp1=imag(sfftI);a=sqrt(RRfdp1.^2+IIfdp1.^2);a=(a-min(min(a)))/(max(max(a))-min(min(a)))*225;figure(2)imshow(real(a));b=angle(fftI);figure(3)imshow(real(b));theta=30;RR1=a*cos(theta);II1=a*sin(theta);fftI1=RR1+i.*II1;C=ifft2(fftI1)*255;figure(4);imshow(real(C));MM=150;RR2=MM*cos(angle(fftI));II2=MM*sin(angle(fftI));fftI2=RR2+i.*II2;D=ifft2(fftI2);figure(5);imshow(real(D));P67例4.10%装入图像RGB=imread('');I=rgb2gray(RGB);%画出图像imshow(I);figure(2);%进行余弦变换J=dct2(I);imshow(log(abs(J),[]);figure(3);J(abs(J)<10)=0;%进行余弦反变换K=idct2(J)/255;imshow(K);p67例4.11RGB=imread('');I=rgb2gray(RGB);I=im2double(I);T=dctmtx(8);B=blkproc(I,[88],'P1*x*P2',T,T');mask=[1111000011100000110000001000000000000000000000000000000000000000];B2=blkproc(B,[88],'P1.*x',mask);I2=blkproc(B2,[88],'P1.*xP2',T',T);figure,imshow(I);figure,imshow(B2);figure,imshow(I2);例5.1Img1=imread()figure,imshow(img1);title('originalimage');img2=imcomplement(img1);figure,imshow(img2);title('negativeimg2');例5.3I=imread('Lena.bmp');I=rgb2gray(I);K=16;H=histeq(I,K);figure,subplot(2,2,1),imshow(I,[])subplot(2,2,1),imshow(I,[]),holdonsubplot(2,2,1),hist(double(I),16),subplot(2,2,4),hist(double(H))例5.5img=rgb2gray(imread('Image.bmp'));figure;imshow(img);img_noise=double(imnoise(img,'salt&pepper',0.06));figure;imshow(img_noise,[]);img_smoothed=imfilter(img_noise,fspecial('average',5));figure;imshow(img_smoothed,[]);例5.7f=lowpassfilter([66],o.3,6);f=f/sum(f(:));I=imread('Image.bmp');I1=uint8(round(filter2(f,double(I))));其中，Butterworth低过滤波器函数为lowpassfilter.functionf=lowpassfilter(sze,cutoff,n)ifcutoff<0|cutoff>0.5error('cutofffrequencymustbebetween0and0.5');endifrem(n,1)~=0|n<1error('nmustbeaninteger>1');endiflength(sze)==1row=sze;cols=sze;elserow=sze(1);cols=sze(2);endifmod(cols,2)xrange=[-(cols-1)/2:(cols-1)/2]/(cols-1);elsexrange=[-cols/2;(cols/2-1)]/cols;endifmod(rows,2)yrange=[-(row-1)/2:(row-1)/2]/(rows-1);elseyrange=[-row/2:(row/2-1)]/rows;end[x,y]=meshgrid(xrange,yrange);radius=sqrt(x.^2+y.^2);f=ifftshift(1.0/(1.0+(radius./cutoff).^(2*n)));例5.8f=highpassfilter([77],0.3,6);f1=zeros(7);f1(4,4)=1;f=f1-f/sum(f(:));I=imread('Image.bmp');I1=filter2(f,double(I));figure,imshow(I,[])figure,imshow(I1,[])例5.9[image_0,map]=imread('Image.bmp');image_1=log(double(image_0)+1);image_2=ff