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