基于MODIS的GF-4/PMS遥感器交叉定标-以巴丹吉林沙漠为参考目标

基于MODIS的GF-4/PMS遥感器交叉定标-以巴丹吉林沙漠为参考目标 Cross calibration of GF-4/PMS Based on MODIS over Badain Jaran Desert 1 first-author 张浩，主要从事光学遥感辐射定标及大气校正； 张浩，主要从事光学遥感辐射定标及大气校正； 1 2 3 corresp 闫东川，研究方向海量遥感数据智能化处理与分析，E-mail： yandc@radi.ac.cn. 闫东川，研究方向海量遥感数据智能化处理与分析，E-mail： yandc@radi.ac.cn. yandc@radi.ac.cn 4 1 5 中国科学院空天信息创新研究院, 北京 100094 中国科学院空天信息创新研究院, 北京 100094 Aerospace Information Research Institute(AIR) Chinese Academy of Sciences(CAS), Beijing 100094, China Aerospace Information Research Institute(AIR) Chinese Academy of Sciences(CAS), Beijing 100094, China 中科星图空间技术有限公司, 西安 710100 中科星图空间技术有限公司, 西安 710100 Zhongke Xingtu Space Technology Co., Ltd, Xi'An 710100, China Zhongke Xingtu Space Technology Co., Ltd, Xi'An 710100, China 中国冶金地质总局矿产资源研究院, 北京 101300 中国冶金地质总局矿产资源研究院, 北京 101300 Institute of Mineral Resources Research, China Metallurgical Geology Bureau, Beijing101300, China Institute of Mineral Resources Research, China Metallurgical Geology Bureau, Beijing101300, China 中国矿业大学(北京), 北京 100083 中国矿业大学(北京), 北京 100083 China University of mining and technology(Beijing), Beijing 100101, China China University of mining and technology(Beijing), Beijing 100101, China 河南理工大学 测绘与国土信息工程学院, 河南 焦作 454003 河南理工大学 测绘与国土信息工程学院, 河南 焦作 454003 School of Surveying and Land Information Engineering, Henan Polytechnic University, Jiaozuo 454003, China School of Surveying and Land Information Engineering, Henan Polytechnic University, Jiaozuo 454003, China GF-4是我国首颗高分辨率静止轨道光学遥感卫星，搭载了一台可见/中波红外凝视相机，由于缺少星上定标系统导致无法提供高频次的星上辐射定标系数。为了更好的实现GF-4/PMS传感器辐射性能监测和数据定量化应用，本文选择了巴丹吉林沙漠南端均匀沙地，利用MODIS作为参考传感器对GF-4/PMS可见近红外波段进行交叉定标。为尽量减少地表不均匀性造成的交叉定标误差，使用较高分辨率的影像（Landsat-8/OLI）搜索巴丹吉林沙漠地区的均匀区域，确定为交叉定标试验区。为进一步减少遥感器成像角度和大气变化的影响，进行了MODIS成像角与GF-4成像角角度差约束（小于20°）、成像时刻约束（小于2小时）、云和成像质量问题，共选择2016年至2018年之间有效图像对13组，通过利用MODTRAN模拟计算二者光谱匹配因子进而实现交叉定标系数计算。结果表明：1） 本文计算的辐射定标系数在朗伯体假定和Ross-Li模型假定情况下具有较高一致性，忽略BRDF效应情况造成的计算误差在0.8%~4%之间，定标结果的不确定度在7.4%以内；2） 2016年至2018年GF-4/PMS遥感器辐射性能呈现缓慢下降特征，年衰减率小于1%。本文的方法可以有效地提高GF-4/PMS的辐射定标频率和定标精度，实现对GF-4/PMS传感器整个生命周期的辐射性能监测。 GF-4 is the first high-resolution geostationary orbit optical remote sensing satellite in China, equipped with a staring camera covering visible to mid-infrared spectrum. It is pretty to provide on-oribit radiometric calibration coefficients frequently due to lacking on-board calibration device. In order to better monitor the radiometric performance of GF-4/PMS and provide the solid basis for quantitative applications, the cross-calibration method was applied to retrieve the calibration coefficients over a select calibration site located in the north of Badan Jaran desert. The cross-calibration site was identified by searching the uniform area in the relative high resolution image (i.e., Landsat-8 /OLI) to minimize the uniform impact on the cross calibration results. To reduce the influence of imaging angles of sensor and atmospheric changing, the restrain was set to the imaging angle difference between MODIS and GF-4 (less than 20°), the imaging time difference (less than 2h), the cloud, and the image quality. Totally, 13 image pairs were available from 2016 to 2018 and were used to calculate the cross calibration coefficients after compensating the spectral matching factors, which were simulated by MODTRAN. The result shows that: 1) the cross calibration coefficients calculated under the Lambertian assumption are highly consistent with that under the Ross-Li BRDF assumption, with uncertainty less than 7.4%, and the relative difference ranges between 0.8% and 4% if the BRDF neglected; and 2) the radiometric performance of GF-4/PMS is slowly decreasing from 2016 to 2018, 1% per year. Thus, the proposed method can effectively improve the radiometric calibration frequency of GF-4/PMS with an acceptable accuracy, and can serve the radiometric performance monitoring over the whole life cycle of the sensor. 交叉定标 GF-4/PMS MODIS 不确定性分析 cross calibration GF-4/PMS MODIS uncertainty analysis 1　引 言 1 高分四号（GF-4）卫星于2015年12月29日在西昌卫星发射中心发射成功，是我国第一颗地球同步轨道遥感卫星，搭载了一台可见光50米/中波红外400米分辨率、大于400KM幅宽的凝视相机，采用面阵凝视方式成像，具备可见光、多光谱和红外成像能力，设计寿命8年，重访周期达20s，通过指向控制，实现对中国及周边地区的观测。 GF-4为减灾、林业、地震和气象等各种应用提供快速、可靠、稳定的光学遥感数据，为灾害风险预报、森林火灾监测、地震构造信息增添新的技术方法，在灾害监测等方面具有巨大的潜力和广阔的应用空间（ 聂娟 等，2018 聂娟 等，2018 ； 张磊 等，2018 张磊 等，2018 ； 吴玮 等，2019 吴玮 等，2019 ）。遗憾的是GF-4没有星上定标系统，传感器的定标完全依赖于场地定标，而场地定标受限于场地、设备、成本等因素，无法满足高频率定标需求，这限制了它的应用（ 高海亮 等，2010 高海亮 等，2010 ； Zhang 等， 2018 Zhang 等， 2018 ）。交叉定标是利用定标精度较高的传感器来标定待标定的卫星传感器的方法，该方法不需要高成本的野外同步试验，也不需要很精确的大气参数测量；同时交叉定标方法可以对历史数据进行标定，是目前最具有广泛应用前景的定标方法之一（ 赵维宁 等，2015 赵维宁 等，2015 ）。 基于交叉定标的优势，国内外学者先后将交叉定标方法应用于NOAA、北京一号、HJ-1、GF-1等卫星遥感器。 Teillet等（1990） Teillet等（1990） 基于美国White Sands场，针对NOAA-10/AVHRR，利用Landsat-5和SPOT/HVR影像，进行了AVHRR的交叉定标。 Teillet等（2007） Teillet等（2007） 分析了不同地物光谱对交叉定标光谱匹配因子的影响，得出利用Railroad试验场交叉定标精度优于利用草地交叉定标。 Rao等（2001） Rao等（2001） 基于Sonoran沙漠试验场，以NOAA-14/AVHRR为参考传感器，对GOES-8传感器进行交叉定标。 Cao（2005） Cao（2005） 提出SNO（Simultaneous Nadir Overpass）交叉定标方法，该方法利用极地地区的卫星影像对完成了NOAA系列极轨卫星15、16、17号搭载的HIRS传感器之间的相互交叉定标。 杨忠东等（2004） 杨忠东等（2004） 使用Landsat-7/ETM+为参考传感器，通过辐射传输模拟和统计分析，实现了对CBERS-01（中巴地球资源卫星）CCD相机的交叉定标。 陈正超（2008） 陈正超（2008） 等在缺少光谱响应函数的情况下，利用SPOT4/HRVIR2、Landst-5/TM和Terra/MODIS三种传感对北京一号小卫星进行交叉辐射定标，得到比较可信的定标系数。 马晓红等（2011） 马晓红等（2011） 选取阿拉伯沙漠作为定标试验场，以MODIS为参考传感器对HJ-1/CCD相机进行交叉辐射定标，结果表明交叉辐射定标的系数具有较高的精度。Zhong和Yang等（2014、2015、2017）等人选取巴丹吉林沙漠，基于分辨率较高的Landsat-7/ETM+、Landsat-8/OL和DEM数据构建地表BRDF模型，实现对HJ-1/CCD、GF-1/WFV和GF-4/PMS的交叉辐射定标。 张玉环等（2016） 张玉环等（2016） 基于反射率低、中、高（植被区、沙漠和雪）场景，使用MODIS作为参考传感器对GOCI进行交叉辐射定标。 Chen等（2017） Chen等（2017） 将交叉辐射定标作为最优逼近问题，以Landsat-8/OLI作为参考传感器，采用SCE-UA算法，通过迭代方程找到最优定标系数和BRDF调整因子，结果表明交叉辐射定标系数计算的地表反射率误差小于5%。 为减少观测角度差异造成的定标误差，通常交叉定标选取两个传感器观测角较小且较为接近的影像进行交叉定标。GF-4/PMS具有较宽的视角，由于纬度原因在中国的大部分地区都有较大的观测角，这对GF-4/PMS交叉定标带来较大困难。本文通过影像搜索巴丹吉林沙漠均匀场地作为交叉定标试验场，以MODIS传感器为参考，利用MODTRAN辐射传输模型，在约束MODIS与GF-4成像角度差异和二者过境时间差异条件下，对GF-4/PMS时序数据进行交叉定标，获取时序的交叉定标系数，实现对GF-4/PMS传感器性能的时序监测和评估。 2　研究方法 1 本文通过较高分辨率影像（Landsat-8/OLI）数据搜索均匀区域对大角度观测的中高分辨率数据（GF-4/PMS）与Terra（Aqua）/MODIS数据进行交叉定标。首先，使用几何和辐射校准较好的Landsat-8/OLI数据寻找合适的均匀区域作为定标试验场；然后，根据选取的试验场选取MODIS和GF-4/PMS时间序列数据；接着，时序数据筛选，除去云和成像质量较差的数据后，利用如下3个条件进一步筛选：1）550nm气溶胶光学厚度小于0.3；2）MODIS成像散射角度（观测方向和太阳入射方向夹角）和GF-4成像散射角度差异小于20度；3）MODIS和GF-4成像时刻小于2小时；最后，对时间序列数据进行预处理，提取试验场时间序列数据的观测几何等信息，使用辐射传输模型计算MODIS和GF-4/PMS模拟的表观辐亮度，计算光谱匹配因子；最后对GF-4/PMS进行交叉定标，总体流程如 图1 图1 所示。 2.1　定标原理 2 交叉定标是利用定标精度较高的传感器作为参考，对待定标的传感器进行定标；其原理是选取对同一目标成像的同步或近似同步的影像对，在分析两个传感器光谱响应、观测几何、大气参数等匹配的基础上，建立两个传感器图像数字计数值之间的关系，利用参考传感器已知的辐射定标系数求解待标定传感器的定标系数（ 高海亮 等，2010 高海亮 等，2010 ；吕文博，2014）。通常情况下，卫星遥感器的DN值与入瞳辐亮度存在线性关系： <math display="block" id="M1"><mi>L</mi><mo>=</mo><mi>g</mi><mi>a</mi><mi>i</mi><mi>n</mi><mi mathvariant="normal">*</mi><mi>D</mi><mi>N</mi><mo>+</mo><mi>o</mi><mi>f</mi><mi>f</mi><mi>s</mi><mi>e</mi><mi>t</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463095&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463103&type= 31.58066559 3.47133350 式中： <math id="M2"><mi>g</mi><mi>a</mi><mi>i</mi><mi>n</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463105&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463104&type= 6.60400009 3.47133350 和 <math id="M3"><mi>o</mi><mi>f</mi><mi>f</mi><mi>s</mi><mi>e</mi><mi>t</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463107&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463106&type= 7.36600018 3.47133350 分别是定标系数的增益和截距， <math id="M4"><mi>L</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463116&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463108&type= 1.86266661 2.62466669 表示传感器的入瞳辐亮度，DN表示图像的数字计数值。 通过两个传感器的光谱匹配因子，计算得到待标定传感器的表观辐亮度的模拟值： <math display="block" id="M5"><msubsup><mrow><mi>L</mi></mrow><mrow><mi>P</mi><mi>M</mi><mi>S</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msubsup><mo>=</mo><mi>Q</mi><mi mathvariant="normal">*</mi><msub><mrow><mi>L</mi></mrow><mrow><mi>M</mi><mi>O</mi><mi>D</mi><mi>I</mi><mi>S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463109&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463118&type= 22.35199928 4.82600021 式中， <math id="M6"><msubsup><mrow><mi>L</mi></mrow><mrow><mi>P</mi><mi>M</mi><mi>S</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463121&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463120&type= 5.75733376 4.06400013 为通过光谱匹配因子和MODIS入瞳辐亮度计算的GF-4/PMS的入瞳辐亮度， <math id="M7"><msub><mrow><mi>L</mi></mrow><mrow><mi>M</mi><mi>O</mi><mi>D</mi><mi>I</mi><mi>S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463125&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463123&type= 7.78933382 3.72533321 为从传感器测量得到的MODIS入瞳辐亮度， <math id="M8"><mi>Q</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463128&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463127&type= 2.53999996 3.38666677 为两个传感器的光谱匹配因子，光谱匹配因子可有下式计算： <math display="block" id="M9"><mi>Q</mi><mo>=</mo><mfrac><mrow><msubsup><mrow><mi>L</mi></mrow><mrow><mi>P</mi><mi>M</mi><mi>S</mi></mrow><mrow><mi>s</mi><mi>i</mi><mi>m</mi><mi>u</mi></mrow></msubsup><mo stretchy="false">(</mo><msub><mrow><mi>θ</mi></mrow><mrow><mi>s</mi><mo>-</mo><mi>P</mi><mi>M</mi><mi>S</mi></mrow></msub><mo>,</mo><msub><mrow><mi>θ</mi></mrow><mrow><mi>v</mi><mo>-</mo><mi>P</mi><mi>M</mi><mi>S</mi></mrow></msub><mo>,</mo><msub><mrow><mi>φ</mi></mrow><mrow><mi>P</mi><mi>M</mi><mi>S</mi></mrow></msub><mo stretchy="false">)</mo></mrow><mrow><msubsup><mrow><mi>L</mi></mrow><mrow><mi>M</mi><mi>O</mi><mi>D</mi><mi>I</mi><mi>S</mi></mrow><mrow><mi>s</mi><mi>i</mi><mi>m</mi><mi>u</mi></mrow></msubsup><mo stretchy="false">(</mo><msub><mrow><mi>θ</mi></mrow><mrow><mi>s</mi><mo>-</mo><mi>M</mi><mi>O</mi><mi>D</mi><mi>I</mi><mi>S</mi></mrow></msub><mo>,</mo><msub><mrow><mi>θ</mi></mrow><mrow><mi>v</mi><mo>-</mo><mi>M</mi><mi>O</mi><mi>D</mi><mi>I</mi><mi>S</mi></mrow></msub><mo>,</mo><msub><mrow><mi>φ</mi></mrow><mrow><mi>M</mi><mi>O</mi><mi>D</mi><mi>I</mi><mi>S</mi></mrow></msub><mo stretchy="false">)</mo></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463143&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463112&type= 50.88466644 9.22866726 式中， <math id="M10"><msubsup><mrow><mi>L</mi></mrow><mrow><mi>P</mi><mi>M</mi><mi>S</mi></mrow><mrow><mi>s</mi><mi>i</mi><mi>m</mi><mi>u</mi></mrow></msubsup><mi mathvariant="normal">、</mi><msubsup><mrow><mi>L</mi></mrow><mrow><mi>M</mi><mi>O</mi><mi>D</mi><mi>I</mi><mi>S</mi></mrow><mrow><mi>s</mi><mi>i</mi><mi>m</mi><mi>u</mi></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463145&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463144&type= 15.66333294 4.40266657 分别为MODTRAN模拟的PMS和MODIS各波段在实际成像角度下的入瞳辐亮度，其中 <math id="M11"><msub><mrow><mi>θ</mi></mrow><mrow><mi mathvariant="normal">s</mi><mo>-</mo><mi>P</mi><mi>M</mi><mi>S</mi></mrow></msub><mi mathvariant="normal">、</mi><msub><mrow><mi>θ</mi></mrow><mrow><mi>s</mi><mo>-</mo><mi>P</mi><mi>M</mi><mi>S</mi></mrow></msub><mi mathvariant="normal">、</mi><msub><mrow><mi>φ</mi></mrow><mrow><mi>P</mi><mi>M</mi><mi>S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463146&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463133&type= 27.60133362 4.48733330 分别对应GF-4的太阳天顶角、观测天顶角和相对方位角， <math id="M12"><msub><mrow><mi>θ</mi></mrow><mrow><mi mathvariant="normal">s</mi><mo>-</mo><mi>M</mi><mi>O</mi><mi>D</mi><mi>I</mi><mi>S</mi></mrow></msub><mi mathvariant="normal">、</mi><msub><mrow><mi>θ</mi></mrow><mrow><mi>v</mi><mo>-</mo><mi>M</mi><mi>O</mi><mi>D</mi><mi>I</mi><mi>S</mi></mrow></msub><mi mathvariant="normal">、</mi><msub><mrow><mi>φ</mi></mrow><mrow><mi>M</mi><mi>O</mi><mi>D</mi><mi>I</mi><mi>S</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463150&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463148&type= 33.95133591 4.48733330 分别对应MODIS的太阳天顶角、观测天顶角和相对方位角。利用 公式（1） 公式（1） 获得待标定传感器的辐亮度之后，便可结合待标定影像的DN值，通过下式拟合计算待定标传感器的定标系数，实现待定标传感器的标定。 <math display="block" id="M13"><msubsup><mrow><mi>L</mi></mrow><mrow><mi>P</mi><mi>M</mi><mi>S</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msubsup><mo>=</mo><mi>g</mi><mi>a</mi><mi>i</mi><msup><mrow><mi>n</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup><mi mathvariant="normal">*</mi><mi>D</mi><msub><mrow><mi>N</mi></mrow><mrow><mi>P</mi><mi>M</mi><mi>S</mi></mrow></msub><mo>+</mo><mi>o</mi><mi>f</mi><mi>f</mi><mi>s</mi><mi>e</mi><msup><mrow><mi>t</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463153&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463151&type= 41.65599823 4.82600021 2.2　试验场 2 为了选取合适的实验场，本文使用巴丹吉林沙漠地区无云的Landsat-8/OLI影像，位置如 图2 图2 所示，搜索DN值的均值和方差最小的500m×500m大小区域作为试验场，最终选取了巴丹吉林沙漠南部边缘的500m×500m的均匀区域（ 图2 图2 ） 图2（a） 图2（a） 为试验场的位置，其中右侧部分为Landsat-8/OLI真彩色图像； 图2（b） 图2（b） 为GF-4/PMS影像上试验场的位置； 图2（c） 图2（c） 为Landsat-8/OLI影像上试验场的位置。本文选取目标试验场主要考虑以下几点：（1）巴丹吉林沙漠随时间在亮度、空间均匀性、季节变化、长时期的稳定，且地表主要地物为沙子（ Yang等，2017 Yang等，2017 ），符合作为交叉定标试验场的条件（ Scott等，1996 Scott等，1996 ）；（2）马晓红等使用影像从阿拉伯沙漠选取均匀场地实现了对环境星的交叉定标（ 马晓红，2011 马晓红，2011 ），结果表明使用影像搜索沙漠地区均匀区域作为交叉辐射定标的试验场是可行的且有较高的精度；（3）在此区域能够获取一定数量连续的影像对数据，可实现时序交叉定标。 1;2,3; 2.3　数据 2 交叉辐射定标的精度依赖于参考传感器的精度，搭载在EOS-Terra/Aqua上的MODIS传感器，配备星上太阳辐射校正系统，使其绝对辐射和相对辐射定标精度分别在5%和1%（ Xiong 等，2006 Xiong 等，2006 ）；另外MODIS传感器重访周期短，覆盖范围广、拥有丰富的数据的特点使得其经常被用来作为参考传感器对其他卫星传感器进行交叉定标研究。HJ-1/CCD、CBERS-02/CCD、北京一号小卫星、NOAA-16、GOCI、ETM+等都以MODIS为参考传感器进行过交叉定标研究（ Hu 等，2001 Hu 等，2001 ； 李小英 等，2005 李小英 等，2005 ； Vermote 等，2006 Vermote 等，2006 ； 陈正超 等，2008 陈正超 等，2008 ； 马晓红，2011 马晓红，2011 ； 张玉环 等，2016 张玉环 等，2016 ）。考虑稳定的辐射性能、较高的辐射定标精度和丰富的数据，因此本文选择MODIS作为参考传感器对GF-4/PMS进行交叉辐射定标研究。 2.4　光谱匹配 2 在对GF-4/PMS进行交叉定前需要考虑与参考传感器（MODIS）光谱响应函数（SRF）的差异， 图3 图3 绘制了GF-4/PMS和MODIS可见光和近红外波段的光谱响应函数。其中，MODIS的光谱响应函数和波段大气层顶太阳辐照度数据来自NASA官网MCST团队 https：//mcst.gsfc.nasa.gov/calibration/parameters https://mcst.gsfc.nasa.gov/calibration/parameters ，GF-4/PMS的光谱响应函数和波段大气层顶太阳辐照度数据来自中国资源卫星中心（ http：//www.cresda.com/CN/Downloads/dbcs/ http://www.cresda.com/CN/Downloads/dbcs/ ），二者在可见光和近红外波段的波段信息如 表1 表1 所示。将影像中提取的试验场观测几何和地表反射光谱等信息输入到辐射传输模型MODTRAN中，计算出两个传感器模拟的表观辐亮度，得到两个传感器表观辐亮度的光谱匹配因子。为评估地表二向性的影响，也同时利用Ross-Li模型进行计算了光谱匹配因子（ Li 等 1992 Li 等 1992 ），相关输入参数来自于MCD43产品（ Schaaf 等，2011 Schaaf 等，2011 ），对应如下形式： <math display="block" id="M14"><mi>R</mi><mo stretchy="false">(</mo><msub><mrow><mi>θ</mi></mrow><mrow><mi>s</mi></mrow></msub><mo>,</mo><msub><mrow><mi>θ</mi></mrow><mrow><mi>v</mi></mrow></msub><mo>,</mo><mi>φ</mi><mo>,</mo><mi>λ</mi><mo stretchy="false">)</mo><mo>=</mo><msub><mrow><mi>f</mi></mrow><mrow><mi>i</mi><mi>s</mi><mi>o</mi></mrow></msub><mo stretchy="false">(</mo><mi>λ</mi><mo stretchy="false">)</mo><mo>+</mo><msub><mrow><mi>f</mi></mrow><mrow><mi>v</mi><mi>o</mi><mi>l</mi></mrow></msub><mo stretchy="false">(</mo><mi>λ</mi><mo stretchy="false">)</mo><msub><mrow><mi>K</mi></mrow><mrow><mi>v</mi><mi>o</mi><mi>l</mi></mrow></msub><mo stretchy="false">(</mo><msub><mrow><mi>θ</mi></mrow><mrow><mi>s</mi></mrow></msub><mo>,</mo><msub><mrow><mi>θ</mi></mrow><mrow><mi>v</mi></mrow></msub><mo>,</mo><mi>φ</mi><mo>,</mo><mi>λ</mi><mo stretchy="false">)</mo><mo>+</mo><msub><mrow><mi>f</mi></mrow><mrow><mi>g</mi><mi>e</mi><mi>o</mi></mrow></msub><mo stretchy="false">(</mo><mi>λ</mi><mo stretchy="false">)</mo><msub><mrow><mi>K</mi></mrow><mrow><mi>g</mi><mi>e</mi><mi>o</mi></mrow></msub><mo stretchy="false">(</mo><msub><mrow><mi>θ</mi></mrow><mrow><mi>s</mi></mrow></msub><mo>,</mo><msub><mrow><mi>θ</mi></mrow><mrow><mi>v</mi></mrow></msub><mo>,</mo><mi>φ</mi><mo>,</mo><mi>λ</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463137&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463136&type= 100.49933624 4.74133301 其中 <math id="M15"><msub><mrow><mi>f</mi></mrow><mrow><mi>i</mi><mi>s</mi><mi>o</mi></mrow></msub><mo stretchy="false">(</mo><mi>λ</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463197&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463195&type= 8.89000034 4.74133301 、 <math id="M16"><msub><mrow><mi>f</mi></mrow><mrow><mi>v</mi><mi>o</mi><mi>l</mi></mrow></msub><mo stretchy="false">(</mo><mi>λ</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463185&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463184&type= 9.05933285 4.74133301 与 <math id="M17"><msub><mrow><mi>f</mi></mrow><mrow><mi>g</mi><mi>e</mi><mi>o</mi></mrow></msub><mo stretchy="false">(</mo><mi>λ</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463202&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463200&type= 9.48266602 4.74133301 分别为各向同性散射、体散射与几何光学散射对应的权重系数， <math id="M18"><msub><mrow><mi>K</mi></mrow><mrow><mi>v</mi><mi>o</mi><mi>l</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463182&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463180&type= 4.65666676 3.72533321 与 <math id="M19"><msub><mrow><mi>K</mi></mrow><mrow><mi>g</mi><mi>e</mi><mi>o</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463193&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463192&type= 5.07999992 3.72533321 分别为体散射核与几何光学核。MCD43 A1产品提供了1~7个波段的 <math id="M20"><msub><mrow><mi>f</mi></mrow><mrow><mi>i</mi><mi>s</mi><mi>o</mi></mrow></msub><mo stretchy="false">(</mo><mi>λ</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463197&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463195&type= 8.89000034 4.74133301 、 <math id="M21"><msub><mrow><mi>f</mi></mrow><mrow><mi>v</mi><mi>o</mi><mi>l</mi></mrow></msub><mo stretchy="false">(</mo><mi>λ</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463185&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463184&type= 9.05933285 4.74133301 与 <math id="M22"><msub><mrow><mi>f</mi></mrow><mrow><mi>g</mi><mi>e</mi><mi>o</mi></mrow></msub><mo stretchy="false">(</mo><mi>λ</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463202&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463200&type= 9.48266602 4.74133301 ，按照MODIS 3、4、1、2波段与GF-4/PMS的2~5波段对应关系从MCD43 A1提取权重因子，然后按照MODTRAN运行要求输入对应参数，利用 公式（3） 公式（3） 计算光谱匹配因子。 2 GF-4/PMS 2 MODIS 波段号 波段范围（nm） 波段号 波段范围（nm） Band2 450 ~ 520 Band3 433 ~ 453 Band3 520 ~ 600 Band4 525 ~ 600 Band4 630 ~ 690 Band1 630 ~ 680 Band5 760 ~ 900 Band2 845 ~ 885 其中，光谱匹配因子计算使用的地面光谱数据来自（ Yang 等，2017 Yang 等，2017 ） 2012年7月在巴丹吉林沙漠地区实测，光谱曲线如 图4 图4 所示；根据（ Lacherade 等，2013 Lacherade 等，2013 ）的研究，在缺少可靠数据的情况下沙漠地区550nm气溶胶光学厚度可以采用默认值0.2，其余相关的大气参数使用MODTRAN模型的默认参数，由参数假设带来光谱配因子的计算误差分析见4.1节。为尽可能减少误差，这里气溶胶光学厚度数据采用了MAIAC算法生产的MCD19产品（ Lyapustin 等，2021 Lyapustin 等，2021 ）。 3　数据处理 1 数据：为获取尽量多的数据，选择MODIS和GF-4在试验区域同日过境的无云、清晰的影像对，在大气、成像角度和成像时刻约束下选取了2016年5月至2018年9月过境的13对GF-4和MODIS影像作为本次交叉定标的数据，各影像对的成像时间、时间差、试验场位置的成像几何和成像散射角度差异信息如 表2 表2 所示。 注： *：成像散射角是指太阳入射方向 <math id="M23"><mo stretchy="false">(</mo><msub><mrow><mi>θ</mi></mrow><mrow><mi>s</mi></mrow></msub><mo>,</mo><msub><mrow><mi>φ</mi></mrow><mrow><mi>s</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463217&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463215&type= 8.04333401 3.55599999 到观测方向 <math id="M24"><mo stretchy="false">(</mo><msub><mrow><mi>θ</mi></mrow><mrow><mi>v</mi></mrow></msub><mo>,</mo><msub><mrow><mi>φ</mi></mrow><mrow><mi>v</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463236&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463234&type= 8.29733276 3.55599999 的角度 <math id="M25"><mi>θ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463238&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463220&type= 1.35466671 2.03200006 （即 <math id="M26"><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mo stretchy="false">(</mo><mi>θ</mi><mo stretchy="false">)</mo><mo>=</mo><mo>-</mo><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mo stretchy="false">(</mo><msub><mrow><mi>θ</mi></mrow><mrow><mi>s</mi></mrow></msub><mo stretchy="false">)</mo><mi mathvariant="normal">*</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mo stretchy="false">(</mo><msub><mrow><mi>θ</mi></mrow><mrow><mi>v</mi></mrow></msub><mo stretchy="false">)</mo><mo>-</mo><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mo stretchy="false">(</mo><msub><mrow><mi>φ</mi></mrow><mrow><mi>s</mi></mrow></msub><mo>-</mo><msub><mrow><mi>φ</mi></mrow><mrow><mi>v</mi></mrow></msub><mo stretchy="false">)</mo><mi mathvariant="normal">*</mi><mo stretchy="false">|</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mo stretchy="false">(</mo><msub><mrow><mi>θ</mi></mrow><mrow><mi>v</mi></mrow></msub><mo stretchy="false">)</mo><mi mathvariant="normal">*</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mo stretchy="false">(</mo><msub><mrow><mi>θ</mi></mrow><mrow><mi>v</mi></mrow></msub><mo stretchy="false">)</mo><mo stretchy="false">|</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463221&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28463240&type= 70.18866730 3.55599999 ），GF-4和MODIS成像散射角差异反映了二者观测几何的相似程度。 2 成像日期 2 成像时间（北京时间） 2 成像时间差（分钟） 2 观测天顶角 2 太阳天顶角 2 相对方位角 2 成像散射角度差* GF-4 MODIS GF-4 MODIS GF-4 MODIS GF-4 MODIS 2016/6/2 12：16 13：05 46 43.36 57.62 21.31 27.93 -31.15 27.13 7.37 2016/9/2 11：30 11：55 25 43.00 31.87 39.55 36.62 -36.04 47.27 2.66 2016/9/25 11：00 12：00 60 43.02 22.53 49.82 43.42 -36.09 55.68 8.22 2016/10/4 9：51 11：50 119 43.02 32.06 62.43 46.99 -49.75 57.65 4.98 2017/6/11 11：24 12：30 66 43.63 31.33 28.57 19.01 -56.52 -134.93 11.39 2017/8/12 13：03 14：20 77 44.58 0.82 25.00 28.33 -1.72 136.62 9.33 2017/8/13 14：05 15：05 60 44.59 57.58 26.90 34.13 32.65 -35.99 8.70 2017/9/16 11：30 11：35 5 46.96 52.73 43.33 42.77 -28.52 48.93 16.49 2018/3/23 12：11 12：00 11 43.18 22.62 42.04 42.90 -19.01 50.26 19.59 2018/5/30 11：32 11：35 3 47.01 52.31 27.84 27.10 -47.74 29.53 2.71 2018/5/30 14：01 14：50 49 43.82 46.70 20.80 28.05 42.08 -26.47 6.75 2018/8/25 11：31 11：40 9 46.95 46.14 37.23 35.64 -37.33 41.37 2.07 2018/9/10 12：00 11：40 20 46.97 46.24 38.31 40.11 -20.87 47.77 16.19 影像对处理：在选取定标试验场之后GF-4影像处理步骤：（1）使用ENVI进行正射校正，裁剪出实验区域2397*1955大小（对应Landsat-8/OLI：3000*2000大小）；（2）使用Landsat-8/OLI影像为基准对GF-4/PMS影像进行几何配准，配准精度优于1个像元；（3）提取GF-4/PMS试验场的观测几何信息（观测天顶角、观测方位角、太阳天顶角和太阳方位角）及各波段的DN值。MODIS影像的处理步骤：（1）使用ENVI进行几何校正；（2）提取MODIS和GF-4在试验场位置的观测几何以及各波段的DN值；（3）根据交叉定标区域位置提取MCD19A2中气溶胶光学厚度、MCD43A1中BRDF数值，其中MCD43A1数值只采用MCD43A2中质量标记为0或者1的数值。将提取的观测及大气参数等信息输入MODTRAN辐射传输模型，计算出GF-4/PMS和MODIS各波段的模拟的表观辐亮度（TOA），求得参考传感器与待标定传感器的光谱匹配因子。利用光谱匹配因子和MODIS观测得到的表观辐亮度计算出GF-4/PMS模拟的入瞳辐亮度；将GF-4/PMS的DN值与模拟出的入瞳辐亮度进行线性拟合，得到交叉辐射定标系数。为进行BRDF效应的对比分析，分别基于实测光谱的朗伯体假设情况和基于MCD43产品的BRDF模型情况两种方式计算光谱匹配因子。 4　结 果 1 4.1　实验结果与分析 2 本文得到的朗伯体假设和BRDF模型计算的GF-4/PMS交叉辐射定标系数如 表3 表3 所示，中国资源卫星应用中心公布的2016至2018年场地定标系数见 表4 表4 。为了更直观的分析交叉定标结果，本文绘制了朗伯体假设和BRDF模型两种情况下计算的交叉定标系数和中国资源卫星应用对比散点图（ 图5 图5 ），其中虚线代表交叉定标系数随时间的变化趋势。结果表明：（1）由于限定了过境时间差和成像角度差，朗伯体假设和BRDF模型假设条件下计算的交叉定标系数具有较高一致性，绿、蓝、红和近红外波段二者的相对差异分别为3.86±2.36%、3.36±1.87%、3.88±2.21%和3.89±2.00%；（2）2016至2018年间GF-4/PMS4各个波段的定标系数整体呈现缓慢增加趋势，表明传感器辐射性能出现下降趋势，朗伯体假设和BRDF模型假设条件下4个波段的年衰减率分别为1.00%、0.69%、0.43%、0.28%和0.90%、0.39%、0.18%、0.06%；（3）每年8-10月附近交叉定标结果平均值与中国资源卫星应用中心场地定标公布的定标系数的相对误差，2017年偏差明显（超过9%），这一方面可能与遥感器的辐射稳定性有关，另一方面与单次定标出现的随机误差有关（ Chen 等， 2014 Chen 等， 2014 ）；（4）相比朗伯体假定条件下计算的交叉定标系数，BRDF模型假定下计算结果与场地定标系数偏差略高，这可能与后者采用朗伯体假定计算的定标方法有关，该问题需要后续更多场地实测BRDF支持情况下进一步验证。 2 影像日期 4 朗伯体假定情况的交叉定标系数 4 BRDF假定情况下的交叉定标系数 Band1 Band2 Band3 Band4 Band1 Band2 Band3 Band4 2016/6/2 0.1794 0.2053 0.1780 0.1335 0.1796 0.1937 0.1722 0.1260 2016/9/2 0.1778 0.1926 0.1596 0.1172 0.1890 0.1975 0.1685 0.1211 2016/9/25 0.1843 0.1953 0.1597 0.1186 0.1973 0.2038 0.1730 0.1266 2016/10/4 0.1895 0.1967 0.1607 0.1193 0.2002 0.2025 0.1710 0.1249 2017/6/11 0.1810 0.2010 0.1649 0.1227 0.1957 0.2110 0.1782 0.1301 2017/8/12 0.1847 0.1998 0.1682 0.1221 0.1978 0.2055 0.1785 0.1261 2017/8/13 0.1929 0.2209 0.1881 0.1376 0.1931 0.2089 0.1835 0.1289 2017/9/16 0.1964 0.2050 0.1644 0.1233 0.2013 0.2002 0.1655 0.1204 2018/3/23 0.1944 0.1984 0.1590 0.1201 0.2046 0.2001 0.1674 0.1242 2018/5/30 0.2097 0.2172 0.1826 0.1325 0.2139 0.2085 0.1795 0.1269 2018/5/30 0.1995 0.2231 0.1822 0.1401 0.2031 0.2130 0.1781 0.1333 2018/8/25 0.1994 0.2054 0.1701 0.1245 0.2069 0.2023 0.1720 0.1232 2018/9/10 0.1997 0.2099 0.1727 0.1242 0.2073 0.2087 0.1768 0.1247 注： *：中国资源卫星应用中心场地辐射定标工作一般在每年9月份前后。 定标年份 * Band1 Band2 Band3 Band4 2016 0.1784 0.1878 0.1515 0.108 2017 0.1749 0.1903 0.1532 0.1073 2018 0.2053 0.2209 0.1664 0.115 1;2;3;4; 从 图5 图5 中进一步发现2018年5月30日两幅GF-4分别与MODIS Terra和Aqua进行交叉定标，二者在4个波段差异小于6%，说明本文方法具有较高精度。为进一步检验该结果，我们放宽角度和成像时间差的差异，发现同一天GF-4对MODIS不同时刻交叉定标的结果保持了较好的一致性，相对误差均小于6%（如 表5 表5 ）。 2 日期 2 成像时刻 4 朗伯体假定情况的交叉定标系数 4 BRDF假定情况下的交叉定标系数 GF-4 MODIS Band1 Band2 Band3 Band4 Band1 Band2 Band3 Band4 2017/6/11 11：24 12：30 0.1810 0.2010 0.1649 0.1227 0.1957 0.2110 0.1782 0.1301 2017/6/11 11：24 14：10 0.1823 0.1988 0.1648 0.1217 0.1975 0.2087 0.1779 0.1286 2017/8/12 13：03 12：20 0.1734 0.1936 0.1608 0.1195 0.1870 0.2041 0.1763 0.1282 2017/8/12 13：03 14：20 0.1847 0.1998 0.1682 0.1221 0.1978 0.2055 0.1785 0.1261 2017/8/13 14：05 11：45 0.1943 0.2139 0.1809 0.1333 0.1985 0.2075 0.1809 0.1283 2017/8/13 14：05 15：05 0.1929 0.2209 0.1881 0.1376 0.1931 0.2089 0.1835 0.1289 2018/3/23 12：11 12：00 0.1944 0.1984 0.1590 0.1201 0.2046 0.2001 0.1674 0.1242 2018/3/23 12：11 13：40 0.1854 0.1979 0.1626 0.1213 0.1952 0.2012 0.1726 0.1271 2018/5/30 11：32 11：35 0.2097 0.2172 0.1826 0.1325 0.2139 0.2085 0.1795 0.1269 2018/5/30 11：32 14：50 0.2041 0.2150 0.1807 0.1302 0.2088 0.2072 0.1784 0.1251 2018/5/30 14：01 11：35 0.2050 0.2254 0.1840 0.1425 0.2080 0.2143 0.1792 0.1352 2018/5/30 14：01 14：50 0.1995 0.2231 0.1822 0.1401 0.2031 0.2130 0.1781 0.1333 2018/8/25 11：31 11：40 0.1994 0.2054 0.1701 0.1245 0.2069 0.2023 0.1720 0.1232 2018/8/25 11：31 15：00 0.1978 0.2094 0.1753 0.1275 0.2023 0.2014 0.1728 0.1227 2018/9/10 12：00 11：40 0.1997 0.2099 0.1727 0.1242 0.2073 0.2087 0.1768 0.1247 2018/9/10 12：00 15：00 0.2024 0.2176 0.1793 0.1290 0.2073 0.2111 0.1784 0.1254 4.2　不确定性分析 2 通过对GF-4/PMS与Terra（Aqua）/MODIS的交叉定标过程分析，影响定标精度的因素主要有以下几方面：（a）MODIS自身定标的不确定性；（b）影像空间配准带来的影响；（c）大气参数误差的影响；（d）地表反射光谱的误差；（e）大气辐射传输模型本身计算的不确定性；（f）地表二向性的不确定性；（f）其它因素引起的误差。 （1）MODIS定标系数的不确定性为±3%（ Chang等，2017 Chang等，2017 ）。 巩慧等（2010） 巩慧等（2010） 用2007年二连浩特对MODIS的同步实验数据，采用反射率基法对MODIS进行了辐射定标和真实性检验研究，证明了MODIS具有较高的定标精度。 （2）为了分析空间像元匹配误差对交叉定标结果的影响，本文将所有选用影像的试验场扩大一倍（1050m*1050m）的平均DN值与试验场的平均DN值做比较最大差异分别为：蓝波段0.28%、绿波段0.24%、红波段0.4%和近红外波段0.19%。因此，即使几何配准误差在1个像元甚至更多像元情况下，DN值的差异也很小。 （3）大气条件对光谱匹配因子的计算的影响主要为：气溶胶类型、气溶胶光学厚度和水汽含量。由于本文交叉定标缺乏同步测量数据，我们使用MODTRAN中所提供的标准大气成分组成，550nm气溶胶光学厚度设置为0.2。在其它变量不变的情况下，分别改变气溶胶类型（沙漠型、乡村型和城市型）、550nm气溶胶光学厚度（0.05至0.4）和水汽含量（0.4至4.0）的值计算光谱匹配因子的差异，计算结果如 表5 表5 所示； 图6 图6 列出了2016年8月25日不同条件对光谱匹配因子的影响；从 表5 表5 和 图6（b） 图6（b） 中均可看出水汽含量对近红外波段光谱匹配因子的影响较大。 1,2;3,4; （4）本文所使用的地面光谱为巴丹吉林沙漠中心地带的光谱曲线，并非试验场实测的光谱，为了分析误差我们使用敦煌定标场的地面光谱来计算光谱匹配因子，分析在没有实测光谱的情况下，使用同类地物光谱替代对交叉定标结果所带来的影响； 4个波段计算的光谱匹配因子差异分别为：5.04%、4.02%、0.02%和2.94%，说明使用同类地物光谱计算光谱匹配因子进行交叉定标在可接受的范围之内；本文所使用的验证地物的光谱曲线如 图7 图7 所示。 （5）前面对比了采用朗伯体假定和Ross-Li BRDF模型假定得到交叉定标差异，蓝、绿、红、近红外波段的相对差异为3.86±2.36%、3.36±1.87%、3.88±2.21%和3.89±2.00%。 综合上述分析，通过对GF-4/PMS与MODIS交叉定标误差来源分析，得到此次交叉定标总的不确定度如 表5 表5 所示，总体误差在7.4%以内，高于垂直观测卫星的交叉定标误差估算结果（ Chen 等，2014 Chen 等，2014 ）。 误差来源 蓝 绿 红 近红外 MODIS定标系数（%） 3.00 3.00 3.00 3.00 MODTRAN精度（%） 2.00 2.00 2.00 2.00 像元几何配准误差（%） 0.28 0.24 0.40 0.19 气溶胶类型（%） 0.28 0.13 0.07 0.14 气溶胶光学厚度（%） 0.69 0.60 0.30 0.80 水汽含量（%） 0.12 0.81 0.07 2.93 地表光谱替代（%） 5.04 4.02 0.02 2.94 地表二向性误差（%） 3.86 3.36 3.88 3.89 其它不确定因素（%） 1.00 1.00 1.00 1.00 总的不确定度（%） 7.41 6.52 5.41 6.86 5　结 论 1 受到太空环境、遥感器自身光电性能下降等因素的影响，遥感器辐射性能不可避免发生变化，因此，一年一次的场地定标频次难以监测遥感器性能的变化情况，无法满足卫星传感器数据长期定量化的需求。本文以辐射定标精度较高的MODIS为参考传感器，对2016年至2018年GF-4/PMS可见光近红外波段进行时序定标，为GF-4/PMS长期辐射性能评估和监测提供一种技术手段。结果表明： （1）GF-4在中国区域具有较大观测角度，通过约束与MODIS成像角度差异和成像时间差异，能较好减少因为大气变化和BRDF差异带来的影响，忽略和考虑沙地BRDF效应带来的差异约4%，两种情况下得到的交叉定标系数具有较高的一致性，定标不确定度小于7.4%。 （2）与中国资源卫星中心的场地定标系数相比，2016和2018年交叉定标结果与场地定标系数较为接近，但2017年偏离较大，这可能与场地定标朗伯体假定、交叉定标自身不确定度有关。 （3）高频次的交叉定标能够有效检查单次定标可能存在的误差，通过交叉定标和场地定标结果相互检查能够剔除存在的粗差点，提高结果的可靠性。 （4）通过时序的交叉定标，我们发现GF-4/PMS各个波段的辐射性能从2016年至2018年呈现缓慢下降趋势，年衰减率小于1%。 由于GF-4卫星静止轨道特征，与极轨卫星有更多交叉成像机会，在约束成像角度和时间情况下，继续监测GF-4的辐射性能变化情况。 参考文献 （References） Chen Z C , Liu X , Li J X , Lou W F , Zang L , Zhang H , Zhang X and Zhang B . 2008 . The Cross Calibration of Beijing-1 Microsatellite Multispectral Sensors. Journal of Astronautics, (02): 637 - 643 The Cross Calibration of Beijing-1 Microsatellite Multispectral Sensors. Journal of Astronautics, (02): 陈正超 , 刘翔 , 李俊生 , 罗文斐 , 张靓 , 张浩 , 张霞 , 张兵 . 2008 . 北京一号小卫星多光谱遥感器交叉定标 . 宇航学报 , ( 02 ): 637 - 643 [ DOI: CNKI:SUN:YHXB.0.2008-02-048 http://dx.doi.org/CNKI:SUN:YHXB.0.2008-02-048 ] 北京一号小卫星多光谱遥感器交叉定标 Gao H L , Gu X F , Yu T , Li X Y , Gong H and Li J G . 2010 . The Research Overview on Visible and Near-infrared Channels Radiometric Calibration of Space-borne Optical Remote Sensors. Remote Sensing Information, (04): 117 - 128 The Research Overview on Visible and Near-infrared Channels Radiometric Calibration of Space-borne Optical Remote Sensors. Remote Sensing Information, (04): 高海亮 , 顾行发 , 余涛 , 李小英 , 巩慧 , 李家国 . 2010 . 星载光学遥感器可见近红外通道辐射定标研究进 . 遥感信息 , ( 04 ): 117 - 128 [ DOI: CNKI:SUN:YGXX.0.2010-04-029 http://dx.doi.org/CNKI:SUN:YGXX.0.2010-04-029 ] 星载光学遥感器可见近红外通道辐射定标研究进 Gong H , Tian G L , Yu T , Zhang Y X , Gu X F and Gao H L . 2010 . Irradiance-based calibration and validation of MODIS visible and near-infrared channels . Journal of Remote Sensing , 14 ( 02 ): 207 - 218 Irradiance-based calibration and validation of MODIS visible and near-infrared channels 巩慧 , 田国良 , 余涛 , 张玉香 , 顾行发 , 高海亮 . 2010 . MODIS辐照度法定标试验研究 . 遥感学报 , 14 ( 02 ): 207 - 218 [ DOI: CNKI:SUN:YGXB.0.2010-02-003 http://dx.doi.org/CNKI:SUN:YGXB.0.2010-02-003 ] MODIS辐照度法定标试验研究 Li X Y , Gu X F , Min X J , Yu T , Fu Q Y , Zhang Y and Li X W . 2005 . Radiometric cross-calibration of the CBERS-02 CCD camera with the TERRA MODIS. Science in China Ser, E Information Sciences, (S1): 41 - 58 Radiometric cross-calibration of the CBERS-02 CCD camera with the TERRA MODIS. Science in China Ser, E Information Sciences, (S1): 李小英 , 顾行发 , 闵祥军 , 余涛 , 傅俏燕 , 张勇 , 李小文 . 2005 . 利用MODIS对CBERS-02卫星CCD相机进行辐射交叉定标 . 中国科学E辑 : 信息科学 , ( S1 ): 41 - 58 [ DOI: CNKI:SUN:JEXK.0.2005-S1-004 http://dx.doi.org/CNKI:SUN:JEXK.0.2005-S1-004 ] 利用MODIS对CBERS-02卫星CCD相机进行辐射交叉定标 Lv W B . 2014 . Wenbo . Research on the method and system of cross-calibration for domestic moderate and high-resolution remote sensing data in reflection bands. Chengdu : Southwest Jiaotong University (吕文博 . 2014. 国产中—高分辨率遥感数据反射波段交叉辐射定标方法与体系研究, 西南交通大学) Wenbo Nie J , Deng L , Hao X L , Liu M and He Y . 2018 . Application of GF-4 satellite in drought remote sensing monitoring: A case study of Southeastern Inner Mongolia . Journal of Remote Sensing , 22 ( 03 ): 400 - 407 Application of GF-4 satellite in drought remote sensing monitoring: A case study of Southeastern Inner Mongolia 聂娟 , 邓磊 , 郝向磊 , 刘明 , 贺英 . 2018 . 高分四号卫星在干旱遥感监测中的应用 . 遥感学报 , 22 ( 03 ): 400 - 407 [ DOI: CNKI:SUN:YGXB.0.2018-03-004 http://dx.doi.org/CNKI:SUN:YGXB.0.2018-03-004 ] 高分四号卫星在干旱遥感监测中的应用 Ma X H . 2011 . Study on the Cross Calibration and Validation of CCD Cameras on HJ-1Satellites . Henan Polytechnic University , 2011 Study on the Cross Calibration and Validation of CCD Cameras on HJ-1Satellites 马晓红 . 2011 . HJ-1星CCD相机交叉定标与真实性检验研究, 河南理工大学 HJ-1星CCD相机交叉定标与真实性检验研究, 河南理工大学 Wu W . 2019 . Application of GF-4 satellite in outburst flood disaster monitoring and assessment . Spacecraft Engineering , 28 ( 2 ): 134 - 140 Application of GF-4 satellite in outburst flood disaster monitoring and assessment 吴玮 . 2019 . 高分四号卫星在溃决型洪水灾害监测评估中的应用 . 航天器工程 , 28 ( 2 ): 134 - 140 [ DOI: CNKI:SUN:HTGC.0.2019-02-020 http://dx.doi.org/CNKI:SUN:HTGC.0.2019-02-020 ] 高分四号卫星在溃决型洪水灾害监测评估中的应用 Yang Z D , Gu S Y , Qiu H , Huang Q , Fan T X . 2004 . CBERS-1 's CCD Image Quality Evaluating and Cross Calibrating Study. Journal of Remote Sensing, (02): 113 - 120 CBERS-1 's CCD Image Quality Evaluating and Cross Calibrating Study. Journal of Remote Sensing, (02): 杨忠东 , 谷松岩 , 邱红 , 黄签 , 范天锡 . 2004 . 中巴地球资源一号卫星CCD图像质量评价和交叉定标研究 . 遥感学报 , ( 02 ): 113 - 120 [ DOI: CNKI:SUN:YGXB.0.2004-02-003 http://dx.doi.org/CNKI:SUN:YGXB.0.2004-02-003 ] 中巴地球资源一号卫星CCD图像质量评价和交叉定标研究 Zhang L . 2018 . Application of GF-4satellite in changes monitoring for water body extent of Hongze Lake in flood season . Spacecraft Engineering , 27 ( 02 ): 129 - 134 Application of GF-4satellite in changes monitoring for water body extent of Hongze Lake in flood season 张磊 . 2018 . 高分四号卫星在洪泽湖汛期水体范围变化监测中的应用 . 航天器工程 , 27 ( 02 ): 129 - 134 [ DOI: CNKI:SUN:HTGC.0.2018-02-020 http://dx.doi.org/CNKI:SUN:HTGC.0.2018-02-020 ] 高分四号卫星在洪泽湖汛期水体范围变化监测中的应用 Zhang Y H , Mao H Q , Wang Z T , Guo W , Li Q , Li Z Q , Chen X F and Chen H . 2016 . Cross Calibration Between GOCI and MODIS Based on Spectrum and Geometry Matching . Journal of Atmospheric and Environmental Optics , 11 ( 6 ): 412 - 422 Cross Calibration Between GOCI and MODIS Based on Spectrum and Geometry Matching 张玉环 , 毛慧琴 , 王中挺 , 郭伟 , 厉青 , 李正强 , 陈兴峰 , 陈辉 . 2016 . 基于光谱与几何匹配的GOCI与MODIS交叉辐射定标 . 大气与环境光学学报 , 11 ( 06 ): 412 - 422 [ DOI: CNKI:SUN:GDJY.0.2016-06-003 http://dx.doi.org/CNKI:SUN:GDJY.0.2016-06-003 ] 基于光谱与几何匹配的GOCI与MODIS交叉辐射定标 Zhao W N , Hu X Q , Fang W , Wang Y P and Xu N . 2015 . Development and applications of intercalibration for satellite optical instruments . Optics and Precision Engineering , 23 ( 7 ): 1921 - 1931 Development and applications of intercalibration for satellite optical instruments 赵维宁 , 胡秀清 , 方伟 , 王玉鹏 , 徐娜 . 2015 . 卫星光学仪器辐射交互定标方法的应用和发展 . 光学精密工程 , 23 ( 07 ): 1921 - 1931 [ DOI: CNKI:SUN:GXJM.0.2015-07-016 http://dx.doi.org/CNKI:SUN:GXJM.0.2015-07-016 ] 卫星光学仪器辐射交互定标方法的应用和发展 Cao C Y , Xu H , Sullivan J , McMillin L , Ciren P and Hou YT . 2005 . Intersatellite radiance biases for the High-Resolution Infrared Radiation Sounders (HIRS) on board NOAA-15,-16, and-17 from simultaneous nadir observations . Journal of Atmospheric and Oceanic Technology , 22 ( 4 ): 381 - 395 [ DOI: 10.1175/JTECH1713.1 http://dx.doi.org/10.1175/JTECH1713.1 ] Intersatellite radiance biases for the High-Resolution Infrared Radiation Sounders (HIRS) on board NOAA-15,-16, and-17 from simultaneous nadir observations Chen YP , Sun KM , Li DR , Bai T and Huang CQ . 2017 . Radiometric cross-calibration of gf-4 pms sensor based on assimilation of landsat-8 oli images . IEEE Transactions on Geoscience and Remote sensing , 55 ( 4 ): 2288 - 2298 [ DOI: 10.3390/rs9080811 http://dx.doi.org/10.3390/rs9080811 ] Radiometric cross-calibration of gf-4 pms sensor based on assimilation of landsat-8 oli images Chen ZC , Zhang B , Zhang H and Zhang W J . 2014 . Vicarious Calibration of Beijing-1 Multispectral Imagers . Remote Sensing , 6 : 1432 - 1450 [ DOI: 10.3390/rs6021432 http://dx.doi.org/10.3390/rs6021432 ] Vicarious Calibration of Beijing-1 Multispectral Imagers Chang T J , Xiong X X , Angal A , Wu A S and Geng X . 2017 . Aqua and Terra MODIS RSB calibration comparison using BRDF modeled reflectance . IEEE Transactions on Geoscience and Remote sensing , 55 ( 4 ): 2288 - 2298 [ DOI: 10.1109/TGRS.2016.2641258 http://dx.doi.org/10.1109/TGRS.2016.2641258 ] Aqua and Terra MODIS RSB calibration comparison using BRDF modeled reflectance Hu C M , Muller-Karger F E , Andrefouet S and Carder K L . 2001 . Atmospheric correction and cross-calibration of LANDSAT-7/ETM+ imagery over aquatic environments: A multiplatform approach using SeaWiFS/MODIS . Remote Sensing of Environment , 78 ( 1-2 ): 99 - 107 [ DOI: 10.1016/S0034-4257(01)00252-8 http://dx.doi.org/10.1016/S0034-4257(01)00252-8 ] Atmospheric correction and cross-calibration of LANDSAT-7/ETM+ imagery over aquatic environments: A multiplatform approach using SeaWiFS/MODIS Lacherade S , Fougnie B , Henry P and Gamet P . 2013 . Cross calibration over desert sites: Description, methodology, and operational implementation . IEEE Transactions on Geoscience and Remote sensing , 51 ( 3 ): 1098 - 1113 [ DOI: 10.1109/TGRS.2012.2227061 http://dx.doi.org/10.1109/TGRS.2012.2227061 ] Cross calibration over desert sites: Description, methodology, and operational implementation Li X and Strahler A H . 1992 . Geometric-optical bidirectional reflectance modeling of the discrete crown vegetation canopy: Effect of crown shape and mutual shadowing . IEEE Trans. Geosci. Remote Sens. , 30 : 276 - 292 . Geometric-optical bidirectional reflectance modeling of the discrete crown vegetation canopy: Effect of crown shape and mutual shadowing Lyapustin A , Wang Y , Go S , et al . 2021 . Atmospheric Correction of DSCOVR EPIC: Version 2 MAIAC Algorithm. Frontiers in Remote Sensing, 2: [ 10.3389/frsen.2021.748362 http://dx.doi.org/10.3389/frsen.2021.748362 ] Atmospheric Correction of DSCOVR EPIC: Version 2 MAIAC Algorithm. Frontiers in Remote Sensing, 2: [ Rao CRN , Zhang N and Sullivan JT . 2001 . Inter-calibration of meteorological satellite sensors in the visible and near-infrared . Advances in Space Research , 28 ( 1 ): 3 - 10 [ DOI: 10.1016/S0273-1177(01)00262-9 http://dx.doi.org/10.1016/S0273-1177(01)00262-9 ] Inter-calibration of meteorological satellite sensors in the visible and near-infrared Schaaf C L B , Liu J , Gao F and Strahler A H . 2011 . MODIS Albedo and Reflectance Anisotropy Products from Aqua and Terra, In Land Remote Sensing and Global Environmental Change: NASA's Earth Observing System and the Science of ASTER and MODIS, Remote Sensing and Digital Image Processing Series, Vol. 11 , B. Ramachandran, C. Justice, M. Abrams, Eds, Springer-Cerlag, 873 pp. MODIS Albedo and Reflectance Anisotropy Products from Aqua and Terra, In Land Remote Sensing and Global Environmental Change: NASA's Earth Observing System and the Science of ASTER and MODIS, Remote Sensing and Digital Image Processing Series, Vol. Scott K P , Thome KJ and Brownlee MR . 1996 . Evaluation of the Railroad Valley playa for use in vicarious calibration. Optics & Photonics, Conference on Multispectral Imaging for Terrestrial Applications, Denver . Evaluation of the Railroad Valley playa for use in vicarious calibration. Optics & Photonics, Conference on Multispectral Imaging for Terrestrial Applications, Denver Teillet P M , Fedosejevs G , Thome K J and Barker J L . 2007 . Impacts of spectral band difference effects on radiometric cross-calibration between satellite sensors in the solar-reflective spectral domain . Remote Sensing of Environment , 110 ( 3 ): 393 - 409 [ DOI: 10.1016/j.rse.2007.03.003 http://dx.doi.org/10.1016/j.rse.2007.03.003 ] Impacts of spectral band difference effects on radiometric cross-calibration between satellite sensors in the solar-reflective spectral domain Teillet P M , Slater P N , Ding Y , Santer RP , Jackson RD and Moran MS . 1990 . Three methods for the absolute calibration of the NOAA AVHRR sensors in-flight . Remote Sensing of Environment , 31 ( 2 ): 105 - 120 [ DOI: 10.1016/0034-4257(90)90060-Y http://dx.doi.org/10.1016/0034-4257(90)90060-Y ] Three methods for the absolute calibration of the NOAA AVHRR sensors in-flight Vermote E F , Saleous N Z . 2006 . Calibration of NOAA16 AVHRR over a desert site using MODIS data . Remote Sensing of Environment , 105 ( 3 ): 214 - 220 [ DOI: 10.1016/j.rse.2006.06.015 http://dx.doi.org/10.1016/j.rse.2006.06.015 ] Calibration of NOAA16 AVHRR over a desert site using MODIS data Xiong X X , Barnes W . 2006 . An overview of MODIS radiometric calibration and characterization . Advances in Atmospheric Sciences , 23 ( 1 ): 69 - 79 [ DOI: 10.1007/s00376-006-0008-3 http://dx.doi.org/10.1007/s00376-006-0008-3 ] An overview of MODIS radiometric calibration and characterization Yang A X , Zhong B , Lv W B , Wu S L and Liu Q H . 2015 . Cross-calibration of GF-1/WFV over a desert site using Landsat-8/OLI imagery and ZY-3/TLC data . Remote Sensing , 7 ( 8 ): 10763 - 10787 [ DOI: 10.3390/rs70810763 http://dx.doi.org/10.3390/rs70810763 ] Cross-calibration of GF-1/WFV over a desert site using Landsat-8/OLI imagery and ZY-3/TLC data Yang A X , Zhong B , Wu S L and Liu Q H . 2017 . Radiometric cross-calibration of GF-4 in multispectral bands . Remote Sensing , 9 ( 3 ): 232 [ DOI: 10.3390/rs9030232 http://dx.doi.org/10.3390/rs9030232 ] Radiometric cross-calibration of GF-4 in multispectral bands Zhong B , Zhang Y H , Du T T , Yang A X , Lv W B and Liu Q H . 2014 . Cross-Calibration of HJ-1/CCD over a desert site using landsat ETM plus imagery and ASTER GDEM product . IEEE Transactions on Geoscience and Remote sensing , 52 ( 11 ): 7247 - 7263 [ DOI: 10.1109/TGRS.2014.2310233 http://dx.doi.org/10.1109/TGRS.2014.2310233 ] Cross-Calibration of HJ-1/CCD over a desert site using landsat ETM plus imagery and ASTER GDEM product Zhang H , Zhang B , Chen Z C , Huang Z H . 2018 . Vicarious radiometric calibration of the hyperspectral imaging microsatellite SPARK-01 and -02 over Dunhuang, China , Remote Sensing , 10 , 120 [ DOI: 10.3390/rs10010120 http://dx.doi.org/10.3390/rs10010120 ] Vicarious radiometric calibration of the hyperspectral imaging microsatellite SPARK-01 and -02 over Dunhuang, China 致谢:本文所使用的巴丹吉林沙漠的地表实测光谱曲线由中国科学院空天信息研究院的杨爱霞博士提供，在此表示衷心的感谢；感谢中国资源卫星中心提供了大量的GF-4数据；感谢匿名审稿人细致的评审。 10.11834/jrs.20221774 TP701 A 张浩，刘涛，闫东川，阎跃观，崔珍珍.XXXX.基于MODIS的GF-4/PMS遥感器交叉定标-以巴丹吉林沙漠为参考目标.遥感学报，XX（XX）： 1-14 张浩，刘涛，闫东川，阎跃观，崔珍珍.XXXX.基于MODIS的GF-4/PMS遥感器交叉定标-以巴丹吉林沙漠为参考目标.遥感学报，XX（XX）： 1-14 ZHANG Hao，LIU Tao，YAN Dongchuan，YAN Yueguan，CUI Zhenzhen. XXXX. Cross calibration of GF-4/PMS Based on MODIS over Badain Jaran Desert. National Remote Sensing Bulletin， XX（XX）：1-14 ZHANG Hao，LIU Tao，YAN Dongchuan，YAN Yueguan，CUI Zhenzhen. XXXX. Cross calibration of GF-4/PMS Based on MODIS over Badain Jaran Desert. National Remote Sensing Bulletin， XX（XX）：1-14 海南省重大科技计划项目—基于天基大数据的海南省生态资源监管关键技术与应用（编号：ZDKJ2019006）资助； 国家自然科学基金（编号：41771397） Supported by Supported by Hainan Provincial Department of Science and Technology under Grant (No. ZDKJ2019006), National Natural Science Foundation of China (No. 41771397) 2022-05-26 遥感学报