首页 > 编程语言 >Java中使用JTS对空间几何计算(读取WKT、距离、点在面内、长度、面积、相交等)

Java中使用JTS对空间几何计算(读取WKT、距离、点在面内、长度、面积、相交等)

时间:2023-08-29 14:48:09浏览次数:56  
标签:面内 Java 10 System JTS 20 println Coordinate new

场景

基于GIS相关的集成系统,需要对空间数据做一些判断处理。比如读取WKT数据、点到点、点到线、点到面的距离,

线的长度、面的面积、点是否在面内等处理。

JTS

(Java Topology Suite) Java拓扑套件,是Java的处理地理数据的API。

github地址:

https://github.com/locationtech/jts

API文档地址:

https://locationtech.github.io/jts/javadoc/

Maven中央仓库地址:

 https://mvnrepository.com/artifact/org.locationtech.jts/jts-core

特点

实现了OGC关于简单要素SQL查询规范定义的空间数据模型
一个完整的、一致的、基本的二维空间算法的实现,包括二元运算(例如touch和overlap)和空间分析方法(例如intersection和buffer)
一个显示的精确模型,用算法优雅的解决导致dimensional collapse(尺度坍塌–专业名词不知道对不对,暂时这样译)的情况。
健壮的实现了关键计算几何操作
提供著名文本格式的I/O接口
JTS是完全100%由Java写的

JTS支持一套完整的二元谓词操作。二元谓词方法将两个几何图形作为参数,
返回一个布尔值来表示几何图形是否有指定的空间关系。它支持的空间关系有:
相等(equals)、分离(disjoint)、相交(intersect)、相接(touches)、
交叉(crosses)、包含于(within)、包含(contains)、覆盖/覆盖于(overlaps)。
同时,也支持一般的关系(relate)操作符。
relate可以被用来确定维度扩展的九交模型(DE-9IM),它可以完全的描述两个几何图形的关系。

空间数据模型

JTS提供了以下空间数据模型

 

 

图形可视化WKT数据

在jts的bin下的testbuilder.bat,双击运行

 

 

即可运行WKT数据可视化界面

 

 

可以在页面上绘制图形并下方生成wkt数据,以及输入wkt数据,点击右边按钮,图形化显示。


实现

1、项目中引入jts的依赖

<!-- https://mvnrepository.com/artifact/org.locationtech.jts/jts-core -->
        <dependency>
            <groupId>org.locationtech.jts</groupId>
            <artifactId>jts-core</artifactId>
            <version>1.18.2</version>
        </dependency>

 

2、从WKT字符串中读取几何图形,读取点、线、面

//read a geometry from a WKT string (using the default geometry factory)
        //从WKT字符串读取几何图形
        Geometry g1 = null;
        try {
            //读取线
            //g1 = new WKTReader().read("LINESTRING (0 0, 10 10, 20 20)");
            //读取点
            //g1 = new WKTReader().read("POINT (2 2)");
            //读取面
            g1 = new WKTReader().read("POLYGON((40 100, 40 20, 120 20, 120 100, 40 100))");
        } catch (ParseException e) {
            e.printStackTrace();
        }
//        Arrays.stream(g1.getCoordinates()).forEach(point -> System.out.println("x:"+point.x+" y:"+point.y));
        //输出结果:
//        x:0.0 y:0.0
//        x:10.0 y:10.0
//        x:20.0 y:20.0

        //Arrays.stream(g1.getCoordinates()).forEach(point ->System.out.println("x:"+point.x+" y:"+point.y));
        //输出:x:2.0 y:2.0

        //Arrays.stream(g1.getCoordinates()).forEach(point ->System.out.println("x:"+point.x+" y:"+point.y));
        //输出:
//        x:40.0 y:100.0
//        x:40.0 y:20.0
//        x:120.0 y:20.0
//        x:120.0 y:100.0
//        x:40.0 y:100.0

 

3、创建点、线

//创建点
        Point point = new GeometryFactory().createPoint(new Coordinate(1, 1));

        // create a geometry by specifying the coordinates directly
        //通过指定坐标创建线
        Coordinate[] coordinates = new Coordinate[]{new Coordinate(0, 0),
                new Coordinate(10, 10), new Coordinate(20, 20)};
        // use the default factory, which gives full double-precision
        Geometry g2 = new GeometryFactory().createLineString(coordinates);
        //System.out.println("Geometry 2: " + g2);
        //输出结果:Geometry 2: LINESTRING (0 0, 10 10, 20 20)

 

4、计算点是否在线上、点是否在面内

//创建点
        Point point = new GeometryFactory().createPoint(new Coordinate(1, 1));
        //输出结果:POINT (1 1)
        //计算点是否在线上
        //System.out.println(g1.contains(point));
        //输出结果:true

        //计算点是否在面内
        Point point2 = new GeometryFactory().createPoint(new Coordinate(70, 70));
        //System.out.println(g1.contains(point2));
        //输出结果: true
        Point point3 = new GeometryFactory().createPoint(new Coordinate(20, 10));
        //System.out.println(g1.contains(point3));
        //输出结果: false

 

5、计算两个几何图形的交点

        // compute the intersection of the two geometries
        //计算两个几何图形的交点
        Geometry g3 = g1.intersection(g2);
        //System.out.println("G1 intersection G2: " + g3);
        //输出结果:G1 intersection G2: MULTILINESTRING ((0 0, 10 10), (10 10, 20 20))

 

6、创建一个MultiPoint多点

// create a factory using default values (e.g. floating precision)
        //创建一个MultiPoint多点
        GeometryFactory fact = new GeometryFactory();

//        Point p1 = fact.createPoint(new Coordinate(0,0));
//        System.out.println(p1);
//
//        Point p2 = fact.createPoint(new Coordinate(1,1));
//        System.out.println(p2);
//
//        MultiPoint mpt = fact.createMultiPointFromCoords(new Coordinate[] { new Coordinate(0,0), new Coordinate(1,1) } );
//        System.out.println(mpt);

        //输出结果:
//        POINT (0 0)
//        POINT (1 1)
//        MULTIPOINT ((0 0), (1 1))

 

7、创建闭合线LinearRing

 //创建闭合线-LinearRing
        LinearRing lr = new GeometryFactory().createLinearRing(new Coordinate[]{new Coordinate(0, 0), new Coordinate(0, 10), new Coordinate(10, 10), new Coordinate(10, 0), new Coordinate(0, 0)});
        //System.out.println(lr);
        //输出结果:LINEARRING (0 0, 0 10, 10 10, 10 0, 0 0)

 

8、创建几何集合列表

//创建几何集合列表
        Geometry[] garray = new Geometry[]{g1,g2};
        GeometryCollection gc = fact.createGeometryCollection(garray);
        //System.out.println(gc.toString());
        //输出结果:GEOMETRYCOLLECTION (POLYGON ((40 100, 40 20, 120 20, 120 100, 40 100)), LINESTRING (0 0, 10 10, 20 20))

 

9、几何关系判断

//几何关系判断,是否相交intersection
        //其他方法类似
//        相等(Equals): 几何形状拓扑上相等。
//        不相交(Disjoint): 几何形状没有共有的点。
//        相交(Intersects): 几何形状至少有一个共有点(区别于脱节)
//        接触(Touches): 几何形状有至少一个公共的边界点,但是没有内部点。
//        交叉(Crosses): 几何形状共享一些但不是所有的内部点。
//        内含(Within): 几何形状A的线都在几何形状B内部。
//        包含(Contains): 几何形状B的线都在几何形状A内部(区别于内含)
//        重叠(Overlaps): 几何形状共享一部分但不是所有的公共点,而且相交处有他们自己相同的区域。

        WKTReader reader = new WKTReader(fact);
        LineString geometry1 = null;
        try {
            geometry1 = (LineString) reader.read("LINESTRING(0 0, 2 0, 5 0)");
        } catch (ParseException e) {
            e.printStackTrace();
        }
        LineString geometry2 = null;
        try {
            geometry2 = (LineString) reader.read("LINESTRING(0 0, 0 2)");
        } catch (ParseException e) {
            e.printStackTrace();
        }
        Geometry interPoint = geometry1.intersection(geometry2);//相交点
        //System.out.println(interPoint.toText());
        //输出结果: POINT (0 0)

 

10、计算距离distance

//计算距离distance
        Point p1 = fact.createPoint(new Coordinate(0,0));
        //System.out.println(p1);

        Point p2 = fact.createPoint(new Coordinate(3,4));
        ///System.out.println(p2);

        //System.out.println(p1.distance(p2));
        //输出结果
//        POINT (0 0)
//        POINT (3 4)
//        5.0

 

11、计算长度和面积

Geometry g5 = null;
        Geometry g6 = null;
        try {
            //读取面
            g5 = new WKTReader().read("POLYGON((40 100, 40 20, 120 20, 120 100, 40 100))");
            g6 = new WKTReader().read("LINESTRING(0 0, 0 2)");
            //计算面积getArea()
            //System.out.println(g5.getArea());
            //输出结果:6400.0
            //计算长度getLength()
            //System.out.println(g6.getLength());
            //输出结果:2.0
        } catch (ParseException e) {
            e.printStackTrace();
        }

 

12、求点到线、点到面的最近距离

GeometryFactory gf = new GeometryFactory();
        WKTReader reader2 = new WKTReader(gf);
        Geometry line2 = null;
        Geometry g7 = null;
        try {
            line2 = reader2.read("LINESTRING(0 0, 10 0, 10 10, 20 10)");
            g7 = new WKTReader().read("POLYGON((40 100, 40 20, 120 20, 120 100, 40 100))");
        } catch (ParseException e) {
            e.printStackTrace();
        }
        Coordinate c = new Coordinate(5, 5);
        PointPairDistance ppd = new PointPairDistance();
        //求点到线的最近距离
        //DistanceToPoint.computeDistance(line2,c,ppd);
        //输出结果:5.0
        //求点到面的最近距离
        DistanceToPoint.computeDistance(g7,c,ppd);
        System.out.println(ppd.getDistance());
        //输出结果38.07886552931954

 

标签:面内,Java,10,System,JTS,20,println,Coordinate,new
From: https://www.cnblogs.com/Fooo/p/17664702.html

相关文章

  • [java基础知识复习] Java基础知识总结分享一
    写代码:1,明确需求。我要做什么?2,分析思路。我要怎么做?1,2,3。3,确定步骤。每一个思路部分用到哪些语句,方法,和对象。4,代码实现。用具体的java语言代码把思路体现出来。学习新技术的四点:1,该技术是什么?2,该技术有什么特点(使用注意):3,该技术怎么使用。demo4,该技术什么时候用?test。————......
  • jts learning
    JTS简介JTS提供了一套操作几何向量的java类库。早期版本com.vividsolutions,已废弃不在维护。现在版本com.locationtech.jts由eclipse开源基金会托管使用说明入门指导GIS开发入门指导jts-core核心库使用说明jts-core核心类库使用说明具体结合示例代码详细介绍JTS下面核心......
  • redis分布式锁,setnx+lua脚本的java实现
    1前言在现在工作中,为保障服务的高可用,应对单点故障、负载量过大等单机部署带来的问题,生产环境常用多机部署。为解决多机房部署导致的数据不一致问题,我们常会选择用分布式锁。目前其他比较常见的实现方案我列举在下面:基于缓存实现分布式锁(本文主要使用redis实现)基于数据库实......
  • 怎么样去学习javaweb
    据Java官方记载,只有JavaSE和JavaEE这类表述,似乎并没有JavaWeb这么一个说法,所以,我更愿意认为,JavaWeb是各位Java前辈们对于Java在Web方向应用的相关技术栈,正如现在的PythonWeb,golangweb等等,只不过是通过某种语言来实现web应用而已。让我们先看看百度百科对于JavaWeb的定义,JavaW......
  • redis分布式锁,setnx+lua脚本的java实现 | 京东物流技术团队
    1前言在现在工作中,为保障服务的高可用,应对单点故障、负载量过大等单机部署带来的问题,生产环境常用多机部署。为解决多机房部署导致的数据不一致问题,我们常会选择用分布式锁。目前其他比较常见的实现方案我列举在下面:基于缓存实现分布式锁(本文主要使用redis实现)基于数据库实现分布......
  • JTS-IntersectionMatrix 使用说明
    参考:https://blog.csdn.net/weixin_40294332/article/details/124124928参考2:https://vimsky.com/examples/detail/java-method-com.vividsolutions.jts.geom.IntersectionMatrix.set.html......
  • Java 15 JSTL实现登录退出
     jstl.jsp<%@pagecontentType="text/html;charset=UTF-8"language="java"%><%@taglibprefix="c"uri="http://java.sun.com/jsp/jstl/core"%><%--if--%><%@taglibprefix="fmt"uri=&......
  • Java 14 EL表达式
     eljstl.jsp<%@pagecontentType="text/html;charset=UTF-8"language="java"%><html><head><title>Title</title></head><body><h1>首页</h1><%--指令包含<%@includefile=......
  • JAVA SE基础《十》 ---- 面对对象高级一
    目录一、static  1、static修饰成员变量2、static修饰成员变量的应用场景:在线人数统计3、static修饰成员方法4、static修饰成员方法的应用场景:工具类5、static的注意事项6、static的应用知识:代码块7、static的应用知识:单例设计模式二、面向......
  • Java代码审计之目录穿越
    一、目录穿越漏洞1、什么是目录穿越所谓的目录穿越指利用操作系统中的文件系统对目录的表示。在文件系统路径中,".."表示上一级目录,当你使用"../"时,你正在引用当前目录的上一级目录。如果你使用"../../",你实际上在两次".."的基础上,再次引用上一级目录,从而返回到上两级目录。......