gpzq.m

来自「此程序为本人编写的关于数字图像处理方面的matlab程序」· M 代码 · 共 54 行

clc
clear
close all
a0=imread('D:\Mary.bmp');
subplot(221);
imshow(a0);
title('原图像');
a=a0;
a=double(a);
%给图像矩阵数据的最上和最下加入两行0
b=zeros(1,256);
c=[b;a;b];
%给图像矩阵数据的最左和最右加入两列0
d=zeros(258,1);
e=[d c d];
%建立一个3×3的平滑滤波器(模板),即3行3列的1的矩阵
panel=[-1 -1 -1;-1 8 -1;-1 -1 -1];
%进行模板运算
for i=2:257
    for j=2:257
        %取出图像中的3×3的矩阵
        image33=[e(i-1,j-1) e(i-1,j) e(i-1,j+1);
                 e(i,j-1)   e(i,j)   e(i,j+1);
                 e(i+1,j-1) e(i+1,j) e(i+1,j+1)];
             temp=0;
             for m=1:3
                 for n=1:3  
                  temp= image33(m,n)*panel(m,n)+temp;
                 end
             end
             panelimage(i,j)=temp/9;
             
    end
end
for  i=2:257
    for j=2:257
        c(i,j)=9*a(i-1,j-1)+panelimage(i,j);
        d(i,j)=.1*a(i-1,j-1)+.1*panelimage(i,j);
    end
end
%d=a+.1*panelimage;
panelimage=mat2gray(panelimage);
subplot(222);
imshow(panelimage);
title('锐化滤波图像');
c=mat2gray(c);
subplot(223);
imshow(c);
title('大系数高平增强图像');
d=mat2gray(d);
subplot(224);
imshow(d);
title('小系数高平增强图像');