[发明专利]一种基于3S技术的森林火灾监测方法

摘要

本发明属森林火灾监测方法技术领域，尤其涉及结合中分辨率遥感数据和多光谱高分辨率遥感影像、GPS定位采样信息对森林火灾进行监测，能够提高火灾信息识别的准确性和效率的一种基于3S技术的森林火灾监测方法，包括受灾区域的中低分辨率MODIS数据和灾前以及灾后多光谱高分辨率哨兵2号遥感监测卫星影像并进行处理，外业采集，火灾区域的识别和生成监测服务报告的步骤。本发明提供的一种外业采样提高遥感解译效率的方法具有能够使技术人员能通过实时获取影像数据，用中等分辨率的影像快速进行火灾区域粗定位，结合GPS定位采样信息，最终准确快速的优点。

主权项

1.一种基于3S技术的森林火灾监测方法，其特征在于，包括以下步骤：步骤1）、当收到火灾初步信息后，立即在NASA官网和欧空局数据分发网站上下载灾后的包括受灾区域的中低分辨率MODIS数据和灾前以及灾后多光谱高分辨率哨兵2号遥感监测卫星影像并进行处理，为后期的火灾识别做数据准备；1.1 MODIS数据的处理步骤：a.预处理，包括影像拼接、重投影、格式的转换、辐射定标和大气校正；所述影像拼接、重投影、格式的转换为利用针对MODIS数据的处理工具MRT完成MODIS数据的影像拼接、重投影和格式的转换，所述MRT是一种针对MODIS数据的处理工具；所述辐射定标为将记录的原始DN值转换为大气外层表面反射率；所述大气校正为将辐射亮度或者表面反射率转换为地表实际反射率；b.数据反演：根据近红外和红外波段反演出归一化植被指数NDVI，反演归一化公式为，其中：：红光反射率，：近红外反射率；然后计算出热红外波段的辐射率，辐射率计算为，其中：：波段增益，取值：0.00084002；：辐射亮度，单位：；：波段偏移量，取值：1577.33972168；：灰度值，根据黑体辐射定律反变换，计算得到地表的亮温值，地表的亮温值计算公式为：其中：T：亮度温度，单位：K，h：普朗克常数，k：波尔兹曼常数，C：光速，：波长，L：辐射亮度，单位：；1.2对哨兵2号遥感监测卫星的L1C级多光谱数据MSI进行辐射定标和大气校正处理；具体步骤为：a.利用ESA发布的用于生产L2A级数据的插件Sen2cor进行辐射定标和大气校正的预处理，其中对哨兵2号遥感监测卫星原始数据的10m和20m空间分辨率波段进行处理；b.数据重采样：将L2A级数据中的B5‑Vegetation Red Edge、B6‑Vegetation Red Edge、B7‑Vegetation Red Edge、B8A‑Vegetation Red Edge、B11‑SWIR、B12‑SWIR六个20m分辨率的波段利用Cubic Convolution的方法进行重采样为10m分辨率；c.波段合成：将L2A级数据中的B2‑Blue、B3‑Green、B4‑Red、B8‑NIR，和经过数据重采样得到的B5‑Vegetation Red Edge、B6‑Vegetation Red Edge、B7‑Vegetation Red Edge、B8A‑Vegetation Red Edge、B11‑SWIR、B12‑SWIR，共计十个10m分辨率的影像波段合成为一景多光谱高分辨率影像；步骤2）、外业采集：派出外业调查人员到受灾现场使用RTK进行精准定位采样，以及利用拍照设备和使用移动端的绘制软件进行外业采样，对受灾现场进行拍照和地块信息的采集，从而获得受灾地点的GPS坐标、受灾区域情况信息等，同时结合土地利用现状图对拍摄的受灾地区记录其他火灾基本信息，其他火灾基本信息包括受灾时间和受灾情况；步骤3）、火灾区域的识别：3.1利用MODIS数据，根据反演的NDVI值和亮温值的相应阈值识别出火灾区域的大概位置；3.2根据得到的火灾大概位置对火灾前后的哨兵2号遥感监测卫星数据进行裁剪，缩小数据的大小和范围3.3基于预处理后的哨兵2号影像数据运用最小噪声分离变换的方法MNF对数据进行降维处理后，利用PPI纯净像元指数通过反复迭代图像中的像素点，从而找到图像中相对纯净的像元作为图像波谱的基本组分，确定终端像元的范围，将终端像元运用混合调谐匹配滤波MTMF混合像元分解，得到终端像元类型的丰度图，从丰度图上提取不同组成比例的像元，从而降低混合像元的影响；3.4基于以上处理的哨兵2号遥感监测卫星数据，根据上述确定的火灾和非火灾的纯净端元反射率的阈值，利用决策树分类建立模型，依据阈值的变化范围进行火灾区域的识别；3.5结合实测样方火灾的采样信息，验证火灾区域的识别精度；步骤4）、生成监测服务报告：对受灾的情况以专题地图的形式进行展示，形成完整的监测服务报告。