geotools读取tiff文件中的高程范围

你可以使用 GeoTools 库来读取 TIFF 文件中的高程数据，并获取其范围。以下是一个示例代码：

import java.io.File;
import org.geotools.coverage.grid.GridCoverage2D;
import org.geotools.coverage.grid.GridEnvelope2D;
import org.geotools.coverage.grid.io.AbstractGridFormat;
import org.geotools.coverage.grid.io.GridFormatFinder;
import org.geotools.data.DataStore;
import org.geotools.data.DataStoreFinder;
import org.geotools.data.simple.SimpleFeatureCollection;
import org.geotools.data.simple.SimpleFeatureIterator;
import org.geotools.factory.CommonFactoryFinder;
import org.geotools.feature.FeatureCollection;
import org.geotools.feature.simple.SimpleFeatureBuilder;
import org.geotools.filter.text.cql2.CQL;
import org.geotools.geometry.jts.JTSFactoryFinder;
import org.geotools.referencing.CRS;
import org.geotools.referencing.GeodeticCalculator;
import org.geotools.referencing.crs.DefaultGeographicCRS;
import org.geotools.referencing.crs.DefaultVerticalCRS;
import org.geotools.referencing.operation.DefaultCoordinateOperationFactory;
import org.opengis.coverage.grid.GridCoverage;
import org.opengis.coverage.grid.GridEnvelope;
import org.opengis.coverage.grid.GridGeometry;
import org.opengis.coverage.grid.GridRange;
import org.opengis.feature.simple.SimpleFeature;
import org.opengis.filter.Filter;
import org.opengis.filter.FilterFactory2;
import org.opengis.filter.expression.Expression;
import org.opengis.geometry.DirectPosition;
import org.opengis.geometry.Envelope;
import org.opengis.geometry.MismatchedDimensionException;
import org.opengis.geometry.coordinate.Position;
import org.opengis.referencing.NoSuchAuthorityCodeException;
import org.opengis.referencing.crs.CoordinateReferenceSystem;

public class TiffReader {

    public static void main(String[] args) throws Exception {
        // 加载 TIFF 文件
        File file = new File("path/to/tiff/file.tif");
        AbstractGridFormat format = GridFormatFinder.findFormat(file);
        DataStore dataStore = format.getDataStore(file.toURI().toURL());
        String[] names = dataStore.getTypeNames();
        String typeName = names[0];

        // 打开图层，获取高程范围
        SimpleFeatureCollection collection = dataStore.getFeatureSource(typeName).getFeatures();
        SimpleFeatureIterator iterator = collection.features();
        Envelope envelope = null;
        while (iterator.hasNext()) {
            SimpleFeature feature = iterator.next();
            if (envelope == null) {
                envelope = feature.getBounds();
            } else {
                envelope.expandToInclude(feature.getBounds());
            }
        }
        iterator.close();

        // 获取高程范围的坐标系
        CoordinateReferenceSystem crs = envelope.getCoordinateReferenceSystem();

        // 打印高程范围
        System.out.println("Min X: " + envelope.getMinX());
        System.out.println("Min Y: " + envelope.getMinY());
        System.out.println("Max X: " + envelope.getMaxX());
        System.out.println("Max Y: " + envelope.getMaxY());

        // 计算高程范围的高度和宽度
        GeodeticCalculator calculator = new GeodeticCalculator();
        calculator.setStartingGeographicPoint(envelope.getMinX(), envelope.getMinY());
        calculator.setDestinationGeographicPoint(envelope.getMaxX(), envelope.getMinY());
        double width = calculator.getOrthodromicDistance();

        calculator.setStartingGeographicPoint(envelope.getMinX(), envelope.getMinY());
        calculator.setDestinationGeographicPoint(envelope.getMinX(), envelope.getMaxY());
        double height = calculator.getOrthodromicDistance();

        System.out.println("Width: " + width);
        System.out.println("Height: " + height);
    }
}

在上面的代码中，我们首先加载 TIFF 文件，然后打开图层，读取每个要素的边界框，并使用 `Envelope` 对象计算高程范围。然后，我们可以获取高程范围的坐标系，计算高程范围的宽度和高度，并将它们打印出来。