第〇部分 前言
一幅带有地理空间参考、且包含海陆分布的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、豆包为本文方法、文案提供的参考和帮助