ch451.m

来自「灰度图象的直方图均衡化；梯度；直方图规定化；灰度图象的直方图数据；灰度直接变换:」· M 代码 · 共 32 行

f = imread('cameraman.tif');
[m,n] = size(f);
f1=double(f);
px=zeros(m,n);
py=zeros(m,n);
for x=1:m-1
for y=1:n-1
 px(x,y)= f1(x+1,y)-f1(x,y);   %对x的偏导
 py(x,y)= f1(x,y+1)-f1(x,y);   %对y的偏导
end;end
Dx=abs(px);           %绝对值
Dy=abs(py);   
Te=(px.^2+py.^2).^0.5; %欧氏距离表示的梯度  
T4=Dx+Dy;              %城区距离表示的梯度
T8=max(Dx,Dy);         %棋盘距离表示的梯度
Te=uint8(Te);
T4=uint8(T4);
T8=uint8(T8);
figure(1);
subplot(231),imshow(f); %
title('原图');
subplot(232),imshow(uint8(px)); %
title('x方向的梯度图');
subplot(233),imshow(uint8(py)); %
title('y方向的梯度图');
subplot(234),imshow(Te); %
title('梯度图：欧氏距离');
subplot(235),imshow(T4); %
title('梯度图：城区距离');
subplot(236),imshow(T8); %
title('梯度图：棋盘距离');