[QGIS]绘制和建模海平面上升

根据美国海洋大气管理局（NOAA）公布的数据显示，2018年全球平均海平面比1993年的平均值高出了81毫米。根据夏威夷大学Philip Thompson相关卫星数据分析显示，自1993年以来，海平面上升速率为平均每年约3.1毫米。

数据来源：

NOAA(https://www.ngdc.noaa.gov/mgg/coastal/coastal.html)

USGS(https://earthexplorer.usgs.gov/)

软件版本：

QGIS 3.10

要求：
1. 使用高分辨率数字高程模型执行海平面变化分析

2. 创建一张清晰直观地显示潜在淹没区域的地图(至少三个不同严重程度的海平面变化场景)

一、获取数字高程模型数据

数字高程模型数据可以从NOAA或者USGS网址获取，此处从NOAA(https://www.ngdc.noaa.gov/mgg/coastal/coastal.html)获取夏威夷的哈纳莱伊数字高程数据，在网页上搜索Hanalei，点击链接然后下载其中的NetCDF文件(.nc).

二、导入图层

点击Layer->Data Source Manager，在Raster标签页的Source栏选择所下载的NetCDF文件(.nc)，添加沿海高程模型图层.

导入栅格图层

三、栅格数据单位转换为英尺

点击Raster->Raster Calculator...打开栅格计算界面，表达式为"hanalei_13_mhw_2011@1"*3.28084（1米=3.28084英尺），注意点击输出文件名称、格式和坐标系。点击ok，将会生成以英尺为单位的新的沿海高程模型图层(hanalei_ft)。

栅格数据转换

四、采用栅格计算的方法生成hanalei_0ft、hanalei_3ft、hanalei_5ft和hanalei_10ft四个图层

这四个图层将海平面高度区分为两类（如hanalei_3ft图层分为海平面高度3英尺以上和海平面高度3英尺及其以下），它们的栅格计算表达式分别为："hanalei_ft@1">0 、"hanalei_ft@1" > 3 、"hanalei_ft@1">5 、"hanalei_ft@1" > 10

五、采用栅格计算的方法生成hanalei_3ft_change、hanalei_5ft_change和hanalei_10ft_change三个图层

这三个图层将海平面高度区分为三类（hanalei_3ft_change图层分为海平面以下、海平面0-3英尺和海平面3英尺以上），它们的栅格计算表达式分别为："hanalei_0ft@1"+"hanalei_3ft@1" 、"hanalei_0ft@1"+"hanalei_5ft@1" 和 "hanalei_0ft@1"+"hanalei_10ft@1" 。

六、安装QuickMapServices插件

点击Plugins->Manage and Install Plugins打开插件界面，搜素QuickMapServices并安装该插件。安装后可以在Web菜单找到QuickMapServices菜单项，如果其中可供选择的底图不够多，可以点击Settings获取更多的底图服务：

底图服务设置入口

获取更多的底图服务

七、添加底图和顶图

分别点击QuickMapServices的ESRI Satellite和ESRI Boundaries&Places，将ESRI Satellite作为底图，ESRI Boundaries&Places作为顶层视图

八、图层着色

hanalei_3ft_change、hanalei_5ft_change和hanalei_10ft_change三个图层的渲染类型均选择为Paleted/Unique values,由于我们仅关心变化，可以将0和2的值的颜色透明度调为0

着色配置

九、生成打印图层并加上标题和图例

Project->New print Layout...生成打印图层，并添加地图、标题和图例。底图透明度调为50，这样可以将海平面上升后的颜色调为蓝色。

hanalei简易效果图

Savannah效果图