clc;
clear;
I=imread（‘1.jpg‘）;????
[mn]=size（I）;???%寬m，長n
%?利用差分算子（像素（cr）在上下左右四個方向的灰度差分絕對值）提取初選點。
g=double（I）;
T=30;???????????????%閾值
Initpointcoord=[];???%初選點坐標
for?i=2:m-1
????for?j=2:n-1
????????dg1=abs（g（ij）-g（i-1j））;???%上側
????????dg2=abs（g（ij）-g（i+1j））;???%下側
????????dg3=abs（g（ij）-g（ij-1））;???%左側
????????dg4=abs（g（ij）-g（ij+1））;???%右側
????????dg=[dg1?dg2?dg3?dg4];
????????temp=sort（dg）;
????????if?temp（3）>T
????????????Initpointcoord=[Initpointcoord;i?j];
????????end
????end
end
%??以初選點（cr）為中心的w*w窗口中，計算協方差矩陣N與誤差橢圓的圓度q確定備選點。
w=3;??????????????????????????%窗口寬的值
Initlen=length（Initpointcoord）;???%初選點個數
k=floor（w/2）;
Tq=0.51;???????????????????????%閾值Tq
candidatepoint=[];?????????????%候選點
for?no=1:Initlen
????Initpoint=Initpointcoord（no:）;
????c=Initpoint（1）;
????r=Initpoint（2）;
????h=c-k;
????l=c+k-1;
????p=r-k;
????w=r+k-1;
????gu2=0;gv2=0;guv=0;
????for?e=h:l
????????for?f=p:w
????????????gu2=gu2+（g（e+1f+1）-g（ef））^2;
????????????gv2=gv2+（g（ef+1）-g（e+1f））^2;
????????????guv=guv+（g（e+1f+1）-g（ef））*（g（ef+1）-g（e+1f））;
????????end
????end
????DetN=gu2*gv2-guv^2;
????trN=gu2+gv2;
????q=4*DetN/（trN*trN）;
????if?q>Tq
????????candidatepoint=[candidatepoint;c?r];
????end
end
clear?h?p?q?j?c?r
w2=10;
cadipotlen=length（candidatepoint）;
k2=floor（5/w2）;
featurepoint=[];
h=floor（m/w2）;???%窗口行數
d=floor（n/w2）;???%窗口列數
for?p=1:h
????for?q=1:d
????????tp=（p-1）*w2;
????????tq=（q