Cesium Terrain Builder
输出瓦片默认是zib压缩后的,在业务中如果传输不是问题,反而增加浏览器的解压处理,希望能支持输出非压缩瓦片。
针对此需求,修改代码并重新编译。
一、代码分析
1、输出数据对象
文件格式:主要为heightmap-1.0、quantized-mesh-1.0,其类图结构如下:
- MeshTile为quantized-mesh-1.0格式数据对象。
- TerrainTile为heightmap-1.0格式数据对象。
2、瓦片生成器
- Meshtiler 用来生成MeshTile对象。
- TerrainTiler 用来生成TerrainTile对象。
3、迭代器和序列化器
4、ctb-tile.exe的主要流程
二、功能扩展
1、增加非压缩支持。
- ctb-tile.cpp中TerrainBuild类,增加参数
// 内部变量
bool gzib;
// 构造函数设置默认值
gzib(false)
// 设置压缩参数
static void setGzib(command_t* command) {
static_cast<TerrainBuild *>(Command::self(command))->gzib = true;
}
- ctb-tile.cpp中main函数解析参数
command.option("-G", "--gzip", "use zib compress file(defalut uncompress)", TerrainBuild::setGzib);
- CTBFileTileSerializer.cpp中,serializeTile函数增加zibFlag参数。根据gzibFlag决定是否使用CTBZFileOutputStream(压缩)或者CTBFileOutputStream(非压缩)
/**
* @details
* Serialize a MeshTile to the Directory store
* 增加zibFlag,是否输出压缩格式文件
*/
bool
ctb::CTBFileTileSerializer::serializeTile(const ctb::MeshTile *tile, bool writeVertexNormals, bool gzibFlag) {
const TileCoordinate *coordinate = tile;
const string filename = getTileFilename(coordinate, moutputDir, "terrain");
const string temp_filename = concat(filename, ".tmp");
if (gzibFlag) {
CTBZFileOutputStream ostream(temp_filename.c_str());
tile->writeFile(ostream, writeVertexNormals);
ostream.close();
} else {
FILE *fp = fopen(temp_filename.c_str(), "wb");
CTBFileOutputStream ostream(fp);
tile->writeFile(ostream, writeVertexNormals);
fclose(fp);
}
if (VSIRename(temp_filename.c_str(), filename.c_str()) != 0) {
throw new CTBException("Could not rename temporary file");
}
return true;
}
- ctb-tile.cpp中buildMesh函数,修改为传递gzlib函数
//serializer.serializeTile(tile, writeVertexNormals);
serializer.serializeTile(tile, writeVertexNormals, command->gzib);
2、输出指定范围的瓦片
- ctb-tile.cpp中TerrainBuild类,增加参数
// 内部变量
int layerStartX, layerStartY, layerEndX, layerEndY;
// 构造函数设置默认值
layerStartX(-1),
layerStartY(-1),
layerEndX(-1),
layerEndY(-1)
//设置参数函数
static void setLayerStartX(command_t *command) {
static_cast<TerrainBuild *>(Command::self(command))->layerStartX = atoi(command->arg);
}
static void setLayerStartY(command_t *command) {
static_cast<TerrainBuild *>(Command::self(command))->layerStartY = atoi(command->arg);
}
static void setLayerEndX(command_t *command) {
static_cast<TerrainBuild *>(Command::self(command))->layerEndX = atoi(command->arg);
}
static void setLayerEndY(command_t *command) {
static_cast<TerrainBuild *>(Command::self(command))->layerEndY = atoi(command->arg);
}
- ctb-tile.cpp中main函数解析参数
command.option("-A", "--start-x <sx>", "tile’s start x ", TerrainBuild::setLayerStartX);
command.option("-D", "--end-x <ex>", "tile’s end x", TerrainBuild::setLayerEndX);
command.option("-W", "--start-y <sy>", "tile’s start y", TerrainBuild::setLayerStartY);
command.option("-S", "--end-y <ey>", "tile’s end y", TerrainBuild::setLayerEndY);
- ctb-tile.cpp修改buildMesh方法
/// Output mesh tiles represented by a tiler to a directory
static void buildMesh(MeshSerializer &serializer, const MeshTiler &tiler, TerrainBuild *command, TerrainMetadata *metadata, bool writeVertexNormals = false) {
i_zoom startZoom = (command->startZoom < 0) ? tiler.maxZoomLevel() : command->startZoom,
endZoom = (command->endZoom < 0) ? 0 : command->endZoom;
// 取得参数中瓦片坐标
int layerStartX = (command->layerStartX < 0)? 0: command->layerStartX,
layerStartY = (command->layerStartY < 0) ? 0 : command->layerStartY,
layerEndX = (command->layerEndX < 0) ? 0 : command->layerEndX,
layerEndY = (command->layerEndY < 0) ? 0 : command->layerEndY;
MeshIterator iter(tiler, startZoom, endZoom);
int currentIndex = incrementIterator(iter, 0);
setIteratorSize(iter);
GDALDatasetReaderWithOverviews reader(tiler);
while (!iter.exhausted()) {
const TileCoordinate *coordinate = iter.GridIterator::operator*();
if (metadata) metadata->add(tiler.grid(), coordinate);
// 如果没有指定瓦片坐标或zoom级别不是一个,使用原来逻辑
if (layerStartX == 0 && layerEndX == 0 && startZoom != endZoom) {
if (serializer.mustSerializeCoordinate(coordinate)) {
MeshTile *tile = iter.operator*(&reader);
serializer.serializeTile(tile, writeVertexNormals, command->gzib);
delete tile;
}
}
else {
// 如果有瓦片坐标,只生成瓦片坐标范围内的文件
if (coordinate->x >= layerStartX && coordinate->x <= layerEndX && coordinate->y >= layerStartY && coordinate->y <= layerEndY) {
if (serializer.mustSerializeCoordinate(coordinate)) {
MeshTile *tile = iter.operator*(&reader);
serializer.serializeTile(tile, writeVertexNormals, command->gzib);
delete tile;
}
}
}
currentIndex = incrementIterator(iter, currentIndex);
showProgress(currentIndex);
}
}
三、命令行说明
编译成功后,就支持非压缩切片和生成指定范围的切片功能了。
- 调用说明
ctb-tile.exe -o D:/test/tile D:/test/dem.tif -f Mesh -s 6 -e 6 -A 105 -D 107 -W 45 -S 47
- 完整参数说明
参数 | 说明 |
---|---|
-V | 输出程序版本。 |
-h | 输出帮助信息。 |
-o | 切片输出路径。 |
-p | 切片规则。默认geodetic,可以指定mercator |
-f | 切片格式。默认Terrain,可以指定Mesh |
-c | 使用线程数。默认同cpu数 |
-t | 瓦片像素大小。默认65*65 |
-s | 瓦片开始级别。 |
-e | 瓦片结束级别。注意:开始级别 > 结束级别 |
-r | 重采样方法。默认平均值,可以指定min,max,average,nearest,bilinear,cubic等。 |
-R | 不覆盖现有文件。 |
-G | 生成zib压缩切片。 |
-A | 切片范围,x开始索引。 |
-D | 切片范围,x结束索引。 |
-W | 切片范围,y开始索引。 |
-S | 切片范围,y结束索引。 |