1.本发明涉及利用地面测量的下行热辐射亮度估算大气水汽含量的技术方法和系统，属于定量遥感参数反演的领域。本发明所反演的大气水汽含量具有较高的精度，可以应用于红外遥感数据的大气订正等。
背景技术：
：：2.大气水汽含量(precipitablewatervapor,pwv)是指从地面到大气顶界的单位面积垂直大气柱中所含水汽全部凝结并降落到地面可产生的降水量，又被称为大气可降水量。大气水汽是引起温室效应的主要气体之一，作为能量循环和水循环的重要因素，大气水汽含量是目前气候变化和数值模式等应用主要监测的因子之一。3.当前常规的大气水汽含量监测手段包括气象遥感卫星和雷达数据的遥感反演、无线电探空监测、地基水汽辐射计或gps监测和利用太阳光度计观测的太阳直射辐射进行估算等方法。其中太阳光度计法利用位于940nm的水汽吸收通道进行大气水汽含量估算，是目前应用最广泛的大气水汽含量的获取方法。但现有方法受限于仪器的便携性和反演算法精度，难以快速、准确地获取大气水汽含量。技术实现要素：4.本发明要解决的技术问题是：如何利用便携式地基热辐射计，便捷地实现大气水汽的精确估算。5.目前在地面获取精确大气水汽含量的方法主要基于太阳分光光度计，该方法虽然可以精确获取大气水汽含量，但是由于太阳分光光度计正常工作需要复杂的参数设置和相对稳定的运行环境，且其机动性差，难以快速、准确地获取大气水汽含量。针对现有技术的不足，本发明公开一种基于地基热辐射计测量的多角度下行辐射亮度便捷地间接估算大气水汽含量的新方法和系统。6.本发明创新点在于：(1)首次提出利用地基红外辐射计多角度观测数据估算大气水汽含量。(2)提供一种便捷的大气水汽地基测量方法，提高野外综合遥感实验中大气水汽含量的测量效率和测量精度。7.本发明采用的技术方案如下：8.一种基于地面红外辐射计的大气水汽含量估算方法，包括以下步骤：9.基于大气辐射传输软件模拟不同大气条件和不同角度情况下大气下行热辐射比值与大气水汽含量的关系，建立利用大气下行热辐射比值估算大气水汽含量的回归模型；10.利用热辐射计测量不同角度下的大气下行热辐射，并计算不同角度下的大气下行热辐射比值；11.将计算得到的不同角度下的大气下行热辐射比值输入所述大气下行热辐射比值与大气水汽含量的回归模型中，获得大气水汽含量估算值。12.进一步地，所述建立利用大气下行热辐射比值估算大气水汽含量的回归模型，包括：13.对热红外辐射传输过程进行分析，推导大气下行辐射与大气水汽的关系，基于理论推导的不同角度的两组大气下行热辐射的比值与大气水汽含量关系，使用模拟的不同大气条件、不同观测角度的组合情况下的大气下行热辐射和大气水汽含量数据集，建立估算大气水汽含量的回归模型查找表。14.进一步地，所述建立利用大气下行热辐射比值估算大气水汽含量的回归模型，包括：通过大气辐射传输软件(如中等光谱分辨率大气透过率及辐射传输算法和计算模型modtran等)模拟不同大气条件下两个独立观测角度(为简化计算，可以设置其中一个观测天顶角为垂直向上0°观测)的下行热辐射与水汽的关系，构建在不同大气条件下两组不同角度下行热辐射的比值与对应大气水汽含量的回归模型。该回归模型设置的三种大气条件包括热带大气、中纬度大气、极地大气等，在设定的五个高斯角度下建立不同角度下行热辐射亮度比值与大气水汽含量的回归模型，形成不同大气条件和不同角度组合的条件下下行热辐射比值与大气水汽的估算关系。15.进一步地，所述利用热辐射计测量不同角度下的大气下行热辐射，包括：利用经黑体校正的便携式热辐射计测量不同角度下的大气下行辐射，记录测量的下行辐射亮度及对应的测量角度、测量地点的纬度和测量日期(用以评估大气条件)。为提高大气水汽的估算精度，可以测量超过三组的不同角度的大气下行热辐射亮度，取大气水汽含量的估算均值来减少随机误差。16.进一步地，基于测量数据和不同大气、角度组合条件下回归模型估算大气水汽含量，具体包括：17.根据使用便携式热辐射计测量的多组大气下行辐射和辐射信息，计算各组数据在不同角度下的大气下行辐射比值，并基于不同大气和不同角度的组合关系，根据测量的角度和大气状况选择对应的非线性回归模型来估算观测时刻对应的大气水汽含量估算值。多组数据取算术平均以降低随机误差：[0018][0019]其中，wi表示第i组大气水汽含量的估算值，w表示n组大气水汽含量的估算值的均值。[0020]一种基于地面红外辐射计的大气水汽含量估算系统，其包括：[0021]估算模型建立模块，用于基于大气辐射传输软件模拟不同大气条件和不同角度情况下大气下行热辐射比值与大气水汽含量的关系，建立利用大气下行热辐射比值估算大气水汽含量的回归模型；[0022]测量模块，用于利用热辐射计测量不同角度下的大气下行热辐射，并计算不同角度下的大气下行热辐射比值；[0023]大气水汽含量估算模块，用于将计算得到的不同角度下的大气下行热辐射比值输入所述大气下行热辐射比值与大气水汽含量的回归模型中，获得大气水汽含量估算值。[0024]本发明的有益效果如下：[0025]本发明基于地面下行热辐射建模的结果，发展了一种新型的基于地基热辐射计的大气水汽含量的估算方法，从而为卫星、航空红外遥感的大气纠正提供大气水汽含量参数，同时也为通过卫星、航空红外遥感技术反演的大气水汽含量产品的地面验证提供数据支撑。附图说明[0026]图1为本发明方法的流程图。其中vza表示观测方位角度。[0027]图2为不同角度下行辐射测量示意图。其中，表示垂直观测角度下红外辐射计测量的辐射亮度，表示观测天顶角为θ时测量的辐射亮度。具体实施方式[0028]为使本发明的上述目的、特征和优点能够更加明显易懂，下面通过具体实施例和附图，对本发明做进一步详细说明：[0029]步骤一、构建大气水汽估算回归模型。[0030](1)大气水汽含量与多角度下行辐射亮度理论建模[0031]根据热红外辐射传输理论，在晴朗无云的大气条件下，因大气的热辐射作用到达地表的天顶角为θ时的下行热辐射latm,↓(θ)为：[0032][0033]其中θ为天顶角，z为大气层高度，ta为大气等效温度，τa(z,θ)为高度z到地面整层大气的透过率，b为普朗克黑体辐射函数。[0034]假定整个大气层的平均温度为利用积分中值定理，上式可以修改为：[0035][0036]根据beer–lambert定律，高度z到地面大气的透过率τa(θ,z)满足指数关系：[0037]τa(θ,z)＝exp[-γa(z)δ(θ)]ꢀꢀꢀ(4)[0038]其中γa(z)为消光系数，与高度γa(z)下大气的性质有关，而δ(θ)为光程，满足关系δ(θ)＝δ(0)/cosθ。所以有：[0039]τa(θ)＝τa(0)sec(θ)ꢀꢀꢀ(5)[0040]其中，τa(θ)表示观测天顶角为θ时的大气透过率，τa(0)表示垂直观测时的大气透过率。[0041]由于大气水汽在热红外吸收波段的主导作用，τa(θ,z)可以分为大气分子透过率τm(θ,z)和大气水汽透过率τw(θ,z)两部分，满足：[0042]τa(θ,z)＝τm(θ,z)τw(θ,z)ꢀꢀꢀ(6)[0043]由于大气分子吸收作用比较稳定，τm(θ,z)可以认为变化较小，为常系数；大气水汽垂直透过率τw(0,z)则与大气水汽含量在水汽吸收谱段近似满足平方根关系：[0044][0045]其中a，b为常数，m为分子质量，而w为大气水汽含量。[0046]考虑到大气平均温度的变化在小角度条件下可以忽略，即[0047][0048]此时，大气下行热辐射的比值与大气水汽含量的透过率满足：[0049][0050]上述超越方程无法直接求解，可以通过构建透过率、角度与大气下行热辐射比值三者的查找表求解得到大气垂直透过率。此时有：[0051][0052]其中a，b，c表示常数系数，由观测波段的吸收特性决定。[0053]两边取对数可以求得：[0054][0055]改写为：[0056]w＝a(logτa(0,z))2+blogτa(0,z)+cꢀꢀꢀ(12)[0057]其中a,b,c为常数系数，可以通过辐射传输方程模拟回归得到。根据上式即可求解得到大气水汽含量估算值w。[0058](2)大气水汽含量回归模型构建[0059]由于不同大气条件和不同角度下，大气水汽含量及其对下行辐射的吸收作用存在差异，所以需要针对不同大气条件和不同角度条件组合情况下构建利用下行辐射亮度比值估算大气水汽含量的回归模型。为简化模型构建复杂度，可以固定一组测量中其中一个角度为垂直向上0°，另一个角度选择对下行辐射具有代表性的五个高斯角度，分别为向上观测天顶角度为11.6°、26.1°、40.3°、53.7°和65°。[0060]在不同的高斯观测角度下，选择热力学初始分析资料(thermodynamicinitialguessretrieval,tigr)大气廓线数据库中的946条晴空廓线(包括热带廓线236条，中纬度廓线258条，亚极地廓线452条，大气水汽含量范围为0.1-6.3cm)作为大气辐射传输软件modtran的输入，得到不同角度、不同大气组合条件下的大气下行热辐射亮度。同时计算946条晴空廓线的大气水汽含量作为真值，构建不同大气条件和不同角度条件组合情况下的模型查找表，见表1。[0061]表1.不同大气条件和角度条件组合形成的大气水汽含量估算模型查找表[0062][0063]通过上述模型查找表，结合测量时刻的大气状况和观测角度，可以选择更精确的大气水汽估算模型。[0064]为说明大气水汽含量回归模型的估算误差，表2中列举了在53°角度下的大气水汽回归模型在不同大气水汽条件的估算误差。表2说明大气水汽含量回归模型在低水汽效果较好，高水汽条件下较差，但总体误差仍小于0.3cm。[0065]表2.不同大气水汽条件下大气水汽含量的估算误差[0066]水汽含量(cm)决定系数(r2)rmse(cm)[0.0,2.5]0.980.089[2.0,3.5]0.870.16[3.0,4.5]0.730.21[4.0,5.5]0.560.26[5.0,6.3]0.580.24[0.0,4]0.980.125[2.5,6.3]0.900.28[0.0,6.3]0.990.14[0067]步骤二、不同角度下行热辐射亮度测量。[0068]需要利用便携式热辐射计精确测量不同角度下大气的下行热辐射亮度及其辅助信息以减小测量误差。根据测量环境的气温确保测量过程中热辐射计处于仪器标准测量环境温度范围内，在测量前需要使用标准黑体对仪器进行绝对辐射定标以提高测量精度；测量场地需要保证空旷、相对平坦且周围无明显遮挡，且天空保持为晴空状态；在测量过程中需要确定仪器视场内无测量人员遮挡或干扰。使用不同角度测量时可以使用经调整水平后的观测架预先设定仪器观测角度或下行辐射亮度测量后使用角度测量工具测量观测角度，注意测量角度不宜偏大，应小于模型查找表表1中角度的最大值65°，以避免大角度情况下的斜路径大气水汽增大带来的不确定性。一组(一组包括一个垂直向上0°测量值和一个倾斜测量的测量值，见图2)测量结束后，记录倾斜测量的角度、测量的纬度、时间等辅助信息，且确保一组数据中两个不同角度的下行辐射测量日期在三分钟内，其测量过程测量视场内大气无明显扰动，否则该组数据作废。在大气比较稳定的情况下，可以连续测量多组有效数据取均值降低随机误差。[0069]步骤三、选择大气水汽估算回归模型估算大气水汽含量。[0070](1)测量时大气廓线类型的确定[0071]步骤二中测量记录的测量点纬度和测量日期等信息可以辅助确定测量的大气廓线类型，进而与观测角度组合选择最适合的大气水汽含量估算回归模型。使用季节-纬度表面温度模型(seasonal-latitudesurfacetemperaturemodel)，见表3，结合观测纬度和观测日期选择表中的类型作为测量时的大气廓线类型。[0072]表3.季节-纬度表面温度模型[0073]纬度(°)*一月三月五月七月九月十一月80亚极地亚极地亚极地中纬度中纬度亚极地70亚极地亚极地中纬度中纬度中纬度亚极地60中纬度中纬度中纬度亚极地亚极地中纬度50中纬度中纬度亚极地亚极地亚极地亚极地40亚极地亚极地亚极地中纬度中纬度亚极地30中纬度中纬度中纬度热带热带中纬度20热带热带热带热带热带热带10热带热带热带热带热带热带0热带热带热带热带热带热带-10热带热带热带热带热带热带-20热带热带热带中纬度中纬度热带-30中纬度中纬度中纬度中纬度中纬度中纬度-40亚极地亚极地亚极地亚极地亚极地亚极地-50亚极地亚极地亚极地中纬度中纬度亚极地-60中纬度中纬度中纬度中纬度中纬度中纬度-70中纬度中纬度中纬度中纬度中纬度中纬度-80中纬度中纬度中纬度中纬度中纬度中纬度[0074]*北纬为正值，南纬为负值[0075](2)大气水汽含量估算回归模型的选择与大气水汽含量估算[0076]根据步骤三的步骤(1)中获取的大气廓线类型信息和测量时记录的角度信息，可以从步骤一中的大气水汽含量估算模型查找表选择最优的大气水汽回归模型，以降低估算的不确定性。具体做法为：通过大气水汽含量估算模型查找表确定适合于该测量角度下的两个相邻的预设高斯角度对应的大气水汽含量估算模型，记为模型a和模型b。使用两个模型分别估算得到大气水汽含量wa和wb。为平衡两个模型误差，使用各自模型拟合时误差δ2的倒数作为权重，得到加权的大气水汽含量的估算值：[0077][0078]其中，表示模型a的拟合误差，表示模型b的拟合误差。[0079]得到各组大气水汽含量的估算值后，利用式(1)得到均值作为测量时空内大气水汽含量的最终估算值。[0080]基于同一发明构思，本发明的另一实施例提供一种基于地面红外辐射计的大气水汽含量估算系统，其包括：[0081]估算模型建立模块，用于基于大气辐射传输软件模拟不同大气条件和不同角度情况下大气下行热辐射比值与大气水汽含量的关系，建立利用大气下行热辐射比值估算大气水汽含量的回归模型；[0082]测量模块，用于利用热辐射计测量不同角度下的大气下行热辐射，并计算不同角度下的大气下行热辐射比值；[0083]大气水汽含量估算模块，用于将计算得到的不同角度下的大气下行热辐射比值输入所述大气下行热辐射比值与大气水汽含量的回归模型中，获得大气水汽含量估算值。[0084]以上实施方式仅用于说明本发明，而非对本发明的限制。尽管参照实施例对发明进行了详细说明，本领域的普通技术人员应当理解，对本发明的技术方案进行各种组合、修改或者等同替换，都不脱离本发明技术方案的精神和范围，均应涵盖在本发明的权利要求范围当中。当前第1页12当前第1页12