角点检测算子的代码描述

最简单的一种角点检测算子，但是由于采用了非最大化抑制，效果和OpenCV相当。

TPointInfo = record

Info: TPoint;

w: extended;

Color: TLabColor;

end;

TSinglePointInfoArray = array of TPointInfo;

procedure CornerDetect(Width, Height: longint);

var

i, j, fi, fj, sum: longint;

begin

PointCount := 0;

for i := 7 to Width - 8 do

for j := 7 to Height - 8 do begin

sum := 0;

for fi := i - 7 to i + 7 do

for fj := j - 7 to j + 7 do

sum := sum + abs(ImageGray[i, j] - ImageGray[fi, fj]);

ImagePoint[i, j] := sum div $100;

end;

{标准角点检测算子部分}

for i := 7 to Width - 8 do

for j := 7 to Height - 8 do begin

sum := ImagePoint[i, j];

if sum > $20 then begin

WBPoint[i, j] := true;

Inc(PointCount);

for fi := i - 7 to i + 7 do begin

for fj := j - 7 to j + 7 do

if ImagePoint[fi, fj] > sum then begin

WBPoint[i, j] := false;

Dec(PointCount);

break;

end;

if not WBPoint[i, j] then break;

end;

end else

WBPoint[i, j] := false;

end;

{用非最大化抑制来抑制假角点}

setlength(CornerPoint, PointCount); fi := 0;

for i := 7 to Width - 8 do

for j := 7 to Height - 8 do

if WBPoint[i, j] then begin

CornerPoint[fi].Info.X := i;

CornerPoint[fi].Info.Y := j;

Inc(fi);

end;

{输出为一个点序列}

end;

输入的ImageGray为图像的灰度描述，WBPoint为Boolean数组，ImagePoint为标准角点检测算子运行后的检测值。