边缘检测（8）Wallis微分算子

来源：互联网发布：mac商店编辑：程序博客网时间：2024/06/05 09:27

Wallis微分算子

模板为L=[ 0, -log, 0;

-log, 4log,-log;

0, -log, 0]

该模板是将Laplacian微分算子和对数算子结合的一种锐化算子。因为人眼对画面信号的处理过程中有一个近似的对数运算环节，解决了Laplacian微分算子对画面比较暗的部分锐化比较弱的缺陷。

%Wallis微分算子clear all;I=imread('lena.bmp');%figure%imshow(I);I=im2double(I);[height width R]=size(I);for i=2:height-1    for j=2:width-1        II(i,j)=log10(I(i,j)+1)-0.25*(log10(I(i-1,j)+1)+log10(I(i+1,j)+1)+log10(I(i,j-1)+1)+log10(I(i,j+1)+1));    endendmin1=min(II);min2=min(min1);for i=2:height-1   for j=2:width-1      II(i,j)=46*II(i,j)-min2+0.4;  endend%figure%imshow(II,[]);for i=1:height-1    for j=1:width-1        if (II(i,j)<-0.035)           II(i,j)=0;        else II(i,j)=1;        end    endendfigure;imshow(II,[]);

%Wallis微分算子clear all;I=imread('lena.bmp');%figure%imshow(I);I=im2double(I);[height width R]=size(I);for i=2:height-1    for j=2:width-1        II(i,j)=log10(I(i,j)+1)-0.25*(log10(I(i-1,j)+1)+log10(I(i+1,j)+1)+log10(I(i,j-1)+1)+log10(I(i,j+1)+1));    endendmin1=min(II);min2=min(min1);for i=2:height-1   for j=2:width-1      II(i,j)=46*II(i,j)-min2+0.4;  endend%figure%imshow(II,[]);for i=1:height-1    for j=1:width-1        if (II(i,j)<-0.035)           II(i,j)=0;        else II(i,j)=1;        end    endendfigure;imshow(II,[]);

0 0