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

adminCesium 生成terrain地形数据—-CTB方式及步骤已关闭评论条评论 45 次浏览

背景:项目前端使用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值,处理公式为:

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控制台,输入指令如:

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

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

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

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