esri的geometry-api学习

<dependency>
            <groupId>com.esri.geometry</groupId>
            <artifactId>esri-geometry-api</artifactId>
            <version>2.2.3</version>
        </dependency>
           

1、2个geojson是否相交

import com.alibaba.fastjson.JSONObject;
import com.esri.core.geometry.*;
import org.geotools.data.shapefile.ShapefileDataStore;
import org.geotools.data.simple.SimpleFeatureIterator;
import org.geotools.data.simple.SimpleFeatureSource;
import org.geotools.geojson.feature.FeatureJSON;
import org.opengis.feature.simple.SimpleFeature;

import java.io.File;
import java.io.StringWriter;
import java.nio.charset.Charset;
import java.util.ArrayList;
import java.util.List;

    		File[] files = new File("E:\\test").listFiles();
		for (File shapeFile : files) {
			if (shapeFile.getName().endsWith(".shp")) {
				// 使用GeoTools读取ShapeFile文件
				ShapefileDataStore store = new ShapefileDataStore(shapeFile.toURI().toURL());
				//设置编码
				Charset charset = Charset.forName("GBK");
				store.setCharset(charset);
				SimpleFeatureSource sfSource = store.getFeatureSource();
				SimpleFeatureIterator sfIter = sfSource.getFeatures().features();
				// 从ShapeFile文件中遍历每一个Feature，然后将Feature转为GeoJSON字符串
				List<String> ls = new ArrayList<>();
				while (sfIter.hasNext()) {
					SimpleFeature feature = sfIter.next();
					// Feature转GeoJSON
					FeatureJSON fjson = new FeatureJSON();
					StringWriter writer = new StringWriter();
					fjson.writeFeature(feature, writer);
					String sjson = writer.toString();
					JSONObject featureStr = JSONObject.parseObject(sjson);
					String geojson = featureStr.get("geometry").toString();
					ls.add(geojson);
					MapGeometry mapGeometry = GeometryEngine.geoJsonToGeometry(geojson,
							WktImportFlags.wktImportDefaults, com.esri.core.geometry.Geometry.Type.Polygon);
					com.esri.core.geometry.Geometry polygon = mapGeometry.getGeometry();
//					double l = GeometryEngine.distance(polygon,new Point(108.41, 30.81),SpatialReference.create(4326));
					boolean b = isWithin(108.49, 30.77,polygon,0.2);
					System.out.println("l = " + b);
				}
				System.out.println("数据导入完成，共耗时" + (System.currentTimeMillis() - start) + "ms");
			}
		}

	// Disjoint	判断两个图形间是否没有相同点
	public static boolean isWithin(double lon, double lat, com.esri.core.geometry.Geometry p, double dis) throws Exception {//dis=20000大约0.2度
		String geojson = "{\"type\":\"Polygon\",\"coordinates\":[[["+(lon-dis)+","+(lat-dis)+"],["+(lon+dis)+","+(lat+dis)+"],[" +
				""+(lon-dis)+","+(lat+dis)+"],["+(lon+dis)+","+(lat-dis)+"]]]}";
		MapGeometry mapGeometry = GeometryEngine.geoJsonToGeometry(geojson,
				WktImportFlags.wktImportDefaults, com.esri.core.geometry.Geometry.Type.Polygon);
		return GeometryEngine.disjoint(mapGeometry.getGeometry(), p, SpatialReference.create(4326));
	}
           

相等(Equals)：	几何形状拓扑上相等。
脱节(Disjoint)：	å‡ ä½•å½¢çŠ¶æ²¡æœ‰å ±æœ‰çš„ç‚¹ã€‚
相交(Intersects)：	å‡ ä½•å½¢çŠ¶è‡³å°‘æœ‰ä¸€ä¸ªå ±æœ‰ç‚¹ï¼ˆåŒºåˆ«äºŽè„±èŠ‚ï¼‰
接触(Touches)：	å‡ ä½•å½¢çŠ¶æœ‰è‡³å°‘ä¸€ä¸ªå ¬å ±çš„è¾¹ç•Œç‚¹ï¼Œä½†æ˜¯æ²¡æœ‰å† éƒ¨ç‚¹ã€‚
交叉(Crosses)：	å‡ ä½•å½¢çŠ¶å ±äº«ä¸€äº›ä½†ä¸æ˜¯æ‰€æœ‰çš„å† éƒ¨ç‚¹ã€‚
å† å«(Within)：	几何形状A的线都在几何形状Bå† éƒ¨ã€‚
åŒ å«(Contains)：	几何形状B的线都在几何形状Aå† éƒ¨ï¼ˆåŒºåˆ«äºŽå† å«ï¼‰
重叠(Overlaps)：	å‡ ä½•å½¢çŠ¶å ±äº«ä¸€éƒ¨åˆ†ä½†ä¸æ˜¯æ‰€æœ‰çš„å ¬å ±ç‚¹ï¼Œè€Œä¸”ç›¸äº¤å¤„æœ‰ä»–ä»¬è‡ªå·±ç›¸åŒçš„åŒºåŸŸã€‚

缓冲区分析（Buffer）	åŒ å«æ‰€æœ‰çš„ç‚¹åœ¨ä¸€ä¸ªæŒ‡å®šè·ç¦»å† çš„å¤šè¾¹å½¢å’Œå¤šå¤šè¾¹å½¢
凸壳分析（ConvexHull）	åŒ å«å‡ ä½•å½¢ä½“çš„æ‰€æœ‰ç‚¹çš„æœ€å°å‡¸å£³å¤šè¾¹å½¢ï¼ˆå¤–åŒ å¤šè¾¹å½¢ï¼‰
交叉分析（Intersection）	交叉操作就是多边形ABä¸­æ‰€æœ‰å ±åŒç‚¹çš„é›†åˆã€‚
联合分析（Union）	AB的联合操作就是AB所有点的集合。
差异分析（Difference）	AB形状的差异分析就是A里有B里没有的所有点的集合。
对称差异分析（SymDifference）	AB形状的对称差异分析就是位于Aä¸­æˆ–è€ B中但不同时在AB中的所有点的集合

参考：https://blog.csdn.net/objectively/article/details/35219605?utm_source=blogxgwz3