%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%?
%???Harris角點(diǎn)提取算法和改進(jìn)的Harris角點(diǎn)提取算法??????????????????????????????%?
%???2005-11-25????????????????????????????????????????????????????????????????%?
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%?
clear;?
FileInfo?=?imfinfo（‘E:\Image1\MagImg1.bmp‘）;?
[ImageDatamap]?=?imread（‘E:\Image1\r4.bmp‘）;?
%??
%?if（strcmp（‘truecolor‘FileInfo.ColorType）?==?1）?
%????ImageData?=?im2uint8（rgb2gray（ImageData））;?
%?elseif（strcmp（‘indexed‘FileInfo.ColorType）?==?1）?
%????ImageData?=?im2uint8（ind2gray（ImageDatamap））;???
%?end?????
ori_im=rgb2gray（ImageData）;?
%%ori_im=histeq（ori_im）;
%ori_im?=?imread（‘a(chǎn)rroyo-r.tiff‘）;?????%?讀取圖像?
?
%?fx?=?[5?0?-5;8?0?-8;5?0?-5];??????????%?高斯函數(shù)一階微分，x方向（用于改進(jìn)的Harris角點(diǎn)提取算法）?
fx?=?[-2?-1?0?1?2];?????????????????%?x方向梯度算子（用于Harris角點(diǎn)提取算法）?
Ix?=?filter2（fxori_im）;??????????????%?x方向?yàn)V波?
%?fy?=?[5?8?5;0?0?0;-5?-8?-5];??????????%?高斯函數(shù)一階微分，y方向（用于改進(jìn)的Harris角點(diǎn)提取算法）?
fy?=?[-2;-1;0;1;2];?????????????????%?y方向梯度算子（用于Harris角點(diǎn)提取算法）?
Iy?=?filter2（fyori_im）;??????????????%?y方向?yàn)V波?
Ix2?=?Ix.^2;?
Iy2?=?Iy.^2;?
Ixy?=?Ix.*Iy;?
clear?Ix;?
clear?Iy;?
?
h=?fspecial（‘gaussian‘[7?7]2）;??????%?產(chǎn)生7*7的高斯窗函數(shù)，sigma=2?
?
Ix2?=?filter2（hIx2）;?
Iy2?=?filter2（hIy2）;?
Ixy?=?filter2（hIxy）;?
?
height?=?size（ori_im1）;?
width?=?size（ori_im2）;?
result?=?zeros（heightwidth）;?????????%?紀(jì)錄角點(diǎn)位置，角點(diǎn)處值為1?
?
R?=?zeros（heightwidth）;?
?
Rmax?=?0;??????????????????????????????%?圖像中最大的R值?
for?i?=?1:height?
????for?j?=?1:width?
????????M?=?[Ix2（ij）?Ixy（ij）;Ixy（ij）?Iy2（ij）];?????????????%?auto?correlation?matrix?
????????R（ij）?=?det（M）-0.06*（trace（M））^2;?????????????????????%?計(jì)算R?
????????if?R（ij）?>?Rmax?
????????????Rmax?=?R（ij）;?
????????end;?
????end;?
end;?
?
cnt?=?0;?
for?i?=?2:height-1?
????for?j?=?2:width-1?
????????%?進(jìn)行非極大抑制，窗口大小3*3?
????????if?R（ij）?>?0.01*Rmax?&&?R（ij）?>?R（i-1j-1）?&&?R（ij）?>?R（i-1j）?&&?R（ij）?>?R（i-1j+1）?&&?R（ij）?>?R（ij-1）?&&?R（ij）?>?R（ij+1）?&&?R（ij）?>?R（i+1j-1）?&&?R（ij）?>?R（i+1j）?&&?R（ij）?>?R（i+1j+1）?
????????????result（ij）?=?1;?
????????????cnt?=?cnt+1;?
????????end;?
????end;?
end;?
?
[posc?posr]?=?find（result?==?1）;?
cnt??????????????????????????????????????%?角點(diǎn)個(gè)數(shù)?
imshow（ori_im）;?
hold?on;?
plot（posrposc‘r+‘）;?

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件?????631648??2010-06-30?16:11??角點(diǎn)檢測\2_2_CSS角點(diǎn)檢測.pdf

?????文件???????2139??2000-01-22?12:00??角點(diǎn)檢測\Corner\a.gif

?????文件???????1707??2000-01-22?12:00??角點(diǎn)檢測\Corner\b.gif

?????文件???????1745??2000-01-22?12:00??角點(diǎn)檢測\Corner\c.gif

?????文件??????10900??2010-07-09?10:29??角點(diǎn)檢測\Corner\corner.jpg

?????文件??????13922??2006-01-24?14:24??角點(diǎn)檢測\Corner\corner.m

?????文件????????338??2010-07-09?10:29??角點(diǎn)檢測\Corner\cout.mat

?????文件???????1426??2000-01-22?12:00??角點(diǎn)檢測\Corner\d.gif

?????文件???????1759??2000-01-22?12:00??角點(diǎn)檢測\Corner\e.gif

?????文件???????9603??2010-07-09?10:29??角點(diǎn)檢測\Corner\edge.jpg

?????文件???????1879??2000-01-22?12:00??角點(diǎn)檢測\Corner\f.gif

?????文件???????1916??2000-01-22?12:00??角點(diǎn)檢測\Corner\g.gif

?????文件??????19024??2005-05-10?16:46??角點(diǎn)檢測\Corner\gui_corner.fig

?????文件???????8565??2005-05-10?16:54??角點(diǎn)檢測\Corner\gui_corner.m

?????文件???????2182??2006-01-17?17:45??角點(diǎn)檢測\Corner\h.gif

?????文件??????66568??2003-10-17?05:43??角點(diǎn)檢測\Corner\house.bmp

?????文件??????21406??2010-07-07?16:21??角點(diǎn)檢測\FART關(guān)鍵點(diǎn)C和MATLab程序\fast-C-src-1.2.zip

?????文件?????174649??2010-07-07?16:21??角點(diǎn)檢測\FART關(guān)鍵點(diǎn)C和MATLab程序\fast-matlab-src.zip

?????文件??????23782??2004-03-07?17:05??角點(diǎn)檢測\Forstner?角點(diǎn)檢測?自編\Forstner?角點(diǎn)檢測?自編\dibapi.cpp

?????文件???????1876??2000-02-03?07:29??角點(diǎn)檢測\Forstner?角點(diǎn)檢測?自編\Forstner?角點(diǎn)檢測?自編\dibapi.h

?????文件??????33304??2010-07-08?11:38??角點(diǎn)檢測\Forstner?角點(diǎn)檢測?自編\Forstner?角點(diǎn)檢測?自編\imagetest.aps

?????文件???????4096??2010-07-08?15:18??角點(diǎn)檢測\Forstner?角點(diǎn)檢測?自編\Forstner?角點(diǎn)檢測?自編\imagetest.clw

?????文件???????4230??2005-09-22?21:51??角點(diǎn)檢測\Forstner?角點(diǎn)檢測?自編\Forstner?角點(diǎn)檢測?自編\imagetest.cpp

?????文件???????4999??2007-09-11?21:15??角點(diǎn)檢測\Forstner?角點(diǎn)檢測?自編\Forstner?角點(diǎn)檢測?自編\imagetest.dsp

?????文件????????541??2005-09-22?21:51??角點(diǎn)檢測\Forstner?角點(diǎn)檢測?自編\Forstner?角點(diǎn)檢測?自編\imagetest.dsw

?????文件???????1389??2005-09-22?21:51??角點(diǎn)檢測\Forstner?角點(diǎn)檢測?自編\Forstner?角點(diǎn)檢測?自編\imagetest.h

?????文件?????181248??2010-07-08?15:18??角點(diǎn)檢測\Forstner?角點(diǎn)檢測?自編\Forstner?角點(diǎn)檢測?自編\imagetest.ncb

?????文件?????386048??2010-07-08?15:18??角點(diǎn)檢測\Forstner?角點(diǎn)檢測?自編\Forstner?角點(diǎn)檢測?自編\imagetest.opt

?????文件???????1587??2007-09-27?15:09??角點(diǎn)檢測\Forstner?角點(diǎn)檢測?自編\Forstner?角點(diǎn)檢測?自編\imagetest.plg

?????文件??????14718??2007-09-11?14:54??角點(diǎn)檢測\Forstner?角點(diǎn)檢測?自編\Forstner?角點(diǎn)檢測?自編\imagetest.rc

............此處省略91個(gè)文件信息