CN103247029A - 一种用于拼接式探测器生成的高光谱图像几何配准方法 - Google Patents
一种用于拼接式探测器生成的高光谱图像几何配准方法 Download PDFInfo
- Publication number
- CN103247029A CN103247029A CN2013100983889A CN201310098388A CN103247029A CN 103247029 A CN103247029 A CN 103247029A CN 2013100983889 A CN2013100983889 A CN 2013100983889A CN 201310098388 A CN201310098388 A CN 201310098388A CN 103247029 A CN103247029 A CN 103247029A
- Authority
- CN
- China
- Prior art keywords
- image
- gray
- value
- point
- geometric
- 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 25
- VMXUWOKSQNHOCA-UKTHLTGXSA-N ranitidine Chemical compound [O-][N+](=O)\C=C(/NC)NCCSCC1=CC=C(CN(C)C)O1 VMXUWOKSQNHOCA-UKTHLTGXSA-N 0.000 claims description 14
- 238000001228 spectrum Methods 0.000 claims description 9
- 238000012417 linear regression Methods 0.000 claims description 7
- 238000012952 Resampling Methods 0.000 claims description 4
- 238000004458 analytical method Methods 0.000 claims description 3
- 238000012216 screening Methods 0.000 claims description 3
- 230000004069 differentiation Effects 0.000 claims description 2
- 230000007704 transition Effects 0.000 claims description 2
- 239000004744 fabric Substances 0.000 description 5
- 238000013461 design Methods 0.000 description 4
- 230000000875 corresponding effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- DGJPPCSCQOIWCP-UHFFFAOYSA-N cadmium mercury Chemical compound [Cd].[Hg] DGJPPCSCQOIWCP-UHFFFAOYSA-N 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Images
Landscapes
- Image Processing (AREA)
Abstract
本发明公开了一种用于拼接式探测器生成的高光谱图像几何配准方法,它应用于品字型排列探测器,其排列方式导致的图像几何错位情况。其方法首先寻找高光谱图像上的狭长地物目标,通过求算质心位置进行目标地物的边缘提取,分别对错位两边图像中目标地物的边缘做线性拟合,比较错位两边图像的拟合偏置,从而得到错位图像的亚象元个数。根据结果对错位图像进行重采样实现几何校正。通过本发明的几何错位配准方法能够高精度探测错位值且不引入其他误差,是该领域的一项关键突破,方法可靠且有实用价值。同时可为同类型探测器的遥感图像几何配准提供参考依据。
Description
技术领域
本发明涉及图像处理领域,具体是指一种用于拼接式探测器生成的高光谱图像的几何配准方法,它用于校正图像中的几何错位现象。
背景技术
为了获得较大视场地物图像,需要增加探测器空间维度。受仪器地面分辨率、幅宽等指标需求,以及红外碲镉汞探测器设计工艺的限制,新研制的高光谱成像仪在短波红外(1000nm-2500nm)范围内,焦平面排列采用四块探测器阵列分列两行错开拼接的设计方案,并采用两个光路***进行逐行的探测器阵列成像,最终将四块探测器阵列所形成的像进行拼接,见图1。然而由于探测器阵列焦平面类似于多元排列,并在横向有一定的错位,表现在图像上为相邻探测器阵列扫描带的错位现象。因此成像后需进行多视场配准校正,以确保图像的几何定位精度能满足设计和应用要求。
传统的几何配准方法大致分为基于频率域、基于灰度相关以及基于特征的方法。基于频率域的方法是利用傅里叶变换将两幅待匹配的图像变换到频域,然后利用互功率谱直接计算出两幅图像间的平移矢量,但图像的多地物多信息性使得计算得到的误差较大;基于灰度相关的方法是寻找两幅图像重叠部分,利用它们所对应的颜色***中灰度级的相似性为准则寻找图像的配准位置,这需要两幅图像之间具有较大面积的重叠区域;比较常用的方法是基于特征的方法,具体算法为从待配准图像中提取特征集,利用特征的对应关系进行匹配。然而,本研究的高光谱成像仪不同探测器阵列之间的重叠区域小于10个像元, 特征难提取,配准精度低,不适用于使用基于灰度相关或者基于特征的方法,会造成很大的匹配误差。因此考虑通过计算错位图像中线性地物的线性拟合偏置差得到错位的亚像元个数,根据求得结果对错位图像进行重采样的方法来完成几何配准。
发明内容
基于上述已有技术存在的问题,本发明的目的是提出一种几何空间图像错位的配准方法,以修正由于光机***和探测器排列限制而导致的图像空间非匹配问题。
该匹配方法首先寻找特殊目标地物,然后根据质心法进行统计数值计算得到地物边界信息,进行边缘提取后对待配准图像进行边缘线性拟合,从而计算两条拟合直线的偏置差,即为偏移的亚像元值,最后根据偏移量对图像进行三次样条重采样以实现校正。
其具体步骤是:
1)筛选高光谱图像。对高光谱成像仪采集到的高光谱图像进行人工目视判读筛选,挑选河流、桥梁、公路、麦田边缘或海岸线这些狭长地物目标,所选图像应是覆盖不同探测器拼接的图像。
2)对地物目标进行边缘提取。具体方法为:
a)记录探测器的固定拼接位置和狭长目标地物的交界点为分割点(x0,y0),以分割点(x0,y0)为基准分别选择错位图像左右两边狭长地物的纵坐标值yj,j=0,1,2...m,选择纵坐标点个数>8个,以保证线性拟合的精度。
b)计算横坐标值x’。假设拼接图像两边灰度值均值分别为对于步骤a)中每个yj,在边界区域并非严格的区分,会存在过渡像元且灰度值介于和之间。定义检测到的边缘为像元被分割为两个灰度值面积一致的部分,寻找边界点即求像元灰度值为时的横坐标位置x’,见图2。设采样点为正方形像元,从左向右水平移动(yj,j=0,1,2...m,xi,i=0,1,2...n),选 择灰度值最接近1/2的像元,设其灰度值为DNij,计算其权重Wij分布为:
x’=xi-(1-Wij) (7)
式中:分别表示边界左右两边的灰度值均值;xi表示yj对应的横坐标值;DNij为x’表示边缘提取出的边界点的横坐标值。
3)对特征边缘进行一元线性回归分析,求解表达式。当得到错位图像两边的一组坐标点后,建立一元线性回归模型并求解线性表达式,具体算法为:首先建立一元线性回归模型表达式:
yt=a+b*xt+et (8)
式中:a和b为待定参数;xt和yt分别代表各个横纵坐标点;t=1,2,...,n为各组观测数据的下标;et为随机变量。
取得最小值
根据极值的必要条件,有:
解方程组,求参数a,b的拟合值为:
4)根据拟合结果计算几何错位亚像元个数。根据步骤3)分别得到待配准图像的地物的边缘信息,考虑到错位两边图像的地物是同一地物,因此拟合得到的线性方程斜率相同,两个方程分别为y1=a1+b*x和y2=a2+b*x,其中a1、a2分别代表线性方程偏置,b代表增益,即斜率。则错位的亚像元个数为Δy=a1-a2。
5)对错位图像重采样以实现几何配准。根据4)中结果,Δy表示左边图像相对右边图像错位量,对右边图像进行三次样条插值,最终得到几何配准后的图像。
本发明有如下有益效果:
1)提出了一种新的几何配准方法,特别对于重叠区域小的待匹配图像,显著提高了配准精度。
2)方法的适用性较高,不仅适用于高光谱图像,对于其他类型的图像也适用,有实用价值。
3)本方法先通过边缘提取计算图像偏移量,再根据偏移量重采样图像,在精度允许方位内可以忽略由于计算引入的误差。
附图说明
图1为高光谱成像仪短波红外探测器排列示意图。
图2为求边界点横坐标示意图。
图3仿真的原始待配准图像。
图4发明流程图。
具体实施方式
根据本发明对一幅空间维错位的图像进行几何配准。考虑用仿真图像可以更精确地考察此方法的配准精度,根据图1的探测器排列方式仿真采样后的效果图像。其中设计图像大小为单波段800*1024,其中在空间维512处位置图像左右存在错位,分界两边的灰度值分别为255和0,效果图见图3。
详述如下:
1)确定图像中错位位置为(511,562),在图像的左右区域设置不同的待计算点的纵坐标y,其中左区域确定的纵坐标值为506、513、520、527、534、538、543、546、550;右边区域纵坐标为570、574、578、581、585、588、595、602、613。
2)根据公式(1)和(2)计算出以上各个y值对应的边界处的x对应值,结果见下表:
左区域横坐标值x | 左区域纵坐标值y | 右区域横坐标值x | 右区域纵坐标值y |
720.2549 | 506 | 812.028 | 570 |
730.2549 | 513 | 817.349 | 574 |
740.2549 | 520 | 823.1765 | 578 |
750.2549 | 527 | 827.349 | 581 |
760.2549 | 534 | 833.1765 | 585 |
766.1137 | 538 | 837.349 | 588 |
773.1765 | 543 | 847.349 | 595 |
777.349 | 546 | 857.349 | 602 |
783.1765 | 550 | 873.1765 | 613 |
3)根据步骤2检测出的边缘点坐标,分别做一元线性回归分析,分别得到表达式为:
y1=1.42518*x-0.77718
y2=1.42518*x-0.57934
计算Δy=0.77718-0.57934=0.19784,而模拟的原始图像的错位值为0.2个像元,其计算的相对误差为
4)根据步骤3)计算得到的图像错位量对右边图像进行剪裁后重新采样,以得到几何配准后的图像。
Claims (1)
1.一种用于拼接式探测器生成的高光谱图像的几何配准方法,其特征在于包括的如下步骤:
1)筛选高光谱图像,从探测器生成的高光谱图像中挑选河流、桥梁、公路、麦田边缘或海岸线这些狭长地物目标,所选图像应是覆盖不同探测器拼接的图像;
2)对地物目标进行边缘提取,具体方法为:
a)记录探测器的固定拼接位置和狭长目标地物的交界点为分割点(x0,y0),以分割点(x0,y0)为基准分别选择错位图像左右两边狭长地物的纵坐标值yj,j=0,1,2...m,选择纵坐标点个数>8个,以保证线性拟合的精度;
b)计算横坐标值x’,假设拼接图像两边灰度值均值分别为对于步骤a)中每个yj,在边界区域并非严格的区分,会存在过渡像元且灰度值介于和之间,定义检测到的边缘为像元被分割为两个灰度值,,面积一致的部分,寻找边界点即求像元灰度值为时的横坐标位置x’,设采样点为正方形像元,从左向右水平移动(yj,j=0,1,2...m,xi,i=0,1,2...n),选择灰度值最接近的像元,设其灰度值为DNij,计算其权重Wij分布为:
x’=xi-(1-Wij) (2)
3)对特征边缘进行一元线性回归分析,求解表达式,当得到错位图像两边的一组横纵坐标点后,建立一元线性回归模型并求解线性表达式:
yt=a+b*xt+et (3)
式中:a和b为待定参数;xt和yt分别代表各个横纵坐标点;t=1,2,...,n为各组观测数据的下标;et为随机变量;
4)根据拟合结果计算几何错位亚像元个数。根据步骤3)分别得到待配准图像的地物的边缘信息,考虑到错位两边图像的地物是同一地物,因此拟合得到的线性方程斜率相同,两个方程分别为y1=a1+b*x和y2=a2+b*x,其中a1、a2分别代表线性方程偏置,b代表增益,即斜率,则错位的亚像元个数为Δy=a1-a2;
5)对错位图像重采样以实现几何配准,根据步骤4)中结果,Δy表示左边图像相对右边图像错位量,对右边图像进行三次样条插值,最终得到几何配准后的图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310098388.9A CN103247029B (zh) | 2013-03-26 | 2013-03-26 | 一种用于拼接式探测器生成的高光谱图像几何配准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310098388.9A CN103247029B (zh) | 2013-03-26 | 2013-03-26 | 一种用于拼接式探测器生成的高光谱图像几何配准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103247029A true CN103247029A (zh) | 2013-08-14 |
CN103247029B CN103247029B (zh) | 2015-11-25 |
Family
ID=48926534
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310098388.9A Expired - Fee Related CN103247029B (zh) | 2013-03-26 | 2013-03-26 | 一种用于拼接式探测器生成的高光谱图像几何配准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103247029B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103971334A (zh) * | 2014-04-25 | 2014-08-06 | 中国科学院遥感与数字地球研究所 | 高光谱遥感图像校正方法及装置 |
CN105069743A (zh) * | 2015-07-28 | 2015-11-18 | 中国科学院长春光学精密机械与物理研究所 | 探测器拼接实时图像配准的方法 |
CN105427244A (zh) * | 2015-11-03 | 2016-03-23 | 中南大学 | 一种遥感影像的拼接方法和装置 |
CN106062815A (zh) * | 2014-01-10 | 2016-10-26 | 珀金埃尔默细胞科技德国公司 | 使用准周期栅格的图像校正的方法及*** |
CN103971334B (zh) * | 2014-04-25 | 2016-11-30 | 中国科学院遥感与数字地球研究所 | 高光谱遥感图像校正方法及装置 |
CN106204601A (zh) * | 2016-07-15 | 2016-12-07 | 华东师范大学 | 一种基于波段扫描形式的活体高光谱序列图像并行配准方法 |
CN107367716A (zh) * | 2017-07-04 | 2017-11-21 | 武汉大学 | 一种高精度星载sar几何定标方法 |
CN107945113A (zh) * | 2017-11-17 | 2018-04-20 | 北京天睿空间科技股份有限公司 | 局部图像拼接错位的矫正方法 |
CN109215005A (zh) * | 2018-08-29 | 2019-01-15 | 北京空间机电研究所 | 一种品字形探测器扫描相机图像融合误差修正方法 |
CN110473235A (zh) * | 2019-09-12 | 2019-11-19 | 杭州依图医疗技术有限公司 | 一种对齐两侧***图像的方法及装置 |
CN112070663A (zh) * | 2020-07-24 | 2020-12-11 | 北京空间机电研究所 | 一种多探测器拼接型相机辐射定标联合数据处理方法 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107220939A (zh) * | 2016-03-21 | 2017-09-29 | 吴健奇 | 一种基于三角形定标法的岩芯高光谱数据几何畸变校正方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110013841A1 (en) * | 2009-01-21 | 2011-01-20 | The Government Of The United States Of America, As Represented By The Secretary Of The Navy | Boundary Extraction Method |
CN101968885A (zh) * | 2010-09-25 | 2011-02-09 | 西北工业大学 | 一种基于边缘和灰度的遥感图像变化检测方法 |
-
2013
- 2013-03-26 CN CN201310098388.9A patent/CN103247029B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110013841A1 (en) * | 2009-01-21 | 2011-01-20 | The Government Of The United States Of America, As Represented By The Secretary Of The Navy | Boundary Extraction Method |
CN101968885A (zh) * | 2010-09-25 | 2011-02-09 | 西北工业大学 | 一种基于边缘和灰度的遥感图像变化检测方法 |
Non-Patent Citations (3)
Title |
---|
CHANGNO LEE等: "Georegistration of Airborne Hyperspectral Image Data", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》, vol. 39, no. 7, 31 July 2001 (2001-07-31), pages 1347 - 1351 * |
汪金花等: "高分辨率遥感影像线性基元的提取", 《科技情报开发与经济》, vol. 17, no. 22, 31 December 2007 (2007-12-31), pages 160 - 161 * |
韦燕凤等: "基于特征的遥感图像自动配准算法", 《电子学报》, vol. 33, no. 1, 31 January 2005 (2005-01-31), pages 161 - 165 * |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106062815A (zh) * | 2014-01-10 | 2016-10-26 | 珀金埃尔默细胞科技德国公司 | 使用准周期栅格的图像校正的方法及*** |
CN106062815B (zh) * | 2014-01-10 | 2019-03-15 | 珀金埃尔默细胞科技德国公司 | 使用准周期栅格的图像校正的方法及*** |
CN103971334B (zh) * | 2014-04-25 | 2016-11-30 | 中国科学院遥感与数字地球研究所 | 高光谱遥感图像校正方法及装置 |
CN103971334A (zh) * | 2014-04-25 | 2014-08-06 | 中国科学院遥感与数字地球研究所 | 高光谱遥感图像校正方法及装置 |
CN105069743B (zh) * | 2015-07-28 | 2018-06-26 | 中国科学院长春光学精密机械与物理研究所 | 探测器拼接实时图像配准的方法 |
CN105069743A (zh) * | 2015-07-28 | 2015-11-18 | 中国科学院长春光学精密机械与物理研究所 | 探测器拼接实时图像配准的方法 |
CN105427244B (zh) * | 2015-11-03 | 2018-09-04 | 中南大学 | 一种遥感影像的拼接方法和装置 |
CN105427244A (zh) * | 2015-11-03 | 2016-03-23 | 中南大学 | 一种遥感影像的拼接方法和装置 |
CN106204601B (zh) * | 2016-07-15 | 2018-09-28 | 华东师范大学 | 一种基于波段扫描形式的活体高光谱序列图像并行配准方法 |
CN106204601A (zh) * | 2016-07-15 | 2016-12-07 | 华东师范大学 | 一种基于波段扫描形式的活体高光谱序列图像并行配准方法 |
CN107367716A (zh) * | 2017-07-04 | 2017-11-21 | 武汉大学 | 一种高精度星载sar几何定标方法 |
CN107367716B (zh) * | 2017-07-04 | 2019-09-10 | 武汉大学 | 一种高精度星载sar几何定标方法 |
CN107945113A (zh) * | 2017-11-17 | 2018-04-20 | 北京天睿空间科技股份有限公司 | 局部图像拼接错位的矫正方法 |
CN109215005A (zh) * | 2018-08-29 | 2019-01-15 | 北京空间机电研究所 | 一种品字形探测器扫描相机图像融合误差修正方法 |
CN109215005B (zh) * | 2018-08-29 | 2021-10-01 | 北京空间机电研究所 | 一种品字形探测器扫描相机图像融合误差修正方法 |
CN110473235A (zh) * | 2019-09-12 | 2019-11-19 | 杭州依图医疗技术有限公司 | 一种对齐两侧***图像的方法及装置 |
CN112070663A (zh) * | 2020-07-24 | 2020-12-11 | 北京空间机电研究所 | 一种多探测器拼接型相机辐射定标联合数据处理方法 |
CN112070663B (zh) * | 2020-07-24 | 2024-03-29 | 北京空间机电研究所 | 一种多探测器拼接型相机辐射定标联合数据处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103247029B (zh) | 2015-11-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103247029A (zh) | 一种用于拼接式探测器生成的高光谱图像几何配准方法 | |
Misra et al. | An automatic satellite image registration technique based on Harris corner detection and Random Sample Consensus (RANSAC) outlier rejection model | |
Puniach et al. | Application of UAV-based orthomosaics for determination of horizontal displacement caused by underground mining | |
Soliman et al. | Effects of vertical accuracy of digital elevation model (DEM) data on automatic lineaments extraction from shaded DEM | |
CN102968788B (zh) | 一种基于规则格网面元的波段配准方法 | |
CN107909018B (zh) | 一种稳健的多模态遥感影像匹配方法和*** | |
US20180130217A1 (en) | Method and apparatus for performing background image registration | |
CN102855628B (zh) | 多源多时相高分辨率卫星遥感影像自动匹配方法 | |
CN104077760A (zh) | 一种航空摄影测量的快速拼接***及其实现方法 | |
Muñoz et al. | Environmental applications of camera images calibrated by means of the Levenberg–Marquardt method | |
CN103337068A (zh) | 空间关系约束的多子区匹配方法 | |
CN104268880A (zh) | 基于特征和区域匹配相结合的深度信息获取方法 | |
Scaioni et al. | Image-based deformation measurement | |
Mathavan et al. | Detection of pavement cracks using tiled fuzzy Hough transform | |
Chen et al. | Improving image-based characterization of road junctions, widths, and connectivity by leveraging OpenStreetMap vector map | |
Liu et al. | Phase correlation pixel‐to‐pixel image co‐registration based on optical flow and median shift propagation | |
Maiwald et al. | Solving photogrammetric cold cases using AI-based image matching: New potential for monitoring the past with historical aerial images | |
Müller et al. | Robust image registration for fusion | |
CN104166977A (zh) | 一种图像匹配相似性测度方法及其图像匹配方法 | |
CN114241022B (zh) | 一种无人机图像自动配准方法及*** | |
Merkle et al. | A new approach for optical and sar satellite image registration | |
RU2525601C1 (ru) | Комплексная корреляционно-экстремальная навигационная система | |
Chen et al. | A wavelet transform based method for road centerline extraction | |
Meza et al. | A structure-from-motion pipeline for topographic reconstructions using unmanned aerial vehicles and open source software | |
CN111222544B (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20151125 |
|
CF01 | Termination of patent right due to non-payment of annual fee |