之前工作中处理过风暴潮的数据，也获取到了不同等级台风水淹的DEM数据，不过还是很好奇是怎么推演出来的，最近一段时间看QGIS比较多，加上看到了一篇文章《GIS软件进行风暴潮沿海洪水建模》 ，于是简单尝试了一下，我也东施效颦了一把。
主要目的是看一下海口市在风暴潮1米、2米、5米、10米情况下的水淹情况，当然这个肯定不精准，未考虑到承灾体、城市建设等情况。
前置条件：还是承接上文，把海口高程地图从海南的shp文件中提取出来后，并叠加海口shp文件，这个是基础。
DEM（Digital Elevation Model）数据格式是用于描述地球表面高度信息的数字地形数据格式。DEM格式最初是由美国联邦地理局（USGS）推出的，是包括高程、坡度、坡向等参数的数字高程模型，广泛应用于水文、地形分析、遥感图像处理、地质勘探和三维可视化等领域。DEM数据格式通常是由等高线制作而成的，它通过一定精度采集和计算，生成高程节点数据，然后根据节点数据生成高程模型网格，以展示出地表上各点的高程值。
第一步，光栅计算 Raster Calculator
分别修改公式中的阈值 "Clipped (mask)@1" <= 1，保存1米、2米、5米、10米高程数据
第二步、定义 NoData 值
Raster >Covertion>Translate（convert format），将NoData 值必须指定为 0
第三步、矢量化处理
为分类图层定义 NoData 值后，将使用 Polygonize 工具将文件转换为矢量文件。通过 Raster >Covertion> > Polygonize（raster to vector）打开工具后，定义输入和输出文件。
第四步、过滤内陆水域
为了过滤满足高程阈值但不与海洋连接的内陆水域区域，应用 QGIS 中的按位置选择工具（vector > Research Tools > Select by location），但其实效果看不大出来
第五步，保存为shp文件
在layer处点击右键，Export ->Save Features As..
不过导出过程中也有一些问题存在，待后面仔细检查。
最后输出的水淹情况如下：
最后欢迎关注公众号：python与大数据分析
