數(shù)字圖像處理上機(jī)大作業(yè)

1、解：（1）源代碼：close all;clear all;clc;%產(chǎn)生如圖所示二值圖像%image1=zeros(256,256);for x=64:192 for y=112:144 image1(x,y)=1; endendfigure(1);subplot(1,2,1);imshow(image1);title('題圖所示圖像');I=fft2(image1);Image1=log(abs(I);subplot(1,2,2)imshow(Image1,-3,5);title('圖像幅度譜');（2）源代碼：%for x=1:256 for y=1:256

2、 f2(x,y)=(-1)(x+y)*image1(x,y); endendfigure(2);subplot(1,2,1);imshow(f2);title('改變后圖像');I2=fft2(f2);Image2=log(abs(I2);subplot(1,2,2)imshow(Image2,-3,5);title('圖像幅度譜');（3）源代碼：%f3=imrotate(image1,90);figure(3)subplot(1,3,1);imshow(f3);title('改變后的圖像');I3=fft2(f3);Image3=log(ab

3、s(I3);subplot(1,3,2)imshow(Image1,-3,5);title('原圖像幅度譜');subplot(1,3,3)imshow(Image3,-3,5);title('改變后圖像幅度譜');（4）源代碼：%f4=f3;for x=1:256 for y=1:256 f5(x,y)=image1(x,y)+f3(x,y); endendfigure(4)subplot(1,3,1);I5=fft2(f5);Image5=log(abs(I5);imshow(Image5,-3,5);title('F5圖像幅度譜');sub

4、plot(1,3,2)imshow(Image1,-3,5);title('F1圖像幅度譜');subplot(1,3,3)imshow(Image3,-3,5);title('F4圖像幅度譜');（5）源代碼：%for x=1:256 for y=1:256 f6(x,y)=f2(x,y)+f3(x,y); endendfigure(5)subplot(1,3,1);I6=fft2(f6);Image6=log(abs(I6);imshow(Image6,-3,5);title('F6圖像幅度譜');subplot(1,3,2)imshow(I

5、mage2,-3,5);title('F2圖像幅度譜');subplot(1,3,3)imshow(Image3,-3,5);title('F3圖像幅度譜');figure(6)subplot(1,2,1);imshow(Image5,-3,5);title('F5圖像幅度譜');subplot(1,2,2)imshow(Image6,-3,5);title('F6圖像幅度譜');解：源代碼：close all;clear all;clc;%產(chǎn)生P104頁(yè)圖ͼ%x=zeros(16,16);x(1:8,1:8

6、)=1;x(9:16,9:16)=1;y=repmat(x,4,4);figure(1);subplot(1,3,1)imshow(y);title('原圖像');%3*3濾波%subplot(1,3,2)h1=0.125*1,1,1;1,0,1;1,1,1;K1=conv2(y,h1);%3*3平均濾波imshow(K1);title('3*3平均濾波');K2=medfilt2(y,3 3,'symmetric');%3*3中值濾波¨subplot(1,2,2)imshow(K2);title('3*3中值濾波')解

7、：源代碼：close all;clear all;clc;%產(chǎn)生P104頁(yè)圖%image=zeros(256,256);for i=1:9 image(29:228,(27+i+24*(i-1):27+(8*i)+17*(i-1)=255;endfigure(1);imshow(image);subplot(2,4,1)imshow(image);title('原圖像');%加入噪聲%J1=imnoise(image,'gaussian',0.05);%加入高斯白噪聲J2=imnoise(image,'salt & pepper',0.0

8、5);%加入椒鹽噪聲subplot(2,4,2)imshow(J1);title('加入高斯白噪聲');subplot(2,4,6)imshow(J2);title('加入椒鹽噪聲');%3*3濾波%subplot(2,4,3)h1=0.125*1,1,1;1,0,1;1,1,1;KJ11=conv2(J1,h1);%3*3平均濾波imshow(KJ11);title('3*3平均濾波');KJ12=medfilt2(J1,3 3,'symmetric');%3*3椒鹽濾波subplot(2,4,4)imshow(KJ12,);t

9、itle('3*3椒鹽濾波')subplot(2,4,7)KJ21=conv2(J2,h1);%3*3平均濾波imshow(KJ21);title('3*3平均濾波');KJ22=medfilt2(J2,3 3,'symmetric');%3*3椒鹽濾波subplot(2,4,8)imshow(KJ22,);title('3*3椒鹽濾波')解：源代碼：close all;clear all;clc;image=imread('E:lesson數(shù)字圖像處理image3.jpg');%讀入圖像I=rgb2gray(im

10、age);%轉(zhuǎn)化為灰度圖像figure(1)subplot(2,3,1)imshow(I);title('原圖像');IR1=edge(I,'roberts');%roberts邊緣檢測(cè)subplot(2,4,2)imshow(IR1);title('roberts邊緣檢測(cè)');IP1=edge(I,'prewitt');%prewitt邊緣檢測(cè)subplot(2,4,3)imshow(IP1);title('prewitt邊緣檢測(cè)');IS1=edge(I,'sobel');%sobel邊緣檢測(cè)s

11、ubplot(2,4,4)imshow(IS1);title('sobel邊緣檢測(cè)');H2=fspecial('prewitt');IP2=filter2(H2,I);%prewitt銳化圖像subplot(2,4,7)imshow(IP2,);title('prewitt銳化圖像');H3=fspecial('sobel');IS2=filter2(H3,I);%sobel銳化圖像subplot(2,4,8)imshow(IS2,);title('sobel銳化圖像');解：源代碼：close all;clear all;clc;I=imread('E:lesson數(shù)字圖像吃力image1.jpg');%讀入圖像image=rgb2gray(I);%圖像灰度化figure(1)subplot(2,2,1)imshow(I);%顯示原圖像title('原圖像');subplot(2,2,2)imshow(image);%顯示灰度圖像title('灰