数字图像处理实验三

数字图像处理—实验三一．实验内容：⑴图像的锐化：使用Sobel，Laplacian，unsharpmasking，high-boostfiltering算子分别对图像进行运算，观察并体会运算结果。⑵综合练习：对需要进行处理的图像分析，正确运用所学的知识，采用正确的步骤，对图像进行各类处理，以得到令人满意的图像效果。二．实验目的：学会用Matlab中的以下函数对输入图像按实验内容进行运算；感受各种不同的图像处理方法对最终图像效果的影响。imfilter；fspecial；imadjust；imadd；immultiply.三．实验步骤：仔细阅读Matlab帮助文件中有关以上函数的使用说明，能充分理解其使用方法并能运用它们完成实验内容。将Fig3.41(c).jpg图像文件读入Matlab，使用imfilter函数分别采用不同的算子对其作锐化运算，显示运算前后的图像。算子的输入可采用直接输入法。其中Sobel，Laplacian，unsharpmasking也可用fspecial函数产生。各类算子如下： High-boostfilteringUnsharpmasking对角线LaplacianSobel垂直梯度High-boostfilteringUnsharpmasking对角线LaplacianSobel垂直梯度将Fig3.46(a).jpg图像文件读入Matlab，按照以下步骤对其进行处理：用带对角线的Laplacian对其处理，以增强边缘。用imadd函数叠加原始图像。可以看出噪声增强了，应想法降低。获取Sobel图像并用imfilter对其进行5×5邻域平均，以减少噪声。用immultiply函数获取Sobel，Laplacian增强图像，噪声得以减少。最后用imadjust函数做power-law灰度变换，以增强dynamicrange。四．实验报告要求： 用imshow函数显示经每一步处理后的图像，打印后提交。本次实验后要求认真写出三次实验的心得体会及提出改良实验的一些建议。实验报告上的其他内容，按常规实验报告要求办。%2i2=imread('Fig3.41(c).jpg');i2=im2double(i2);subplot(3,2,1);imshow(i2);title('orgimg');h1=[-1-2-1000121];h2=[1111-81111];h3=[1111-91111];A=-1;h4=[1111-A-81111];subplot(3,2,2);g1=imfilter(i2,h1);g1=imadd(i2,g1);imshow(g1);title('g1');subplot(3,2,3);g2=imfilter(i2,h2);g2=imadd(i2,g2);imshow(g2);title('g2');subplot(3,2,4);g3=imfilter(i2,h3);g3=imadd(i2,g3);imshow(g3);title('g3');subplot(3,2,5);g4=imfilter(i2,h4);g4=imadd(i2,g4);imshow(g4);title('g4');%5i5=imread('Fig3.46(a).jpg');i5=im2double(i5);subplot(3,2,6);imshow(i5);title('5orgimg');%51figure;subplot(3,2,1);g51=imfilter(i5,h2);imshow(g51);title('51');%52subplot(3,2,2);i52=imadd(i5,g51);imshow(i52);title('52');subplot(3,2,3);h52=fspecial('average',[33]);%生成一3×3邻域平均窗函数g52=filter2(h52,i51);%求邻域平均imshow(g52);title('g52');%53subplot(3,2,4);H=fspecial('sobel');SOBEL=imfilter(i5,H);i53=imadd(i5,SOBEL);imshow(i53);title('53');subplot(3,2,5);H=fspecial('average',[55]);%生成一5×5邻域平均窗函数g53=filter2(H,i53);%求邻域平均imshow(g53);title('g53');%54subplot(3,2,6);i54l=immultiply(i52,i52);imshow(i54);title('54l');figure;subplot(3,2,1);i54s=immultiply(i53,i53);imshow(i54);title('54s');%55subplot(3,2,2);i55l=imadjust(i52,[],[],0.2);imshow(i55l);title('55l');subplot(3,2,3);i55=imadjust(i54l,[],[],0.1);imshow(i55);title('g55l');subplot(3,2,4);i55s=imadjust(i53,[],[],0.1);imshow(