CesiumLab 地形数据处理

最近接连有用户反应地形数据处理的各种问题，我也是各种测试，想想还是整理一个文档彻底说明一下。

地形栅格数据格式，一般是tif ，也有dem或者img，但是我个人强烈建议使用tif格式，因为cesiumlab都是用这个格式测试的。

我们先说一下地形栅格数据的几个关键参数：

强烈建议处理之前，先用arcmap打开，确保以下参数正常

1，空间参考（spatial reference）

spatial reference 定义了 栅格数据的坐标空间，下图为正常参数，如果arcmap 这里显示undefined，那表示空间参考未知，必然无法正常处理。

CesiumLab 地形数据处理

栅格的空间参考信息

如果tif中不含空间参考，但是空间参考数据处理人员又明确知道，那么可以在cesiumlab中指定。反过来说如果不确定，那就不要乱尝试，只会浪费处理时间。

CesiumLab 地形数据处理

cesiumlab中指定空间索引

2，无数据值（nodata）

nodata定义了删格中的无效数据数据值，也就是说如果栅格里某个像素的值 = nodata，意味着该像素没有数据。

对于nodata的处理，cesiumlab把该位置使用高程0填充。

但是原始数据中这个值很乱，完全是用户设置的，栅格里可能没有这个nodata信息，也可能为任意值。一般用户会设置一个超大的负数值，比如 -10000  ，-99999等。这里强烈建议，不要设置数据极值，例如-1.797693e+308 ，这种极值会导致中间计算出现非数字，导致地形处理异常。

CesiumLab 地形数据处理

无数据值

如果原始栅格中没有nodata值，只有一种情况处理地形高程会正常（数据的空间参考是wgs84，地形数据边界为矩形），除此之外都必须设置此值。在cesiumlab选中对应tif，如下设置。

CesiumLab 地形数据处理

这个-9999.005859是怎么来的？ 我们还是必须通过arcmap打开tif，使用indentify工具，点击我们明确知道是无数据区域，那么获得他的值。

CesiumLab 地形数据处理

arcgis里人工判定nodata

3，地理范围（extent)

栅格数据实际上是一张图片，必须要一个转换，把像素坐标转到地理坐标，这就是extent，在gdal里是geotransform。在arcmap里如下图所示

CesiumLab 地形数据处理

栅格的地理范围

正常这里应该显示的地理坐标，如果你发现这个extent的四个值里有0，那么基本意味着，你的tif的地理范围信息缺失。

tif 一般把这个信息存储在文件内，还有另一种方式存储为 同名的tfw 文件。

所以我们处理或者拷贝数据的时候，如果有tfw，请一并拷贝。

CesiumLab 地形数据处理

tfw文件

确保上面三个关键参数正确，我们再分析cesiumlab处理地形的异常情况原因

1， 瞬间结束      空间参考缺失  或者  地理范围缺失

2，处理完成之后半边黑球      nodata 过大

3，地形高程异常      nodata 缺失或者不正确

4，位置不正确或伴随高程不正常      一般是空间参考不正确根据异常原因，请在arcmap使用相关工具处理你的栅格数据。

©著作权归作者所有：来自51CTO博客作者cesium实验的原创作品，如需转载，请注明出处，否则将追究法律责任