图像处理技术实验指导书(201603)程序代码.docVIP

实验一 图像变换、相关运算及直方图均衡化（2学时） 实验目的与要求 熟练掌握数字图像的读取和计算机显示，熟练掌握彩色数字图像灰度化； 熟练掌握数字图像的傅里叶变换、余弦变换； 熟练掌握数字图像卷积运算 掌握数字图像相关运算定位方法 熟练掌握图像直方图的概念，熟练掌握图像直方图均衡化方法； 实验内容 读取一彩色图像并显示图像（imread、imshow）（如:“\matlab7\toolbox\images\imdemos\” 下的图像），将彩色图像转为灰度图像（rgb2gray） 将上述灰度图像进行二维傅里叶变换（fft2），并显示傅里叶谱，并把频谱原点移到中心（fftshift）。 从原图像中剪切一部分图像，并另存为子图像，把原图像和子图像转为灰度图像，进行卷积运算（imfilter），查看卷积效果。 试用相关运算法把题3中的子图像从原图像中找到位置，并在原图像中把子图像的区域用白色边框圈起来。（选作） 编程实现图像直方图求取显示（imhist），并用直方图均衡化对原图像进行灰度处理（histeq）； a=imread(11.bmp); figure(1),subplot(1,2,1),imshow(a),title(原始图像); b=rgb2gray(a); figure(2),subplot(1,2,2),imshow(b),title(灰度图); fb=fft2(b); afb=abs(fb); max_afb=max(max(afb)); afb_1=afb*255/max_afb; figure(2),imshow(afb_1),title(傅里叶幅值谱); afb_shift=fftshift(afb_1); figure(3),imshow(afb_shift),title(半波平移中心); cb=dct2(b); max_cb=max(max(cb)); cb_1=cb*255/max_cb; figure(4),imshow(cb_1),title(余弦变换); data_11=imread(11.bmp); gray_11=rgb2gray(data_11); figure(5),imshow(gray_11); data_1=imread(1.bmp); gray_1=rgb2gray(data_1); figure(6),imshow(gray_1); C = real(ifft2(fft2(gray_11,1024,1024) .* fft2(gray_1,1024,1024))); C = (C/max(max(C))*255); figure(7),imshow(C,[]); title(卷积运算); CC = real(ifft2(fft2(C,1024,1024) ./ fft2(gray_1,1024,1024))); CC2=(CC/max(max(CC))*255); figure(8),imshow(CC2,[]); title(反卷积运算); f=rgb2gray(imread(11.bmp)); figure(9),imshow(f); figure(10),imhist(f); ha=histeq(f); figure(11),imshow(ha); 实验二 图像增强（2学时） 一、实验目的与要求 熟练掌握图像的中值滤波、邻域平均、低通滤波等平滑方法； 熟练掌握图像的空域、频域锐化方法 熟练掌握图像的直方图的概念； 二、实验内容 读取一幅图像，并加入噪声(imnoise)，采用中值滤波（medfilt2）、邻域平均（imfilter）、低通滤波、等方法进行滤波处理，并比较滤波效果； 读取一幅图像，采用进行8邻域、4邻域拉普拉斯微分算子对该图像进行图像锐化并显示锐化图像，采用高通滤波器对该图像进行锐化并显示锐化图像。 读取一图像，编程实现直方图求取并进行直方图显示； 不采用系统函数实现图像的直方图统计并显示直方图； clear RGB = imread(1.bmp); figure(1),subplot(2,3,1),imshow(RGB),title(原始图像); I1=rgb2gray(RGB); subplot(2,3,2),imshow(I1),title(灰度图像); GN1=imnoise(I1,gaussian); subplot(2,3,3),imshow(GN1),title(加高斯噪声); I2=medfilt2(GN1,[3 3]); subplot(2,3,4),imshow(I2),title(中值滤波); h=ones(3,3)/9; I3=imfilter(GN1,h); subplot(2,3,5),imshow(I3),title(邻域平均)