昼夜遥感影像短波红外火山温度-发射率联合反演方法

摘要

本发明涉及一种昼夜遥感影像短波红外火山温度‑发射率联合反演方法，包括Landsat8遥感影像数据获取，日间遥感影像预处理，夜间遥感影像预处理，确定火山区域，建立感兴趣区，两景影像对应相同区域，大气上界太阳辐照度计算，建立昼‑夜火山温度联合反演方程以及温度‑发射率联合反演求解。本发明通过基于Landsat8短波红外昼夜遥感影像像元尺度火山温度‑发射率联合反演方法，获得像元尺度的火山温度及对应发射率，解决了数学算法求解缺乏物理意义、热红外影像存在混合像元以及点尺度与像元尺度对应关系未知的问题；在不需要实地观测的条件下，通过多时相遥感影像即可获取火山温度及发射率。

权利要求

1.一种昼夜遥感影像短波红外火山温度‑发射率联合反演方法，其特征在于，包括以下步骤：A、Landsat8遥感影像数据获取：选择相近时相的昼、夜影像数据；

B、日间遥感影像预处理：对Landsat8影像进行辐射定标、大气校正处理；

C、夜间遥感影像预处理：对Landsat8影像进行辐射定标后，计算大气透过率，并通过Landsat8波段响应函数得到各波段大气透过率，其中，大气透过率参数由火山地区经纬度及影像时相确定，夜间影像数据大气校正处理为：其中，LN为大气校正后夜间影像波段λ的辐射亮度，LN0为辐射定标后夜间影像波段λ的辐射亮度，τ为波段λ的夜间大气透过率；

D、确定火山区域，建立感兴趣区，两景影像对应相同区域；

E、大气上界太阳辐照度E0的计算：大气上界太阳辐照度E0由Landsat8影像数据头文件中参量系数间接获取：其中，Lmax为探测器可接收到的最大辐射亮度，Lmin为探测器可接收到的最小辐射亮度，ρmax为探测器可接收到的最大反射率，ρmin为探测器可接收到的最小反射率；

F、建立昼‑夜火山温度联合反演方程：利用昼‑夜两幅影像对应像元的6，7波段分别建立温度反演模型，其中各波段发射率不同但不受温度变化影响，即日、夜间方程中对应波段发射率相同，受大气条件影响，大气校正存在误差，因此在日、夜间影像各波段间温度反演存在差异，故反演温度不为相同值，以此建立联合反演方程：其中， 为中心波长为λi波段的日间影像辐射亮度， 为中心波长为λi波段的夜间影像辐射亮度， 为中心波长为λi波段的发射率， 为中心波长为λi波段的日间温度，‑34为中心波长为λi波段的夜间温度，h为普朗克函数，取值为6.626×10 J·s，k为玻尔‑23 ‑8兹曼常数，取值为1.3806×10 J/K，c为光速，取值为2.998×10 m/s；

G、温度‑发射率联合反演求解：将各波段日‑夜间温度‑发射率联合反演方程，转化为温度反演方程，假设不同波段间温度反演差异最小，且日、夜间两幅影像反演温度差别最小，并建立日、夜间各波段温度反演差值最小约束条件，进而采用遗传算法求解发射率、温度最优解：其中，T0、T1分别为日、夜间各波段温度反演均值。

2.根据权利要求1所述的一种昼夜遥感影像短波红外火山温度‑发射率联合反演方法，其特征在于：步骤B，采用FLAASH大气校正方法进行大气校正。

3.根据权利要求1所述的一种昼夜遥感影像短波红外火山温度‑发射率联合反演方法，其特征在于：步骤C，采用MORTRAN4.0大气辐射传输模型计算大气透过率。

说明书

昼夜遥感影像短波红外火山温度‑发射率联合反演方法

技术领域

[0001] 本发明属于定量遥感观测技术领域，具体涉及一种具体涉及一种遥感影像火山温度‑发射率联合反演方法，特别涉及一种基于Landsat8短波红外昼夜遥感影像火山温度‑发射率联合反演方法。

背景技术

[0002] 火山喷发活动具有突发性和剧烈性特征，对地表环境及人类生命财产安全带来巨大损失，利用遥感卫星实现火山温度反演是定量遥感中的重要议题。火山温度的反演对于对火山识别及火山活动动态监测具有现实意义。

[0003] 由于火山活动中难以进行实地勘察，遥感影像成为其活动监测的有效手段。目前，火山监测主要采用NOAA/AVHRR和MODIS影像或成品数据，通常采用中红外或热红外波段进行火山目标的识别。温度与发射率分离算法(Temperature and emissivity separate,TES)在红外波段已有广泛研究。从算法角度来看，基于ADE、牛顿迭代等方法进行温度‑发射率分离解算，其依据经验统计或数学方法进行求解，难以得到具有明确遥感物理意义的精确解。从数据分辨率角度来看，虽然NOAA/AVHRR、MODIS影像等数据时间分辨率较高，但其空间分辨率较低(1km、500m，250m)，导致反演结果通常为混合像元尺度的温度与发射率。尤其对于小型火山，难以准确圈定火山温度异常区域范围大小和温度范围。

[0004] 已有研究表明短波红外是火山等高温目标探测的敏感波段，而在短波红外波段高温目标温度反演研究还较少研究。Landsat8短波红外数据相较NOAA/AVHRR、MODIS等热红外分辨率较高，有学者利用Landsat短波红外数据对火山熔岩流、岩浆湖、火山湖进行热异常识别或温度反演，其中发射率作为反演参数，有学者采用实地点观测值或实验观测经验值进行替代，但存在数学算法求解缺乏物理意义、热红外影像存在混合像元以及点尺度与像元尺度对应关系未知的问题，故点尺度与像元尺度的对应关系依然是有待研究的科学问题。

发明内容

[0005] 本发明的目的就在于提供一种基于Landsat8短波红外昼夜遥感影像像元尺度火山温度‑发射率联合反演方法，以解决数学算法求解缺乏物理意义、热红外影像存在混合像元以及点尺度与像元尺度对应关系未知的问题，通过建立多波段昼‑夜火山温度联合反演方程，以昼夜两影像中对应像元温度差最小作为约束条件，对预处理后的昼夜遥感影像中火山纯像元进行火山温度及发射率反演，获得像元尺度的火山温度及对应发射率，在不需要实地观测的条件下，通过多时相遥感影像即可获取火山温度及发射率。

[0006] 本发明的目的是通过以下技术方案实现的：

[0007] 一种昼夜遥感影像短波红外火山温度‑发射率联合反演方法，包括以下步骤：

[0008] A、Landsat8遥感影像数据获取：选择相近时相的昼、夜影像数据；

[0009] B、日间遥感影像预处理：对Landsat8影像进行辐射定标、大气校正处理；

[0010] C、夜间遥感影像预处理：对Landsat8影像进行辐射定标后，计算大气透过率，并通过Landsat8波段响应函数得到各波段大气透过率，其中，大气透过率参数由火山地区经纬度及影像时相确定，夜间影像数据大气校正可处理为：

[0011]

[0012] 其中，LN为大气校正后夜间影像波段λ的辐射亮度，LN0为辐射定标后夜间影像波段λ的辐射亮度，τ为波段λ的夜间大气透过率；

[0013] D、确定火山区域，建立感兴趣区，两景影像对应相同区域；

[0014] E、大气上界太阳辐照度E0的计算：大气上界太阳辐照度E0可由Landsat8影像数据头文件中参量系数间接获取：

[0015]

[0016] 其中，Lmax为探测器可接收到的最大辐射亮度，Lmin为探测器可接收到的最小辐射亮度，ρmax为探测器可接收到的最大反射率，ρmin为探测器可接收到的最小反射率；

[0017] F、建立昼‑夜火山温度联合反演方程：利用昼‑夜两幅影像对应像元的6，7波段分别建立温度反演模型，其中各波段发射率不同但不受温度变化影响，即日、夜间方程中对应波段发射率相同，受大气条件等影响，大气校正存在误差，因此在日、夜间影像各波段间温度反演存在差异，故反演温度不为相同值，以此建立联合反演方程：

[0018]

[0019] 其中， 为中心波长为λi波段的日间影像辐射亮度， 为中心波长为λi波段的夜间影像辐射亮度， 为中心波长为λi波段的发射率， 为中心波长为λi波段的日间‑34温度， 为中心波长为λi波段的夜间温度，h为普朗克函数，取值为6.626×10 J·s，k为‑23 ‑8
玻尔兹曼常数，取值为1.3806×10 J/K，c为光速，取值为2.998×10 m/s；

[0020] G、温度‑发射率联合反演求解：将各波段日‑夜间温度‑发射率联合反演方程，转化为温度反演方程，假设不同波段间温度反演差异最小，且日、夜间两幅影像反演温度差别最小，并建立日、夜间各波段温度反演差值最小约束条件，进而采用遗传算法求解发射率、温度最优解：

[0021]

[0022] 其中，T0、T1分别为日、夜间n个波段温度反演均值。

[0023] 进一步地，步骤B，采用FLAASH大气校正方法进行大气校正。

[0024] 进一步地，步骤C，采用MORTRAN4.0大气辐射传输模型计算大气透过率。

[0025] 与现有技术相比，本发明的有益效果是：

[0026] 本发明通过基于Landsat8短波红外昼夜遥感影像像元尺度火山温度‑发射率联合反演方法，获得像元尺度的火山温度及对应发射率，解决了数学算法求解缺乏物理意义、热红外影像存在混合像元，避免了点尺度与像元尺度对应关系未知的问题；在不需要实地观测的条件下，通过多时相遥感影像即可获取火山温度及发射率。

附图说明

[0027] 为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。

[0028] 图1昼夜遥感影像短波红外火山温度‑发射率联合反演方法流程图；

[0029] 图2本申请实施例中基拉韦厄火山2014年4月24日Landsat8日间遥感影像图；

[0030] 图3本申请实施例中基拉韦厄火山2014年4月26日Landsat8夜间遥感影像图；

[0031] 图4本申请实施例中基于MORTRAN4.0大气透过率的计算结果；

[0032] 图5本申请实施例中火山感兴趣区，图5a日间火山感兴趣区，图5b夜间火山感兴趣区；

[0033] 图6本申请实施例中火山感兴趣区发射率反演结果，图6a第五波段部分像元发射率结果，图6b第六波段部分像元发射率结果，图6c第七波段部分像元发射率结果；

[0034] 图7本申请实施例中火山感兴趣区温度反演结果。

具体实施方式

[0035] 下面结合具体实施案例对本发明作进一步的说明，本实施案例在以本发明技术为前提下进行实施，给出了详细的实施方式，但本发明的保护范围不限于以下的实施例。

[0036] 一种昼夜遥感影像短波红外火山温度‑发射率联合反演方法，包括以下步骤：

[0037] A、Landsat8遥感影像数据获取：选择相近时相的昼、夜影像数据；

[0038] B、日间遥感影像预处理：对Landsat8影像进行辐射定标，采用FLAASH大气校正方法进行大气校正处理；

[0039] C、夜间遥感影像预处理：对Landsat8影像进行辐射定标后，采用MORTRAN4.0大气辐射传输模型计算大气透过率，并通过Landsat8波段响应函数得到各波段大气透过率，其中，大气透过率参数由火山地区经纬度及影像时相确定，由于月球对地辐照度与火山目标辐射亮度相比可以忽略，因此夜间影像数据大气校正可处理为：

[0040]

[0041] 其中，LN为大气校正后夜间影像辐射亮度，LN0为辐射定标后夜间影像辐射亮度，τ为夜间大气透过率；

[0042] D、确定火山区域，建立感兴趣区，两景影像对应相同区域；

[0043] E、大气上界太阳辐照度计算：大气上界太阳辐照度可由Landsat8影像数据头文件中参量系数间接获取：

[0044]

[0045] 其中，Lmax为探测器可接收到的最大辐射亮度，Lmin为探测器可接收到的最小辐射亮度，ρmax为探测器可接收到的最大反射率，ρmin为探测器可接收到的最小反射率；

[0046] F、建立昼‑夜火山温度联合反演方程：利用昼‑夜两幅影像对应像元的6，7波段分别建立温度反演模型，其中各波段发射率不同但不受温度变化影响，即日、夜间方程中对应波段发射率相同，受大气条件等影响，大气校正存在误差，因此在日、夜间影像各波段间温度反演存在差异，故反演温度不为相同值，以此建立联合反演方程：

[0047]

[0048] 其中， 为中心波长为λi波段的日间影像辐射亮度， 为中心波长为λi波段的夜间影像辐射亮度， 为中心波长为λi波段的发射率， 为中心波长为λi波段的日间‑34温度， 为中心波长为λi波段的夜间温度，h为普朗克函数，取值为6.626×10 J·s，k为‑23 ‑8
玻尔兹曼常数，取值为1.3806×10 J/K，c为光速，取值为2.998×10 m/s；

[0049] G、温度‑发射率联合反演求解。将各波段日‑夜间温度‑发射率联合反演方程，转化为温度反演方程，假设不同波段间温度反演差异最小，且日、夜间两幅影像反演温度差别最小，并建立日、夜间各波段温度反演差值最小约束条件，进而采用遗传算法求解发射率、温度最优解：

[0050]

[0051] 其中，T0、T1分别为日、夜间n个波段反演温度均值。

[0052] 实施例

[0053] 本申请实施例按照图1所示流程图进行昼夜遥感影像短波红外火山温度‑发射率联合反演。

[0054] 如图1所示，一种昼夜遥感影像短波红外火山温度‑发射率联合反演方法，包括以下步骤：

[0055] A、获取Landsat8遥感影像数据。Landsat8影像选择相近时相的昼、夜影像数据。图2为基拉韦厄火山2014年4月24日Landsat8日间遥感影像图，图3表示基拉韦厄火山2014年4月26日Landsat8夜间遥感影像图。

[0056] B、对日间遥感影像进行预处理。包括Landsat8影像辐射定标和大气校正，采用FLAASH大气校正方法进行大气校正。

[0057] C、对夜间遥感影像进行预处理。对Landsat8影像进行辐射定标后，采用MORTRAN4.0大气辐射传输模型计算大气透过率，并通过Landsat8波段响应函数得到各波段大气透过率。其中大气透过率参数由火山地区经纬度及影像时相确定。图4为基于MORTRAN4.0大气透过率的计算结果。

[0058] D、确定火山区域，建立感兴趣区。图5(a)为日间影像感兴趣区范围，图5(b)为夜间影像感兴趣区范围。

[0059] E、大气上界太阳辐照度计算。大气上界太阳辐照度可由Landsat8影像数据头文件中参量系数间接获取：

[0060]

[0061] 分别计算得到近红外‑短波红外(第5、6、7波段)的大气上界太阳辐照度。

[0062] F、建立昼‑夜火山温度‑发射率联合反演方程。对于7波段未饱和的波段，采用6，7波段联立温度‑发射率反演方程，对于7波段饱和像元，采用5，6波段联立温度‑发射率反演方程，对于6，7波段饱和像元，采用5波段进行温度‑发射率联合反演。方程建立思想为利用昼‑夜两幅影像对应像元的6，7波段分别建立温度反演模型：

[0063]

[0064] G、温度‑发射率联合反演求解。将各波段日‑夜间温度‑发射率联合反演方程，转化为温度反演方程，假设不同波段间温度反演差异最小，且日、夜间两幅影像反演温度差别最小，并建立日、夜间各波段温度反演差值最小约束条件，进而采用遗传算法求解发射率、温度最优解：

[0065]

[0066] 其中，T0、T1分别为日、夜间n个波段反演温度均值；h为普朗克函数，取值为6.626×‑34 ‑23 ‑810 J·s；k为玻尔兹曼常数，取值为1.3806×10 J/K；c为光速，取值为2.998×10 m/s。反演发射率及温度结果由图6、图7所示。

[0067] 注意，上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解，本发明不限于这里所述的特定实施例，对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此，虽然通过以上实施例对本发明进行了较为详细的说明，但是本发明不仅仅限于以上实施例，在不脱离本发明构思的情况下，还可以包括更多其他等效实施例，而本发明的范围由所附的权利要求范围决定。