Java地理工具如何创建覆盖格网

<dependency>
    <groupId>org.geotools</groupId>
    <artifactId>gt-shapefile</artifactId>
    <version>xx.xx</version>
</dependency>

import org.geotools.coverage.grid.GridCoverage2D;
import org.geotools.coverage.grid.GridGeometry2D;
import org.geotools.geometry.Envelope2D;
import org.geotools.referencing.crs.DefaultGeographicCRS;
import org.opengis.geometry.Envelope;

// 定义格网属性
int numColumns = 10; // 列数
int numRows = 10; // 行数
double minX = -180.0; // 最小经度
double minY = -90.0; // 最小纬度
double maxX = 180.0; // 最大经度
double maxY = 90.0; // 最大纬度

// 创建格网生成器
Envelope envelope = new Envelope2D(DefaultGeographicCRS.WGS84, minX, minY, maxX - minX, maxY - minY);
GridGeometry2D gridGeometry = new GridGeometry2D(envelope, numColumns, numRows);
GridCoverage2D gridCoverage = new GridCoverage2D(gridGeometry);

// 遍历格网的每个单元格
for (int row = 0; row < numRows; row++) {
    for (int col = 0; col < numColumns; col++) {
        // 计算单元格的中心点坐标
        double cellCenterX = minX + (col + 0.5) * (maxX - minX) / numColumns;
        double cellCenterY = minY + (row + 0.5) * (maxY - minY) / numRows;
        
        // 在此处可以进行其他操作，例如计算属性值等
        // ...
    }
}