数字图像实验报告.docx

实验报告专业名称班级学号学生姓名实验一 图像的基本运算实验内容1、程序：I=imread(lena8.jpg);figure;subplot(2,3,1);imshow(I);title(原图);J=imadjust(I,0.3;0.6,0.1;0.9); %设置灰度变换的范围subplot(2,3,2);imshow(J);title(线性扩展);I1=double(I); %将图像转换为double类型I2=I1/255; %归一化此图像C=2;K=C*log(1+I2); %求图像的对数变换subplot(2,3,3);imshow(K);title(非线性扩展);M=im2bw(I,0.5);M=M;%M=255-I; %将此图像取反%Figuresubplot(2,3,4);imshow(M);title(灰度倒置);N1=im2bw(I,0.4); %将此图像二值化，阈值为0.4N2=im2bw(I,0.7); %将此图像二值化，阈值为0.7subplot(2,3,5);imshow(N1);title(二值化阈值0.4);subplot(2,3,6);imshow(N2);title(二值化阈值0.7);结果：2、程序：I=imread(hough.bmp);%I=rgb2gray(I);J=imread(rice.bmp);I=im2double(I); %将图像转换成double型J=im2double(J);K=I+0.3*J; %两幅图像相加subplot(1,3,1);imshow(I);title(物图);subplot(1,3,2);imshow(J);title(背景图);subplot(1,3,3);imshow(K);title(相加后的图);imwrite(K,lena1.jpg);结果；3、程序：A=imread(lena1.jpg);B=imread(rice.bmp);A=im2double(A); B=im2double(B);C=A-0.3*B;subplot(1,3,1);imshow(A);title(混合图);subplot(1,3,2);imshow(B);title(背景图);subplot(1,3,3);imshow(C);title(分离后的图);结果：4、程序：A=imread(rice.bmp);%A=rgb2gray(A);A=im2double(A);subplot(1,2,1);imshow(A);title(原图);B=zeros(256,256);B(40:200,40:200)=1;K=A.*B;subplot(1,2,2);imshow(K);title(局部图);结果：5、程序：A=imread(lena8.jpg);Bl=imresize(A,1.5); %比例放大1.5倍,默认的采用的是最近邻法进行线性插值B2=imresize (A, 420 384); %非比例放大至420:384 Cl=imresize (A, 0 . 7) ; %比例缩小0.7倍C1=imresize(A, 0.7) ;C2=imresize(A, 150 180) ; %非比例缩小到150:180figure;subplot(2,2,1);imshow(Bl);title(比例放大图);subplot(2,2,2);imshow(B2);title(非比例放大图);subplot(2,2,3);imshow(C1);title(比例缩小图);subplot(2,2,4);imshow(C2);title(非比例缩小图);结果：实验二 图像的变换实验内容1、 程序：I=imread(1.bmp);%I=imread(LENA.JPG);%imshow(I);figure();subplot(3,2,1);imshow(real(I);I=I(:,:,3);fftI=fft2(I);sfftI=fftshift(fftI); %求离散傅里叶频谱%对原始图像进行二维离散傅里叶变换，并将其坐标原点移到频谱中央位置RRfdpl=real(sfftI);IIfdpl=imag(sfftI);a=sqrt(RRfdpl.2+IIfdpl.2);a=(a-min(min(a)/(max(max(a)-min(min(a)*225;subplot(3,2,2);imshow(real(a);I=imread(2.bmp);%I=imread(LENA.JPG);%imshow(I);subplot(3,2,3);imshow(real(I);I=I(:,:,3);fftI=fft2(I);sfftI=fftshift(fftI); %求离散傅里叶频谱%对原始图像进行二维离散傅里叶变换，并将其坐标原点移到频谱中央位置RRfdpl=real(sfftI);IIfdpl=imag(sfftI);a=sqrt(RRfdpl.2+IIfdpl.2);a=(a-min(min(a)/(max(max(a)-min(min(a)*225;subplot(3,2,4);imshow(real(a);I=imread(3.bmp);%I=imread(LENA.JPG);%imshow(I);subplot(3,2,5);imshow(real(I);I=I(:,:,3);fftI=fft2(I);sfftI=fftshift(fftI); %求离散傅里叶频谱%对原始图像进行二维离散傅里叶变换，并将其坐标原点移到频谱中央位置RRfdpl=real(sfftI);IIfdpl=imag(sfftI);a=sqrt(RRfdpl.2+IIfdpl.2);a=(a-min(min(a)/(max(max(a)-min(min(a)*225;subplot(3,2,6);imshow(real(a);结果：2、程序：I=zeros(256,256);I(88:168,124:132)=1; %图像范围是256*256，前一值是纵向比，后一值是横向比figure(),subplot(2,2,1);imshow(I); %求原始图像的傅里叶变换J=fft2(I);F=abs(J);J1=fftshift(F);subplot(2,2,2);imshow(J1,5 50);%对原始图像进行旋转J=imrotate(I,90,bilinear,crop);subplot(2,2,3);imshow(J);%求旋转后图像的傅里叶频谱J1=fft2(J);F=abs(J1);J2=fftshift(F);subplot(2,2,4);imshow(J2,5 50);结果：3、程序：%对 lena24.jpg文件计算二维DCT变换RGB = imread(lena24.jpg);figure(1),subplot(1,3,1);imshow(RGB);I = rgb2gray(RGB);%真彩色图像转换成灰度图像J = dct2(I);%计算二维DCT变换subplot(1,3,2);imshow(log(abs(J),);%对图像大部分能量集中在上右角处subplot(1,3,3);J(abs(J) =Step*(k-1)&IStep*k); Image(Index)=Projectm(k); end figure, imshow(Image,)end结果：3、程序：clear allclose allI=double(imread(lena8.jpg);figure,subplot(2,4,1);imshow(I,);% 1.均值低通滤波 H=fspecial(average,5); F1=double(filter2(H,I); subplot(2,4,2);,imshow(F1,);% 2 . gaussian低通滤波 H=fspecial(gaussian,7,3); F2=double(filter2(H,I); subplot(2,4,3);,imshow(F2,);% 3.增强图像-原图-均值低通滤波F3=2*I-F1;subplot(2,4,4);,imshow(uint8 (F3),);% 4.增强图像=原图-高斯低通滤波F4=2*I-F2;subplot(2,4,5);,imshow(uint8 (F4),);%5. prewitt边缘算子增强 H=fspecial(prewitt);F 5=uint8(I+filter2(H,I); subplot(2,4,6);,imshow(F5,);%6. sobel边缘算子增强 H=fspecial(sobel);F6=uint8(I + filter2(H,I); subplot(2,4,7);,imshow(F6,);结果：实验四 图像的复原实验内容1、程序：clear;close all;%1.生成含噪图像img = imread(lena8.bmp);figure,subplot(2,3,1); imshow(img);img =double(imnoise(img,salt & pepper, 0.01);subplot(2,3,2); ,imshow(img,);%2.采用均值滤波N=5; %滤波模板大小h=fspecial(average,N);I=filter2(h,img);subplot(2,3,3); ,imshow(I,)%3.中值滤波I=medfilt2(img,N N);subplot(2,3,4); ,imshow(I,)%4.最大值滤波I=ordfilt2(img,N*N,true(N);subplot(2,3,5); ,imshow(I,)%5.最小值滤波I=ordfilt2(img,1,true(N); subplot(2,3,6); ,imshow(I,)结果：2、程序：close allclear all%1.生成波纹噪声图像img = double(imread(lena8.bmp);figure,subplot(2,3,1); imshow(img,);sizec=size(img);w=0.4*2*pi;%噪声的数字频率N=2*pi/w; %噪声每一周期的采样点数img_noise=img+20*ones(sizec(1),1)*sin(w*1:sizec(2); subplot(2,3,2); ,imshow(img_noise,);%图像频谱 F0=fft2(img);F0=fftshift(F0);subplot(2,3,3); ,imshow(log(abs(F0),);F=fft2(img_noise);F=fftshift(F);subplot(2,3,4); ,imshow(log(abs(F),);%2.设计理