技术领域本发明涉及一种自动确定天气情况的方法，尤其涉及一种自动确定全天无云和部分有云天气情况的方法。
背景技术：
蒸散发是水循环与地表能量平衡的重要组成部分，在水文学、气象学、农学、地学等学科研究中，都发挥着重要的作用。遥感技术从技术性、经济性和实效性等方面来看，被认为是获得地球表面区域尺度蒸散发分布唯一可行的方法。研究表明，局地、区域和全球尺度的蒸散发分布可以依据遥感反演得到的地表参数和地面观测数据进行估算。但是遥感获取蒸散发所需地面参数大多数是在卫星过境时刻为晴天的情况下进行。目前大量的对遥感反演地表蒸散发的研究，都存在对地表能量辐射各组分、水汽压、阻抗、潜在蒸发、参考蒸发等蒸散发的控制/表征因子受云影响的日变化过程及相互作用机理的认识严重不足的问题，使得蒸散发遥感估算精度受到影响。因此，基于地面入射太阳辐射和短波大气透射率数据的日内变化特征，开展全天无云、部分有云等天气情况的判断，研究蒸散发各控制/表征因子受云影响的情况，对于提高遥感反演地表蒸散发精度具有重要意义。
技术实现要素：
为了解决上述技术所存在的不足之处，本发明提供了一种自动确定全天无云和部分有云天气情况的方法。为了解决以上技术问题，本发明采用的技术方案是：一种自动确定全天无云和部分有云天气情况的方法，包括以下步骤：步骤一：数据获取及其预处理，包括以下步骤：(1)首先获得半小时尺度地面入射太阳短波辐射数据，其为地面直接测量数据，记为RS；(2)计算对应半小时时间尺度的短波大气透射率数据；短波大气透射率为入射地面太阳短波辐射数据与大气层顶太阳辐射数据的比率，记为τ,通过公式Ⅰ计算得出：τ=RsRp---I]]>其中，RS为入射地面的太阳短波辐射；Rp为大气层顶太阳辐射；大气层顶太阳辐射Rp通过公式Ⅱ计算得出：Rp＝S0cosθⅡ其中，S0为太阳常数，取1365Wm-2；θ为太阳天顶角；步骤二：确定一天中的日出日落时刻：日出时刻为一天中首次地面入射太阳短波辐射大于5W/m2的时刻，记为tr；日落时刻为一天结束前最后一次地面入射太阳短波辐射大于5W/m2的时刻，记为ts；步骤三：全天无云天气的确定：在给定的一天中，如果半小时尺度地面入射太阳辐射RS和对应的短波大气透射率数据τ从日出时刻tr到正午时刻tm单调增加，同时满足从正午时刻tm到日落时刻ts单调减少，则将这一天选中为晴天，即晴天条件满足公式III：其中，RS为地面入射短波太阳辐射；ti为第i个半小时对应的时刻；ti+1为第i+1个半小时对应的时刻；tr为日出时刻；tm为正午时刻；ts为日落时刻；然后结合遥感中分辨率成像光谱仪MODIS(moderate-resolutionimagespectroradiometer)云掩膜产品MOD35的质量控制数据，即质量标记首两位值为11(可信度高的晴空)的数据，从所选晴天中排除持续有云天气，则所剩晴天为全天无云天；步骤四：部分有云天气的确定，具体步骤为：(1)在给定的一天中，如果半小时尺度地面入射太阳短波辐射RS和相应的短波大气透射率数据τ从日出时刻tr到正午时刻tm不是单调增加或者从正午时刻tm到日落时刻ts不是单调减小，则这一天为部分有云天；(2)在部分有云天中继续选取云出现间隔有晴天出现的部分；先从白天时间选择9:00-10:00、10:00-11:00、11:00-12:00、12:00-13:00、13:00-14:00、14:00-15:00的时间间隔；由于短波大气透射率τ随着太阳天顶角和总体积含水量变化而变化，为了判断部分有云天各时间间隔中是否有晴天出现，首先将τ按照公式IV转换为垂直方向的值：τ0＝τcosθIV其中，τ0为垂直方向的短波大气透射率；θ为太阳天顶角；若部分有云天给定时间间隔的垂直方向短波大气透射率τ0大于两个相邻晴天对应时间间隔的τ0，即满足公式V，则该部分有云天的该时间段有晴天出现，被选择用于进一步分析；否则，该部分有云天被舍弃；τ0i＞τ0i(before)&τ0i＞τ0i(after)V其中，τ0i为该有云天给定时间间隔垂直方向的短波大气透射率；τ0i(before)为该有云天前一相邻晴天对应时间间隔垂直方向的短波大气透射率，τ0i(after)为该有云天后一相邻晴天对应时间间隔垂直方向的短波大气透射率。根据以上步骤，可快速得到完全晴天和部分有云天的天气判断。本发明与现有技术相比，仅通过入射地面太阳辐射和短波大气透射率在一天中的变化情况，结合MODIS云掩膜产品即可进行全天无云和部分有云天气的判断，避免了获取大量天气状况数据的问题。本发明原理简单易操作、确定精度高，能快速得到全天无云和部分有云天天气情况，并得到一天中各时间段间隔有晴天的部分有云天。具体实施方式下面结合具体实施方式对本发明作进一步详细的说明。本发明包括以下步骤：步骤一：数据获取及其预处理，包括以下步骤：(1)首先获得半小时尺度地面入射太阳短波辐射数据，其为地面直接测量数据，记为RS；(2)计算对应半小时时间尺度的短波大气透射率数据；短波大气透射率为入射地面太阳短波辐射数据与大气层顶太阳辐射数据的比率，记为τ,通过公式Ⅰ计算得出：τ=RsRp---I]]>其中，RS为入射地面的太阳短波辐射；Rp为大气层顶太阳辐射；大气层顶太阳辐射Rp通过公式Ⅱ计算得出：Rp＝S0cosθⅡ其中，S0为太阳常数，取1365Wm-2；θ为太阳天顶角；步骤二：确定一天中的日出日落时刻：日出时刻为一天中首次地面入射太阳短波辐射大于5W/m2的时刻，记为tr；日落时刻为一天结束前最后一次地面入射太阳短波辐射大于5W/m2的时刻，记为ts；步骤三：全天无云天气的确定：在给定的一天中，如果半小时尺度地面入射太阳辐射RS和对应的短波大气透射率数据τ从日出时刻tr到正午时刻tm单调增加，同时从正午时刻tm到日落时刻ts单调减少，则将这一天选中为晴天，即晴天条件满足公式III：其中，RS为地面入射短波太阳辐射；ti为第i个半小时对应的时刻；ti+1为第i+1个半小时对应的时刻；tr为日出时刻；tm为正午时刻；ts为日落时刻；然后结合遥感中分辨率成像光谱仪MODIS(moderate-resolutionimagespectroradiometer)云掩膜产品MOD35的质量控制数据，即质量标记首两位值为11(可信度高的晴空)的数据，从所选晴天中排除持续有云天气，则所剩晴天为全天无云天；步骤四：部分有云天气的确定，具体步骤为：(1)在给定的一天中，如果半小时尺度地面入射太阳短波辐射RS和相应的短波大气透射率数据τ从日出时刻tr到正午时刻tm不是单调增加或者从正午时刻tm到日落时刻ts不是单调减小，则这一天为部分有云天；(2)在部分有云天中继续选取云出现间隔有晴天出现的部分；先从白天时间选择9:00-10:00、10:00-11:00、11:00-12:00、12:00-13:00、13:00-14:00、14:00-15:00的时间间隔；由于短波大气透射率τ随着太阳天顶角和总体积含水量变化而变化，为了判断部分有云天各时间间隔中是否有晴天出现，首先将τ按照公式IV转换为垂直方向的值：τ0＝τcosθIV其中，τ0为垂直方向的短波大气透射率；θ为太阳天顶角；若部分有云天给定时间间隔的垂直方向短波大气透射率τ0大于两个相邻晴天对应时间间隔的τ0，即满足公式V，则该部分有云天的该时间段有晴天出现，被选择用于进一步分析；否则，该部分有云天被舍弃；τ0i＞τ0i(before)&τ0i＞τ0i(after)V其中，τ0i为该有云天给定时间间隔垂直方向的短波大气透射率；τ0i(before)为该有云天前一相邻晴天对应时间间隔垂直方向的短波大气透射率，τ0i(after)为该有云天后一相邻晴天对应时间间隔垂直方向的短波大气透射率。根据以上步骤，可快速得到完全晴天和部分有云天的天气判断。本发明充分利用给定一天中半小时尺度的地面入射太阳短波辐射数据和相应短波大气透射率数据在一天的变化情况，结合遥感MODIS云掩膜产品得到完全晴天和部分有云天的天气情况判断，同时对部分有云天各时间段短波大气透射率数据进行进一步分析，得出各时间段有晴天出现的部分有云天。本发明只需地面入射太阳辐射数据，即可确定日出日落时刻；利用半小时尺度地面入射太阳辐射数据、短波大气透射率数据，结合遥感MODIS云掩膜产品，即可得到全天无云和部分有云天的天气判断。本发明能简单快速地确定完全晴天和部分有云天气状况，避免了需要大量收集天气状况数据的问题，是研究蒸散发各控制/表征因子受云影响情况的前提，对于提高遥感反演地表蒸散发精度具有重要作用。上述实施方式并非是对本发明的限制，本发明也并不仅限于上述举例，本技术领域的技术人员在本发明的技术方案范围内所做出的变化、改型、添加或替换，也均属于本发明的保护范围。