留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

交叉定标中的不确定度分析及定标系数计算改进

高帅 李元 白廷柱 张玉香 郑小兵

高帅, 李元, 白廷柱, 张玉香, 郑小兵. 交叉定标中的不确定度分析及定标系数计算改进[J]. 中国光学. doi: 10.3788/CO.2019-0215
引用本文: 高帅, 李元, 白廷柱, 张玉香, 郑小兵. 交叉定标中的不确定度分析及定标系数计算改进[J]. 中国光学. doi: 10.3788/CO.2019-0215
GAO Shuai, LI Yuan, BAI Ting-zhu, ZHANG Yu-xiang, ZHENG Xiao-bing. Uncertainty Analysis in Cross-calibration and Optimization Calculation of Calibration Coefficients[J]. Chinese Optics. doi: 10.3788/CO.2019-0215
Citation: GAO Shuai, LI Yuan, BAI Ting-zhu, ZHANG Yu-xiang, ZHENG Xiao-bing. Uncertainty Analysis in Cross-calibration and Optimization Calculation of Calibration Coefficients[J]. Chinese Optics. doi: 10.3788/CO.2019-0215

交叉定标中的不确定度分析及定标系数计算改进

doi: 10.3788/CO.2019-0215
基金项目: 国家重点研发计划课题(No.2018YFB0504601);国家自然基金委员会资助项目(No.41271373)
详细信息
    作者简介:

    高帅(1995—),女,山西吕梁人,北京理工大学硕士研究生在读,研究方向为遥感器交叉定标。E-mail:arcbeacon@foxmail.com

    李元(1978—),女,北京丰台人,毕业于北京理工大学,工学博士,国家卫星气象中心副研究员,主要从事在轨遥感器定标与产品反演等方面的研究。E-mail:liyuan@cma.gov.cn

  • 中图分类号: TP732

Uncertainty Analysis in Cross-calibration and Optimization Calculation of Calibration Coefficients

Funds: Supported by National Key R&D Program of China (No.2018YFB0504601) and National Natural Science Foundation of China (No.41271373)
  • 摘要: 针对我国在研的太阳反射波段高光谱基准载荷对高精度辐射基准传递方法的需求,考虑一般交叉定标将时间、空间、观测几何、光谱匹配得到的近一致观测数据点对以普通最小二乘法的回归方式来计算定标系数,忽略了获取的各数据点对的质量差异,降低了定标系数的有效性问题,提出了基于不确定度计算的定标改进方法,利用不确定度分析方法计算数据点对中各辐射基准值的不确定度并给出权重系数,采用加权最小二乘法回归定标系数。选取与在研的基准载荷参数最为接近的HYPERION作为辐射基准替代载荷,分别利用普通最小二乘法和加权最小二乘法对MODIS CH1~7进行了定标,采用MODIS官方定标系数作为真值对定标结果进行验证。结果表明,采用加权最小二乘法对MODIS的1、2、4、5、6、7通道回归的定标系数更接近真值,定标结果的最大相对误差与传统方法相比降幅在3~5%范围内,平均相对误差降幅在0.5~1.5%范围内,说明本文的加权最小二乘法可进一步提升交叉定标的精度。
  • 图  1  交叉定标流程图

    Figure  1.  Flowchart of cross-calibration approach

    图  2  2012年7月10日HYPERION(左)和MODIS(右)过境Liyba4场地时的影像

    Figure  2.  Image of HYPERION(left) and MODIS(right) when they passed the Liyba4 site on 10th July, 2012

    图  3  Kref的概率分布图

    Figure  3.  The probability density distribution of Kref

    图  4  Kspace的概率分布图

    Figure  4.  The probability density distribution of Kspace

    图  5  Kspectral的概率分布图

    Figure  5.  The probability density distribution of Kspectral

    表  1  与基准载荷指标相近的高光谱成像仪参数

    Table  1.   Parameters of partial hyperspectral imagers

    载荷名称光谱覆盖
    范围(nm)
    光谱分辨率(nm)空间分辨率(m)幅宽(km)通道数
    HICO380~9605.7100
    HYPERION400~2 50010307.7242
    LAC930~1 5802~6250185256
    HIS450~950410050128
    下载: 导出CSV

    表  2  WLS和OLS对MODIS CH1~7通道的定标系数回归结果

    Table  2.   The calibration result of MODIS CH1~7 by WLS and OLS

    MODIS CH1234567
    β00.000 00.000 00.000 00.000 00.000 00.000 00.000 0
    β10.027 20.010 20.036 10.024 90.005 90.002 70.000 8
    βOLS 0−1.902 5−1.951 6−4.023 8−0.900 5−2.885 8−0.746 6−0.496 0
    βOLS 10.026 10.009 60.031 80.023 40.005 30.002 40.000 8
    βWLS 0−0.325 3−0.288 3−5.891 5−2.501 9−2.379 2−0.730 2−0.364 4
    βWLS 10.026 10.009 60.032 50.024 00.005 30.002 40.000 8
    ${{\rm{\epsilon}}_{{{\max}}\;{\rm{OLS}}}}$0.091 20.110 10.215 40.062 90.213 10.142 80.229 2
    ${{\rm{\epsilon}}_{{{\max}}\;{\rm{WLS}}}}$0.048 60.065 40.238 60.042 00.160 90.106 00.170 2
    εOLS0.058 90.076 10.164 70.068 90.149 90.126 50.071 6
    εWLS0.043 00.060 40.164 30.064 90.138 80.119 40.054 4
    RMSEOLS6.904 67.364 816.573 16.934 68.633 23.893 70.466 3
    RMSEWLS5.455 26.254 116.095 95.896 08.061 93.667 80.359 7
    下载: 导出CSV
  • [1] CHANDER G, HEWISON T J, FOX N, et al. Overview of intercalibration of satellite instruments[J]. IEEE Transactions on Geoscience and Remote Sensing, 2013, 51(3): 1056-1080. doi:  10.1109/TGRS.2012.2228654
    [2] LUKASHIN C, WIELICKI B A, YOUNG D F, et al. Uncertainty estimates for imager reference inter-calibration with CLARREO reflected solar spectrometer[J]. IEEE Transactions on Geoscience and Remote Sensing, 2013, 51(3): 1425-1436. doi:  10.1109/TGRS.2012.2233480
    [3] National Research Council. Earth Science and Applications from Space: National Imperatives for the Next Decade and Beyond[M]. Washington, DC, USA: National Academies Press, 2007.
    [4] CEOS, WMO, GSICS, et al.. SI-traceable space-based climate observing system workshop[R]. London, UK: NPL, 2019. (请核对修改是否正确)
    [5] GORROÑO J, BANKS A C, FOX N P, et al. Radiometric inter-sensor cross-calibration uncertainty using a traceable high accuracy reference hyperspectral imager[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2017, 130: 393-417. doi:  10.1016/j.isprsjprs.2017.07.002
    [6] ROITHMAYR C M, LUKASHIN C, SPETH P W, et al. CLARREO approach for reference intercalibration of reflected solar sensors: on-orbit data matching and sampling[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(10): 6762-6774. doi:  10.1109/TGRS.2014.2302397
    [7] 王玉鹏, 胡秀清, 王红睿, 等. 可在轨溯源的太阳反射波段光学遥感仪器辐射定标基准传递链路[J]. 光学 精密工程,2015,23(7):1807-1812. doi:  10.3788/OPE.20152307.1807

    WANG Y P, HU X Q, WANG H R, et al. Standard transfer chain for radiometric calibration of optical sensing instruments with traceability[J]. Optics and Precision Engineering, 2015, 23(7): 1807-1812. (in Chinese) doi:  10.3788/OPE.20152307.1807
    [8] 陈申玮, 徐娜, 戴铁, 等. 空间辐射基准传递不确定性的光谱敏感性分析[J]. 光学学报,2018,38(1):0128004. doi:  10.3788/AOS201838.0128004

    CHEN SH W, XU N, DAI T, et al. Sensitivity of intercalibration uncertainty on spectral sampling of space-based radiance standard[J]. Acta Optica Sinica, 2018, 38(1): 0128004. (in Chinese) doi:  10.3788/AOS201838.0128004
    [9] 赵维宁, 胡秀清, 方伟, 等. 卫星光学仪器辐射交互定标方法的应用和发展[J]. 光学 精密工程,2015,23(7):1921-1931. doi:  10.3788/OPE.20152307.1921

    ZHAO W N, HU X Q, FANG W, et al. Development and applications of intercalibration for satellite optical instruments[J]. Optics and Precision Engineering, 2015, 23(7): 1921-1931. (in Chinese) doi:  10.3788/OPE.20152307.1921
    [10] XIONG X X, ANGAL A, BUTLER J, et al. Global space-based inter-calibration system reflective solar calibration reference: from Aqua MODIS to S-NPP VIIRS[J]. Proceedings of SPIE, 2016, 9881: 98811D.
    [11] 徐文斌, 郑小兵, 易维宁. 基于超光谱成像仪Hyperion的交叉定标方法[J]. 光学学报,2013,33(5):0528002. doi:  10.3788/AOS201333.0528002

    XU W B, ZHENG X B, YI W N. Cross-calibration method based on hyperspectral imager Hyperion[J]. Acta Optica Sinica, 2013, 33(5): 0528002. (in Chinese) doi:  10.3788/AOS201333.0528002
    [12] HUNG K W, SIU W C. Improved image interpolation using bilateral filter for weighted least square estimation[C]. Proceedings of 2010 IEEE International Conference on Image Processing, IEEE, 2010.
    [13] PETRAKOS G, ARTELARIS P. European regional convergence revisited: a weighted least squares approach[J]. Growth and Change, 2009, 40(2): 314-331. doi:  10.1111/j.1468-2257.2009.00477.x
    [14] USGS. EO-1 HYPERION L1GST product[EB/OL]. USGS. https://earthexplorer.usgs.gov/. (请核对标题信息并且补充引用日期)
    [15] NASA. TERRA MODIS Level 1B Calibrated Radiances[EB/OL]. NASA EOSDIS LP DAAC. https://ladsweb.modaps.eosdis.nasa.gov/search/. (未找到本条文献信息, 请核对)
    [16] CHANDER G, HELDER D L, AARON D, et al. Assessment of spectral, misregistration, and spatial uncertainties inherent in the cross-calibration study[J]. IEEE Transactions on Geoscience and Remote Sensing, 2013, 51(3): 1282-1296. doi:  10.1109/TGRS.2012.2228008
    [17] JCGM. JCGM 101: 2008 Evaluation of measurement data – supplement 1 to the “guide to the expression of uncertainty in measurement” – propagation of distributions using a Monte Carlo method[S]. BIPM Joint Committee for Guides in Metrology, 2008. (未找到出版地信息, 请补充)
    [18] XIONG X, CHE N, XIE Y, et al. Four-years of on-orbit spectral characterization results for Aqua MODIS reflective solar bands[J]. Proceedings of SPIE, 2006, 6361: 63610S. doi:  10.1117/12.687163
  • [1] 迟明波, 韩欣欣, 徐阳, 舒风风, 吴一辉.  宽谱段高分辨扫描光谱定标技术 . 中国光学, doi: 10.3788/CO.20201302.0249
    [2] 孙德贝, 李志刚, 李福田.  用于太阳光谱仪的光电探测系统线性度测试装置 . 中国光学, doi: 10.3788/CO.20191202.0294
    [3] 毛靖华, 王咏梅, 石恩涛, 张仲谋, 江芳.  基于中阶梯光栅的波长定标方法研究 . 中国光学, doi: 10.3788/CO.20171003.0376
    [4] 穆治亚, 艾华, 樊孝贺, 何昕, 何丁龙, 韩冬松, 于国栋.  采用整体最小二乘法的条纹图配准方法 . 中国光学, doi: 10.3788/CO.20160906.0625
    [5] 孙景旭, 刘洪兴, 许艳军, 刘则洵, 李葆勇, 任建伟.  4m直径均匀扩展定标光源 . 中国光学, doi: 10.3788/CO.20150805.0823
    [6] 于海, 万秋华, 王树洁, 卢新然, 杜颖财.  编码器动态检测系统高实时性高精度角度基准设计 . 中国光学, doi: 10.3788/CO.20150803.0447
    [7] 李浩, 宋玲玲, 张立钧, 王焕然, 李亮, 王菲, 衣云骥, 张大明.  4通道交叉型二氧化硅光波导延迟线阵列的设计与制备 . 中国光学, doi: 10.3788/CO.20140703.0435
    [8] 王龙, 蔺超, 郑玉权.  CO2探测仪星上定标铝漫反射板的制备与试验 . 中国光学, doi: 10.3788/CO.20130604.0591
    [9] 薛旭成, 傅瑶, 韩诚山.  TDI CCD相机的卫星姿态稳定度确定 . 中国光学, doi: 10.3788/CO.20130605.0767
    [10] 陈健, 王伟国, 高慧斌, 刘廷霞, 吉桐伯, 于洪君.  紫外探测器的辐射定标及标准传递 . 中国光学, doi: 10.3788/CO.20120504.0423
    [11] 赵亚辉, 范长江.  组合式发光二极管路灯反射器的设计 . 中国光学, doi: 10.3788/CO.20120505.0520
    [12] 安岩, 孙强, 刘英, 李淳.  交叉型消像散Czerny-Turner结构光谱仪设计 . 中国光学, doi: 10.3788/CO.20120505.0470
    [13] 沈宏海, 黄猛, 李嘉全, 刘晶红, 戴明, 贾平.  国外先进航空光电载荷的进展与关键技术分析 . 中国光学, doi: 10.3788/CO.20120501.0020
    [14] 朱时雨, 张新, 李威.  计算机辅助装调与传统基准传递技术相结合实现三镜消像散系统的装调 . 中国光学,
    [15] 鲁剑锋.  无人机光电载荷图像处理器的设计 . 中国光学,
    [16] 曲艺.  真空紫外波段光栅二级衍射效率测量装置 . 中国光学,
    [17] 于杰.  用于相移点衍射干涉仪的加权最小二乘相位提取算法 . 中国光学,
    [18] 李东景, 于平, 李先锋, 赵昶宇.  星上定标系统探测器高精度信号处理电路的实现 . 中国光学,
    [19] 张军强, 颜昌翔.  碳/环氧复合材料在航天有效载荷支撑结构中的应用 . 中国光学,
    [20] 李晓晖, 颜昌翔.  成像光谱仪星上定标技术 . 中国光学,
  • 加载中
计量
  • 文章访问数:  231
  • HTML全文浏览量:  109
  • 被引次数: 0
出版历程
  • 网络出版日期:  2020-04-03

交叉定标中的不确定度分析及定标系数计算改进

doi: 10.3788/CO.2019-0215
    基金项目:  国家重点研发计划课题(No.2018YFB0504601);国家自然基金委员会资助项目(No.41271373)
    作者简介:

    高帅(1995—),女,山西吕梁人,北京理工大学硕士研究生在读,研究方向为遥感器交叉定标。E-mail:arcbeacon@foxmail.com

    李元(1978—),女,北京丰台人,毕业于北京理工大学,工学博士,国家卫星气象中心副研究员,主要从事在轨遥感器定标与产品反演等方面的研究。E-mail:liyuan@cma.gov.cn

  • 中图分类号: TP732

摘要: 针对我国在研的太阳反射波段高光谱基准载荷对高精度辐射基准传递方法的需求,考虑一般交叉定标将时间、空间、观测几何、光谱匹配得到的近一致观测数据点对以普通最小二乘法的回归方式来计算定标系数,忽略了获取的各数据点对的质量差异,降低了定标系数的有效性问题,提出了基于不确定度计算的定标改进方法,利用不确定度分析方法计算数据点对中各辐射基准值的不确定度并给出权重系数,采用加权最小二乘法回归定标系数。选取与在研的基准载荷参数最为接近的HYPERION作为辐射基准替代载荷,分别利用普通最小二乘法和加权最小二乘法对MODIS CH1~7进行了定标,采用MODIS官方定标系数作为真值对定标结果进行验证。结果表明,采用加权最小二乘法对MODIS的1、2、4、5、6、7通道回归的定标系数更接近真值,定标结果的最大相对误差与传统方法相比降幅在3~5%范围内,平均相对误差降幅在0.5~1.5%范围内,说明本文的加权最小二乘法可进一步提升交叉定标的精度。

English Abstract

高帅, 李元, 白廷柱, 张玉香, 郑小兵. 交叉定标中的不确定度分析及定标系数计算改进[J]. 中国光学. doi: 10.3788/CO.2019-0215
引用本文: 高帅, 李元, 白廷柱, 张玉香, 郑小兵. 交叉定标中的不确定度分析及定标系数计算改进[J]. 中国光学. doi: 10.3788/CO.2019-0215
GAO Shuai, LI Yuan, BAI Ting-zhu, ZHANG Yu-xiang, ZHENG Xiao-bing. Uncertainty Analysis in Cross-calibration and Optimization Calculation of Calibration Coefficients[J]. Chinese Optics. doi: 10.3788/CO.2019-0215
Citation: GAO Shuai, LI Yuan, BAI Ting-zhu, ZHANG Yu-xiang, ZHENG Xiao-bing. Uncertainty Analysis in Cross-calibration and Optimization Calculation of Calibration Coefficients[J]. Chinese Optics. doi: 10.3788/CO.2019-0215
    • 卫星遥感器在发射过程中周围环境变化剧烈导致其物理性能改变,在轨后随着时间的推移也会出现不同程度的衰减,因此需要对遥感器重新进行在轨定标以保证其观测数据的科学准确性[1]

      目前太阳反射波段(RSB)卫星载荷在轨定标精度在10−2量级[2],与观测气候变化所需的10−3量级的精度水平还存在一定差距[3],考虑到进一步提高现有定标方法精度的难度较大,各国均将发展空基溯源标准的计划列入日程,通过在太空建立辐射基准并采用交叉定标的方法将辐射基准进行传递,从而提高在轨定标精度水平。国际上已经有NASA主导的CLARREO计划和NPL主导的THRUTH计划[4][5][6]。我国也于近年来开展了国产高精度在轨辐射基准的相关研制计划[4][7][8]

      交叉定标方法是将在轨辐射基准传递至其他卫星的重要方法之一,是一种利用已有的定标精度较高的卫星作为参考对其他卫星进行定标的方法,其不受定标场景的限制,定标成本低,定标频次高[1][9]。目前国际上关于交叉定标的研究已经取得了诸多成果[1][9],针对红外波段(IR)已经建立起基于METOP IASI和AQUA AIRS的交叉定标系统[9],在RSB波段普遍以定标手段丰富且定标精度较高的AQUA MODIS作为辐射基准传递至其他卫星,但AQUA MODIS已在轨运行了16年,因此S-NPP VIIRS正成为RSB波段交叉定标基准的新选择[10]。综合可知目前RSB波段交叉定标能够达到的平均精度水平在5%左右[1][2][11]

      一般的交叉定标方法是将经过近一致观测筛选匹配后的数据点对代入回归模型中,然后使用普通最小二乘法(OLS,Ordinary Least Square)得出定标系数。而根据实际的交叉定标流程,各项匹配环节对基准辐射值影响的差异使得各数据点对的质量不同,因此若此时继续使用OLS回归定标系数将会降低结果的有效性[12][13]。针对以上问题,本文提出一种基于不确定度分析的交叉定标改进方法,利用加权最小二乘法(WLS,Weighted Least Square)回归定标系数,从而提高定标结果的准确性。

    • 考虑RSB波段高光谱成像仪作为参考仪器对宽通道成像仪进行交叉定标时的情况,分析交叉定标过程中的辐射误差传递源,通过蒙特卡罗法量化计算各数据点对中辐射基准值的不确定度,从而解算出权重系数矩阵,利用WLS方法回归定标系数,并通过与待定标仪器官方定标系数进行比对来检验WLS方法的有效性。

    • 因现阶段的基准载荷尚在研制中,缺乏可用数据,为提供具体的算法实现以及验证本文所提出WLS方法的有效性,本文通过比较现有的常见RSB波段高光谱载荷参数(表1),拟选用与基准载荷预定指标(光谱分辨率优于10 nm、光谱范围在380~2 350 nm内)最接近的EO-1卫星上搭载的HYPERION作为参考载荷实现具体的定标算法。

      表 1  与基准载荷指标相近的高光谱成像仪参数

      Table 1.  Parameters of partial hyperspectral imagers

      载荷名称光谱覆盖
      范围(nm)
      光谱分辨率(nm)空间分辨率(m)幅宽(km)通道数
      HICO380~9605.7100
      HYPERION400~2 50010307.7242
      LAC930~1 5802~6250185256
      HIS450~950410050128

      由于TERRA MODIS与HYPERION的轨道交点较多,易于满足交叉定标的时空匹配条件,且其定标精度高定标系数可以用于本文的算法验证,同时还具有产品数据易获取的优点,因此本文利用HYPERION作为参考仪器对TERRA MODIS 1~7通道进行定标,选取2012年7月4日起连续10天内HYPERION和MODIS的观测数据作为基础数据集[14][15],通过交叉定标各匹配环节获取满足近一致观测数据点对,从而建立起MODIS 1~7通道的回归模型。

    • 交叉定标流程如图1所示,通过时间、空间、观测几何、光谱匹配环节选取多个满足近一致观测的数据点对,建立起回归模型,通过线性回归获得待定标仪器定标系数。值得注意的是在使用高光谱遥感器作为参考仪器对宽通道式仪器进行定标时,直接采用光谱卷积的方式对其他宽通道式成像仪进行光谱匹配。

      图  1  交叉定标流程图

      Figure 1.  Flowchart of cross-calibration approach

      设筛选后满足近一致观测的数据点对共N个,将来自参考仪器和待定标仪器的数据分别用N×1的列向量LREFDNCAL表示,$\left({{L_{REF\;i}},D{N_{CAL\;i}}} \right)$代表一个数据点对,β0和β1为定标系数,将定标系数用2×1的列向量β表示,回归模型可用矩阵表示为:

      $$\begin{split} & {L_{REF}} = D{N_{CAL}}\beta \\ & \left[ {\begin{array}{*{20}{c}} {\begin{array}{*{20}{c}} {{L_{REF\;1}}}\\ {{L_{REF\;2}}}\\ \vdots \end{array}}\\ {{L_{REF\;N}}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {\begin{array}{*{20}{c}} 1&{D{N_{CAL\;1}}}\\ 1&{D{N_{CAL\;2}}}\\ \vdots & \vdots \end{array}}\\ {\begin{array}{*{20}{c}} 1&{D{N_{CAL\;N}}} \end{array}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {{\beta _0}}\\ {{\beta _1}} \end{array}} \right] \end{split}$$ (1)

      一般采用OLS对定标系数矩阵进行求解,得出待定标遥感器的定标系数。OLS通过寻找系数矩阵使得残差平方和最小来对模型参数进行估计[13],可以得出β的解为:

      $$\beta = {\left( {DN_{cal}^TD{N_{cal}}} \right)^{ - 1}}DN_{cal}^T{L_{ref}}$$ (2)

      考虑交叉定标回归模型,由于定标过程中各项匹配环节对不同数据点对中基准辐射值的影响不同,使得获取的各数据点对中辐射基准值的不确定度存在差异,即数据点对的质量存在差异,此时若继续使用OLS来回归定标系数将会降低定标系数的有效性,影响定标精度,由于交叉定标时各数据点对的获取过程相互独立,各数据点对相互之间没有影响,因此本文采用WLS法来计算定标系数,提高定标系数的有效性。

    • 加权最小二乘法通过给获取的每个交叉定标数据点对加上权重系数提高回归结果的有效性,该权重系数可以反映该定标数据点对的质量以及对定标曲线回归的影响大小,数据点对对应的权重系数一般用数据点对中自变量不确定度平方的倒数表示[12][13]。结合交叉定标回归模型,权重系数用数据点对中辐射基准值的不确定度平方的倒数表示,即获取的数据点对中辐射基准值的不确定度越大,说明数据点对的质量越差,对应权重系数越小,则该数据点对对定标系数回归的影响就越小,第i个定标数据点对对应的权重系数wi可表示为$1/{u^2}\left({{L_{ref\;i}}} \right)$。WLS法通过最小化广义残差平方和求出模型系数[12][13],可表示为:

      $$\beta = {\left( {DN_{cal}^TWD{N_{cal}}} \right)^{ - 1}}DN_{cal}^TW{L_{ref}}$$ (3)

      W为权重系数矩阵,结合定标模型,因此权重系数矩阵可以表示为:

      $$W = \left[ {\begin{array}{*{20}{c}} {\dfrac{1}{{{u^2}\left( {{L_{ref\;1}}} \right)}}}&0& \cdots &0\\ 0&{\dfrac{1}{{{u^2}\left( {{L_{ref\;2}}} \right)}}}& \vdots & \vdots \\ 0&0& \ddots &0\\ 0&0& \cdots &{\dfrac{1}{{{u^2}\left( {{L_{ref\;N}}} \right)}}} \end{array}} \right]$$ (4)

      在交叉定标时,要使用WLS计算定标系数,则需要计算定标数据点对中辐射基准值的不确定度,具体的计算方法如下文。

    • 在交叉定标数据点对的筛选过程中,Lref的误差主要来源于以下几点:参考仪器自身观测的不确定度、时相匹配中过境时间的差异、空间匹配时空间分辨率不一致和空间位置匹配不当产生的误差、地表方向性即观测几何匹配不当引起的误差、光谱匹配时光谱响应函数的衰减和波长漂移等引起的误差[2][16]

      在待定标仪器单个像元观测视场内,将与待定标仪器满足时间,空间位置,观测角度一致的像元输出值设为L0,经过近一致观测后配准的数据点对中辐射基准值表示为Lref,并将各环节不确定度对最终获取的辐射基准值的影响看作乘性影响。则可在L0基础上添加各项乘性因子(参考仪器自身Kref、时间匹配Ktime、空间匹配Kspace、观测几何匹配Kangle、光谱匹配Kspectral),Lref表示为:

      $${L_{ref}} = {L_0} \cdot {K_{ref}} \cdot {K_{time}} \cdot {K_{space}} \cdot {K_{angle}} \cdot {K_{spectral}}$$ (5)

      则根据不确定度计算的表达式法[17],辐射基准值的不确定度可以表示为:

      $$\begin{split} {{{\left( {\frac{{u\left( {{L_{ref}}} \right)}}{{{L_{ref}}}}} \right)}^2}}=& {{{\left( {\frac{{u\left( {{K_{ref}}} \right)}}{{{K_{ref}}}}} \right)}^2} + {{\left( {\frac{{u\left( {{K_{time}}} \right)}}{{{K_{time}}}}} \right)}^2} + {{\left( {\frac{{u\left( {{K_{space}}} \right)}}{{{K_{space}}}}} \right)}^2}}\\ & { + {{\left( {\frac{{u\left( {{K_{angle}}} \right)}}{{{K_{angle}}}}} \right)}^2} + {{\left( {\frac{{u\left( {{K_{spectral}}} \right)}}{{{K_{spectral}}}}} \right)}^2}} \end{split}$$ (6)

      公式(6)中$u\left({{K_{ref}}} \right)$$u\left({{K_{time}}} \right)$$u\left({{K_{space}}} \right)$$u\left({{K_{angle}}} \right)$$u\left({{K_{spectral}}} \right)$为对应各项K因子的不确定度,采用蒙特卡罗法(MCM)计算各K因子的不确定度。MCM法在原有输入数据的基础上经过添加扰动项重复试验得出输出量的概率分布后,给定置信概率根据概率分布进行统计得出置信区间从而得出输出量的不确定度[17]。利用MCM法各项K因子不确定度计算的具体描述如下:

      参考仪器自身观测不确定度引起Kref的不确定度,采用MCM法结合参考仪器自身观测值的不确定度在参考仪器观测值的基础上加适当的扰动,输入交叉定标流程,重复试验,最终得出Kref的分布,给定置信率从而得出Kref的不确定度$u\left({{K_{ref}}} \right)$

      对于时相匹配不当引起的Ktime因子的变化,主要考虑当地气象条件变化的影响,利用大气辐射传输模型,结合过境时间差内当地气象参数的变化,利用MCM法得出Ktime的分布,结合置信率给出置信区间,得出其不确定度$u\left({{K_{time}}} \right)$

      空间位置配准不当则采用滑动窗口法[16]进行分析,使得两卫星图像在原有的配准基础上发生相对位移,分析辐亮度值的变化,利用MCM法得出Kspace分布及其不确定度$u\left({{K_{space}}} \right)$

      观测几何匹配不当的影响主要受地表方向的影响,结合当地的BRDF模型和实际的观测几何角度的变化,输入辐射传输模型,利用MCM方法得出Kangle的分布及其不确定度$u\left({{K_{angle}}} \right)$

      光谱匹配主要考虑遥感器光谱响应的变化,通过在已知光谱响应函数的基础上适当地加扰动,MCM方法得出Kspectral的分布及其不确定度$u\left({{K_{spectral}}} \right)$

      在得出各项K因子的不确定度后,利用公式(6)得到辐射基准值的合成不确定度,针对每个数据点对都计算其辐射基准值的不确定度,从而得出权重系数矩阵用于WLS法回归模型系数。

    • 以HYPERION作为参考仪器对TERRA MODIS的1~7通道进行定标,验证WLS方法的有效性。为得出WLS法的权重系数矩阵,需要量化定标数据点对中辐射基准值的不确定度。以2012年7月10日卫星过境Libya4定标场地经纬度坐标(28.55°N,23.67°E)对应的数据点对为例,展示单个数据点对中辐射基准值的不确定度的获取过程。图2(a)(b)分别为HYPERION和MODIS过境Liyba4场地时的影像。

      图  2  2012年7月10日HYPERION(左)和MODIS(右)过境Liyba4场地时的影像

      Figure 2.  Image of HYPERION(left) and MODIS(right) when they passed the Liyba4 site on 10th July, 2012

      在前文的不确定度分析方法中,关于时间和观测几何匹配引入的不确定度计算需要获取地面配合测量的气象参数输入至大气传输模型,本文通过选取晴朗无云场景的观测图像以及设置尽量小的过境时间差和观测几何角度差阈值,减小时相匹配和观测几何匹配对数据点对中基准辐射值的影响,针对所有的数据点对中的辐射基准值忽略了时相和观测几何匹配引入的误差,使KtimeKangle的不确定度近似为零。

      已知HYPERION自身观测相对不确定度为3.4%[11];认为MODIS和HYPERION的空间位置匹配误差不超过单个MODIS像元大小即空间匹配的误差在1 km内;MODIS在轨的光谱响应的波长漂移在−2~2 nm范围内,响应度变化在5%范围内[18];(对获取的所有数据点对均使用以上条件)采用蒙特卡罗法分别计算KrefKspaceKspectral的分布,求出置信率95%时对应的置信区间,求得各因子对应的不确定度,忽略时间和观测几何引入的不确定度计算合成不确定度。针对MODIS CH01通道Libya4定标场获取的数据点对的辐射基准值计算其不确定度时,获得的各因子的分布如下图3(a)(b)(c)所示。

      图  3  Kref的概率分布图

      Figure 3.  The probability density distribution of Kref

      图  4  Kspace的概率分布图

      Figure 4.  The probability density distribution of Kspace

      图  5  Kspectral的概率分布图

      Figure 5.  The probability density distribution of Kspectral

      在该数据点对不确定度的量化过程中,由空间匹配不当引入的不确定度占比最大,Kspace的相对不确定度为1.87%,其次是遥感器自身观测精度的影响,Kref的相对不确定度为1.58%,由光谱匹配不当引入的不确定度占比最小,Kspectral的相对不确定度为0.067%,图3(c)可以看出由光谱匹配不当引起的Kspectral概率分布呈偏斜状,是由MODIS CH1的光谱响应函数和HYPERION观测输出的高光谱辐亮度曲线线型共同作用导致的。针对不同场景获得的高光谱的辐亮度曲线不同,在光谱匹配不确定度分析时获得的Kspectral概率分布图也不尽相同。最终计算得出MODIS CH01通道在Libya4定标场获取的数据点对中辐射基准值总的相对不确定度为2.45%。

      对所有的数据点对重复上述计算得出每个数据点对中辐射基准值的不确定度,从而得出权重系数矩阵,利用WLS和OLS分别计算HYPERION对MODIS CH 1~7通道的定标系数βWLSβOLS,并与MODIS产品数据中的定标系数β进行比对。为评价两种方法计算出的定标系数的定标效果,本文利用WLS和OLS定标系数对同一DN值序列进行定标并与官方定标系数的定标结果进行比对,分别计算定标结果与官方定标结果的平均相对误差ε以及均方根误差RMSE。设DN值序列共包含N个点,对DN值序列中第i个DN值采用官方系数定标后的辐亮度值为${L_{0\;i}}$,经计算获取的系数定标后的辐亮度值为Li,则相对差异计算公式为:

      $$\epsilon_i = \frac{{\left| {{L_i} - {L_{0\;i}}} \right|}}{{{L_{0\;i}}}}$$ (7)

      平均相对误差计算公式为:

      $$\varepsilon = \frac{1}{N}\mathop \sum \nolimits{\epsilon_i}$$ (8)

      均方根误差计算公式为:

      $$RMSE = \sqrt {\frac{1}{N}\sum {{\left( {{L_i} - {L_{0\;i}}} \right)}^2}} $$ (9)

      以HYPERION对MODIS CH01通道定标为例,采用WLS和OLS的回归曲线如图4所示,WLS和OLS计算的定标系数分别为(-1.9025,0.0261)、(-0.3253,0.0261),官方定标系数为(0.0,0.0272),可以看出WLS的定标系数要更接近官方定标系数,同时可以发现两种方法回归得到的定标系数的差异主要体现在截距β0上。将WLS和OLS得到的定标系数对同一DN值序列进行定标并计算与官方定标结果的相对差异如图5所示。从两种回归方法获得的定标系数与官方定标系数定标比较结果来看,采用WLS计算的定标系数的定标相对误差在低端明显减少,最大由9.12%降低到4.86%,平均相对误差由5.89%降低到4.30%,均方根误差由6.905 W/(m2∙sr∙μm)降低至5.455 W/(m2∙sr∙μm),说明WLS方法有效提高了该通道的定标曲线拟合精度。

      HYPERION对MODIS CH1~7通道的WLS和OLS计算定标系数结果βWLSβOLS、官方定标系数β、定标评价参数最大相对误差${\epsilon_{{\rm{max}}}}$、平均相对误差ε以及均方根误差RMSE的计算结果如下表2所示。

      表 2  WLS和OLS对MODIS CH1~7通道的定标系数回归结果

      Table 2.  The calibration result of MODIS CH1~7 by WLS and OLS

      MODIS CH1234567
      β00.000 00.000 00.000 00.000 00.000 00.000 00.000 0
      β10.027 20.010 20.036 10.024 90.005 90.002 70.000 8
      βOLS 0−1.902 5−1.951 6−4.023 8−0.900 5−2.885 8−0.746 6−0.496 0
      βOLS 10.026 10.009 60.031 80.023 40.005 30.002 40.000 8
      βWLS 0−0.325 3−0.288 3−5.891 5−2.501 9−2.379 2−0.730 2−0.364 4
      βWLS 10.026 10.009 60.032 50.024 00.005 30.002 40.000 8
      ${{\rm{\epsilon}}_{{{\max}}\;{\rm{OLS}}}}$0.091 20.110 10.215 40.062 90.213 10.142 80.229 2
      ${{\rm{\epsilon}}_{{{\max}}\;{\rm{WLS}}}}$0.048 60.065 40.238 60.042 00.160 90.106 00.170 2
      εOLS0.058 90.076 10.164 70.068 90.149 90.126 50.071 6
      εWLS0.043 00.060 40.164 30.064 90.138 80.119 40.054 4
      RMSEOLS6.904 67.364 816.573 16.934 68.633 23.893 70.466 3
      RMSEWLS5.455 26.254 116.095 95.896 08.061 93.667 80.359 7

      表2数据中可以看出除MODIS CH3通道WLS和OLS的定标评价参数相近外,其余通道WLS的定标效果均要好于OLS的定标效果。不同的MODIS通道在进行光谱匹配后获取的定标数据点对不同,各项K因子的不确定度计算结果也存在差异,因此WLS方法获取定标结果也不尽相同。对于MODIS的1、2、4、5、6、7通道WLS方法回归得到的定标系数与OLS相比,其定标平均相对误差的降幅在0.5~1.5%范围内,在观测低亮度目标时的定标相对误差降幅在3~5%范围内。针对MODIS CH3通道WLS方法定标效果不理想的分析原因如下:MODIS CH3通道带宽较窄(19 nm左右)仅对应2个高光谱通道,相对于其它通道可对应4~9个高光谱通道而言,CH3通道采用光谱卷积得出辐射基准值的准确性较差,各项K因子的不确定度计算准确度下降权重系数计算不理想,导致最终WLS回归结果准确性下降。同时可以发现采用WLS和OLS回归的系数与官方系数相比,斜率的差异不明显,主要体现在截距的差异上。亮度较低的目标一般存在较大的不确定度,WLS方法的引进可以有效的修正低亮度目标对截距的干扰,使定标结果更接近真实值。还需要指出的是HYPERION载荷自身的定标精度低于MODIS载荷的定标精度,因此WLS计算所得的MODIS定标系数与官方定标系数还存在一定差异,待我国自研的高精度太阳反射波段高光谱基准载荷在轨后使用其作为参考仪器向MODIS传递溯源基准,将有望消减该差异从而更准确地检验WLS方法的优势。

    • 本文提出的基于不确定度分析的加权最小二乘法交叉定标方法,在得出近一致观测的数据点对后,利用蒙特卡罗法对各项匹配环节对辐射基准值的影响进行量化,针对每个数据点对都可算出其辐射基准值的不确定度,利用不确定度平方的倒数作为权重系数,采用WLS方法回归了定标系数。以HYPERION作为参考基准载荷对MODIS CH1~7通道的定标结果表明,1、2、4、5、6、7通道采用WLS法回归的定标系数的评价参数要优于OLS法,其中定标系数的定标平均相对误差的降幅在0.5~1.5%范围内,在观测低亮度目标时的定标相对误差降幅在3~5%范围内,说明WLS方法有效地提高了定标系数的精度,验证了基于不确定度的WLS方法的有效性。

      在后续的工作中将进一步评估通道带宽较窄,对应高光谱通道有限时,基于高光谱地物谱型库的光谱匹配方法对应的不确定度,同时在计算不确定度时进一步分析时间匹配和观测几何匹配引入的不确定度,从而进一步提高权重系数和定标数据点对的计算精度,提升WLS的回归效果。

      致谢:本文由国家重点研发计划课题(2018YFB0504601)以及国家自然基金委员会面上项目(41271373)资助。感谢ladsweb.modaps.eosdis.nasa.gov提供MODIS数据及earthexplorer.usgs.gov提供HYPERION数据。

WeChat 关注分享

返回顶部

目录

    /

    返回文章
    返回