一种利用Sentinel-2日间影像探测火炬的方法技术

本发明专利技术公开一种利用Sentinel

【技术实现步骤摘要】
一种利用Sentinel
‑
2日间影像探测火炬的方法


[0001]本专利技术涉及一种利用Sentinel
‑
2日间影像探测火炬的方法，属于遥感地学应用


技术介绍

[0002]在石油生产过程中，会产生伴生天然气，因为天然气大量堆积会造成的安全隐患以及天然气回收利用过高的经济成本，往往会采用燃烧的方式处理伴生天然气，燃除的天然气被称为火炬。火炬会排放大量的二氧化碳和黑炭加剧全球变暖，还会排放SO2、NOX、PM、CO，等污染气体严重影响燃烧区域周边区域的空气质量，使长期暴露在燃除区域工作的工人和居民产生呼吸道疾病，目前，已有80个政府和石油天然气公司做出了“零常规燃除”的承诺，但是由于天然气燃除行为通常位于偏远地区，很少有国家会有详细的燃除地点和燃除量的公开报告，因此，需要采用遥感探测技术来跟踪全球的天然气燃除情况。
[0003]由于夜间没有来自背景的太阳辐射的干扰，早期大部分研究使用的都是夜间传感器来探测火炬，但是夜间传感器的空间分辨率往往都较低，一些燃烧比较微弱的火炬常常无法被探测出来，或者距离较近的火炬无法被区分开。为了解决这一问题，目前已经开始有研究在日间高空间分辨率影像上使用多波段阈值算法或热异常指数等进行火炬的探测，但是日间影像探测结果往往会受到太阳辐射的干扰，将一些高反射率地物如：建筑物，生物质燃烧等误认为是火炬，有研究采用全球居住建成数据对建筑物进行剔除，但是该数据集具有一定的误差，无法完全剔除建筑物的干扰，有时还可能将邻近建成区的火炬误判为建筑物进行剔除。

技术实现思路

[0004]鉴于上述问题，本专利技术提出一种首先使用日间影像进行火炬的初步探测，再根据潜在火炬的夜间温度进行筛选，在保证火炬的高精度的探测下，尽最大可能减少其他高反射率地物的干扰，以获取精准无误的火炬分布结果的利用Sentinel
‑
2日间影像探测火炬的方法。
[0005]技术方案：为了解决上述技术问题，本专利技术采用的技术方案如下：
[0006]一种利用Sentinel
‑
2日间影像探测火炬的方法，包括如下步骤：
[0007]步骤1：下载研究区域Sentinel
‑
2 MSI影像，再将各波段重采样为20m空间分辨率，预处理得到大气表观反射率图像；
[0008]步骤2：使用近红外（Narrow NIR）波段、短波红外1（SWIR1） 波段、短波红外2（SWIR2）波段，计算热异常指数（Thermal Anomaly Index, TAI）,并结合火炬样本确定提取火炬的TAI阈值；
[0009]步骤3：基于下载的全年的Sentinel
‑
2 影像，根据TAI阈值和一年被探测到的次数确定潜在火炬；
[0010]步骤4：下载研究区域NOAA
‑
20 VIIRS夜间影像，基于普朗克定律计算出潜在火炬
一年的夜间温度，并结合火炬样本得到火炬在燃烧月份出现温度高于1600K的平均次数作为在潜在火炬中筛选出真实火炬的阈值；
[0011]步骤5：对潜在火炬使用步骤4得到的阈值进行筛选，最终得到研究区域的火炬分布结果。
[0012]进一步地，步骤1中，下载的Sentinel
‑
2 MSI 影像为Level 1c的大气表观反射率数据，再使用SNAP软件将各个波段的空间分辨率都重采样为20m，之后在ENVI中用波段计算器将影像记录的DN值除以1000转换为真实的大气表观反射率。
[0013]进一步地，步骤2的具体实施方法为：
[0014]步骤21：首先使用处理好的Sentinel
‑
2的遥感影像的近红外（Narrow NIR）波段、短波红外1（SWIR1） 波段、短波红外2（SWIR2）波段，计算热异常指数（Thermal Anomaly Index, TAI），所述TAI指数计算公式如下：
[0015][0016]式中：, , 分别为Sentinel
‑
2 遥感影像的短波红外2波段（SWIR2）、短波红外1波段（SWIR1）和近红外波段（Narrow NIR）的大气表观反射率数据。
[0017]步骤22：再统计所有火炬样本在Sentinel
‑
2影像对应位置像元的TAI值，有90%以上的火炬点的TAI值是高于0.45的，所以将提取火炬的TAI阈值设定为0.45。
[0018]进一步地，步骤3中基于下载的全年的Sentinel
‑
2 影像，根据TAI阈值和一年被探测到的次数确定潜在火炬的具体判断条件是当影像中某一处像元位置的TAI计算值在一年内大于0.45的次数大于两次，或者一年中最大的TAI值超过了1，就将其判定为潜在火炬。
[0019]进一步地，步骤4的具体实施方法为：
[0020]步骤41：下载研究区域NOAA
‑
20 VIIRS夜间SDR影像，提取出M7、 M8、 M10这三个波段的夜间辐亮度数据，基于普朗克定律计算出潜在火炬一年的夜间温度，所述普朗克定律公式如下：
[0021][0022]式中： λ 是M7、M8、M10这三个波段的波长，值分别为0.87μm、1.24μm、1.6μm, 是影像中记录的M7、M8、M10这三个波段的夜间辐亮度，h是普朗克常数，值为6.6310
‑
34 Js， c是光速，值为3.010
8 m/s，K是玻尔兹曼常数，值为1.3810
‑
23 J/K，T 是燃烧火炬的绝对温度， 是燃烧火炬占据的像元比例,
ꢀꢀ
是燃烧火炬的比辐射率。
[0023]将存在潜在火炬的像元位置的M7、M8、M10波段的夜间辐亮度值和各个波段的波长作为普朗克定律拟合的输入条件，拟合程序的初始T为1000K，初始为1，约束条件为温度介于600K
‑
6000K，使用Matlab进行函数最优解的拟合，最终得到潜在火炬在一年的夜间拟合温度。
[0024]步骤42：汇总火炬样本的夜间拟合温度，当温度大于1600K时，认为火炬的状态为燃烧，统计火炬在燃烧月份出现燃烧状态的平均次数作为在潜在火炬中筛选出真实火炬的阈值，平均次数的计算公式为：
[0025][0026]式中：是火炬样本在一年燃烧月份中夜间温度大于1600K的总次数，是火炬样本在一年燃烧月份中夜间温度大于600K的总次数，31为一个月的总天数。
[0027]进一步地，步骤5中，对潜在火炬使用步骤4得到的阈值进行筛选，判定条件为：当潜在火炬在一年燃烧月份中至少有一个月的夜间温度大于1600K的天数大于N，则为真实火炬，或者当潜在火炬在一年燃烧月份中至少有两个月的夜间温度大于1600K的天数大于N/2，则为真实火炬，最终得到研究区域的火炬分布结果。
[0028]有益效果：与现有技术相比，本专利技术具有以下优点：
[0029]本发本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种利用Sentinel
‑
2日间影像探测火炬的方法，其特征在于：包括如下步骤：步骤1：获取研究区域Sentinel
‑
2 MSI影像，再将各波段重采样，然后预处理得到大气表观反射率图像；步骤2：使用近红外波段、短波红外1波段和短波红外2波段，计算热异常指数TAI，并结合火炬样本确定提取火炬的TAI阈值；步骤3：基于获取的全年的Sentinel
‑
2影像，根据TAI阈值和一年被探测到的次数确定潜在火炬；步骤4：下载研究区域NOAA
‑
20 VIIRS夜间影像，基于普朗克定律计算出潜在火炬一年的夜间温度，并结合火炬样本得到火炬在燃烧月份出现温度高于1600K的平均次数作为在潜在火炬中筛选出真实火炬的阈值；步骤5：对潜在火炬使用步骤4得到的阈值进行筛选，最终得到研究区域的火炬分布结果。2.根据权利要求1所述的利用Sentinel
‑
2日间影像探测火炬的方法，其特征在于：步骤1中，下载的Sentinel
‑
2 MSI 影像为Level 1c的大气表观反射率数据，再使用SNAP软件将各个波段的空间分辨率都重采样为20m，之后在ENVI中用波段计算器将影像记录的DN值除以1000转换为真实的大气表观反射率。3.根据权利要求1所述的利用Sentinel
‑
2日间影像探测火炬的方法，其特征在于：步骤2的具体实施方法为：步骤21：首先使用处理好的Sentinel
‑
2的遥感影像的近红外波段、短波红外1波段和短波红外2波段，计算热异常指数，热异常指数计算公式如下：；式中：, , 分别为Sentinel
‑
2遥感影像的短波红外2波段、短波红外1波段和近红外波段的大气表观反射率数据；步骤22：再统计所有火炬样本在Sentinel
‑
2影像对应位置像元的TAI值，有90%以上的火炬点的TAI值是高于0.45的，所以将提取火炬的TAI阈值设定为0.45。4.根据权利要求1所述的利用Sentinel
‑<...

【专利技术属性】
技术研发人员：张秀英，胡成林，高倩，
申请(专利权)人：南京大学，
类型：发明
国别省市：