• CesiumLab 地形数据处理


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

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

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

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

    1,空间参考(spatial reference)

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

     
    栅格的空间参考信息

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

     
    cesiumlab中指定空间索引

    2,无数据值(nodata)

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

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

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

     
    无数据值

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

     
     

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

     
    arcgis里人工判定nodata

    3,地理范围(extent)

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

     
    栅格的地理范围

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

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

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

     
    tfw文件

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

    1, 瞬间结束

          空间参考缺失  或者  地理范围缺失

    2,处理完成之后半边黑球

          nodata 过大

    3,地形高程异常

          nodata 缺失或者不正确

    4,位置不正确或伴随高程不正常

          一般是空间参考不正确

    根据异常原因,请在arcmap使用相关工具处理你的栅格数据。



  • 相关阅读:
    课程总结
    第十四周课程总结&实验报告
    第十三周课程总结
    第十二周课程总结
    第十一周课程总结
    第十周课程总结
    第七次java实验报告
    第六次java实验报告
    第五次java实验报告
    第四次java实验报告
  • 原文地址:https://www.cnblogs.com/cesium1/p/10062892.html
Copyright © 2020-2023  润新知