一种基于光子点云辅助的卫星摄影测量对地定位方法技术

一种基于光子点云辅助的卫星摄影测量对地定位方法，它属于卫星摄影测量与遥感领域。本发明专利技术解决了现有的联合处理方法的定位精度低的问题。本发明专利技术方法的主要步骤为：获取目标区域的立体影像数据，并生成地物的摄影测量点云；获取光子点云并滤波，根据滤波后的光子点云生成高程采样点；将高程采样点云与摄影测量点云进行配准获得变换矩阵；根据变换矩阵获得控制点；根据获得的控制点对摄影测量点云的系统误差进行修正。本发明专利技术方法可以应用于卫星摄影测量与遥感领域。影测量与遥感领域。影测量与遥感领域。

【技术实现步骤摘要】
一种基于光子点云辅助的卫星摄影测量对地定位方法


[0001]本专利技术属于卫星摄影测量与遥感领域，具体涉及一种基于光子点云辅助的卫星摄影测量对地定位方法。

技术介绍

[0002]随着光学遥感卫星和数字摄影测量技术的发展，卫星摄影测量已成为一项被广泛应用于国防经济建设领域的遥感测绘技术。现今光学遥感卫星如资源三号、高分七号、wordview3/4、Pleiades和Cartosat
‑
3等，其地面分辨率均已全面迈入分米级水平。各类高分辨率遥感卫星陆续发射运行，为各行业领域提供数量丰富、质量优良的影像数据。同时，国内外学者对利用遥感影像进行地形测绘的理论方法进行了大量研究，形成了较为成熟的卫星摄影测量技术体系。
[0003]近些年，新兴的卫星激光测高技术上取得重大进展，对地观测定位能力大幅增强，为地形测绘等领域提供了强力的数据支撑。ICESat
‑
2卫星采用先进地形激光测高系统(Advanced Topographic Laser Altimeter，ATLAS)，具备强大的激光探测能力，能够采集观测频率更高、定位精度更高的地表激光点云数据。
[0004]由于光学卫星在成像的过程中，受到姿轨观测误差、平台颤振、大气折光和内部拼接等各种因素影响，其定位结果中含有复杂的系统误差，需要采集一定数量均匀分布的地面控制点，通过区域网平差提高定位精度。但对于如山林、沙漠和高原等区域，获取足够的高精度地面控制点非常困难，作业成本高昂。而卫星激光测高技术可以不受地形地貌的影响，在轨自动采集高精度的测高点数据。因此对ICESat
‑
2光子点云和光学遥感影像联合处理，利用激光测高点代替传统人工外业测量点，可大幅提高卫星摄影测量技术提取地表点云的效率。对于星载激光雷达数据与卫星遥感影像的联合处理，已有的研究表明，融合激光测高数据可以显著提高卫星摄影测量技术的定位精度。当前常用的联合处理方案是先从点云中筛选出位于平坦区域的激光点作为控制点，然后在像方坐标中，将控制点作为带权观测值进行区域网平差。该方案可有效提高无外业控制点下的卫星摄影测量定位精度，但面临以下待进一步解决的问题：
[0005](1)与传统脉冲式激光雷达不同，当前的激光测高卫星采用新型的光子计数激光雷达，其观测数据为点云形式。由于激光雷达对植被具有一定的穿透性，山林地区光子点云中含有来自树冠反射和地表反射的辐射信息，而光学遥感影像所见主要为树冠目标。
[0006](2)建立激光点与影像的准确对应关系是进行两种数据联合处理的前提条件，对此目前通常采用的是像方配准方案，即先选取平坦地区的激光点，利用原始的或经过自由网平差处理后的影像定位模型，计算激光点的像方投影点作为影像对应点，由于定位模型不可避免地存在误差，导致激光点物像对应关系不准确，不适用于对于地势崎岖的观测区域。
[0007]综上所述，由于光子点云中包含地表反射的辐射信息以及像方配准方案中影像定位模型存在误差，导致现有的联合处理方法的定位精度低。因此对卫星影像数据和光子雷
达数据进行联合处理时，需要从光子点云中分类提取出树冠高程，以保证目标高程的一致性。那么研究光子点云滤波分类算法以提取准确地表高程点，对于联合处理方案具有重要意义。并且鉴于现有联合处理方案的不足，有必要深入研究激光点与遥感影像的物方配准方案，为卫星影像数据和光子雷达数据的联合处理提供必要条件。

技术实现思路

[0008]本专利技术的目的是为解决现有的联合处理方法的定位精度低的问题，而提出的一种基于光子点云辅助的卫星摄影测量对地定位方法。
[0009]本专利技术为解决上述技术问题所采取的技术方案是：
[0010]一种基于光子点云辅助的卫星摄影测量对地定位方法，所述方法具体包括以下步骤：
[0011]步骤一、获取目标区域的立体影像数据；
[0012]步骤二、利用立体影像数据生成地物的摄影测量点云；
[0013]步骤三、获取立体像对区域内的ICESat
‑
2光子点云；
[0014]步骤四、对步骤三中获取的光子点云进行栅格滤波，滤除噪声点后，得到滤波后的光子点云；
[0015]步骤五、根据滤波后的光子点云等间隔生成高程采样点，再分别计算出每个高程采样点的地面高程值；
[0016]步骤六、根据各个高程采样点的地面高程值，将高程采样点云和步骤二生成的摄影测量点云进行配准，获得变换矩阵；再根据变换矩阵获得作为控制点的高程采样点；
[0017]步骤七、根据获得的控制点对摄影测量点云的系统误差进行修正。
[0018]本专利技术的有益效果是：
[0019]1、本专利技术提出基于栅格连续性约束的光子点云滤波算法，光子点云滤波及高程提取效果优良。对ICESat
‑
2/ATL03数据的试验结果表明，该算法能高效滤除噪声点并提取地表高程点，滤波成功率高于94％，高程提取精度优于5.4m，对不同地形变化、地物类型和激光能量的光子点云数据均可取得稳定处理效果。
[0020]2、本专利技术提出的基于线面状点云配准的联合处理算法，对摄影测量点云定位精度提高显著。试验结果表明引入ICESat
‑
2光子点云辅助后，资源三号摄影测点云的X、Y、Z方向精度分别达到0.46m、2.43m和2.13m，高分七号摄影测点云的X、Y、Z方向精度分别达到0.36m、1.34m和0.59m，定位精度得到显著提高。
附图说明
[0021]图1为本专利技术方法的流程图；
[0022]图2为ICESat
‑
2光子点云滤波及高程提取流程图；
[0023]图3为栅格连续性判断示意图；
[0024]图4为点云降趋势处理示意图；
[0025]图5为高程采样点示意图；
[0026]图6a为获取高程采样点的邻近点的示意图；
[0027]图6b为邻域点云高程频数统计示意图；
[0028]图7为光子点云与摄影测量点云配准流程图；
[0029]图8为高程相似度计算示意图；
[0030]图9a为资源三号摄影测量点云与ICESat
‑
2点云的示意图；
[0031]图9b为资源三号摄影测量点云系统误差修正前的示意图；
[0032]图9c为资源三号摄影测量点云系统误差修正后的示意图；
[0033]图10a为高分七号摄影测量点云与ICESat
‑
2点云的示意图；
[0034]图10b为高分七号摄影测量点云系统误差修正前的示意图；
[0035]图10c为高分七号摄影测量点云系统误差修正后的示意图。
具体实施方式
[0036]具体实施方式一、结合图1说明本实施方式。本实施方式所述的一种基于光子点云辅助的卫星摄影测量对地定位方法，所述方法具体包括以下步骤：
[0037]步骤一、获取目标区域的立体影像数据；
[0038]步骤二、利用立体影像数本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于光子点云辅助的卫星摄影测量对地定位方法，其特征在于，所述方法具体包括以下步骤：步骤一、获取目标区域的立体影像数据；步骤二、利用立体影像数据生成地物的摄影测量点云；步骤三、获取立体像对区域内的ICESat
‑
2光子点云；步骤四、对步骤三中获取的光子点云进行栅格滤波，滤除噪声点后，得到滤波后的光子点云；步骤五、根据滤波后的光子点云等间隔生成高程采样点，再分别计算出每个高程采样点的地面高程值；步骤六、根据各个高程采样点的地面高程值，将高程采样点云和步骤二生成的摄影测量点云进行配准，获得变换矩阵；再根据变换矩阵获得作为控制点的高程采样点；步骤七、根据获得的控制点对摄影测量点云的系统误差进行修正。2.根据权利要求1所述的一种基于光子点云辅助的卫星摄影测量对地定位方法，其特征在于，所述步骤四的具体过程为：步骤四一、栅格划分对整个光子点云剖面范围进行栅格划分，形成大小为w
×
h的格网；步骤四二、分别统计格网中每个栅格内包含的点数量，对于任意的一列栅格，从该列的全部栅格中选取出包含点数量最多的前t个栅格作为该列的候选栅格：其中，G
j
代表提取出的第j列中的候选栅格，为第j列中第m个栅格内包含的点数量，m为该列中的栅格总个数，代表从第j列中选取出包含点数量最多的t个栅格；同理，分别获得每列的候选栅格；步骤四三、栅格连续性检测步骤四三一、将第j列中的某个候选栅格作为待检测栅格，将第j
‑
1列、第j
‑
2列、
…
、第j
‑
k
’
列中的候选栅格以及第j+1列、第j+2列、
…
、第j+k
’
列中的候选栅格作为待检测栅格的邻域范围，待检测栅格与第j
‑
1列、第j
‑
2列、
…
、第j
‑
k
’
列中的候选栅格形成左连通域，待检测栅格与第j+1列、第j+2列、
…
、第j+k
’
列中的候选栅格形成右连通域；列中的候选栅格形成右连通域；式中，为左连通域内的栅格集合，为右连通域内的栅格集合，(i,j)为待检测栅格在格网中的行列序号，为左连通域内的栅格集合中的第m1个栅格，为右连通域内的栅格集合中的第n1个栅格；步骤四三二、计算待检测栅格的连续性得分Continue_score：
式中，为左连通域内的栅格数量，为右连通域内的栅格数量；步骤四三三、遍历计算出格网的第j列中每个候选栅格的连续性得分，将连续性得分最高的候选栅格以及连续性得分最高候选栅格的上、下各相邻的t'个栅格作为第j列中的信号栅格；步骤四三四、重复步骤四三一至步骤四三三三的过程，分别获得格网中每列的信号栅格；缩小单位栅格的宽度和高度，再基于缩小后单位栅格的宽度和高度对获得的信号栅格点云剖面范围进行栅格划分；步骤四四、根据步骤四三四获得的栅格划分结果，重复执行步骤四二和步骤四三的过程；直至缩小后单位栅格的宽度和高度均小于给定的最小值时停止，最后一次迭代所获得的信号栅格中包含的点作为滤波后的点云。3.根据权利要求2所述的一种基于光子点云辅助的卫星摄影测量对地定位方法，其特征在于，所述t的取值为3。4.根据权利要求3所述的一种基于光子点云辅助的卫星摄影测量对地定位方法，其特征在于，所述t'的取值为1。5.根据权利要求4所述的一种基于光子点云辅助的卫星摄影测量对地定位方法，其特征在于，所述缩小单位栅格的宽度和高度的具体方法为：width
i
’
+1
＝width
i
’
/r
w
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)height
i
’
+1
＝height
i
’
/r
h
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)式中，width
i
’
为第i
’
次滤波的单位栅格宽度，height
i
’
为第i
’
次滤波的单位栅格高度，r
w
为单位栅格宽度的缩放比例系数，r
h
为单位栅格高度的缩放比例系数。6.根据权利要求5所述的一种基于光子点云辅助的卫星摄影测量对地定位方法，其特征在于，所述步骤五的具体过程为：步骤五一、点云降趋势处理对滤波后的光子点云进行栅格划分，形成大小为w0
×
h0的格网；对于格网中的第j列栅格，计算该列中所有点的高程值的中位数，将中位数作为该列的趋势常数项；再分别利用该列中每个点的高程值减去该列的趋势常数项，得到该列中每个点对应的降趋势后的高程值：值：其中，C_detrend
j
为第j列的趋势常数项，Grid
j
为w0
×
h0的格网中第j列内包含的点的高程值的集合，代表格网中第j列栅格内的第i”个点的高程值，i”＝1,2，
…
，I，I为第j列栅格内的点的总个数，为第j列栅格内的第i”个点降趋势后的高程值；同理，分别获得格网中每个点降趋势后的高程值；步骤五二、生成等间隔高程采样点在降趋势后点云的沿轨距离方向上，根据降趋势后点云沿轨距离的最小值和最大值，等间隔内插形成高程采样点；
式中，为降趋势后点云沿轨距离的最大值，为降趋势后点云沿轨距离的最小值，d
sample
为相邻的两个高程采样点之间的间距，N
sample
为高程采样点的总数量，为第k个高程采样点的沿轨距离；步骤五三、获取邻近点对于第k个高程采样点，按给定的邻域距离，选取第k个高程采样点的邻近点集合：其中，为第k个高程采样点的邻近点集合，表示集合内的第n个点，为邻近点集合内的点的沿轨距离，d
ε
为邻域距离；步骤五四、统计邻近点集合的频数分布根据集合内降趋势后点云高程的最大值、最小值，并设置组距，统计邻近点集合的频数分布直方图；步骤五五、计算点云分类阈值其中，bin
j
...

【专利技术属性】
技术研发人员：余俊鹏，吴伟东，
申请(专利权)人：广东工业大学，
类型：发明
国别省市：