CN102436652A - 一种多源遥感图像自动配准方法 - Google Patents
一种多源遥感图像自动配准方法 Download PDFInfo
- Publication number
- CN102436652A CN102436652A CN201110254958XA CN201110254958A CN102436652A CN 102436652 A CN102436652 A CN 102436652A CN 201110254958X A CN201110254958X A CN 201110254958XA CN 201110254958 A CN201110254958 A CN 201110254958A CN 102436652 A CN102436652 A CN 102436652A
- Authority
- CN
- China
- Prior art keywords
- image
- same place
- registration
- cutting
- coupling
- 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
Images
Landscapes
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种多源遥感图像自动配准方法,包括对基准图像和待配准图像重叠区进行裁剪和对待配准图像进行小波尺度变换,使其分辨率与基准图像一致的步骤;对基准图像进行切片并进行特征点检测,提取均匀分布的特征点并获得初始匹配同名点对的步骤;利用初始匹配同名点对统计求得待配准图像经纬度方向整体偏移量的步骤;根据经纬度方向整体偏移量修正初始匹配同名点对,并对修正后特征点进行双向互相关特征点匹配获取亚像素定位精度的精匹配同名点,最终重采样得到配准图像的步骤。采用本发明有效克服了HJ星多源遥感图像大视场、大跨度、异特性等特点,可快速检测大量分布均匀、稳定、可靠的高精度匹配同名点对,并进行基于密集三角网的微分纠正。
Description
技术领域
本发明涉及一种多源遥感图像自动配准方法,属于卫星图像处理技术领域,可用于HJ星红外、CCD和高光谱同类传感器和异类传感器图像的自动配准。
背景技术
图像配准是图像处理的基本任务之一,是对取自不同时间、不同传感器或不同视角的同一景物的两幅图像或多幅图像进行匹配、叠加的过程。已经被广泛应用于图像分析、变化检测和信息融合等领域。设参考图像和待配准图像分别用I1(x1,y1)和I2(x2,y2)表示,则图像配准在数学上可以表示为
I1(x1,y1)=g(I2(T(x2,y2)))
其中,g(·)是一维灰度变换,T(·)是二维空间位置变换。由于在大多数情况下,配准后的两幅图像一般不要求灰度一致,而且配准也并不能通过一个或多个变换来达到两幅图像的完全等同,而只能做到某种程度上的近似。
HJ星多源遥感图像配准具有其自身的特殊性,主要表现在HJ星CCD相机和红外相机分别通过双相机拼接(二台±30度视场)和双面镜扫描(±29度扫描范围)实现大视场观测;高光谱成像仪幅宽50km,侧向可视范围±30度。HJ星多源遥感图像空间分辨率从30m、100m、150m到300m不等,跨度大。将低分辨率图像向与之分辨率差异较大的较高分辨率转换,要利用插值算法,增加带宽等方式,填充大量虚假信息,会引起图像失真,影响配准精度。
典型的图像配准算法大致可以分为基于像素的、基于灰度的、基于特征的和基于物理模型的方法。目前,对于多源遥感图像的配准,主要采用的是基于特征的配准方法。在基于特征的图像配准方法中,特征的选择与提取无疑对图像配准的质量有着很大的影响,其将极大地影响图像配准成功率。传统的图像配准技术要达到较高的精度,通常需要引入一定的人工干预,手工方式不但降低了***的效率,而且影响了***的应用性。目前,还没有应用于业务化***的HJ数据快速自配准方法。
如何实现特征点的自动提取与自动匹配、如何提高宽幅遥感影像匹配效率、如何采用考虑成像几何的自动配准更好满足大视场、大跨度、异特性特点,是本发明重点解决的问题。
发明内容
本发明的技术解决问题是:克服现有技术的不足,提供了一种多源遥感图像自动配准方法,该方法有效克服了HJ星多源遥感图像大视场、大跨度、异特性等特点,可自动快速检测大量分布均匀、稳定、可靠的高精度匹配同名点对,并进行基于密集三角网的微分纠正,支持HJ星红外、CCD和高光谱同类传感器和异类传感器的多源遥感图像的配准。本发明解决了传统多源影像配准过程中需要人为干预的问题,在处理初始阶段通过预设参数就可实现全流程自动配准。
本发明的技术解决方案是:一种多源遥感图像自动配准方法,包括下列步骤:
1.根据遥感数据自带的经纬度信息,对基准图I1和待配准图像I2的重叠区进行裁剪,得到基准图裁剪结果影像I1 *和待配准图裁剪结果影像I2 *。
2.根据基准图和待配准图像分辨率信息R1和R2,判断分辨率是否一致,一致则进行下一步;不一致则对待配准图裁剪结果影像I2w进行尺度变换:
(1)对待配准图裁剪结果影像I2 *进行Mallat小波分解;得到一个低频图像信号以及不同分辨率上的高频图像信号。
(2)图像经过小波分解后,待配准图裁剪结果影像I2 *被分解成沿水平、垂直和对角三个方向的一系列高频信号图像与一个低频图像,其中同一方向的各高频图像信号之间具有相似性。小波尺度变换就是利用这种相似性,进行高频外推插值;
(3)利用小波反变换重构出比原始图像分辨率更高的插值图像WD,本方法在高频外推插值时采用的插值方法为双三次插值。
3.对基准图裁剪结果影像I1 *并进行FAST特征点检测,提取均匀分布的特征点;依据2级图像的地理信息对应关系在待校正图像上定位匹配点,得到初始匹配同名点对:
(1)将基准图裁剪结果影像I1 *按照预设大小进行切片,得到图像片I1i(i=1,2,...,M,M是图像片总数)。一个线程一次对一个切片的数据I1i进行FAST特征点检测,实现特征点检测阶段的并行处理:
(2)记基准图裁剪结果影像I1 *检测的特征点为PixI1i(i=1,2,3,...,N)。GDAL遥感图像库的OpenGIS坐标转换方法提供了用来描绘坐标***(投影和基准面)以及坐标***相互之间转换的服务,本发明通过这种坐标转换结合地理定位信息在待配准裁剪图像的尺度变换图像WD上获取初始匹配同名点PixWi(i=1,2,3,...,N),N为特征点总数目。
4.选取初始匹配同名点对{PixI1i,PixWi}(i=1,2,3,...,N)的一个子集{PixI1j,PixWj}(j=1,2,3,...,M)。基于金字塔和相关系数进行特征点匹配(匹配搜索范围优于HJ影像***几何校正的精度),一个线程一次匹配一个点,达到并行计算目的,经粗差剔除后获得更准确的同名点对,统计求得待配准图像经纬度方向整体偏移量:
(1)将基准图裁剪结果影像I1 *和待配准裁剪图像的尺度变换图像WD建立3层高斯金字塔,图像金字塔从下到上对应I1 *、I11 *、I12 *和WD、WD1、WD2,从下往上每层图像缩小1/3。
(2)第1层:
(2.1)提取初始匹配同名点子集{PixI1j,PixWj}(j=1,2,3,...,M)的一个初始匹配同名点{PixI1j,PixWj}。
(2.2)在顶层影像I12 *上以PixI1j/9为中心获取127*127图像块,WD2上以PixWj/9为中心在周围5个像素范围滑窗,进行相关匹配,获取相关系数矩阵C及其峰值对应位置(p,q)和峰值处相关系数Cmax,当Cmax>0.5时,进行下一步;否则,回到步骤(2.1)计算下1个匹配点。
(2.3)在WD2上以(p,q)为中心获取127*127图像块,顶层影像I12 *上以PixI1j/9为中心在其周围滑窗,进行相关匹配,获取相关系数矩阵及其峰值对应位置(m,n),当(m,n)与PixI1j/9对应的点位置相等时,认为PixI1j/9与(p,q)为1对匹配同名点;否则,回到步骤(2.1)计算下1个匹配点。
(2.4)在顶层影像上遍历匹配同名点子集{PixI1j,PixWj}(j=1,2,3,...,M)的每一个点,进行RANSAC(Random Sample Consensus)剔除,获取可靠匹配同名点对{LvI2t,LvW2t}(t=1,2,...,T)。
(3)第2层:
(3.1)根据顶层第1层匹配同名点{LvI2t,LvW2t}(t=1,2,...,T),各点位置乘以金子塔放大倍数3,计算第2层待配准图像到基准图像的仿射变换关系;
(3.2)提取初始匹配同名点子集{PixI1j,PixWj}(j=1,2,3,...,M)的一个初始匹配同名点{PixI1j,PixWj}。
(3.3)在第2层影像I11 *上,根据放射关系计算PixI1j/3在WD1上对应的位置PixW* j,计算WD1上PixWj/3与PixW* j在x和y方向上的距离,若均小于15则进行下一步;否则,回到步骤(3.2)计算下1个匹配点。
(3.4)在第2层影像I11 *上,以PixI1j/3为中心获取127*127图像块,WD1上以PixW* j为中心在周围5个像素范围滑窗,进行相关匹配,获取相关系数矩阵C及其峰值对应位置(p,q)。
(3.5)在WD1上以(p,q)为中心获取127*127图像块,第2层影像I11 *上,以PixI1j/3为中心在其周围滑窗,进行相关匹配,获取相关系数矩阵及其峰值对应位置(m,n),当(m,n)与PixI1j/3对应的点位置相等时,认为PixI1j/3与(p,q)为1对匹配同名点;否则,回到步骤(3.2)计算下1个匹配点。
(3.6)在第二影像上遍历匹配同名点子集{PixI1j,PixWj}(j=1,2,3,...,M)的每一个点,获取匹配同名点对{LvI1t,LvW1t}(t=1,2,...,T)。
(4)第3层(最底层):
(4.1)根据顶层第2层匹配同名点{LvI1t,LvW1t}(t=1,2,...,T),各点位置乘以金子塔放大倍数3,计算第3层待配准图像到基准图像的仿射变换关系;
(4.2)提取初始匹配同名点子集{PixI1j,PixWj}(j=1,2,3,...,M)的一个初始匹配同名点{PixI1j,PixWj}。
(4.3)在第3层影像I1 *上,根据放射关系计算PixI1j在WD上对应的位置PixW* j,计算WD1上PixWj与PixW* j在x和y方向上的距离,若均小于15则进行下一步;否则,回到步骤(4.2)计算下1个匹配点。
(4.4)在第3层影像I1 *上,以PixI1j中心获取127*127图像块,WD上以PixW* j为中心在周围5个像素范围滑窗,进行相关匹配,获取相关系数矩阵C及其峰值对应位置(p,q)。
(4.5)在WD上以(p,q)为中心获取127*127图像块,在第3层影像I1 *上,以PixI1j为中心在其周围滑窗,进行相关匹配,获取相关系数矩阵及其峰值对应位置(m,n),当(m,n)与PixI1j对应的点位置相等时,认为PixI1j与(p,q)为1对匹配同名点;否则,回到步骤(3.2)计算下1个匹配点。
(4.6)在第3影像上遍历匹配同名点子集{PixI1j,PixWj}(j=1,2,3,...,M)的每一个点,获取匹配同名点对{LvIt,LvWt}(t=1,2,...,T)。
(5)统计经纬度方向整体偏移量Δx,Δy:由匹配过程可知,最底层匹配同名点{LvIt,LvWt}(t=1,2,...,T)中的LvIt(t=1,2,...,T)为始匹配同名点子集{PixI1j,PixWj}(j=1,2,3,...,M)中PixI1j的一个子集,可得到LvIt(t=1,2,...,T)对应的初始匹配点PixWt(t=1,2,...,T),计算LvWt(t=1,2,...,T)与PixWt(t=1,2,...,T)在经纬度方向(对应图像横轴x和纵轴y方向)的均值偏移量Δx和Δy。
5.根据经纬度方向整体偏移量Δx,Δy修正初始匹配同名点对{PixI1i,PixWi}(i=1,2,3,...,N)中待配准图对应匹配点PixWi(i=1,2,3,...,N),其x坐标减去Δx,y坐标减去Δy;得到修正后匹配同名点对{PixI1i,PixW* i}(i=1,2,3,...,N)。
6.对进行双向互相关特征点匹配,匹配峰值处采用二元三点插值获取亚像素定位精度的精匹配同名点;同时对精匹配同名点对进行多项式拟合验证,剔除显著误匹配点:
(1)提取修正后匹配同名点对{PixI1i,PixW* i}(i=1,2,3,...,N)的一个匹配同名点{PixI1i,PixW* i}
(2)在基准图裁剪结果影像I1 *上以PixI1j为中心获取127*127图像块,待配准裁剪图像的尺度变换图像WD上以PixW* i为中心在周围5个像素范围滑窗,进行相关匹配,获取相关系数矩阵C及其峰值对应位置(p,q)和峰值处相关系数Cmax,当Cmax>0.5时,进行下一步;否则,回到步骤(1)计算下1个匹配点。
(3)在WD上以(p,q)为中心获取127*127图像块,基准图裁剪结果影像I1 *上以PixI1j/9为中心在其周围滑窗,进行相关匹配,获取相关系数矩阵及其峰值对应位置(m,n),当(m,n)与PixI1j对应的点位置相等时,认为PixI1j与(p,q)为1对匹配同名点,进行下一步;否则,回到步骤(1)计算下1个匹配点。
(4)对相关系数矩阵C峰值对应位置(p,q)周围3*3邻域内采用二元三点插值算法对相关系数进行插值,插值间隔选取1/20,获取插值后3*3邻域相关系数矩阵Z,获取Z的对应的峰值位置(p*,q*),此时PixI1j与(p*,q*)为一对高精度匹配同名点。
(5)回到步骤(1),遍历修正后匹配同名点对{PixI1i,PixW* i}(i=1,2,3,...,N)的每一个同名点,得到亚像素匹配同名点对{PixID1i,PixWDi}(i=1,2,3,...,M)。
(6)对精匹配同名点对{PixID1i,PixWDi}(i=1,2,3,...,M)进行RANSAC剔除,获取可靠精匹配同名点对{PixID1i,PixWDi}(i=1,2,3,...,D)。
7.结合获取的大量精匹配同名点对{PixID1i,PixWDi}(i=1,2,3,...,D),构建密集Delaunay三角网,对各小面元分别进行纠正,得到高精度配准图像。
本发明与现有技术相比的有益效果在于:
(1)本发明对HJ星多源遥感图像特性进行了深入分析,设计了包括尺度变换、基于成像几何的初始匹配同名点对获取、经纬度整体偏移消除、同名点对精匹配、密集三角网纠正等步骤的全自动配准流程;在整个配准过程中考虑了HJ星配准数据量和计算性能,在特征点检测和匹配阶段均采用并行设计,提高了配准效率,满足业务化运行需求。
(2)HJ星多源遥感图像配准误差主要是经纬度方向误差,本发明结合成像几何对配准误差的影响,针对经纬度方向的整体误差采用了金字塔与互相关双向匹配的方法,提高了初始匹配精度,减小了后续匹配计算量。
(3)本发明针对同名点精匹配,采用互相关双向匹配与二元三点插值相结合的方法,将匹配精度提高到了子像元级。
(4)本发明结合HJ星几何纠正数据配准误差主要是平移误差的特点,具有同名点对分布均匀、稳定、可靠,定位精度高等优势,可检测大量匹配同名点,支持HJ星红外、CCD和高光谱同类传感器和异类传感器图像的配准,大量的试验结果表明,配准精度可达到平地优于2个像素,山地优于5个像素。
附图说明
图1为本发明的整体流程图;
图2为本发明基于金字塔和相关系数定位经纬度方向偏移量流程图;
图3为本发明特征点精匹配流程图。
具体实施方式
下面结合附图1、图2和图3对本发明的具体实施方式进行进一步的详细描述:
1.根据遥感数据自带的经纬度信息,对基准图I1和待配准图像I2的重叠区进行裁剪,得到基准图裁剪结果影像I1 *和待配准图裁剪结果影像I2 *。
2.根据基准图和待配准图像分辨率信息R1和R2,判断分辨率是否一致,一致则进行下一步;不一致则对待配准图裁剪结果影像I2 *进行尺度变换:
(1)小波分解过程表述为:令待配准图裁剪结果影像I2 *表示大小为M×N的原始图像,l(i)表示相对于分析小波的低通滤波器系数,i=0,1,2,…Nl-1,Nl表示滤波器L的支撑长度;h(i)表示相对于分析小波的高通滤波器系数,i=0,1,2,…Nh-1,Nh表示滤波器H的支撑长度,则小波分解为:
对逼近子图重复此过程,直到确定的分解水平,这样就能实现图像从高分辨率到低分辨率的小波分解。低频信号能反映出原始信号的基本波形,它是原始信号在低一级分辨率上的表示,而另一方面,高频信号则包含着由高分辨率变换到低分辨率所丢失的细节信息。如果输出的低频信号被当作新的信号输入到同一个滤波器组,再次进行离散小波变换,并且如此循环若干次,最终会得到一个低频信号以及不同分辨率上的高频细节信息。
(2)图像经过小波分解后,待配准图裁剪结果影像I2 *被分解成沿水平、垂直和对角三个方向的一系列高频信号图像与一个低频图像,其中同一方向的各高频图像信号之间具有相似性。小波尺度变换就是利用这种相似性,进行高频外推插值,然后利用小波反变换重构出比原始图像分辨率更高的插值图像WD。
不妨设某图像为一低通滤波后得到的图像IL,先通过小波变换把IL分解成IH1,IV1,ID1以及低频图像IL1,然后采用相似变换外推得到高频细节IH,IV,ID,再用小波反变换得到比原始图像分辨率更高的图像I。
I=DWT-1(IL,IH,IV,ID)
本方法在高频外推插值时采用的插值方法为双三次插值。
3.对基准图裁剪结果影像I1 *并进行FAST特征点检测,提取均匀分布的特征点;依据2级图像的地理信息对应关系在待校正图像上定位匹配点,得到初始匹配同名点对:
(1)将基准图裁剪结果影像I1 *按照预设大小进行切片,得到图像片I1i(i=1,2,...,M,M是图像片总数)。一个线程一次对一个切片的数据I1i进行FAST特征点检测,实现特征点检测阶段的并行处理:
对切片进行FAST特征点检测,表述为:检查待检测点c周围的圆,寻找其中最长的圆弧,如果圆弧中所有的点的灰度值都大于I(c)或都小于I(c),则被判定为角点,其中I(c)为c的灰度值。在离散情况下,圆大小为3×3区域,弧长为离散点数目。灰度判定公式为:
|I(c′)-I(c)|≤t
其中I(c′)为3×3圆弧上的点的灰度值,t为阈值。只要满足上述公式的连续弧长的数目大于等于9,则该中心点为角点。
在对切片I1i计算FAST特征点提取时,为了得到稀疏且分布均匀的特征点,对图像切片I1i进行20*20个像素的细粒度划分,计算每个20*20像素粒度里面的特征点响应Rj(j=1,2,...,400)。取角点响应的最大值Rmax=max(Rj),如果Rmax>0,则标记该角点为特征点。通过上述计算,可以保证基准图上成像区每20*20个像素区域内包含一个特征点。
(2)记基准图裁剪结果影像I1 *检测的特征点为PixI1i(i=1,2,3,...,N)。GDAL遥感图像库的OpenGIS坐标转换方法提供了用来描绘坐标***(投影和基准面)以及坐标***相互之间转换的服务,本发明通过这种坐标转换结合地理定位信息在待配准裁剪图像的尺度变换图像WD上获取初始匹配同名点PixWi(i=1,2,3,...,N),N为特征点总数目。具体过程表述如下:
遍历PixI1i(i=1,2,3,...,N)的每一个点,得到对应的初始匹配点PixWi(i=1,2,3,...,N)。
4.选取初始匹配同名点对{PixI1i,PixWi}(i=1,2,3,...,N)的一个子集{PixI1j,PixWj}(j=1,2,3,...,M<N),本发明中M=50。基于金字塔和相关系数进行特征点匹配(匹配搜索范围优于HJ影像***几何校正的精度),一个线程一次匹配一个点,达到并行计算目的,经粗差剔除后获得更准确的同名点对,统计求得待配准图像经纬度方向整体偏移量:
(1)将基准图裁剪结果影像I1 *和待配准裁剪图像的尺度变换图像WD建立3层高斯金字塔,图像金字塔从下到上对应I1 *、I11 *、I12 *和WD、WD1、WD2,从下往上每层图像缩小1/3。
(2)第1层:
(2.1)提取初始匹配同名点子集{PixI1j,PixWj}(j=1,2,3,...,M)的一个初始匹配同名点{PixI1j,PixWj}。
(2.2)在顶层影像I12 *上以PixI1j/9为中心获取127*127图像块,WD2上以PixWj/9为中心在周围5个像素范围滑窗,进行相关匹配,获取相关系数矩阵C及其峰值对应位置(p,q)和峰值处相关系数Cmax,当Cmax>0.5时,进行下一步;否则,回到步骤(2.1)计算下1个匹配点。
相关系数匹配表述为:它是基于灰度相关匹配算法中最常用的算法,经过去均值和归一化处理后的相关系数对图像对比度、明暗变得不敏感,使得归一化互相关系数更可靠、适应性更强。归一化互相关系数的定义如下:
由归一化互相关的定义可知,相关系数C(A,B)最大值时对应的位置即为两幅图像的匹配位置。
(2.3)在WD2上以(p,q)为中心获取127*127图像块,顶层影像I12 *上以PixI1j/9为中心在其周围滑窗,进行相关匹配,获取相关系数矩阵及其峰值对应位置(m,n),当(m,n)与PixI1j/9对应的点位置相等时,认为PixI1j/9与(p,q)为1对匹配同名点;否则,回到步骤(2.1)计算下1个匹配点。
(2.4)在顶层影像上遍历匹配同名点子集{PixI1j,PixWj}(j=1,2,3,...,M)的每一个点,进行RANSAC(Random Sample Consensus)剔除,获取可靠匹配同名点对{LvI2t,LvW2t}(t=1,2,...,T)。
RANSAC剔除表述为:迭代地在输入数据{PixI1j,PixWj}(j=1,2,3,...,M)中采样所谓的最小点集(3对匹配同名点)。并利用每次采样所得到的最小点集(3对匹配同名点)估计出所要确定的参数(仿射变换参数),同时根据阈值则来判别输入数据中哪些是与该参数相一致的,即内点,哪些是不一致的,即外点。如此迭代一定次数后,将对应输入数据中内点比例最高的所估计参数值以及所筛选出来的内点作为RANSAC最后解。将此解作为其它方法的初始值进一步优化计算,从而得到最后的估计参数值。抽样次数N的计算公式如下:
其中σ为预期的原始数据错误率,s代表至少有一个最小子集包含所有内点的概率。
(3)第2层:
(3.1)根据顶层第1层匹配同名点{LvI2t,LvW2t}(t=1,2,...,T),各点位置乘以金子塔放大倍数3,计算第2层待配准图像到基准图像的仿射变换关系;
(3.2)提取初始匹配同名点子集{PixI1j,PixWj}(j=1,2,3,...,M)的一个初始匹配同名点{PixI1j,PixWj}。
(3.3)在第2层影像I11 *上,根据放射关系计算PixI1j/3在WD1上对应的位置PixW* j,计算WD1上PixWj/3与PixW* j在x和y方向上的距离,若均小于15则进行下一步;否则,回到步骤(3.2)计算下1个匹配点。
(3.4)在第2层影像I11 *上,以PixI1j/3为中心获取127*127图像块,WD1上以PixW* j为中心在周围5个像素范围滑窗,进行相关匹配,获取相关系数矩阵C及其峰值对应位置(p,q)。
(3.5)在WD1上以(p,q)为中心获取127*127图像块,第2层影像I11 *上,以PixI1j/3为中心在其周围滑窗,进行相关匹配,获取相关系数矩阵及其峰值对应位置(m,n),当(m,n)与PixI1j/3对应的点位置相等时,认为PixI1j/3与(p,q)为1对匹配同名点;否则,回到步骤(3.2)计算下1个匹配点。
(3.6)在第二影像上遍历匹配同名点子集{PixI1j,PixWj}(j=1,2,3,...,M)的每一个点,获取匹配同名点对{LvI1t,LvW1t}(t=1,2,...,T)。
(4)第3层(最底层):
(4.1)根据顶层第2层匹配同名点{LvI1t,LvW1t}(t=1,2,...,T),各点位置乘以金子塔放大倍数3,计算第3层待配准图像到基准图像的仿射变换关系;
(4.2)提取初始匹配同名点子集{PixI1j,PixWj}(j=1,2,3,...,M)的一个初始匹配同名点{PixI1j,PixWj}。
(4.3)在第3层影像I1 *上,根据放射关系计算PixI1j在WD上对应的位置PixW* j,计算WD1上PixWj与PixW* j在x和y方向上的距离,若均小于15则进行下一步;否则,回到步骤(4.2)计算下1个匹配点。
(4.4)在第3层影像I1 *上,以PixI1j中心获取127*127图像块,WD上以PixW* j为中心在周围5个像素范围滑窗,进行相关匹配,获取相关系数矩阵C及其峰值对应位置(p,q)。
(4.5)在WD上以(p,q)为中心获取127*127图像块,在第3层影像I1 *上,以PixI1j为中心在其周围滑窗,进行相关匹配,获取相关系数矩阵及其峰值对应位置(m,n),当(m,n)与PixI1j对应的点位置相等时,认为PixI1j与(P,q)为1对匹配同名点;否则,回到步骤(3.2)计算下1个匹配点。
(4.6)在第3影像上遍历匹配同名点子集{PixI1j,PixWj}(j=1,2,3,...,M)的每一个点,获取匹配同名点对{LvIt,LvWt}(t=1,2,...,T)。
(5)统计经纬度方向整体偏移量Δx,Δy:由匹配过程可知,最底层匹配同名点{LvIt,LvWt}(t=1,2,...,T)中的LvIt(t=1,2,...,T)为始匹配同名点子集{PixI1j,PixWj}(j=1,2,3,...,M)中PixI1j的一个子集,可得到LvIt(t=1,2,...,T)对应的初始匹配点PixWt(t=1,2,...,T),计算LvWt(t=1,2,...,T)与PixWt(t=1,2,...,T)在经纬度方向(对应图像横轴x和纵轴y方向)的均值偏移量Δx和Δy。
5.根据经纬度方向整体偏移量Δx,Δy修正初始匹配同名点对{PixI1i,PixWi}(i=1,2,3,...,N)中待配准图对应匹配点PixWi(i=1,2,3,...,N),其x坐标减去Δx,y坐标减去Δy;得到修正后匹配同名点对{PixI1i,PixW* i}(i=1,2,3,...,N)。
6.对进行双向互相关特征点匹配,匹配峰值处采用二元三点插值获取亚像素定位精度的精匹配同名点;同时对精匹配同名点对进行多项式拟合验证,剔除显著误匹配点:
(1)提取修正后匹配同名点对{PixI1i,PixW* i}(i=1,2,3,...,N)的一个匹配同名点{PixI1i,PixW* i}
(2)在基准图裁剪结果影像I1 *上以PixI1j为中心获取127*127图像块,待配准裁剪图像的尺度变换图像WD上以PixW* i为中心在周围5个像素范围滑窗,进行相关匹配,获取相关系数矩阵C及其峰值对应位置(p,q)和峰值处相关系数Cmax,当Cmax>0.5时,进行下一步;否则,回到步骤(1)计算下1个匹配点。
(3)在WD上以(p,q)为中心获取127*127图像块,基准图裁剪结果影像I1 *上以PixI1j/9为中心在其周围滑窗,进行相关匹配,获取相关系数矩阵及其峰值对应位置(m,n),当(m,n)与PixI1j对应的点位置相等时,认为PixI1j与(p,q)为1对匹配同名点,进行下一步;否则,回到步骤(1)计算下1个匹配点。
(4)对相关系数矩阵C峰值对应位置(p,q)周围3*3邻域内采用二元三点插值算法对相关系数进行插值,插值间隔选取1/20,获取插值后3*3邻域相关系数矩阵Z,获取Z的对应的峰值位置(p*,q*),此时PixI1j与(P*,q*)为一对高精度匹配同名点。
(5)回到步骤(1),遍历修正后匹配同名点对{PixI1i,PixW* i}(i=1,2,3,...,N)的每一个同名点,得到亚像素匹配同名点对{PixID1i,PixWDi}(i=1,2,3,...,M)。
(6)对精匹配同名点对{PixID1i,PixWDi}(i=1,2,3,...,M)进行RANSAC剔除,获取可靠精匹配同名点对{PixID1i,PixWDi}(i=1,2,3,...,D)。
7.结合获取的大量精匹配同名点对{PixID1i,PixWDi}(i=1,2,3,...,D),构建密集Delaunay三角网,对各小面元分别进行纠正,得到高精度配准图像。
本发明未详细阐述的部分属于本领域公知技术。
Claims (5)
1.一种多源遥感图像自动配准方法,其特征在于该方法包含以下步骤:
A.对基准图像和待配准图像的经纬度重叠区进行裁剪,若基准图裁剪结果影像和待配准裁剪结果影像的分辨率不一致,则对待配准裁剪结果影像进行Mallat小波尺度变换;
B.对基准图裁剪结果影像进行FAST特征点检测,提取均匀分布的特征点;依据基准图裁剪结果影像和待配准裁剪结果影像的地理信息对应关系在待配准裁剪结果影像上定位特征点的匹配点,得到初始匹配同名点对;
C.提取初始匹配同名点对的一个子集,在子集内基于金字塔和相关系数进行特征点匹配,剔除显著错误点,统计计算得到经纬度方向整体偏移量;
D.通过经纬度方向整体偏移量对初始匹配同名点对进行修正,剔除平移误差;对每个匹配同名点进行双向互相关匹配并在峰值处进行二元三点插值,获取亚像素级匹配精度的匹配同名点对,利用亚像素级匹配精度的匹配同名点对获得配准图像。
2.如权利要求1中所述的一种多源遥感图像自动配准方法,其特征在于:在所述步骤A中所述的对待配准裁剪结果影像进行Mallat小波尺度变换的方法为:
A1.对原始待配准图裁剪结果影像进行Mallat小波分解,得到一个低频图像信号以及不同分辨率上的高频图像信号;
A2.根据其中同一方向的各高频图像信号之间的相似性,进行高频外推插值;
A3.对插值后的低频图像信号以及不同分辨率上的高频图像信号利用小波反变换重构出比原始待配准图裁剪结果影像分辨率更高的待配准图裁剪结果影像。
3.如权利要求1中所述的一种多源遥感图像自动配准方法,其特征在于:在步骤B中获得初始匹配同名点对的步骤中:
利用GDAL遥感图像库的OpenGIS坐标转换方法以及坐标***相互之间转换的服务,通过地理定位信息确定待配准裁剪图像的尺度变换图像上对应的初始匹配同名点。
4.如权利要求1中所述的一种多源遥感图像自动配准方法,其特征在于:在所述步骤C中获得经纬度方向整体偏移量的步骤中:
对基准图像和待配准图像裁剪结果图像建立3层高斯金字塔,初始匹配同名点的子集进行变换后先在金子塔顶层进行双向互相关匹配,匹配关系逐步映射到下层金字塔。
5.如权利要求1中所述的一种多源遥感图像自动配准方法,其特征在于:在所述步骤D中获取亚像素级匹配精度的匹配同名点对的步骤为:
D1.通过经纬度方向整体偏移量对初始匹配同名点对进行修正,剔除平移误差后,对每个匹配同名点进行相关匹配;
D2.相关系数峰值处采用采用二元三点插值获取亚像素定位精度。所有修正后匹配同名点遍历完后,对精匹配同名点对进行RANSAC剔除显著误匹配点,获取匹配同名点对。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110254958.XA CN102436652B (zh) | 2011-08-31 | 2011-08-31 | 一种多源遥感图像自动配准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110254958.XA CN102436652B (zh) | 2011-08-31 | 2011-08-31 | 一种多源遥感图像自动配准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102436652A true CN102436652A (zh) | 2012-05-02 |
CN102436652B CN102436652B (zh) | 2014-08-27 |
Family
ID=45984695
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110254958.XA Active CN102436652B (zh) | 2011-08-31 | 2011-08-31 | 一种多源遥感图像自动配准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102436652B (zh) |
Cited By (33)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102750668A (zh) * | 2012-07-04 | 2012-10-24 | 西南交通大学 | 结合局部方向特征的数字图像三倍插值放大方法 |
CN103115614A (zh) * | 2013-01-21 | 2013-05-22 | 武汉大学 | 多源多轨长条带卫星遥感影像联合并行匹配方法 |
CN103337052A (zh) * | 2013-04-17 | 2013-10-02 | 国家测绘地理信息局卫星测绘应用中心 | 面向宽幅遥感影像的自动几何纠正方法 |
CN103679675A (zh) * | 2013-11-29 | 2014-03-26 | 航天恒星科技有限公司 | 一种面向水质定量遥感应用的遥感影像融合方法 |
CN103905746A (zh) * | 2012-12-28 | 2014-07-02 | 清华大学 | 亚像素级图像偏移定位及叠加方法和装置以及视频设备 |
CN104021556A (zh) * | 2014-06-13 | 2014-09-03 | 西南交通大学 | 一种基于几何结构相似性的异源遥感影像配准方法 |
CN104809724A (zh) * | 2015-04-21 | 2015-07-29 | 电子科技大学 | 多波段遥感影像的自动精配准方法 |
CN104966283A (zh) * | 2015-05-22 | 2015-10-07 | 北京邮电大学 | 图像分层配准方法 |
CN105136164A (zh) * | 2015-08-13 | 2015-12-09 | 航天恒星科技有限公司 | 考虑卫星综合运动的凝视成像仿真及质量评估方法及装置 |
CN105160624A (zh) * | 2015-08-20 | 2015-12-16 | 中电科海洋信息技术研究院有限公司 | 地理信息图片自动配准方法及装置 |
CN105205812A (zh) * | 2015-09-01 | 2015-12-30 | 哈尔滨工业大学 | 一种基于微小卫星星座的多帧图像重构方法 |
CN105913435A (zh) * | 2016-04-13 | 2016-08-31 | 西安航天天绘数据技术有限公司 | 一种适用于大区域的多尺度遥感影像匹配方法及*** |
CN106056625A (zh) * | 2016-05-25 | 2016-10-26 | 中国民航大学 | 一种基于地理同名点配准的机载红外运动目标检测方法 |
CN106228593A (zh) * | 2015-05-28 | 2016-12-14 | 长沙维纳斯克信息技术有限公司 | 一种图像密集匹配方法 |
CN103823889B (zh) * | 2014-03-10 | 2017-02-01 | 北京大学 | 基于l1范数全局几何一致性检验的错误匹配检测方法 |
CN106447613A (zh) * | 2016-09-27 | 2017-02-22 | 西安蒜泥电子科技有限责任公司 | 一种基于图像局部配准的全景去虚影方法及*** |
CN107016695A (zh) * | 2017-04-13 | 2017-08-04 | 首都师范大学 | 一种亚像素级影像配准方法及*** |
CN107609183A (zh) * | 2017-09-29 | 2018-01-19 | 浙江科澜信息技术有限公司 | 将原始坐标数据转换为球面坐标数据的方法、装置及设备 |
CN107705244A (zh) * | 2017-09-11 | 2018-02-16 | 中国国土资源航空物探遥感中心 | 一种适用于大区域多幅遥感影像的接边纠正方法 |
CN107945216A (zh) * | 2017-11-10 | 2018-04-20 | 西安电子科技大学 | 基于最小二乘估计的多图像联合配准方法 |
CN108629798A (zh) * | 2018-04-28 | 2018-10-09 | 安徽大学 | 基于gpu的图像快速配准方法 |
CN109144095A (zh) * | 2018-04-03 | 2019-01-04 | 奥瞳***科技有限公司 | 用于无人驾驶飞行器的基于嵌入式立体视觉的障碍躲避*** |
CN109493298A (zh) * | 2018-11-13 | 2019-03-19 | 中国国土资源航空物探遥感中心 | 一种机载摆扫式高光谱数据快速几何校正方法 |
CN110009670A (zh) * | 2019-03-28 | 2019-07-12 | 上海交通大学 | 基于fast特征提取和piifd特征描述的异源图像配准方法 |
CN110163896A (zh) * | 2019-03-29 | 2019-08-23 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | Sar图像匹配方法 |
CN110310308A (zh) * | 2019-06-18 | 2019-10-08 | 中南林业科技大学 | 一种基于子图的图像配准方法 |
CN110543886A (zh) * | 2018-05-28 | 2019-12-06 | 天津理工大学 | 一种基于月面ccd影像的金字塔成像匹配算法 |
CN111062976A (zh) * | 2019-12-25 | 2020-04-24 | 中国科学院长春光学精密机械与物理研究所 | 基于fmt的低轨卫星太阳望远镜遥感图像配准方法 |
CN111242006A (zh) * | 2020-01-10 | 2020-06-05 | 长江水利委员会长江科学院 | 一种将基于Mask R-CNN的遥感影像地物检测实现为地理WPS服务的方法 |
CN113223065A (zh) * | 2021-03-30 | 2021-08-06 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | Sar卫星图像与光学图像自动匹配方法 |
CN113344937A (zh) * | 2021-04-25 | 2021-09-03 | 中国科学院空天信息创新研究院 | 一种自动识别和去除遥感影像黑边的方法 |
CN113781529A (zh) * | 2021-09-24 | 2021-12-10 | 中国科学院精密测量科学与技术创新研究院 | 一种两次分块策略的广域sar复图像序列快速配准方法 |
CN115265424A (zh) * | 2022-09-27 | 2022-11-01 | 威海晶合数字矿山技术有限公司 | 一种基于合成孔径雷达技术的地质灾害边坡位移监测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6681057B1 (en) * | 2000-02-22 | 2004-01-20 | National Instruments Corporation | Image registration system and method implementing PID control techniques |
CN102117483A (zh) * | 2009-12-31 | 2011-07-06 | 核工业北京地质研究院 | 不同空间分辨率的多光谱遥感图像融合方法 |
-
2011
- 2011-08-31 CN CN201110254958.XA patent/CN102436652B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6681057B1 (en) * | 2000-02-22 | 2004-01-20 | National Instruments Corporation | Image registration system and method implementing PID control techniques |
CN102117483A (zh) * | 2009-12-31 | 2011-07-06 | 核工业北京地质研究院 | 不同空间分辨率的多光谱遥感图像融合方法 |
Non-Patent Citations (4)
Title |
---|
屈有山: "基于小波双三次插值提高光学遥感图像空间分辨率的研究", 《光子学报》, vol. 33, no. 5, 31 May 2004 (2004-05-31), pages 601 - 604 * |
张继贤: "多源遥感影像高精度自动配准的方法研究", 《遥感学报》, vol. 9, no. 1, 31 January 2005 (2005-01-31), pages 73 - 77 * |
成英燕: "用InSAR技术进行形变监测的研究", 《测绘科学》, vol. 31, no. 3, 31 May 2006 (2006-05-31), pages 56 - 59 * |
蔡志凌: "MODIS近红外遥感反演水汽含量及水汽含量对气象数据的回归分析", 《中国科技论文在线》, 12 May 2008 (2008-05-12), pages 1 - 6 * |
Cited By (48)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102750668B (zh) * | 2012-07-04 | 2014-07-02 | 西南交通大学 | 结合局部方向特征的数字图像三倍插值放大方法 |
CN102750668A (zh) * | 2012-07-04 | 2012-10-24 | 西南交通大学 | 结合局部方向特征的数字图像三倍插值放大方法 |
CN103905746A (zh) * | 2012-12-28 | 2014-07-02 | 清华大学 | 亚像素级图像偏移定位及叠加方法和装置以及视频设备 |
CN103905746B (zh) * | 2012-12-28 | 2017-02-22 | 清华大学 | 亚像素级图像偏移定位及叠加方法和装置以及视频设备 |
CN103115614A (zh) * | 2013-01-21 | 2013-05-22 | 武汉大学 | 多源多轨长条带卫星遥感影像联合并行匹配方法 |
CN103115614B (zh) * | 2013-01-21 | 2014-12-17 | 武汉大学 | 多源多轨长条带卫星遥感影像联合并行匹配方法 |
CN103337052A (zh) * | 2013-04-17 | 2013-10-02 | 国家测绘地理信息局卫星测绘应用中心 | 面向宽幅遥感影像的自动几何纠正方法 |
CN103337052B (zh) * | 2013-04-17 | 2016-07-13 | 国家测绘地理信息局卫星测绘应用中心 | 面向宽幅遥感影像的自动几何纠正方法 |
CN103679675A (zh) * | 2013-11-29 | 2014-03-26 | 航天恒星科技有限公司 | 一种面向水质定量遥感应用的遥感影像融合方法 |
CN103679675B (zh) * | 2013-11-29 | 2017-01-11 | 航天恒星科技有限公司 | 一种面向水质定量遥感应用的遥感影像融合方法 |
CN103823889B (zh) * | 2014-03-10 | 2017-02-01 | 北京大学 | 基于l1范数全局几何一致性检验的错误匹配检测方法 |
CN104021556A (zh) * | 2014-06-13 | 2014-09-03 | 西南交通大学 | 一种基于几何结构相似性的异源遥感影像配准方法 |
CN104809724A (zh) * | 2015-04-21 | 2015-07-29 | 电子科技大学 | 多波段遥感影像的自动精配准方法 |
CN104966283A (zh) * | 2015-05-22 | 2015-10-07 | 北京邮电大学 | 图像分层配准方法 |
CN106228593B (zh) * | 2015-05-28 | 2019-05-17 | 长沙维纳斯克信息技术有限公司 | 一种图像密集匹配方法 |
CN106228593A (zh) * | 2015-05-28 | 2016-12-14 | 长沙维纳斯克信息技术有限公司 | 一种图像密集匹配方法 |
CN105136164B (zh) * | 2015-08-13 | 2019-04-05 | 航天恒星科技有限公司 | 考虑卫星综合运动的凝视成像仿真及质量评估方法及装置 |
CN105136164A (zh) * | 2015-08-13 | 2015-12-09 | 航天恒星科技有限公司 | 考虑卫星综合运动的凝视成像仿真及质量评估方法及装置 |
CN105160624A (zh) * | 2015-08-20 | 2015-12-16 | 中电科海洋信息技术研究院有限公司 | 地理信息图片自动配准方法及装置 |
CN105205812A (zh) * | 2015-09-01 | 2015-12-30 | 哈尔滨工业大学 | 一种基于微小卫星星座的多帧图像重构方法 |
CN105913435A (zh) * | 2016-04-13 | 2016-08-31 | 西安航天天绘数据技术有限公司 | 一种适用于大区域的多尺度遥感影像匹配方法及*** |
CN105913435B (zh) * | 2016-04-13 | 2019-05-28 | 西安航天天绘数据技术有限公司 | 一种适用于大区域的多尺度遥感影像匹配方法及*** |
CN106056625B (zh) * | 2016-05-25 | 2018-11-27 | 中国民航大学 | 一种基于地理同名点配准的机载红外运动目标检测方法 |
CN106056625A (zh) * | 2016-05-25 | 2016-10-26 | 中国民航大学 | 一种基于地理同名点配准的机载红外运动目标检测方法 |
CN106447613A (zh) * | 2016-09-27 | 2017-02-22 | 西安蒜泥电子科技有限责任公司 | 一种基于图像局部配准的全景去虚影方法及*** |
CN107016695A (zh) * | 2017-04-13 | 2017-08-04 | 首都师范大学 | 一种亚像素级影像配准方法及*** |
CN107016695B (zh) * | 2017-04-13 | 2019-09-17 | 首都师范大学 | 一种亚像素级影像配准方法及*** |
CN107705244A (zh) * | 2017-09-11 | 2018-02-16 | 中国国土资源航空物探遥感中心 | 一种适用于大区域多幅遥感影像的接边纠正方法 |
CN107609183A (zh) * | 2017-09-29 | 2018-01-19 | 浙江科澜信息技术有限公司 | 将原始坐标数据转换为球面坐标数据的方法、装置及设备 |
CN107945216A (zh) * | 2017-11-10 | 2018-04-20 | 西安电子科技大学 | 基于最小二乘估计的多图像联合配准方法 |
CN107945216B (zh) * | 2017-11-10 | 2019-10-11 | 西安电子科技大学 | 基于最小二乘估计的多图像联合配准方法 |
CN109144095A (zh) * | 2018-04-03 | 2019-01-04 | 奥瞳***科技有限公司 | 用于无人驾驶飞行器的基于嵌入式立体视觉的障碍躲避*** |
CN108629798A (zh) * | 2018-04-28 | 2018-10-09 | 安徽大学 | 基于gpu的图像快速配准方法 |
CN110543886A (zh) * | 2018-05-28 | 2019-12-06 | 天津理工大学 | 一种基于月面ccd影像的金字塔成像匹配算法 |
CN109493298A (zh) * | 2018-11-13 | 2019-03-19 | 中国国土资源航空物探遥感中心 | 一种机载摆扫式高光谱数据快速几何校正方法 |
CN110009670A (zh) * | 2019-03-28 | 2019-07-12 | 上海交通大学 | 基于fast特征提取和piifd特征描述的异源图像配准方法 |
CN110163896B (zh) * | 2019-03-29 | 2023-02-03 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | Sar图像匹配方法 |
CN110163896A (zh) * | 2019-03-29 | 2019-08-23 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | Sar图像匹配方法 |
CN110310308A (zh) * | 2019-06-18 | 2019-10-08 | 中南林业科技大学 | 一种基于子图的图像配准方法 |
CN111062976A (zh) * | 2019-12-25 | 2020-04-24 | 中国科学院长春光学精密机械与物理研究所 | 基于fmt的低轨卫星太阳望远镜遥感图像配准方法 |
CN111062976B (zh) * | 2019-12-25 | 2023-02-28 | 中国科学院长春光学精密机械与物理研究所 | 基于fmt的低轨卫星太阳望远镜遥感图像配准方法 |
CN111242006B (zh) * | 2020-01-10 | 2021-04-09 | 长江水利委员会长江科学院 | 一种将基于Mask R-CNN的遥感影像地物检测实现为地理WPS服务的方法 |
CN111242006A (zh) * | 2020-01-10 | 2020-06-05 | 长江水利委员会长江科学院 | 一种将基于Mask R-CNN的遥感影像地物检测实现为地理WPS服务的方法 |
CN113223065A (zh) * | 2021-03-30 | 2021-08-06 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | Sar卫星图像与光学图像自动匹配方法 |
CN113223065B (zh) * | 2021-03-30 | 2023-02-03 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | Sar卫星图像与光学图像自动匹配方法 |
CN113344937A (zh) * | 2021-04-25 | 2021-09-03 | 中国科学院空天信息创新研究院 | 一种自动识别和去除遥感影像黑边的方法 |
CN113781529A (zh) * | 2021-09-24 | 2021-12-10 | 中国科学院精密测量科学与技术创新研究院 | 一种两次分块策略的广域sar复图像序列快速配准方法 |
CN115265424A (zh) * | 2022-09-27 | 2022-11-01 | 威海晶合数字矿山技术有限公司 | 一种基于合成孔径雷达技术的地质灾害边坡位移监测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102436652B (zh) | 2014-08-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102436652B (zh) | 一种多源遥感图像自动配准方法 | |
US11244197B2 (en) | Fast and robust multimodal remote sensing image matching method and system | |
CN104574347B (zh) | 基于多源遥感数据的在轨卫星图像几何定位精度评价方法 | |
CN103679714B (zh) | 一种基于梯度互相关的光学和sar图像自动配准方法 | |
CN103400151B (zh) | 一体化的光学遥感影像与gis自动配准与水体提取方法 | |
CN102567979B (zh) | 车载红外夜视***及其多源图像融合方法 | |
CN107610164B (zh) | 一种基于多特征混合的高分四号影像配准方法 | |
CN104021556A (zh) | 一种基于几何结构相似性的异源遥感影像配准方法 | |
CN103839265A (zh) | 基于sift和归一化互信息的sar图像配准方法 | |
CN105184801A (zh) | 一种基于多层次策略的光学与sar影像高精度配准方法 | |
CN103226826B (zh) | 基于局部熵视觉注意模型的遥感图像变化检测方法 | |
CN111008664B (zh) | 一种基于空谱联合特征的高光谱海冰检测方法 | |
CN103295239A (zh) | 一种基于平面基准影像的激光点云数据的自动配准方法 | |
CN103914678A (zh) | 基于纹理与植被指数的撂荒地遥感识别方法 | |
CN105651263A (zh) | 浅海水深多源遥感融合反演方法 | |
CN103839243B (zh) | 基于Shearlet变换的多通道卫星云图融合方法 | |
CN103093459A (zh) | 利用机载LiDAR点云数据辅助影像匹配的方法 | |
CN104809724A (zh) | 多波段遥感影像的自动精配准方法 | |
CN103489178A (zh) | 一种图像配准方法和*** | |
CN107688777A (zh) | 一种协同多源遥感影像的城市绿地提取方法 | |
CN103218811A (zh) | 一种基于统计分布的卫星多光谱图像波段配准方法 | |
CN104318583A (zh) | 一种可见光宽带光谱图像配准方法 | |
CN107688776A (zh) | 一种城市水体提取方法 | |
CN116878748A (zh) | 一种激光与图像融合的智能气体泄漏定位方法及装置 | |
Hou et al. | A pointer meter reading recognition method based on YOLOX and semantic segmentation technology |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |