基于 Java 和 GeoTools 根据矢量 BBox 自动生成格网文件

2、生成经纬网要素集合

在获取了矢量数据的 BBox 后，我们就可以根据 BBox 的最大最小值来生成对应的经纬线。因此这里就成了文章的核心，如何根据 bbox 来生成经纬线数据。由于最后我们需要生成一个完整标准的 shp 文件，因此我们首先就需要创建一个要素信息，创建要素类型的代码如下：

/** 创建要素类型 */
private static SimpleFeatureType createFeatureType() throws SchemaException {
    SimpleFeatureTypeBuilder builder = new SimpleFeatureTypeBuilder();
    builder.setName("Graticule");
    builder.setCRS(WGS84);
    builder.add("geometry", LineString.class); // 几何字段
    builder.add("type", String.class); // 类型（经线/纬线）
    builder.add("degree", Double.class); // 度数
    return builder.buildFeatureType();
}

上面的代码就好像是在关系型数据库中，我们已经创建了一个表结构。里面包含了三个字段，第一个字段是空间属性数据，使用 geometry 来进行存储，除此之外还有两个属性字段，即类型（纬线还是经线），最后是一个标准的度数数据，这里采用的是默认的 WGS84 的地理坐标。有了表结构之后，接下来就是往表中写入数据。根据 bbox 来自动生成格网的关键代码如下：

/** 创建线要素 */
private static SimpleFeature createFeature(SimpleFeatureType type, LineString line, String typeStr, double degree) {
    SimpleFeatureBuilder featureBuilder = new SimpleFeatureBuilder(type);
    featureBuilder.set("geometry", line);
    featureBuilder.set("type", typeStr);
    featureBuilder.set("degree", degree);
    return featureBuilder.buildFeature(null);
}

/** 创建线段 */
private static LineString createLine(double x1, double y1, double x2, double y2) {
    org.locationtech.jts.geom.Coordinate[] coordinates = {
        new org.locationtech.jts.geom.Coordinate(x1, y1),
        new org.locationtech.jts.geom.Coordinate(x2, y2)
    };
    return JTSFactoryFinder.getGeometryFactory().createLineString(coordinates);
}

接下来是生成格网的核心算法，这里采用循环的方式，自动按照矢量的边界范围进行生成。示例代码如下：

/**
 * 生成经纬网要素集合
 * @param minLon 最小经度
 * @param maxLon 最大经度
 * @param minLat 最小纬度
 * @param maxLat 最大纬度
 * @param interval 间隔（度）
 */
private static SimpleFeatureCollection generateGraticule(double minLon, double maxLon, double minLat, double maxLat, double interval) throws SchemaException {
    // 创建要素类型
    SimpleFeatureType TYPE = createFeatureType();
    // 存储所有要素的集合
    List<SimpleFeature> features = new ArrayList<>();
    // 生成经线（固定经度，纬度从 minLat 到 maxLat）
    for (double lon = minLon; lon <= maxLon; lon += interval) {
        LineString line = createLine(lon, minLat, lon, maxLat);
        features.add(createFeature(TYPE, line, "LON", lon));
    }
    // 生成纬线（固定纬度，经度从 minLon 到 maxLon）
    for (double lat = minLat; lat <= maxLat; lat += interval) {
        LineString line = createLine(minLon, lat, maxLon, lat);
        features.add(createFeature(TYPE, line, "LAT", lat));
    }
    return new ListFeatureCollection(TYPE, features);
}

3、写入 Shapefile

这里我们将 shp 文件准备完成之后，接下来就是最重要的一步，即往 shp 文件中写入要素信息，这样就实现了生成的格网 shapefile 的落盘处理。写出 shapefile 的实现比较简单，示例代码如下：

private static void writeToShapefile(SimpleFeatureCollection features, File file) throws IOException {
    // 方法 1：直接使用 ShapefileDataStoreFactory
    ShapefileDataStore dataStore = new ShapefileDataStore(file.toURI().toURL());
    dataStore.createSchema(features.getSchema());
    // 写入数据
    try (FeatureWriter<SimpleFeatureType, SimpleFeature> writer = dataStore.getFeatureWriterAppend(null)) {
        SimpleFeatureIterator iterator = features.features();
        while (iterator.hasNext()) {
            SimpleFeature feature = iterator.next();
            SimpleFeature newFeature = writer.next();
            newFeature.setAttributes(feature.getAttributes());
            writer.write();
        }
        iterator.close();
    }
    dataStore.dispose();
}

在需要调用的逻辑如下，需要定义一个本地保存的矢量文件。

// 3、写入 Shapefile
// 定义输出路径
File shapefile = new File("D:/shp_test/graticule250508.shp");
writeToShapefile(features, shapefile);
System.out.println("完成!!!");

三、结果输出与验证

本节将对经过前面的步骤生成的格网成果进行输出和验证。

1、格网文件输出格式

经过前面的步骤在目标磁盘可以看到已经成功的生成我们需要的 shp 文件，打开磁盘目录如下图所示： 表名需要的文件已经成功的生成到我们的目标磁盘中。

2、GIS 工具验证

将生成 shp 文件使用 QGIS 这样的 GIS 软件中进行查看，可以看到以下的效果： 为了验证当前生成的经纬网线数据是否是符合目标的文件 bbox 范围的，我们再叠加原始的数据图层来校验一下。 可以看到生成的格网数据是完全按照我们的预期来生成，需求得到了满足。

四、总结

以上就是本文的主要内容，本实践旨在探索利用 Java 和 GeoTools 的强大功能，构建一个高效、准确的矢量文件转格网文件的系统。从矢量文件的读取与解析开始，深入研究如何根据矢量数据的空间分布和属性特征，自动划分格网并生成相应的格网文件。在实践过程中，我们将面临诸多挑战，如矢量数据的复杂性、格网划分的精度控制、数据转换的效率提升等，需要通过不断尝试和优化算法来解决。本文将详细阐述基于 Java 和 GeoTools 进行矢量文件转格网文件实践的过程、方法和技术细节，包括环境搭建、矢量数据处理、格网生成算法设计以及结果验证等方面。希望本实践研究能够为地理空间数据处理领域提供一种新的方法和思路，提高矢量文件转格网文件的效率和质量，为相关领域的研究和应用提供有力的支持和参考。当然在生成的结果中可以看到，还是存在一些不足，即我们生成的格网的末端不是很平滑，因此在未来我们可以优化以下末端的处理。同时对于经纬线的标注可以将值进行格式化一下，这样看起来会更加符合人类的习惯。