一种红外与激光联合的地基测云方法
【专利摘要】红外与激光联合的地基测云方法，包括以下步骤：(1)分别利用非致冷红外焦平面阵列传感器获得大气向下红外辐射数据，利用激光传感器获得天顶后向消光系数廓线数据；获得上述两个数据的时间同步；(2)联合上述数据对云下水汽和气溶胶辐射进行估计，并利用辐射传输模式计算的晴空阈值进行初步的云检测；并假定云为黑体，反演得到云底高；(3)对高时间分辨率的红外辐射图像进行序列分析，结合晴空阈值做进一步的云检测，并计算云量；(4)拟合红外辐射反演的云底高与激光测量的云底高之间的比例系数；(5)进行全视场云底高修正，并计算得出每10分钟的典型云底高度。
【专利说明】一种红外与激光联合的地基测云方法
【技术领域】
[0001]本发明涉及一种地基测云方法，尤其是一种联合地基红外与激光进行云检测和全视场云底高修正的方法，其主要用于地面气象观测自动化测量中。
【背景技术】
[0002]云通常覆盖地球大气约50%左右，是地气系统辐射收支的主要调节者。云的观测工作一直是气候研究、天气分析与预报工作的基础。长期以来，目测是气象业务中常用的估计云量、确定云高的方法，但由于人为观测误差大，加之夜间光照条件差，云的观测资料的准确性和昼夜一致性难以保证，同时对实际天空的描述也较为困难。目前主要的地基测云技术体制有可见光被动遥感体制、红外被动遥感体制和激光主动遥感体制。可见光测云仪由于其夜间无法测量，很难应用于实际观测业务中。红外测云仪和激光云高仪能昼夜连续观测且各有特点。
[0003]地基红外被动遥感测云技术国外以美国ICI测云仪研究团队为代表，国内以中科院大气物理所和解放军理工大学气象海洋学院为代表，均取得大量创新性成果。研究表明，8-14 μ m波段向下红外辐射主要受到水汽和气溶胶的影响，单一红外测云仪由于对云下水汽和气溶胶的估计不足会造成云检测出现偏差。另外，红外辐射反演的云底高只适用近似为黑体的云，其他情况下由于无法得知云的光学厚度，反演的云底高存在很大误差。
[0004]地基激光主动遥感测云技术，国外以Vaisala公司为代表，其研制的激光云高仪在20世纪90年代应用于美国地面自动观测系统(ASOS)中替代人工目测。国内也有多家公司研发出了激光云高仪。一些学者研究发现，激光云高仪会在一些明显无云的情况下出现错误的有云信号，这可能是由于路径上存在一个很薄的湿气层或气溶胶层造成的。此外噪声的影响使其测量卷云这一弱回波信号时的数据获取率极低，这一问题在我国雾霾多发地体现的尤为明显。更重要的是，激光云高仪采用时间积分计算的云量存在严重的空间代表性问题。
[0005]综上所述，红外被动遥感测云和激光主动遥感测云这两种最有可能实现测云业务化的技术体制目前仍然在不同程度上存在问题，尚不能完全替代人工目测。
【发明内容】
[0006]本发明的目的是:提出一种将地基红外被动遥感和激光主动遥感联合进行云检测和全视场云底高修正的方法，从而最大程度地克服地基红外被动遥感测云和激光主动遥感测云的使用场合的不足，提高地基昼夜连续测云的准确性。
[0007]本发明的技术方案为:一种红外与激光联合的地基测云方法。该方法包括以下步骤:
[0008](I)分别利用非致冷红外焦平面阵列传感器获得大气向下红外辐射数据，利用激光传感器获得天顶后向消光系数廓线数据；获取数据的时间同步；(如通过数据采集系统每I分钟同步进行数据采集，并将数据存储。)
[0009](2)联合上述数据对云下水汽辐射和气溶胶辐射进行估计，并利用辐射传输模式计算的晴空阈值进行初步的云检测；并假定云为黑体，反演得到云底高；
[0010](3)对高时间分辨率的红外辐射图像进行序列分析，结合晴空阈值做进一步的云检测，并计算云量；
[0011](4)拟合红外辐射反演的云底高与激光测量的云底高之间的比例系数；
[0012](5)进行全视场云底高修正，并计算得出每10分钟的典型云底高度。
[0013]进一步，通过数据采集系统每I分钟同步进行非致冷红外焦平面阵列传感器数据采集和激光传感器的数据采集。
[0014]进一步，联合非致冷红外焦平面阵列传感器获取的红外辐射数据和激光传感器获取的天顶回波数据对云下水汽和气溶胶辐射进行估计，并计算出晴空辐射阈值。
[0015](2)的具体化:联合上述数据对云下水汽和气溶胶辐射进行估计，并利用辐射传输模式计算的晴空阈值进行初步的云检测；并假定云为黑体，反演得到云底高:
[0016]步骤1:云下气溶胶福射估计；将消光廓线转换为光学厚度廓线，代入Libradtran辐射传输模式，估算垂直路径的气溶胶向下辐射值I.。-;
[0017]步骤2:云下水汽辐射估计:利用激光无云时的天顶向下红外辐射Iiky和地面水汽密度Pwat?—vap■估算垂直路径的可降水含量PWV(PWV = f(IIE clear, pwat?—vap.))，其中函数f利用一段时间(例如不少于I年的时间)的探空廓线、天空红外辐射和地面水汽密度统计得出，具体统计模型可使用二元线性回归模型，并利用最小二乘法得到回归系数。将PffV代入Libradtran福射传输模式(但不限于Libradtran),即可估算垂直路径的水汽向下辐射值
Iwater_vapour °
[0018]步骤3:晴空福射阈值确定；利用Libradtran福射传输模式计算大气廓线以及臭氧等其他稳定少变气体的固定福射值I。，结合步骤1、2得到的IaMOSaL和IwatOT—vap.，计算出晴空辐射阈值Itl =
laerosol + Iwater—vapour+Ic ?
[0019]步骤4:初步的云检测；遍历红外辐射图像，判断实际测量的大气向下红外辐射I与Itl的大小，若实际测量的辐射I大于晴空辐射阈值Itl，则判为有云，否则无云；此步骤得到初步估计云量记为CCl ;并进一步得到修正后的云体向下红外辐射I1 = 1-1aerosol_Iwater_
-T
vapour 丄 c 0
[0020]步骤5:利用红外辐射反演云底高；根据步骤4中得到的I1，假设云为黑体，求出云层底亮温，根据该地区该时间段内的典型温度递减率，得到红外辐射反演的云底高ΗΙΚ。
[0021]对高时间分辨率的红外辐射图像进行序列分析，结合晴空阈值做进一步的云检测，并计算云量；具体步骤为:
[0022]步骤1:求取各像素点的序列辐射标准差。以I分钟为迭代周期，按10分钟为间隔进行图像的序列分析，将每个像素点的10次辐射数据求标准差σ ;利用辐射标准差判断每个像素随时间的均匀性；
[0023]步骤2:标记可疑像素云点。如步骤I求取各像素点的序列辐射标准差，遍历红外辐射标准差图像，得σ，若σ >0 ^，则此像素点标记为可疑云点，否则标记为晴空点，其中σ ^通过大量晴空样本统计得出。此步骤得到可疑云量记为CC2。
[0024]步骤3:对初步云检测进行修正:
[0025]对于⑵中判断为无云的区域，若(3)中步骤2得到的是可疑云区，则晴空辐射阈值Itl按0.2ff/ (m2.sr)向下调整，直到CCl与CC2的差值绝对值小于I ;
[0026]对于⑵中判断为有云的区域，且根据红外辐射反演的云底高在6000m以上，若
(3)中步骤2得到的是无云区，则晴空辐射阈值1(|按0.21/(!112.81.)向上调整，直到CCl与CC2的差值绝对值小于I ;对于其他情况不做修正。
[0027]通过此步骤，进一步提高了对高云的检测能力。
[0028]步骤4:按照调整后的晴空阈值进行云检测，并得到修正后的云量CC。
[0029]拟合红外辐射反演的云底高与激光测量的云底高之间的比例系数:具体方法是，统计10分钟内，天顶5°以内红外辐射反演的云底高Hik和最低层天顶激光云底高的中间值(10组数据)，若两类数据均存在，且大于等于3组，则拟合得到比例系数；若其中任一类不存在，则天顶扩大至10°，统计时间向前扩展至15分钟。若其中任一类仍不存在或者两类数据均存在的情况不足3组数据，则放弃拟合，利用前一次拟合结果作为本次拟合系数。将拟合系数记入拟合系数数据库，将有效的HIK、^asCT记入云底高数据库。
[0030]进行全视场云底高修正，并计算得出每10分钟的典型云底高度:具体步骤为；
[0031]步骤1:天空均匀性判断。利用大视场的红外辐射数据的标准差判断天空的均匀性，分为均匀天空和不均匀天空两类。
[0032]步骤2:全视场云底高修正。根据步骤I的结果，若为均匀天空，则直接利用拟合系数进行全视场云底高修正，得到全天空云底高分布图；若为不均匀天空，则根据待修正的红外辐射反演的云底高查找云底高数据库，将之前2小时内相差在500m以内的数据对找出，拟合得到修正比例因子，并将云底高乘以修正比例因子进行修正。
[0033]步骤3:典型云底高度的求取。按照低云(0-2500m)、中云(2500_6000m)和高云(6000m以上)三档分别进行统计。将每一档的云底高以10分钟为统计范围，由低到高排序，然后求得前2%所对应的云底高，作为每一档的典型云底高度。
[0034]与现有技术相比，本发明的有益效果在于:本发明公开了一种将地基红外被动遥感和激光主动遥感联合进行云检测和全视场云底高修正的方法。
[0035]本能更加准确地估计云下气溶胶和水汽的向下红外辐射，这将使得辐射传输模式计算的晴空阈值更加接近真实值，从而提高了云检测效果。进一步地，利用高时间分辨率的序列分析可以提闻对卷Z?的检测能力。结合激光Z?底闻可以修正红外反演的Z?底闻，不但能得到更加准确的全天空云底高分布，还可统计得出全天空典型云底高度信息。本发明通过充分联合地基测量的大气向下红外辐射数据、天顶激光回波数据、高时间分辨率的红外辐射序列图像，采用综合观测的理念，为新一代自动气象站中云的自动观测提供一种联合探测的技术途径，较好地克服了单一测云手段的不足。
【专利附图】
【附图说明】
[0036]图1红外与激光联合的地基测云方法流程图；
[0037]图2本发明得到的云量与只使用红外测云仪的云量结果对比；
[0038]图3本发明得到的修正后天顶云底高与只使用红外测云仪的天顶云底高对比。
【具体实施方式】
[0039]为了使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本发明进一步详细说明。
[0040]图1是本发明提出的红外与激光联合的地基测云方法流程图，所述方法的【具体实施方式】如下:
[0041](I)分别利用非致冷红外焦平面阵列传感器获得大气向下红外辐射数据，利用激光传感器获得天顶后向消光系数廓线数据。这一步骤中，关键在于获取数据的时间同步，通过数据采集系统每I分钟同步进行数据采集，并将数据存储。
[0042](2)联合上述数据对云下水汽和气溶胶辐射进行估计，并利用辐射传输模式计算的晴空阈值进行初步的云检测；并假定云为黑体，反演得到云底高。具体步骤有:
[0043]步骤1:云下气溶胶福射估计。将消光廓线转换为光学厚度廓线,代入Libradtran辐射传输模式，估算垂直路径的气溶胶向下辐射Ia_S()1。
[0044]步骤2:云下水汽辐射估计:利用激光无云时的天顶向下红外辐射Iiky和地面水汽密度Pwat?—vap■估算垂直路径的可降水含量PWV(PWV = f(IIE clear, Pwat?—vap.))，其中函数f利用一段时间(例如不少于I年的时间)的探空廓线、天空红外辐射和地面水汽密度统计得出，具体统计模型可使用二元线性回归模型，并利用最小二乘法得到回归系数。将PWV代入Libradtran福射传输模式(但不限于Libradtran),即可估算垂直路径的水汽向下辐射值
Iwater_vapour °
[0045]步骤3:晴空福射阈值确定。利用Libradtran福射传输模式计算大气廓线以及臭氧等其他稳定少变气体的固定福射值I。，结合步骤1、2得到的IaMOS()1和IwatOT—vap.，计算出晴空辐射阈值Itl =
laerosol + Iwater—vapour+Ic °
[0046]步骤4:初步的云检测。遍历红外辐射图像，判断实际测量的大气向下红外辐射I与Itl的大小，若实际测量的辐射I大于晴空辐射阈值Itl，则判为有云，否则无云。此步骤得到初步估计云量记为CC1。并进一步得到修正后的云体向下红外辐射I1 = 1-1aerosol_Iwater_
-T
vapour 丄 c 0
[0047]步骤5:利用红外辐射反演云底高。根据步骤4中得到的I1，假设云为黑体，求出云层底亮温，根据该地区该时间段内的典型温度递减率，得到红外辐射反演的云底高Hik。此步骤由于大多数利用红外辐射反演云底高均采用此方法，不做详细介绍。
[0048](3)对高时间分辨率的红外辐射图像进行序列分析，结合晴空阈值做进一步的云检测，并计算云量。具体步骤为:
[0049]步骤1:求取各像素点的序列辐射标准差。以I分钟为迭代周期，按10分钟为间隔进行图像的序列分析，将每个像素点的10次辐射数据求标准差O。利用辐射标准差判断每个像素随时间的均匀性。
[0050]步骤2:标记可疑云点。遍历红外辐射标准差图像，若σ>σ ^，则此像素点标记为可疑云点，否则标记为晴空点，其中Oci通过大量晴空样本统计得出。此步骤得到可疑云量记为CC2。
[0051]步骤3:对初步云检测进行修正。对于(2)中判断为无云的区域，若(3)中步骤2得到的是可疑云区，则晴空辐射阈值Itl按0.2ff/(m2.Sr)向下调整，直到CCl与CC2的差值绝对值小于I。对于(2)中判断为有云的区域，且根据红外辐射反演的云底高在6000m以上，若(3)中步骤2得到的是无云区，则晴空辐射阈值Itl按0.2W/(m2.sr)向上调整，直到CCl与CC2的差值绝对值小于I。对于其他情况不做修正。通过此步骤，进一步提高了对高云的检测能力。
[0052]步骤4:按照调整后的晴空阈值进行云检测，并得到修正后的云量CC。图2为使用本发明获得的云量与不使用本发明获得的云量差异。
[0053](4)拟合红外辐射反演的云底高与激光测量的云底高之间的比例系数。具体方法是:统计10分钟内，天顶5°以内红外福射反演的云底高Hik和最低层天顶激光云底高Huisot的中间值(10组数据)，若两类数据均存在，且大于等于3组，则拟合得到比例系数；若其中任一类不存在，则天顶扩大至10°，统计时间向前扩展至15分钟。若其中任一类仍不存在或者两类数据均存在的情况不足3组数据，则放弃拟合，利用前一次拟合结果作为本次拟合系数。将拟合系数记入拟合系数数据库，将有效的ΗΙΚ、Η^记入云底高数据库。
[0054](5)进行全视场云底高修正，并计算得出每10分钟的典型云底高度。图3为使用本发明获得的天顶云底高与不使用本发明获得的天顶云底高差异。具体步骤为:
[0055]步骤1:天空均匀性判断。利用大视场的红外辐射数据的标准差判断天空的均匀性，分为均匀天空和不均匀天空两类。
[0056]步骤2:全视场云底高修正。根据步骤I的结果，若为均匀天空，则直接利用拟合系数进行全视场云底高修正，得到全天空云底高分布图；若为不均匀天空，则根据待修正的红外辐射反演的云底高查找云底高数据库，将之前2小时内相差在500m以内的数据对找出，拟合得到修正比例因子，并将云底高乘以修正比例因子进行修正。
[0057]步骤3:典型云底高度的求取。按照低云(0-2500m)、中云(2500_6000m)和高云(6000m以上)三档分别进行统计。将每一档的云底高以10分钟为统计范围，由低到高排序，然后求得前2 %所对应的云底高，作为每一档的典型云底高度。
[0058]本发明未涉及部分均与现有技术相同或可采用现有技术加以实现。虽然本发明已以较佳实施例揭露如上，然其并非用以限定本发明。本发明所属【技术领域】中具有通常知识者，在不脱离本发明的精神和范围内，当可作各种的更动与润饰。因此，本发明的保护范围当视权利要求书所界定者为准。
【权利要求】
1.一种红外与激光联合的地基测云方法，其特征在于，包括以下步骤: (1)分别利用非致冷红外焦平面阵列传感器获得大气向下红外辐射数据，利用激光传感器获得天顶后向消光系数廓线数据；获得上述两个数据的时间同步； (2)联合上述数据对云下水汽和气溶胶辐射进行估计，并利用辐射传输模式计算的晴空阈值进行初步的云检测；并假定云为黑体，反演得到云底高； (3)对高时间分辨率的红外辐射图像进行序列分析，结合晴空阈值做进一步的云检测，并计算云量； (4)拟合红外辐射反演的云底高与激光测量的云底高之间的比例系数； (5)进行全视场云底高修正，并计算得出每10分钟的典型云底高度。
2.根据权利要求1所述的方法，其特征在于:通过数据采集系统每I分钟同步进行非致冷红外焦平面阵列传感器数据采集和激光传感器的数据采集。
3.根据权利要求1所述的方法，其特征在于:联合非致冷红外焦平面阵列传感器获取的红外辐射数据和激光传感器获取的天顶回波数据对云下水汽和气溶胶辐射进行估计，并计算出晴空辐射阈值。
4.根据权利要求1所述的方法，其特征在于:(2)的具体化:根据(I)数据对云下水汽和气溶胶辐射进行估计，并利用辐射传输模式计算的晴空阈值进行初步的云检测；并假定云为黑体，反演得到云底高: 步骤1:云下气溶胶福射估计；将消光廓线转换为光学厚度廓线，代入Libradtran福射传输模式，估算垂直路径的气溶胶向下辐射值I.。-; 步骤2:云下水汽辐射估计:利用激光无云时的天顶向下红外辐射Iik d■和地面水汽密度P water_vapour估算垂直路径的可降水含量PWV，PffV = f(IIE clear, P water_vapour)，其中函数f利用不少于I年的时间的探空廓线、天空红外辐射和地面水汽密度统计得出，具体统计模型可使用二元线性回归模型，并利用最小二乘法得到回归系数；将PWV代入包括Libradtran的辐射传输模式，即可估算垂直路径的水汽向下辐射值IwatCT—vap.; 步骤3:晴空福射阈值确定；利用Libradtran福射传输模式计算大气廓线以及臭氧等其他稳定少变气体的固定辐射值I。，结合步骤1、2得到的I.。-和Iwatw vap■，计算出晴空辐射阈值Itl = laerosol + Iwater—vapour+Ic ? 步骤4:初步的云检测；遍历红外辐射图像，判断实际测量的大气向下红外辐射I与Itl的大小，若实际测量的辐射I大于晴空辐射阈值Itl,则判为有云，否则无云；此步骤得到初步估计云量记为CCl ;并进一步得到修正后的云体向下红外辐射I1 = 1-1aerosol_Iwater_
-Tvapour 丄 c 0 步骤5:利用红外辐射反演云底高；根据步骤4中得到的I1，假设云为黑体，求出云层底亮温，根据该地区该时间段内的典型温度递减率，得到红外辐射反演的云底高HIK。
5.根据权利要求4所述的方法，其特征在于:对高时间分辨率的红外辐射图像进行序列分析，结合晴空阈值做进一步的云检测，并计算云量；具体步骤为: 步骤1:求取各像素点的序列辐射标准差。以I分钟为迭代周期，按10分钟为间隔进行图像的序列分析，将每个像素点的10次辐射数据求标准差σ ;利用辐射标准差判断每个像素随时间的均匀性； 步骤2:标记可疑像素云点。如步骤I求取各像素点的序列辐射标准差，遍历红外辐射标准差图像，得σ，若σ > σ ^，则此像素点标记为可疑云点，否则标记为晴空点，其中O。通过大量晴空样本统计得出。此步骤得到可疑云量记为CC2。 步骤3:对初步云检测进行修正: 对于(2)中判断为无云的区域，若(3)中步骤2得到的是可疑云区，则晴空辐射阈值Itl按0.2ff/ (m2.sr)向下调整，直到CCl与CC2的差值绝对值小于I ; 对于(2)中判断为有云的区域，且根据红外辐射反演的云底高在6000m以上，若(3)中步骤2得到的是无云区，则晴空辐射阈值Itl按0.2ff/(m2.sr)向上调整，直到CCl与CC2的差值绝对值小于I ;对于其他情况不做修正； 步骤4:按照调整后的晴空阈值进行云检测，并得到修正后的云量CC。
6.根据权利要求1所述的方法，其特征在于:拟合红外辐射反演的云底高与激光测量的云底高之间的比例系数:具体方法是，统计10分钟内，天顶5°以内红外辐射反演的云底高Hik和最低层天顶激光云底高的中间值(10组数据)，若两类数据均存在，且大于等于3组，则拟合得到比例系数；若其中任一类不存在，则天顶扩大至10°，统计时间向前扩展至15分钟。若其中任一类仍不存在或者两类数据均存在的情况不足3组数据，则放弃拟合，利用前一次拟合结果作为本次拟合系数。将拟合系数记入拟合系数数据库，将有效的^IE Λ ^Laser
记入云底高数据库。
7.根据权利要求1所述的方法，其特征在于:进行全视场云底高修正，并计算得出每10分钟的典型云底高度:具体步骤为； 步骤1:天空均匀性判断:利用大视场的红外辐射数据的标准差判断天空的均匀性，分为均匀天空和不均匀天空两类。 步骤2:全视场云底高修正:根据步骤I的结果，若为均匀天空，则直接利用拟合系数进行全视场云底高修正，得到全天空云底高分布图；若为不均匀天空，则根据待修正的红外辐射反演的云底高查找云底高数据库，将之前2小时内相差在500m以内的数据对找出，拟合得到修正比例因子，并将云底高乘以修正比例因子进行修正； 步骤3:典型云底高度的求取:按照0-2500m低云、2500-6000m中云和6000m以上的高云三档分别进行统计；将每一档的云底高以10分钟为统计范围，由低到高排序，然后求得前2%所对应的云底高，作为每一档的典型云底高度。
【文档编号】G01W1/00GK104181612SQ201410398747
【公开日】2014年12月3日 申请日期:2014年8月13日 优先权日:2014年8月13日
【发明者】刘磊, 高太长, 张婷, 赵世军, 韩文宇 申请人:中国人民解放军理工大学