开篇

第一次写博客，有点紧张和激动。GIS群里的小伙伴都在猜我是谁，先留个悬念吧，哈哈。

最近正好在工作换挡期，又想把三维捡起来倒腾一下，就记录一下如何用Cesium构建三维场景。这个有很多前辈都已经实现过，但是自己尝试很多次，这次算是成功，也算是一个小小的里程碑。

首先是数据的准备

三维场景需要地形数据和叠加正摄影像，而Cesium支持的地形数据有两种，一种为STK World Terrain ，另一种是 Small Terrain。具体这两种的区别可参见前辈们的文章。

在Cesium群中，有群友Wonder贡献出一个getTerrain.py，用Python2.x来运行，可以下载.terrain的文件，但是用gallery中的Terrain示例却无法调用。

查找原因后无果，准备用Python3添加请求头再试试。

对请求头添加请求头，其中包括User-Agent、Accept、Accept-Encoding、Accept-Language等，请求地址为

https://assets.agi.com/stk-terrain/world/0/0/0.terrain?v=1.31376.0&f=TerrainTile

注意，下载下来的文件大小为3896字节，而请求页面上显示为2035字节。原来下载下来的文件经过gzip解压了，所以必须在python代码中添加gzip解压这段代码进去，否则下载下来的文件只有2035字节并且Cesium中无法使用

下载完terrain后，还需要下载layer.json这个可以找个STKTerrain的例子，获得该文件。

打开该layer.json并格式化（在线格式化工具http://tool.oschina.net/codeformat/js/），可以看出，现在的数据版本变化了

Small Terrain 是中等高分辨率,基于heightmap的地形，渲染出的地形效果不如quantized mesh的地形，但也基本能接受。

网上已经有一些开源的生成工具可以由DEM 数据生成这种规范的.terrain文件,因此我们下面重点说明这种类型的地形生成。

在易至天工开源的tiff转terrain源码中，有这样一段话

//前7层建议自己到到CES上去下，下载地址  string url = "http://cesiumjs.org/smallterrain/" + layer.ToString() + "/" + x.ToString() + "/" + y.ToString() + ".terrain";

所以前0-6级我们也可以通过python进行爬取

将上面的代码地址换一下即可，同样它也是需要gzip解压的。

接下来可以到SRTM（http://srtm.csi.cgiar.org/index.asp）上下载90m分辨率的DEM数据，并通过易至天工开源的源码进行转换成terrain文件，注意输入等级时候要按照以下格式

其次，通过Cesium加载三维场景

在Cesium 上加载本地地形Provider