[发明专利]一种基于夜间灯光数据提取地震烈度的方法

摘要

本发明提出的一种基于夜间灯光数据提取地震烈度的方法，属于夜间灯光遥感数据领域。本方法首先利用下载得到的全球夜间灯光数据裁剪出研究区地震前一年份和研究区地震年份的夜间灯光数据；提取地震前一年份和地震年份的夜间灯光数据差值，生成研究区各个区域的灯光亮度变化值的栅格数据；将栅格数据转换为点状矢量数据后，剔除0值，并进行空间插值，生成研究区灯光亮度变化值的栅格数据；然后对研究区灯光亮度变化值的栅格数据重分类，将栅格数据转换为面状矢量数据，并进行烈度赋值，进一步平滑处理后，最终生成地震烈度分布的面状矢量数据。本发明通过夜间灯光数据，准确提取地震烈度的分布情况，及时有效反映地震破坏情况。

主权项

一种基于夜间灯光数据提取地震烈度的方法，其特征在于，包括以下步骤：1)裁剪夜间灯光数据；从夜间灯光数据官方网站下载研究区的地震前一年份和地震年份的全球夜间灯光数据，根据研究区的矢量数据确定裁剪的区域,利用ArcGIS软件对下载得到的数据进行裁剪处理，生成裁剪后的研究区地震前一年份和研究区地震年份的夜间灯光数据；所述研究区为需要提取地震烈度的区域；2)提取研究区地震前一年份和研究区地震年份的夜间灯光数据差值；通过ArcGIS软件中ArcToolbox工具箱的地图代数工具对步骤1)得到的裁剪后的研究区地震前一年份和研究区地震年份的夜间灯光数据进行差值计算，生成研究区各个区域的灯光亮度变化值的栅格数据；3)利用ArcGIS软件将步骤2)得到的栅格数据，生成研究区的各个区域的灯光亮度变化值的点状矢量数据；4)剔除步骤3)生成的研究区的各个区域的灯光亮度变化值的点状矢量数据中的0值；5)对剔除0值后的研究区的的灯光亮度变化值的点状矢量数据进行空间插值，得到空间插值后的研究区灯光亮度变化值的栅格数据；6)根据灯光亮度变化程度，将空间插值后的研究区灯光亮度变化值的栅格数据进行重分类，生成重分类后的研究区灯光亮度变化值的栅格数据；所述灯光亮度变化程度分为1‑9共9个级别；7)通过ArcGIS软件将步骤6)得到的重分类后的研究区灯光亮度变化值的栅格数据转换为面状矢量数据；8)对步骤7)得到的研究区灯光亮度变化值的面状矢量数据进行地震烈度赋值，地震烈度随灯光亮度变化值的变化而变化：令灯光亮度变化值7及以上的为XI度，4‑6为X度，1‑3为IX度，并赋予不同的烈度不同的颜色以示区别；9)利用ArcGIS软件对步骤7)得到的面状矢量数据进行平滑处理，生成平滑后的地震烈度分布的面状矢量数据。