「GIS」Cesium生成terrain地形数据-CTB方式及步骤

背景：项目前端使用Cesium，地形服务一直使用外网的，常常因为翻墙访问的问题，导致地形数据取不到，进而导致地球不能加载，故决定搭建自己的地形服务，彻底解决这个问题。

转载地址: liubf.com

背景：项目前端使用Cesium，地形服务一直使用外网的，常常因为翻墙访问的问题，导致地形数据取不到，进而导致地球不能加载，故决定搭建自己的地形服务，彻底解决这个问题。博文包含以下几个过程：

1. 下载原始地形数据，格式为.tif。
2. 处理地形数据，将零散的地形文件整合成一个地形文件。
3. 配置CTB环境
4. 使用ctb-tile指令将地形文件(.tif)加工成.terrain文件。
5. 发布地形服务并使用Cesium调用。

一、下载原始地形数据

• 下载地址(精度90m)：
  http://srtm.csi.cgiar.org/SELECTION/inputCoord.asp
• 如果你刚好需要的也是全国地形，可以从楼主网盘下载，更快更方便(密码zyt9)：https://pan.baidu.com/s/1YYC25b48QAVb64mW77vCEw

网盘包含两部分内容：

• dem.rar：原始数据，需要自行使用ArcMap工具整合。
• 合并后的数据：经过ArcMap工具处理的数据，已经整合成了一个地形文件，可以直接跳到第三步。

在网站下载流程：

1. 打开网页http://srtm.csi.cgiar.org/SELECTION/inputCoord.asp。如图:
2. 根据需要选择方便的选取方式：
   
   MutilpleSelection:点击哪个选哪个，可以选择多个，双击取消。
   EnableMouseDrag:框选，框选了什么范围就是什么范围。
   InputCoordinates:输入最西，最东。最南，最北的坐标来下载该范围内的所有数据。
3. 选择好范围后，点击“Click here to Begin Search”,进入下一步。
4. 可以看到刚刚选择了几个数据，点击每个数据的DownLoad，下载数据（这个地方比较繁琐，需要耐心）。
5. 下载完数据后解压所有的压缩包，并将tifwen文件拷贝到统一的文件夹内，这里的内容和我网盘中的dem.rar是一样的：

二、处理地形数据

由于CTB工具不支持DEM为NoData值和float的数据，所以需要对数据进行处理。

• 多个tif必须县合并镶嵌成一张tif
• pixeltype从float转为int
• NoData值处理为0

我使用的是ArcGis10.5，打开ArcMap:

打开 ArcToolbox->Data Management Tools->Raster->Raster Dataset ->Mosaic。

这个处理过程很久，可以在Result面板里查看转换过程。成功后，处理NoData值，处理公式为:

Con(IsNull("xxxx.tif"),0,"xxxx.tif")//xxxx.tif对应你的合并后的文件

三、配置CTB环境

1. 下载CTB工具包：资源
2. 将ctb和gdal-data拷贝至合适的目录下如：D:\soft\CTB
3. 配置环境变量：

四、使用CTB生成.terrain文件：

注意：执行前请确保输出文件夹存在且是空的，同时保证磁盘有足够的空间，全球数据大约800G。生成瓦片的等级不要太高，中国范围内单18级就需要大约700G。

打开cmd控制台,输入指令如：

ctb-tile -o D:\terrain\china\terrain -s 14 -e 0 -r nearest -c 4 C:\Users\Administrator\Documents\ArcGIS\dem_result\noData\new.tif

这个过程更加费时间，需要等待很久，直到执行结束。

五、发布并使用Cesium调用

1. 将生成的文件拷贝到Tomcat服务器下。

2. 拷贝“覆盖至生成tiles结果”目录中的layer.json文件和0文件夹至生成结果目录下：

3. 调用代码：

   1. var terrainLayer = new Cesium.CesiumTerrainProvider({
   2. url: “http://localhost:8080/terrain”, // 默认立体地表
   3. });
   4. scene.terrainProvider = terrainLayer;

   至此，地形发布完成。如有错误之处欢迎指正。

后记：有朋友反馈，太大的地形处理出来之后有问题，使用的时候会报错，这里我在给出一个简单的工具：CesiumLab

操作方式仍然需要执行完第二步，也就是整合地形文件为一个tif，然后使用CesiumLab的地形处理工具处理，亲测可用，0-12级，需要约2天时间。