[发明专利]一种气溶胶厚度与激光雷达测距延迟的自动标定方法

说明书

技术领域

本发明属于遥感对地观测技术领域，具体涉及一种气溶胶厚度与激光雷达测距延迟的自动标定方法。

背景技术

为了提高立体测图卫星的精度，亟需解决激光测高仪测高测距数据不准确的这一关键技术问题。其中最重要内容就是定期对激光测高数据的水汽延迟量的大小进行标定，并以此进行水汽延迟改正，从而获取精确的测距数据，从而提高激光测距的精度。

现有的激光测距改正技术需要依靠NCEP提供的准确大气数据，进而对测距进行精确的改正。这种方法不适用于国产激光测高卫星，一是由于NCEP并不及时对中国用户开放，造成国产激光测高卫星由于缺乏大气辅助数据而不进行精确大气改正的情况，从而无法获得准确的测距数据，最终导致国产立体测图卫星的整个工程任务的失败；第二个因素是NPEC数据空间分辨率有限，其1*1°的分辨率很难区分出激光传输路径上的各种小云层等等，使用NPEC进行大气改正的精度有限。

通过各种国产数据获取大气气溶胶含量的分布，并改正国产立体测图卫星的激光测距值成为测距精确修正的最可行手段。新一代立体测图卫星含有双线阵和多光谱影像，与激光测距仪基本同时工作，这些影像上准确捕获了激光测距时的瞬时气溶胶分布情况。从这些影像上提取出的气溶胶数据，能够最准确反映出激光发射/反射路径上的真实气溶胶分布。然而，通过图像处理方法识别出的气溶胶是一个无量纲数值，不代表任何物理含义，需要进行标定。

由于激光落脚点的不可见性质，且1064nm波长的光斑到达地面的位置难以预知，在陆地表面进行气溶胶的方法具有极大的技术挑战，然而海平面的实时高程是可以通过潮汐模型精确计算出来的，因此能够通过海面的高程和激光的发射位置精确计算出实际的激光路径长度，从而计算出实际大气延迟量，进而完成标定工作。这一过程完全不依赖于地面实测值，具有高效、高精度、价廉的多重优点。

通过提取的气溶胶厚度，与海表计算得到的实际大气延迟数值进行回归分析，就可以得到不同气溶胶厚度情况下的气溶胶像元值与激光测距值改正数的关系。从而在轨工作期间计算出各个测量数据的实际测距改正值，以此为国产高分辨率立体测图卫星系统的激光器提供精确的测距值，最终获取高精度的控制点数据。这对于提升我国对地观测卫星数据几何处理精度和产品质量具有极为重要的意义。

发明内容

本发明所要解决的问题是：在不依赖于地面测量的情况下，提供一种高精度、便捷的自动标定方法，实时标定出气溶胶对大气测距延迟的影响大小，从而大幅度提高陆地表面激光测高数据的准确性。

本发明提供的技术方案采用如下步骤：

(1)从国产立体测图卫星在海洋测区的足印与多光谱图像上提取气溶胶含量分布图；分布图提取方法步骤如下：

(1-1)将待识别气溶胶厚度的多光谱图像与辅助影像进行配准；

(1-2)对多光谱图像上的任意一个点在整副图像A上找到与其同类的像元，根据辅助影像筛选出类别不变的点；

(1-3)建立每一个像点多光谱图像与辅助影像回归关系，用辅助影像拟合出多光谱图像；

(1-4)拟合的亮度值除以真实亮度值，得到每一个像元的气溶胶含量；

(1-5)对各个波段得到的气溶胶含量图像低通滤波，对各波段进行加权平均得到最终的气溶胶含量分布图；

(2)对测区内测高数据进行潮汐改正；方法步骤如下：

(2-1)从NASA服务器中获取行星的三维笛卡儿太阳系质心坐标；

(2-2)利用美国海军天文台提供的C版本矢量天文子程序NOVASVer2.0，计算得到太阳和月亮任意时刻在协议地球坐标中的坐标；

(2-3)使用精密星历，通过拉格朗日多项式内插,就可以求得任意时刻卫星在协议地球坐标系中的坐标，逐步计算即可得到地球固体潮改正；

(3)使用全球大气分布模型确定大气成分和厚度，对测距值进行大气干项延迟改正；方法步骤如下：

(3-1)对在时间和空间上都不连续的气象数据进行内插，得到激光发射时刻和位置的气象数据；

(3-2)对经过内插计算得到的与测高仪时空间隔相同的气象数据，通过未经大气延迟修正的测高仪高程转换的位势高度寻找位势高度相邻的两个标准大气压力层，计算地表大气压数值积分；

(3-3)使用地表大气压数值对测高数据进行修正；

(4)通过海表高程数据和大气干项改正后的测高数据，计算气溶胶延迟量；方法步骤如下：

(4-1)根据激光发射时刻，用TOPEX/Poseidon数据计算出海平面在发射时刻的海拔高度；