第〇部分 前言

一幅帶有地理空間參考、且包含海陸分布的DEM圖像在研究區的繪制中非常常見，本文將實現以下圖像的繪制

關鍵步驟：

（1）NOAA-NCEI官方下載最新的ETOPO-2022 DEM數據

（2）在ArcGIS（ArcMap）中使用“柵格計算器”分別獲取海洋、陸地圖層，并分別施加色帶

第一部分 ETOPO-2022 DEM數據下載

ETOPO是一個全球高分辨率的地形和海洋深度數據集，提供了地表的詳細地形和海洋深度。特點是“全球覆蓋”、“高分辨率”。當前最新版本是“ETOPO 2022”

ETOPO-2022 DEM數據下載網址：https://www.ncei.noaa.gov/products/etopo-global-relief-model

對于“ETOPO 2022”的DEM數據，官網提供了“GEOTIFF”和“NC”兩種格式的下載。本文以GEOTIFF為例，GEOTIFF分成了三種分辨率：15弧秒分辨率（部分地區）、30弧秒分辨率（全球）、60弧秒分辨率（全球）。弧秒數越大，空間分辨率越低，數據量越小。本文以60弧秒的全球范圍為例。

第二部分 在ArcGIS的ArcMap中操作

（1）將下載好的原始圖層導入ArcGIS新建的工程中（請無視ETOPO_2022_OCEAN.tif 和 ETOPO_2022_LAND.tif 兩個圖層）

（2）在作者已知的知識范圍內，沒有一種方法能夠實現同一張TIF圖像使用兩個不同的色帶（colorbar）。因此便嘗試以“0”值作為分界，分別基于原始TIF生成兩個圖層，一個是高程高于0的陸地圖層，一個是高程低于等于0的海洋圖層。這樣、兩個圖層可以分別擁有不同的色帶

以上想法基于“柵格計算器”實現：ArcToolBox—空間分析工具（Spatial Analysis Tool）—地圖代數（Map Algebra）—柵格計算器（Raster Caculator）（圖中“空間的（Spatial）”拼寫錯誤請忽略）

（3.1）先來獲取陸地的圖層：

表達式中輸入：

Con("ETOPO_2022_v1_60s_N90W180_bed.tif"< 0,SetNull("ETOPO_2022_v1_60s_N90W180_bed.tif"<0,"ETOPO_2022_v1_60s_N90W180_bed.tif"),"ETOPO_2022_v1_60s_N90W180_bed.tif")

輸出路徑選擇合適的位置，盡量保證路徑上不存在中文：

C:\Users\--o\Desktop\新建文件夾\ETOPO_2022_LAND.tif

?（3.2）然后獲取海洋的圖層：

表達式中輸入：

Con("ETOPO_2022_v1_60s_N90W180_bed.tif" >= 0, SetNull("ETOPO_2022_v1_60s_N90W180_bed.tif" >= 0, "ETOPO_2022_v1_60s_N90W180_bed.tif"), "ETOPO_2022_v1_60s_N90W180_bed.tif")

輸出路徑選擇合適的位置，盡量保證路徑上不存在中文：

C:\Users\--o\Desktop\新建文件夾\ETOPO_2022_OCEAN.tif

（4）如此一來，就分別得到了陸地圖層（ETOPO_2022_LAND.tif ）和海洋圖層（ETOPO_2022_OCEAN）：

（5）之后，分別對陸地圖層和海洋圖層選擇合適的色帶即可，此處陸地圖層應用了“紅-綠漸變（反向）”、海洋圖層應用了“藍色漸變（反向）”。此時已經可以實現和前言及thumbnail相同的結果

參考內容：

[1] https://zhuanlan.zhihu.com/p/628414578（ETOPO-2022 DEM數據下載）

鳴謝：

感謝人工智能助手kimi、豆包為本文方法、文案提供的參考和幫助