CN102519879A - 大气臭氧柱总量探测方法 - Google Patents

大气臭氧柱总量探测方法 Download PDF

Info

Publication number
CN102519879A
CN102519879A CN2011104003522A CN201110400352A CN102519879A CN 102519879 A CN102519879 A CN 102519879A CN 2011104003522 A CN2011104003522 A CN 2011104003522A CN 201110400352 A CN201110400352 A CN 201110400352A CN 102519879 A CN102519879 A CN 102519879A
Authority
CN
China
Prior art keywords
passage
calibration coefficient
shortwave
ozone
low side
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2011104003522A
Other languages
English (en)
Other versions
CN102519879B (zh
Inventor
王维和
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
STATE SATELLITE METEROLOGICAL CENTER
National Satellite Meteorological Center
Original Assignee
STATE SATELLITE METEROLOGICAL CENTER
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by STATE SATELLITE METEROLOGICAL CENTER filed Critical STATE SATELLITE METEROLOGICAL CENTER
Priority to CN201110400352.2A priority Critical patent/CN102519879B/zh
Publication of CN102519879A publication Critical patent/CN102519879A/zh
Application granted granted Critical
Publication of CN102519879B publication Critical patent/CN102519879B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明涉及一种大气臭氧柱总量探测方法,包括:获取星载紫外臭氧总量探测仪的辐亮度探测数据；根据臭氧总量基准数据和辐射传输方程，逐级确定各个通道的在轨定标系数;从所述辐亮度探测数据中提取臭氧光谱信息，根据所述在轨定标系数对所述臭氧光谱信息进行反演,得到大气臭氧柱总量。本发明利用定标效果较好的卫星观测的臭氧总量数据，结合辐射传输方程，逐级对各个通道进行定标系数校正，解决了传统的交叉定标方法的局限，能够在星载紫外臭氧探测仪发射前实验室定标发生重大失误，或者卫星在运行过程中性能发生巨大改变导致星上定标设备失败的情况下，比较准确地测定紫外臭氧总量。

Description

大气臭氧柱总量探测方法
技术领域
本发明涉及卫星探测领域，尤其涉及一种大气臭氧柱总量探测方法。
背景技术
卫星探测器的定标是指确定卫星探测信号(模拟信号或数字信号)与物理量(例如温度、亮度等)之间的对应关系的过程。一般情况下，在卫星发射前对卫星探测器在实验室进行测试，确定用于将卫星探测信号转化成物理量的定标系数。遥感卫星探测器在轨道运行过程中，经常因卫星发射前实验室确定的定标系数不准确而导致被探测的物理量发生偏差。卫星运行过程中传感器的衰减也会导致被探测的物理量产生偏差。在轨定标系数订正是指对卫星在轨道运行过程中定标系数的订正。现有的在轨定标系数订正方法有三种：一是依赖卫星上定标设备的测量结果进行计算，二是通过地面同步观测配合卫星测量完成在轨定标系数的订正，三是通过不同卫星之间同步比对实现在轨定标系数的订正。一般情况下，通过上述三种方法可以追踪卫星在轨运行中性能的变化从而实现在轨定标系数的订正，但是，当卫星发射前在实验室定标过程中发生重大失误，或者卫星在运行过程中性能发生巨大改变时，利用上述方法进行在轨定标系数的订正会遇到困难。
臭氧是对生态、气候和环境有重要影响的大气微量气体。1985年Farman等首次发现了南极臭氧洞，之后在中国的青藏高原以及世界上其他高原地区又相继发现了很多微型的臭氧洞。研究臭氧洞的成因和臭氧变化的机理需要高时空分辨率的臭氧监测资料。早在1924年，英国物理学家和气象学家G.M.Dobson就发明了利用直射太阳光测量大气臭氧含量的紫外光谱仪，即Dobosn仪器，目前国际上很多地基臭氧观测站仍采用Dobson仪器观测大气臭氧总量，臭氧总量的单位也以他的名字命名，即Dobson(简写为DU)单位。和地基臭氧观测仪器相比，卫星臭氧探测仪(即搭载在卫星上的紫外臭氧探测仪)在时空方面具有明显的优势，使卫星臭氧探测成为监测全球臭氧变化的重要手段。利用大气对太阳紫外线的后向散射反演大气臭氧总量是由Dave等人在20世纪60年代提出的，之后美国、欧洲相继发射了十几个紫外臭氧探测仪，其中最早的是美国的TOMS系列仪器，分别搭载在美国的NIMBUS-7(1978-1993)和Earth Probe(简称EP，1996-2004)、前苏联的Meteor 3(1991-1994)卫星平台上。我国于2008年5月27日发射的“风云三号”极轨气象卫星也搭载了我国的第一台紫外臭氧总量探测仪(Total OzoneUnit，缩写为TOU)，该仪器与TOMS仪器具有类似的工作原理。
我国的星载紫外臭氧总量探测仪TOU的探测原理和定标方法与TOMS仪器是类似的，因此发射前的实验室定标以及发射后在轨定标方法与TOMS仪器也是相近的。就星载紫外臭氧总量探测仪而言，一般来说，卫星发射后，如果卫星状态没有大的变化，均采用发射前实验室测定的定标系数进行卫星探测信号到物理量的转换，当卫星状态发生变化时，会利用卫星上定标设备的观测值对卫星的性能或定标系数随时间的变化进行重新计算，或者也会采用一些替代方法例如不同卫星之间辐射量进行比较，从而进行相互定标。
通过卫星之间相互比对从而实现定标只能在卫星同步观测同一区域进行，即进行同步测量时才能发生作用，并且两个卫星观测点的几何参数要完全一致才能够互相比对。一般来说，采用星下点比对才能实现，即两颗卫星星下点同时观测同一地面目标。这种定标方法在NOAA卫星搭载的SBUV/2仪器的定标中获得了成功，被称为SNO(Simultaneous Nadir Overpass)方法。但是，当卫星发射前实验室定标发生重大错误，或者卫星进入轨道后探测器发生重大改变时如何对定标系数进行校正，这种方法并没有得到验证，该方法是直接比较两个卫星的探测量，因此最大的局限性就是两个卫星对同一地区的观测必须完全同步才能进行。这也是其他类型的卫星探测器常用的交叉定标技术。
我国风云三号卫星首次搭载的紫外臭氧总量探测仪与2008年5月发射升空，在分析数据时发现，反演的大气臭氧总量与国外同类卫星反演的结果以及地面观测结果相比有非常大的偏差，相对偏差甚至大于50％以上，这与仪器的设计指标(相对偏差为10％)有很大的差距。我们首先利用星上设备进行了分析，发现星上的定标设备即太阳辐照度观测结果也同样与正常的结果偏差很大，因此无法使用星上的定标设备进行误差订正。利用普遍采用的不同卫星之间交叉比对，发现由于条件的苛刻，即要求两个卫星观测必须完全同步的情况下才能实现，因此可以利用的比对样本数目不够，一般只能在高纬度地区找到少量满足要求的样本，但由于高纬度地区因太阳天顶角较大，大气散射的太阳紫外辐射比较弱，无法分析在整个动态范围内我国卫星观测结果与国外卫星观测结果之间存在的统计规律。初步的结论是，风云三号卫星紫外臭氧总量探测仪的实验室定标系数发生了重大失误，而星上定标设备也失败了，无法通过常规的方法解决定标系数的订正问题。如果无法解决，意味着卫星发射彻底失败。因此，必须找到仪器误差产生的原因，在此基础上找到订正方法。
由上可见，如何在星载紫外臭氧探测仪发射前实验室定标发生重大失误，或者卫星在运行过程中性能发生巨大改变导致星上定标设备失败的情况下，比较准确地测定紫外臭氧总量，是臭氧监测领域一个亟待解决的重要问题。
发明内容
本发明所要解决的技术问题是提供一种大气臭氧柱总量探测方法，能够在星载紫外臭氧探测仪发射前实验室定标发生重大失误，或者卫星在运行过程中性能发生巨大改变导致星上定标设备失败的情况下，比较准确地测定紫外臭氧总量，是臭氧监测领域一个亟待解决的重要问题。
为解决上述技术问题，本发明提出了一种大气臭氧柱总量探测方法，包括：
获取星载紫外臭氧总量探测仪的辐亮度探测数据；
根据臭氧总量基准数据和辐射传输方程，逐级确定各个通道的在轨定标系数；
从所述辐亮度探测数据中提取臭氧光谱信息，根据所述在轨定标系数对所述臭氧光谱信息进行反演，得到大气臭氧柱总量。
进一步地，上述方法还可具有以下特点，所述根据臭氧总量基准数据和辐射传输方程，逐级确定各个波段通道的在轨定标系数包括：
根据长波通道星下点在晴空海洋的辐亮度探测数据和辐射传输方程，确定长波通道能量低端星下点方向的在轨定标系数，所述长波通道的中心波长为360nm，宽度为1nm；
根据所述长波通道除星下点外的各个扫描方向在晴空海洋的辐亮度探测数据和辐射传输方程，确定所述长波通道除星下点外的各个扫描方向的能量低端的在轨定标系数；
根据长波通道各个扫描方向的能量低端的在轨定标系数和辐射传输方程，确定第一短波通道各个扫描方向能量低端的在轨定标系数，所述第一短波通道的中心波长为322nm，宽度为1nm；
根据所述长波通道能量低端的在轨定标系数和第一短波通道能量低端各个扫描方向的在轨定标系数、臭氧总量基准数据和辐射传输方程，确定第二短波通道各个扫描方向能量低端的在轨定标系数，所述第二短波通道的中心波长为312nm，宽度为1nm；
根据所述长波通道能量低端、第一短波通道能量低端以及第二短波通道整个动态范围各个扫描方向的在轨定标系数、臭氧总量基准数据和辐射传输方程，确定长波通道能量高端，第一短波通道能量高端以及其他通道各个扫描方向的在轨定标系数。
进一步地，上述方法还可具有以下特点，根据长波通道星下点在晴空海洋的辐亮度探测数据和辐射传输方程，确定长波通道能量低端星下点方向的在轨定标系数包括：
利用辐射传输方程计算星下点晴空海洋的紫外后向散射辐亮度Icorrected，
获取星载紫外臭氧总量探测仪星下点在晴空海洋的紫外后向散射辐亮度的探测数据Imeasured；
根据Icorrected和Imeasured，通过第一拟合公式确定拟合系数c0和c1，所述第一拟合公式如下：
Icorrected＝c0+c1×Imeasured
拟合系数c0和c1即为在轨定标系数；
根据每个扫描方向长波通道探测的辐亮度数据与辐射传输方程计算的晴空海洋的后向散射辐亮度数据，通过第一拟合公式得到长波通道各个扫描方向在轨定标系数。
进一步地，上述方法还可具有以下特点，根据长波通道各个扫描方向的能量低端的在轨定标系数和辐射传输方程，确定第一短波通道各个扫描方向能量低端的在轨定标系数包括：
根据长波通道各个扫描方向的能量低端的在轨定标系数，计算各个像元对应的第一下垫面反射率；
根据所述第一下垫面反射率和辐射传输方程得到所述第一短波通道的辐亮度计算数据；
对所述第一短波通道的辐亮度计算数据和所述第一短波通道的辐亮度探测数据用第二拟合公式进行拟合，所述获得第一短波通道能量低端各个扫描方向的在轨定标系数，其中，第二拟合公式如下：
Icorrected＝c0+c1×Imeasured+c2×Imeasured 2+c3×Imeasured 3。
进一步地，上述方法还可具有以下特点，根据所述长波通道能量低端的在轨定标系数和第一短波通道能量低端各个扫描方向的在轨定标系数、臭氧总量基准数据和辐射传输方程，确定第二短波通道各个扫描方向能量低端的在轨定标系数包括：
根据所述长波通道能量低端和第一短波通道能量低端各个扫描方向的在轨定标系数、臭氧总量基准数据，用辐射传输方程计算得到第二短波通道能量低端的辐亮度计算数据；
对所述第二短波通道能量低端的辐亮度计算数据和第二短波通道能量低端的辐亮度探测数据用第一拟合公式进行拟合，获得所述第二短波通道能量低端各个扫描方向的在轨定标系数；
将所述第二短波通道能量低端各个扫描方向的在轨定标系数外推到所述第二短波通道的整个动态范围，获得所述第二短波通道各个扫描方向的整个动态范围内的在轨定标系数。
进一步地，上述方法还可具有以下特点，根据所述长波通道能量低端、第一短波通道能量低端以及第二短波通道整个动态范围各个扫描方向的在轨定标系数、臭氧总量基准数据和辐射传输方程，确定长波通道能量高端，第一短波通道能量高端以及其他通道各个扫描方向的在轨定标系数：
根据所述长波通道能量低端和第一短波通道能量低端以及第二短波通道整个动态范围各个扫描方向的在轨定标系数，计算第二下垫面反射率；
根据所述第二下垫面反射率、臭氧总量基准数据，利用辐射传输方程，确定长波通道能量高端，第一短波通道能量高端以及除第二短波通道以外的其他通道的辐亮度计算数据；
对所述长波通道能量高端、第一短波通道能量高端以及其他通道辐亮度计算数据和辐亮度探测数据用第二拟合公式进行拟合，确定长波通道能量高端、第一短波通道能量高端以及除第二短波通道以外的其他通道各个扫描方向的在轨定标系数。
本发明利用其他卫星(定标效果较好的卫星)观测的臭氧总量数据，结合辐射传输方程，逐级对各个通道进行定标系数校正，解决了传统的交叉定标方法的局限，能够在星载紫外臭氧探测仪发射前实验室定标发生重大失误，或者卫星在运行过程中性能发生巨大改变导致星上定标设备失败的情况下，比较准确地测定紫外臭氧总量。
附图说明
具体实施方式
以下对本发明的原理和特征进行描述，所举实例只用于解释本发明，并非用于限定本发明的范围。
本发明的大气臭氧柱总量探测方法包括如下步骤：
步骤a，获取星载紫外臭氧总量探测仪的辐亮度探测数据；
星载紫外臭氧总量探测仪的辐亮度探测数据可以从卫星探测信号中获得。
步骤b，根据臭氧总量基准数据和辐射传输方程，逐级确定各个波段通道的在轨定标系数；
这里，臭氧总量基准数据可以是实验室定标准确的其他卫星的臭氧总量探测数据，也可以是地面观测得到的臭氧总量探测数据。
步骤c，从步骤a获取的辐亮度探测数据中提取臭氧光谱信息，根据步骤b确定的在轨定标系数对臭氧光谱信息进行反演，得到大气中臭氧柱总量，简称臭氧总量。
反演方法是本领域人员的公知常识，此处不再详述。
下面通过实施例，说明本发明的臭氧总量的测定方法中步骤b的实现过程。
根据辐亮度基准数据和辐射传输方程，逐级确定各个波段通道的在轨定标系数具体可以通过如下步骤实现：
步骤一，根据长波通道(中心波长为360nm，宽度为1nm)星下点在晴空海洋的辐亮度探测数据和辐射传输方程，确定长波通道能量低端的在轨定标系数，也即利用长波通道星下点在晴空海洋测量的结果，对长波通道动态范围的能量低端进行订正；
将卫星的观测几何参数作为输入，利用辐射传输方程计算长波通道的辐亮度(称为辐亮度计算数据)，再与星载紫外臭氧总量探测仪的辐亮度探测数据进行对比，即可确定长波通道能量低端的在轨定标系数(也可以称为辐射值订正系数)。其中，卫星的观测几何参数可以是卫星天顶角、卫星方位角、太阳天顶角、太阳方位角、经纬度等。
步骤一可以通过以下具体的子步骤实现：
步骤11，利用辐射传输方程计算星下点晴空海洋的紫外后向散射辐亮度Icorrected；
星下点是指卫星与地心的连线与地面交叉处对应的观测点。以我国风云三号卫星搭载的紫外臭氧总量探测仪TOU为例，太阳紫外线中的一部分辐射(波长区间为300nm-360nm)在大气的传输过程中经大气中臭氧分子的吸收、氧气和氮气的散射以及地球表面的反射后到达卫星高度，这部分紫外辐射也叫紫外后向散射，紫外臭氧总量探测仪TOU接收了紫外后向散射后，计算被臭氧吸收的紫外辐射所占的比例，从而可以获得大气中臭氧总含量。紫外臭氧总量探测仪TOU并不是连续测量300nm-360nm内全部的紫外后向散射，而是不连续地测量6个波长位置上宽度很窄的紫外波段内的紫外后向散射，这6个波段的宽度大约为1nm(纳米)左右，也可以称这6个波段为6个通道。由于太阳辐射的强度不一样，每个波段内太阳紫外辐射在大气中传播后在到卫星仪器高度上，其能量的分布也发生了变化，因此每个通道测量的紫外后向散射能量都具有不同的动态范围，一般来说波长越长，该通道的能量范围越大，将某一通道能量较低的部分称为能量低端，其余部分称为能量高端。对于紫外臭氧总量探测仪TOU，能量低端对应的模拟信号电压值范围大致为0～2V。由于长波通道臭氧吸收可以忽略，海洋对于紫外通道的反射率很低，在5％左右，因此可以首先确定长波通道能量低端的在轨定标系数。
卫星接收到的信号一般用计数值DC表示，如果实验室定标准确，可以通过实验室定标系数将计数值转化成物理量即辐亮度I，辐亮度I与DC通过实验室定标系数建立如公式(1)的线性关系：
I＝a×DC+b                公式(1)
公式(1)中，a、b为实验室定标系数。
紫外臭氧总量探测仪TOU在发射后，利用实验室定标系数得到的辐亮度进行臭氧总量的反演(即从卫星测量的辐亮度中提取臭氧光谱信息，从而计算出大气臭氧总量)，发现与其他卫星以及地面观测得到的大气臭氧总量几乎完全不一致，因此得到结论，紫外臭氧总量探测仪TOU发射前实验室定标系数存在问题，需要重新进行在轨定标。
紫外臭氧总量探测仪TOU在中心波长为360nm的通道6观测的紫外后向散射是不受大气中臭氧的影响的，当卫星在晴空海洋上进行观测时，卫星实际观测到的紫外后向散射I可以利用辐射传输方程计算。辐射传输方程如公式(2)所示：
I(λ，θ，θ0，Ω，P0，R)＝Ia(λ，θ，θ0，φ，Ω，P0)+Is(λ，θ，θ0，φ，Ω，P0，R)公式(2)
公式(2)中，各参数的含义如下：
I(λ，θ，θ0，Ω，P0，R)表示卫星观测到的后向散射辐亮度；
Ia(λ，θ，θ0，Φ，Ω，P0)表示卫星观测到的后向散射辐亮度中来自大气的贡献；
Is(λ，θ，θ0，Φ，Ω，P0，R)表示卫星观测到的后向散射辐亮度中来自地面的贡献；
λ表示波长，θ表示卫星天顶角，θ0表示太阳天顶角，Φ表示相对方位角，Ω表示臭氧总量，P0表示表面气压，R表示表面等效反射率。
表面反射贡献项可以表述公式(3)：
I s ( λ , θ , θ 0 , Ω , P 0 , R ) = RT ( λ , θ , θ 0 , Ω , P 0 ) 1 - RS b ( λ , Ω , P 0 ) 公式(3)
T(λ，θ，θ0，Ω，P0)＝Id(λ，θ，θ0，Ω，P0)f(λ，θ，Ω，P0)公式(4)
公式(3)和公式(4)中，Sb表示大气球反射率，Id表示到达P0处的直射漫射辐射的和，f表示表面方向(即卫星观测方向)反射比。
来自大气的辐射贡献可表述为公式(5)：
I a = I 0 + I 1 cos φ + I 2 cos 2 φ + R S T 1 - R S S b 公式(5)
公式(5)中，Φ为太阳矢量和卫星扫描平面之间的夹角，即相对方位角，RS为地表等效反射率。
在晴空海洋上空，对于长波通道(中心波长为360nm)，通常表面反射率很低，一般在5％左右，臭氧的吸收可以忽略，而观测几何参数如卫星天顶角、卫星方位角、太阳天顶角、太阳方位角、卫星和太阳之间的相对方位角都是通过卫星定位设备获取的，因此可以直接计算长波通道的后向散射辐亮度。
晴空的认定是这样的，在没有云的情况下，后向散射能量偏低。通常情况下，全球每天的晴空面积至少占总面积的10％以上，所以对于相同的观测几何条件，能量较低的少部分像元被视为晴空像元。
在计算出后向散射辐亮度I后，就可以通过拟合过程确定长波通道在晴空的情况下，也就是能量低端的在轨定标系数。
步骤12，获取星载紫外臭氧总量探测仪星下点在晴空海洋的紫外后向散射辐亮度的探测数据Imeasured；
步骤13，根据Icorrected和Imeasured，通过拟合公式确定拟合系数c0和c1，拟合公式如公式(6)。
Icorrected＝c0+c1×Imeasured    公式(6)
其中Icorrected代表辐射方程计算出的辐射值，Imeasured代表卫星观测值(即星载紫外臭氧总量探测仪的辐亮度探测数据)，c0和c1为拟合系数，也即在轨定标系数。
步骤14，根据每个扫描方向长波通道探测的辐亮度数据与辐射传输方程计算的晴空海洋的后向散射辐亮度数据，通过第一拟合公式得到长波通道各个扫描方向在轨定标系数。
步骤二，根据长波通道除星下点外的各个扫描方向在晴空海洋的辐亮度探测数据和辐射传输方程，确定长波通道除星下点外的各个扫描方向的能量低端的在轨定标系数；
紫外臭氧总量探测仪TOU在观测过程中不断自右向左进行扫描，每次扫描除了星下点外，左右各有15个观测方向。对每个观测方向，重复步骤1过程，可以获得长波通道在每个扫描方向的辐射订正系数，也就相当于得到了卫星观测随角度的变化规律，即卫星观测对扫描角的依赖关系。
步骤三，根据长波通道各个扫描方向的能量低端的在轨定标系数和辐射传输方程，确定第一短波通道(中心波长为322nm)能量低端各个扫描方向的在轨定标系数；
步骤三可以包括：
步骤31，根据长波通道各个扫描方向的能量低端的在轨定标系数，计算各个像元对应的第一下垫面反射率；
每个像元对应的下垫面反射率利用公式(7)计算：
R = ( I m - I a ) [ T - S b ( I m - I a ) ] 公式(7)
公式(7)中，Im为TOU的长波通道的观测辐亮度，Ia为公式(5)定义的来自大气的辐射贡献。
步骤32，根据第一下垫面反射率和辐射传输方程得到第一短波通道(中心波长为322nm)的辐亮度计算数据；
步骤33，对第一短波通道(中心波长为322nm)的辐亮度计算数据和第一短波通道(中心波长为322nm)的辐亮度观测数据进行拟合，获得第一短波通道(中心波长为322nm)能量低端各个扫描方向的在轨定标系数，其中，拟合公式为公式(8)。
Icorrected＝c0+c1×Imeasured+c2×Imeasured 2+c3×Imeasured 3公式(8)
公式(8)与公式(6)类似，增加了非线性订正系数c2和c3。
步骤四，根据长波通道(中心波长为360nm)能量低端的在轨定标系数和第一短波通道(中心波长为322nm)能量低端各个扫描方向的在轨定标系数、臭氧总量基准数据和辐射传输方程，确定第二短波通道(312nm)能量低端各个扫描方向的在轨定标系数；
步骤四可以包括：
步骤41，根据长波通道能量低端和第一短波通道能量低端各个扫描方向的在轨定标系数、臭氧总量基准数据，用辐射传输方程计算得到第二短波通道(中心波长为312nm)能量低端的辐亮度计算数据；
由于第二短波通道在整个动态范围内变化是线性的，因此第二短波通道能量低端在轨定标系数适用于第二短波通道整个动态范围。
步骤42，对第二短波通道能量低端的辐亮度计算数据和第二短波通道能量低端的辐亮度探测数据进行拟合，拟合公式为公式(6)，获得第二短波通道能量低端各个扫描方向的在轨定标系数；
步骤43，将第二短波通道能量低端各个扫描方向的在轨定标系数外推到第二短波通道的整个动态范围，获得第二短波通道的整个动态范围内各个扫描方向的在轨定标系数。
步骤五，根据长波通道能量低端、第一短波通道能量低端以及第二短波通道整个动态范围各个扫描方向的在轨定标系数、臭氧总量基准数据和辐射传输方程，确定长波通道能量高端，第一短波通道能量高端以及其他通道各个扫描方向的在轨定标系数。
步骤五可以包括：
步骤51，根据长波通道能量低端和第一短波通道能量低端以及第二短波通道整个动态范围各个扫描方向的在轨定标系数，计算第二下垫面反射率；
步骤52，根据第二下垫面反射率、臭氧总量基准数据，利用辐射传输方程，确定确定长波通道能量高端、第一短波通道能量高端以及除第二短波通道以外的其他通道的辐亮度计算数据；
步骤53，对长波通道能量高端、第一短波通道能量高端以及其他通道辐亮度计算数据和辐亮度探测数据进行拟合，拟合公式为公式(8)，确定长波通道能量高端、第一短波通道能量高端以及除第二短波通道以外的其他通道的在轨定标系数。
由上可见，本发明利用其它卫星观测的臭氧总量数据和辐射传输方程，从长波开始，逐级计算各个通道的在轨定标系数，最后给出所有通道的在轨定标系数。
表2～11是根据本发明方法，对表1中6个通道所获得的在轨定标系数结果。
表1通道序号及其对应的中心波长数据表
  通道   中心波长(nm)
  1   308.68
  2   312.59
  3   317.61
  4   322.40
  5   331.31
  6   360.11
表2通道1在轨定标系数数据表
  扫描序号   C0   C1   C2   C3
  1   0.0359033   0.6775104   0.2131372   -0.0940822
  2   0.0252710   0.7199647   0.1600021   -0.0596947
  3   0.0185551   0.7555540   0.1017969   -0.0301795
  4   0.0153165   0.7720801   0.0794389   -0.0207171
  5   0.0109418   0.7989720   0.0417998   -0.0074007
  6   0.0096461   0.8080102   0.0395559   -0.0079444
  7   0.0076898   0.8278320   0.0221054   -0.0021331
  8   0.0064840   0.8444843   0.0169983   -0.0010210
  9   0.0055192   0.8596525   0.0173554   -0.0019305
  10   0.0047153   0.8769161   0.0198828   -0.0031116
  11   0.0043342   0.8968122   0.0206546   -0.0034745
  12   0.0026803   0.9246721   0.0172966   -0.0022581
  13   0.0012419   0.9536421   0.0184009   -0.0030102
  14   0.0109568   0.9278868   0.0850720   -0.0237314
  15   -0.0010934   1.0213622   0.0199448   -0.0029662
  16   -0.0032501   1.0627820   0.0183565   -0.0026086
  17   -0.0090825   1.1298427   -0.0103031   0.0056987
  18   -0.0107194   1.1705892   -0.0031878   0.0007986
  19   -0.0144362   1.2224208   0.0002768   -0.0064032
  20   -0.0156030   1.2617406   0.0282390   -0.0264729
  21   -0.0195585   1.3270188   0.0063555   -0.0228630
  22   -0.0225820   1.3855111   -0.0078025   -0.0218295
  23   -0.0297946   1.4670979   -0.0682518   -0.0012608
  24   -0.0304361   1.4967075   -0.0452498   -0.0183590
  25   -0.0409466   1.6018071   -0.1680849   0.0401728
  26   -0.0443106   1.6321200   -0.1500327   0.0247877
  27   -0.0515302   1.6951851   -0.2359937   0.0797442
  28   -0.0522168   1.6976664   -0.2362135   0.0950354
  29   -0.0374187   1.5945346   -0.1052385   0.0569727
  30   -0.0498965   1.6555638   -0.2422923   0.1492256
  31   -0.0453921   1.6028749   -0.2590242   0.2029319
表3通道2在轨定标系数数据表
  扫描序号   C0   C1
  1   0.8154903   0.0206316
  2   0.8275054   0.0104101
  3   0.8283755   0.0070976
  4   0.8250939   0.0078477
  5   0.8196939   0.0100395
  6   0.8193689   0.0123852
  7   0.8227852   0.0118101
  8   0.8297912   0.0112010
  9   0.8392259   0.0109319
  10   0.8537546   0.0082781
  11   0.8702817   0.0054194
  12   0.8911720   0.0014864
  13   0.9153658   -0.0028190
  14   0.9443767   -0.0082448
  15   0.9771625   -0.0138943
  16   1.0120161   -0.0207612
  17   1.0445996   -0.0204913
  18   1.0812056   -0.0246719
  19   1.1212193   -0.0290828
  20   1.1619139   -0.0333946
  21   1.2027173   -0.0381964
  22   1.2445063   -0.0450735
  23   1.2748307   -0.0420634
  24   1.3068073   -0.0434295
  25   1.3320191   -0.0415200
  26   1.3505832   -0.0364045
  27   1.3717854   -0.0394490
  28   1.3770808   -0.0398079
  29   1.3730069   -0.0437393
  30   1.3590504   -0.0477537
  31   1.3064378   -0.0353390
表4通道3能量低端在轨定标系数数据表
  扫描序号   C0   C1   C2   C3
  1   0.1015035   0.8408877   0.0220264   -0.0025657
  2   0.0725485   0.8616217   0.0170489   -0.0020572
  3   0.0588935   0.8747478   0.0101021   -0.0012879
  4   0.0642179   0.8666194   0.0091118   -0.0011576
  5   0.0683435   0.8644880   0.0054063   -0.0007251
  6   0.0675216   0.8650547   0.0039703   -0.0005943
  7   0.0716849   0.8672447   0.0021824   -0.0003644
  8   0.0692047   0.8750809   0.0008946   -0.0002409
  9   0.0703852   0.8760163   0.0035158   -0.0004954
  10   0.0692308   0.8819324   0.0057538   -0.0006902
  11   0.0651522   0.8953096   0.0064509   -0.0007570
  12   0.0574840   0.9133390   0.0074886   -0.0008770
  13   0.0481547   0.9353966   0.0083270   -0.0009590
  14   0.0367829   0.9623843   0.0087014   -0.0009639
  15   0.0285531   0.9905674   0.0099574   -0.0010538
  16   0.0140792   1.0234101   0.0114971   -0.0012519
  17   0.0159967   1.0372425   0.0205763   -0.0022176
  18   0.0048857   1.0698562   0.0234575   -0.0025654
  19   -0.0012569   1.1013674   0.0267825   -0.0028784
  20   -0.0037419   1.1290983   0.0323871   -0.0035031
  21   -0.0148018   1.1715717   0.0321283   -0.0034497
  22   -0.0160343   1.2005875   0.0367749   -0.0039768
  23   -0.0108438   1.2179711   0.0452950   -0.0050603
  24   -0.0034403   1.2344565   0.0521665   -0.0059073
  25   0.0019383   1.2491233   0.0578312   -0.0065940
  26   0.0016243   1.2659084   0.0607158   -0.0068802
  27   -0.0034022   1.2879734   0.0576674   -0.0064389
  28   -0.0400285   1.3469356   0.0367482   -0.0043283
  29   -0.0404876   1.3633407   0.0209752   -0.0025153
  30   -0.0686327   1.3700891   0.0133294   -0.0019163
  31   -0.0146560   1.2812453   0.0238393   -0.0032292
表5通道3能量高端在轨定标系数数据表
  扫描序号   C0   C1   C2   C3
  1   0.8867492   0.7146997   -0.0049710   0.0017622
  2   -4.9790075   3.3617248   -0.3973862   0.0209284
  3   8.3966531   -3.0425430   0.6176057   -0.0323852
  4   10.6856790   -4.0971691   0.7761099   -0.0403408
  5   9.5456171   -3.5706789   0.6922652   -0.0359259
  6   7.0439930   -2.4492899   0.5250098   -0.0277050
  7   5.8598219   -1.9158016   0.4453821   -0.0237473
  8   4.3943071   -1.2579944   0.3490634   -0.0190541
  9   3.6237487   -0.9221577   0.3031999   -0.0169524
  10   2.5764858   -0.4534018   0.2379000   -0.0138977
  11   2.1183862   -0.2220007   0.2039229   -0.0121381
  12   2.4489002   -0.3375462   0.2215911   -0.0128874
  13   2.0464265   -0.1597049   0.2021339   -0.0120984
  14   -1.0255090   1.1466362   0.0279209   -0.0043427
  15   1.7475645   -0.0143832   0.1959299   -0.0120967
  16   6.1908674   -1.8739624   0.4596028   -0.0240940
  17   8.0735795   -2.6683882   0.5805688   -0.0299964
  18   10.7234520   -3.9248534   0.7856392   -0.0405999
  19   14.6755920   -5.7744551   1.0815657   -0.0557870
  20   26.4107740   -11.2555710   1.9438540   -0.1005173
  21   30.4419710   -13.1596050   2.2520771   -0.1165352
  22   -52.6848350   27.5399250   -4.3557640   0.2402669
  23   40.6225240   -20.1283610   3.7300285   -0.2143268
  24   37.4359470   -17.5900780   3.1694638   -0.1759914
  25   5.2751910   -2.1567507   0.7221737   -0.0470891
  26   50.0400220   -22.4542290   3.7708403   -0.1982070
  27   254.8921800   -122.4624900   20.0006330   -1.0731531
  28   25.2156390   -10.6327130   1.9056208   -0.1001459
  29   17.7863350   -7.1766804   1.3697614   -0.0727249
  30   35.2139390   -15.6468630   2.7310534   -0.1456289
  31   28.3103490   -12.0235520   2.0895919   -0.1091992
表6通道4能量低端在轨定标系数数据表
  扫描序号   C0   C1
  1   0.9223665   0.0991571
  2   0.9262605   0.0802458
  3   0.9208407   0.0760403
  4   0.9071272   0.0887456
  5   0.8937509   0.1031730
  6   0.8867096   0.1117618
  7   0.8845692   0.1143139
  8   0.8868585   0.1139396
  9   0.8947053   0.1066570
  10   0.9067596   0.0974699
  11   0.9215302   0.0874597
  12   0.9415537   0.0729339
  13   0.9653577   0.0574393
  14   0.9938876   0.0396949
  15   1.0269352   0.0204767
  16   1.0625540   0.0005445
  17   1.1010019   -0.0200707
  18   1.1405954   -0.0386427
  19   1.1810891   -0.0541650
  20   1.2235477   -0.0704812
  21   1.2655789   -0.0841175
  22   1.3063062   -0.0969309
  23   1.3428900   -0.1034473
  24   1.3762936   -0.1088127
  25   1.4054329   -0.1147424
  26   1.4278483   -0.1171278
  27   1.4440783   -0.1226279
  28   1.4355886   -0.0931265
  29   1.3999033   -0.0255106
  30   1.3681026   0.0031532
  31   1.2909452   0.0851826
表7通道4能量高端在轨定标系数数据表
  扫描序号   C0   C1   C2   C3
  1   -17.6240580   7.4622103   -0.7738206   0.0285888
  2   -20.0920900   8.3928856   -0.8887887   0.0332496
  3   -15.9700610   6.7970058   -0.6873394   0.0248928
  4   -10.3645860   4.7705972   -0.4483975   0.0155572
  5   -3.9563210   2.3808527   -0.1585344   0.0040056
  6   -8.9712340   4.2146315   -0.3808702   0.0128561
  7   -10.1641310   4.6327230   -0.4298309   0.0147698
  8   -2.6043536   1.9591743   -0.1190211   0.0029266
  9   -7.9364512   3.8597397   -0.3399282   0.0113721
  10   -9.8015822   4.5313786   -0.4164526   0.0142548
  11   -8.5027096   4.1102470   -0.3688735   0.0125405
  12   -7.7539079   3.8691214   -0.3396568   0.0114479
  13   -8.9579450   4.3259729   -0.3909522   0.0133995
  14   -9.7661686   4.6118773   -0.4175851   0.0142572
  15   -0.3860137   1.3793467   -0.0454584   0.0003324
  16   -10.9344820   5.1152142   -0.4716709   0.0163743
  17   -13.7959630   6.2187911   -0.6013621   0.0214075
  18   -17.3532960   7.4579198   -0.7358172   0.0263531
  19   -17.5600080   7.4994325   -0.7263527   0.0255459
  20   -29.7222450   12.0471690   -1.2741921   0.0472665
  21   -63.1757720   25.0090080   -2.9253556   0.1170128
  22   -112.4609800   45.7123710   -5.7898556   0.2482439
  23   15.9663100   -6.1456718   1.1653109   -0.0611309
  24   40.4065290   -16.9977210   2.7705718   -0.1396986
  25   103.8926100   -41.1738580   5.8136675   -0.2657640
  26   23.2706280   -8.7431933   1.4799494   -0.0728848
  27   -71.6183050   28.7397570   -3.4244760   0.1399893
  28   -17.1337870   7.5258402   -0.6878297   0.0230908
  29   -15.1381360   6.7980494   -0.6027742   0.0197080
  30   -14.4475130   6.4074191   -0.5454898   0.0169023
  31   -18.4796680   8.2515094   -0.8318841   0.0305892
表8通道5能量低端在轨定标系数数据表
  扫描序号   C0   C1   C2   C3
  1   2.3691183   -0.5120454   0.2748754   -0.0183541
  2   1.6466773   -0.1325703   0.2056645   -0.0141775
  3   0.4802968   0.5857385   0.0576614   -0.0043025
  4   0.6312558   0.4657858   0.0859563   -0.0065609
  5   0.0520221   0.8348535   0.0094193   -0.0015897
  6   -0.0266176   0.8868530   -0.0017071   -0.0009178
  7   0.0647731   0.8099327   0.0192695   -0.0027151
  8   -0.0165992   0.8656729   0.0081541   -0.0019882
  9   0.1350258   0.7577083   0.0341257   -0.0038870
  10   0.1167329   0.7766850   0.0305838   -0.0035444
  11   0.2288328   0.7004961   0.0500393   -0.0048924
  12   0.2615744   0.6847387   0.0554209   -0.0051441
  13   0.2412759   0.7152930   0.0499964   -0.0045440
  14   0.3570654   0.6352138   0.0733063   -0.0061827
  15   0.2688481   0.7160593   0.0587308   -0.0049701
  16   0.3092071   0.7050895   0.0652116   -0.0052333
  17   0.3358937   0.7008710   0.0720565   -0.0056438
  18   0.3369952   0.7266532   0.0707715   -0.0053512
  19   0.3258842   0.7628031   0.0678712   -0.0049413
  20   0.4161581   0.7182175   0.0860899   -0.0062351
  21   0.4925421   0.6872258   0.1013331   -0.0073224
  22   0.5538694   0.6651980   0.1148779   -0.0083022
  23   0.6412183   0.6340187   0.1284614   -0.0092217
  24   0.6218802   0.6816637   0.1230970   -0.0087794
  25   0.5709840   0.7517297   0.1118674   -0.0079645
  26   0.5758327   0.7798066   0.1095100   -0.0078507
  27   0.5072878   0.8409070   0.1005273   -0.0073662
  28   0.3350175   0.9782107   0.0732581   -0.0058052
  29   0.4480879   0.8759415   0.1035262   -0.0086963
  30   0.0682541   1.1511858   0.0371428   -0.0041423
  31   -0.1112393   1.3066127   -0.0025346   -0.0021281
表9通道5能量高端在轨定标系数数据表
  扫描序号   C0   C1   C2   C3
  1   2.2370897   0.5139213   0.0011845   -0.0000917
  2   1.9819164   0.5684713   -0.0019705   -0.0000386
  3   1.6392560   0.6478466   -0.0080359   0.0000973
  4   1.6013284   0.6498084   -0.0088955   0.0001177
  5   1.6102111   0.6492417   -0.0101173   0.0001600
  6   1.6379207   0.6419073   -0.0102034   0.0001654
  7   1.5943224   0.6544018   -0.0115747   0.0002030
  8   1.6906223   0.6386363   -0.0107808   0.0001914
  9   1.7560909   0.6304957   -0.0100794   0.0001765
  10   1.8330348   0.6228944   -0.0092440   0.0001591
  11   1.8581815   0.6297368   -0.0093587   0.0001645
  12   1.8981929   0.6379446   -0.0093500   0.0001644
  13   1.9841204   0.6400456   -0.0088357   0.0001563
  14   2.0981545   0.6368481   -0.0076398   0.0001333
  15   2.1779081   0.6485389   -0.0074797   0.0001341
  16   2.2699492   0.6586377   -0.0069854   0.0001221
  17   2.4068329   0.6529016   -0.0044931   0.0000466
  18   2.5476244   0.6530052   -0.0028647   0.0000043
  19   2.8247412   0.6245551   0.0009916   -0.0000899
  20   3.4391298   0.5035327   0.0128660   -0.0004023
  21   3.5843084   0.4990426   0.0158731   -0.0004999
  22   4.1951140   0.3736517   0.0286718   -0.0008518
  23   4.3246168   0.3729798   0.0314988   -0.0009732
  24   4.6167742   0.3251103   0.0377105   -0.0011609
  25   4.7906357   0.3076873   0.0409099   -0.0012638
  26   4.2542973   0.4910837   0.0250853   -0.0008060
  27   4.1585975   0.5433144   0.0204231   -0.0006549
  28   3.8085929   0.6527147   0.0092898   -0.0003200
  29   4.1970951   0.5656113   0.0121726   -0.0003239
  30   4.4161326   0.4992916   0.0146866   -0.0003631
  31   4.8627719   0.3500106   0.0217718   -0.0005204
表10通道6能量低端在轨定标系数数据表
  扫描序号   C0   C1
  1   0   0.8684080
  2   0   0.8630570
  3   0   0.8589060
  4   0   0.8561410
  5   0   0.8549200
  6   0   0.8553790
  7   0   0.8576390
  8   0   0.8618070
  9   0   0.8679860
  10   0   0.8762710
  11   0   0.8867510
  12   0   0.8995140
  13   0   0.9146350
  14   0   0.9321820
  15   0   0.9522000
  16   0   0.9747040
  17   0   0.9996680
  18   0   1.0269980
  19   0   1.0565110
  20   0   1.0879030
  21   0   1.1207090
  22   0   1.1542650
  23   0   1.1876600
  24   0   1.2197030
  25   0   1.2489070
  26   0   1.2735090
  27   0   1.2915360
  28   0   1.3009450
  29   0   1.2998290
  30   0   1.2866700
  31   0   1.2606020
表11通道6能量高端在轨定标系数数据表
  扫描序号   C0   C1   C2   C3
  1   0.9988047   0.7579074   -0.0096199   0.0000885
  2   -0.0225412   0.9727928   -0.0212594   0.0002741
  3   -0.3953856   1.0794272   -0.0289221   0.0004234
  4   -0.0368804   1.0303416   -0.0278741   0.0004182
  5   0.3036928   0.9791600   -0.0268795   0.0004233
  6   0.7444024   0.9020927   -0.0236350   0.0003766
  7   1.0411747   0.8541061   -0.0219664   0.0003615
  8   1.4259679   0.7820873   -0.0182421   0.0003014
  9   1.6992849   0.7310034   -0.0153458   0.0002532
  10   1.9314219   0.6891667   -0.0126666   0.0002060
  11   2.0883346   0.6637498   -0.0108864   0.0001772
  12   2.2480016   0.6403806   -0.0089025   0.0001406
  13   2.3653011   0.6283937   -0.0074627   0.0001150
  14   2.5769215   0.5996943   -0.0048458   0.0000689
  15   2.6096403   0.6129241   -0.0046781   0.0000675
  16   2.6491435   0.6278139   -0.0041958   0.0000531
  17   2.6693343   0.6436419   -0.0030550   0.0000120
  18   2.6608823   0.6731230   -0.0030309   0.0000026
  19   2.7455452   0.6866575   -0.0020138   -0.0000257
  20   3.1193737   0.6311935   0.0040241   -0.0001714
  21   3.2271128   0.6404805   0.0055343   -0.0002187
  22   3.6706764   0.5714301   0.0125472   -0.0003867
  23   3.7329516   0.5976308   0.0125345   -0.0004109
  24   4.0315614   0.5632692   0.0164054   -0.0005050
  25   4.1178604   0.5738854   0.0170849   -0.0005352
  26   3.9761451   0.6531118   0.0106443   -0.0003600
  27   4.5227685   0.5648966   0.0157713   -0.0004384
  28   4.8349322   0.5246662   0.0152332   -0.0003651
  29   5.7295358   0.2498095   0.0224354   -0.0004479
  30   6.3103705   0.2356152   0.0264354   -0.0004944
  31   7.5649342   -0.0253438   0.0368444   -0.0006560
通过本发明方法获得的定标系数后再进行大气臭氧总量反演，与地面观测结果相比相对偏差平方根为3％左右，而相对偏差平均值小于1％，远远优于设计指标(10％)，达到了国际先进水平。并且，通过使用本发明方法对风云三号星载紫外臭氧总量探测仪进行定标系数订正实验，证实了本发明方法是有效的。
由上可见，本发明利用其他卫星(定标效果较好的卫星)观测的臭氧总量数据，结合辐射传输方程，逐级对各个通道进行定标系数校正，解决了传统的交叉定标方法的局限(传统交叉定标要求两个卫星必须同步观测，并且其观测几何参数严格一致，才能直接对观测量进行对比，这样就无法获取足够的样本，在整个动态范围内进行定标系数的校正)，能够在星载紫外臭氧探测仪发射前实验室定标发生重大失误，或者卫星在运行过程中性能发生巨大改变导致星上定标设备失败的情况下，比较准确地测定大气臭氧柱总量探。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (6)

1.一种大气臭氧柱总量探测方法,其特征在于,包括：
获取星载紫外臭氧总量探测仪的辐亮度探测数据；
根据臭氧总量基准数据和辐射传输方程，逐级确定各个通道的在轨定标系数;
从所述辐亮度探测数据中提取臭氧光谱信息，根据所述在轨定标系数对所述臭氧光谱信息进行反演,得到大气臭氧柱总量。
2.根据权利要求1所述的大气臭氧柱总量探测方法, 其特征在于,所述根据臭氧总量基准数据和辐射传输方程，逐级确定各个波段通道的在轨定标系数包括：
根据长波通道星下点在晴空海洋的辐亮度探测数据和辐射传输方程，确定长波通道能量低端星下点方向的在轨定标系数，所述长波通道的中心波长为360nm，宽度为1 nm；
根据所述长波通道除星下点外的各个扫描方向在晴空海洋的辐亮度探测数据和辐射传输方程，确定所述长波通道除星下点外的各个扫描方向的能量低端的在轨定标系数；
根据长波通道各个扫描方向的能量低端的在轨定标系数和辐射传输方程，确定第一短波通道各个扫描方向能量低端的在轨定标系数，所述第一短波通道的中心波长为322nm，宽度为1 nm；
根据所述长波通道能量低端的在轨定标系数和第一短波通道能量低端各个扫描方向的在轨定标系数、臭氧总量基准数据和辐射传输方程，确定第二短波通道各个扫描方向能量低端的在轨定标系数，所述第二短波通道的中心波长为312nm，宽度为1 nm；
根据所述长波通道能量低端、第一短波通道能量低端以及第二短波通道整个动态范围各个扫描方向的在轨定标系数、臭氧总量基准数据和辐射传输方程，确定长波通道能量高端，第一短波通道能量高端以及其他通道各个扫描方向的在轨定标系数。
3.根据权利要求2所述的大气臭氧柱总量探测方法, 其特征在于, 根据长波通道星下点在晴空海洋的辐亮度探测数据和辐射传输方程，确定长波通道能量低端星下点方向的在轨定标系数包括：
利用辐射传输方程计算星下点晴空海洋的紫外后向散射辐亮度Icorrected；
获取星载紫外臭氧总量探测仪星下点在晴空海洋的紫外后向散射辐亮度的探测数据Imeasured；
根据Icorrected和Imeasured，通过第一拟合公式确定拟合系数c0和c1，所述第一拟合公式如下：
Figure 2011104003522100001DEST_PATH_IMAGE001
拟合系数c0和c1即为在轨定标系数；
根据每个扫描方向长波通道探测的辐亮度数据与辐射传输方程计算的晴空海洋的后向散射辐亮度数据，通过第一拟合公式得到长波通道各个扫描方向在轨定标系数。
4.根据权利要求2所述的大气臭氧柱总量探测方法, 其特征在于, 根据长波通道各个扫描方向的能量低端的在轨定标系数和辐射传输方程，确定第一短波通道各个扫描方向能量低端的在轨定标系数包括：
根据长波通道各个扫描方向的能量低端的在轨定标系数，计算各个像元对应的第一下垫面反射率；
根据所述第一下垫面反射率和辐射传输方程得到所述第一短波通道的辐亮度计算数据；
对所述第一短波通道的辐亮度计算数据和所述第一短波通道的辐亮度探测数据用第二拟合公式进行拟合，所述获得第一短波通道能量低端各个扫描方向的在轨定标系数，其中，第二拟合公式如下：
。
5.根据权利要求2所述的大气臭氧柱总量探测方法, 其特征在于,根据所述长波通道能量低端的在轨定标系数和第一短波通道能量低端各个扫描方向的在轨定标系数、臭氧总量基准数据和辐射传输方程，确定第二短波通道各个扫描方向能量低端的在轨定标系数包括：
    根据所述长波通道能量低端和第一短波通道能量低端各个扫描方向的在轨定标系数、臭氧总量基准数据，用辐射传输方程计算得到第二短波通道能量低端的辐亮度计算数据；
对所述第二短波通道能量低端的辐亮度计算数据和第二短波通道能量低端的辐亮度探测数据用第一拟合公式进行拟合，获得所述第二短波通道能量低端各个扫描方向的在轨定标系数；
    将所述第二短波通道能量低端各个扫描方向的在轨定标系数外推到所述第二短波通道的整个动态范围，获得所述第二短波通道各个扫描方向的整个动态范围内的在轨定标系数。
6.根据权利要求2所述的大气臭氧柱总量探测方法, 其特征在于,根据所述长波通道能量低端、第一短波通道能量低端以及第二短波通道整个动态范围各个扫描方向的在轨定标系数、臭氧总量基准数据和辐射传输方程，确定长波通道能量高端，第一短波通道能量高端以及其他通道各个扫描方向的在轨定标系数：
    根据所述长波通道能量低端和第一短波通道能量低端以及第二短波通道整个动态范围各个扫描方向的在轨定标系数，计算第二下垫面反射率；
根据所述第二下垫面反射率、臭氧总量基准数据，利用辐射传输方程，确定长波通道能量高端，第一短波通道能量高端以及除第二短波通道以外的其他通道的辐亮度计算数据；
    对所述长波通道能量高端、第一短波通道能量高端以及其他通道辐亮度计算数据和辐亮度探测数据用第二拟合公式进行拟合，确定长波通道能量高端、第一短波通道能量高端以及除第二短波通道以外的其他通道各个扫描方向的在轨定标系数。
CN201110400352.2A 2011-12-05 2011-12-05 大气臭氧柱总量探测方法 Expired - Fee Related CN102519879B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201110400352.2A CN102519879B (zh) 2011-12-05 2011-12-05 大气臭氧柱总量探测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110400352.2A CN102519879B (zh) 2011-12-05 2011-12-05 大气臭氧柱总量探测方法

Publications (2)

Publication Number Publication Date
CN102519879A true CN102519879A (zh) 2012-06-27
CN102519879B CN102519879B (zh) 2014-12-17

Family

ID=46290875

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110400352.2A Expired - Fee Related CN102519879B (zh) 2011-12-05 2011-12-05 大气臭氧柱总量探测方法

Country Status (1)

Country Link
CN (1) CN102519879B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105737980A (zh) * 2016-02-14 2016-07-06 国家卫星气象中心 辐射定标方法及装置
CN115964677A (zh) * 2023-01-03 2023-04-14 国家卫星气象中心(国家空间天气监测预警中心) 一种遥感臭氧数据集的生成方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101504353A (zh) * 2008-02-06 2009-08-12 香港科技大学 用于提供近地面悬浮颗粒物分布的方法和相应的***
CN101598543A (zh) * 2009-07-29 2009-12-09 中国科学院对地观测与数字地球科学中心 一种实用的遥感影像大气校正方法
CN101936877A (zh) * 2010-03-05 2011-01-05 中国农业科学院农业资源与农业区划研究所 从modis数据反演大气水汽含量方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101504353A (zh) * 2008-02-06 2009-08-12 香港科技大学 用于提供近地面悬浮颗粒物分布的方法和相应的***
CN101598543A (zh) * 2009-07-29 2009-12-09 中国科学院对地观测与数字地球科学中心 一种实用的遥感影像大气校正方法
CN101936877A (zh) * 2010-03-05 2011-01-05 中国农业科学院农业资源与农业区划研究所 从modis数据反演大气水汽含量方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105737980A (zh) * 2016-02-14 2016-07-06 国家卫星气象中心 辐射定标方法及装置
CN115964677A (zh) * 2023-01-03 2023-04-14 国家卫星气象中心(国家空间天气监测预警中心) 一种遥感臭氧数据集的生成方法

Also Published As

Publication number Publication date
CN102519879B (zh) 2014-12-17

Similar Documents

Publication Publication Date Title
Washenfelder et al. Carbon dioxide column abundances at the Wisconsin Tall Tower site
Chander et al. Overview of intercalibration of satellite instruments
Nassar et al. Validation of Tropospheric Emission Spectrometer (TES) nadir ozone profiles using ozonesonde measurements
Buchwitz et al. Carbon monoxide, methane and carbon dioxide columns retrieved from SCIAMACHY by WFM-DOAS: year 2003 initial data set
CN106569186B (zh) 一种星载微波辐射计的偏差校正方法
Bock et al. Accuracy assessment of water vapour measurements from in-situ and remote sensing techniques during the DEMEVAP 2011 campaign at OHP
Hong et al. Retrieval of global carbon dioxide from TanSat satellite and comprehensive validation with TCCON measurements and satellite observations
CN110716185A (zh) 一种用于星载微波辐射计的在轨辐射校正方法
CN108827878B (zh) 一种地表气压的被动微波遥感探测方法
CN112798013B (zh) 一种对光学载荷在轨绝对辐射定标结果进行验证的方法
Li et al. Comprehensive calibration and validation site for information remote sensing
Wang et al. Water vapor retrievals from near-infrared channels of the advanced Medium Resolution Spectral Imager instrument onboard the Fengyun-3D satellite
Wang et al. Analysis for retrieval and validation results of FY-3 Total Ozone Unit (TOU)
Antón et al. Validation of GOME-2/MetOp-A total water vapour column using reference radiosonde data from the GRUAN network
Liu et al. A review of low Earth orbit occultation using microwave and infrared-laser signals for monitoring the atmosphere and climate
Noël et al. Preliminary results of GOME-2 water vapour retrievals and first applications in polar regions
CN102519879B (zh) 大气臭氧柱总量探测方法
Martins et al. Ozone correlations between mid-tropospheric partial columns and the near-surface at two mid-atlantic sites during the DISCOVER-AQ campaign in July 2011
Jakobson et al. Diurnal variability of precipitable water in the Baltic region, impact on transmittance of the direct solar radiation
Zhou et al. Ground-Based Atmospheric CO2, CH4, and CO Column Measurements at Golmud in the Qinghai-Tibetan Plateau and Comparisons with TROPOMI/S5P Satellite Observations
Thompson et al. Ozonesondes: Instrumentation and Data Applications
Zhang et al. Comparison analysis of global methane concentration derived from SCIAMACHY, AIRS, and GOSAT with surface station measurements
Rault et al. Validation of ozone profiles retrieved from Stratospheric Aerosol and Gas Experiment III limb scatter measurements
Mzé et al. Climatology and comparison of ozone from ENVISAT/GOMOS and SHADOZ/balloon-sonde observations in the southern tropics
Hendrick et al. Multi-year comparison of stratospheric BrO vertical profiles retrieved from SCIAMACHY limb and ground-based UV-visible measurements

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20141217

Termination date: 20151205

EXPY Termination of patent right or utility model