2021/01/11/publish-dem-data-with-mapbox-terrain-rgb

[beginor]

2021/01/11/publish-dem-data-with-mapbox-terrain-rgb

[chenchaoqdu]

您好，我通过上面的文档发布了一个高程服务，通过mapbox引入，但是这块高程区域比其他图层高出一块来是怎么回事，能解答一下么

[beginor]

@chenchaoqdu 主要和你的栅格数据有关， 有几个可能出问题的地方：

• 栅格数据的 0 值是多少；
• 栅格数据的分辨率和它的缩放不匹配，也可能导致这样的问题；
• 这种 rgb 格式的高层， 精确度不高(分米级)，需要对栅格数据进行差值，尽量适配；

[chenchaoqdu]

谢谢解答，我再找一下原因

[felix0917]

@chenchaoqdu 我也遇到了相同的问题，请问问题原因是什么，解决了吗？感谢~

[Sun404]

@chenchaoqdu 请问解决了吗， 我也遇到类似问题，高程区域边界有明显的楞线

[beginor]

@chenchaoqdu @undefinedGIS @Sun404 能否将 dem 的 rio info 结果贴出来看看？

[chenchaoqdu]

{"blockxsize": 256, "blockysize": 256, "bounds": [13303711.89897836, 4242401.635796383, 13476483.548124677, 4459833.289155098], "colorinterp": ["gray"], "compress": "deflate", "count": 1, "crs": "EPSG:3857", "descriptions": [null], "driver": "GTiff", "dtype": "float32", "height": 5039, "indexes": [1], "interleave": "band", "lnglat": [120.28529440920339, 36.36612646280445], "mask_flags": [["all_valid"]], "nodata": null, "res": [43.14976252405522, 43.14976252405522], "shape": [5039, 4004], "tiled": true, "transform": [43.14976252405522, 0.0, 13303711.89897836, 0.0, -43.14976252405522, 4459833.289155098, 0.0, 0.0, 1.0], "units": [null], "width": 4004}

[beginor]

@chenchaoqdu 分辨率是 43.14976252405522 ， 不能对应到 web 墨卡托坐标系中对应图层的标准分辨率， 可以尝试对影像进行插值， 插值到最接近图层的标准分辨率， 然后再进行发布。

[chenchaoqdu]

这个影像如何插值啊，能不能讲一下

[beginor]

@chenchaoqdu 准确的说法应该是 栅格数据重采样 ，自行搜索， arcgis， qgis，gdal 都可以实现

[chenchaoqdu]

是不是用gdalwarp处理时，加上-r参数进行栅格数据重采样也可以？

[beginor]

试一下不就知道了， gdal 的文档比我要权威吧。

[wangsihao19981113]

您好，我最近通过你这个文章成功发布了类似的服务，进而使用了maptalks做出了地形三维可视化，我想写一篇类似的博客，先是阐述你这个方法，再阐述我使用maptalks做地形三维可视化的方法，借此想问一下能不能复述你的方法呢，当然我会附上参考的连接。

[beginor]

@wangsihao19981113 没问题的，不过maptalks似乎很少人使用啊

[wangsihao19981113]

maptalks确实用的人不是很多，但各方面拓展还算可以

[beginor]

请问你发表的文章的地址是什么，方便看一下么？