一种根据激光雷达数据自动反演水平能见度的方法
【技术领域】
[0001] 本发明涉及目标光学特征研究技术，特别涉及一种根据激光雷达数据自动反演水 平能见度的方法。
【背景技术】
[0002] 能见度是在一定的大气透明度下，正常人眼能发现水平天空为背景的黑色目标物 (视场角大于30° )的最大距离。能见度是大气探测的一个重要因素，可反映大气层稳定 程度、研究大气污染等。
[0003] 能见度测量方法主要经历了由人工目测到仪器测量的发展过程。最初，能见度一 般采用人工目测法测量；随着器测研发的进展，先后出现了数学摄像式、透射式、散射式能 见度仪和激光雷达大气能见度测量设备。
[0004] 激光雷达是传统雷达技术和现代激光技术相结合的产物，它以激光作为载波，可 以用振幅、频率、相位和偏振来搭载信息，是一种极为有效的主动式现代光学遥感设备，被 广泛地应用于大气科学、环境保护、气象与气候和大气遥感等领域。
[0005] 1961年激光器问世不久，以测距为主要功能的激光雷达便诞生了，60年代后期， 研制重点转向单脉冲激光雷达和〇)2激光雷达。随着二极管激光栗浦的发展，美国开发了 以DPSS微脉冲激光雷达，并在此基础上又开发了差分吸收、拉曼散射、多普勒激光雷达。目 前，我国正紧跟国际先进技术，积极开展大气探测激光雷达的研究，逐步缩小与国外高水平 国家的差距。
[0006] 激光雷达通过激光测量大气消光系数，进而推算能见度，相对而言是目前最为客 观和准确的能见度及其他光学参量检测设备。
[0007] 因此，在本领域中，寻找一种有效、稳定、精度高的反演方法一直都是大气探测激 光雷达研究中的一项重要目标。迄今为止，在利用激光雷达的回波信号反演大气消光系数 的方法中，比较有效的有Collis斜率法以及Klett数据分析反演方法。
[0008] 激光雷达作为一种新型的大气观测工具，可以通过直接探测激光与大气相互作用 的光辐射信号来定量确定大气能见度。在利用激光雷达反演水平大气能见度时，一般可采 用斜率法。在所述斜率法中，通过拟合直线斜率可导出大气的消光系数，利用消光系数可获 得水平大气能见度。
[0009] 在激光雷达进行水平测量时，近地面线性区域的选择对计算能见度有比较大的影 响，在现有技术中，一般都是采用手动处理的方式来利用激光雷达数据获取水平大气能见 度，因此无法有效、稳定地进行水平能见度的自动反演。
【发明内容】
[0010] 有鉴于此，本发明提供一种根据激光雷达数据自动反演水平能见度的方法，从而 可以有效、稳定地进行水平能见度的自动反演。
[0011] 本发明的技术方案具体是这样实现的：
[0012] -种根据激光雷达数据自动反演水平能见度的方法，该方法包括：
[0013] 使用激光雷达进行水平大气光学特性测量，获取激光雷达的测量数据；
[0014] 对激光雷达的测量数据进行数据处理，获得相对应的距离校准对数函数；
[0015] 从距离校准对数函数中选取符合预设线性特性的线性区域；
[0016] 对所选取的线性区域进行消光系数反演，得到该线性区域的线性拟合斜率，并根 据所述线性拟合斜率计算得到水平能见度。
[0017] 可选的，所述从距离校准对数函数中选取符合预设线性特性的线性区域包括：
[0018] 从距离校准对数函数中选取一个初始区域；
[0019] 将所述初始区域中的数据划分为多个数组，每个数组中包括N个数据；
[0020] 采用最小二乘法进行线性拟合，计算每个数组的斜率；
[0021 ] 根据各个数组的斜率，根据预设起始条件确定线性区域的起始位置restart ;
[0022] 从起始位置开始，计算数据组（n_start :n_start+N)的斜率；
[0023] 依次增加 N个数据，分别计算数据组（n_start :n_start+NX (j-Ι))的斜率，其中 j = 1，2，...，50 ;
[0024] 将计算得到的各个斜率中的最小值所对应的数组的数据终止位置设置为线性区 域的终止位置。
[0025] 可选的，所述N的值为50。
[0026] 可选的，所述起始条件为：
[0027] 起始位置的第一个斜率为负值，且起始位置的后三个斜率中只有一个为正值。
[0028] 可选的，所述激光雷达的测量数据为距离r处的大气散射回波信号能量P (r):
[0030] 其中，EO是发射的激光脉冲能量；c是光速;Ar是接收望远镜的有效面积；Y(r) 是系统的几何重叠因子，β (r)是距离r处的大气后向散射系数；τ SyS是系统的总透过率； xa(r)是大气透过率；其中，
[0032] 所述σ (r)是r处的大气消光系数。
[0033] 可选的，所述距离校准对数函数为：
[0034] S (r) = In [r2P (r)],
[0035] 其中，S(r)为距离校准对数函数，
[0036] 可选的，所述距离校准对数函数S (r)满足如下的方程：
[0038] 其中，β (r)、σ (r)分别为气溶胶的后向散射系数和消光系数，ra为参考距离；所 述方程的微分形式为：
[0040] 可选的，水平方向上大气均匀分布时的消光系数〇为：
[0043] 可选的，根据如下公式计算得到水平能见度：
[0045] 其中，Vis为水平能见度。
[0046] 如上可见，在本发明的根据激光雷达数据自动反演水平能见度的方法中，由于可 以对激光雷达的测量数据进行数据处理，获得相对应的距离校准对数函数，然后从距离校 准对数函数中选取符合预设线性特性的线性区域，再对所选取的线性区域进行消光系数反 演，得到该线性区域的线性拟合斜率，并根据所述线性拟合斜率计算得到水平能见度，从而 实现了激光雷达测量数据距离校准对数函数线性区域的自动选取，避免了人工干扰的不确 定性，实现了激光雷达数据自动反演能见度，可以有效、稳定、快速、高效地进行水平能见度 的自动反演。
【附图说明】
[0047] 图1为本发明实施例中的根据激光雷达数据自动反演水平能见度的方法的流程 示意图。
[0048] 图2为本发明实施例中激光雷达的距离校准对数函数分布图的示意图。
[0049] 图3为本发明实施例中水平能见度的自动反演的计算结果的示意图。
【具体实施方式】
[0050] 为使本发明的目的、技术方案及优点更加清楚明白，以下参照附图并举实施例，对 本发明进一步详细说明。
[0051] 本实施例提供了一种根据激光雷达数据自动反演水平能见度的方法。
[0052] 在本发明的技术方案中，可以对激光雷达的水平测量数据的特点进行分析。
[0053] 当激光雷达发射的激光光束通过大气时，会与大气中的分子及气溶胶发生各种作 用，其中散射的作用最为明显。激光雷达探测器接收到距离r处的大气散射回波信号能量 P (r)可由激光雷达方程确定：
[0055] 其中，EO是发射的激光脉冲能量；c是光速；Ar是接收望远镜的有效面积；Y(r) 是系统的几何重叠因子，β (r)是距离r处的大气后向散射系数；τ SyS是系统的总透过率； xa(r)是大气透过率，定义为
[0057] 根据上述公式（2)可知，激光在传输过程中受到大气分子和气溶胶粒子的衰减程 度。其中σ (r)是r处的大气消光系数。
[0058] 在现有技术中，Collis针对气溶胶散射很强且分布均匀的情况，提出了一种求解 激光雷达方程的简单方法，即Collis斜率法。
[0059] 在Collis斜率法中，将对经距离平方校准后的激光雷达测量数据取对数：
[0060] S (r) = In [r2P (r) ] (3)
[0061] 其中，S(r)为距离校准对数函数。S(r)满足如下所述的方程：
[0063] 其中，β (r)、σ (r)分别为气溶胶的后向散射系数和消光系数，ra为参考距离。上 式的微分形式为：
[0065] 当水平方向上的同层大气处于均匀时，β (r)、〇 (r)均为常数，不随距离的变化而 改变，即存在
，由此可以得出水平方向上大气均匀分布时的消光系数〇 :
[0067] 消光系数σ与水平能见度Vis之间的关系可用下面的公式表述：
[0071] 此种方法原理明晰、简单，计算量小，因而在大气均匀分布时，是求解大气水平消 光系数和能见度的有效方法。
[0072] 图1为本发明实施例中的根据激光雷达数据自动反演水平能见度的方法的流程 示意图。如图1所示，本发