简介CTB地形切片生成器是一款面向Cesium三维GIS开发者的实用工具可将TIF高程影像切割为Cesium可直接加载的.terrain地形数据解决三维场景中地形瓦片制作流程繁琐的问题适合具备一定GIS基础、从事Web三维可视化开发的技术人员使用。资源包共108个文件约20.53MB以40个dll动态库和29个csv投影参数文件为核心辅以gfs、wkt、xsd等坐标与投影配置以及exe可执行程序、svg图标、txt说明和少量xml、json、ini等配置文件整体构成一套完整可运行的切片工具链。目前已有3493人学习下载说明其在Cesium地形处理场景中具有较高参考价值。借助该工具读者可省去自行编写切片脚本的环节直接完成从TIF到terrain的转换并结合配套博文理解投影参数与切片流程快速搭建三维地形底图提升项目开发效率。1. CTB地形切片生成器为什么你的三维地形加载总在关键时刻卡住做过三维 GIS 或者数字孪生项目的同行大概都遇到过这种场景地形数据一加载浏览器内存直接飙到几个 G稍微转动视角就掉帧移动端更是直接白屏。问题往往不在模型本身而在于地形数据的组织方式——你喂给渲染引擎的是一整块几十公里见方的 DEM 高程栅格它不卡谁卡。CTB 地形切片生成器要解决的就是这件事把大范围、高精度的地形高程数据按照瓦片金字塔结构切成一层层、一块块的小切片让前端按需加载、逐级细化。CTB 是 Cesium Terrain Builder 的缩写它输出的切片格式直接兼容 Cesium 的CesiumTerrainProvider是目前三维地球场景里最常用的离线地形方案之一。这篇文章面向的是需要自己动手把原始 DEM 转成可发布地形服务的工程师不管你是第一次接触 CTB还是已经用过但总在参数上翻车下面这些内容都能直接拿去复现。2. CTB 切片到底切了什么从 DEM 到瓦片金字塔的完整链路2.1 为什么不能直接把 DEM 丢给 Cesium很多人第一次做离线地形直觉就是把 GeoTIFF 往 Cesium 里一塞结果发现要么加载不出来要么加载出来是一整块平板。原因在于 Cesium 的地形渲染依赖的是分层分块的瓦片索引而不是一张完整的高程图。它需要知道当前视角下应该请求哪一级、哪一块瓦片每块瓦片对应地理范围是什么瓦片内部的高程值怎么排列。CTB 生成器做的就是这套转换工作——读取原始 DEM重投影到目标坐标系按 TMS 规则切分再编码成 Cesium 能识别的.terrain格式。这里有一个容易被忽略的点CTB 输出的瓦片不是图片而是二进制高程数据。每个.terrain文件里存的是该瓦片范围内按网格排列的高程值经过量化压缩后体积很小但精度足够渲染用。这也是为什么同样范围的地形切片后加载流畅直接加载原始 DEM 就卡死——前者是按需请求后者是全量加载。2.2 瓦片金字塔的层级与编号规则CTB 遵循的是 TMSTile Map Service规范层级从 0 开始第 0 层是整个世界的一张瓦片第 1 层切成 2×2第 n 层切成 2^n × 2^n。每个瓦片有唯一的z/x/y编号z是层级x是列号y是行号。Cesium 请求地形时会根据当前相机高度计算出需要的层级再根据视锥体范围算出需要的x/y范围只请求这些瓦片。CTB 生成器在切分时会同时生成一个layer.json文件里面记录了可用层级范围、瓦片格式、投影方式、可用范围等元信息。Cesium 加载地形服务时首先读的就是这个文件如果layer.json里的available数组和实际瓦片对不上就会出现某些区域加载不出来或者请求 404 的情况。2.3 用 ctb 命令行走通一次最小切片流程CTB 的常见用法是通过命令行工具ctb来执行。假设你手头有一个dem.tif想切成 Cesium 可用的地形瓦片最简流程如下# 第一步查看 DEM 基本信息确认坐标系和范围 gdalinfo dem.tif # 第二步如果 DEM 不是 WGS84 地理坐标系先重投影 gdalwarp -t_srs EPSG:4326 -r bilinear dem.tif dem_wgs84.tif # 第三步用 ctb 生成地形瓦片 ctb-tile -f Terrain -C -N -o ./terrain_output dem_wgs84.tif # 第四步生成 layer.json 索引文件 ctb-tile -f Terrain -C -N -l -o ./terrain_output dem_wgs84.tif这几条命令里-f Terrain指定输出格式为 Cesium 地形-C表示输出为紧凑格式减少文件体积-N表示不生成法线贴图除非你需要光照效果-o指定输出目录-l表示只生成layer.json而不重新切片。逻辑说明gdalinfo和gdalwarp是 GDAL 工具链的命令前者用来确认源数据的坐标系和范围后者用来做重投影。CTB 要求输入 DEM 必须是 WGS84 地理坐标系EPSG:4326如果源数据是投影坐标系比如 UTM不重投影直接切片会导致瓦片位置偏移甚至切片失败。ctb-tile是 CTB 的核心命令第一次执行切片第二次执行只更新索引。参数说明-f Terrain是固定写法表示输出 Cesium 地形格式-C开启紧凑模式每个瓦片文件更小适合网络传输-N跳过法线生成如果你不需要地形光照可以加上能省不少时间-o后面跟输出目录目录不存在会自动创建-l只生成layer.json不重新切片适合在已有瓦片基础上补索引。提示如果你的 DEM 范围很大切片时间可能从几十分钟到几个小时不等。建议先用gdal_translate裁一小块测试确认流程走通再跑全量。2.4 切片输出目录结构长什么样切片完成后输出目录的结构大致如下terrain_output/ ├── layer.json ├── 0/ │ └── 0/ │ └── 0.terrain ├── 1/ │ ├── 0/ │ │ ├── 0.terrain │ │ └── 1.terrain │ └── 1/ │ ├── 0.terrain │ └── 1.terrain └── ...每一级目录对应一个层级层级目录下再按x/y分目录最终文件是.terrain后缀。layer.json在根目录记录了所有可用层级和范围。部署时把这个目录整体放到静态服务器上Cesium 里用CesiumTerrainProvider指向这个目录的 URL 就能加载。3. 参数怎么调影响切片质量和性能的五个关键设置3.1 层级范围不是越多越好CTB 默认会根据 DEM 分辨率自动计算最大层级但你可以通过-z参数手动指定。层级越多地形细节越丰富但瓦片数量呈指数增长。第 10 层有 100 多万张瓦片第 12 层就超过 1600 万张。对于大多数项目10 到 12 层已经足够再高不仅切片慢前端请求也会成为瓶颈。我一般会先算一下DEM 的分辨率是 30 米那么第 12 层的瓦片分辨率大约是 38 米已经接近原始精度再切更高层级就是插值放大没有实际意义。所以最大层级应该由 DEM 分辨率决定而不是拍脑袋定。3.2 重采样方法bilinear 还是 cubic在gdalwarp重投影阶段-r参数控制重采样方法。常见的有near最近邻、bilinear双线性、cubic三次卷积。对于地形高程数据bilinear是比较稳妥的选择它在平滑度和计算量之间取得了平衡。near会保留原始值但可能产生锯齿cubic更平滑但计算更慢而且在边缘处可能产生过冲。如果你的 DEM 本身质量很高用bilinear就够了。如果发现切片后地形有明显的阶梯感可以试试cubic但要注意检查边缘区域是否有异常值。3.3 紧凑模式与法线贴图什么时候该开什么时候该关-C紧凑模式建议始终开启它会把瓦片数据用 gzip 压缩后存储体积能减少 60% 以上前端加载也更快。-N法线贴图则要看需求如果你需要在 Cesium 里开启地形光照enableLighting就需要法线数据如果只是普通渲染关掉能省不少切片时间和存储空间。实际项目中我见过有人为了省事关掉了法线结果后来想加光照效果只能重新切一遍。所以如果项目后期可能用到光照建议一开始就保留法线。3.4 并行切片用 -p 参数加速CTB 支持多线程切片-p参数指定并行进程数。默认是单线程对于大范围 DEM 来说太慢。一般设置为 CPU 核心数即可比如 8 核机器用-p 8。但要注意并行切片对内存消耗也更大如果 DEM 很大而内存有限反而会拖慢速度甚至崩溃。# 8 线程并行切片适合大范围 DEM ctb-tile -f Terrain -C -N -p 8 -o ./terrain_output dem_wgs84.tif3.5 输出格式Terrain 还是 Quantized-MeshCTB 支持两种输出格式Terrain和Quantized-Mesh。前者是 Cesium 旧版地形格式兼容性好但体积较大后者是新版格式支持更高效的数据压缩和更丰富的元信息但需要 Cesium 1.40 以上版本。新项目建议直接用Quantized-Mesh老项目如果已经用了Terrain格式迁移成本不高但也没必要强行换。参数作用推荐值注意事项-f输出格式Quantized-Mesh新项目首选旧版 Cesium 用 Terrain-C紧凑模式开启减少体积建议始终开启-N跳过法线按需需要光照则不加此参数-p并行进程数CPU 核心数内存不足时降低-z层级范围根据 DEM 分辨率定不是越多越好4. 避坑指南CTB 切片最常见的五个翻车现场4.1 切片成功但 Cesium 加载后地形是平的现象切片命令跑完了目录里也有.terrain文件但 Cesium 加载后地形没有任何起伏像一张平面。原因最常见的情况是 DEM 的坐标系不是 WGS84切片时没有重投影导致瓦片位置和实际地理范围对不上。另一种可能是layer.json里的available数组为空或不完整Cesium 找不到可用瓦片。解决先用gdalinfo确认 DEM 坐标系如果不是 EPSG:4326必须先用gdalwarp重投影。然后检查layer.json里的available字段确保每个层级的startX/startY/endX/endY范围正确。如果available为空重新执行带-l参数的命令生成索引。4.2 部分区域加载不出来控制台报 404现象地形整体能加载但某些区域始终是空白浏览器控制台能看到.terrain文件请求返回 404。原因通常是layer.json里记录的可用范围和实际瓦片文件不匹配。比如切片时用了-z限制了层级范围但layer.json没有更新Cesium 请求了不存在的瓦片。解决删除旧的layer.json重新执行ctb-tile -l生成索引。如果问题依旧检查切片输出目录里对应层级的x/y目录是否存在确认瓦片文件确实生成成功。4.3 切片速度极慢几个小时跑不完现象一个中等范围的 DEM 切片跑了几个小时还没结束CPU 占用率很低。原因默认单线程切片没有加-p参数。另外如果 DEM 分辨率很高、范围很大切片本身就是计算密集型任务单线程跑确实慢。解决加上-p参数设置为 CPU 核心数。如果内存充足可以适当提高并行数。另外建议先用gdal_translate裁剪出感兴趣区域减少不必要的切片范围。4.4 前端加载地形后内存暴涨现象Cesium 加载地形后浏览器内存持续上升转动视角一段时间后页面卡死。原因通常是切片层级过高前端请求了过多瓦片。Cesium 默认会根据视角动态请求但如果最大层级设得太高相机稍微拉远就会请求大量瓦片。解决降低最大层级或者在 Cesium 里设置maximumScreenSpaceError参数控制瓦片请求的激进程度。另外确保服务器开启了 gzip 压缩减少传输体积。4.5 layer.json 里的 projection 字段写错了现象Cesium 加载地形时报错提示投影不匹配或无法解析。原因layer.json里的projection字段必须是EPSG:4326如果切片时源数据坐标系不对这个字段可能写成其他值。解决手动检查layer.json确保projection为EPSG:4326format为quantized-mesh-1.0或heightmap-1.0scheme为tms。如果不对重新切片或手动修正。5. 进阶技巧用 Python 脚本批量切片并自动校验5.1 批量切片脚本的骨架实际项目中往往需要处理多个 DEM 文件手动一条条跑命令效率太低。我一般会写一个 Python 脚本自动完成重投影、切片、索引生成和校验。下面是一个简化版的骨架import subprocess import os import json def reproject_dem(src, dst): 重投影到 WGS84 cmd [ gdalwarp, -t_srs, EPSG:4326, -r, bilinear, src, dst ] subprocess.run(cmd, checkTrue) def generate_terrain(dem_path, output_dir, max_zoom12, parallel8): 生成地形瓦片 cmd [ ctb-tile, -f, Quantized-Mesh, -C, -N, -p, str(parallel), -z, str(max_zoom), -o, output_dir, dem_path ] subprocess.run(cmd, checkTrue) # 生成 layer.json cmd_index [ ctb-tile, -f, Quantized-Mesh, -C, -N, -l, -o, output_dir, dem_path ] subprocess.run(cmd_index, checkTrue) def validate_layer_json(output_dir): 校验 layer.json 关键字段 layer_path os.path.join(output_dir, layer.json) with open(layer_path, r) as f: data json.load(f) assert data[projection] EPSG:4326, 投影字段错误 assert data[scheme] tms, 瓦片方案错误 assert len(data[available]) 0, available 数组为空 print(layer.json 校验通过) if __name__ __main__: src_dem dem.tif wgs84_dem dem_wgs84.tif out_dir ./terrain_output reproject_dem(src_dem, wgs84_dem) generate_terrain(wgs84_dem, out_dir) validate_layer_json(out_dir)逻辑说明脚本分三步——重投影、切片、校验。reproject_dem调用gdalwarp把源 DEM 转到 WGS84generate_terrain调用ctb-tile切片并生成索引validate_layer_json读取layer.json检查关键字段。这样每次处理新数据只需要改几个路径参数不用重复敲命令。参数说明max_zoom控制最大层级根据 DEM 分辨率调整parallel是并行进程数建议设为 CPU 核心数validate_layer_json里的断言可以根据项目需求增减比如还可以检查available数组的层级范围是否连续。5.2 切片完成后的自动化校验切片跑完不代表万事大吉我习惯加一步自动校验随机抽取几个瓦片文件检查文件大小是否在合理范围内太小可能是空瓦片太大可能是数据异常再检查layer.json里的层级范围是否和实际目录一致。import random def spot_check_tiles(output_dir, sample_count10): 随机抽查瓦片文件大小 terrain_files [] for root, dirs, files in os.walk(output_dir): for f in files: if f.endswith(.terrain): terrain_files.append(os.path.join(root, f)) if len(terrain_files) sample_count: sample_count len(terrain_files) samples random.sample(terrain_files, sample_count) for path in samples: size os.path.getsize(path) if size 100: print(f警告{path} 文件过小可能为空瓦片) elif size 1024 * 1024: print(f警告{path} 文件过大检查数据是否异常) else: print(f正常{path} 大小 {size} 字节)这个抽查逻辑很简单但能提前发现大部分切片异常。比如某个区域 DEM 有空洞切片后对应瓦片会特别小或者数据范围计算错误导致某些瓦片异常大。5.3 部署时的几个实用建议切片输出目录部署到静态服务器时有几个细节值得注意。第一确保服务器开启了 gzip 压缩.terrain文件本身已经是压缩格式但layer.json是纯文本压缩后传输更快。第二设置正确的 MIME 类型.terrain文件建议用application/octet-stream避免浏览器尝试解析。第三如果瓦片数量很大建议用 CDN 或者对象存储避免单台服务器磁盘 IO 成为瓶颈。另外Cesium 加载地形时默认会请求layer.json如果这个文件更新了但浏览器有缓存可能会导致加载异常。部署时可以在 URL 后面加版本号参数强制刷新缓存。5.4 一个容易被忽略的细节NoData 值的处理DEM 数据里经常有 NoData 值比如 -9999表示该区域没有有效高程。CTB 切片时如果不对 NoData 做处理这些值会被当成真实高程导致地形出现异常凹陷或尖刺。常见做法是在gdalwarp阶段用-dstnodata参数指定 NoData 值或者在切片前用gdal_fillnodata填充空洞。我一般会在重投影之后加一步检查# 查看 DEM 的 NoData 值 gdalinfo -stats dem_wgs84.tif | grep NoData # 如果 NoData 值不合理用 gdal_translate 重新指定 gdal_translate -a_nodata -9999 dem_wgs84.tif dem_fixed.tif这个细节看起来小但实际项目中因为 NoData 导致地形异常的案例不少。提前处理能省掉很多排查时间。5.5 关于 CTB 版本选择的一点经验CTB 有几个不同的发行版本功能和支持的格式略有差异。新版本通常支持 Quantized-Mesh 格式和更好的并行切片但可能对某些旧版 Cesium 兼容性稍差。我的习惯是新项目直接用最新稳定版老项目如果已经跑通就不要轻易升级除非有明确的需求比如需要 Quantized-Mesh 格式。升级前建议先用一小块 DEM 做对比测试确认切片结果和加载效果没有异常再全量迁移。毕竟切片一次不容易跑几个小时甚至一整天的都有因为版本问题翻车就太不划算了。希望帮到你。本文还有配套的精品资源点击获取
简介CTB地形切片生成器是一款面向Cesium三维GIS开发者的实用工具可将TIF高程影像切割为Cesium可直接加载的.terrain地形数据解决三维场景中地形瓦片制作流程繁琐的问题适合具备一定GIS基础、从事Web三维可视化开发的技术人员使用。资源包共108个文件约20.53MB以40个dll动态库和29个csv投影参数文件为核心辅以gfs、wkt、xsd等坐标与投影配置以及exe可执行程序、svg图标、txt说明和少量xml、json、ini等配置文件整体构成一套完整可运行的切片工具链。目前已有3493人学习下载说明其在Cesium地形处理场景中具有较高参考价值。借助该工具读者可省去自行编写切片脚本的环节直接完成从TIF到terrain的转换并结合配套博文理解投影参数与切片流程快速搭建三维地形底图提升项目开发效率。1. CTB地形切片生成器为什么你的三维地形加载总在关键时刻卡住做过三维 GIS 或者数字孪生项目的同行大概都遇到过这种场景地形数据一加载浏览器内存直接飙到几个 G稍微转动视角就掉帧移动端更是直接白屏。问题往往不在模型本身而在于地形数据的组织方式——你喂给渲染引擎的是一整块几十公里见方的 DEM 高程栅格它不卡谁卡。CTB 地形切片生成器要解决的就是这件事把大范围、高精度的地形高程数据按照瓦片金字塔结构切成一层层、一块块的小切片让前端按需加载、逐级细化。CTB 是 Cesium Terrain Builder 的缩写它输出的切片格式直接兼容 Cesium 的CesiumTerrainProvider是目前三维地球场景里最常用的离线地形方案之一。这篇文章面向的是需要自己动手把原始 DEM 转成可发布地形服务的工程师不管你是第一次接触 CTB还是已经用过但总在参数上翻车下面这些内容都能直接拿去复现。2. CTB 切片到底切了什么从 DEM 到瓦片金字塔的完整链路2.1 为什么不能直接把 DEM 丢给 Cesium很多人第一次做离线地形直觉就是把 GeoTIFF 往 Cesium 里一塞结果发现要么加载不出来要么加载出来是一整块平板。原因在于 Cesium 的地形渲染依赖的是分层分块的瓦片索引而不是一张完整的高程图。它需要知道当前视角下应该请求哪一级、哪一块瓦片每块瓦片对应地理范围是什么瓦片内部的高程值怎么排列。CTB 生成器做的就是这套转换工作——读取原始 DEM重投影到目标坐标系按 TMS 规则切分再编码成 Cesium 能识别的.terrain格式。这里有一个容易被忽略的点CTB 输出的瓦片不是图片而是二进制高程数据。每个.terrain文件里存的是该瓦片范围内按网格排列的高程值经过量化压缩后体积很小但精度足够渲染用。这也是为什么同样范围的地形切片后加载流畅直接加载原始 DEM 就卡死——前者是按需请求后者是全量加载。2.2 瓦片金字塔的层级与编号规则CTB 遵循的是 TMSTile Map Service规范层级从 0 开始第 0 层是整个世界的一张瓦片第 1 层切成 2×2第 n 层切成 2^n × 2^n。每个瓦片有唯一的z/x/y编号z是层级x是列号y是行号。Cesium 请求地形时会根据当前相机高度计算出需要的层级再根据视锥体范围算出需要的x/y范围只请求这些瓦片。CTB 生成器在切分时会同时生成一个layer.json文件里面记录了可用层级范围、瓦片格式、投影方式、可用范围等元信息。Cesium 加载地形服务时首先读的就是这个文件如果layer.json里的available数组和实际瓦片对不上就会出现某些区域加载不出来或者请求 404 的情况。2.3 用 ctb 命令行走通一次最小切片流程CTB 的常见用法是通过命令行工具ctb来执行。假设你手头有一个dem.tif想切成 Cesium 可用的地形瓦片最简流程如下# 第一步查看 DEM 基本信息确认坐标系和范围 gdalinfo dem.tif # 第二步如果 DEM 不是 WGS84 地理坐标系先重投影 gdalwarp -t_srs EPSG:4326 -r bilinear dem.tif dem_wgs84.tif # 第三步用 ctb 生成地形瓦片 ctb-tile -f Terrain -C -N -o ./terrain_output dem_wgs84.tif # 第四步生成 layer.json 索引文件 ctb-tile -f Terrain -C -N -l -o ./terrain_output dem_wgs84.tif这几条命令里-f Terrain指定输出格式为 Cesium 地形-C表示输出为紧凑格式减少文件体积-N表示不生成法线贴图除非你需要光照效果-o指定输出目录-l表示只生成layer.json而不重新切片。逻辑说明gdalinfo和gdalwarp是 GDAL 工具链的命令前者用来确认源数据的坐标系和范围后者用来做重投影。CTB 要求输入 DEM 必须是 WGS84 地理坐标系EPSG:4326如果源数据是投影坐标系比如 UTM不重投影直接切片会导致瓦片位置偏移甚至切片失败。ctb-tile是 CTB 的核心命令第一次执行切片第二次执行只更新索引。参数说明-f Terrain是固定写法表示输出 Cesium 地形格式-C开启紧凑模式每个瓦片文件更小适合网络传输-N跳过法线生成如果你不需要地形光照可以加上能省不少时间-o后面跟输出目录目录不存在会自动创建-l只生成layer.json不重新切片适合在已有瓦片基础上补索引。提示如果你的 DEM 范围很大切片时间可能从几十分钟到几个小时不等。建议先用gdal_translate裁一小块测试确认流程走通再跑全量。2.4 切片输出目录结构长什么样切片完成后输出目录的结构大致如下terrain_output/ ├── layer.json ├── 0/ │ └── 0/ │ └── 0.terrain ├── 1/ │ ├── 0/ │ │ ├── 0.terrain │ │ └── 1.terrain │ └── 1/ │ ├── 0.terrain │ └── 1.terrain └── ...每一级目录对应一个层级层级目录下再按x/y分目录最终文件是.terrain后缀。layer.json在根目录记录了所有可用层级和范围。部署时把这个目录整体放到静态服务器上Cesium 里用CesiumTerrainProvider指向这个目录的 URL 就能加载。3. 参数怎么调影响切片质量和性能的五个关键设置3.1 层级范围不是越多越好CTB 默认会根据 DEM 分辨率自动计算最大层级但你可以通过-z参数手动指定。层级越多地形细节越丰富但瓦片数量呈指数增长。第 10 层有 100 多万张瓦片第 12 层就超过 1600 万张。对于大多数项目10 到 12 层已经足够再高不仅切片慢前端请求也会成为瓶颈。我一般会先算一下DEM 的分辨率是 30 米那么第 12 层的瓦片分辨率大约是 38 米已经接近原始精度再切更高层级就是插值放大没有实际意义。所以最大层级应该由 DEM 分辨率决定而不是拍脑袋定。3.2 重采样方法bilinear 还是 cubic在gdalwarp重投影阶段-r参数控制重采样方法。常见的有near最近邻、bilinear双线性、cubic三次卷积。对于地形高程数据bilinear是比较稳妥的选择它在平滑度和计算量之间取得了平衡。near会保留原始值但可能产生锯齿cubic更平滑但计算更慢而且在边缘处可能产生过冲。如果你的 DEM 本身质量很高用bilinear就够了。如果发现切片后地形有明显的阶梯感可以试试cubic但要注意检查边缘区域是否有异常值。3.3 紧凑模式与法线贴图什么时候该开什么时候该关-C紧凑模式建议始终开启它会把瓦片数据用 gzip 压缩后存储体积能减少 60% 以上前端加载也更快。-N法线贴图则要看需求如果你需要在 Cesium 里开启地形光照enableLighting就需要法线数据如果只是普通渲染关掉能省不少切片时间和存储空间。实际项目中我见过有人为了省事关掉了法线结果后来想加光照效果只能重新切一遍。所以如果项目后期可能用到光照建议一开始就保留法线。3.4 并行切片用 -p 参数加速CTB 支持多线程切片-p参数指定并行进程数。默认是单线程对于大范围 DEM 来说太慢。一般设置为 CPU 核心数即可比如 8 核机器用-p 8。但要注意并行切片对内存消耗也更大如果 DEM 很大而内存有限反而会拖慢速度甚至崩溃。# 8 线程并行切片适合大范围 DEM ctb-tile -f Terrain -C -N -p 8 -o ./terrain_output dem_wgs84.tif3.5 输出格式Terrain 还是 Quantized-MeshCTB 支持两种输出格式Terrain和Quantized-Mesh。前者是 Cesium 旧版地形格式兼容性好但体积较大后者是新版格式支持更高效的数据压缩和更丰富的元信息但需要 Cesium 1.40 以上版本。新项目建议直接用Quantized-Mesh老项目如果已经用了Terrain格式迁移成本不高但也没必要强行换。参数作用推荐值注意事项-f输出格式Quantized-Mesh新项目首选旧版 Cesium 用 Terrain-C紧凑模式开启减少体积建议始终开启-N跳过法线按需需要光照则不加此参数-p并行进程数CPU 核心数内存不足时降低-z层级范围根据 DEM 分辨率定不是越多越好4. 避坑指南CTB 切片最常见的五个翻车现场4.1 切片成功但 Cesium 加载后地形是平的现象切片命令跑完了目录里也有.terrain文件但 Cesium 加载后地形没有任何起伏像一张平面。原因最常见的情况是 DEM 的坐标系不是 WGS84切片时没有重投影导致瓦片位置和实际地理范围对不上。另一种可能是layer.json里的available数组为空或不完整Cesium 找不到可用瓦片。解决先用gdalinfo确认 DEM 坐标系如果不是 EPSG:4326必须先用gdalwarp重投影。然后检查layer.json里的available字段确保每个层级的startX/startY/endX/endY范围正确。如果available为空重新执行带-l参数的命令生成索引。4.2 部分区域加载不出来控制台报 404现象地形整体能加载但某些区域始终是空白浏览器控制台能看到.terrain文件请求返回 404。原因通常是layer.json里记录的可用范围和实际瓦片文件不匹配。比如切片时用了-z限制了层级范围但layer.json没有更新Cesium 请求了不存在的瓦片。解决删除旧的layer.json重新执行ctb-tile -l生成索引。如果问题依旧检查切片输出目录里对应层级的x/y目录是否存在确认瓦片文件确实生成成功。4.3 切片速度极慢几个小时跑不完现象一个中等范围的 DEM 切片跑了几个小时还没结束CPU 占用率很低。原因默认单线程切片没有加-p参数。另外如果 DEM 分辨率很高、范围很大切片本身就是计算密集型任务单线程跑确实慢。解决加上-p参数设置为 CPU 核心数。如果内存充足可以适当提高并行数。另外建议先用gdal_translate裁剪出感兴趣区域减少不必要的切片范围。4.4 前端加载地形后内存暴涨现象Cesium 加载地形后浏览器内存持续上升转动视角一段时间后页面卡死。原因通常是切片层级过高前端请求了过多瓦片。Cesium 默认会根据视角动态请求但如果最大层级设得太高相机稍微拉远就会请求大量瓦片。解决降低最大层级或者在 Cesium 里设置maximumScreenSpaceError参数控制瓦片请求的激进程度。另外确保服务器开启了 gzip 压缩减少传输体积。4.5 layer.json 里的 projection 字段写错了现象Cesium 加载地形时报错提示投影不匹配或无法解析。原因layer.json里的projection字段必须是EPSG:4326如果切片时源数据坐标系不对这个字段可能写成其他值。解决手动检查layer.json确保projection为EPSG:4326format为quantized-mesh-1.0或heightmap-1.0scheme为tms。如果不对重新切片或手动修正。5. 进阶技巧用 Python 脚本批量切片并自动校验5.1 批量切片脚本的骨架实际项目中往往需要处理多个 DEM 文件手动一条条跑命令效率太低。我一般会写一个 Python 脚本自动完成重投影、切片、索引生成和校验。下面是一个简化版的骨架import subprocess import os import json def reproject_dem(src, dst): 重投影到 WGS84 cmd [ gdalwarp, -t_srs, EPSG:4326, -r, bilinear, src, dst ] subprocess.run(cmd, checkTrue) def generate_terrain(dem_path, output_dir, max_zoom12, parallel8): 生成地形瓦片 cmd [ ctb-tile, -f, Quantized-Mesh, -C, -N, -p, str(parallel), -z, str(max_zoom), -o, output_dir, dem_path ] subprocess.run(cmd, checkTrue) # 生成 layer.json cmd_index [ ctb-tile, -f, Quantized-Mesh, -C, -N, -l, -o, output_dir, dem_path ] subprocess.run(cmd_index, checkTrue) def validate_layer_json(output_dir): 校验 layer.json 关键字段 layer_path os.path.join(output_dir, layer.json) with open(layer_path, r) as f: data json.load(f) assert data[projection] EPSG:4326, 投影字段错误 assert data[scheme] tms, 瓦片方案错误 assert len(data[available]) 0, available 数组为空 print(layer.json 校验通过) if __name__ __main__: src_dem dem.tif wgs84_dem dem_wgs84.tif out_dir ./terrain_output reproject_dem(src_dem, wgs84_dem) generate_terrain(wgs84_dem, out_dir) validate_layer_json(out_dir)逻辑说明脚本分三步——重投影、切片、校验。reproject_dem调用gdalwarp把源 DEM 转到 WGS84generate_terrain调用ctb-tile切片并生成索引validate_layer_json读取layer.json检查关键字段。这样每次处理新数据只需要改几个路径参数不用重复敲命令。参数说明max_zoom控制最大层级根据 DEM 分辨率调整parallel是并行进程数建议设为 CPU 核心数validate_layer_json里的断言可以根据项目需求增减比如还可以检查available数组的层级范围是否连续。5.2 切片完成后的自动化校验切片跑完不代表万事大吉我习惯加一步自动校验随机抽取几个瓦片文件检查文件大小是否在合理范围内太小可能是空瓦片太大可能是数据异常再检查layer.json里的层级范围是否和实际目录一致。import random def spot_check_tiles(output_dir, sample_count10): 随机抽查瓦片文件大小 terrain_files [] for root, dirs, files in os.walk(output_dir): for f in files: if f.endswith(.terrain): terrain_files.append(os.path.join(root, f)) if len(terrain_files) sample_count: sample_count len(terrain_files) samples random.sample(terrain_files, sample_count) for path in samples: size os.path.getsize(path) if size 100: print(f警告{path} 文件过小可能为空瓦片) elif size 1024 * 1024: print(f警告{path} 文件过大检查数据是否异常) else: print(f正常{path} 大小 {size} 字节)这个抽查逻辑很简单但能提前发现大部分切片异常。比如某个区域 DEM 有空洞切片后对应瓦片会特别小或者数据范围计算错误导致某些瓦片异常大。5.3 部署时的几个实用建议切片输出目录部署到静态服务器时有几个细节值得注意。第一确保服务器开启了 gzip 压缩.terrain文件本身已经是压缩格式但layer.json是纯文本压缩后传输更快。第二设置正确的 MIME 类型.terrain文件建议用application/octet-stream避免浏览器尝试解析。第三如果瓦片数量很大建议用 CDN 或者对象存储避免单台服务器磁盘 IO 成为瓶颈。另外Cesium 加载地形时默认会请求layer.json如果这个文件更新了但浏览器有缓存可能会导致加载异常。部署时可以在 URL 后面加版本号参数强制刷新缓存。5.4 一个容易被忽略的细节NoData 值的处理DEM 数据里经常有 NoData 值比如 -9999表示该区域没有有效高程。CTB 切片时如果不对 NoData 做处理这些值会被当成真实高程导致地形出现异常凹陷或尖刺。常见做法是在gdalwarp阶段用-dstnodata参数指定 NoData 值或者在切片前用gdal_fillnodata填充空洞。我一般会在重投影之后加一步检查# 查看 DEM 的 NoData 值 gdalinfo -stats dem_wgs84.tif | grep NoData # 如果 NoData 值不合理用 gdal_translate 重新指定 gdal_translate -a_nodata -9999 dem_wgs84.tif dem_fixed.tif这个细节看起来小但实际项目中因为 NoData 导致地形异常的案例不少。提前处理能省掉很多排查时间。5.5 关于 CTB 版本选择的一点经验CTB 有几个不同的发行版本功能和支持的格式略有差异。新版本通常支持 Quantized-Mesh 格式和更好的并行切片但可能对某些旧版 Cesium 兼容性稍差。我的习惯是新项目直接用最新稳定版老项目如果已经跑通就不要轻易升级除非有明确的需求比如需要 Quantized-Mesh 格式。升级前建议先用一小块 DEM 做对比测试确认切片结果和加载效果没有异常再全量迁移。毕竟切片一次不容易跑几个小时甚至一整天的都有因为版本问题翻车就太不划算了。希望帮到你。本文还有配套的精品资源点击获取
