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

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

Info

Publication number
CN102519879B
CN102519879B CN201110400352.2A CN201110400352A CN102519879B CN 102519879 B CN102519879 B CN 102519879B CN 201110400352 A CN201110400352 A CN 201110400352A CN 102519879 B CN102519879 B CN 102519879B
Authority
CN
China
Prior art keywords
shortwave
calibration coefficient
passage
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.)
Expired - Fee Related
Application number
CN201110400352.2A
Other languages
English (en)
Other versions
CN102519879A (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)、前苏联的Meteor3(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 - R S 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.3498095 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 (5)

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

Families Citing this family (2)

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

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数据反演大气水汽含量方法

Also Published As

Publication number Publication date
CN102519879A (zh) 2012-06-27

Similar Documents

Publication Publication Date Title
Chander et al. Overview of intercalibration of satellite instruments
Washenfelder et al. Carbon dioxide column abundances at the Wisconsin Tall Tower site
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
CN110716185A (zh) 一种用于星载微波辐射计的在轨辐射校正方法
CN108827878B (zh) 一种地表气压的被动微波遥感探测方法
CN112798013B (zh) 一种对光学载荷在轨绝对辐射定标结果进行验证的方法
Stauffer et al. Propagation of radiosonde pressure sensor errors to ozonesonde measurements
CN113408111B (zh) 大气可降水量反演方法及系统、电子设备和存储介质
Wang et al. Water vapor retrievals from near-infrared channels of the advanced Medium Resolution Spectral Imager instrument onboard the Fengyun-3D satellite
Li et al. Comprehensive calibration and validation site for information remote sensing
CN113552080A (zh) 一种超光谱遥感地表大气相对湿度的实时反演算法
CN112596046A (zh) 一种用于星载大气探测激光雷达的高精度定标处理方法
Liu et al. Retrieval of atmospheric thermodynamic state from synergistic use of radio occultation and hyperspectral infrared radiances observations
CN102519879B (zh) 大气臭氧柱总量探测方法
Shimada et al. Calibration and validation of the ocean color version-3 product from ADEOS OCTS
Noël et al. Preliminary results of GOME-2 water vapour retrievals and first applications in polar regions
CN112162158B (zh) 在轨太赫兹探测仪天线指向失配的评估方法及装置
CN112345151B (zh) 一种基于自然大气的mwts-ii对海面气压的灵敏性测试方法
Jakobson et al. Diurnal variability of precipitable water in the Baltic region, impact on transmittance of the direct solar radiation
Thompson et al. Ozonesondes: Instrumentation and Data Applications
Mzé et al. Climatology and comparison of ozone from ENVISAT/GOMOS and SHADOZ/balloon-sonde observations in the southern tropics
Chanin The intercomparison ozone campaign held in France in June 1981: Description of the campaign
Boesch et al. Monitoring Greenhouses Gases over China Using Space-Based Observations
Ren et al. A Low-cost UAV Coordinated Carbon observation Network (LUCCN): an analysis of environment impact on ground base measurement node

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