maketip.m

来自「matlab处理图像的一些基本方法。其中有一部分mex程序需要安装编译」· M 代码 · 共 39 行

function T=maketip(x, y)
%MAKETIP Make the tip matrix from a profile

%
% Copyright (c) 1995 by Claudio Rivetti and Mark Young
% [e-mail],    [e-mail]
%

A=1;
B=length(x);
C=find(y==min(min(y)));

T=zeros(B);


for col=A:C
  for row=A:B
    l=sqrt((C-col)^2 + (row-B/2)^2);
      if l < C-A
        T(row,col)=y(C-l);
      else
        T(row,col)=y(A);
      end
  end
end

for col=C+1:B
  for row=A:B
    l=sqrt((C-col)^2 + (row-B/2)^2);
      if l < B-C
        T(row,col)=y(C+l);
      else
        T(row,col)=y(B);
      end
  end
end

return