使用geotools求出四至

Geotools是一个基于Java的开源GIS库，可以用于处理地理空间数据。要求出四至，需要先将要素转换为地理坐标系，然后获取要素的边界范围。

以下是一个简单的示例代码：

// 加载要素
File file = new File("path/to/shapefile.shp");
FileDataStore store = FileDataStoreFinder.getDataStore(file);
SimpleFeatureSource featureSource = store.getFeatureSource();

// 获取地理坐标系
CoordinateReferenceSystem crs = featureSource.getSchema().getCoordinateReferenceSystem();

// 创建转换器，将要素转换为地理坐标系
MathTransform transform = CRS.findMathTransform(DefaultGeographicCRS.WGS84, crs);
SimpleFeatureCollection featureCollection = featureSource.getFeatures();
SimpleFeatureIterator iterator = featureCollection.features();

// 计算要素的边界范围
ReferencedEnvelope envelope = new ReferencedEnvelope(crs);
while (iterator.hasNext()) {
    SimpleFeature feature = iterator.next();
    Geometry geometry = (Geometry) feature.getDefaultGeometry();
    Geometry transformedGeometry = JTS.transform(geometry, transform);
    envelope.expandToInclude(transformedGeometry.getEnvelopeInternal());
}

// 输出边界范围的四至
System.out.println("左下角经度：" + envelope.getMinX());
System.out.println("左下角纬度：" + envelope.getMinY());
System.out.println("右上角经度：" + envelope.getMaxX());
System.out.println("右上角纬度：" + envelope.getMaxY());

需要注意的是，在使用Geotools时，要素的坐标系需要与计算四至时使用的坐标系一致，否则计算结果可能会有误。