matlab锐化代码

x=imread('02.bmp'); [m,n]=size(x); figure,imshow(x); x=double(x); b=zeros(m,n); c=zeros(m,n); d=zeros(m,n); e=zeros(m,n); f=zeros(m,n); for i=1:m-2 for j=1:n-2 c(i+1,j+1)=abs(8*x(i+1,j+1)-x(i,j)-x(i,j+1)-x(i,j+2)-x(i+1,j)-x(i+1,j+2)-x(i+2,j)-x(i+2,j+1)-x(i+2,j+2)); %采用模板a=[-1 -1 -1;-1 8 -1;-1 -1 -1] b(i+1,j+1)=abs(4*x(i+1,j+1)-x(i,j+1)-x(i+1,j)-x(i+1,j+2)-x(i+2,j+1)); %采用模板a=[0 -1 0;-1 4-1;0-1 0] d(i+1,j+1)= abs(5*x(i+1,j+1)-x(i,j+1)-x(i+1,j)-x(i+1,j+2)-x(i+2,j+1)); %采用模板a=[0 -1 0;-1 5-1;0 -1 0] e(i+1,j+1)=abs(4*x(i+1,j+1)+x(i,j)-2*x(i,j+1)+x(i,j+2)-2*x(i+1,j)-2*x(i+1,j+2)+x(i+2,j)-2*x(i+2,j+1)+x(i+2,j+2)); %采用模板a=[1 -2 -1;-2 4 -2;1 -2 1] f(i+1,j+1)=abs(9*x(i+1,j+1)-x(i,j)-x(i,j+1)-x(i,j+2)-x(i+1,j)-x(i+1,j+2)-x(i+2,j)-x(i+2,j+1)-x(i+2,j+2)); %采用模板a=[-1 -1 -1;-1 9 -1;-1 -1 -1] %这儿处理有两种方式，绝对值和统一加一个数,本题采用的是取绝对值 %if b(i+1,j+1)<0 这部分为无用代码 % b(i+1,j+1)=-b(i+1,j+1); %end %if c(i+1,j+1)<0 % c(i+1,j+1)=-c(i+1,j+1) %end end end figure;imshow(uint8(b)) figure;imshow(uint8(c)) figure;imshow(uint8(d)) figure,imshow(uint8(e)) figure,imshow(uint8(f)) %这就是二阶微分及其变形算法的效果