CN106548462A - 基于薄板样条插值的非线性sar图像几何校正方法 - Google Patents

基于薄板样条插值的非线性sar图像几何校正方法 Download PDF

Info

Publication number
CN106548462A
CN106548462A CN201610948742.6A CN201610948742A CN106548462A CN 106548462 A CN106548462 A CN 106548462A CN 201610948742 A CN201610948742 A CN 201610948742A CN 106548462 A CN106548462 A CN 106548462A
Authority
CN
China
Prior art keywords
point
transformation matrix
affine transformation
match
affine
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201610948742.6A
Other languages
English (en)
Other versions
CN106548462B (zh
Inventor
李亚超
项宇泽
邹倩雅
全英汇
邓欢
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xidian University
Original Assignee
Xidian University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xidian University filed Critical Xidian University
Priority to CN201610948742.6A priority Critical patent/CN106548462B/zh
Publication of CN106548462A publication Critical patent/CN106548462A/zh
Application granted granted Critical
Publication of CN106548462B publication Critical patent/CN106548462B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/80Geometric correction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/14Transformations for image registration, e.g. adjusting or mapping for alignment of images
    • G06T3/147Transformations for image registration, e.g. adjusting or mapping for alignment of images using affine transformations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10032Satellite or aerial image; Remote sensing
    • G06T2207/10044Radar image

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Image Analysis (AREA)

Abstract

本发明公开了一种基于薄板样条插值的非线性SAR图像几何校正方法,主要解决现有技术无法精确校正图像非线性失真的问题。其实现步骤是:1.利用SIFT‑Like算法分别提取基准图和实时图的特征点;2.利用相似性度量法在基准图和实时图的特征点中选取匹配点;3.利用RANSAC算法对匹配点进行提纯并求解仿射变换矩阵,得到精确的仿射变换模型,即图像的线性映射关系;4.利用薄板样条插值函数在仿射变换模型上对仿射变换矩阵与实际情况的误差进行局部非线性拟合,得到图像的非线性映射关系,实现对实时图的精确几何校正。本发明能消除非线性几何畸变的影响,提高了校正精度,可用于对雷达接收图像信号非线性失真的精确几何校正。

Description

基于薄板样条插值的非线性SAR图像几何校正方法
技术领域
本发明属于雷达技术领域,主要涉及一种非线性SAR图像校正方法,可用于对雷达接收图像信号的精确几何校正。
背景技术
作为SAR系统成像的重要评估手段,针对SAR图像几何校正的研究是成像算法设计以及基于图像的目标识别的理论依据。合成孔径雷达SAR在成像过程中,由于地形起伏和平台运动状态的变化等原因,使得图像存在一定的非线性几何失真,因此所获取图像在几何位置上发生畸变。
在SAR图像中存在一定程度的几何形变,这中间包括线性映射和非线性映射。目前传统方法是将图像间的非线性映射变换关系近似为一种线性映射关系。一般采用RANSAC算法求取其仿射变换矩阵,该算法通过对获取的匹配点进行提纯,删除错误匹配点,利用最终正确匹配点集估计图像之间的线性变换关系,如旋转变换、偏移变换以及伸缩变换等,最终求得基准图与待匹配图之间的仿射变换矩阵。基于RANSANC算法求解的基准图与实时图之间的映射关系是线性变换关系,匹配线性图像时,空间变换通常选为仿射变换,该变换为全局的变换。然而由于实时图中存在一定的非线性几何失真,使得各匹配点的对应关系会有一定的非线性误差,仿射变换并不能反映局部的和全部的非线性差异,因此会带来匹配的误差,导致图像畸变。
发明内容
本发明的目的在于针对上述现有技术的不足,提出一种基于薄板样条插值的非线性SAR图像几何校正方法,获得更高的校正精度。
为实现上述目的,本发明的技术方案包括如下:
(1)将通过星载数据获得的无失真图像作为基准图,利用SIFT-Like算法分别提取基准图与机载的实时图特征点;
(2)根据相似性度量法在基准图特征点和机载实时图特征点中选取出两者的匹配点作为仿射变换模型的控制点,再利用RANSAC算法对控制点进行提纯,求解基准图和实时图这两幅图像各点之间的仿射变换关系,即仿射变换矩阵,得到精确的仿射变换模型;
(3)利用薄板样条插值函数在仿射变换模型上对仿射矩阵与实际情况的误差进行局部非线性拟合,实现对实时图的精确几何校正。
本发明与现有技术相比较,有如下优点:
1.获得更高的校正精度
经典的模型参数估计方法是最小二乘法,但是它没有检测和排除错误数据的方法,用最小二乘法拟合得到的参数准确性不高,甚至出现错误。本发明利用对匹配集采样和对匹配点验证的策略,在包含异常数据的匹配点集中,通过逐步的迭代和筛选,可剔除错误匹配点,最终根据有效样本数据拟合仿射变换矩阵,消除错误数据对仿射变换矩阵的干扰,提高了拟合精确度,获得更高的校正精度。
2.提高SAR图像的质量
传统方法是将图像间的非线性映射变换关系近似为一种线性映射关系,忽略了图像的非线性几何失真,并不能全面反映实时图与基准图之间的映射关系。本发明在仿射变换模型的基础上,运用薄板样条插值函数对仿射矩阵与实际情况的误差进行局部非线性拟合,得到实时图和基准图之间的非线性映射关系,对实时图像进行非线性几何校正,减少图像畸变,提高了SAR图像的质量。
附图说明
图1是本发明的实现流程图;
图2是原始数据三维坐标图;
图3是用本发明对原始数据随机抽样后的三维坐标图;
图4是现有薄板样条插值函数的三维坐标图;
图5是用本发明对随机抽样后数据进行薄板样条插值后的数据三维坐标图;
图6是对本发明插值后数据与原始数据与作差后的失真校正误差三维坐标图;
具体实施方式
以下参照附图对本发明作进一步作详细描述:
参照图1,本发明的实现步骤如下:
步骤1,SIFT-Like算法提取特征点。
图像的特征主要有三种,即点、线、面,其中点特征是最为常见的特征之一,主要有房屋的角点,图像的边缘点,河流、道路等的交汇点等,点特征提取方法有Moravec算法、Forstner算法、SIFT算法及其改进算法等。本发明利用SIFT算法的一种改进算法SIFT-Like分别提取基准图与机载的实时图特征点,其步骤如下:
(1a)将通过星载数据获得的无失真图像作为基准图,在基准图和实时图上分别生成SAR-Harris尺度空间,采用指数加权均值比梯度法得到SAR-Harris函数;
(1b)在SAR-Harris尺度空间,计算每个点的角点响应值,利用其SAR-Harris函数得到初步角点,对初步角点进行非极大值抑制,检测出二维坐标空间中的极值点;
(1c)计算每个坐标空间中极值点的拉普拉斯响应值,并与其相邻两个尺度空间中对应位置的拉普拉斯响应值比较,若该点为拉普拉斯极值点,则该点为尺度空间极值点,即图像的特征点。
步骤2,利用相似性度量法提取匹配点。
利用相似性的度量方法有很多,在图像匹配方面可以通过选取不同的特征描述子来提取特征点,特征描述子是对图像特征的定量化描述,它以向量形式刻画了特征点邻域局部区域的形状特征以及纹理结构特征。理想的特征描述子应该具备高的稳健性、独特性和匹配速度。特征描述子主要分为三类:基于滤波的描述子、基于分布的描述子和基于微分的描述子。常见的描述子包括:GLOH描述子、不变矩描述子、形状描述子、Fourier描述子、线特征描述子、SIFT描述子、形状上下文描述子。
本实例运用但不限于SIFT描述子对每个提取出的特征点进行特征描述,即先将每个点用一个128×1维的向量表示;再采用遍历搜索算法,分别计算实时图中每一个点与基准图中各点间的特征向量欧式距离,当特征点对的最近邻与次近邻的欧式距离比值小于0.8时,则认为该特征点对为匹配点对。
步骤3,提纯匹配点并获得线性映射关系。
按照不同检测方法使用的技术原理,常用的误匹配点检测技术可以分为三类:基于函数拟合的方法、基于统计模型的方法、基于图的方法。基于函数拟合的方法是假设所有正确匹配点满足某一个函数模型,再使用最小二乘拟合方法求解函数模型的系数,然后通过计算匹配点和函数模型的误差来检测误匹配点;基于统计模型的方法是统计满足给定条件的最多数量的正确匹配点,其典型代表是RANSAC方法;基于图的方法是使用原始图像和基准图像上匹配点的分布和邻域关系来检测误匹配点。
本实例运用但不限于RANSAC算法检测出错误匹配点,提纯出正确的匹配点集,求解基准图和实时图之间的仿射变换矩阵获得基准图和实时图之间的线性映射关系,其步骤如下:
(2a)对一个匹配点集P,假设初始化模型参数需要的样本数为n,从匹配点集P中随机抽取n对匹配点组成匹配点集P的一个子集S,并根据抽取的这n对匹配点初始化仿射变换矩阵H;
(2b)设置一个误差容忍度,判断子集S的补集中匹配点是否在该误差容忍度范围内符合仿射变换矩阵H,将补集中的符合仿射变换矩阵H匹配点加入到子集S中,补集中符合仿射变换矩阵H的匹配点称为内点,该子集S称为仿射变换矩阵H的一致集,一致集S中的匹配点个数为t;
(2c)将匹配点集P中匹配点数目的90%设置为特征点数目阈值T,根据当前一致集S中的匹配点个数t判断是否大于特征点数目阈值T:
若t≥T,则执行步骤(2d);
若t<T,则执行步骤(2e)。
(2d)设置一个仿射变换矩阵H的最小错误概率Pmin,判断当前仿射变换矩阵H的错误概率PH是否小于最小错误概率Pmin
若PH≤Pmin,根据一致集S中的所有匹配点,采用最小二乘法对当前仿射变换矩阵H进行重新拟合,得到最优的仿射变换矩阵H′;
若PH>Pmin,则执行步骤(2e);
(2e)通过迭代来求取最优的仿射变换矩阵H′:
(2e1)计算最大的迭代次数N:
其中ω为匹配点集P中任意一点为内点的概率,p为在N次迭代中至少有一次抽样的n个样本均为内点的概率;
(2e2)重复步骤(2a)~(2d),不断寻找最优的仿射变换矩阵H′,记重复次数为m:
若m<N时循环结束,则找到了最优的仿射变换矩阵H′;
若m=N时还未找到最优的仿射变换矩阵H′,则选取当前仿射变换矩阵H作为最优的仿射变换矩阵H′。
步骤4,获得非线性映射关系。
径向基函数是一种兼顾全局和局部单元的变换,可以将两幅图像上的匹配点一一对应,图像上其它点的位置通过匹配点插值得到,用于图像的非线性变换插值与统计形状分析,反映图像间的局部非线性差异。常用的径向基函数有TPS函数、高斯插值函数、multiquadric函数、inverse quadric函数、inverse multi-quadric函数、reflectedsigmoid函数。
本实例运用但不限于TPS插值函数对仿射变换矩阵与实际情况的误差进行局部非线性拟合,获得基准图和实时图之间的非线性映射关系,其步骤如下:
(3a)建立基于薄板样条函数TPS图像之间映射关系数学模型:
其中,X′,Y′分别表示实时图上匹配点(X′,Y′)的横坐标和纵坐标,X,Y分别表示基准图上与(X′,Y′)对应的匹配点(X,Y)的横坐标和纵坐标,f1表示实时图匹配点横坐标X′与基准图匹配点(X,Y)之间的映射关系,f2表示实时图匹配点纵坐标Y′与基准图匹配点(X,Y)之间的映射关系;A1+A2X+A3Y,C1+C2X+C3Y分别表示匹配点对(X′,Y′)与(X,Y)之间在横坐标方向和纵坐标方向上的线性映射关系,A1表示横坐标方向上线性映射常系数,A2表示横坐标方向上与X有关的线性映射系数,A3表示横坐标方向上与Y有关的线性映射系数,C1表示纵坐标方向上线性映射常系数,C2表示纵坐标方向上与X有关的线性映射系数,C3表示纵坐标方向上与Y有关的线性映射系数; 分别表示匹配点对(X′,Y′)与(X,Y)之间在横坐标方向和纵坐标方向上的非线性映射关系,为核函数,rij表示(Xj,Yj),(Xi,Yi)两点间的距离;Bi为横坐标方向上的非线性映射系数,Di为纵坐标方向上的非线性映射系数;
(3b)求解(3a)中基于薄板样条函数TPS图像之间映射关系数学模型,得到横坐标方向非线性映射系数Bi和纵坐标方向上非线性映射系数Di,其步骤如下:
(3b1)将最优的仿射变换矩阵H'与基准图匹配点坐标(X,Y)相乘得到横坐标方向上的三个线性映射系数A1,A2,A3和纵坐标方向上的三个线性映射系数C1,C2,C3
(3b2)生成核函数的矩阵:
其中rij=||(Xj,Yj)-(Xi,Yi)||2+d2,rij表示实时图上第i个点到第j个点的距离,d2为调节插值平面的弯曲程度的参数,通过实验确定,当d2增加时,插值平面更加均匀;
(3b3)将横坐标方向上的三个线性映射系数A1,A2,A3的矩阵带入等式求解得横坐标方向上的非线性映射系数Bi
(3b4)将纵坐标方向上的三个线性映射系数C1,C2,C3的矩阵带入等式求解得纵坐标方向上的非线性映射系数Di
本发明的效果可以通过仿真实验进一步说明:
一.仿真条件
SAR实时图像的非线性失真误差值会随图像中各点到场景中心点距离的增大而增大,而二次函数刚好可以很好的表现出这一特性。本仿真实验用一个二次函数作为原始数据,模拟出SAR实时图像上各点的失真误差随与中心点距离的变化关系。以场景中心点为原点建立三维坐标系,X轴、Y轴分别表示实际图像中各点的横、纵坐标,Z轴表示SAR实时图像与基准图像的同名点间的失真误差大小。
仿真实验采用图2所示的原始数据三维坐标图。
二.仿真内容
仿真1,根据实际工程中基于特征匹配的算法思想,用本发明对图2所示的原始数据随机抽样,用随机抽样获得的点模拟图像中的特征点,随机抽样后的数据点三维坐标图如图3所示。
仿真2,运用如图4所示现有的薄板样条函数对仿真1得到的数据点进行薄板样条插值拟合,拟合后各数据点的三维坐标图如图5所示。
仿真3,将仿真2得到的数据与原始数据作差得到各点的失真校正误差,失真校正误差三维坐标图如图6所示。
图6表明,经薄板样条插值拟合后的失真校正误差与原始数据相比较明显减小,越靠近场景中心点薄板样条函数插值拟合的效果越好。
仿真4,提取仿真3得到的数据分别在X轴,Y轴方向上与场景中心点距离为-50m,-30m,-20m,0m,20m,30m,50m的数据点失真校正误差,如表1所示。
表1薄板样条插值拟合后数据点失真校正误差
从表1可见,距离场景中心点越近的特征点失真误差值越小,同时整体的拟合误差均不大于1m,符合实际应用过程中精确几何校正需求。

Claims (5)

1.基于薄板样条插值的非线性SAR图像几何校正方法,包括有如下步骤:
(1)将通过星载数据获得的无失真图像作为基准图,利用SIFT-Like算法分别提取基准图与机载的实时图特征点;
(2)根据相似性度量法在基准图特征点和机载实时图特征点中选取出两者的匹配点作为仿射变换模型的控制点,再利用RANSAC算法对控制点进行提纯,求解基准图和实时图这两幅图像各点之间的仿射变换关系,即仿射变换矩阵,得到精确的仿射变换模型;
(3)利用薄板样条插值函数在仿射变换模型上对仿射矩阵与实际情况的误差进行局部非线性拟合,实现对实时图的精确几何校正。
2.根据权利要求1中所述方法,其中步骤(1)中利用SIFT-Like算法分别提取基准图与机载的实时图特征点,按如下步骤进行:
(1a)在基准图和实时图上分别生成SAR-Harris尺度空间,采用指数加权均值比梯度法得到SAR-Harris函数;
(1b)在SAR-Harris尺度空间,计算每个点的角点响应值,利用其SAR-Harris函数得到初步角点,对初步角点进行非极大值抑制,检测出二维坐标空间中的极值点;
(1c)计算每个坐标空间中极值点的拉普拉斯响应值,并与其相邻两个尺度空间中对应位置的拉普拉斯响应值比较,若该点为拉普拉斯极值点,则该点为尺度空间极值点,即图像的特征点。
3.根据权利要求1中所述方法,其中步骤(2)中利用RANSAC算法对控制点进行提纯,求解基准图和实时图这两幅图像各点之间的仿射变换关系,按如下步骤进行:
(2a)对一个匹配点集P,假设初始化模型参数需要的样本数为n,从匹配点集P中随机抽取n对匹配点组成匹配点集P的一个子集S,并根据抽取的这n对匹配点初始化仿射变换矩阵H;
(2b)设置一个误差容忍度,判断子集S的补集中匹配点是否在该误差容忍度范围内符合仿射变换矩阵H,将补集中的符合仿射变换矩阵H匹配点加入到子集S中,补集中符合仿射变换矩阵H的匹配点称为内点,该子集S称为仿射变换矩阵H的一致集,一致集S中的匹配点个数为t;
(2c)将匹配点集P中匹配点数目的90%设置为特征点数目阈值T,根据当前一致集S中的匹配点个数t判断是否大于特征点数目阈值T:
若t≥T,则执行步骤(2d);
若t<T,则执行步骤(2e)。
(2d)设置一个仿射变换矩阵H的最小错误概率Pmin,判断当前仿射变换矩阵H的错误概率PH是否小于最小错误概率Pmin
若PH≤Pmin,根据一致集S中的所有匹配点,采用最小二乘法对当前仿射变换矩阵H进行重新拟合,得到最优的仿射变换矩阵H′;
若PH>Pmin,则执行步骤(2e);
(2e)通过迭代来求取最优的仿射变换矩阵H′:
(2e1)计算最大的迭代次数N:
N = l o g ( 1 - p ) log ( 1 - &omega; n )
其中ω为匹配点集P中任意一点为内点的概率,p为在N次迭代中至少有一次抽样的n个样本均为内点的概率;
(2e2)重复步骤(2a)~(2d),不断寻找最优的仿射变换矩阵H′,记重复次数为m:
若m<N时循环结束,则找到了最优的仿射变换矩阵H′;
若m=N时还未找到最优的仿射变换矩阵H′,则选取当前仿射变换矩阵H作为最优的仿射变换矩阵H′。
4.根据权利要求1中所述方法,其中步骤(3)中利用薄板样条插值函数在仿射变换模型上对仿射矩阵与实际情况的误差进行局部非线性拟合,按如下步骤进行:
(3a)建立基于薄板样条函数TPS图像之间映射关系数学模型:
其中,X′,Y′分别表示实时图上匹配点(X′,Y′)的横坐标和纵坐标,X,Y分别表示基准图上与(X′,Y′)对应的匹配点(X,Y)的横坐标和纵坐标,f1表示实时图匹配点横坐标X′与基准图匹配点(X,Y)之间的映射关系,f2表示实时图匹配点纵坐标Y′与基准图匹配点(X,Y)之间的映射关系;A1+A2X+A3Y,C1+C2X+C3Y分别表示匹配点对(X′,Y′)与(X,Y)之间在横坐标方向和纵坐标方向上的线性映射关系,A1表示横坐标方向上线性映射常系数,A2表示横坐标方向上与X有关的线性映射系数,A3表示横坐标方向上与Y有关的线性映射系数,C1表示纵坐标方向上线性映射常系数,C2表示纵坐标方向上与X有关的线性映射系数,C3表示纵坐标方向上与Y有关的线性映射系数;
分别表示匹配点对(X′,Y′)与(X,Y)之间在横坐标方向和纵坐标方向上的非线性映射关系,为核函数,rij表示(Xj,Yj),(Xi,Yi)两点间的距离;Bi为横坐标方向上的非线性映射系数,Di为纵坐标方向上的非线性映射系数;
(3b)求解(3a)中基于薄板样条函数TPS图像之间映射关系数学模型,得到横坐标方向非线性映射系数Bi和纵坐标方向上非线性映射系数Di
5.根据权利要求4或3所述的方法,其中步骤(3b)中求解(3a)中基于薄板样条函数TPS图像之间映射关系数学模型的非线性映射系数Bi,Di,按如下步骤进行:
(3b1)将最优的仿射变换矩阵H'与基准图匹配点坐标(X,Y)相乘得到横坐标方向上的三个线性映射系数A1,A2,A3和纵坐标方向上的三个线性映射系数C1,C2,C3
(3b2)生成核函数的矩阵:
其中rij=||(Xj,Yj)-(Xi,Yi)||2+d2,rij表示实时图上第i个点到第j个点的距离,d2为调节插值平面的弯曲程度的参数,通过实验确定,当d2增加时,插值平面更加均匀;
(3b3)将横坐标方向上的三个线性映射系数A1,A2,A3的矩阵带入等式求解得横坐标方向上的非线性映射系数Bi
(3b4)将纵坐标方向上的三个线性映射系数C1,C2,C3的矩阵带入等式求解得纵坐标方向上的非线性映射系数Di
CN201610948742.6A 2016-11-02 2016-11-02 基于薄板样条插值的非线性sar图像几何校正方法 Active CN106548462B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610948742.6A CN106548462B (zh) 2016-11-02 2016-11-02 基于薄板样条插值的非线性sar图像几何校正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610948742.6A CN106548462B (zh) 2016-11-02 2016-11-02 基于薄板样条插值的非线性sar图像几何校正方法

Publications (2)

Publication Number Publication Date
CN106548462A true CN106548462A (zh) 2017-03-29
CN106548462B CN106548462B (zh) 2019-02-19

Family

ID=58393526

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610948742.6A Active CN106548462B (zh) 2016-11-02 2016-11-02 基于薄板样条插值的非线性sar图像几何校正方法

Country Status (1)

Country Link
CN (1) CN106548462B (zh)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107330856A (zh) * 2017-06-17 2017-11-07 南京理工大学 一种基于投影变换和薄板样条的全景成像方法
CN107798327A (zh) * 2017-10-31 2018-03-13 北京小米移动软件有限公司 字符识别方法及装置
CN108346152A (zh) * 2018-03-19 2018-07-31 北京大学口腔医院 基于牙根尖周片自动评价牙科临床诊疗效果的方法
CN108871185A (zh) * 2018-05-10 2018-11-23 苏州大学 零件检测的方法、装置、设备以及计算机可读存储介质
CN109615584A (zh) * 2018-12-17 2019-04-12 辽宁工程技术大学 一种基于单应性约束的sar图像序列map超分辨率重建方法
CN109903352A (zh) * 2018-12-24 2019-06-18 中国科学院遥感与数字地球研究所 一种卫星遥感影像大区域无缝正射影像制作方法
CN110031898A (zh) * 2019-05-14 2019-07-19 中国科学院地质与地球物理研究所 数据优化方法及积分法叠前深度偏移方法
CN110047032A (zh) * 2019-03-01 2019-07-23 中国人民解放军军事科学院国防科技创新研究院 一种基于径向基函数拟合的局部自适应误匹配点剔除方法
CN110390338A (zh) * 2019-07-10 2019-10-29 武汉大学 一种基于非线性引导滤波与比率梯度的sar高精度匹配方法
CN110738693A (zh) * 2019-09-27 2020-01-31 北京理工大学 一种地基成像雷达多角度图像配准方法
CN111160098A (zh) * 2019-11-21 2020-05-15 长春理工大学 一种基于sift特征的表情变化人脸识别方法
CN112149442A (zh) * 2020-09-15 2020-12-29 浙江大华技术股份有限公司 畸变二维码的识别方法和装置、存储介质及电子装置
CN112200731A (zh) * 2019-07-08 2021-01-08 上海隽珑信息技术有限公司 一种图像像素的校正方法
CN112330727A (zh) * 2020-11-02 2021-02-05 珠海大横琴科技发展有限公司 图像匹配方法、装置、计算机设备及存储介质
CN113658080A (zh) * 2021-08-23 2021-11-16 宁波棱镜空间智能科技有限公司 一种基于特征点匹配的线扫圆柱体几何校正方法及装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103839253A (zh) * 2013-11-21 2014-06-04 苏州盛景空间信息技术有限公司 一种基于局部仿射变换的任意点匹配方法
CN104517124A (zh) * 2014-12-25 2015-04-15 西安电子科技大学 基于sift特征点的sar图像变化检测方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103839253A (zh) * 2013-11-21 2014-06-04 苏州盛景空间信息技术有限公司 一种基于局部仿射变换的任意点匹配方法
CN104517124A (zh) * 2014-12-25 2015-04-15 西安电子科技大学 基于sift特征点的sar图像变化检测方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
FLORA DELLINGER ETC.: ""SAR-SIFT:A SIFT-Like Algorithm for SAR Images"", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》 *
YUHAO ZHOU ETC.: ""AN AUTOMATIC GLOBAL-TO-LOCAL IMAGE REGISTRATION BASED ON SIFT AND THIN-PLATE SPLINE"", 《2013 IEEE INTERNATIONAL GEOSCIENCE AND REMOTE SENSING SYMPOSIUM》 *
万国挺: ""图像拼接技术与质量评价方法研究"", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *

Cited By (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107330856B (zh) * 2017-06-17 2020-11-13 南京理工大学 一种基于投影变换和薄板样条的全景成像方法
CN107330856A (zh) * 2017-06-17 2017-11-07 南京理工大学 一种基于投影变换和薄板样条的全景成像方法
CN107798327A (zh) * 2017-10-31 2018-03-13 北京小米移动软件有限公司 字符识别方法及装置
CN108346152A (zh) * 2018-03-19 2018-07-31 北京大学口腔医院 基于牙根尖周片自动评价牙科临床诊疗效果的方法
CN108871185A (zh) * 2018-05-10 2018-11-23 苏州大学 零件检测的方法、装置、设备以及计算机可读存储介质
CN108871185B (zh) * 2018-05-10 2020-12-29 苏州大学 零件检测的方法、装置、设备以及计算机可读存储介质
CN109615584A (zh) * 2018-12-17 2019-04-12 辽宁工程技术大学 一种基于单应性约束的sar图像序列map超分辨率重建方法
CN109615584B (zh) * 2018-12-17 2022-12-06 辽宁工程技术大学 一种基于单应性约束的sar图像序列map超分辨率重建方法
CN109903352A (zh) * 2018-12-24 2019-06-18 中国科学院遥感与数字地球研究所 一种卫星遥感影像大区域无缝正射影像制作方法
CN110047032A (zh) * 2019-03-01 2019-07-23 中国人民解放军军事科学院国防科技创新研究院 一种基于径向基函数拟合的局部自适应误匹配点剔除方法
CN110031898A (zh) * 2019-05-14 2019-07-19 中国科学院地质与地球物理研究所 数据优化方法及积分法叠前深度偏移方法
CN110031898B (zh) * 2019-05-14 2019-12-13 中国科学院地质与地球物理研究所 数据优化方法及积分法叠前深度偏移方法
CN112200731A (zh) * 2019-07-08 2021-01-08 上海隽珑信息技术有限公司 一种图像像素的校正方法
CN110390338B (zh) * 2019-07-10 2022-08-05 武汉大学 一种基于非线性引导滤波与比率梯度的sar高精度匹配方法
CN110390338A (zh) * 2019-07-10 2019-10-29 武汉大学 一种基于非线性引导滤波与比率梯度的sar高精度匹配方法
CN110738693A (zh) * 2019-09-27 2020-01-31 北京理工大学 一种地基成像雷达多角度图像配准方法
CN110738693B (zh) * 2019-09-27 2022-05-20 北京理工大学 一种地基成像雷达多角度图像配准方法
CN111160098A (zh) * 2019-11-21 2020-05-15 长春理工大学 一种基于sift特征的表情变化人脸识别方法
CN112149442A (zh) * 2020-09-15 2020-12-29 浙江大华技术股份有限公司 畸变二维码的识别方法和装置、存储介质及电子装置
CN112330727A (zh) * 2020-11-02 2021-02-05 珠海大横琴科技发展有限公司 图像匹配方法、装置、计算机设备及存储介质
CN113658080A (zh) * 2021-08-23 2021-11-16 宁波棱镜空间智能科技有限公司 一种基于特征点匹配的线扫圆柱体几何校正方法及装置
CN113658080B (zh) * 2021-08-23 2023-12-22 宁波棱镜空间智能科技有限公司 一种基于特征点匹配的线扫圆柱体几何校正方法及装置

Also Published As

Publication number Publication date
CN106548462B (zh) 2019-02-19

Similar Documents

Publication Publication Date Title
CN106548462A (zh) 基于薄板样条插值的非线性sar图像几何校正方法
CN107301661B (zh) 基于边缘点特征的高分辨率遥感图像配准方法
CN104318548B (zh) 一种基于空间稀疏度和sift特征提取的快速图像配准实现方法
CN105744256B (zh) 基于图谱视觉显著的立体图像质量客观评价方法
CN110992263B (zh) 一种图像拼接方法及系统
CN102800096B (zh) 一种摄像机参数的鲁棒性估计算法
CN106709948A (zh) 一种基于超像素分割的快速双目立体匹配方法
CN107560592B (zh) 一种用于光电跟踪仪联动目标的精确测距方法
CN104732546B (zh) 区域相似性和局部空间约束的非刚性sar图像配准方法
CN105354841B (zh) 一种快速遥感影像匹配方法及系统
CN101639947A (zh) 基于图像的植物三维形态测量与重建方法及系统
CN103679714A (zh) 一种基于梯度互相关的光学和sar图像自动配准方法
CN103971378A (zh) 一种混合视觉系统中全景图像的三维重建方法
CN108647580B (zh) 基于改进sift引导isar图像特征点提取匹配方法
CN105866776B (zh) 一种地基sar的高质量动态ps点选择方法
CN103886606B (zh) 基于联合广义伽玛分布参数的sar图像分割方法
CN105957096A (zh) 一种用于三维数字图像相关的相机外参标定方法
CN104103047B (zh) 一种心电图像倾斜度校正方法
CN108765476A (zh) 一种偏振图像配准方法
CN107341824B (zh) 一种图像配准的综合评价指标生成方法
CN105488541A (zh) 增强现实系统中基于机器学习的自然特征点识别方法
CN112329726B (zh) 人脸识别方法及装置
CN109461132A (zh) 基于特征点几何拓扑关系的sar图像自动配准方法
CN103646395A (zh) 一种基于格网法的高精度影像配准方法
CN105631872A (zh) 基于多特征点的遥感图像配准方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant