Abstract
  1. Cesium 生成 terrain 地形数据 ----CTB 方式及步骤
Keywords Cesium  Cesium 
Citation Yao Qing-sheng.Cesium 生成 Terrain 地形数据 ----CTB 方式及步骤.FUTURE & CIVILIZATION Natural/Social Philosophy & Infomation Sciences,20240808. https://yaoqs.github.io/20240808/cesium-sheng-cheng-terrain-di-xing-shu-ju-ctb-fang-shi-ji-bu-zou/

转载自 Cesium 生成 terrain 地形数据 ----CTB 方式及步骤

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

后记：如果你只是需要将 tif 处理成 terrain, 而不是学习处理过程，可以直接使用这个简单的工具：CesiumLab。不需要拼接成一个 TIF，也不需要处理 No-Data，地形功能免费。如果对你有用，点个赞加个关注吧๑乛◡乛๑

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

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

一、下载原始地形数据

网盘包含两部分内容：

  • 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 值，ArcToolbox->spatial Analyst Tool->Map Algebra->Raster Calculator 处理公式为:

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

三、配置 CTB 环境

  1. 下载 CTB 工具包：资源 （没有积分的可以加我 QQ）
  2. 将 ctb 和 gdal-data 拷贝至合适的目录下如：D:\soft\CTB
  3. 配置环境变量：                                                             

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

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

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

1
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
    2
    3
    4
    var terrainLayer = new Cesium.CesiumTerrainProvider({
            url: "http://localhost:8080/terrain", // 默认立体地表
         });
    scene.terrainProvider = terrainLayer;

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

References

*Corresponding Author: Yao Qing-sheng <Email: 350788415@qq.com> BY-NC-SA .

  1. Yao Qing-sheng.Cesium 生成 Terrain 地形数据 ----CTB 方式及步骤.FUTURE & CIVILIZATION.Natural/Social Philosophy & Infomation Sciences,20240808.https://yaoqs.github.io/20240808/cesium-sheng-cheng-terrain-di-xing-shu-ju-ctb-fang-shi-ji-bu-zou/
  2. 版权声明：本文为「LordYao」的原创文章，遵循 CC 4.0 BY-NC-SA 版权协议，转载请附上原文出处链接及本声明。
    原文链接： https://yaoqs.github.io/20240808/cesium-sheng-cheng-terrain-di-xing-shu-ju-ctb-fang-shi-ji-bu-zou/