[发明专利]面向非线性辐射畸变的稳健遥感影像特征点对提取方法

权利要求书

1.一种面向非线性辐射畸变的稳健遥感影像特征点对提取方法，其特征在于，其具体步骤如下：

S1，对获取的两幅多模态遥感影像进行影像预处理；

S2，对主从影像分别提取其精确角点与边缘特征点，获取每个影像的特征点：

S3，生成特征点的特征向量；

S4，对多模态影像特征点的特征向量进行配准；

S5，对所有候选匹配点对，利用RANSAC算法作进一步筛选，最终得到精准、均匀且对非线性辐射畸变鲁棒的特征匹配点对。

2.一种如权利要求1所述的面向非线性辐射畸变的稳健遥感影像特征点对提取方法，其特征在于，所述的步骤S1，其具体步骤为：

S11，判断所获取的多模态遥感影像是否为全色黑白，如果是，则转入S12，否则将多模态遥感影像加权成为黑白影像；

S12，设置多模态遥感影像中的一幅影像为主影像，另一幅影像为从影像；

S13，将主影像和从影像分别进行重采样，使其分辨率大小相同；

S14，分别提取主影像和从影像的相位一致性最大矩图与最小矩图；

设置滤波器尺度数N_s＝4，方向数N_o＝6；

对主影像的每个像素点(x,y)，计算其在所有不同尺度s、不同方向o下的相位一致性测度PC(x,y)，计算公式为：

式中，w_o(x,y)为像素点(x,y)处的基于频响范围的权值函数，T为噪声阈值，ξ为偏移量，运算符表示当且仅当其表达式为正时取其本身，否则取0，像素点(x,y)处的振幅分量A_so(x,y)、相位偏差函数ΔΦ_so(x,y)的计算公式分别为：

其中，与的计算公式为：

式中，[E_so(x,y),o_so(x,y)]为像素点(x,y)处的滤波后的响应分量，其表达式为：

[E_so(x,y),O_so(x,y)]

＝[I(x,y)*L^even(x,y,s,o),I(x,y)*L^odd(x,y,s,o)]，

式中,L^even(x,y,s,o)与L^odd(x,y,s,o)分别为二维log-Gabor滤波函数在空间域内像素点(x,y)处的实部与虚部，*表示卷积运算，I(x,y)为像素点(x,y)处的灰度值，其中二维log-Gabor滤波函数的表达式为：

式中，(ρ,θ)代表对数极坐标，(s,o)为滤波器的尺度与方向，ρ_s和θ_so为中心频率对应参数，σ_ρ和σ_θ分别为ρ和θ的带宽；

取方向o的不同尺度数下的PC(x,y)的最大值，作为方向o的相位一致性测度，记作PC(θ_o)，θ_o表示PC(x,y)取最大值时方向o的角度值；

对主影像和从影像，分别计算其所包含的像素点(x,y)的主轴ψ，最大矩M_ψ，最小矩m_ψ，其具体过程包括，

首先计算中间变量a、b、c，其计算公式为，

根据中间变量，计算像素点(x,y)的主轴ψ、最大矩M_ψ、最小矩m_ψ，其计算公式为，

影像中所有像素点的最大矩M_ψ与最小矩m_ψ的集合，分别构成该影像的相位一致性最大矩图与最小矩图。