首页 / 专利库 / 资料储存系统 / 数据库 / 一种基于地面高温点的宽覆盖相机几何定位方法

一种基于地面高温点的宽覆盖相机几何定位方法

阅读:116发布:2021-04-14

专利汇可以提供一种基于地面高温点的宽覆盖相机几何定位方法专利检索,专利查询,专利分析的服务。并且本 发明 公开一种基于地面高温点的宽 覆盖 相机几何 定位 方法,该方法步骤为:(1)构建地面高温点 数据库 ;(2)对 云 层遮挡的遥感影像进行云层消除处理;(3)根据高温点的 位置 信息,在高 分辨率 高 精度 遥感图像上进行标注;(4)在视场内,采用基于海岸线、河流和山脉等特征的方法,实现热红外和可见谱段高分辨率高精度遥感图像间配准;(5)计算热红外图像中目标点和标记点之间的欧式距离,满足 阈值 要求时识别成功;(6)以识别后的地面高温点集为参考,解算有理函数模型参数矩阵,实现对地目标高精度定位。本发明方法能够更精确地解算模型参数,定位精度高,且有效避免了严格成像模型内外方位元素测量时系统误差带来的影响。,下面是一种基于地面高温点的宽覆盖相机几何定位方法专利的具体信息内容。

1.一种基于地面高温点的宽覆盖相机几何定位方法,其特征在于包括以下步骤:
1)构建地面高温点数据库,在全球范围内统计地面的油井,高温烟囱,核电站或活火山等高温点位置的经纬度和高程信息,并对其进行编号,构建地面高温点数据库;
2)对层遮挡的遥感影像进行云层消除处理,具体步骤为:
2-1)采用同态滤波的方法消除薄云影响;
2-2)对于厚云覆盖的区域,联合利用多源遥感影像进行去云处理,具体步骤如下:
2-2-1)对多源遥感影像,采用概率松弛匹配法进行配准;
2-2-2)采用仿射变换和双线性内插方法对配准图像进行纠正;
2-2-3)采用直方图匹配和灰度加权内插方法对纠正图像进行辐射改正,消除云层影响;
3)根据高温点的位置信息,在高分辨率精度遥感图像上进行标注;
4)在相机视场内,采用基于海岸线、河流和山脉轮廓特征的方法,实现热红外和可见谱段高分辨率高精度遥感图像间配准,具体步骤为:
4-1)采用维纳滤波器分别对可见与热红外图像进行降噪处理;
4-2)通过Canny算子检测图像边缘、轮廓特征;
4-3)基于Hausdorff距离进行轮廓线匹配,实现图像配准
5)计算热红外图像中目标点和标记点之间的欧式距离D,满足阈值要求:D≤3δ时,其中δ为谱段间配准精度,取1个像元,识别成功;
6)以识别后的地面高温点集为参考,构建有理函数模型求解其参数矩阵,实现对地目标高精度定位;具体步骤如下:
6-1)构建有理函数模型:
式中,(Xn,Yn,Zn),(xn,yn)分别为地面点坐标(X,Y,Z),像点坐标(x,y)经平移和缩放后的正则化坐标,取值为[-1,1];Q1(Xn,Yn,Zn)表示如下:
Q1(Xn,Yn,Zn)=a0+a1Zn+a2Yn+a3Xn+a4ZnYn+a5ZnXn+a6YnXn+a7Zn2+a8Yn2+a9Xn2+a10ZnYnXn+a11Zn2Yn+a12Zn2Xn+a13Yn2Zn+a14Yn2Xn+a15ZnXn2+a16YnXn2+a17Zn3+a18Yn3+a19Xn3    (2)其中,将aj(j=0~19)替换为bj(j=0~19)可得Q2(Xn,Yn,Zn),同理,将aj(j=0~19)替换为cj(j=0~19)可得Q3(Xn,Yn,Zn),将aj(j=0~19)替换为dj(j=0~19)可得Q4(Xn,Yn,Zn),取b0,d0值为1;
6-2)由(1)式可得:
可得误差方程为:
V=Bx-l,P                          (4)
式中,
x=[a0 a1 … a19 b1 … b19 c0 c1 … c19 d1 … d19]T为有理函数模型的系数矩阵,l=[Fx0,Fy0]T为利用物方坐标当前值求得的误差向量,P为虚拟观测值Fx,Fy的权矩阵,这里设为单位矩阵;由此得到法方程:
BTPBx=BTPl                           (5)
6-3)采用步骤(5)中识别的参考点,根据最小二乘平差原理求解该函数模型的参数估计值为:
x=(BTPB)-1BTPl                        (6)
式中x即为模型系数矩阵;由于参考点分布不均匀或者模型的过度参数化会引起法方程的病态,导致模型参数求解不稳定,这里采用岭估计的方法修正:
x=(BTPB+kE)-1BTPl                    (7)
式中,x为模型系数矩阵,E为单位矩阵,k值可由L曲线法求得;求得模型参数后利用该有理函数模型实现高精度对地目标定位。

说明书全文

一种基于地面高温点的宽覆盖相机几何定位方法

技术领域

[0001] 本发明属于空间卫星定位技术领域,尤其涉及一种地球静止轨道宽覆盖遥感探测热红外相机几何定位方法。

背景技术

[0002] 随着卫星技术的发展,越来越多的卫星影像数据在商用、民用领域逐渐推广开来。同时,几十米量级的遥感卫星影像对地目标定位精度已难以满足要求。人们对高精度甚至甚高精度定位技术的要求日益强烈。
[0003] 目前,卫星遥感影像的目标定位方法主要包括两种:一是根据卫星定轨测姿系统提供的影像定向参数,采用严格处理模型或者通用几何处理模型,直接将像点坐标转换为对应地面点的物方空间坐标,称为直接对地目标定位。受限于卫星定轨测姿系统所提供观测数据的精度,卫星在发射、在轨时受加速度、冲击、振动、温度、气压等环境因素影响导致的传感器参数变化,国产遥感卫星影像的直接对地目标定位精度与国外的差距还在一个数量级以上。二是利用同一区域的多重覆盖影像,通过影像间的连接点所满足的构像几何关系,联合求解影像定向参数及其对应地面点的物方空间坐标,称为间接对地目标定位。该方法定位精度一方面依赖于地面控制点的精度和数量,另一方面在模型参数强相关时,物点坐标精度与影像定向精度并不一定能完全对应。综上,目前定位方法难以满足高精度定位的要求。

发明内容

[0004] 本发明的目的在于克服现有技术的不足,提供一种基于地面高温点的宽覆盖相机几何定位方法,解决遥感卫星对地目标定位的精度问题。
[0005] 本发明解决其技术问题是采用以下技术方案实现的:
[0006] 1.一种基于地面高温点的宽覆盖相机几何定位方法,包括以下步骤:
[0007] 1)构建地面高温点数据库,在全球范围内统计地面的油井,高温烟囱,核电站或活火山等高温点位置的经纬度和高程信息,并对其进行编号,构建地面高温点数据库;
[0008] 2)对层遮挡的遥感影像进行云层消除处理,具体步骤为:
[0009] 2-1)采用同态滤波的方法消除薄云影响;
[0010] 2-2)对于厚云覆盖的区域,联合利用多源遥感影像进行去云处理,具体步骤如下:
[0011] 2-2-1)对多源遥感影像,采用概率松弛匹配法进行配准;
[0012] 2-2-2)采用仿射变换和双线性内插方法对配准图像进行纠正;
[0013] 2-2-3)采用直方图匹配和灰度加权内插方法对纠正图像进行辐射改正,消除云层影响;
[0014] 3)根据高温点的位置信息,在高分辨率高精度遥感图像上进行标注;
[0015] 4)在相机视场内,采用基于海岸线、河流和山脉等轮廓特征的方法,实现热红外和可见谱段高分辨率高精度遥感图像间配准,具体步骤为:
[0016] 4-1)采用维纳滤波器分别对可见与热红外图像进行降噪处理;
[0017] 4-2)通过Canny算子检测图像边缘、轮廓特征;
[0018] 4-3)基于Hausdorff距离进行轮廓线匹配,实现图像配准
[0019] 5)计算热红外图像中目标点和标记点之间的欧式距离D,满足阈值要求:D≤3δ时,其中δ为谱段间配准精度,取1个像元,识别成功;
[0020] 6)以识别后的地面高温点为参考点,构建有理函数模型并求解其参数矩阵,实现对地目标高精度定位。具体步骤如下:
[0021] 6-1)构建有理函数模型:
[0022]
[0023] 式中,(Xn,Yn,Zn),(xn,yn)分别为地面点坐标(X,Y,Z),像点坐标(x,y)经平移和缩放后的正则化坐标,取值为[-1,1],具体求解方法如下:
[0024]
[0025]
[0026] 式中,X0,Y0,Z0,x0,y0为正则化的平移参数,Xs,Ys,Zs,xs,ys为正则化的比例系数,具体求解方法如下:
[0027]
[0028]
[0029] 式中m为参考点数。
[0030] Q1(Xn,Yn,Zn)表示如下:
[0031] Q1(Xn,Yn,Zn)=a0+a1Zn+a2Yn+a3Xn+a4ZnYn+a5ZnXn+a6YnXn
[0032] +a7Zn2+a8Yn2+a9Xn2+a10ZnYnXn+a11Zn2Yn+a12Zn2Xn+a13Yn2Zn
[0033] +a14Yn2Xn+a15ZnXn2+a16YnXn2+a17Zn3+a18Yn3+a19Xn3  (6)
[0034] 其中,将aj(j=0~19)替换为bj(j=0~19)可得Q2(Xn,Yn,Zn),同理,将aj(j=0~19)替换为cj(j=0~19)可得Q3(Xn,Yn,Zn),将aj(j=0~19)替换为dj(j=0~19)可得Q4(Xn,Yn,Zn),取b0,d0值为1;
[0035] 6-2)由(1)式可得:
[0036]
[0037] 可得误差方程为:
[0038] V=Bx-l,P  (8)
[0039] 式中,T
x=[a0 a1 … a19 b1 … b19 c0 c1 … c19 d1 … d19] 为有理函数模型的系数矩阵,l=[Fx0,Fy0]T为利用物方坐标当前值求得的误差向量,P为虚拟观测值Fx,Fy的权矩阵,这里设为单位矩阵。由此得到法方程:
[0040] BTPBx=BTPl  (9)
[0041] 6-3)采用步骤(5)中识别的参考点,根据最小二乘平差原理求解该函数模型的参数估计值为:
[0042] x=(BTPB)-1BTPl  (10)
[0043] 式中x为函数模型的系数矩阵,由于参考点分布不均匀或者模型的过度参数化会引起法方程的病态,导致模型参数求解不稳定,这里采用岭估计的方法修正:
[0044] x=(BTPB+kE)-1BTPl  (11)
[0045] 式中,x为模型系数矩阵,E为单位矩阵,k值可由L曲线法求得,具体方法为:首先选择一组不同的k值,按照式(11)求解出系数矩阵x;然后以lg‖Bx-l‖k为横坐标,以lg‖x‖k为纵坐标作图,得到以k为参变量的曲线;最后选取曲线上曲率最大的点所对应的k值即为所求的最优岭参数值。经过修正后即可得到准确的模型参数,求得模型参数后利用该有理函数模型实现高精度对地目标定位。
[0046] 本发明的优点和积极效果是:
[0047] 1.本发明构建了高精度的地面高温点数据库,并以其为参考点,精确解算有理函数模型的参数矩阵,有效避免了通过卫星定轨测姿系统及内外方位元素求解严格几何成像模型时带来的系统误差,同时克服了卫星在发射、在轨时受加速度、冲击力、振动、温度、真空等环境因素导致的传感器参数变化的影响,显著提高对地目标的定位精度,全局定位精度可达像元级,靠近参考点的局部定位精度可达1个像元。
[0048] 2.本发明设计合理,具有精度高、使用方便等特点,可以在高分辨率遥感卫星对地目标定位的普通民用领域应用推广。附图说明
[0049] 图1为本发明的总体流程图
[0050] 图2为有云层遮挡的遥感图像。
[0051] 图3为消除云层遮挡的遥感图像。

具体实施方式

[0052] 下面以视场为15°×15°的静止轨道多谱段凝视相机为例,结合附图对本发明的实施方式作进一步详述:
[0053] 本发明是通过构建高精度的地面高温点数据库,并以其为参考点实现对有理函数成像模型参数矩阵的高精度结算,以此实现遥感卫星高精度对地目标定位。主要包括以下步骤:
[0054] 1.构建地面高温点数据库;首先,在全球范围内统计地面高温点信息,如油井,高温烟囱,核电站和活火山等热源点,并对其进行编号。然后,利用高温点位置的经纬度和高程信息构建精确的地面高温点数据库;
[0055] 2.对云层遮挡的遥感影像进行云层消除处理;采用同态滤波的方法消除薄云影响;对于厚云覆盖的区域,联合利用多源遥感影像进行去云处理,具体步骤为:首先,对多源遥感影像,采用概率松弛匹配法进行配准;然后,采用仿射变换和双线性内插方法对配准图像进行纠正;最后,采用直方图匹配和灰度加权内插方法对纠正图像进行辐射改正,消除云层影响,云层消除效果示意图如图2,图3所示。
[0056] 3.根据数据库高温点的位置信息,在高分辨率高精度遥感图像上进行标注,图2,图3中编号各点为地面的油井,核电站等,即高温点;
[0057] 4.在相机视场内,采用基于海岸线、河流和山脉等轮廓特征的方法,实现热红外和可见谱段高分辨率高精度遥感图像间配准;首先,采用维纳滤波器分别对可见与热红外图像进行降噪处理;然后,通过Canny算子检测图像中边缘、轮廓等特征信息;最后采用基于Hausdorff距离的方法利用轮廓线等特征进行图像匹配,实现图像谱段间配准,配准精度为1个像元;
[0058] 5.通过步骤4进行谱段间图像配准、融合后,计算融合图像中目标点和标记点之间的欧式距离D,当满足D小于3倍的配准精度时,认为目标点识别成功;
[0059] 6.采用步骤5获取的地面高温点集,构建有理函数模型并求解其参数矩阵,实现对地目标高精度定位。具体步骤包括:
[0060] 1)构建有理函数模型:
[0061]
[0062] 式中,(Xn,Yn,Zn),(xn,yn)分别为地面点坐标(X,Y,Z),像点坐标(x,y)经平移和缩放后的正则化坐标,取值为[-1,1],具体求解方法如下:
[0063]
[0064]
[0065] 式中,X0,Y0,Z0,x0,y0为正则化的平移参数,Xs,Ys,Zs,xs,ys为正则化的比例系数,具体求解方法如下:
[0066]
[0067]
[0068] 式中m为参考点数。
[0069] Q1(Xn,Yn,Zn)表示如下:
[0070] Q1(Xn,Yn,Zn)=a0+a1Zn+a2Yn+a3Xn+a4ZnYn+a5ZnXn+a6YnXn
[0071] +a7Zn2+a8Yn2+a9Xn2+a10ZnYnXn+a11Zn2Yn+a12Zn2Xn+a13Yn2Zn
[0072] +a14Yn2Xn+a15ZnXn2+a16YnXn2+a17Zn3+a18Yn3+a19Xn3  (6)
[0073] 其中,将aj(j=0~19)替换为bj(j=0~19)可得Q2(Xn,Yn,Zn),同理,将aj(j=0~19)替换为cj(j=0~19)可得Q3(Xn,Yn,Zn),将aj(j=0~19)替换为dj(j=0~19)可得Q4(Xn,Yn,Zn),取b0,d0值为1。
[0074] 2)由(1)式可得:
[0075]
[0076] 可得误差方程为:
[0077] V=Bx-l,P  (8)
[0078] 式中,x=[a0 a1 … a19 b1 … b19 c0 c1 … c19 d1 … d19]T为有理函数模型的系数矩阵,l=[Fx0,Fy0]T为利用物方坐标当前值求得的误差向量,P为虚拟观测值Fx,Fy的权矩阵,这里设为单位矩阵。由此得到法方程:
[0079] BTPBx=BTPl  (9)
[0080] 3)采用步骤(5)中识别的参考点,根据最小二乘平差原理求解该函数模型的参数估计值为:
[0081] x=(BTPB)-1BTPl  (10)
[0082] 式中x为模型系数矩阵。由于参考点分布不均匀或者模型的过度参数化会引起法方程的病态,导致模型参数求解不稳定,这里采用岭估计的方法修正:
[0083] x=(BTPB+kE)-1BTPl  (11)
[0084] 式中,x为模型系数矩阵,E为单位矩阵,k值可由L曲线法求得,具体方法为:首先选择一组不同的k值,按照式(11)求解出未知数x;然后,以lg‖Bx-l‖k为横坐标,以lg‖x‖k为纵坐标作图,得到以k为参变量的曲线;最后选取曲线上曲率最大的点所对应的k值即为所求的最优岭参数值。经过修正后即可得到准确的模型参数,求得模型参数后利用该有理函数模型实现高精度对地目标定位。
[0085] 通过以上方法,至少选取80个地面点即可求解函数模型参数矩阵,适当增加参考点的数量可实现对地目标定位精度达1个像元。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈