From e2a7b8aff862419dc4b868fc7bb2759012be6409 Mon Sep 17 00:00:00 2001 From: liuchengqian Date: Thu, 13 Apr 2023 17:21:33 +0800 Subject: [PATCH] =?UTF-8?q?=E7=81=AB=E7=82=B9=E9=87=8D=E6=9E=84-=E5=AE=8C?= =?UTF-8?q?=E5=96=84=E5=9C=B0=E7=89=A9=E7=B1=BB=E5=9E=8B=E7=99=BE=E5=88=86?= =?UTF-8?q?=E6=AF=94?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- pom.xml | 90 +++++++++++++ .../xkrs/utilsnew/BufferAnalysisUtils.java | 123 ++++++++++++++++++ .../java/com/xkrs/utilsnew/ShpFileUtils.java | 50 +++++++ 3 files changed, 263 insertions(+) create mode 100644 src/main/java/com/xkrs/utilsnew/BufferAnalysisUtils.java create mode 100644 src/main/java/com/xkrs/utilsnew/ShpFileUtils.java diff --git a/pom.xml b/pom.xml index b9a7e15..cdf70de 100644 --- a/pom.xml +++ b/pom.xml @@ -35,6 +35,7 @@ 2.0.4 4.4.3 0.4.0 + 28.2 @@ -255,8 +256,97 @@ 1.13.1 + + org.geotools + gt-xml + ${geotools.version} + + + + org.geotools + gt-shapefile + ${geotools.version} + + + + org.geotools + gt-main + ${geotools.version} + + + jts-core + org.locationtech.jts + + + + + + org.geotools.xsd + gt-xsd-sld + ${geotools.version} + + + + org.geotools + gt-swing + ${geotools.version} + + + jts-core + org.locationtech.jts + + + guava + com.google.guava + + + + + + org.geotools + gt-geojson + ${geotools.version} + + + + org.geotools + gt-epsg-hsql + ${geotools.version} + + + + org.geotools.jdbc + gt-jdbc-postgis + ${geotools.version} + + + + + osgeo + OSGeo Release Repository + https://repo.osgeo.org/repository/release/ + + false + + + true + + + + osgeo-snapshot + OSGeo Snapshot Repository + https://repo.osgeo.org/repository/snapshot/ + + true + + + false + + + + diff --git a/src/main/java/com/xkrs/utilsnew/BufferAnalysisUtils.java b/src/main/java/com/xkrs/utilsnew/BufferAnalysisUtils.java new file mode 100644 index 0000000..01e4472 --- /dev/null +++ b/src/main/java/com/xkrs/utilsnew/BufferAnalysisUtils.java @@ -0,0 +1,123 @@ +package com.xkrs.utilsnew; + +import org.apache.hc.core5.util.TextUtils; +import org.geotools.data.simple.SimpleFeatureCollection; +import org.geotools.data.simple.SimpleFeatureIterator; +import org.geotools.filter.text.cql2.CQL; +import org.geotools.filter.text.cql2.CQLException; +import org.locationtech.jts.geom.*; +import org.locationtech.jts.operation.buffer.BufferParameters; +import org.opengis.feature.simple.SimpleFeature; +import org.opengis.filter.Filter; + +import java.text.DecimalFormat; +import java.util.ArrayList; +import java.util.HashMap; +import java.util.List; +import java.util.Map; + +public class BufferAnalysisUtils { + + + private BufferAnalysisUtils() { + } + + public static void main(String[] args) throws CQLException { + String shapeFilePath = "/Users/liuchengqian/Desktop/data/shandong_land_type_2022.shp"; + SimpleFeatureCollection simpleFeatureCollection = ShpFileUtils.readShpFileFeatureCollection(shapeFilePath); + double longitude = 116.218;//经度 + double latitude = 36.0852;//纬度 + double distance = getBufferDistance(1500L);// buffer 宽度 1.5公里 + GeometryFactory geometryFactory = new GeometryFactory(new PrecisionModel(PrecisionModel.FLOATING), 4326); + Coordinate[] coordinates = new Coordinate[]{new Coordinate(longitude, latitude)}; + Point point = new Point(geometryFactory.getCoordinateSequenceFactory().create(coordinates), geometryFactory); + Geometry pointBuffer = point.buffer(distance, 40, BufferParameters.CAP_ROUND); + String pointBufferWkt = pointBuffer.toText(); + double pointBufferArea = pointBuffer.getArea();//总面积 + // 相交过滤 + Filter intersectsFilter = CQL.toFilter(String.format("INTERSECTS(the_geom, %s)", pointBufferWkt)); + SimpleFeatureCollection filterSimpleFeatureCollection = simpleFeatureCollection.subCollection(intersectsFilter); + //地物类型对应的字段 + String landTypeProperty = "tblxmc"; + String defaultLandType = "其它"; + List targetLandTypeList = new ArrayList<>(); + targetLandTypeList.add("林地"); + targetLandTypeList.add("农作物"); + targetLandTypeList.add("水域"); + targetLandTypeList.add("建筑物"); + targetLandTypeList.add("水产养殖设施农用地"); +// targetLandTypeList.add("其他土地"); + Map targetLandTypeAreaMap = new HashMap<>(); + for (String targetLandType : targetLandTypeList) { + targetLandTypeAreaMap.put(targetLandType, 0D); + } + SimpleFeatureIterator featureIterator = filterSimpleFeatureCollection.features(); + while (featureIterator.hasNext()) { + try { + SimpleFeature simpleFeature = featureIterator.next(); + if (simpleFeature == null) { + System.out.println("simpleFeature == null"); + continue; + } + Geometry geometry = (Geometry) simpleFeature.getDefaultGeometry(); + if (geometry == null) { + System.out.println("geometry == null"); + continue; + } + String landType = simpleFeature.getProperty(landTypeProperty).getValue().toString().trim(); + if (TextUtils.isEmpty(landType) || (!targetLandTypeAreaMap.containsKey(landType))) { + System.out.println("!targetLandTypeMap.containsKey(landType)， landType = " + landType); + continue; + } + // 检查多边形是否存在自相交或重叠的情况，如果有则修复 + Geometry geometryBuffer = geometry.buffer(0D); + if (!pointBuffer.intersects(geometryBuffer)) { + System.out.println("!pointBuffer.intersects(geometryBuffer)"); + continue; + } + Geometry intersectGeometry = pointBuffer.intersection(geometryBuffer); + if (intersectGeometry == null) { + System.out.println("intersectGeometry == null"); + continue; + } + double intersectArea = intersectGeometry.getArea(); + if (intersectArea < 0D) { + System.out.println("intersectArea < 0"); + continue; + } + targetLandTypeAreaMap.put(landType, targetLandTypeAreaMap.get(landType) + intersectArea); + } catch (Exception e) { + e.printStackTrace(); + } + } + double targetLandAreaSum = 0D; + StringBuilder builder = new StringBuilder(); + for (String targetLandType : targetLandTypeList) { + double currentArea = targetLandTypeAreaMap.get(targetLandType); + targetLandAreaSum += currentArea; + builder.append(targetLandType); + builder.append("="); + builder.append(formatDoublePer(currentArea / pointBufferArea)); + builder.append(", "); + } + builder.append(defaultLandType); + builder.append("="); + builder.append(formatDoublePer((pointBufferArea - targetLandAreaSum) / pointBufferArea)); + String result = builder.toString(); + System.out.println(result); + } + + + public static final int RADIUS = 6377830; + + private static double getBufferDistance(long meter) { + return meter * 180.D / Math.PI / RADIUS; + } + + private static final DecimalFormat decimalFormat = new DecimalFormat("0.00"); + + private static String formatDoublePer(double value) { + return decimalFormat.format(Math.abs(value) * 100D) + "%"; + } + +} diff --git a/src/main/java/com/xkrs/utilsnew/ShpFileUtils.java b/src/main/java/com/xkrs/utilsnew/ShpFileUtils.java new file mode 100644 index 0000000..717b574 --- /dev/null +++ b/src/main/java/com/xkrs/utilsnew/ShpFileUtils.java @@ -0,0 +1,50 @@ +package com.xkrs.utilsnew; + +import org.geotools.data.FeatureSource; +import org.geotools.data.Query; +import org.geotools.data.shapefile.ShapefileDataStore; +import org.geotools.data.simple.SimpleFeatureCollection; +import org.geotools.geometry.jts.ReferencedEnvelope; +import org.opengis.feature.type.AttributeDescriptor; +import org.opengis.feature.type.GeometryType; +import org.opengis.referencing.crs.CoordinateReferenceSystem; + +import java.io.File; +import java.io.IOException; +import java.nio.charset.StandardCharsets; +import java.util.List; + +public class ShpFileUtils { + + private ShpFileUtils() { + } + + public static SimpleFeatureCollection readShpFileFeatureCollection(String shpPath) { + File shpFile = new File(shpPath); + SimpleFeatureCollection simpleFeatureCollection = null; + try { + ShapefileDataStore shapefileDataStore = new ShapefileDataStore(shpFile.toURI().toURL()); + // 设置编码,防止属性的中文字符出现乱码 + shapefileDataStore.setCharset(StandardCharsets.UTF_8); + // 这个typeNamae不传递，默认是文件名称 + FeatureSource featuresource = shapefileDataStore.getFeatureSource(shapefileDataStore.getTypeNames()[0]); + // 读取bbox + ReferencedEnvelope bbox = featuresource.getBounds(); + // 读取投影 + CoordinateReferenceSystem crs = featuresource.getSchema().getCoordinateReferenceSystem(); + // 特征总数 + int count = featuresource.getCount(Query.ALL); + // 获取当前数据的geometry类型（点、线、面） + GeometryType geometryType = featuresource.getSchema().getGeometryDescriptor().getType(); + // 读取要素 + simpleFeatureCollection = (SimpleFeatureCollection) featuresource.getFeatures(); + // 获取当前矢量数据有哪些属性字段值 + List attributes = simpleFeatureCollection.getSchema().getAttributeDescriptors(); + + } catch (IOException e) { + e.printStackTrace(); + } + System.out.println("读取完成！"); + return simpleFeatureCollection; + } +}