根據(jù)美國海洋大氣管理局(NOAA)公布的數(shù)據(jù)顯示,2018年全球平均海平面比1993年的平均值高出了81毫米。根據(jù)夏威夷大學(xué)Philip Thompson相關(guān)衛(wèi)星數(shù)據(jù)分析顯示,自1993年以來,海平面上升速率為平均每年約3.1毫米。
數(shù)據(jù)來源:
NOAA(https://www.ngdc.noaa.gov/mgg/coastal/coastal.html)
USGS(https://earthexplorer.usgs.gov/)
軟件版本:
QGIS 3.10
要求:
1. 使用高分辨率數(shù)字高程模型執(zhí)行海平面變化分析
2. 創(chuàng)建一張清晰直觀地顯示潛在淹沒區(qū)域的地圖(至少三個不同嚴重程度的海平面變化場景)
一、獲取數(shù)字高程模型數(shù)據(jù)
數(shù)字高程模型數(shù)據(jù)可以從NOAA或者USGS網(wǎng)址獲取,此處從NOAA(https://www.ngdc.noaa.gov/mgg/coastal/coastal.html)獲取夏威夷的哈納萊伊數(shù)字高程數(shù)據(jù),在網(wǎng)頁上搜索Hanalei,點擊鏈接然后下載其中的NetCDF文件(.nc).
二、導(dǎo)入圖層
點擊Layer->Data Source Manager,在Raster標(biāo)簽頁的Source欄選擇所下載的NetCDF文件(.nc),添加沿海高程模型圖層.

三、柵格數(shù)據(jù)單位轉(zhuǎn)換為英尺
點擊Raster->Raster Calculator...打開柵格計算界面,表達式為"hanalei_13_mhw_2011@1"*3.28084(1米=3.28084英尺),注意點擊輸出文件名稱、格式和坐標(biāo)系。點擊ok,將會生成以英尺為單位的新的沿海高程模型圖層(hanalei_ft)。

四、采用柵格計算的方法生成hanalei_0ft、hanalei_3ft、hanalei_5ft和hanalei_10ft四個圖層
這四個圖層將海平面高度區(qū)分為兩類(如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三個圖層
這三個圖層將海平面高度區(qū)分為三類(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獲取更多的底圖服務(wù):


七、添加底圖和頂圖
分別點擊QuickMapServices的ESRI Satellite和ESRI Boundaries&Places,將ESRI Satellite作為底圖,ESRI Boundaries&Places作為頂層視圖
八、圖層著色
hanalei_3ft_change、hanalei_5ft_change和hanalei_10ft_change三個圖層的渲染類型均選擇為Paleted/Unique values,由于我們僅關(guān)心變化,可以將0和2的值的顏色透明度調(diào)為0

九、生成打印圖層并加上標(biāo)題和圖例
Project->New print Layout...生成打印圖層,并添加地圖、標(biāo)題和圖例。底圖透明度調(diào)為50,這樣可以將海平面上升后的顏色調(diào)為藍色。

