CN104535979B - 一种陆地云光学厚度的遥感反演方法及系统 - Google Patents

一种陆地云光学厚度的遥感反演方法及系统 Download PDF

Info

Publication number
CN104535979B
CN104535979B CN201410811094.0A CN201410811094A CN104535979B CN 104535979 B CN104535979 B CN 104535979B CN 201410811094 A CN201410811094 A CN 201410811094A CN 104535979 B CN104535979 B CN 104535979B
Authority
CN
China
Prior art keywords
cloud
visible image
look
image
optical thickness
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.)
Expired - Fee Related
Application number
CN201410811094.0A
Other languages
English (en)
Other versions
CN104535979A (zh
Inventor
张海龙
李丽
辛晓洲
余珊珊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Remote Sensing and Digital Earth of CAS
Original Assignee
Institute of Remote Sensing and Digital Earth of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Institute of Remote Sensing and Digital Earth of CAS filed Critical Institute of Remote Sensing and Digital Earth of CAS
Priority to CN201410811094.0A priority Critical patent/CN104535979B/zh
Publication of CN104535979A publication Critical patent/CN104535979A/zh
Application granted granted Critical
Publication of CN104535979B publication Critical patent/CN104535979B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
    • G01S7/4802—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/88—Lidar systems specially adapted for specific applications
    • G01S17/89—Lidar systems specially adapted for specific applications for mapping or imaging
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/88—Lidar systems specially adapted for specific applications
    • G01S17/95—Lidar systems specially adapted for specific applications for meteorological use
    • 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

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Electromagnetism (AREA)
  • Image Processing (AREA)

Abstract

本发明提供了一种陆地云光学厚度的遥感反演方法及系统，该方法包括：获取陆地的可见光影像，生成所述可见光影像的地表反射率；根据预设参数利用辐射传输模型建立查找表；对所述可见光影像进行辐射定标，转换为可见光影像的星上辐照度；结合多光谱影像数据，将所述可见光影像生成云掩膜图像；根据所述可见光影像的地表反射率和星上辐照度，利用所述查找表，将所述云掩膜图像中的云覆盖像元的星上辐照度反演为云光学厚度。本发明提出的陆地云光学厚度的反演方法及系统，基于云掩膜影像的提取精度，可以快速准确的实现日间云光学厚度的提取，且独立于特定的传感器和卫星，在高亮地表，如冰、雪等下垫面情形下，取得较好的反演结果。

Description

一种陆地云光学厚度的遥感反演方法及系统
技术领域
本发明涉及大气遥感技术领域，尤其涉及一种陆地云光学厚度的遥感反演方法及系统。
背景技术
云性质及其在时空尺度上的变化对全球气候变化研究十分重要。云光学厚度是计算云辐射效应和光学作用的重要参量之一，是地气辐射系统的主要调节参量。遥感技术为大范围实时云光学参数反演提供了有力手段。在云光学厚度大的情况下，云的反射函数可写为符合渐进理论的表达形式，其中反射辐射在非吸收的大气波段是云光学厚度、地面反射率和非对称因子的函数。当云光学厚度大于9后，无论是非守恒大气还是守恒大气，反射函数的渐进理论表达式均可以适用。Ou等人在辐射传输模式和参数化的基础上，选取了AVHRR的3.7μm和10.9μm两个通道数据，提出了反演卷云光学厚度的算法，并将该算法应用于第一次国际卫星云气候计划区域试验卷云加密观测中，结果表明，反演得到的云参数与地基和机载探测结果较为一致。Rosenfeld等人利用NOAA-AVHRR的0.65通道数据反演了云顶附近的云光学厚度信息，并与雷达回波的结果进行了对比分析。第一个全球范围的光学厚度产品由ISCCP提供，该方法是针对有效粒子半径为10μm的水云，通过查找表的方式得到，其光学厚度值域为0.5～100。
20世纪80年代以来，我国利用飞机、雷达和气象卫星遥感资料，开展了一系列云辐射特性方面的研究工作，取得了众多研究成果。利用卫星可见光通道资料反演云光学厚度的基础理论已趋于成熟，即在可见光波段，云的反射函数主要依赖于云的光学厚度，但受多种因素的共同影响，云的辐射特性复杂，在空间和时间上具有高度的可变性。现有MODIS产品的业务算法也采用近红外和可见光两个通道同时实现光学厚度和有效粒子半径的反演，其理论基础依然是平面平行理论，其反演方法会高估向阳方向的云光学厚度，低估背光方向的云光学厚度。
在实现本发明过程中，发明人发现现有技术中至少存在如下问题：
目前而言，卫星遥感云光学厚度的精度还存在着一定的问题，其中最主要的问题是云的宏观和微观不均匀性对云光学厚度的影响的程度至今尚无准确结论。对冰云而言，球形粒子假定和Mie理论已不再适用，加上冰云的时空分布变化大和其水平取向的多样性，使得冰云的反演难度更大。Nakajima指出，对于光学厚度大于4，有效粒子半径大于6μm的层状水云来说，可以用0.75μm反演光学厚度，而对于光学薄云层来说，反演结果是不确定的，解不是唯一的，由0.75反演的云光学厚度代表整层云光学厚度的20％～40％。
发明内容
(一)要解决的技术问题
本发明提出一种陆地云光学厚度的遥感反演方法及系统，可以快速准确的实现日间云光学厚度的提取，对于高亮地表，如冰、雪等下垫面情形，得到较好的反演结果。
(二)技术方案
为解决上述技术问题，本发明提供了一种陆地云光学厚度的遥感反演方法，所述方法包括：
获取陆地的可见光影像，生成所述可见光影像的地表反射率；
根据预设参数利用辐射传输模型建立查找表；
对所述可见光影像进行辐射定标，转换为可见光影像的星上辐照度；
结合多光谱影像数据，将所述可见光影像生成云掩膜图像；
根据所述可见光影像的地表反射率和星上辐照度，利用所述查找表，将所述云掩膜图像中的云覆盖像元的星上辐照度反演为云光学厚度。
优选地，所述可见光影像为极轨卫星采集的可见光波段影像数据或静止卫星采集的可见光波段影像数据。
优选地，所述生成所述可见光影像的地表反射率数据具体包括：
在预设时间段内，选取逐日同一成像时刻的可见光影像；
对所述预设时间段内选取的所述可见光影像进行大气校正；
根据校正后的可见光影像，提取地表反射率的最小值作为地表反射率。
优选地，所述预设时间段的时间步长为一个月，在所述预设时间段内地表反射率特性变化最小，且该预设时间段内至少有一天是晴空状况。
优选地，所述根据预设参数利用辐射传输模型建立查找表，具体包括：
根据预设参数，所述预设参数包括多个太阳天顶角、多个观测天顶角、相对方位角、多个地表高程、多个气溶胶光学厚度和多个云光学厚度，采用大气辐射传输模型SBDART计算多组不同地表、大气和观测几何下的云光学厚度与星上辐照度的值，生成查找表。
优选地，根据云相态的不同，所述查找表包括：水云查找表和冰云查找表。
优选地，结合多光谱影像数据，生成云掩膜图像具体包括：
通过多光谱影像数据不同波段的波段反射率或亮温与背景值的差值进行比较，预设不同的差值阈值，利用光谱阈值法生成第一层云掩膜；
通过对多光谱影像划分不同的像元矩阵，在划分出的每块像元矩阵内设定不同的动态阈值，高于所述动态阈值的像元矩阵判定为有云像元，生成第二层云掩膜；
利用云的高反射率和低亮温特性，针对可见光波段和热红外波段设定不同的覆盖阈值，高于所述覆盖阈值的像元判定为云覆盖像元，生成第三层云掩膜；
将上述三层云掩膜进行图像融合，生成云掩膜图像。
优选地，根据所述可见光影像的地表反射率和星上辐照度，利用所述查找表，将所述云掩膜图像中的云覆盖像元的星上辐照度反演为云光学厚度，具体包括：
对所述可见光影像进行地形识别，识别出下垫面为陆地的区域；
获取所述区域的地表高程和气溶胶光学厚度；
计算当前卫星的观测天顶角与观测方位角以及太阳天顶角与太阳方位角，根据观测方位角和太阳方位角计算相对方位角；
根据3.7μm通道的亮温进行云相态识别，选择相应的查找表；
根据所述可见光影像的地表反射率、星上辐照度、地表高程、观测天顶角、太阳天顶角、相对方位角以及气溶胶光学厚度，在所述查找表中选取对应的光学厚度；
利用线性插值法获取所述云掩膜图像中的云覆盖像元的光学厚度。
优选地，所述根据3.7μm通道的亮温进行云相态识别，选择相应的查找表，具体包括：
如果所述3.7μm通道的亮温高于238K，则选用水云查找表进行云参数反演；如果所述3.7μm通道的亮温低于或等于238K，则选用冰云查找表进行云参数反演。
相应的，本发明还提出了一种陆地云光学厚度的遥感反演系统，所述系统包括：
地表反射率生成模块，用于获取陆地的可见光影像，生成所述可见光影像的地表反射率；
查找表建立模块，用于根据预设参数利用辐射传输模型建立查找表；
辐射定标模块，用于对所述可见光影像进行辐射定标，转换为可见光影像的星上辐照度；
云掩膜图像生成模块，用于结合多光谱影像数据，将所述可见光影像生成云掩膜图像；
反演模块，用于根据所述可见光影像的地表反射率和星上辐照度，利用所述查找表，将所述云掩膜图像中的云覆盖像元的星上辐照度反演为云光学厚度。
(三)有益效果
本发明提供了一种陆地云光学厚度的遥感反演方法及系统，基于云掩膜影像的提取，可以快速准确的实现日间云光学厚度的提取，且独立于特定的传感器和卫星，在高亮地表，如冰、雪等下垫面情形下，取得较好的反演结果。
附图说明
通过参考附图会更加清楚的理解本发明的特征和优点，附图是示意性的而不应理解为对本发明进行任何限制，在附图中：
图1为本发明提出的一种陆地云光学厚度的遥感反演方法流程图；
图2为本发明提出的一种陆地云光学厚度的遥感反演系统模块图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
图1为本发明提出的一种陆地云光学厚度的遥感反演方法流程图，如图1所示，该方法包括：
S101，获取陆地的可见光影像，生成所述可见光影像的地表反射率；所述可见光影像为极轨卫星采集的可见光波段影像数据或静止卫星采集的可见光波段影像数据。
其中，生成所述可见光影像的地表反射率数据具体包括：在预设时间段内，选取逐日同一成像时刻的可见光影像；对所述预设时间段内选取的所述可见光影像进行大气校正；根据校正后的可见光影像，提取地表反射率的最小值作为地表反射率，所述预设时间段的时间步长为一个月，在所述预设时间段内地表反射率特性变化最小，且该预设时间段内至少有一天是晴空状况。
本发明实施例中，结合遥感影像可见光通道以及辐射传输模型，综合考虑地表反射、大气层的反射、日-地-星间的几何关系等参数，求算出可见光波段的地表反射率。以一个月的时间窗口作为步长，将连续一年每个月高于0.005的最低地表反射率作为晴空背景值，并录入数据库。其中，高反射地表(大于0.4)的像元标识为冰雪覆盖区。
S102，根据预设参数利用辐射传输模型建立查找表，具体包括：
根据预设参数，所述预设参数包括多个太阳天顶角、多个观测天顶角、相对方位角、多个地表高程、多个气溶胶光学厚度和多个云光学厚度，采用大气辐射传输模型SBDART计算多组不同地表、大气和观测几何下的云光学厚度与星上辐照度的值，生成查找表，根据云相态的不同，所述查找表包括：水云查找表和冰云查找表。
本发明实施例中，通过设置不同日-地-星间的几何参数(包括19个太阳天顶角、19个观测天顶角和19个相对方位角)，7个地表高程、11个地表反射率、7个气溶胶光学厚度、13个云光学厚度值，采用SBDART辐射传输模型计算多组不同地表、大气和观测几何下的云光学厚度与星上辐照度的值而成N×M维的查找表。根据云相态的不同，分别针对水云和冰云构建不同的查找表。
S103，对所述可见光影像进行辐射定标，转换为可见光影像的星上辐照度；
本发明实施例中，利用辅助文件中提供的定标系数，将影像DN值转换为星上辐照度R(W/m2.sr)，具体如下：
R＝DN×gain+bias
其中，gain为增益量，bias为偏移量，定标系数根据静止或极轨卫星以及仪器灵敏度的不同而会不定期的进行更新。
S104，结合多光谱影像数据，将所述可见光影像生成云掩膜图像，具体包括：
通过多光谱影像数据不同波段的波段反射率或亮温与背景值的差值进行比较，预设不同的差值阈值，利用光谱阈值法生成第一层云掩膜；通过对多光谱影像划分不同的像元矩阵，在划分出的每块像元矩阵内设定不同的动态阈值，高于所述动态阈值的像元矩阵判定为有云像元，生成第二层云掩膜；利用云的高反射率和低亮温特性，针对可见光波段和热红外波段设定不同的覆盖阈值，高于所述覆盖阈值的像元判定为云覆盖像元，生成第三层云掩膜；将上述三层云掩膜进行图像融合，生成云掩膜图像。
本发明实施例中，为保证云掩膜图像的精度，结合可见光、热红外波段3.7μm和12.0μm波段影像进行云识别。首先根据不同波段，通过该波段反射率或亮温与背景值的差值进行比较，根据不同的阈值，利用光谱阈值法生成第一层云掩膜；通过对影像划分不同的像元矩阵，针对每块子区域内设定不同的动态阈值，高于该阈值的像元判定为有云像元，生成第二层云掩膜；利用云的高反射率和低亮温特性，针对可见光波段和热红外波段设定不同的阈值，高于该阈值的像元判定为云覆盖像元，生成第三层云掩膜。三层云掩膜图像融合即为区域的云掩膜影像。
S105，根据所述可见光影像的地表反射率和星上辐照度，利用所述查找表，将所述云掩膜图像中的云覆盖像元的星上辐照度反演为云光学厚度，具体包括：
对所述可见光影像进行地形识别，识别出下垫面为陆地的区域；
获取所述区域的地表高程和气溶胶光学厚度；可见光影像对应的像元区域的地表高程和气溶胶光学厚度为现有数据。
计算当前卫星的观测天顶角与观测方位角以及太阳天顶角与太阳方位角，根据观测方位角和太阳方位角计算相对方位角；
根据3.7μm通道的亮温并进行云相态识别，选择相应的查找表；
根据所述可见光影像的地表反射率、星上辐照度、地表高程、观测天顶角、太阳天顶角、相对方位角以及气溶胶光学厚度，在所述查找表中选取对应的光学厚度；
利用线性插值法获取所述云掩膜图像中的云覆盖像元的光学厚度。
其中，根据3.7μm通道的亮温进行云相态识别，选择相应的查找表，具体包括：如果所述3.7μm通道的亮温高于238K，则选用水云查找表进行云参数反演；如果所述3.7μm通道的亮温低于或等于238K，则选用冰云查找表进行云参数反演。
本发明实施例中，采用的星上辐照度转光学厚度(τc)的算法如Lambertain面上的光学厚度与反射函数的表达式：
式中，K为逃逸函数，Ag为地面反射率，g为非对称因子，μ为观测天顶角余弦，μ0为太阳天顶角余弦，φ为相对方位角，k为散射指数(特征向量)，A*为半无限大气的球面反照率，m、n、l为常数。上述k、A*、m、n以及l等五个参数与单次散射反照率(ω0)紧密相关，与非对称因子(g)关系较小，并且可以由相似参数s表示，具体如下：
云层的反射特性主要由光学厚度指数和相似参数确定，同时还与地表反射率有关。而在波长小于1μm的可见光波段，相似参数s近乎为0，因此云光学厚度可以由这一波段的反射函数获得。
本发明实施例中，根据地形识别结果，判断下垫面情形，区分陆地、内陆水体和海洋，只对标记为陆地的区域进行像元计算，提取地表高程。
本发明实施例中，以静止卫星为例，详细介绍计算卫星天顶角与方位角，具体如下：
静止卫星的天顶角Zenv可以表示为：
Zenv＝π/2-tan-1[(rs/re-cosγ)/sinγ]-γ
cos(γ)＝cos(Le)cos(Ls)cos(ls-le)+sin(Le)sin(Ls)
式中，rs为卫星高度，re为地球半径，Ls为静止卫星所处的纬度，Le为像元某点所处的纬度，ls为静止卫星所处经度，le为像元某点所处经度。
静止卫星的方位角可以表示为：
卫星方位角定义为从正北方向起，顺时针旋转，则各个象限的方位角为：
静止卫星的天顶角和方位角是固定的，因此只需要计算一次即可。
本发明实施例中，计算太阳天顶角与方位角，具体如下：
根据太阳视轨道方程，太阳高度角Zens可以表示为：
式中，为地理纬度，δ为太阳赤纬，ω为太阳时角。
太阳方位角Azis的计算为：
由于有云天气下水汽、臭氧等参数对可见光波段的影响远低于云的影响，在本算法中都取默认值。根据读取的太阳天顶角、太阳方位角、观测天顶角、观测方位角、地表高程和气溶胶光学厚度，结合得到的云相态以及得到的地表反射率，在查找表中选取相应的数据，通过对查找表不断降维处理的方法，采用线性插值得到整幅图像云覆盖像元的光学厚度。
图2为本发明提出的一种陆地云光学厚度的遥感反演系统模块图，如图2所示，该系统包括：
地表反射率生成模块201，用于获取陆地的可见光影像，生成所述可见光影像的地表反射率；
查找表建立模块202，用于根据预设参数利用辐射传输模型建立查找表；
辐射定标模块203，用于对所述可见光影像进行辐射定标，转换为可见光影像的星上辐照度；
云掩膜图像生成模块204，用于结合多光谱影像数据，将所述可见光影像生成云掩膜图像；
反演模块205，用于根据所述可见光影像的地表反射率和星上辐照度，利用所述查找表，将所述云掩膜图像中的云覆盖像元的星上辐照度反演为云光学厚度。
采用本发明提供了陆地云光学厚度的遥感反演方法及系统，可以快速准确的实现日间云光学厚度的提取，且独立于特定的传感器和卫星，对于高亮地表，如冰、雪等下垫面情形下的反演结果较好。其反演精度主要依赖于云掩膜影像的提取精度，可以实现太阳天顶角在75°以下像元云光学厚度的自动反演提取。
以上实施方式仅用于说明本发明，而并非对本发明的限制，有关技术领域的普通技术人员，在不脱离本发明的精神和范围的情况下，还可以做出各种变化和变型，因此所有等同的技术方案也属于本发明的范畴，本发明的发明保护范围应由权利要求限定。

Claims (10)

1.一种陆地云光学厚度的遥感反演方法，其特征在于，所述方法包括：
获取陆地的可见光影像，生成所述可见光影像的地表反射率；
根据预设参数利用辐射传输模型建立查找表；
对所述可见光影像进行辐射定标，转换为可见光影像的星上辐照度；
结合多光谱影像数据，将所述可见光影像生成云掩膜图像；
根据所述可见光影像的地表反射率和星上辐照度，利用所述查找表，将所述云掩膜图像中的云覆盖像元的星上辐照度反演为云光学厚度。
2.根据权利要求1所述的方法，其特征在于，所述可见光影像为极轨卫星采集的可见光波段影像数据或静止卫星采集的可见光波段影像数据。
3.根据权利要求1或2所述的方法，其特征在于，所述生成所述可见光影像的地表反射率具体包括：
在预设时间段内，选取逐日同一成像时刻的可见光影像；
对所述预设时间段内选取的所述可见光影像进行大气校正；
根据校正后的可见光影像，提取地表反射率的最小值作为地表反射率。
4.根据权利要求3所述的方法，其特征在于，所述预设时间段的时间步长为一个月，在所述预设时间段内地表反射率特性变化最小，且该预设时间段内至少有一天是晴空状况。
5.根据权利要求1所述的方法，其特征在于，所述根据预设参数利用辐射传输模型建立查找表，具体包括：
根据预设参数，所述预设参数包括多个太阳天顶角、多个观测天顶角、相对方位角、多个地表高程、多个气溶胶光学厚度和多个云光学厚度，采用大气辐射传输模型SBDART计算多组不同地表、大气和观测几何下的云光学厚度与星上辐照度的值，生成查找表。
6.根据权利要求5所述的方法，其特征在于，根据云相态的不同，所述查找表包括：水云查找表和冰云查找表。
7.根据权利要求1所述的方法，其特征在于，结合多光谱影像数据，生成云掩膜图像具体包括：
通过多光谱影像数据不同波段的波段反射率或亮温与背景值的差值进行比较，预设不同的差值阈值，利用光谱阈值法生成第一层云掩膜；
通过对多光谱影像划分不同的像元矩阵，在划分出的每块像元矩阵内设定不同的动态阈值，高于所述动态阈值的像元矩阵判定为有云像元，生成第二层云掩膜；
利用云的高反射率和低亮温特性，针对可见光波段和热红外波段设定不同的覆盖阈值，高于所述覆盖阈值的像元判定为云覆盖像元，生成第三层云掩膜；
将上述三层云掩膜进行图像融合，生成云掩膜图像。
8.根据权利要求1所述的方法，其特征在于，根据所述可见光影像的地表反射率和星上辐照度，利用所述查找表，将所述云掩膜图像中的云覆盖像元的星上辐照度反演为云光学厚度，具体包括：
对所述可见光影像进行地形识别，识别出下垫面为陆地的区域；
获取所述区域的地表高程和气溶胶光学厚度；
计算当前卫星的观测天顶角与观测方位角以及太阳天顶角与太阳方位角，根据观测方位角和太阳方位角计算相对方位角；
根据3.7μm通道的亮温进行云相态识别，选择相应的查找表；
根据所述可见光影像的地表反射率、星上辐照度、地表高程、观测天顶角、太阳天顶角、相对方位角以及气溶胶光学厚度，在所述查找表中选取对应的光学厚度；
利用线性插值法获取所述云掩膜图像中的云覆盖像元的光学厚度。
9.根据权利要求8所述的方法，其特征在于，所述根据3.7μm通道的亮温进行云相态识别，选择相应的查找表，具体包括：
如果所述3.7μm通道的亮温高于238K，则选用水云查找表进行云参数反演；如果所述3.7μm通道的亮温低于或等于238K，则选用冰云查找表进行云参数反演。
10.一种陆地云光学厚度的遥感反演系统，其特征在于，所述系统包括：
地表反射率生成模块，用于获取陆地的可见光影像，生成所述可见光影像的地表反射率；
查找表建立模块，用于根据预设参数利用辐射传输模型建立查找表；
辐射定标模块，用于对所述可见光影像进行辐射定标，转换为可见光影像的星上辐照度；
云掩膜图像生成模块，用于结合多光谱影像数据，将所述可见光影像生成云掩膜图像；
反演模块，用于根据所述可见光影像的地表反射率和星上辐照度，利用所述查找表，将所述云掩膜图像中的云覆盖像元的星上辐照度反演为云光学厚度。
CN201410811094.0A 2014-12-23 2014-12-23 一种陆地云光学厚度的遥感反演方法及系统 Expired - Fee Related CN104535979B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410811094.0A CN104535979B (zh) 2014-12-23 2014-12-23 一种陆地云光学厚度的遥感反演方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410811094.0A CN104535979B (zh) 2014-12-23 2014-12-23 一种陆地云光学厚度的遥感反演方法及系统

Publications (2)

Publication Number Publication Date
CN104535979A CN104535979A (zh) 2015-04-22
CN104535979B true CN104535979B (zh) 2017-03-29

Family

ID=52851541

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410811094.0A Expired - Fee Related CN104535979B (zh) 2014-12-23 2014-12-23 一种陆地云光学厚度的遥感反演方法及系统

Country Status (1)

Country Link
CN (1) CN104535979B (zh)

Families Citing this family (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105911613B (zh) * 2016-07-01 2018-07-24 中国科学院地理科学与资源研究所 一种自动确定全天无云和部分有云天气情况的方法
FR3066301B1 (fr) * 2017-05-15 2022-02-25 Airbus Ds Geo Sa Procede de traitement d’une image optique acquise en presence d’aerosols et/ou de nuages dans l’atmosphere
CN110163035A (zh) * 2018-02-11 2019-08-23 青岛星科瑞升信息科技有限公司 一种先验数据支持的云阴影识别方法
CN109323992B (zh) * 2018-09-19 2021-03-09 中国科学院遥感与数字地球研究所 一种大气乙二醛遥感监测方法
CN109543721B (zh) * 2018-11-05 2023-07-04 中国科学院西北生态环境资源研究院 一种晴间多云天气情况下的太阳辐照度超短期预报方法
CN110186822B (zh) * 2019-05-13 2020-09-11 中国科学院遥感与数字地球研究所 一种气溶胶光学厚度遥感反演方法
CN110287587B (zh) * 2019-06-24 2020-10-13 中国科学院遥感与数字地球研究所 一种确定污染云微物理属性的方法
CN110455689B (zh) * 2019-08-01 2020-12-29 中国科学院遥感与数字地球研究所 一种模拟冰晶粒子的光散射特性的方法
CN110411927B (zh) * 2019-08-02 2020-12-04 中国科学院遥感与数字地球研究所 一种大气细粒子aod和地表偏振反射率协同反演方法
CN110543837B (zh) * 2019-08-16 2022-03-15 北京航空航天大学 一种基于潜在目标点的可见光机场飞机检测方法
CN110688606B (zh) * 2019-08-30 2021-10-12 中国科学院遥感与数字地球研究所 一种热红外遥感反演云微物理参数的方法
CN110726653B (zh) * 2019-09-25 2022-03-04 中国电子科技集团公司第二十七研究所 基于天地一体化信息的pm2.5浓度监测方法
CN110633693A (zh) * 2019-09-26 2019-12-31 南京信息工程大学 一种基于云反演资料的深对流系统自动识别追踪方法
CN111123382B (zh) * 2019-12-25 2020-11-13 中国科学院遥感与数字地球研究所 一种气溶胶和地表参数联合反演方法
CN111598802B (zh) * 2020-05-12 2023-04-25 中国科学院合肥物质科学研究院 一种地基全天空云参数反演系统和方法
CN113093183B (zh) * 2021-04-02 2022-03-08 国家卫星海洋应用中心 阈值确定方法、海冰厚度反演方法、装置、设备及介质
CN114112070B (zh) * 2021-12-29 2022-06-24 国家卫星海洋应用中心 一种星载微波辐射计的外定标方法及装置
CN115267786A (zh) * 2022-07-29 2022-11-01 北京彩彻区明科技有限公司 融合卫星观测亮温和高程的ResUNet-GAN全球雷达反演方法及装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101915914A (zh) * 2010-07-30 2010-12-15 南京信息工程大学 一种基于查找表的遥感影像逐像元大气校正方法
CN102103204A (zh) * 2011-01-26 2011-06-22 环境保护部卫星环境应用中心 基于环境一号卫星的陆地气溶胶光学厚度反演方法
CN102636143A (zh) * 2012-04-20 2012-08-15 中国科学院遥感应用研究所 一种气溶胶光学厚度遥感反演方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101915914A (zh) * 2010-07-30 2010-12-15 南京信息工程大学 一种基于查找表的遥感影像逐像元大气校正方法
CN102103204A (zh) * 2011-01-26 2011-06-22 环境保护部卫星环境应用中心 基于环境一号卫星的陆地气溶胶光学厚度反演方法
CN102636143A (zh) * 2012-04-20 2012-08-15 中国科学院遥感应用研究所 一种气溶胶光学厚度遥感反演方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
利用MODIS数据反演多层云光学厚度和有效粒子半径;叶晶 等;《气象学报》;20091231(第4期);613-622 *
利用透射太阳辐射反演云光学厚度及有效粒子半径：方法研究;王越 等;《自然科学进展》;20060731;第16卷(第7期);850-858 *

Also Published As

Publication number Publication date
CN104535979A (zh) 2015-04-22

Similar Documents

Publication Publication Date Title
CN104535979B (zh) 一种陆地云光学厚度的遥感反演方法及系统
Li et al. Temperature-based and radiance-based validation of the collection 6 MYD11 and MYD21 land surface temperature products over barren surfaces in northwestern China
Li et al. Watershed allied telemetry experimental research
Li et al. Comparison of the MuSyQ and MODIS Collection 6 land surface temperature products over barren surfaces in the Heihe River Basin, China
West et al. Cassini Imaging Science Subsystem observations of Titan’s south polar cloud
Walden et al. Atmospheric components of the surface energy budget over young sea ice: Results from the N‐ICE2015 campaign
CN107436193B (zh) 一种基于遥感的林下光照强度估测方法
Yu et al. Estimating the land-surface temperature of pixels covered by clouds in MODIS products
Zhu et al. Retrieval of land surface temperature with topographic effect correction from Landsat 8 thermal infrared data in mountainous areas
CN107389617A (zh) 基于高分四号卫星的气溶胶光学厚度的反演方法及设备
Li et al. Land Surface Temperature Retrieval from Sentinel-3A SLSTR data: Comparison among Split-window, Dual-window, Three-channel and Dual-angle algorithms
Diallo-Dudek et al. Land surface temperature in the urban area of Lyon metropolis: A comparative study of remote sensing data and MesoNH model simulation
Parol et al. Review of capabilities of multi-angle and polarization cloud measurements from POLDER
Nauss et al. The intercomparison of selected cloud retrieval algorithms
Wang et al. Estimation of surface longwave radiation over the Tibetan plateau region using MODIS data for cloud-free skies
Al Fahdawi et al. Measurement Albedo Coefficient for Land Cover (LC) and Land Use (LU), using remote sensing techniques, a study case: Fallujah City
Yin et al. Estimate of ocean wind vectors inside tropical cyclones from polarimetric radiometer
Cao et al. Attitude measurement method of spinning projectile based on infrared focal plane array
Curtis et al. Gulf Stream thermal gradients from satellite, ship, and aircraft observations
Burrage et al. Deriving sea surface salinity and density variations from satellite and aircraft microwave radiometer measurements: application to coastal plumes using STARRS
Lyu et al. Solar-blind ultraviolet upwelling radiance diurnal variation led by observation geometry factors on geostationary attitude sensor limb viewing
Kelly et al. Assimilation of nontraditional datasets to improve atmospheric compensation
Cao et al. Optimizing the protocol of near-surface remote sensing experiments over heterogeneous canopy using DART simulated images
Lei et al. Simulation of future geostationary ocean color images
Peng et al. Estimating of the total atmospheric precipitable water vapor amount from the Chinese new generation polar orbit FengYun meteorological satellite (FY-3) data

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170329

Termination date: 20181223