geotools获取给定点的DEM高程值

概述

本文讲述结合geotools实现获取给定点的DEM(tif格式)高程值。

需求描述

1、在web端绘制一条曲线; 2、获取各节点处的高程值; 3、根据高程值绘制高程堆积图。

实现效果

实现思路及代码

1、绘制曲线,将绘制的结果传给后台

            var draw = new ol.interaction.Draw({
                source: source,
                type: "LineString",
                freehand: false//是否手绘
            });
            map.addInteraction(draw);
            draw.setActive(false);
            draw.on("drawend", function (evt) {
                draw.setActive(false);
                var coords = evt.feature.getGeometry().getCoordinates(),
                    points = [];
                for(var i=0;i<coords.length;i++){
                    points.push(coords[i].join(","));
                }
                $.get("dem/query", {points: points.join(";")}, function (result) {
                   console.info(result);
                })
            });

2、后台解析数据并获取点的高程数据

        String points_data = request.getParameter("points");
        String demPath = "/Users/lzugis/Documents/ncdata/bj_dem.tif";
        File file = new File(demPath);
        GeoTiffReader tifReader = new GeoTiffReader(file);
        GridCoverage2D coverage = tifReader.read(null);

        CoordinateReferenceSystem crs = coverage.getCoordinateReferenceSystem2D();
        String[] points = points_data.split(";");

        List list = new ArrayList();

        for(int i=0;i<points.length;i++) {
            String strLonlat = points[i];
            String[] strLonlats = strLonlat.split(",");

            double lon = Double.parseDouble(strLonlats[0]),
                    lat = Double.parseDouble(strLonlats[1]);

            DirectPosition position = new DirectPosition2D(crs, lon, lat);
            int[] results = (int[]) coverage.evaluate(position);
            results = coverage.evaluate(position, results);
            Map map = new HashMap();
            map.put("lon", lon);
            map.put("lat", lon);
            map.put("dem", results[0]);
            list.add(JSONObject.toJSONString(map));
            response.getWriter().println(JSONArray.toJSONString(list));
        }

3、前段接受数据并展示

                var dom = document.getElementById("chart");
                var myChart = echarts.init(dom);
                var option = {
                    title: {
                        text: '高程图'
                    },
                    backgroundColor:"rgba(255, 255, 255, .8)",
                    grid: {
                        left: '5%',
                        right: '4%',
                        bottom: '1%',
                        containLabel: true
                    },
                    xAxis : [
                        {
                            show:false,
                            type : 'category',
                            boundaryGap : false,
                            data : ["点1", "点2", "点3", "点4"]
                        }
                    ],
                    yAxis : [
                        {
                            type : 'value',
                            offset:15
                        }
                    ],
                    series : [
                        {
                            name: 'DEM',
                            type: 'line',
                            stack: '总量',
                            label: {
                                normal: {
                                    show: true,
                                    position: 'top'
                                }
                            },
                            areaStyle: {normal: {}},
                            data: [714, 1000, 1063, 864]
                        }
                    ]
                };
                myChart.setOption(option);
说明:

1、dem数据来源于地理空间数据云

2、展示用echart实现的。

本文参与腾讯云自媒体分享计划,欢迎正在阅读的你也加入,一起分享。

发表于

我来说两句

0 条评论
登录 后参与评论

相关文章

来自专栏码匠的流水账

FluxSink实例及解析

reactor-core-3.1.3.RELEASE-sources.jar!/reactor/core/publisher/FluxSink.java

1272
来自专栏曾大稳的博客

Android ClassLoader流程解读并简单方式实现热更新

ClassLoader在启动Activity的时候会调用loadClass方法,我们就从这里入手:

2632
来自专栏菩提树下的杨过

温故而知新:设计模式之原型模式(Prototype)

原型模式个人以为最适合的场景:参照现有的某一个对象实例,快速得到多个完整的实例副本。(通常是深拷贝的副本) 深拷贝在c#中实现的最简单方式莫过于通过反序列化得到...

2355
来自专栏码匠的流水账

聊聊storm的IEventLogger

storm-2.0.0/storm-client/src/jvm/org/apache/storm/metric/IEventLogger.java

1343
来自专栏码匠的流水账

聊聊kafka 0.8 ConsumerFetcherManager的MaxLag指标

本文主要研究一下kafka0.8.2.2版本中ConsumerFetcherManager的MaxLag指标的统计。

861
来自专栏码匠的流水账

聊聊flink的RichParallelSourceFunction

flink-streaming-java_2.11-1.6.2-sources.jar!/org/apache/flink/streaming/api/func...

1991
来自专栏高性能服务器开发

(五)sparkline微线图

sparkline这个单词,我第一次看的时候,也不知道这什么意思啊,以前根本没听过啊,但是这真真实实的出现在了redis的代码中了,刚刚开始以为这也是属于普通...

37012
来自专栏前端新视界

由移动端级联选择器所引发的对于数据结构的思考

GitHub:https://github.com/nzbin/Framework7-CityPicker Demo:https://nzbin.githu...

3808
来自专栏码匠的流水账

bloomfilter的简单实现

布隆过滤器(英语:Bloom Filter)是1970年由布隆提出的,可以用于检索一个元素是否在一个集合中。

681
来自专栏码匠的流水账

聊聊storm的IEventLogger

storm-2.0.0/storm-client/src/jvm/org/apache/storm/metric/IEventLogger.java

842

扫码关注云+社区