数字图像处理实验报告(20211125232055)

1、数字图像处理实验报告姓名：高爱军学号： 20091003056 指导老师：傅华明实验地点：信息楼201 实验一一实验目的掌握 matlab 的基本操作， 在命令窗口或建立m 文件进行简单图像处理， 学会读入图像并对图像矩阵进行操作二实验内容1.计算图像的灰度中值。2.对图像进行直方图规定化。3.对图像施加伽马噪声并进行中值滤波法（3*3 窗口）。三实验步骤1.计算图像的灰度中值。使用 imread 函数读取图像，再用size 函数获取图像矩阵的高度和宽度，用变量height1 和 width1 保存。循环比较找出最大值和最小值，求二者的平均得中值或将全部灰度按从小到大排序求其中间的值2.对图像

2、进行直方图规定化。使用 imread 读入图像并显示原图，用grayxform 函数进行图像的规定化。3.对图像施加伽马噪声并进行中值滤波法（3*3 窗口）。使用 imnoise 对读入的图像添加椒盐噪声，显示出加噪声后的图像用medfilt 进行中值滤波再显示出出来完后的图像。将原图和加了噪声、 处理后的图像显示在同一窗口，便于观察分析结果。四实验代码和结果1.计算图像的灰度中值。x=imread(p09.tif); x=sort(x(:); size(x) x(480*200,1) median(x) 原图：求中值的代码。clear all; clc; close all; i1 = im

3、read(p09.tif); height1,width1,channel1 = size(i1); %求图像高度和宽度min=1; %求图像灰度中值max=0; mid=0; for i=1:height1 for j=1:width1 if i1(i,j)max max=i1(i,j); end if i1(i,j)min min=i1(i,j); end end end max=max min=min mid=(max+min)/2 求得中值为 125. 2.对图像进行直方图规定化clear all i=imread(p10.tif); j=histeq(i,32); counts,x=

4、imhist(j); q=imread(p10.tif); figure; imshow(q); title(原图像 ); figure; imhist(q); title(原图像直方图 ); m=histeq(q,counts); figure; imshow(m); title(直方图规定化后的图像 ); figure; imhist(m); title(规定直方图 ); 3.对图像施加伽马噪声并进行中值滤波法（3*3 窗口）。img=imread(p09.tif); subplot(2,2,1);imshow(img); img_noise=double(imnoise(img,salt

5、 & pepper,0.06); %添加椒盐噪声subplot(2,2,2),imshow(img_noise,); img_mean=imfilter(img_noise,fspecial(average,5); subplot(2,2,3);imshow(img_mean,); img_median=medfilt2(img_noise); %进行中值滤波，默认窗口为3*3 subplot(2,2,4);imshow(img_median),); 运行结果如下：实验二一实验目的学会设计简易图像并用进行傅立叶变换的性质的验证。二实验内容1.设计图象4*30*30/256*256 ，主

6、对角线排列。2. 利用现有的离散傅立叶变换、离散余弦变换、 walsh-hadamard 变换对同一图象实施变换，比较三种变换所得到的频谱3. 根据窗函数法设计一个高通滤波器，对两图象（f1(x,y) 为 30*30/256*256 的图象；f2(x,y)=p2图象）进行低通滤波，观察分析空域图象和频谱分布的变化。三实验步骤，代码与结果1.设计图象4*30*30/256*256 ，主对角线排列。clear all m=zeros(256,256); for i=38:68 for j=38:68 m(i,j)=1; end end for i=88:118 for j=88:118 m(i,j

7、)=1; end end for i=138:168 for j=138:168 m(i,j)=1; end end for i=188:218 for j=188:218 m(i,j)=1; end end n=fft2(m); n=fftshift(n); subplot(1,2,1);imshow(m) title(图像) subplot(1,2,2);imshow(log(abs(n) title(图像) 【运行结果图像】2. 利用现有的离散傅立叶变换、离散余弦变换、 walsh-hadamard 变换对同一图象实施变换，比较三种变换所得到的频谱clear m=zeros(256,256); for i=98:158 for j=98:158 m(i,j)=1; end end n=fft2(m); n=fftshift(n); subplot(1,3,1);imshow(log(abs(n