实验二图像直方图分析.docVIP

实验二 图像直方图分析 一．实验目的及要求 1．了解MATLAB的操作环境和基本功能。 2．掌握MATLAB中图像 3．加深理解图像增强与平滑的算法原理。 二、实验内容 （一）研究以下程序，分析程序功能；输入执行各命令行，认真观察命令执行的结果。熟悉程序中所使用函数的调用方法，改变有关参数，观察试验结果。（可将每段程序保存为一个.m文件） 1．直方图均衡化 clear all; close all % Clear the MATLAB workspace of any variables % and close open figure windows. I = imread(pout.tif); % eads the sample images ‘ pout.tif’, and stores it in imshow(I) % an array named I.display the image figure, imhist(I) % Create a histogram of the image and display it in % a new figure window. [I2,T] = histeq(I); % Histogram equalization. figure, imshow(I2) % Display the new equalized image, I2, in a new figure window. figure, imhist(I2) % Create a histogram of the equalized image I2. figure,plot((0:255)/255,T); % plot the transformation curve. imwrite (I2, pout2.png); % Write the newly adjusted image I2 to a disk file named % ‘pout2.png’. imfinfo(pout2.png) % Check the contents of the newly written file 2．直接灰度变换 clear all; close all I = imread(cameraman.tif); J = imadjust(I,[0 0.2],[0.5 1]); imshow(I) figure, imshow(J) [X,map] = imread(forest.tif); figure,imshow(X,map) I2 = ind2gray(X,map); J2 = imadjust(I2,[],[],0.5); figure,imshow(I2) figure, imshow(J2) J3 = imadjust(I2,[],[],1.5); figure, imshow(J3) help imadjust % Display the imadjust() function information. 3．空域平滑滤波（模糊、去噪） clear all; close all I = imread(eight.tif); h1 = ones(3,3) / 9; h2 = ones(5,5) / 25; I1 = imfilter(I,h1); I2 = imfilter(I,h2); figure(1), imshow(I), title(Original Image); figure(2), imshow(I1), title(Filtered Image With 3*3 ) figure(3), imshow(I2), title(Filtered Image With 5*5 ) % 加入Gaussian 噪声 J1 = imnoise(I,gaussian,0,0.005); % 加入椒盐噪声 J2 = imnoise(I,salt pepper,0.02); % 对J1、J2进行平均值平滑滤波 K1 = imfilter(J1,fspecial(average,3)); K2 = imfilter(J2,fspecial(average,3)); figure(4); subplot(2,2,1), imshow(J1) , title(gaussian); subplot(2,2,2), imshow(J2), title(salt pepper ); subplot(2,2,3), imshow(K1), title(a