一种激光测高数据辅助的卫星线阵影像拼接方法及系统技术方案

本公开实施例公开了一种激光测高数据辅助的卫星线阵影像拼接方法及系统，涉及影像拼接技术领域

【技术实现步骤摘要】
一种激光测高数据辅助的卫星线阵影像拼接方法及系统


[0001]本公开实施例涉及卫星影像拼接
，尤其涉及一种激光测高数据辅助的卫星线阵影像拼接方法及系统
。

技术介绍

[0002]目前，高分辨率线阵推扫式卫星多采用延迟积分电荷耦合器件
(Time Delayed and Integration Charge
‑
Coupled Device,TDI
‑
CCD)
设计，以此获取高质量卫星影像，同时为了增加幅宽，在焦平面上一般会交错排列安装多个电荷耦合器件
(Charge
‑
Coupled Device
，
CCD)。
在实际使用时一般会将多个
CCD
获取的影像进行拼接后提供给用户使用
。
目前影像拼接方式主要有两种：基于像方的影像拼接和基于物方的影像拼接
。
基于像方的影像拼接主要是在影像重叠区进行同名点匹配，以此获得影像间的平移或者仿射关系，进而完成影像间的拼接
。
这种方法严重依赖影像匹配精度，往往拼接精度较低，同时没有考虑卫星平台稳定性，对于影像内部存在的畸变不能很好的消除，影响后续的应用
。
基于物方的影像拼接主要是基于影像严密成像几何模型，通过构建统一的影像成像模型完成各
CCD
影像的对应关系构建，最终实现影像的高精度拼接
。
这种方法可以确保影像几何无缝，并且消除影像内部畸变，但是没有其他辅助数据，拼接影像的几何定位精度依赖卫星量测数据，导致最终的几何精度不甚理想
。

技术实现思路

[0003]本公开实施例的目的在于提供一种激光测高数据辅助的卫星线阵影像拼接方法及系统，从而解决现有技术中存在的前述问题
。
[0004]为了实现上述目的，本公开实施例采用的技术方案如下：
[0005]本公开实施例一方面提供了一种激光测高数据辅助的卫星线阵影像拼接方法，所述方法包括：
[0006]构建相机各分片真实电荷耦合器件
CCD
的严密几何成像模型
。
[0007]在各分片
CCD
影像重叠区进行连接点匹配，获取分布均匀
、
一定数量的连接点
。
[0008]获取影像成像时间范围内的激光测高数据，得到激光测高数据在
CCD
影像上的像点坐标
。
[0009]根据连接点像点坐标和物方坐标，以及激光测高点的像点坐标
、
物方坐标和测距值解算所述真实
CCD
的严密几何成像模型的外参数，更新所述真实
CCD
的严密几何成像模型
。
[0010]根据所述真实
CCD
的严密几何成像模型的内参数，在相机焦平面上进行直线拟合，构建虚拟
CCD
，确定所述虚拟
CCD
的每个探元的严密几何成像模型内参数，结合相机姿态
、
轨道测量值以及解算后的所述真实
CCD
的严密几何成像模型的外参数，构建虚拟
CCD
严密几何成像模型
。
[0011]根据所述虚拟
CCD
严密几何成像模型
、
所述真实
CCD
的严密几何成像模型以及高程
为相应的改正数系数矩阵，
L1和
P1分别是初值计算常量和权矩阵
。
[0028]可选的，所述激光测高点的测距信息的测距方程为：
[0029][0030]其中，
(X,
′
Y
′
,Z
′
)
T
是激光测高点在地面物方坐标系下的物方坐标，
(X
S
,Y
S
,Z
S
)
T
是卫星本体坐标系中心在地面物方坐标系的物方坐标，
ρ
是激光测高点的精确测距值；
[0031]所述根据激光测高点的测距信息，对激光测高点构建距离误差方程
(6)
：
[0032][0033]将公式
(6)
写成矩阵形式为：
[0034]V2＝
B2x
‑
L2,P2ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(7)
[0035]其中，
V2＝
v
las
是激光测高点距离误差值，
x
＝
(
Δ
X
Δ
Y
Δ
Z)
T
是激光测高点地面物方坐标改正数，
B2为相应的改正数系数矩阵，
L2和
P2分别是初值计算常量和权矩阵
。
[0036]可选的，所述对前两步骤构建的误差方程进行法化，利用最小二乘法求解真实
CCD
的严密几何成像模型外参数和连接点
、
激光测高点物方坐标改正数之前，结合所述连接点和激光测高点误差方程公式
(4)
和所述激光测高点的距离误差方程公式
(7)
构建整体误差方程
:
[0037]V
＝
At+Bx
‑
L,P
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(8)
[0038]式中，
V
＝
(V
1 V2)
T
表示观测值残差向量，
A,B
为相应的改正数系数矩阵，
L
和
P
分别是初值计算常量和权矩阵，
x
是连接点和激光测高点物方坐标改正数
。
[0039]可选的，所述获取影像成像时间范围内的激光测高数据，得到激光测高数据在
CCD
影像上的像点坐标，包括：
[0040]将激光测高点物方坐标代入构建的所述真实
CCD
的严密几何模型中，得到激光测高点在各分片
CCD
影像上的概略像点坐标
。
[0041]将所述得到的激光测高点在
CCD
影像上的概略像点坐标与足印影像中的激光测高点像点坐标，在足印大小范围内进行定点匹配，得到激光测高点在
CCD
影像上的精确像点坐标
。
[0042]可选的，所述获取影像成像时间范围内的激光测高数据，得到激光测高数据在
CCD
影像上的像点坐标，还包括：
[0043]如果定点匹配失败，则在足印大小范围内对足印影像和
CCD
影像进行同名点匹配，建立
CCD
影像与足印影像的映射关系，代入激光测高点在足印影像本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.
一种激光测高数据辅助的卫星线阵影像拼接方法，其特征在于，所述方法包括：构建相机各分片真实电荷耦合器件
CCD
的严密几何成像模型；在各分片
CCD
影像重叠区进行连接点匹配，获取分布均匀
、
一定数量的连接点；获取影像成像时间范围内的激光测高数据，得到激光测高数据在
CCD
影像上的像点坐标；根据连接点像点坐标和物方坐标，以及激光测高点的像点坐标
、
物方坐标和测距值解算所述真实
CCD
的严密几何成像模型的外参数，更新所述真实
CCD
的严密几何成像模型；根据所述真实
CCD
的严密几何成像模型的内参数，在相机焦平面上进行直线拟合，构建虚拟
CCD
，确定所述虚拟
CCD
的每个探元的严密几何成像模型内参数，结合相机姿态
、
轨道测量值以及解算后的所述真实
CCD
的严密几何成像模型的外参数，构建虚拟
CCD
严密几何成像模型；根据所述虚拟
CCD
严密几何成像模型
、
所述真实
CCD
的严密几何成像模型以及高程数据，对线阵卫星影像进行拼接
。2.
根据权利要求1所述的方法，其特征在于，所述根据连接点像点坐标和物方坐标，以及激光测高点的像点坐标
、
物方坐标和测距值解算所述真实
CCD
的严密几何成像模型的外参数，更新所述真实
CCD
的严密几何成像模型，包括：分别对所述真实
CCD
的严密几何成像模型外参数
、
连接点和激光测高点物方坐标构建误差方程；根据激光测高点的测距信息，对激光测高点构建距离误差方程；对前两个步骤构建的误差方程进行法化，利用最小二乘法求解真实
CCD
的严密几何成像模型外参数和连接点
、
激光测高点物方坐标改正数；根据改正数更新几何成像模型外参数和连接点以及激光测高点物方坐标，重复以上步骤直至迭代结束，得到最终真实
CCD
的严密几何成像模型
。3.
根据权利要求2所述的方法，其特征在于，所述
CCD
的严密几何成像模型为：其中，
(X,Y,Z)
T
是影像点在地面物方坐标系下的物方坐标，
(X
S
,Y
S
,Z
S
)
T
是卫星本体坐标系中心在地面物方坐标系的物方坐标，
m
是成像比例系数，是卫星本体坐标系到
J2000
惯性坐标系的转换矩阵，是
J2000
惯性坐标系到物方坐标系的转换矩阵；
(
ψ
x
,
ψ
y
)
是影像探元的指向角，为影像成像严密几何模型内参数，代表了影像内部精度，
R
u
是影像成像严密几何模型外参数构成的矩阵，代表了影像的外部精度；严密几何成像模型外参数矩阵
R
u
由公式
(2)
确定，具体形式如下：
其中，
ω
,
κ
分别是围绕
Y,X,Z
轴旋转角；所述分别对所述真实
CCD
的严密几何成像模型外参数
、
连接点和激光测高点物方坐标构建误差方程，得到公式
(3)
：将公式
(3)
写成矩阵形式为：
V1＝
At+B1x
‑
L1,P1ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(4)
其中，
V1＝
(v
x v
y
)
T
是观测值投影误差，是影像严密几何成像模型外参数改正数，
x
＝
(
Δ
X Δ
Y Δ
Z)
T
是激光测高点或连接点地面物方坐标改正数，
A,B1为相应的改正数系数矩阵，
L1和
P1分别是初值计算常量和权矩阵
。4.
根据权利要求2或3所述的方法，其特征在于，所述激光测高点的测距信息的测距方程为：其中，
(X,
′
Y
′
,Z
′
)
T
是激光测高点在地面物方坐标系下的物方坐标，
(X
S
,Y
S
,Z
S
)
T
是卫星本体坐标系中心在地面物方坐标系的物方坐标，
ρ
是激光测高点的精确测距值；所述根据激光测高点的测距信息，对激光测高点构建距离误差方程
(6)
：将公式
(6)
写成矩阵形式为：
V2＝...

【专利技术属性】
技术研发人员：刘昌儒，李鸿洲，郭莉，陈颖，尹烁，王艺颖，岳明宇，
申请(专利权)人：自然资源部国土卫星遥感应用中心，
类型：发明
国别省市：