CN112665829B - 一种光学遥感卫星的波段间定标方法 - Google Patents
一种光学遥感卫星的波段间定标方法 Download PDFInfo
- Publication number
- CN112665829B CN112665829B CN202011492747.5A CN202011492747A CN112665829B CN 112665829 B CN112665829 B CN 112665829B CN 202011492747 A CN202011492747 A CN 202011492747A CN 112665829 B CN112665829 B CN 112665829B
- Authority
- CN
- China
- Prior art keywords
- calibration
- remote sensing
- optical remote
- sensing satellite
- camera
- 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.)
- Active
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Image Processing (AREA)
- Spectrometry And Color Measurement (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
本发明公开了一种光学遥感卫星的波段间定标方法,包括:确定光学遥感卫星各相机各波段间定标精度;在地面定标场测量光学遥感卫星相机的全色波段以及多光谱各波段的地表真实反演辐亮度值;针对每个所述地表真实反演辐亮度值,量测对应光学遥感卫星影像像元的表观辐亮度值,并计算地表真实反演辐亮度值与对应表观辐亮度值的线性关系;基于各所述线性关系确定各波段间的线性斜率误差。本发明公开的光学遥感卫星的波段间定标方法,能够实现光学遥感卫星的波段间定标。
Description
技术领域
本发明属于航天技术领域,尤其涉及一种光学遥感卫星的波段间定标方法。
背景技术
光学遥感卫星在轨期间,针对辐射定标,国内外学者主要是利用卫星-地面定标场同步测量数据,获取卫星相机的全色波段与多光谱各波段的真实反演辐亮度值与表观辐亮度值线性关系,进而求取各波段的绝对定标系数,以进行卫星数据的绝对辐射校正。而针对遥感卫星波段间的定标,尚未查找到有效文献与著作。
遥感卫星波段间的定标旨在获取辐射定标后,各波段绝对定标结果是否精确,能否达到业务系统使用要求,不仅可对单波段辐射定标结果进行评价,还可针对整星相机各波段间辐射标定后关系是否符合相机实验室要求。
发明内容
本发明要解决的技术问题为:提供一种光学遥感卫星的波段间定标方法,以获取不同波段绝对定标测量结果的比率值关系。
为了解决上述技术问题,本发明公开了一种光学遥感卫星的波段间定标方法,其中,所述方法包括:
确定光学遥感卫星各相机各波段间定标精度;
在地面定标场测量光学遥感卫星相机的全色波段以及多光谱各波段的地表真实反演辐亮度值;
针对每个所述地表真实反演辐亮度值,量测对应光学遥感卫星影像像元的表观辐亮度值,并计算地表真实反演辐亮度值与对应表观辐亮度值的线性关系;
基于各所述线性关系确定各波段间的线性斜率误差。
可选的,所述确定光学遥感卫星各相机各波段间定标精度的步骤,包括:
进行大气气象参数同步测量、典型地物及靶标光谱反射率同步测量,得到同步测量数据;
对所述同步测量数据进行预处理;
根据辐射定标流程算法,对预处理后的所述同步测量数据进行计算,得到相机的绝对辐射定标系数;
基于多次定标得到的相机绝对辐射定标系数、MODIS辐亮度比较验证结果,确定相机的绝对辐射定标精度。
可选的,在所述在地面定标场测量光学遥感卫星相机的全色波段以及多光谱各波段的地表真实反演辐亮度值的步骤之前,所述方法还包括:
判断所述相机的绝对辐射定标精度是否满足预设条件;
在满足预设条件的情况下,执行所述在地面定标场测量光学遥感卫星相机的全色波段以及多光谱各波段的地表真实反演辐亮度值的步骤。
可选的,在所述进行大气气象参数同步测量、典型地物及靶标光谱反射率同步测量,得到同步测量数据的步骤之前,所述方法还包括:
根据任务特点及需求制定外场辐射定标试验方案,其中,所述试验方案中包括:成像计划、试验日程、场地要求、仪器设备校准方案、靶标布设方案;
根据所述试验方案对定标场地进行实地勘探考察,选取合适的区域完成灰阶靶标的布设;
对实验所需的太阳光度计、光谱辐射计、参考板进行定标检测。
可选的,在所述对实验所需的太阳光度计、光谱辐射计、参考板进行定标检测的步骤之后,所述方法还包括:
建立光学遥感规划成像、太阳高度数据、地物光谱数据、气象信息、定位信息的同步测量机制,以实现星地同步获取光学遥感卫星数据、地面各测量数据的联动机制。
可选的,地表真实反演辐亮度值与对应表观辐亮度值的线性关系可通过如下公式表示:
L'i=(Li+Lbias)*K*Ki
其中,L'i为i波段的等效辐亮度值,Li为地表反射率反演得到的真实辐亮度值,Lbias为理论线性值的偏差,K为理论线性的斜率,Ki为理论线性斜率偏差。
本发明具有以下优点:
本发明实施例公开了一种光学遥感卫星的波段间定标方法,基于光学遥感卫星-地面定标场同步测量数据,在地面定标场测量光学遥感卫星相机的全色波段以及多光谱各波段的地表真实反演辐亮度值,通过量测对应光学遥感卫星影像像元的表观辐亮度值,建立地表真实反演辐亮度值与相应表观辐亮度值的线性关系,获得各波段间的线性斜率误差,实现光学遥感卫星的波段间定标,通过波段间定标可反应不同波段绝对定标测量结果的比率值关系。
附图说明
图1为本发明实施例的一种光学遥感卫星的波段间定标方法的步骤流程图;
图2为本发明实施例的又一种光学遥感卫星的波段间定标方法的步骤流程图。
具体实施方式
下面根据具体的实施例,结合附图针对本发明进行详细说明。应当理解,此处所述的具体实施例仅用于解释本发明,并不用于限定本发明。
本发明中基于光学遥感卫星-地面定标场同步测量数据,测量卫星相机的全色波段地表真实反演辐亮度值与表观辐亮度值线性关系,多光谱各波段的真实反演辐亮度值与表观辐亮度值线性关系,求取各波段间线性斜率误差。波段间定标反应不同波段绝对定标测量结果的比率值关系,基于此理念设计出如下述各实施例中所示的光学遥感卫星的波段间定标方法。
图1为本发明实施例提供的一种光学遥感卫星的波段间定位方法的步骤流程图。
本发明实施例的光学遥感卫星的波段间定标方法包括如下步骤:
步骤101:确定光学遥感卫星各相机各波段间定标精度。
本步骤中评判光学遥感卫星各相机各波段间定标精度,应用于卫星的在轨辐射标定。
一种可选地确定光学遥感卫星各相机各波段间定标精度的步骤,包括以下子步骤:
子步骤一:进行大气气象参数同步测量、典型地物及靶标光谱反射率同步测量,得到同步测量数据;
子步骤二:对所述同步测量数据进行预处理;
子步骤三:根据辐射定标流程算法,对预处理后的所述同步测量数据进行计算,得到相机的绝对辐射定标系数;
子步骤四:基于多次定标得到的相机绝对辐射定标系数、MODIS辐亮度比较验证结果,确定相机的绝对辐射定标精度。
MODIS(Moderate-resolution Imaging Spectroradiometer,中分辨率成像光谱仪)是搭载在TERRA和AQUA卫星上的一个重要的传感器,是卫星上唯一将实时观测数据通过x波段向全世界直接广播,并可以免费接收数据并无偿使用的星载仪器。
在一种可选地实施例中,在进行大气气象参数同步测量、典型地物及靶标光谱反射率同步测量,得到同步测量数据的步骤之前,还包括如下子步骤:
子步骤五:根据任务特点及需求制定外场辐射定标试验方案,其中,所述试验方案中包括:成像计划、试验日程、场地要求、仪器设备校准方案、靶标布设方案;
子步骤六根据所述试验方案对定标场地进行实地勘探考察,选取合适的区域完成灰阶靶标的布设;
子步骤七:对实验所需的太阳光度计、光谱辐射计、参考板进行定标检测。
在对实验所需的太阳光度计、光谱辐射计、参考板进行定标检测的步骤之后,方法还可以包括:子步骤八:建立光学遥感规划成像、太阳高度数据、地物光谱数据、气象信息、定位信息的同步测量机制,以实现星地同步获取光学遥感卫星数据、地面各测量数据的联动机制。
步骤102:在地面定标场测量光学遥感卫星相机的全色波段以及多光谱各波段的地表真实反演辐亮度值。
本发明实施例中,基于光学遥感卫星在轨期间,测量各相机各波段辐亮度值。
在一种可选地实施例中,在在地面定标场测量光学遥感卫星相机的全色波段以及多光谱各波段的地表真实反演辐亮度值的步骤之前,还包括如下步骤:
判断相机的绝对辐射定标精度是否满足预设条件;在满足预设条件的情况下,执行在地面定标场测量光学遥感卫星相机的全色波段以及多光谱各波段的地表真实反演辐亮度值的步骤。
预设条件可设定为大于预设精度阈值。
步骤103:针对每个地表真实反演辐亮度值,量测对应光学遥感卫星影像像元的表观辐亮度值,并计算地表真实反演辐亮度值与对应表观辐亮度值的线性关系。
其中,地表真实反演辐亮度值与对应表观辐亮度值的线性关系可通过如下公式表示:
L'i=(Li+Lbias)*K*Ki
其中,L'i为i波段的等效辐亮度值,Li为地表反射率反演得到的真实辐亮度值,Lbias为理论线性值的偏差,K为理论线性的斜率,Ki为理论线性斜率偏差。
步骤104:基于各线性关系确定各波段间的线性斜率误差。
所建立的各线性关系,能够反应不同波段绝对定标测量结果的比率值关系,达到光学遥感卫星辐射检校验证的要求。
本发明实施例公开了一种光学遥感卫星的波段间定标方法,基于光学遥感卫星-地面定标场同步测量数据,在地面定标场测量光学遥感卫星相机的全色波段以及多光谱各波段的地表真实反演辐亮度值,通过量测对应光学遥感卫星影像像元的表观辐亮度值,建立地表真实反演辐亮度值与相应表观辐亮度值的线性关系,获得各波段间的线性斜率误差,实现光学遥感卫星的波段间定标,通过波段间定标可反应不同波段绝对定标测量结果的比率值关系。
图2为本发明实施例提供的又一种光学遥感卫星的波段间定位方法的步骤流程图。
下面结合附图2对本发明的具体实施方式进行进一步的详细描述,定位方法整体可分为如下三部分:
第一部分:光学遥感卫星任务规划,同步准备地面定标场测量仪器及布设靶标。
首先根据任务特点及需求制定了外场辐射定标试验方案,明确了成像计划、试验日程、场地要求、仪器设备校准方案、靶标布设方案等内容,根据试验方案对定标场地进行实地勘探考察,选取合适的区域完成灰阶靶标的布设,同时对参试的太阳光度计、光谱辐射计、参考板等定标检测,确保仪器状态正常,以上准备工作完成后开展了卫星规划观测、大气气象参数同步测量、典型地物及靶标光谱反射率同步测量。
第二部分:对同步测量数据质量检查完成地面数据预处理,然后根据辐射定标流程算法计算出相机的绝对辐射定标系数,最后结合多次定标结果统计分析及MODIS辐亮度比较验证等完成绝对辐射定标精度分析评估。
第三部分:建立地表真实反演辐亮度值(真值)与相应表观辐亮度值的数学模型,获得各波段间的线性斜率误差。
定标后地表反射率反演可得到的真实辐亮度值(真值)与等效辐亮度间关系具体可通过以下公式得到:
L'i=(Li+Lbias)*K*Ki
其中,L'i为i波段的等效辐亮度值,Li为地表反射率反演得到的真实辐亮度值(真值),通过测量地表反射率反演得到,Lbias为理论线性值的偏差,K为理论线性的斜率,Ki为理论线性斜率偏差。
本实施例提供的光学遥感卫星的波段间定标方法,评判光学遥感卫星各相机各波段间定标精度,应用于卫星的在轨辐射标定。基于光学遥感卫星-地面定标场同步测量数据,建立地表真实反演辐亮度值(真值)与相应表观辐亮度值的数学模型,获得各波段间的线性斜率误差,该设计可应用于光学遥感卫星在轨辐射标定检校。反应不同波段绝对定标测量结果的比率值关系,达到光学遥感卫星辐射检校的要求。
需要说明的是,以上说明仅是本发明的优选实施方式,应当理解,对于本领域技术人员来说,在不脱离本发明技术构思的前提下还可以做出若干改变和改进,这些都包括在本发明的保护范围内。
本发明说明书中未作详细描述的内容属于本领域技术人员的公知技术。
Claims (5)
1.一种光学遥感卫星的波段间定标方法,其特征在于,所述方法包括:
确定光学遥感卫星各相机各波段间定标精度;
在地面定标场测量光学遥感卫星相机的全色波段以及多光谱各波段的地表真实反演辐亮度值;
针对每个所述地表真实反演辐亮度值,量测对应光学遥感卫星影像像元的表观辐亮度值,并计算地表真实反演辐亮度值与对应表观辐亮度值的线性关系;
基于各所述线性关系确定各波段间的线性斜率误差;
地表真实反演辐亮度值与对应表观辐亮度值的线性关系通过如下公式表示:
L'i=(Li+Lbias)*K*Ki
其中,L'i为i波段的表观辐亮度值,Li为地表反射率反演得到的地表真实反演辐亮度值,Lbias为理论线性值偏差,K为理论线性斜率,Ki为理论线性斜率误差。
2.根据权利要求1所述的方法,其特征在于,所述确定光学遥感卫星各相机各波段间定标精度的步骤,包括:
进行大气气象参数同步测量、典型地物及靶标光谱反射率同步测量,得到同步测量数据;
对所述同步测量数据进行预处理;
根据辐射定标流程算法,对预处理后的所述同步测量数据进行计算,得到相机的绝对辐射定标系数;
基于多次定标得到的相机绝对辐射定标系数、MODIS辐亮度比较验证结果,确定相机的绝对辐射定标精度。
3.根据权利要求2所述的方法,其特征在于,在所述在地面定标场测量光学遥感卫星相机的全色波段以及多光谱各波段的地表真实反演辐亮度值的步骤之前,所述方法还包括:
判断所述相机的绝对辐射定标精度是否满足预设条件;
在满足预设条件的情况下,执行所述在地面定标场测量光学遥感卫星相机的全色波段以及多光谱各波段的地表真实反演辐亮度值的步骤。
4.根据权利要求2所述的方法,其特征在于,在所述进行大气气象参数同步测量、典型地物及靶标光谱反射率同步测量,得到同步测量数据的步骤之前,所述方法还包括:
根据任务特点及需求制定外场辐射定标试验方案,其中,所述试验方案中包括:成像计划、试验日程、场地要求、仪器设备校准方案、靶标布设方案;
根据所述试验方案对定标场地进行实地勘探考察,选取合适的区域完成灰阶靶标的布设;
对实验所需的太阳光度计、光谱辐射计、参考板进行定标检测。
5.根据权利要求4所述的方法,其特征在于,在所述对实验所需的太阳光度计、光谱辐射计、参考板进行定标检测的步骤之后,所述方法还包括:
建立光学遥感规划成像、太阳高度数据、地物光谱数据、气象信息、定位信息的同步测量机制,以实现星地同步获取光学遥感卫星数据、地面各测量数据的联动机制。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011492747.5A CN112665829B (zh) | 2020-12-16 | 2020-12-16 | 一种光学遥感卫星的波段间定标方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011492747.5A CN112665829B (zh) | 2020-12-16 | 2020-12-16 | 一种光学遥感卫星的波段间定标方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112665829A CN112665829A (zh) | 2021-04-16 |
CN112665829B true CN112665829B (zh) | 2022-12-27 |
Family
ID=75404458
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011492747.5A Active CN112665829B (zh) | 2020-12-16 | 2020-12-16 | 一种光学遥感卫星的波段间定标方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112665829B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113447138B (zh) * | 2021-07-27 | 2023-05-23 | 郭强 | 一种星载遥感器热红外波段无辐射参考源定标方法 |
CN114646601A (zh) * | 2022-03-29 | 2022-06-21 | 中国科学技术大学 | 基于多波段耦合的卫星臭氧垂直廓线反演方法及系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11203443A (ja) * | 1998-01-19 | 1999-07-30 | Hitachi Ltd | マルチスペクトル衛星画像の処理方法と処理システム並びに水圏評価方法 |
CN102288956A (zh) * | 2011-05-10 | 2011-12-21 | 中国资源卫星应用中心 | 一种遥感卫星多光谱数据的大气订正方法 |
CN102901516A (zh) * | 2012-09-29 | 2013-01-30 | 航天恒星科技有限公司 | 一种基于绝对辐射定标的多光谱影像辐射校正方法 |
CN107230186A (zh) * | 2016-03-23 | 2017-10-03 | 中国科学院遥感与数字地球研究所 | 一种卫星遥感影像物理匀色方法 |
CN108955883A (zh) * | 2018-07-19 | 2018-12-07 | 国家卫星气象中心 | 基于高光谱数据的遥感器辐射定标方法及定标系统 |
CN110411565A (zh) * | 2019-07-05 | 2019-11-05 | 中国船舶重工集团公司第七一六研究所 | 一种由多光谱地表反射率推演高光谱地表反射率的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109974665B (zh) * | 2019-04-04 | 2021-05-07 | 东北师范大学 | 一种针对缺少短波红外数据的气溶胶遥感反演方法及系统 |
-
2020
- 2020-12-16 CN CN202011492747.5A patent/CN112665829B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11203443A (ja) * | 1998-01-19 | 1999-07-30 | Hitachi Ltd | マルチスペクトル衛星画像の処理方法と処理システム並びに水圏評価方法 |
CN102288956A (zh) * | 2011-05-10 | 2011-12-21 | 中国资源卫星应用中心 | 一种遥感卫星多光谱数据的大气订正方法 |
CN102901516A (zh) * | 2012-09-29 | 2013-01-30 | 航天恒星科技有限公司 | 一种基于绝对辐射定标的多光谱影像辐射校正方法 |
CN107230186A (zh) * | 2016-03-23 | 2017-10-03 | 中国科学院遥感与数字地球研究所 | 一种卫星遥感影像物理匀色方法 |
CN108955883A (zh) * | 2018-07-19 | 2018-12-07 | 国家卫星气象中心 | 基于高光谱数据的遥感器辐射定标方法及定标系统 |
CN110411565A (zh) * | 2019-07-05 | 2019-11-05 | 中国船舶重工集团公司第七一六研究所 | 一种由多光谱地表反射率推演高光谱地表反射率的方法 |
Non-Patent Citations (2)
Title |
---|
卫星光学传感器全过程辐射定标;乔延利等;《遥感学报》;20060930;第10卷(第5期);第3节 * |
遥感25号无场化相对辐射标定;张过等;《测绘学报》;20170831;第46卷(第8期);第1节 * |
Also Published As
Publication number | Publication date |
---|---|
CN112665829A (zh) | 2021-04-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102901516B (zh) | 一种基于绝对辐射定标的多光谱影像辐射校正方法 | |
Lacherade et al. | Cross calibration over desert sites: Description, methodology, and operational implementation | |
Coll et al. | Temperature and emissivity separation from calibrated data of the Digital Airborne Imaging Spectrometer | |
CN112798013B (zh) | 一种对光学载荷在轨绝对辐射定标结果进行验证的方法 | |
CN112665829B (zh) | 一种光学遥感卫星的波段间定标方法 | |
Clevers et al. | MERIS and the red-edge position | |
Sobrino et al. | Accuracy of ASTER level-2 thermal-infrared standard products of an agricultural area in Spain | |
CN102279393A (zh) | 一种基于多光谱传感器对高光谱传感器交叉辐射定标方法 | |
Bendig et al. | Solar-induced chlorophyll fluorescence measured from an unmanned aircraft system: sensor etaloning and platform motion correction | |
CN114219994A (zh) | 一种基于空海协同观测的海洋光学卫星辐射定标方法 | |
CN113066057A (zh) | 一种气溶胶光学厚度监测方法 | |
Valorge et al. | Forty years of experience with SPOT in-flight calibration | |
Markelin et al. | Radiometric calibration and characterization of large-format digital photogrammetric sensors in a test field | |
CN110702228B (zh) | 一种航空高光谱影像的边缘辐射校正方法 | |
CN105737980B (zh) | 辐射定标方法及装置 | |
CN112924391A (zh) | 基于遥感大数据的fy-4a/agri交叉辐射定标方法 | |
CN113836731B (zh) | 陆表稳定目标大气层顶反射率模型的构建方法和装置 | |
CN113029977B (zh) | 一种针对宽视场角多光谱传感器的自动交叉辐射定标方法 | |
CN114296061A (zh) | 基于多元变量检测和不同辐射传输模型的交叉定标方法 | |
Slater et al. | In-flight radiometric calibration of ASTER by reference to well-characterized scenes | |
Babu et al. | An empirical comparison of calibration and validation methodologies for airborne imaging spectroscopy | |
Salaswati et al. | Potential Area for Radiometric Calibration of Satellite Cameras in Indonesia | |
Valor et al. | The Adjusted Normalized Emissivity Method (ANEM) for land surface temperature and emissivity recovery | |
CN116794748A (zh) | 一种遥感卫星的多时相绝对辐射定标方法 | |
Wang et al. | Consistent Transfer Radiometric Calibration Technology for Optical Remotely Sensor and Field Campaign Validation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |