资源描述:
《图像处理技术实验指导书(201603)程序代码》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、实验一图像变换、相关运算及直方图均衡化(2学时)一、实验目的与要求1、熟练掌握数字图像的读取和计算机显示,熟练掌握彩色数字图像灰度化;2、熟练掌握数字图像的傅里叶变换、余弦变换;3、熟练掌握数字图像卷积运算4、掌握数字图像相关运算定位方法5、熟练掌握图像直方图的概念,熟练掌握图像直方图均衡化方法;二、实验内容1、读取一彩色图像并显示图像(imread、imshow)(如:“matlab7toolboximagesimdemos”下的图像),将彩色图像转为灰度图像(rgb2gray)2、
2、将上述灰度图像进行二维傅里叶变换(fft2),并显示傅里叶谱,并把频谱原点移到中心(fftshift)。3、从原图像中剪切一部分图像,并另存为子图像,把原图像和子图像转为灰度图像,进行卷积运算(imfilter),查看卷积效果。4、试用相关运算法把题3中的子图像从原图像中找到位置,并在原图像中把子图像的区域用白色边框圈起来。(选作)5、编程实现图像直方图求取显示(imhist),并用直方图均衡化对原图像进行灰度处理(histeq);a=imread('11.bmp');figure(1),sub
3、plot(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(af
4、b_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),imsh
5、ow(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('反
6、卷积运算');f=rgb2gray(imread('11.bmp'));figure(9),imshow(f);figure(10),imhist(f);ha=histeq(f);figure(11),imshow(ha);实验二图像增强(2学时)一、实验目的与要求1、熟练掌握图像的中值滤波、邻域平均、低通滤波等平滑方法;2、熟练掌握图像的空域、频域锐化方法3、熟练掌握图像的直方图的概念;二、实验内容读取一幅图像,并加入噪声(imnoise),采用中值滤波(medfilt2)、邻域平均(imfi
7、lter)、低通滤波、等方法进行滤波处理,并比较滤波效果;1、读取一幅图像,采用进行8邻域、4邻域拉普拉斯微分算子对该图像进行图像锐化并显示锐化图像,采用高通滤波器对该图像进行锐化并显示锐化图像。2、读取一图像,编程实现直方图求取并进行直方图显示;1、不采用系统函数实现图像的直方图统计并显示直方图;clearRGB=imread('1.bmp');figure(1),subplot(2,3,1),imshow(RGB),title('原始图像');I1=rgb2gray(RGB);subplot
8、(2,3,2),imshow(I1),title('灰度图像');GN1=imnoise(I1,'gaussian');subplot(2,3,3),imshow(GN1),title('加高斯噪声');I2=medfilt2(GN1,[33]);subplot(2,3,4),imshow(I2),title('中值滤波');h=ones(3,3)/9;I3=imfilter(GN1,h);subplot(2,3,5),imshow(I3),title('邻域平均');I4=imfilter(G