GIS算法：JAVA拓扑套件JTS

常用可以用于GIS数据处理和空间计算的java包有geotool和jts。

相对来说，geotool功能更全面，还可以用于数据转换、瓦片地图发布、栅格影像分析等，jts只能进行基本的数据处理和空间计算。

但大多数情况下jts就完全够用了。

geotool的官网：https://www.geotools.org/

本例只讲jts的用法：

maven依赖：

<dependency><groupId>com.vividsolutions</groupId><artifactId>jts</artifactId><version>1.13</version>
</dependency>

import com.vividsolutions.jts.geom.Coordinate;
import com.vividsolutions.jts.geom.Envelope;
import com.vividsolutions.jts.geom.Geometry;
import com.vividsolutions.jts.geom.GeometryFactory;
import com.vividsolutions.jts.index.strtree.STRtree;import java.util.List;public class RTreeDemo {public static void main(String[] args){//声明STRtreeSTRtree strTree = new STRtree();GeometryFactory geometryFactory = new GeometryFactory();for (int i=0;i<50;i++){for (int j=0;j<50;j++){//新建一个点Geometry p=geometryFactory.createPoint(new Coordinate(i,j));//以点为中心，取一个半径为1.0的圆，插入RTree
//insert(Envelope itemEnv, Object item)strTree.insert(p.buffer(1.0).getEnvelopeInternal(),p.buffer(1.0));}}//构建RTreestrTree.build();//查询点[1,1]落在哪些圆中List querys=strTree.query(geometryFactory.createPoint(new Coordinate(1,1)).getEnvelopeInternal());for (Object obj:querys) {System.out.println(obj);}}
}

线裁切面：如图所示，用红色线去裁切蓝色面，结果会生成3个面。

import com.vividsolutions.jts.geom.Geometry;
import com.vividsolutions.jts.io.ParseException;
import com.vividsolutions.jts.io.WKTReader;
import com.vividsolutions.jts.operation.polygonize.Polygonizer;import java.util.ArrayList;
import java.util.Collection;
import java.util.List;public class PolygonDemo {public static void main(String[] args){// wkt工具，将wkt文本转为geometry对象WKTReader wktReader=new WKTReader();try {//读被裁切面Geometry polygon=wktReader.read("POLYGON ((220 350, 400 440, 635 249, 380 80, 174 164, 179 265, 220 350))");//读裁切线Geometry polyline=wktReader.read("LINESTRING (570 400, 392 315, 299 215, 430 140, 530 240, 450 360, 460 480)");//取面的边线Geometry boundary=polygon.getBoundary();//将裁切线与面的边线联合，交点会被打断polyline=polyline.union(boundary);List<Geometry> clipPolygon = new ArrayList<>();List<Geometry> lineList=new ArrayList<>();for(int i=0;i<polyline.getNumGeometries();i++){lineList.add(polyline.getGeometryN(i));}// 构造面生成器Polygonizer p = new Polygonizer();p.add(lineList);//取构面结果Collection<Geometry> polys = p.getPolygons();//取buffer，以免因精度损失，遗漏构面结果Geometry buffer=polygon.buffer(1);for(Geometry geometry:polys){//如果包含在buffer中，则添加if(buffer.contains(geometry)){clipPolygon.add(geometry);}}for (Geometry presult:clipPolygon) {System.out.println(presult);}} catch (ParseException e) {e.printStackTrace();}}
}