JTS基本概念和使用

转x33g5p2x  于2022-05-20 转载在 其他  
字(6.9k)|赞(0)|评价(0)|浏览(586)

简介

  1. JTS是加拿大的 Vivid Solutions公司做的一套开放源码的 Java API。它提供了一套空间数据操作的核心算法。为在兼容OGC标准的空间对象模型中进行基础的几何操作提供2D空间谓词API。

操作

  1. 表示Geometry对象

  2. Geometry类型介绍见另一篇文章：WKT WKB和GeoJSON

  1. package com.alibaba.autonavi;
  3. import com.vividsolutions.jts.geom.Coordinate;
  4. import com.vividsolutions.jts.geom.Geometry;
  5. import com.vividsolutions.jts.geom.GeometryCollection;
  6. import com.vividsolutions.jts.geom.GeometryFactory;
  7. import com.vividsolutions.jts.geom.LineString;
  8. import com.vividsolutions.jts.geom.LinearRing;
  9. import com.vividsolutions.jts.geom.Point;
  10. import com.vividsolutions.jts.geom.Polygon;
  11. import com.vividsolutions.jts.geom.MultiPolygon;
  12. import com.vividsolutions.jts.geom.MultiLineString;
  13. import com.vividsolutions.jts.geom.MultiPoint;
  14. import com.vividsolutions.jts.io.ParseException;
  15. import com.vividsolutions.jts.io.WKTReader;
  17. public class GeometryDemo {
  19.     private GeometryFactory geometryFactory = new GeometryFactory();
  21.     /**
  22.      * create a point
  23.      * @return
  24.      */
  25.     public Point createPoint(){
  26.         Coordinate coord = new Coordinate(109.013388, 32.715519);
  27.         Point point = geometryFactory.createPoint( coord );
  28.         return point;
  29.     }
  30.     
  31.     /**
  32.      * create a point by WKT
  33.      * @return
  34.      * @throws ParseException 
  35.      */
  36.     public Point createPointByWKT() throws ParseException{
  37.         WKTReader reader = new WKTReader( geometryFactory );
  38.         Point point = (Point) reader.read("POINT (109.013388 32.715519)");
  39.         return point;
  40.     }
  41.     
  42.     /**
  43.      * create multiPoint by wkt
  44.      * @return
  45.      */
  46.     public MultiPoint createMulPointByWKT()throws ParseException{
  47.         WKTReader reader = new WKTReader( geometryFactory );
  48.         MultiPoint mpoint = (MultiPoint) reader.read("MULTIPOINT(109.013388 32.715519,119.32488 31.435678)");
  49.         return mpoint;
  50.     }
  51.     /**
  52.      * 
  53.      * create a line
  54.      * @return
  55.      */
  56.     public LineString createLine(){
  57.         Coordinate[] coords  = new Coordinate[] {new Coordinate(2, 2), new Coordinate(2, 2)};
  58.         LineString line = geometryFactory.createLineString(coords);
  59.         return line;
  60.     }
  61.     
  62.     /**
  63.      * create a line by WKT
  64.      * @return
  65.      * @throws ParseException
  66.      */
  67.     public LineString createLineByWKT() throws ParseException{
  68.         WKTReader reader = new WKTReader( geometryFactory );
  69.         LineString line = (LineString) reader.read("LINESTRING(0 0, 2 0)");
  70.         return line;
  71.     }
  72.     
  73.     /**
  74.      * create multiLine 
  75.      * @return
  76.      */
  77.     public MultiLineString createMLine(){
  78.         Coordinate[] coords1  = new Coordinate[] {new Coordinate(2, 2), new Coordinate(2, 2)};
  79.         LineString line1 = geometryFactory.createLineString(coords1);
  80.         Coordinate[] coords2  = new Coordinate[] {new Coordinate(2, 2), new Coordinate(2, 2)};
  81.         LineString line2 = geometryFactory.createLineString(coords2);
  82.         LineString[] lineStrings = new LineString[2];
  83.         lineStrings[0]= line1;
  84.         lineStrings[1] = line2;
  85.         MultiLineString ms = geometryFactory.createMultiLineString(lineStrings);
  86.         return ms;
  87.     }
  88.     
  89.     /**
  90.      * create multiLine by WKT
  91.      * @return
  92.      * @throws ParseException
  93.      */
  94.     public MultiLineString createMLineByWKT()throws ParseException{
  95.         WKTReader reader = new WKTReader( geometryFactory );
  96.         MultiLineString line = (MultiLineString) reader.read("MULTILINESTRING((0 0, 2 0),(1 1,2 2))");
  97.         return line;
  98.     }
  99.     
  100.     /**
  101.      * create a polygon(多边形) by WKT
  102.      * @return
  103.      * @throws ParseException
  104.      */
  105.     public Polygon createPolygonByWKT() throws ParseException{
  106.         WKTReader reader = new WKTReader( geometryFactory );
  107.         Polygon polygon = (Polygon) reader.read("POLYGON((20 10, 30 0, 40 10, 30 20, 20 10))");
  108.         return polygon;
  109.     }
  110.     
  111.     /**
  112.      * create multi polygon by wkt
  113.      * @return
  114.      * @throws ParseException
  115.      */
  116.     public MultiPolygon createMulPolygonByWKT() throws ParseException{
  117.         WKTReader reader = new WKTReader( geometryFactory );
  118.         MultiPolygon mpolygon = (MultiPolygon) reader.read("MULTIPOLYGON(((40 10, 30 0, 40 10, 30 20, 40 10),(30 10, 30 0, 40 10, 30 20, 30 10)))");
  119.         return mpolygon;
  120.     }
  121.     
  122.     /**
  123.      * create GeometryCollection  contain point or multiPoint or line or multiLine or polygon or multiPolygon
  124.      * @return
  125.      * @throws ParseException
  126.      */
  127.     public GeometryCollection createGeoCollect() throws ParseException{
  128.         LineString line = createLine();
  129.         Polygon poly =  createPolygonByWKT();
  130.         Geometry g1 = geometryFactory.createGeometry(line);
  131.         Geometry g2 = geometryFactory.createGeometry(poly);
  132.         Geometry[] garray = new Geometry[]{g1,g2};
  133.         GeometryCollection gc = geometryFactory.createGeometryCollection(garray);
  134.         return gc;
  135.     }
  136.     
  137.     /**
  138.      * create a Circle  创建一个圆，圆心(x,y) 半径RADIUS
  139.      * @param x
  140.      * @param y
  141.      * @param RADIUS
  142.      * @return
  143.      */
  144.     public Polygon createCircle(double x, double y, final double RADIUS){
  145.         final int SIDES = 32;//圆上面的点个数
  146.         Coordinate coords[] = new Coordinate[SIDES+1];
  147.         for( int i = 0; i < SIDES; i++){
  148.             double angle = ((double) i / (double) SIDES) * Math.PI * 2.0;
  149.             double dx = Math.cos( angle ) * RADIUS;
  150.             double dy = Math.sin( angle ) * RADIUS;
  151.             coords[i] = new Coordinate( (double) x + dx, (double) y + dy );
  152.         }
  153.         coords[SIDES] = coords[0];
  154.         LinearRing ring = geometryFactory.createLinearRing( coords );
  155.         Polygon polygon = geometryFactory.createPolygon( ring, null );
  156.         return polygon;
  157.     }
  159.     /**
  160.      * @param args
  161.      * @throws ParseException 
  162.      */
  163.     public static void main(String[] args) throws ParseException {
  164.         GeometryDemo gt = new GeometryDemo();
  165.         Polygon p = gt.createCircle(0, 1, 2);
  166.         //圆上所有的坐标(32个)
  167.         Coordinate coords[] = p.getCoordinates();
  168.         for(Coordinate coord:coords){
  169.             System.out.println(coord.x+","+coord.y);
  170.         }
  171.     }
  172. }

  1. Geometry之间的关系

  2. 支持的空间操作包括

相等(Equals)：几何形状拓扑上相等。
脱节(Disjoint)：几何形状没有共有的点。
相交(Intersects)：几何形状至少有一个共有点（区别于脱节）
接触(Touches)：几何形状有至少一个公共的边界点，但是没有内部点。
交叉(Crosses)：几何形状共享一些但不是所有的内部点。
内含(Within)：几何形状A的线都在几何形状B内部。
包含(Contains)：几何形状B的线都在几何形状A内部（区别于内含）
重叠(Overlaps)：几何形状共享一部分但不是所有的公共点，而且相交处有他们自己相同的区域。
  1. package com.alibaba.autonavi;
  3. import com.vividsolutions.jts.geom.*;
  4. import com.vividsolutions.jts.io.ParseException;
  5. import com.vividsolutions.jts.io.WKTReader;
  7. /**
  8.  * gemotry之间的关系
  9.  *
  10.  */
  11. public class GeometryRelated {
  13.     private GeometryFactory geometryFactory = new GeometryFactory();
  14.     
  15.     /**
  16.      *  两个几何对象是否是重叠的
  17.      * @return
  18.      * @throws ParseException
  19.      */
  20.     public boolean equalsGeo() throws ParseException{
  21.         WKTReader reader = new WKTReader( geometryFactory );
  22.         LineString geometry1 = (LineString) reader.read("LINESTRING(0 0, 2 0, 5 0)");
  23.         LineString geometry2 = (LineString) reader.read("LINESTRING(5 0, 0 0)");
  24.         return geometry1.equals(geometry2);//true
  25.     }
  26.     
  27.     /**
  28.      * 几何对象没有交点(相邻)
  29.      * @return
  30.      * @throws ParseException
  31.      */
  32.     public boolean disjointGeo() throws ParseException{
  33.         WKTReader reader = new WKTReader( geometryFactory );
  34.         LineString geometry1 = (LineString) reader.read("LINESTRING(0 0, 2 0, 5 0)");
  35.         LineString geometry2 = (LineString) reader.read("LINESTRING(0 1, 0 2)");
  36.         return geometry1.disjoint(geometry2);
  37.     }
  38.     
  39.     /**
  40.      * 至少一个公共点(相交)
  41.      * @return
  42.      * @throws ParseException
  43.      */
  44.     public boolean intersectsGeo() throws ParseException{
  45.         WKTReader reader = new WKTReader( geometryFactory );
  46.         LineString geometry1 = (LineString) reader.read("LINESTRING(0 0, 2 0, 5 0)");
  47.         LineString geometry2 = (LineString) reader.read("LINESTRING(0 0, 0 2)");
  48.         Geometry interPoint = geometry1.intersection(geometry2);//相交点
  49.         System.out.println(interPoint.toText());//输出 POINT (0 0)
  50.         return geometry1.intersects(geometry2);
  51.     }
  53.     /**
  54.      * 判断以x,y为坐标的点point(x,y)是否在geometry表示的Polygon中
  55.      * @param x
  56.      * @param y
  57.      * @param geometry wkt格式
  58.      * @return
  59.      */
  60.     public boolean withinGeo(double x,double y,String geometry) throws ParseException {
  62.         Coordinate coord = new Coordinate(x,y);
  63.         Point point = geometryFactory.createPoint( coord );
  65.         WKTReader reader = new WKTReader( geometryFactory );
  66.         Polygon polygon = (Polygon) reader.read(geometry);
  67.         return point.within(polygon);
  68.     }
  69.     /**
  70.      * @param args
  71.      * @throws ParseException 
  72.      */
  73.     public static void main(String[] args) throws ParseException {
  74.         GeometryRelated gr = new GeometryRelated();
  75.         System.out.println(gr.equalsGeo());
  76.         System.out.println(gr.disjointGeo());
  77.         System.out.println(gr.intersectsGeo());
  78.         System.out.println(gr.withinGeo(5,5,"POLYGON((0 0, 10 0, 10 10, 0 10,0 0))"));
  79.     }
  81. }

版权说明 : 本文为转载文章, 版权归原作者所有 版权申明

原文链接 : https://lebron.blog.csdn.net/article/details/124874405

内容来源于网络，如有侵权，请联系作者删除！

Java开发语言JTS图形算法

相关文章

    查看更多

    热门标签

    更多
    JavaquerypythonNode开发语言requestUtil数据库Table后端算法LoggerMessageElementParser

    最新文章

    更多
    • 蜀ICP备13028337号-1 大数据知识库 https://www.saoniuhuo.com © All rights reserved
    • 本站内容来源互联网,如果侵犯您的权益请联系我们删除, 联系方式：448109455@qq.com