《一种红外与激光联合的地基测云方法.pdf》由会员分享，可在线阅读，更多相关《一种红外与激光联合的地基测云方法.pdf（9页完整版）》请在专利查询网上搜索。
1、10申请公布号CN104181612A43申请公布日20141203CN104181612A21申请号201410398747722申请日20140813G01W1/00200601G01S17/95200601G06F17/3020060171申请人中国人民解放军理工大学地址211101江苏省南京市江宁区双龙街60号理工大学72发明人刘磊高太长张婷赵世军韩文宇74专利代理机构南京瑞弘专利商标事务所普通合伙32249代理人陈建和54发明名称一种红外与激光联合的地基测云方法57摘要红外与激光联合的地基测云方法，包括以下步骤1分别利用非致冷红外焦平面阵列传感器获得大气向下红外辐射数据，利用激光传感。
2、器获得天顶后向消光系数廓线数据；获得上述两个数据的时间同步；2联合上述数据对云下水汽和气溶胶辐射进行估计，并利用辐射传输模式计算的晴空阈值进行初步的云检测；并假定云为黑体，反演得到云底高；3对高时间分辨率的红外辐射图像进行序列分析，结合晴空阈值做进一步的云检测，并计算云量；4拟合红外辐射反演的云底高与激光测量的云底高之间的比例系数；5进行全视场云底高修正，并计算得出每10分钟的典型云底高度。51INTCL权利要求书2页说明书5页附图1页19中华人民共和国国家知识产权局12发明专利申请权利要求书2页说明书5页附图1页10申请公布号CN104181612ACN104181612A1/2页21一种红。
3、外与激光联合的地基测云方法，其特征在于，包括以下步骤1分别利用非致冷红外焦平面阵列传感器获得大气向下红外辐射数据，利用激光传感器获得天顶后向消光系数廓线数据；获得上述两个数据的时间同步；2联合上述数据对云下水汽和气溶胶辐射进行估计，并利用辐射传输模式计算的晴空阈值进行初步的云检测；并假定云为黑体，反演得到云底高；3对高时间分辨率的红外辐射图像进行序列分析，结合晴空阈值做进一步的云检测，并计算云量；4拟合红外辐射反演的云底高与激光测量的云底高之间的比例系数；5进行全视场云底高修正，并计算得出每10分钟的典型云底高度。2根据权利要求1所述的方法，其特征在于通过数据采集系统每1分钟同步进行非致冷红外。
4、焦平面阵列传感器数据采集和激光传感器的数据采集。3根据权利要求1所述的方法，其特征在于联合非致冷红外焦平面阵列传感器获取的红外辐射数据和激光传感器获取的天顶回波数据对云下水汽和气溶胶辐射进行估计，并计算出晴空辐射阈值。4根据权利要求1所述的方法，其特征在于2的具体化根据1数据对云下水汽和气溶胶辐射进行估计，并利用辐射传输模式计算的晴空阈值进行初步的云检测；并假定云为黑体，反演得到云底高步骤1云下气溶胶辐射估计；将消光廓线转换为光学厚度廓线，代入LIBRADTRAN辐射传输模式，估算垂直路径的气溶胶向下辐射值IAEROSOL；步骤2云下水汽辐射估计利用激光无云时的天顶向下红外辐射IIR_CLEA。
5、R和地面水汽密度WATER_VAPOUR估算垂直路径的可降水含量PWV，PWVFIIR_CLEAR，WATER_VAPOUR，其中函数F利用不少于1年的时间的探空廓线、天空红外辐射和地面水汽密度统计得出，具体统计模型可使用二元线性回归模型，并利用最小二乘法得到回归系数；将PWV代入包括LIBRADTRAN的辐射传输模式，即可估算垂直路径的水汽向下辐射值IWATER_VAPOUR；步骤3晴空辐射阈值确定；利用LIBRADTRAN辐射传输模式计算大气廓线以及臭氧等其他稳定少变气体的固定辐射值IC，结合步骤1、2得到的IAEROSOL和IWATER_VAPOUR，计算出晴空辐射阈值I0IAEROSO。
6、LIWATER_VAPOURIC；步骤4初步的云检测；遍历红外辐射图像，判断实际测量的大气向下红外辐射I与I0的大小，若实际测量的辐射I大于晴空辐射阈值I0，则判为有云，否则无云；此步骤得到初步估计云量记为CC1；并进一步得到修正后的云体向下红外辐射I1IIAEROSOLIWATER_VAPOURIC。步骤5利用红外辐射反演云底高；根据步骤4中得到的I1，假设云为黑体，求出云层底亮温，根据该地区该时间段内的典型温度递减率，得到红外辐射反演的云底高HIR。5根据权利要求4所述的方法，其特征在于对高时间分辨率的红外辐射图像进行序列分析，结合晴空阈值做进一步的云检测，并计算云量；具体步骤为步骤1求取。
7、各像素点的序列辐射标准差。以1分钟为迭代周期，按10分钟为间隔进行图像的序列分析，将每个像素点的10次辐射数据求标准差；利用辐射标准差判断每个像素随时间的均匀性；步骤2标记可疑像素云点。如步骤1求取各像素点的序列辐射标准差，遍历红外辐射权利要求书CN104181612A2/2页3标准差图像，得，若0，则此像素点标记为可疑云点，否则标记为晴空点，其中0通过大量晴空样本统计得出。此步骤得到可疑云量记为CC2。步骤3对初步云检测进行修正对于2中判断为无云的区域，若3中步骤2得到的是可疑云区，则晴空辐射阈值I0按02W/M2SR向下调整，直到CC1与CC2的差值绝对值小于1；对于2中判断为有云的区域，。
8、且根据红外辐射反演的云底高在6000M以上，若3中步骤2得到的是无云区，则晴空辐射阈值I0按02W/M2SR向上调整，直到CC1与CC2的差值绝对值小于1；对于其他情况不做修正；步骤4按照调整后的晴空阈值进行云检测，并得到修正后的云量CC。6根据权利要求1所述的方法，其特征在于拟合红外辐射反演的云底高与激光测量的云底高之间的比例系数具体方法是，统计10分钟内，天顶5以内红外辐射反演的云底高HIR和最低层天顶激光云底高HLASER的中间值10组数据，若两类数据均存在，且大于等于3组，则拟合得到比例系数；若其中任一类不存在，则天顶扩大至10，统计时间向前扩展至15分钟。若其中任一类仍不存在或者两类。
9、数据均存在的情况不足3组数据，则放弃拟合，利用前一次拟合结果作为本次拟合系数。将拟合系数记入拟合系数数据库，将有效的HIR、HLASER记入云底高数据库。7根据权利要求1所述的方法，其特征在于进行全视场云底高修正，并计算得出每10分钟的典型云底高度具体步骤为；步骤1天空均匀性判断利用大视场的红外辐射数据的标准差判断天空的均匀性，分为均匀天空和不均匀天空两类。步骤2全视场云底高修正根据步骤1的结果，若为均匀天空，则直接利用拟合系数进行全视场云底高修正，得到全天空云底高分布图；若为不均匀天空，则根据待修正的红外辐射反演的云底高查找云底高数据库，将之前2小时内相差在500M以内的数据对找出，拟合得到。
10、修正比例因子，并将云底高乘以修正比例因子进行修正；步骤3典型云底高度的求取按照02500M低云、25006000M中云和6000M以上的高云三档分别进行统计；将每一档的云底高以10分钟为统计范围，由低到高排序，然后求得前2所对应的云底高，作为每一档的典型云底高度。权利要求书CN104181612A1/5页4一种红外与激光联合的地基测云方法技术领域0001本发明涉及一种地基测云方法，尤其是一种联合地基红外与激光进行云检测和全视场云底高修正的方法，其主要用于地面气象观测自动化测量中。背景技术0002云通常覆盖地球大气约50左右，是地气系统辐射收支的主要调节者。云的观测工作一直是气候研究、天气分析与。
11、预报工作的基础。长期以来，目测是气象业务中常用的估计云量、确定云高的方法，但由于人为观测误差大，加之夜间光照条件差，云的观测资料的准确性和昼夜一致性难以保证，同时对实际天空的描述也较为困难。目前主要的地基测云技术体制有可见光被动遥感体制、红外被动遥感体制和激光主动遥感体制。可见光测云仪由于其夜间无法测量，很难应用于实际观测业务中。红外测云仪和激光云高仪能昼夜连续观测且各有特点。0003地基红外被动遥感测云技术国外以美国ICI测云仪研究团队为代表，国内以中科院大气物理所和解放军理工大学气象海洋学院为代表，均取得大量创新性成果。研究表明，814M波段向下红外辐射主要受到水汽和气溶胶的影响，单一红外。
12、测云仪由于对云下水汽和气溶胶的估计不足会造成云检测出现偏差。另外，红外辐射反演的云底高只适用近似为黑体的云，其他情况下由于无法得知云的光学厚度，反演的云底高存在很大误差。0004地基激光主动遥感测云技术，国外以VAISALA公司为代表，其研制的激光云高仪在20世纪90年代应用于美国地面自动观测系统ASOS中替代人工目测。国内也有多家公司研发出了激光云高仪。一些学者研究发现，激光云高仪会在一些明显无云的情况下出现错误的有云信号，这可能是由于路径上存在一个很薄的湿气层或气溶胶层造成的。此外噪声的影响使其测量卷云这一弱回波信号时的数据获取率极低，这一问题在我国雾霾多发地体现的尤为明显。更重要的是，激。
13、光云高仪采用时间积分计算的云量存在严重的空间代表性问题。0005综上所述，红外被动遥感测云和激光主动遥感测云这两种最有可能实现测云业务化的技术体制目前仍然在不同程度上存在问题，尚不能完全替代人工目测。发明内容0006本发明的目的是提出一种将地基红外被动遥感和激光主动遥感联合进行云检测和全视场云底高修正的方法，从而最大程度地克服地基红外被动遥感测云和激光主动遥感测云的使用场合的不足，提高地基昼夜连续测云的准确性。0007本发明的技术方案为一种红外与激光联合的地基测云方法。该方法包括以下步骤00081分别利用非致冷红外焦平面阵列传感器获得大气向下红外辐射数据，利用激光传感器获得天顶后向消光系数廓线。
14、数据；获取数据的时间同步；如通过数据采集系统每1分钟同步进行数据采集，并将数据存储。说明书CN104181612A2/5页500092联合上述数据对云下水汽辐射和气溶胶辐射进行估计，并利用辐射传输模式计算的晴空阈值进行初步的云检测；并假定云为黑体，反演得到云底高；00103对高时间分辨率的红外辐射图像进行序列分析，结合晴空阈值做进一步的云检测，并计算云量；00114拟合红外辐射反演的云底高与激光测量的云底高之间的比例系数；00125进行全视场云底高修正，并计算得出每10分钟的典型云底高度。0013进一步，通过数据采集系统每1分钟同步进行非致冷红外焦平面阵列传感器数据采集和激光传感器的数据采集。。
15、0014进一步，联合非致冷红外焦平面阵列传感器获取的红外辐射数据和激光传感器获取的天顶回波数据对云下水汽和气溶胶辐射进行估计，并计算出晴空辐射阈值。00152的具体化联合上述数据对云下水汽和气溶胶辐射进行估计，并利用辐射传输模式计算的晴空阈值进行初步的云检测；并假定云为黑体，反演得到云底高0016步骤1云下气溶胶辐射估计；将消光廓线转换为光学厚度廓线，代入LIBRADTRAN辐射传输模式，估算垂直路径的气溶胶向下辐射值IAEROSOL；0017步骤2云下水汽辐射估计利用激光无云时的天顶向下红外辐射IIR_CLEAR和地面水汽密度WATER_VAPOUR估算垂直路径的可降水含量PWVPWVFII。
16、R_CLEAR，WATER_VAPOUR，其中函数F利用一段时间例如不少于1年的时间的探空廓线、天空红外辐射和地面水汽密度统计得出，具体统计模型可使用二元线性回归模型，并利用最小二乘法得到回归系数。将PWV代入LIBRADTRAN辐射传输模式但不限于LIBRADTRAN，即可估算垂直路径的水汽向下辐射值IWATER_VAPOUR。0018步骤3晴空辐射阈值确定；利用LIBRADTRAN辐射传输模式计算大气廓线以及臭氧等其他稳定少变气体的固定辐射值IC，结合步骤1、2得到的IAEROSOL和IWATER_VAPOUR，计算出晴空辐射阈值I0IAEROSOLIWATER_VAPOURIC；0019。
17、步骤4初步的云检测；遍历红外辐射图像，判断实际测量的大气向下红外辐射I与I0的大小，若实际测量的辐射I大于晴空辐射阈值I0，则判为有云，否则无云；此步骤得到初步估计云量记为CC1；并进一步得到修正后的云体向下红外辐射I1IIAEROSOLIWATER_VAPOURIC。0020步骤5利用红外辐射反演云底高；根据步骤4中得到的I1，假设云为黑体，求出云层底亮温，根据该地区该时间段内的典型温度递减率，得到红外辐射反演的云底高HIR。0021对高时间分辨率的红外辐射图像进行序列分析，结合晴空阈值做进一步的云检测，并计算云量；具体步骤为0022步骤1求取各像素点的序列辐射标准差。以1分钟为迭代周期，按。
18、10分钟为间隔进行图像的序列分析，将每个像素点的10次辐射数据求标准差；利用辐射标准差判断每个像素随时间的均匀性；0023步骤2标记可疑像素云点。如步骤1求取各像素点的序列辐射标准差，遍历红外辐射标准差图像，得，若0，则此像素点标记为可疑云点，否则标记为晴空点，其中0通过大量晴空样本统计得出。此步骤得到可疑云量记为CC2。0024步骤3对初步云检测进行修正0025对于2中判断为无云的区域，若3中步骤2得到的是可疑云区，则晴空辐射阈说明书CN104181612A3/5页6值I0按02W/M2SR向下调整，直到CC1与CC2的差值绝对值小于1；0026对于2中判断为有云的区域，且根据红外辐射反演的。
19、云底高在6000M以上，若3中步骤2得到的是无云区，则晴空辐射阈值I0按02W/M2SR向上调整，直到CC1与CC2的差值绝对值小于1；对于其他情况不做修正。0027通过此步骤，进一步提高了对高云的检测能力。0028步骤4按照调整后的晴空阈值进行云检测，并得到修正后的云量CC。0029拟合红外辐射反演的云底高与激光测量的云底高之间的比例系数具体方法是，统计10分钟内，天顶5以内红外辐射反演的云底高HIR和最低层天顶激光云底高HLASER的中间值10组数据，若两类数据均存在，且大于等于3组，则拟合得到比例系数；若其中任一类不存在，则天顶扩大至10，统计时间向前扩展至15分钟。若其中任一类仍不存在。
20、或者两类数据均存在的情况不足3组数据，则放弃拟合，利用前一次拟合结果作为本次拟合系数。将拟合系数记入拟合系数数据库，将有效的HIR、HLASER记入云底高数据库。0030进行全视场云底高修正，并计算得出每10分钟的典型云底高度具体步骤为；0031步骤1天空均匀性判断。利用大视场的红外辐射数据的标准差判断天空的均匀性，分为均匀天空和不均匀天空两类。0032步骤2全视场云底高修正。根据步骤1的结果，若为均匀天空，则直接利用拟合系数进行全视场云底高修正，得到全天空云底高分布图；若为不均匀天空，则根据待修正的红外辐射反演的云底高查找云底高数据库，将之前2小时内相差在500M以内的数据对找出，拟合得到修。
21、正比例因子，并将云底高乘以修正比例因子进行修正。0033步骤3典型云底高度的求取。按照低云02500M、中云25006000M和高云6000M以上三档分别进行统计。将每一档的云底高以10分钟为统计范围，由低到高排序，然后求得前2所对应的云底高，作为每一档的典型云底高度。0034与现有技术相比，本发明的有益效果在于本发明公开了一种将地基红外被动遥感和激光主动遥感联合进行云检测和全视场云底高修正的方法。0035本能更加准确地估计云下气溶胶和水汽的向下红外辐射，这将使得辐射传输模式计算的晴空阈值更加接近真实值，从而提高了云检测效果。进一步地，利用高时间分辨率的序列分析可以提高对卷云的检测能力。结合激。
22、光云底高可以修正红外反演的云底高，不但能得到更加准确的全天空云底高分布，还可统计得出全天空典型云底高度信息。本发明通过充分联合地基测量的大气向下红外辐射数据、天顶激光回波数据、高时间分辨率的红外辐射序列图像，采用综合观测的理念，为新一代自动气象站中云的自动观测提供一种联合探测的技术途径，较好地克服了单一测云手段的不足。附图说明0036图1红外与激光联合的地基测云方法流程图；0037图2本发明得到的云量与只使用红外测云仪的云量结果对比；0038图3本发明得到的修正后天顶云底高与只使用红外测云仪的天顶云底高对比。具体实施方式0039为了使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例。
23、，并参说明书CN104181612A4/5页7照附图，对本发明进一步详细说明。0040图1是本发明提出的红外与激光联合的地基测云方法流程图，所述方法的具体实施方式如下00411分别利用非致冷红外焦平面阵列传感器获得大气向下红外辐射数据，利用激光传感器获得天顶后向消光系数廓线数据。这一步骤中，关键在于获取数据的时间同步，通过数据采集系统每1分钟同步进行数据采集，并将数据存储。00422联合上述数据对云下水汽和气溶胶辐射进行估计，并利用辐射传输模式计算的晴空阈值进行初步的云检测；并假定云为黑体，反演得到云底高。具体步骤有0043步骤1云下气溶胶辐射估计。将消光廓线转换为光学厚度廓线，代入LIBRA。
24、DTRAN辐射传输模式，估算垂直路径的气溶胶向下辐射IAEROSOL。0044步骤2云下水汽辐射估计利用激光无云时的天顶向下红外辐射IIR_CLEAR和地面水汽密度WATER_VAPOUR估算垂直路径的可降水含量PWVPWVFIIR_CLEAR，WATER_VAPOUR，其中函数F利用一段时间例如不少于1年的时间的探空廓线、天空红外辐射和地面水汽密度统计得出，具体统计模型可使用二元线性回归模型，并利用最小二乘法得到回归系数。将PWV代入LIBRADTRAN辐射传输模式但不限于LIBRADTRAN，即可估算垂直路径的水汽向下辐射值IWATER_VAPOUR。0045步骤3晴空辐射阈值确定。利用L。
25、IBRADTRAN辐射传输模式计算大气廓线以及臭氧等其他稳定少变气体的固定辐射值IC，结合步骤1、2得到的IAEROSOL和IWATER_VAPOUR，计算出晴空辐射阈值I0IAEROSOLIWATER_VAPOURIC。0046步骤4初步的云检测。遍历红外辐射图像，判断实际测量的大气向下红外辐射I与I0的大小，若实际测量的辐射I大于晴空辐射阈值I0，则判为有云，否则无云。此步骤得到初步估计云量记为CC1。并进一步得到修正后的云体向下红外辐射I1IIAEROSOLIWATER_VAPOURIC。0047步骤5利用红外辐射反演云底高。根据步骤4中得到的I1，假设云为黑体，求出云层底亮温，根据该地。
26、区该时间段内的典型温度递减率，得到红外辐射反演的云底高HIR。此步骤由于大多数利用红外辐射反演云底高均采用此方法，不做详细介绍。00483对高时间分辨率的红外辐射图像进行序列分析，结合晴空阈值做进一步的云检测，并计算云量。具体步骤为0049步骤1求取各像素点的序列辐射标准差。以1分钟为迭代周期，按10分钟为间隔进行图像的序列分析，将每个像素点的10次辐射数据求标准差。利用辐射标准差判断每个像素随时间的均匀性。0050步骤2标记可疑云点。遍历红外辐射标准差图像，若0，则此像素点标记为可疑云点，否则标记为晴空点，其中0通过大量晴空样本统计得出。此步骤得到可疑云量记为CC2。0051步骤3对初步云检。
27、测进行修正。对于2中判断为无云的区域，若3中步骤2得到的是可疑云区，则晴空辐射阈值I0按02W/M2SR向下调整，直到CC1与CC2的差值绝对值小于1。对于2中判断为有云的区域，且根据红外辐射反演的云底高在6000M以上，若3中步骤2得到的是无云区，则晴空辐射阈值I0按02W/M2SR向上调整，直到CC1与CC2的差值绝对值小于1。对于其他情况不做修正。通过此步骤，进一步提高了对高说明书CN104181612A5/5页8云的检测能力。0052步骤4按照调整后的晴空阈值进行云检测，并得到修正后的云量CC。图2为使用本发明获得的云量与不使用本发明获得的云量差异。00534拟合红外辐射反演的云底高与。
28、激光测量的云底高之间的比例系数。具体方法是统计10分钟内，天顶5以内红外辐射反演的云底高HIR和最低层天顶激光云底高HLASER的中间值10组数据，若两类数据均存在，且大于等于3组，则拟合得到比例系数；若其中任一类不存在，则天顶扩大至10，统计时间向前扩展至15分钟。若其中任一类仍不存在或者两类数据均存在的情况不足3组数据，则放弃拟合，利用前一次拟合结果作为本次拟合系数。将拟合系数记入拟合系数数据库，将有效的HIR、HLASER记入云底高数据库。00545进行全视场云底高修正，并计算得出每10分钟的典型云底高度。图3为使用本发明获得的天顶云底高与不使用本发明获得的天顶云底高差异。具体步骤为00。
29、55步骤1天空均匀性判断。利用大视场的红外辐射数据的标准差判断天空的均匀性，分为均匀天空和不均匀天空两类。0056步骤2全视场云底高修正。根据步骤1的结果，若为均匀天空，则直接利用拟合系数进行全视场云底高修正，得到全天空云底高分布图；若为不均匀天空，则根据待修正的红外辐射反演的云底高查找云底高数据库，将之前2小时内相差在500M以内的数据对找出，拟合得到修正比例因子，并将云底高乘以修正比例因子进行修正。0057步骤3典型云底高度的求取。按照低云02500M、中云25006000M和高云6000M以上三档分别进行统计。将每一档的云底高以10分钟为统计范围，由低到高排序，然后求得前2所对应的云底高，作为每一档的典型云底高度。0058本发明未涉及部分均与现有技术相同或可采用现有技术加以实现。虽然本发明已以较佳实施例揭露如上，然其并非用以限定本发明。本发明所属技术领域中具有通常知识者，在不脱离本发明的精神和范围内，当可作各种的更动与润饰。因此，本发明的保护范围当视权利要求书所界定者为准。说明书CN104181612A1/1页9图1图2图3说明书附图CN104181612A。