前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >leaflet显示高程

leaflet显示高程

作者头像
hotqin888
发布2023-02-26 21:09:13
1.3K0
发布2023-02-26 21:09:13
举报
文章被收录于专栏:hotqin888的专栏hotqin888的专栏

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

1 下载高程数据

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

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

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

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

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

2 导出geotif文件

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

3 将tif转换为terrainRGB

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

用命令

代码语言:javascript
复制
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],实现鼠标移动,高程数据不断变化。

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

7 地形切剖面

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

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

再用Leaflet.Elevation生成剖面。

代码语言:javascript
复制
    // 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);
本文参与 腾讯云自媒体分享计划,分享自作者个人站点/博客。
原始发表:2023-02-11,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 作者个人站点/博客 前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体分享计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
目录
  • 1 下载高程数据
  • 2 导出geotif文件
  • 3 将tif转换为terrainRGB
  • 4 重命名
  • 5 将png文件发布
  • 6 用leaflet-topography获得高程
  • 7 地形切剖面
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档