一种基于半导体激光云高仪的云高自动反演方法
【技术领域】
[0001] 本发明涉及云高测量领域，具体为一种基于半导体激光云高仪的云高自动反演方 法。
【背景技术】
[0002] 云是悬浮在大气中的小水滴、过冷水滴、冰晶或它们的混合物组成的可见聚合物； 有时也包含一些较大的雨滴、冰晶和雪晶，其底部不接触地面。云的生成、外形特征、云量、 云的分布及其演变不仅反映了当时大气的运动、稳定程度和水汽状况，同时预示着未来天 气变化。云层厚度对温室效应有着重要影响，云层的性质和天气特征分析对人工增雨抗旱 减灾有着重要指导意义，云层信息尤其云高参数是影响飞机起飞、着陆、飞行安全的重要气 象要素之一。如何方便有效地准确确定云层高度并监测其变化过程，对天气系统发展、降雨 形成和大气福射传输模式研宄具有十分重要意义。
[0003] 目前对云高的测量主要是依赖于目测，要比较准确的判断、辨别各种云的高度，需 要掌握云的外形特征、了解云的物理过程，对观测人员业务知识要求高，并且人为误差较 大。激光云高仪采用激光原理和微弱信号检测技术有效提取云和大气的后向散射信号，能 够实时提取云高的信息，具有智能、便携、实时和动态范围大等优点。基于半导体激光云高 仪的反演主要有微分零交叉法、阈值法。微分零交叉法对大气和云的后向散射信号微分，寻 找微分信号零点位置，零点位置对应云峰的高度，但是，在复杂的大气环境下气溶胶的干扰 将会产生很多过零点，并且白噪声的抖动也会引入额外的过零点，这势必对过零点的有效 提取造成很大的影响；阈值法是首先选择一个适当的阈值，以超过设定阈值作为判断基准， 阈值的选取基于经验的方法，而实际的大气是剧烈变化的，白天和夜晚的噪声干扰差别明 显，因此，基于经验值的阈值法很难适应于复杂的大气环境，白噪声和气溶胶有时也会产生 超过设定阈值的奇异值，阈值法没有对奇异值作进一步的处理和筛选，这也将在云高的反 演中带入一定的误差。
【发明内容】
[0004] 本发明要解决的技术问题是：克服现有技术的不足，提供一种基于半导体激光云 高仪的云高自动反演方法，以解决传统方法中云高观测人为误差较大，受气溶胶和白噪声 干扰大，鉴别能力不高，自适应性不强等问题。
[0005] 为达到上述目的，本发明所采用的技术方案为：
[0006] 一种基于半导体激光云高仪的云高自动反演方法，其特征在于实现步骤如下：
[0007] (1)半导体激光云高仪接收l〇64nm后向散射回波信号；
[0008] ⑵由l〇64nm后向散射回波信号特征获取分段高度；
[0009] (3)对采集到的1064nm后向散射回波信号进行距离平方修正，并将距离平方修正 后的信号取名为PRR信号，分别对l〇64nm后向散射回波信号和PRR信号进行归一化处理；
[0010] (4)根据分段信号分布特征获取滑动窗口积分阈值以及云顶、云底、云峰、云厚度 和云幅度变化阈值；
[0011] (5)米用滑动窗口积分的方法，从地面开始向上逐一滑动窗口，在窗口内对分段信 号进行累加平均，当信号累加平均值大于等于滑动窗口积分阈值时，初步筛选出云层分布 的大致位置；
[0012] (6)由分段信号分布特征确定云底、云峰、云顶的精确高度；
[0013] (7)利用云厚度阈值、云峰阈值、云幅度变化阈值对云层进行进一步的筛选，当经 过（6)步骤确定的云厚、云变化幅度、云峰的幅度其中两项同时大于等于对应的阈值时，即 筛选过后的云高分布信息。
[0014] 本发明与现有技术相比的有益效果：
[0015] (1)本发明通过选取背景信号计算出标准偏差，以此为基准对1064nm后向散射回 波信号斜率进行自动化判别，有效识别分段高度，并利用分段高度对l〇64nm后向散射回波 信号和PRR信号分段处理，有利于阈值的获取和设定。
[0016] (2)本发明可以在复杂的大气环境条件下，有效减小气溶胶和白噪声在云高分布 测量中的干扰，自适应地确定云高识别阈值，高准确性地提取云高分布信息。
[0017] (3)本发明计算便于程序化，可以长时间在线高时间分辨率和高空间分辨率的对 云高分布进行主动遥感测量，这是传统仪器所无法实现的。
【附图说明】
[0018] 图1为本发明方法实现的流程图；
[0019] 图2为1064nm后向散射回波信号和PRR信号；
[0020] 图3为基于云高自动反演方法的云高仪云高自动监测软件；
[0021] 图4黑点代表了云层时空分布演化图。
【具体实施方式】
[0022] 根据本发明中提出的一种基于半导体激光云高仪的云高自动反演方法，以某次实 验数据为例进行分析，所述方法步骤为：
[0023] (1)半导体激光云高仪接收l〇64nm后向散射回波信号，激光云高仪的1064nm后向 散射回波信号可表示为：
[0024]
[0025] 其中，P(A，Z)为波长A，距离z处的激光云高仪l〇64nm后向散射回波信号，C为 系统常数，0 (X，z)为气体分子和云或气溶胶的总的后向散射系数，a (A，z)为气体分子 和云或气溶胶的消光系数。由于云的后向散射系数远大于气溶胶，所以当激光在大气中传 输与云产生后向散射作用时，P(A，z)将会突然增大，出现突变，呈现出周围气溶胶回波信 号截然不同的分布特征，主要表现在：云层回波信号幅度较大；云层回波信号变化幅度大， 因此云层内部标准偏差较气溶胶回波信号标准偏差明显偏大；云层回波信号脉冲宽度宽， 而干扰噪声具有随机性，脉冲很窄；云层内部可能剧烈变化，表现为云峰出现多个，但认为 同一层云。因此，针对云层回波信号独特的特征设定六个阈值，针对云层回波信号幅度，设 定滑动窗口积分阈值和云峰阈值；针对云层回波变化幅度，设定云幅度变化阈值；针对云 层回波信号宽度，设定云厚度阈值；为了有效识别出云底、云顶信息设定云底、云顶阈值； 以六个阈值为判断依据从而由l〇64nm后向散射回波信号中筛选出云层信息。
[0026] (2)由1064nm后向散射回波信号特征获取分段高度，选取12km~15km的1064nm 后向散射回波信号作为背景信号，计算背景信号的标准偏差；由近地面高度向上将l〇64nm 后向散射回波信号与背景信号的标准偏差逐一进行比较，当两者比值第一次小于等于1时 所对用的高度为分段高度；
[0027] (3)受气溶胶干扰和几何因子的影响，1064nm后向散射回波信号在分段高度以下 呈现先上升后下降的变化特征，不利于阈值的选取和设定，而l〇64nm后向散射回波信号经 过距离平方修正后的PRR信号平整形较好，有利于阈值的选取和设定，所以首先对l〇64nm 后向散射回波信号进行距离平方修正、归一化处理、三角滤波后得到PRR信号，如下所示：
[0028]
[0029] 利用PRR信号确定六个阈值；在此高度以上的1064nm后向散射回波信号随距离变 化不明显，而PRR受背景噪声影响，信号抖动随距离增加而加剧，因此直接利用归一化处理 后的1064nm后向散射回波信号确定六个阈值。
[0030] (4)扣除奇异点后，求分段信号的平均值和标准偏差作为滑动窗口积分阈值以及 云顶、云底、云峰、云厚度和云幅度变化阈值的选择基准，以此基准为基础乘以是适当的系 数作为阈值。所述分段信号在分段高度以下为PRR信号，在分段高度以上为归一化处理后 的1064nm后向散射回波信号；所述滑动窗口积分阈值设为分段信号平均值加上标准偏差 的2~5倍；所述云顶、云底阈值相同，并设为分