HJ-1BB08的有效波段宽度计算方法及定标方法专利检索-低温物理学物理专利检索查询-专利查询网

HJ-1BB08的有效波段宽度计算方法及定标方法

阅读：272发布：2020-05-26

专利汇可以提供HJ-1BB08的有效波段宽度计算方法及定标方法专利检索，专利查询，专利分析的服务。并且本发明公开了一种HJ-1B B08的有效波段宽度计算方法及定标方法，具体为：1)卫星在发射前根据实验室里面的标准面源黑体数据建立面源黑体温度与辐照度转换方程，计算在不同面源黑体温度时的辐照度；2)根据地面标准面源黑体下的比辐射率、温度和通道响应函数，构建有效辐亮度与面源黑体温度的方程，并计算在特定面源黑体温度时面源黑体的有效辐亮度；3)根据光学物理学中，辐照度和辐亮度存在的基本关系，构建基于地面实验室数据的有效波段宽度和面源黑体温度的查找表，得到不同黑体温度条件下的有效波段宽度。本发明的有效波段宽度计算方法，不仅考虑了通道响应函数，还充分考虑了其它可能影响到有效波段宽度值的参数，其计算出的有效波段宽度值准确能很好的应用于星上定标。，下面是HJ-1BB08的有效波段宽度计算方法及定标方法专利的具体信息内容。

权利要求

1.一种HJ-1B B08的有效波段宽度计算的方法，具体为：1)HJ-1B B08通道在发射前根据实验室里面的标准面源黑体数据建立面源黑体温度与辐照度转换方程，计算在不同面源黑体温度Tm时的辐照度NTm；2)根据地面标准面源黑体下的比辐射率εm、温度Tm和通道响应函数CRF，利用Planck方程，构建有效辐亮度Lm与面源黑体温度Tm的方程，并计算在特定面源黑体温度Tm时面源黑体的有效辐亮度Lm；3)根据光学物理学中，辐照度和辐亮度存在的基本关系，构建基于地面实验室数据的有效波段宽度和面源黑体温度的查找表，得到不同黑体温度Tm条件下的有效波段宽度Δλ。
2.如权利要求1所述的HJ-1B B08的有效波段宽度的方法，其特征在于，所述HJ-1B卫星的面源黑体温度与辐照度转换方程为 k0，k1，k2，k3是已知
常数。
3.如权利要求1所述的HJ-1B B08的有效波段宽度的方法，其特征在于，
所述有效辐亮度Lm与面源黑体温度Tm的方程为
式中h为Planck常数，k为玻耳兹曼常数，c为光速，f(λ)为传感器通道响应函数CRF，其中方程。
4.如权利要求1所述的HJ-1B B08的有效波段宽度的方法，其特征在于，辐照度和辐亮度存在的基本关系为N＝π×Δλ×L。
5.一种采用如权利要求1-4任一项的有效波段宽度计算方法对HJ-1BB08在轨星上的定标方法，具体为：1)卫星发射后根据面源黑体温度与辐照度转换方程，计算星上定标黑体在特定温度Tc时，定标黑体发射的辐照度NTc；2)根据上述建立的查找表，计算星上定标黑体在温度Tc时对应的有效波段宽度，再根据辐照度和辐亮度存在的基本关系，得到定标黑体的辐亮度LTc；3)利用星上定标黑体在高温点(328±5K)和低温点(293±5K)两点的L定标黑体和DN定标黑体，基于L＝(DN-b)/g回归得到星上定标系数：增益g和偏移b；4)根据得到的增益g和偏移b，利用式L＝(DN-b)/g对遥感影像进行计算，得到卫星观测地面时的辐亮度L。

说明书全文

HJ-1B B08的有效波段宽度计算方法及定标方法

技术领域

[0001] 本发明涉及一种轨星上有效波段宽度计算及定标方法，尤指一种HJ-1B B08在轨星上的有效波段宽度计算及定标方法。

背景技术

[0002] 通道有效波段宽度是星上定标实施过程中的重要参数，现有计算有效波段宽度的方法均以传感器通道响应CRF为基础，其中主要有半高宽法和矩方法。半高宽法以光谱最大响应的一半高度(0.5Rmax)来表示，经验化且简单，如图1所示。矩方法通过中心波长来计算有效波段宽度，但是由于以上两种方法没有充分考虑影响有效波段宽度的因素，以及HJ-1B B08通道响应的特性，因此半高宽法和矩方法计算出的通道有效波段宽度往往不准确，从而影响定标结果，半高宽法过低计算了通道有效波段宽度，定标结果偏高，而矩方法定标结果偏低，半高宽法和矩方法并不适用于HJ-1B B08通道有效波段宽度的计算。发明内容

[0003] 针对现有技术存在的问题，本发明提供一种适用于HJ-1B B08通道有效波段宽度计算，且计算准确，能很好的应用于星上定标的有效波段宽度计算方法及其定标方法。 [0004] 为实现上述目的，本发明的HJ-1B B08在轨星上定标有效波段宽度的方法，具体为：1)HJ-1B B08通道在发射前根据实验室里面的标准面源黑体数据建立面源黑体温度与辐照度转换方程，计算在不同面源黑体温度Tm时的辐照度NTm；2)根据地面标准面源黑体下的比辐射率εm、温度Tm和通道响应函数CRF，利用Planck方程，构建有效辐亮度Lm与面源黑体温度Tm的方程，并计算在特定面源黑体温度Tm时面源黑体的有效辐亮度Lm；3)根据光学物理学中，辐照度和辐亮度存在的基本关系，构建基于地面实验室数据的有效波段宽度和面源黑体温度的查找表，得到不同黑体温度Tm条件下的有效波段宽度Δλ。 [0005] 进一步，所述HJ-1B卫星的面源黑体温度与辐照度转换方程为 k0，k1，k2，k3是已知常数。

[0006] 进一步，所述有效辐亮度 Lm 与面源黑体温度Tm的方程为式中h为Planck常数，k为玻耳兹曼常数，c为光速，f(λ)为传感器通道响应函数CRF，其中为Planck万程。

[0007] 进一步，辐照度和辐亮度存在的基本关系为N＝π×Δλ×L。

[0008] 本发明的HJ-1B B08在轨星上的定标方法，具体为：1)卫星发射后根据面源黑体温度与辐照度转换方程，计算星上定标黑体在特定温度Tc时，定标黑体发射的辐照度NTc；2)根据上述建立的查找表，计算星上定标黑体在温度Tc时对应的有效波段宽度，再根据辐照度和辐亮度存在的基本关系，得到定标黑体的辐亮度LTc；3)利用星上定标黑体在高温点(328±5K)和低温点(293±5K)两点的L定标黑体和DN定标黑体，基于L＝(DN-b)/g回归得到星上定标系数：增益g和偏移b；4)根据得到的增益g和偏移b，利用式L＝(DN-b)/g对遥感影像进行计算，得到卫星观测地面时的辐亮度L。

[0009] 本发明的HJ-1B B08的有效波段宽度计算方法，不仅考虑了HJ-1B B08的通道响应函数CRF，还充分考虑了其它可能影响到有效波段宽度值的参数，得到的有效波段宽度和面源黑体温度的查找表，有效波段宽度Δλ从星上定标低温点293±5K到高温点328±5K，由2.010近似线性变化到2.016，相对差异在0.3％以内，非常稳定。根据本发明的方法计算出的有效波段宽度进行定标，在低温点时，探元本身的记数差异在4DN左右，辐射差异在-2 -1 -10.02W·m ·μm ·sr 左右，探元之间的相对差异较大；在高温点时，探元本身的记数差异-2 -1 -1
在2DN左右，辐射差异在0.02W·m ·μm ·sr 左右，线性变化较为稳定，探元之间也存在一定的相对差异。为了降低记数差异和辐射差异的影响，本发明的定标方法对所有探元的响应作均值处理，分别得到低温点和高温点的辐亮度L及计数值DN，再通过线性回归，得到在不同有效波段宽度下的定标系数，定标准确，能很好的应用于星上定标。附图说明

[0010] 图1为半高宽法计算的有效波段宽度；

[0011] 图2为本发明的查找表法计算的有效波段宽度；

[0012] 图3为2009年8月5日星上定标高温点计数值；

[0013] 图4为2009年8月5日星上定标低温点计数值；

[0014] 图5为2009年8月14日星上定标高温点计数值；

[0015] 图6为2009年8月14日星上定标低温点计数值；

[0016] 图7为校正黑体辐射温度Tc与等效辐射率NTc对应曲线；

[0017] 图8为2009年8月5日星上定标高温点辐射率；

[0018] 图9为2009年8月5日星上定标低温点辐射率；

[0019] 图10为2009年8月14日星上定标高温点辐射率；

[0020] 图11为2009年8月14日星上定标低温点辐射率；

[0021] 图12为青海湖湖区HJ-1B B08遥感影像；

[0022] 图13为不同定标系数得到的水平剖面辐亮度；

[0023] 图14为不同定标系数得到的垂直剖面辐亮度。

具体实施方式

[0024] 以下对本发明的HJ-1B B08的有效波段宽度计算方法，具体说明如下： [0025] 1.式(1)是卫星在发射前根据实验室里面的标准面源黑体数据建立起来的温度与辐照度转换方程，利用式(1)计算在不同面源温度Tm时的辐照度NTm。 [0026] 2.按照式(2)，利用Planck方程，地面标准面源黑体比辐射率εm、温度Tm和通道响应函数CRF，如图1所示，卷积积分计算在特定温度Tm时面源黑体的有效辐亮度Lm。其中，面源黑体比辐射率εm＝0.998，通道响应函数CRF是传感器的属性，由实验室里测量得到。

[0027] 3.在光学物理学中，辐照度和辐亮度存在基本关系，如式(3)。因此，可通过NTm和πLm的比值构建基于地面实验室数据的Δλ(Tm)查找表，得到在不同黑体温度Tm条件下的Δλ。

[0028] 式(1)-2

[0029] 上式(1)中，NTm为辐照度，单位W·m 。Tm为黑体温度，单位K。k0，k1，k2，k3是已知常数，如

[0030] 表1所示。

[0031] 式(2)-2 -1 -1

[0032] 上式(2)中，Lm为辐亮度，单位W·m ·μm ·sr ，h为Planck常数，取值-34 -23 -16.626×10 ，单位J.s，k为玻耳兹曼常数，取值1.3806×10 ，单位J·K ，c为光速，取值
8 -1
2.998×10，单位m·s ，f(λ)为传感器通道响应CRF。其中为Planck方
程。

[0033] N＝π×Δλ×L 式(3)

[0034] 上式(3)中，Δλ为有效波段宽度，单位μm。

[0035] 表1HJ-1B B08通道辐射率N和温度T的反演系数

[0036]

[0037] 查找表法得到的有效波段宽度结果如图2和表2所示。从图2可知，Δλ从星上定标低温点293±5K到高温点328±5K，由2.010近似线性变化到2.016，相对差异在0.3％以内，非常稳定。

[0038] 表2查找表法计算的有效波段宽度

[0039]

[0040]

[0041] 以下介绍本发明的查找表计算星上有效波段宽度的定标应用： [0042] 1.卫星发射后，根据式(1)计算星上定标黑体在特定温度Tc时，定标黑体发射的辐照度NTc。

[0043] 2.根据建立的查找表，计算星上定标黑体在温度Tc时对应的有效波段宽度，再根据式(3)，得到定标黑体的辐亮度LTc。

[0044] 3.利用星上定标黑体在高温点(328±5K)和常温点(293±5K)两点的L定标黑体和DN定标黑体，基于式(4)回归得到星上定标系数：增益g和偏移b。

[0045] L＝(DN-b)/g 式(4)

[0046] 式(4)中DN为传感器计数值，即为影像记录值，b为定标系数的偏移，单位DN，g-1 2为定标系数的增益，单位DN·W ·m·sr·μm。

[0047] 4.根据得到的增益g和偏移b，利用式(4)对遥感影像进行计算，得到卫星观测地面时的辐亮度L。

[0048] 现有的方法均以传感器通道响应CRF为基础，进行有效波段宽度的计算，其中主要有半高宽法和矩方法。半高宽法以光谱最大响应的一半高度(0.5Rmax)来表示，经验化且简单，如图1所示。矩方法通过中心波长来计算有效波段宽度，如式(5)所示。表3为三种方法计算出的Δλ结果，从表中可知，半高宽法和矩方法计算出的结果差异较大，相对差异在20％左右。查找表法得到的Δλ结果介于半高宽法和矩方法得到的结果之间。 [0049] 式(5)

[0050] 上式(4)中，中λmax和λmin是波段的最大值和最小值，在这个波段范围之外传感器的响应函数是0，其它参数与上述参数一致。

[0051] 表3 半高宽法和矩方法计算HJ-1B B08有效波段宽度

[0052]

[0053]

[0054] 以2009年8月5日和2009年8月14日，星上定标黑体的记录数据为对象，利用三种不同的方法得到的有效波段宽度计算星上定标系数，并通过星上定标的精度来说明不同有效波段宽度值的优劣性。

[0055] HJ-1B B08波段共有10个扫描探元，每次星上定标时共进行225 帧测控，测控参数主要记录在黑体状态文件(rcp格式)和传感器状态文件(ca1格式)中。ca1文件记录了传感器对热辐射响应的计数值DN，定标时每帧每个探元分别记录8次DN值，对应于ca1文件的最后8列，计算时取8个DN的均值。图3和图4分别为2009年8月5日星上定标高温点和低温点计数值计算结果，图5和图6分别为2009年8月14日星上定标高温点和低温点计数值计算结果。

[0056] rcp文件包含了定标时校正黑体的温度Tc和双面镜左、右边缘温度TM_L和TM_R，根据由地面标准面源黑体建立的式(1)计算校正黑体在温度Tc时传感器经由内部光路系-2统接受到的等效辐射率NTc(单位：W·m )，其中ki为发射前实验室测量的常数，图7为根据ki计算出来的不同黑体温度Tc对应的辐射率NTc曲线。

[0057] B08波段通过双面镜观测成像，双面镜温度对内部光路系统具有辐射贡献，因此通过式(6)对双面镜的辐射贡献进行纠正，得到纠正后的辐照度Nc，式中aci和bci分别为双面镜温度纠正常数。图8和图9分别为2009年8月5日星上定标高温点和低温点辐照度计算结果，图10和图11分别为2009年8月14日星上定标高温点和低温点辐照度计算结果。

[0058] Nc＝ac+bc×NTc

[0059] ac＝ac0+ac1×TM

[0060] 式(6)

[0061] bc＝ac0+bc1×TM

[0062] TM＝(TM_L+TM_R)/2

[0063] 由图3至图6的计数值和图8至图11的辐射率值可知，HJ-1B B08传感器的10个探元本身及相互之间存在着记数差异和辐射差异。在低温点时，探元本身的记数差异在4DN-2 -1 -1左右，辐射差异在0.02W·m ·μm ·sr 左右，探元之间的相对差异较大；在高温点时，探-2 -1 -1
元本身的记数差异在2DN左右，辐射差异在0.02W·m ·μm ·sr 左右，线性变化较为稳定，探元之间也存在一定的相对差异。为了降低记数差异和辐射差异的影响，对所有探元的响应作均值处理，分别得到低温点和高温点的辐亮度L及计数值DN，再通过线性回归，得到在不同有效波段宽度下的定标系数，如表4所示。

[0064] 表48月5日和8月14日星上绝对辐射定标结果

[0065]

[0066] 由表4可知，半高宽法、矩方法和查找表法得到的定标系数相互之间存在差异，矩方法增益最大，半高宽法与查找表法相近，查找表法偏移最大，半高宽法和矩方法相近。以MODIS数据和星地同步观测实验数据比对不同定标系数得到的HJ-1B B08表观辐亮度，进而分析不同有效波段宽度计算方法的精度及适用性。

[0067] 利用表4中8月5日三组星上定标结果，分别定标2009年8月12日04点22分(UTM)过境于中国定标与真实性检验场地之一的青海湖湖区遥感影像，如图13所示，得到不同定标系数下的表观辐亮度L，并与同天04点18分(UTM)过境的MODIS传感器第31波段和32波段数据进行比较，图14和图15分别为MODIS B31，32数据和HJ-1B B08数据水平方向剖面序列和垂直方向剖面序列。由图14和图15可知，矩方法值低于所有其它值，半高宽法值和查找表法值基本介于MODIS B31，B32之间，但半高宽法在水体地物处明显高于MODIS值。结合图11中HJ-1B B08等效中心波长位于MODIS B31，32中间及热红外辐射传输理论可知，矩方法过高计算了通道有效波段宽度，定标结果明显偏低，误差最大。 [0068] 由于水体在大范围内温度变化不大且近似于黑体，因此是定标与真实性检验的倾向性目标之一。基于2009年8月12日卫星过境时刻在青海湖实验区(如图12所示)获取的星地同步观测水表参数和大气探空廓线，利用热红外辐射传输方程计算HJ-1B B08在卫星入瞳处的表观辐亮度，结果如图13和图14中“实验值”所示。表3为实验值与水体垂直剖面表观辐亮度均值的比对。结合图表可知，半高宽法过低计算了通道有效波段宽度，定标结果偏高，矩方法定标结果偏低，查找表法结果与实验值非常接近，表观辐亮度绝对误差-2 -1 -1为0.04W·m ·sr ·μm 。

[0069] 综上，通道有效波段宽度是星上定标实施过程中的重要参数，基于HJ-1B星上定标原理，利用发射前地面实验数据和Planck方程，构建有效波段宽度计算的查找表法，结果显示查找表法计算的有效波段宽度从低温点的2.010μm变化到高温点的2.016μm，非常稳定，介于半高宽法和矩方法得到的有效波段宽度之间。通过比对分析不同定标系数计算出的HJ-1B B08、MODIS B31，32及星地同步观测实验的表观辐亮度可知，半高宽法、矩方-2 -1 -1法和查找表法产生的表观辐亮度绝对误差分别为0.32、-0.81和0.04W·m ·sr ·μm ，对应的亮温差异分别在3K、-8K和0.4K左右。由于HJ-1B B08通道响应的特性，半高宽法和矩方法并不适用于HJ-1B B08通道有效波段宽度的计算，基于发射前地面实验数据构建的查找表法，能很好的应用于星上定标。

标题	发布/更新时间	阅读量
利用空气能，太阳能发电，海水淡化的方法及装置	2020-06-07	484
一种表征深层碳酸盐岩岩石物理学特征方法	2020-05-12	804
橡胶种子的贮藏技术	2020-05-27	344
利用郁金香生育期有效积温生产种球及控制花期的方法	2020-05-20	192
治疗和生产原料用甲型H1N1流感病毒灭活抗体血浆制备工艺	2020-05-28	1014
有益な効果を奏するアルギニンの送達	2020-05-21	273
Method of depositing a high aspect ratio structures	2020-06-01	913
Mechanism of refractory disease and treatment thereof	2020-06-16	417
Gel formulations and their use	2020-06-05	252
OPTICALLY NON-LINEAR SEMICONDUCTOR MATERIAL AND A METHOD FOR THE PRODUCTION THEREOF	2020-06-15	166

HJ-1BB08的有效波段宽度计算方法及定标方法

HJ-1B B08的有效波段宽度计算方法及定标方法

技术领域

背景技术

具体实施方式

该功能需要专业版企业版VIP权限，您可以：