CN105787943A - 基于多尺度图像块特征和稀疏表示的sar图像配准方法 - Google Patents
基于多尺度图像块特征和稀疏表示的sar图像配准方法 Download PDFInfo
- Publication number
- CN105787943A CN105787943A CN201610118533.9A CN201610118533A CN105787943A CN 105787943 A CN105787943 A CN 105787943A CN 201610118533 A CN201610118533 A CN 201610118533A CN 105787943 A CN105787943 A CN 105787943A
- Authority
- CN
- China
- Prior art keywords
- image
- point
- points
- reference image
- feature
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 50
- 230000009466 transformation Effects 0.000 claims abstract description 7
- 230000002159 abnormal effect Effects 0.000 claims abstract description 5
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 claims abstract description 4
- 238000000354 decomposition reaction Methods 0.000 claims description 17
- 239000013598 vector Substances 0.000 claims description 10
- 238000004422 calculation algorithm Methods 0.000 claims description 6
- 238000013507 mapping Methods 0.000 claims description 3
- 238000012216 screening Methods 0.000 claims description 3
- 230000008859 change Effects 0.000 abstract description 3
- 230000000694 effects Effects 0.000 abstract description 3
- 238000001514 detection method Methods 0.000 abstract description 2
- 230000004927 fusion Effects 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000002708 enhancing effect Effects 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/14—Transformations for image registration, e.g. adjusting or mapping for alignment of images
- G06T3/147—Transformations for image registration, e.g. adjusting or mapping for alignment of images using affine transformations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
- G06T2207/10044—Radar image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20048—Transform domain processing
- G06T2207/20064—Wavelet transform [DWT]
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种基于多尺度图像块特征和稀疏表示的SAR图像配准方法,主要解决现有配准方法应用于SAR图像配准效果不佳的问题,其实现的步骤是:1)输入两幅SAR图像,任选一幅作为参考图像,将另一幅作为待配准图像;2)选取参考图像特征点;3)利用多尺度图像块特征构建参考图像和待配准图像的特征点描述符;4)建立参考图像和待配准图像间的匹配点对;5)去除匹配点对中异常点;6)根据最终得到的匹配点对,建立仿射变换模型,采用最小二乘法获得几何形变参数,得到配准结果。本发明与现有技术相比,增强了对斑点噪声的鲁棒性,提高了匹配点对的准确性以及配准精度,可用于图像融合和变化检测。
Description
技术领域
本发明属于图像处理技术领域,具体涉及雷达图像处理中的图像配准方法,可用于图像融合和变化检测。
背景技术
合成孔径雷达SAR***因其全天时,全天候,具有穿透性等特点而广泛应用于军事和民用邻域。SAR图像配准作为SAR图像应用中的关键环节,它是对取自不同时间、不同视角的同一景物的两幅或多幅SAR图像进行匹配、叠加的过程。
对于图像配准问题,目前提出的方法大致可以分为两类:基于灰度和基于特征的配准方法。基于灰度的配准方法直接利用图像的灰度信息,通过建立图像像素间的某种相似性度量测度来寻找到最优匹配时对应的平移、旋转等配准参数。最常用的基于灰度的配准方法是基于互信息的配准方法。这种方法虽然比较直观,容易实现,但计算复杂度高,容易陷入局部最优解,而且易受噪声的影响。基于特征的配准方法由于不直接作用于图像灰度值,而是作用于图像本身的特征,因而对灰度变化有较强的适应能力,且计算量小,能够处理图像间的配准问题。最常用的基于特征的配准方法是基于尺度不变特征变换SIFT特征的配准方法。然而由于SAR图像中存在斑点噪声,基于特征的配准方法在处理SAR图像配准时,很有可能将斑点噪声检测为特征点,从而会带来大量的错误匹配点,导致错误的配准结果。
发明内容
本发明的目的在于提出一种基于多尺度图像块特征和稀疏表示的SAR图像配准方法,以解决现有技术中进行SAR图像配准时出现大量错误匹配点,导致配准精度不高的问题。
实现本发明目的的技术思路是:利用空间相关性选择可靠性强的特征点,采用多尺度图像块特征形成特征描述符,根据稀疏表示技术计算的最小差异准则获取最佳的匹配点对,有效地增强对斑点噪声的鲁棒性,其实现步骤包括如下:
(1)输入两幅图像,任选一幅作为参考图像I1,将另一幅作为待配准图像I2;
(2)选取参考图像特征点:
(2a)采用SIFT算法提取参考图像I1的特征点,并将I1的所有特征点存放在第一集合R中;
(2b)从参考图像特征点集合R中任意选取一个特征点ri,利用平稳小波变换方法计算每个特征点对应的空间相关性ρ(ri);
(2c)设置阈值E=0.05,若得到的ρ(ri)满足ρ(ri)≥E,则将参考图像特征点ri保留,否则,删除该特征点;
(2d)遍历参考图像所有特征点,重复步骤(2b)-(2c),得到筛选后的参考图像特征点;
(2e)分别计算上述筛选后参考图像特征点集合中任意两个特征点间的欧式距离Ed,若Ed≥15,则保留这两个特征点,否则去除;
(2f)将步骤(2e)得到的特征点集合中的前10个特征点作为最终的参考图像特征点;
(3)利用多尺度图像块特征构建参考图像特征点描述符:
(3a)任意选取一个参考图像特征点a,取该特征点周围邻域15×15的图像块P(a);
(3b)采用平稳小波变换对图像块P(a)进行多尺度分解,得到三个不同分解尺度的图像块Ps(a),s=3,4,5;
(3c)分别计算上述三个不同分解尺度图像块的灰度直方图矢量Hs(a),将其作为特征描述符的灰度特征;
(3d)分别计算上述三个不同分解尺度图像块的梯度方向直方图矢量Gs(a),将其作为特征描述符的梯度特征;
(3e)将不同分解尺度图像块对应的灰度特征和梯度特征串联起来,得到参考图像特征点a对应的特征描述符F(a)={H3(a),H4(a),H5(a),G3(a),G4(a),G5(a)};
(4)对于任意一个待配准图像特征点b,根据步骤(3)同样的操作得到其特征点描述符F(b),利用参考图像特征点描述符F(a)与待配准图像特征点描述符F(b)间的相似性,建立参考图像和待配准图像间的匹配点对;
(5)去除步骤(4)得到的匹配点对中的异常点:
(5a)从步骤(4)得到的匹配点对中任意选取一个参考图像特征点rc,在已匹配的参考图像特征点中取该特征点最近邻的3个邻域点,并将所取的这3个邻域点映射到与参考图像特征点rc成为匹配点的待配准图像特征点tc的邻域内,计算该参考图像特征点rc与其邻域点间的几何代价
其中,表示参考图像特征点rc的第k个最近邻点,tc表示与参考图像特征点rc相匹配的待配准图像特征点,m(·)表示匹配函数,||·||表示欧式距离,c表示匹配点对的索引,其取值范围为1到10,k表示第c个参考图像特征点所取的最近邻点的索引,其取值范围为1到3;
(5b)遍历所有已匹配参考图像特征点,重复步骤(5a),得到已匹配参考图像特征点与其各自的邻域点间的几何代价,将几何代价最小值对应的所有特征点作为基准点集合,表示为:
其中,(rc,tc)表示匹配点对,tm(k)表示与邻域点对应的匹配点;
(5c)采用下式计算剩余匹配点到基准点间的几何代价:
其中,rc′和tc′分别表示已匹配点对中剩余的参考图像特征点和待配准图像特征点,表示基准点集合中第o个参考图像特征点,tm(o)表示与集准点集合中第o个参考图像特征点对应的匹配点,c′表示剩余匹配点对的索引,其取值范围为1到6,o表示基准点的索引,其取值范围为1到4;
(5d)设置阈值Eo=0.03,若得到的几何代价满足则将(rc′,tc′)作为正确的匹配点,否则,将该匹配点对删除;
(5e)重复步骤(5c)-(5d),得到参考图像和待配准图像最终的匹配点对;
(6)根据上述得到的最终匹配点对,建立仿射变换模型,计算待配准图像的几何形变参数,并利用该几何形变参数,将待配准图像进行几何变换,得到配准结果。
本发明与现有技术相比具有以下优点:
第一,由于本发明在对SAR图像配准的过程中,采用空间相关性选择特征点,同时利用多尺度图像块特征构建特征描述符,克服了现有技术仅采用单一尺度图像块信息无法准确描述特征点属性的不足,使得本发明提高了特征点的显著性,增强了对斑点噪声的鲁棒性。
第二,由于本发明利用基于稀疏表示技术的最小差异准则建立匹配点对,同时采用匹配点与其邻域点间的几何约束关系滤除特征异常点,克服了现有技术采用欧式距离比方法在建立匹配点对的过程中容易出现错误匹配点的不足,使得本发明提高了匹配点的准确性。
附图说明
图1是本发明的实现流程图;
图2是本发明的第一组仿真实验结果图;
图3是本发明的第二组仿真实验结果图。
具体实施方式
下面结合附图对本发明做进一步描述:
参照图1,本发明的实施步骤如下:
步骤1,输入两幅图像,任选一幅作为参考图像I1,将另一幅作为待配准图像I2。
输入的两幅图像是在某机载雷达获取的不同极化方式或不同时相两幅SAR图像上分别截取的。
步骤2,选取参考图像特征点。
2.1)采用SIFT算法提取参考图像I1的特征点,将I1的所有特征点存放在第一集合R中;
2.2)从参考图像特征点集合R中任意选取一个特征点ri,利用平稳小波变换方法计算每个特征点对应的空间相关性ρ(ri);
2.2a)利用平稳小波变换对输入图像进行s个尺度分解,得到输入图像在不同尺度上3种不同的细节图像;
2.2b)定义任意图像像素点x在s尺度下的幅值Ms(x)表示为:
其中,分别表示在s尺度下图像在水平方向,垂直方向和对角线方向上的细节信息,|·|表示绝对值操作运算;
2.2c)采用下式计算像素点x的空间相关性ρ(x):
其中,∏(·)表示相乘操作;
2.3)设置阈值E=0.05,若得到的ρ(ri)满足ρ(ri)≥E,则将参考图像特征点ri保留,否则,删除该特征点;
2.4)遍历参考图像所有特征点,重复步骤2.2)-2.3),得到筛选后的参考图像特征点;
2.5)分别计算上述筛选后参考图像特征点集合中任意两个特征点间的欧式距离Ed,若Ed≥15,则保留这两个特征点,否则去除;
2.6)将步骤2.5)得到的特征点集合中的前10个特征点作为最终的参考图像特征点。
步骤3:利用多尺度图像块特征构建参考图像特征点描述符。
3.1)任意选取一个参考图像特征点a,取该特征点周围邻域15×15的图像块P(a);
3.2)采用平稳小波变换对图像块P(a)进行多尺度分解,得到三个不同分解尺度的图像块Ps(a),s=3,4,5;
3.3)分别计算上述三个不同分解尺度图像块的灰度直方图矢量Hs(a),将其作为特征点描述符的灰度特征;
3.4)分别计算上述三个不同分解尺度图像块的梯度方向直方图矢量Gs(a),将其作为特征点描述符的梯度特征;
3.5)将不同分解尺度图像块对应的灰度特征和梯度特征串联起来,得到参考图像特征点a对应的特征描述符F(a)={H3(a),H4(a),H5(a),G3(a),G4(a),G5(a)}。
步骤4:对于任意一个待配准图像特征点b,根据步骤(3)同样的操作得到其特征点描述符F(b),利用参考图像特征点描述符F(a)与待配准图像特征点描述符F(b)间的相似性,建立参考图像和待配准图像间的匹配点对。
形成匹配点对的方法有很多种,常见的有最近邻方法,欧式距离比等,本实例中采用但不限于基于稀疏表示技术的最小差异准则建立匹配点对,其具体实施步骤如下:
4.1)从步骤(2)得到的参考图像特征点中任意选取一个特征点ri,计算其对应的特征描述符F(ri);
4.2)设参考图像和待配准图像间的最大的偏移量为l=100,在待配准图像中选取L×L大小的窗口Wi作为特征点ri的搜索区域,并将该区域内所有像素点{Vq:q=1,2,…,Q,Q=L×L}均作为特征点ri的候选匹配点,其中L=2×l+1;
4.3)对于上述的每一个候选匹配点Vq,取该像素点周围邻域20×20的窗口,并将该窗口内的所有像素点un对应的特征描述符F(un)作为Vq对应的稀疏字典Dq={F(un),n=1,2,…,J,J=20×20};
4.4)利用正交匹配追踪算法求解如下式子,得到特征点ri在稀疏字典Dq下的稀疏向量αq
其中,argmin(·)表示函数达到最小值时自变量的取值,||·||表示欧式距离,||·||0表示零范数,C表示稀疏度;
4.5)遍历窗口Wi内所有像素点,重复步骤4.3)-4.4),得到特征点ri在Q个不同的稀疏字典下的稀疏向量αq;
4.6)按下式计算特征点ri与任意候选匹配点Vq间的重构误差,并将具有最小重构误差的像素点作为特征点ri的匹配点m(ri)
4.7)遍历所有参考图像特征点,重复步骤4.1)-4.6),得到所有参考图像特征点在待配准图像内对应的匹配点,将待配准图像的所有特征点存放在第二集合T中。
步骤5:去除步骤(4)得到的匹配点对中的异常点。
5.1)从步骤(4)得到的匹配点对中任意选取一个参考图像特征点rc,在已匹配的参考图像特征点中取该特征点最近邻的3个邻域点,并将所取的这3个邻域点映射到与参考图像特征点rc成为匹配点的待配准图像特征点tc的邻域内,计算该参考图像特征点rc与其邻域点间的几何代价
其中,表示参考图像特征点rc的第k个最近邻点,tc表示与参考图像特征点rc相匹配的待配准图像特征点,m(·)表示匹配函数,||·||表示欧式距离,c表示匹配点对的索引,其取值范围为1到10,k表示第c个参考图像特征点所取的最近邻点的索引,其取值范围为1到3;
5.2)遍历所有已匹配参考图像特征点,重复步骤5.1),得到已匹配参考图像特征点与其各自的邻域点间的几何代价,将几何代价最小值对应的所有特征点作为基准点集合,表示为:
其中,(rc,tc)表示匹配点对,tm(k)表示与邻域点对应的匹配点;
5.3)采用下式计算剩余匹配点到基准点间的几何代价:
其中,rc′和tc′分别表示已匹配点对中剩余的参考图像特征点和待配准图像特征点,表示基准点集合中第o个参考图像特征点,tm(o)表示与集准点集合中第o个参考图像特征点对应的匹配点,c′表示剩余匹配点对的索引,其取值范围为1到6,o表示基准点的索引,其取值范围为1到4;
5.4)设置阈值Eo=0.03,若得到的几何代价满足则将(rc′,tc′)作为正确的匹配点,否则,将该匹配点对删除;
5.5)重复步骤5.3)-5.4),得到参考图像和待配准图像最终的匹配点对。
步骤6:根据最终匹配点对,得到配准结果。
根据上述得到的最终匹配点对建立仿射变换模型,计算待配准图像的几何形变参数,本实例具体采用最小二乘法计算几何形变参数,并利用该几何形变参数,将待配准图像进行几何变换,得到配准结果。
下面结合实验仿真对本发明效果做进一步说明。
1.仿真条件:
本发明的仿真实验平台采用Intel(R)Pentium(R)CPUG32403.10GHz,内存为4GB,运行Windows7的PC机,编程语言为Matlab2011b。
2.仿真内容与结果分析:
仿真1、分别应用基于SIFT-OCT的图像配准方法、基于BFSIFT的图像配准方法,基于NDSS-SIFT的图像配准方法和本发明对SAR图像进行配准,结果如图2所示,其中图2(a)是基于SIFT-OCT方法的匹配点连线图,图2(b)是基于BFSIFT方法的匹配点连线图,图2(c)是基于NDSS-SIFT方法的匹配点连线图,图2(d)是本发明的匹配点连线图。图2中黄色实线表示正确匹配点对,红色线表示错误匹配点对。
仿真2、分别应用基于SIFT-OCT的图像配准方法、基于BFSIFT的图像配准方法,基于NDSS-SIFT的图像配准方法和本发明对SAR图像进行配准,结果如图3所示,其中图3(a)是基于SIFT-OCT方法的匹配点连线图,图3(b)是基于BFSIFT方法的匹配点连线图,图3(c)是基于NDSS-SIFT方法的匹配点连线图,图3(d)是本发明的匹配点连线图。图3中黄色实线表示正确匹配点对,红色线表示错误匹配点对。
从图2(a)-2(c)和图3(a)-3(c)可以看出,对于2组实测SAR图像对,SIFT-OCT方法,BFSIFT方法和NDSS-SIFT方法3种算法都存在比较多的错误匹配点。其中,SIFT-OCT方法效果最差,包含的错误匹配点最多。BFSIFT方法和NDSS-SIFT方法中的错误匹配点对相对较少。
从图2(d)和图3(d)可以看出,对于2组实测SAR图像对,本发明得到的配准结果更准确,不包含错误匹配点。这源于本发明利用空间相关性选取可靠的特征点,同时利用多尺度图像块特征构建特征描述符,增强了特征描述符的显著性和对斑点噪声的鲁棒性。此外,利用基于稀疏表示技术的最小差异准则建立匹配点对,提高了匹配点对的准确性,从而解决了现有技术应用于SAR图像配准时出现大量错误匹配点对的问题。
Claims (3)
1.一种基于多尺度图像块特征和稀疏表示的SAR图像配准方法,包括如下步骤:
(1)输入两幅图像,任选一幅作为参考图像I1,将另一幅作为待配准图像I2;
(2)选取参考图像特征点:
(2a)采用SIFT算法提取参考图像I1的特征点,并将I1的所有特征点存放在第一集合R中;
(2b)从参考图像特征点集合R中任意选取一个特征点ri,利用平稳小波变换方法计算每个特征点对应的空间相关性ρ(ri);
(2c)设置阈值E=0.05,若得到的ρ(ri)满足ρ(ri)≥E,则将参考图像特征点ri保留,否则,删除该特征点;
(2d)遍历参考图像所有特征点,重复步骤(2b)-(2c),得到筛选后的参考图像特征点;
(2e)分别计算上述筛选后参考图像特征点集合中任意两个特征点间的欧式距离Ed,若Ed≥15,则保留这两个特征点,否则去除;
(2f)将步骤(2e)得到的特征点集合中的前10个特征点作为最终的参考图像特征点;
(3)利用多尺度图像块特征构建参考图像特征点描述符:
(3a)任意选取一个参考图像特征点a,取该特征点周围邻域15×15的图像块P(a);
(3b)采用平稳小波变换对图像块P(a)进行多尺度分解,得到三个不同分解尺度的图像块Ps(a),s=3,4,5;
(3c)分别计算上述三个不同分解尺度图像块的灰度直方图矢量Hs(a),将其作为特征描述符的灰度特征;
(3d)分别计算上述三个不同分解尺度图像块的梯度方向直方图矢量Gs(a),将其作为特征描述符的梯度特征;
(3e)将不同分解尺度图像块对应的灰度特征和梯度特征串联起来,得到参考图像特征点a对应的特征描述符F(a)={H3(a),H4(a),H5(a),G3(a),G4(a),G5(a)};
(4)对于任意一个待配准图像特征点b,根据步骤(3)同样的操作得到其特征点描述符F(b),利用参考图像特征点描述符F(a)与待配准图像特征点描述符F(b)间的相似性,建立参考图像和待配准图像间的匹配点对;
(5)去除步骤(4)得到的匹配点对中的异常点:
(5a)从步骤(4)得到的匹配点对中任意选取一个参考图像特征点rc,在已匹配的参考图像特征点中取该特征点最近邻的3个邻域点,并将所取的这3个邻域点映射到与参考图像特征点rc成为匹配点的待配准图像特征点tc的邻域内,计算该参考图像特征点rc与其邻域点间的几何代价
其中,表示参考图像特征点rc的第k个最近邻点,tc表示与参考图像特征点rc相匹配的待配准图像特征点,m(·)表示匹配函数,||·||表示欧式距离,c表示匹配点对的索引,其取值范围为1到10,k表示第c个参考图像特征点所取的最近邻点的索引,其取值范围为1到3;
(5b)遍历所有已匹配参考图像特征点,重复步骤(5a),得到已匹配参考图像特征点与其各自的邻域点间的几何代价,将几何代价最小值对应的所有特征点作为基准点集合,表示为:
其中,(rc,tc)表示匹配点对,tm(k)表示与邻域点对应的匹配点;
(5c)采用下式计算剩余匹配点到基准点间的几何代价:
其中,rc′和tc′分别表示已匹配点对中剩余的参考图像特征点和待配准图像特征点,表示基准点集合中第o个参考图像特征点,tm(o)表示与集准点集合中第o个参考图像特征点对应的匹配点,c′表示剩余匹配点对的索引,其取值范围为1到6,o表示基准点的索引,其取值范围为1到4。
(5d)设置阈值Eo=0.03,若得到的几何代价满足则将(rc′,tc′)作为正确的匹配点,否则,将该匹配点对删除;
(5e)重复步骤(5c)-(5d),得到参考图像和待配准图像最终的匹配点对;
(6)根据上述得到的最终匹配点对,建立仿射变换模型,计算待配准图像的几何形变参数,并利用该几何形变参数,将待配准图像进行几何变换,得到配准结果。
2.根据权利要求1所述的的SAR图像配准方法,其特征在于步骤(2b)中利用平稳小波变换方法计算每个特征点对应的空间相关性,按如下步骤进行:
2b1)利用平稳小波变换对输入图像进行s个尺度分解,得到输入图像在不同尺度上3种不同的细节图像;
2b2)定义任意图像像素点x在s尺度下的幅值Ms(x)表示为:
其中,分别表示在s尺度下图像在水平方向,垂直方向和对角线方向上的细节信息,|·|表示绝对值操作运算;
2b3)采用下式计算像素点x的空间相关性ρ(x):
其中,∏(·)表示相乘操作。
3.根据权利要求1所述的基于多尺度图像块特征和稀疏表示的SAR图像配准方法,其特征在于所述步骤(4)中建立参考图像和待配准图像间的匹配点对,按如下步骤进行:
4a)从步骤(2)得到的参考图像特征点中任意选取一个特征点ri,计算其对应的特征描述符F(ri);
4b)设参考图像和待配准图像间的最大的偏移量为l=100,在待配准图像中选取L×L大小的窗口Wi作为特征点ri的搜索区域,将该区域内所有像素点{Vq:q=1,2,…,Q,Q=L×L}均作为特征点ri的候选匹配点,其中L=2×l+1;
4c)对于上述的每一个候选匹配点Vq,取该像素点周围邻域20×20的窗口,并将该窗口内的所有像素点un对应的特征描述符F(un)作为Vq对应的稀疏字典Dq={F(un),n=1,2,…,J,J=20×20};
4d)利用正交匹配追踪算法求解如下式子,得到特征点ri在稀疏字典Dq下的稀疏向量αq
其中,argmin(·)表示函数达到最小值时自变量的取值,||·||表示欧式距离,||·||0表示零范数,C表示稀疏度;
4e)遍历窗口Wi内所有像素点,重复步骤4c)-4d),得到特征点ri在Q个不同的稀疏字典下的稀疏向量αq。
4f)按下式计算特征点ri与任意候选匹配点Vq间的重构误差,并将具有最小重构误差的像素点作为特征点ri的匹配点m(ri)
4g)遍历所有参考图像特征点,重复步骤(4a)-(4f),得到所有参考图像特征点在待配准图像内对应的匹配点,将待配准图像的所有特征点存放在第二集合T中。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610118533.9A CN105787943B (zh) | 2016-03-03 | 2016-03-03 | 基于多尺度图像块特征和稀疏表示的sar图像配准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610118533.9A CN105787943B (zh) | 2016-03-03 | 2016-03-03 | 基于多尺度图像块特征和稀疏表示的sar图像配准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105787943A true CN105787943A (zh) | 2016-07-20 |
CN105787943B CN105787943B (zh) | 2018-08-31 |
Family
ID=56386825
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610118533.9A Active CN105787943B (zh) | 2016-03-03 | 2016-03-03 | 基于多尺度图像块特征和稀疏表示的sar图像配准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105787943B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108627835A (zh) * | 2018-06-29 | 2018-10-09 | 中国科学院电子学研究所 | 全极化差分sar层析的目标重构方法 |
CN109166336A (zh) * | 2018-10-19 | 2019-01-08 | 福建工程学院 | 一种基于区块链技术的实时路况信息采集推送方法 |
CN109741275A (zh) * | 2018-12-28 | 2019-05-10 | 济南大学 | 一种mvct图像的增强方法及*** |
CN109816619A (zh) * | 2019-01-28 | 2019-05-28 | 努比亚技术有限公司 | 图像融合方法、装置、终端及计算机可读存储介质 |
CN110033422A (zh) * | 2019-04-10 | 2019-07-19 | 北京科技大学 | 一种眼底oct图像融合方法及装置 |
CN112184785A (zh) * | 2020-09-30 | 2021-01-05 | 西安电子科技大学 | 基于mcd度量和vtm的多模态遥感图像配准方法 |
CN116008911A (zh) * | 2022-12-02 | 2023-04-25 | 南昌工程学院 | 一种基于新型原子匹配准则的正交匹配追踪声源识别方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103177444A (zh) * | 2013-03-08 | 2013-06-26 | 中国电子科技集团公司第十四研究所 | 一种sar图像自动配准方法 |
CN103839265A (zh) * | 2014-02-26 | 2014-06-04 | 西安电子科技大学 | 基于sift和归一化互信息的sar图像配准方法 |
-
2016
- 2016-03-03 CN CN201610118533.9A patent/CN105787943B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103177444A (zh) * | 2013-03-08 | 2013-06-26 | 中国电子科技集团公司第十四研究所 | 一种sar图像自动配准方法 |
CN103839265A (zh) * | 2014-02-26 | 2014-06-04 | 西安电子科技大学 | 基于sift和归一化互信息的sar图像配准方法 |
Non-Patent Citations (2)
Title |
---|
GUITING WANG等: "Change detection based on image segment and fusion in multitemporal SAR images", 《SYNTHETIC APERTURE RADAR, 2009. APSAR 2009. 2ND ASIAN-PACIFIC CONFERENCE ON》 * |
刘爱平等: "基于MAR-MRF的SAR图像分割方法", 《电子与信息学报》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108627835A (zh) * | 2018-06-29 | 2018-10-09 | 中国科学院电子学研究所 | 全极化差分sar层析的目标重构方法 |
CN108627835B (zh) * | 2018-06-29 | 2021-07-27 | 中国科学院电子学研究所 | 全极化差分sar层析的目标重构方法 |
CN109166336B (zh) * | 2018-10-19 | 2020-08-07 | 福建工程学院 | 一种基于区块链技术的实时路况信息采集推送方法 |
CN109166336A (zh) * | 2018-10-19 | 2019-01-08 | 福建工程学院 | 一种基于区块链技术的实时路况信息采集推送方法 |
CN109741275B (zh) * | 2018-12-28 | 2020-06-12 | 济南大学 | 一种mvct图像的增强方法及*** |
CN109741275A (zh) * | 2018-12-28 | 2019-05-10 | 济南大学 | 一种mvct图像的增强方法及*** |
CN109816619A (zh) * | 2019-01-28 | 2019-05-28 | 努比亚技术有限公司 | 图像融合方法、装置、终端及计算机可读存储介质 |
CN110033422A (zh) * | 2019-04-10 | 2019-07-19 | 北京科技大学 | 一种眼底oct图像融合方法及装置 |
CN110033422B (zh) * | 2019-04-10 | 2021-03-23 | 北京科技大学 | 一种眼底oct图像融合方法及装置 |
CN112184785A (zh) * | 2020-09-30 | 2021-01-05 | 西安电子科技大学 | 基于mcd度量和vtm的多模态遥感图像配准方法 |
CN112184785B (zh) * | 2020-09-30 | 2023-03-24 | 西安电子科技大学 | 基于mcd度量和vtm的多模态遥感图像配准方法 |
CN116008911A (zh) * | 2022-12-02 | 2023-04-25 | 南昌工程学院 | 一种基于新型原子匹配准则的正交匹配追踪声源识别方法 |
CN116008911B (zh) * | 2022-12-02 | 2023-08-22 | 南昌工程学院 | 一种基于新型原子匹配准则的正交匹配追踪声源识别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105787943B (zh) | 2018-08-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105787943B (zh) | 基于多尺度图像块特征和稀疏表示的sar图像配准方法 | |
CN107301661B (zh) | 基于边缘点特征的高分辨率遥感图像配准方法 | |
US11244197B2 (en) | Fast and robust multimodal remote sensing image matching method and system | |
CN111028277B (zh) | 基于伪孪生卷积神经网络的sar和光学遥感图像配准方法 | |
CN104318548B (zh) | 一种基于空间稀疏度和sift特征提取的快速图像配准实现方法 | |
CN104376330B (zh) | 基于超像素散射机制的极化sar图像舰船目标检测方法 | |
CN102169545B (zh) | 一种高分辨率遥感图像的变化检测方法 | |
CN107145829B (zh) | 一种融合纹理特征和尺度不变特征的掌静脉识别方法 | |
CN104732546B (zh) | 区域相似性和局部空间约束的非刚性sar图像配准方法 | |
US7912321B1 (en) | Image registration with uncertainty analysis | |
CN108021890B (zh) | 一种基于plsa和bow的高分遥感影像港口检测方法 | |
CN105354841B (zh) | 一种快速遥感影像匹配方法及*** | |
CN105631872B (zh) | 基于多特征点的遥感图像配准方法 | |
CN107862319B (zh) | 一种基于邻域投票的异源高分光学影像匹配误差剔除方法 | |
CN104077782A (zh) | 一种星载遥感图像匹配方法 | |
CN105160344A (zh) | 一种三维点云的局部特征提取方法及装置 | |
CN103729462B (zh) | 一种基于稀疏表示处理遮挡的行人检索方法 | |
CN109509216A (zh) | 一种Landsat与SPOT卫星影像自动配准方法 | |
Huang et al. | SAR and optical images registration using shape context | |
CN114119686A (zh) | 空间布局相似计算的多源遥感影像配准方法 | |
CN104484647B (zh) | 一种高分辨率遥感图像云高度检测方法 | |
Tashlinskii et al. | Synthesis of stochastic algorithms for image registration by the criterion of maximum mutual information | |
Li et al. | A new framework of hyperspectral image classification based on spatial spectral interest point | |
CN103605980A (zh) | 基于局部特征匹配的isar目标部件检测方法 | |
CN107451610A (zh) | 一种提高特征匹配精度的图像检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |