[发明专利]一种自适应门限的NCC配准方法

权利要求书

1.一种自适应门限的NCC配准方法，其特征是：

(1)选取一个初始的NCC门限值对两幅图像中的特征点进行配准；

(2)利用配准的特征点对求得两幅图像间的仿射变换矩阵；

(3)将待配准图像利用求得的仿射变换矩阵做仿射变换；

(4)求出仿射变换后的待配准图像与基准图像间的互信息熵；

(5)利用加入了进退法与黄金分割法的powell寻优迭代算法，将互信息作为搜索准则进行NCC门限的自适应迭代；即选取一个适当步长，改变初始门限再次进行(1)—(4)步得到新的互信息熵值；

(6)利用powell算法不断改变NCC门限值，进行迭代；

(7)迭代结束后，找到最大的互信息熵并且得到此时对应的NCC门限值；

(8)将得到的NCC门限值作为最终的门限值进行特征点对的配准。

2.根据权利要求1所述的自适应门限的NCC配准方法，其特征是所述进退法的主要步骤为：

1)给定初始点x0和步长h,h>0；

2)计算f(x0)和f(x0+h)；

3)比较f(x0)和f(x0+h)，有如下三种情况：

●若f(x0)=f(x0+h),则[x0,x0+h]为搜索区间；

●若f(x0)<f(x0+h),则说明x即极大点在x0的左边，置x0=x0+h,h=-h，转到第4）步；

●若f(x0)＞f(x0+h),则说明x即极大点在x0的右边，转到第4）步；

4)计算f[x0+(2^k-1)*h]，k=1,2,…，直到某一个m，使得

f[x0+(2^m-1-1)h]<=f[x0+(2^m-1)h]>=f[x0+(2^m+1-1)h]成立；

令x1=x0+(2^m+1-1)h，x2=x0+(2^m-1-1)h；

5)令a=min{x1，x2}，b=max{x1,x2}，则区间[a,b]为单谷函数f(x)的极小点的一个搜索区间。

3.根据权利要求1或2所述的自适应门限的NCC配准方法，其特征是所述黄金分割法的主要步骤为：

1）给定初始点x0、步长h和收敛精度e；

2）用进退法确定函数f(x)的极小点的一个搜索区间[a,b]；

3）计算：x1=b-k(b-a)和x2=a+k(b-a)，其中k=0.618；

4）比较f(x1)和f(x2)；

5）若f(x1)＞=f(x2)，置b=x2，x2=x1，求新搜索区间的插入点x1=b-k(b-a)，转到步骤6）；若f(x1)<f(x2)，置a=x1,x1=x2，求新搜索区间的插入点x2=a+k(b-a)；

6）若b-a<e，则置x=(a+b)/2，转到步骤7）；否则，转到步骤4）；

7）输出求得的x和f(x)，停止程序执行，结束。