实验二图像锐化及伪彩色处理..docVIP

实验二 ?图像锐化及伪彩色处理 有的图片为三维的（MATLAB的workspace窗口可以看到），无法生成直方图，可以读入图像后用下面的命令转换为二维图像 I=I（：，：，1）；  1.能够对单色图像进行伪彩色处理，能够分析彩色图像。 2.能够自行评价各主要算子在无噪声条件下和噪声条件下的锐化性能。完成规定图像的处理并要求正确评价处理结果，能够从理论上作出合理的解释。 三、实验仪器 PC 一台，MATLAB软件 。 四、实验内容 (一) 试分别用roberts、sobel 、Prewitt、log四种算子编写MATLAB程序对Miss.bmp图像及添加高斯噪声的Miss图像进行锐化，提取边缘信息。 要求： 1. 显示原图像和对其锐化后的各个图像，并分别为之命名为： “原始图像”、“Robert算子锐化结果”、“sobel算子锐化结果”、“Prewitt算子锐化结果”、“高斯－拉普拉斯算子锐化结果” 2. 显示添加噪声的图像及对其锐化后的各个图像，并分别为之命名为： “添加噪声后的图像”、“Robert算子锐化结果”、“sobel算子锐化结果”、“Prewitt算子锐化结果”、“高斯－拉普拉斯算子锐化结果” 3.比较几种方法的锐化效果，并分析原因 所用到的函数说明如下： 1. 采用边缘算子分割图像函数为： [g,t]=edge(image,method,threshold,direction) image:输入的图像 method：采用的方法类型，有roberts、sobel 、Prewitt、log threshold：阈值，通常采用默认值，表示为[] direction: 所寻找边缘的方向，常用both ； g：返回的二值图像 2.title(图像名字); %此函数可为图像命名 3.添加高斯噪声 I=imnoise(image,gaussian,0.02); （1）I=imread(Miss.bmp); [g1,t]=edge(I,roberts,[],both); [g2,t]=edge(I,sobel,[],both); [g3,t]=edge(I,Prewitt,[],both); [g4,t]=edge(I,log,[],both); subplot(231);imshow(I);title(原始图像); subplot(232);imshow(g1);title(reborts算子锐化结果); subplot(233);imshow(g2);title(sobel算子锐化结果); subplot(234);imshow(g3);title(prewitt算子锐化结果); subplot(235);imshow(g4);title(log算子锐化结果); （2） I=imread(Miss.bmp); J=imnoise(I,gaussian,0.02); [g5,t]=edge(I,roberts,[],both); [g6,t]=edge(I,sobel,[],both); [g7,t]=edge(I,Prewitt,[],both); [g8,t]=edge(I,log,[],both); subplot(231);imshow(J);title(添加高斯噪声图像 ); subplot(232);imshow(g5);title(reborts算子锐化结果‘); subplot(233);imshow(g6);title(sobel算子锐化结果); subplot(234);imshow(g7);title(prewit算子锐化结果); subplot(235);imshow(g8);title(log??算子锐化结果); (二) 运行下列采用8个灰度级密度分割的伪彩色变换程序，观察结果，并分别采用2、16、32个灰度级进行伪彩色处理，编程后运行，观察并比较结果。 所用到的函数说明如下： G2C=grayslice(I,m); %该函数用多重（即m－1个）等间隔阈值将灰度图像I转换为索引图像，即m色图像。 % 伪彩色变换——密度分割 I=imread(Miss.bmp); G2C=grayslice(I,8); figure; imshow(I); figure; subplot(222);imshow(G2C,cool(8)); %采用冷色调显示 figure; subplot(223);imshow(G2C,hot(8)); %采用暖色调显示 (三) 根据下图（b），采用灰度级彩色变换的方法，编写MATLAB程序，实现Miss.bmp图像的灰度级彩色变换。 注： 1