leaflet显示高程

很多地图软件都能随鼠标移动动态显示高程。这里介绍一种方法，我所得出的。

1 下载高程数据

一般有12.5m数据下载，可惜精度根本不够，比如mapbox的免费在线的，或者91卫图提供百度网盘打包下载的，没法用，差距太大。而91卫图自己使用的，精度估计是5m的，是可以的。

所以还是用框选的办法在91中下载，你没有企业版也下载不了，再找其他办法吧。

可以一次下载全省：点广东省——下载……数据太大，不好处理

所以还是放大到8级，把经纬度打开，显示几个经纬度，按框选下载18级的dem高程数据

广东省分成了这么十几个框

2 导出geotif文件

下载后导出数据，每个导出文件为1G多，为tif文件

3 将tif转换为terrainRGB

将tif文件其中一个拷贝到dem2terrainRGB文件夹内

用命令

D:\dem2terrainrgb>python main.py --dem ./liyuzhou.tif --dist ./tiles --zoom 18
filled NODATA value successfully: ./tmp/liyuzhou_without_nodata.tif
rgbified successfully: ./tmp/liyuzhou_RGB.tif
created tileset successfully: ./tiles

生成**_without_nodata.tif（907M）和**_RGB.tif（279M），再切割成瓦片到tiles目录下

这个过程请参考我的文档。

4 重命名

由于它生成的瓦片命名是TMS标准的，我们leaflet需要WMTS标准的，所以要重新命名

代码写好了，不过是golang语言的，编译好的exe也提供了，放到与tiles文件夹同级的目录下运行exe，它自动递归tiles文件夹里所有文件，根据tiles文件夹里目录的数字，对png文件命名进行转换，自己看代码吧。

5 将png文件发布

比如放到engineercms的static目录下某个文件夹里

地址写好放到leaflet框架下的js文件中调用。

6 用leaflet-topography获得高程

用leaflet这个插件[leaflet-topography](https://github.com/slutske22/leaflet-topography)，实现鼠标移动，高程数据不断变化。

原理就是鼠标移动后，取得鼠标所在经纬度，请求这个经纬度所在的png文件，这个插件根据png文件的RGB三色值，计算得出高程数值，显示到页面。

7 地形切剖面

思路：取得lineString对象，用turf.length和turf.lineChunk进行分割，返回分割后的线段对象

用turf.coordEach循环分割后的对象，根据循环返回值其中的坐标，用eaflet-topography取得高程，并将高程赋予坐标中。

再用Leaflet.Elevation生成剖面。

    // let line = turf.lineString(data.data);var line = turf.lineString([[-95, 40],[-93, 45],[-85, 50]]);// var chunk = turf.lineChunk(line, 15, {units: 'miles'});var options = { units: 'kilometers' };var length = turf.length(line, options);var along = turf.lineChunk(line, length / 100, options);var bbox = turf.bbox(line);console.log(along);// var features = turf.featureCollection([//   turf.point([26, 37], {"foo": "bar"}),//   turf.point([36, 53], {"hello": "world"}),//   line,// ]);turf.coordEach(along, function(currentCoord, coordIndex, featureIndex, multiFeatureIndex, geometryIndex) {currentCoord[2] = 333;//=currentCoord//=coordIndex//=featureIndex//=multiFeatureIndex//=geometryIndex});console.log(along);