一种根据激光雷达数据自动反演水平能见度的方法技术

本发明专利技术公开了一种根据激光雷达数据自动反演水平能见度的方法，该方法包括：使用激光雷达进行水平大气光学特性测量，获取激光雷达的测量数据；对激光雷达的测量数据进行数据处理，获得相对应的距离校准对数函数；从距离校准对数函数中选取符合预设线性特性的线性区域；对所选取的线性区域进行消光系数反演，得到该线性区域的线性拟合斜率，并根据所述线性拟合斜率计算得到水平能见度。通过使用本发明专利技术中的根据激光雷达数据自动反演水平能见度的方法，可以有效、稳定地进行水平能见度的自动反演。

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

【技术保护点】
一种根据激光雷达数据自动反演水平能见度的方法，其特征在于，该方法包括：使用激光雷达进行水平大气光学特性测量，获取激光雷达的测量数据；对激光雷达的测量数据进行数据处理，获得相对应的距离校准对数函数；从距离校准对数函数中选取符合预设线性特性的线性区域；对所选取的线性区域进行消光系数反演，得到该线性区域的线性拟合斜率，并根据所述线性拟合斜率计算得到水平能见度。

【技术特征摘要】
1.一种根据激光雷达数据自动反演水平能见度的方法，其特征在于，该方法包括：使用激光雷达进行水平大气光学特性测量，获取激光雷达的测量数据；对激光雷达的测量数据进行数据处理，获得相对应的距离校准对数函数；从距离校准对数函数中选取符合预设线性特性的线性区域，具体包括：从距离校准对数函数中选取一个初始区域；将所述初始区域中的数据划分为多个数组，每个数组中包括N个数据；采用最小二乘法进行线性拟合，计算每个数组的斜率；根据各个数组的斜率，根据预设起始条件确定线性区域的起始位置n_start；从起始位置开始，计算数据组(n_start：n_start+N)的斜率；依次增加N个数据，分别计算数据组(n_start：n_start+N×(j-1))的斜率，其中j＝1，2，…，50；将计算得到的各个斜率中的最小值所对应的数组的数据终止位置设置为线性区域的终止位置；对所选取的线性区域进行消光系数反演，得到该线性区域的线性拟合斜率，并根据所述线性拟合斜率计算得到水平能见度。2.根据权利要求1所述的方法，其特征在于：所述N的值为50。3.根据权利要求2所述的方法，其特征在于，所述起始条件为：起始位置的第一个斜率为负值，且起始位置的后三个斜率中只有一个为正值。4.根据权利要求1所述的方法，其特征在于，所述激光雷达的测量数据为距离r处的大气散射回波信号能量P(r)：

【专利技术属性】
技术研发人员：刘兴润，李霞，刘浩，雷浩，
申请(专利权)人：北京环境特性研究所，
类型：发明
国别省市：北京;11