查找Shapefile中是否存在纬度/经度的Java解决方案

Question

我的任务是尝试查找给定的纬度/经度(以度为单位，例如 -50.55555、78.232222)是否存在于给定的 Shapefile 中。

我从这里下载的 Shapefile:https://www.census.gov/geo/maps-data/data/cbf/cbf_nation.html (它包含美国及其所有领土的形状文件)。

我正在寻找一种基于 Java 的解决方案，在给定随机生成的纬度/经度的情况下，确定该点是否位于美国境内。我发现有用(但我是新手)的一种工具是 GeoTools。

那么几个问题:

• 是否可以使用 GeoTools 来完成此要求(确定某些纬度/经度是否位于美国)？

• 如果 GeoTools(或单独的 GeoTools)无法做到这一点，是否有另一种方法可以为此找到 Java 解决方案(无需求助于 Google map 等网络服务)。

另请注意:我下载了 QGIS 并使用 mmqgis 插件能够导出我上面提到的 Shapefile 的所有形状文件内容(shape_id、lat、long)。也许我可以使用另一个 Java 库来读取这些形状并测试随机纬度/经度是否落在这些形状内？

我是这一切的新手，任何/所有问题请提问!非常感谢。

Answer 1

你可以通过多种方式做到这一点，但我会这样做:

   private SimpleFeatureCollection features;
   private ReferencedEnvelope env;
   private boolean isInShape(Point p) {
        if (!env.contains(p.getCoordinate())) {
            return false;
        }
        Expression propertyName = filterFactory.property(features.getSchema()
                .getGeometryDescriptor().getName());
        Filter filter = filterFactory.contains(propertyName,
                filterFactory.literal(p));
        SimpleFeatureCollection sub = features.subCollection(filter);
        if (sub.size() > 0) {
            return true;
        }
        return false;
    }

然后读入您的 shapefile 并提取特征并向其 throw 点:

public static void main(String[] args) throws IOException {
    File file = null;
    if (args.length == 0) {
        // display a data store file chooser dialog for shapefiles
        file = JFileDataStoreChooser.showOpenFile("shp", null);
        if (file == null) {
            return;
        }
    } else {
        file = new File(args[0]);
        if (!file.exists()) {
            System.err.println(file + " doesn't exist");
            return;
        }
    }
    PointInPolygon tester = new PointInPolygon();
    FileDataStore store = FileDataStoreFinder.getDataStore(file);
    SimpleFeatureSource featureSource = store.getFeatureSource();
    tester.setFeatures(featureSource.getFeatures());
    GeometryFactory fac = new GeometryFactory();
    for (int i = 0; i < 1000; i++) {

        double lat = (Math.random() * 180.0) - 90.0;
        double lon = (Math.random() * 360.0) - 180.0;
        Point p = fac.createPoint(new Coordinate(lat, lon));
        boolean flag = tester.isInShape(p);
        if (flag) {
            System.out.println(p + " is in States ");
        }
    }
}

完整的例子是here .