基于MODIS影像的表观反射率模型构建方法、系统及定标方法技术方案

本发明专利技术涉及基于MODIS影像的表观反射率模型构建方法、系统及定标方法，通过选取有效定标场地影像，得到对应的图像表观反射率及相关角度信息，进一步利用核驱动模型，得到表观反射率校正系数，从而建立表观反射率模型，用于获取表观反射率，在此基础上实现不同传感器卫星影像对的交叉定标，本发明专利技术方法可适用于多光谱卫星的在轨定标，尤其是大角度影像对之间的交叉定标，该方法具有定标精度高，可用于历史数据再定标等优点，本发明专利技术不用开展星地同步实验，取消了场地必须具有朗伯性的假设，不仅可用于小角度影像的交叉定标，而且可用于大角度下影像对的交叉定标，极大提高了定标的频率和精度。

【技术实现步骤摘要】
基于MODIS影像的表观反射率模型构建方法、系统及定标方法
本专利技术涉及一种基于MODIS影像的表观反射率模型构建方法、系统及定标方法，属于定标检校

技术介绍
高精度、高频次的在轨辐射定标是实现定量遥感产品反演的前提和基础。在轨辐射定标是在卫星运行期间，利用地面或星上特定目标作为参考辐射源，建立辐射亮度与图像灰度值的转换关系，得到传感器的定标系数。根据参考辐射源的不同，在轨辐射定标方法又分为星上定标器定标法、太阳定标法、月球定标法、场地定标法、场景定标法、交叉定标法等。其中，星上定标器定标法利用星上内置积分球作为定标源，实现在轨辐射定标。太阳定标法和月球定标法分别选择太阳和月球作为参考辐射源，通过调整卫星姿态，利用星上装置获取太阳和月球影像，实现在轨辐射定标。这三种方法都要求在卫星上配置相应的定标装置，在定标时需要对卫星姿态进行调整，对卫星载荷的控制要求高，可定标的频次有限，目前国产陆地卫星大部分都不具备这三种的定标能力。场地定标法依赖地面星地同步实验数据，利用定标场实测地表和大气参数，实现在轨辐射定标。该方法的定标精度和频率受限于星地同步实验的次数及过境时刻的天气情况，每次定标都需要耗费大量的人力物力，可开展的定标次数有限。场景定标法利用特定的自然场景目标作为参考辐射源，利用模型辐亮度和图像的灰度值，实现在轨辐射定标。根据自然场景的不同，又分为沙漠场景法、云场景法、海洋场景法、极地场景法等。场景定标法的优点在于定标频率高、成本低、可实现历史影像数据的定标。但不足之处在于定标精度受限于定标场景的模型，该方法假设不同地区不同时相的场景具有相同的模型，而实际上不同地区不同时相的场景影像，由于观测几何、大气条件、地表变化的差异，这一假设存在较大的不确定性，导致场景定标法精度较低。交叉定标法选择参考卫星作为参考辐射源，通过建立目标传感器与参考传感器图像之间的转换关系，实现传感器的绝对辐射定标。该方法考虑了不同卫星成像条件、光谱分辨率等之间的差异，其定标精度和选择的参考影像对直接相关。通常情况下，选择相同观测角度、相同观测时间下的影像对具有较高的定标精度。现有的《基于冷云目标的气象卫星太阳反射波段辐射定标方法》(CN105092055B)公开了一种基于冷云目标的气象卫星太阳反射波段辐射定标方法。该方法通过冷云目标物的提取，冷云反射率归一化处理；辐射定标响应的在轨状态监控与日衰减模型建立等处理，利用辐射参考卫星，实现气象卫星太阳反射波段重新定标。该方法属于交叉定标的一种。本专利技术实现了卫星仪器观测的连续性和一致性；不受天气条件影响，省时省力，可以方便快捷、准实时得到卫星响应变化。然而，该方法并未考虑大观测角度下(观测角度大于30度)地表方向反射率的影响，不适用于大观测角度下影像对的交叉定标。现有的交叉定标研究中，大多数研究都假设地表为朗伯体，忽略了不同卫星观测几何的影响，这一假设在卫星观测角度偏小时合适，但当卫星具有较大观测角度(大于30度)时，忽略场地反射率方向性的影响，将会带来较大的误差。现有的方法《考虑地物BRDF特性改进后的CBERS202卫星CCD相机的辐射定标系数》(遥感学报2006，第10卷第5期)虽然考虑了场地方向反射率的影响，但该方法采用的是地面测量的多角度反射率数据，只考虑了地表反射率方向性的影响，而未考虑大气层顶表观反射率方向性的影响。现有的方法《基于深度学习的在轨辐射定标方法研究》(航天返回与遥感，2017，第38卷第2期)采用时间序列MODIS作为参考卫星利用定标场地的大量历史卫星影像、历史大气数据和历史光谱数据，通过对这些数据的学习和筛选，构建和真实场景最接近的定标场地模型。利用这一定标场地模型，模拟出待定标卫星成像时刻对应观测几何下的表观反射率，实现传感器的绝对辐射定标。该方法虽然也取得了比较好的定标结果，但该方法对如何实现有效数据的筛选，并未给出详细的过程，且该方法是否适用于大角度影像对的交叉定标，也未得到有效验证。综上所述，现有交叉定标方法中，较少考虑大观测角度下卫星影像对之间方向反射率的差异，或采用地面实测反射率代替大气层顶表观反射率方向性的影响，导致现有交叉定标方法不适用于大观测角度卫星影像对的交叉定标。而在实际定标过程中，绝大多数交叉定标影像对的观测角度会有较大的差异。因此，提出一种可消除大观测角度差异的表观反射率模型，实现大角度下的高精度交叉定标，具有重要的研究意义和实际应用需求。
技术实现思路
本专利技术的目的在于克服现有技术的上述缺陷，提供一种基于MODIS影像的表观反射率模型构建方法，该方法可消除大观测角度差异，通过建立表观反射率模型获取表观反射率，进一步完成卫星在轨定标，具有定标精度高，可用于历史数据再定标等优点。本专利技术的另外一个目的在于提供一种基于MODIS影像的表观反射率模型构建系统。本专利技术的又一个目的在于提供一种基于MODIS影像的定标方法及定标系统。本专利技术的上述目的主要是通过如下技术方案予以实现的：基于MODIS影像的表观反射率模型构建方法，包括：获取定标场地时间序列MODIS卫星影像；根据定标场地位置，从所述MODIS卫星影像中提取定标场地上空对应的图像灰度值、图像定标系数和图像观测几何信息；根据所述图像定标系数，得到定标场地时间序列表观辐亮度、表观反射率和表观亮温；根据所述表观亮温得到包络线亮温，建立包络线亮温曲线，选取包络线亮温曲线中，包络线亮温和图像亮温的差值小于10K的图像；计算MODIS可见近红外区间第1波段的场地图像灰度值的变异系数，选取变异系数大于4％的图像；将同时满足包络线亮温和图像亮温的差值小于10K、变异系数大于4％的图像作为有效定标场地影像，得到所述有效定标场地影像对应的图像表观反射率和角度信息；根据所述有效定标场地影像对应的图像表观反射率以及角度信息，利用核驱动模型，得到表观反射率校正系数；将所述表观反射率校正系数代入所述核驱动模型，建立表观反射率模型，用于获取表观反射率。在上述基于MODIS影像的表观反射率模型构建方法中，所述定标场地的要求为：(1)场地面积不小于5公里×5公里；(2)定标场地均匀平坦、无植被覆盖；(3)场地晴天数多；(4)场地位于干旱区；获取定标场地1年以上的时间序列MOD021KM卫星影像。在上述基于MODIS影像的表观反射率模型构建方法中，所述图像灰度值包括MODIS可见近红外区间的波段1-7，和热红外波段31，32；所述图像定标系数包括各通道反射率定标系数和辐亮度定标系数；所述图像观测几何信息包括场地位置对应的太阳天顶角、太阳方位角、观测天顶角和观测方位角。在上述基于MODIS影像的表观反射率模型构建方法中，所述表观辐亮度计算公式如下：其中，Li为MODIS第i波段表观辐亮度，和分别为MODIS第i波段的表观辐亮度定标系数的增益和截距；所述表观反射率计算公式如下：其中，ρi为MODIS第i波段表观反射率，和分别为MODIS第i波段的表观反射率定标系数的增益和截距，θs为对应的太阳天顶角；DN为遥感图像的计数值；所述表观亮温计算公式如下：其中，T为表观亮温，h为普朗克常数，K为玻尔兹曼常数，c为光速，λ为波长，L为表观辐亮度。在上述基于MODIS影像的表观反射率模型构建方法中，所述计算MODIS本文档来自技高网...

【技术保护点】
1.基于MODIS影像的表观反射率模型构建方法，其特征在于：包括：获取定标场地时间序列MODIS卫星影像；根据定标场地位置，从所述MODIS卫星影像中提取定标场地上空对应的图像灰度值、图像定标系数和图像观测几何信息；根据所述图像定标系数，得到定标场地时间序列表观辐亮度、表观反射率和表观亮温；根据所述表观亮温得到包络线亮温，建立包络线亮温曲线，选取包络线亮温曲线中，包络线亮温和图像亮温的差值小于10K的图像；计算MODIS可见近红外区间第1波段的场地图像灰度值的变异系数，选取变异系数大于4％的图像；将同时满足包络线亮温和图像亮温的差值小于10K、变异系数大于4％的图像作为有效定标场地影像，得到所述有效定标场地影像对应的图像表观反射率和角度信息；根据所述有效定标场地影像对应的图像表观反射率以及角度信息，利用核驱动模型，得到表观反射率校正系数；将所述表观反射率校正系数代入所述核驱动模型，建立表观反射率模型，用于获取表观反射率。

【技术特征摘要】
1.基于MODIS影像的表观反射率模型构建方法，其特征在于：包括：获取定标场地时间序列MODIS卫星影像；根据定标场地位置，从所述MODIS卫星影像中提取定标场地上空对应的图像灰度值、图像定标系数和图像观测几何信息；根据所述图像定标系数，得到定标场地时间序列表观辐亮度、表观反射率和表观亮温；根据所述表观亮温得到包络线亮温，建立包络线亮温曲线，选取包络线亮温曲线中，包络线亮温和图像亮温的差值小于10K的图像；计算MODIS可见近红外区间第1波段的场地图像灰度值的变异系数，选取变异系数大于4％的图像；将同时满足包络线亮温和图像亮温的差值小于10K、变异系数大于4％的图像作为有效定标场地影像，得到所述有效定标场地影像对应的图像表观反射率和角度信息；根据所述有效定标场地影像对应的图像表观反射率以及角度信息，利用核驱动模型，得到表观反射率校正系数；将所述表观反射率校正系数代入所述核驱动模型，建立表观反射率模型，用于获取表观反射率。2.根据权利要求1所述的基于MODIS影像的表观反射率模型构建方法，其特征在于：所述定标场地的要求为：(1)场地面积不小于5公里×5公里；(2)定标场地均匀平坦、无植被覆盖；(3)场地晴天数多；(4)场地位于干旱区；获取定标场地1年以上的时间序列MOD021KM卫星影像。3.根据权利要求1所述的基于MODIS影像的定标场地表观反射率模型构建方法，其特征在于：所述图像灰度值包括MODIS可见近红外区间的波段1-7，和热红外波段31，32；所述图像定标系数包括各通道反射率定标系数和辐亮度定标系数；所述图像观测几何信息包括场地位置对应的太阳天顶角、太阳方位角、观测天顶角和观测方位角。4.根据权利要求1所述的基于MODIS影像的表观反射率模型构建方法，其特征在于：所述表观辐亮度计算公式如下：其中，Li为MODIS第i波段表观辐亮度，和分别为MODIS第i波段的表观辐亮度定标系数的增益和截距；所述表观反射率计算公式如下：其中，ρi为MODIS第i波段表观反射率，和分别为MODIS第i波段的表观反射率定标系数的增益和截距，θs为对应的太阳天顶角；DN为遥感图像的计数值；所述表观亮温计算公式如下：其中，T为表观亮温，h为普朗克常数，K为玻尔兹曼常数，c为光速，λ为波长，L为表观辐亮度。5.根据权利要求1所述的基于MODIS影像的表观反射率模型构建方法，其特征在于：所述计算MODIS可见近红外区间第1波段的场地图像灰度值的变异系数，即计算场地图像灰度值标准差与均值的比值；所述角度信息包括观测天顶角、观测方位角、太阳天顶角和太阳方位角。6.根据权利要求1所述的基于MODIS影像的表观反射率模型构建方法，其特征在于：所述核驱动模型表示如下：其中，为二向表观反射率；kgeo为几何光学核，kvol为体散射核；θs为太阳天顶角；θv为观测天顶角；为相对方位角，即太阳方位角与观测方位角的相对差值；fiso，fgeo，fvol为校正系数，分别表示各向均匀散射、几何...

【专利技术属性】
技术研发人员：刘李，高海亮，程天海，闵祥军，潘志强，张学文，
申请(专利权)人：中国资源卫星应用中心，
类型：发明
国别省市：北京,11