frierdscripwhbrightness.m

来自「一个简单的傅立叶带亮度描绘子提取图像的形状特征程序」· M 代码 · 共 32 行

function ft=frierdscripwhbrightness(im)

ed=edge(im,'canny');
ethin=bwmorph(ed,'thin',Inf);
[m,n]=size(im);

ft=zeros(1,9);
%boundary pixels
[x,y]=find(ethin==1);
 
%centroid distance
N=length(x);
val=zeros(N,1);

xc=sum(x)/N;
yc=sum(y)/N;

rd=sqrt((x-xc).^2+(y-yc).^2);
for i=1:N
    val(i)=im(x(i),y(i));
end

t=fft(rd.*val);
abst=abs(t);
%nt=abst/abst(1);
ft=abst(1:8);

u11=(x-xc)'*(y-yc);
u20=(x-xc)'*(x-xc);
u02=(y-yc)'*(y-yc);
e=(u20+u02-sqrt((u20-u02)^2+4*u11^2))/(u20+u02+sqrt((u20-u02)^2+4*u11^2));
ft(9)=e;