计算栅格中的波段数

Question

我一直在尝试计算光栅文件中的波段数。虽然下面的代码大部分都有效，但当我在投影栅格上尝试代码时似乎会抛出异常。

我试过以下代码：

AbstractGridFormat format = GridFormatFinder.findFormat(rasterFile);
        // this is a bit hacky but does make more geotiffs work
        Hints hints = new Hints();
        if (format instanceof GeoTiffFormat) {
            hints = new Hints(Hints.FORCE_LONGITUDE_FIRST_AXIS_ORDER, Boolean.TRUE);
        }
        reader = format.getReader(rasterFile, hints);
        GridCoverage2D result = reader.read(null);
        int bands = result.getNumSampleDimensions();
        System.out.println(bands)

我收到以下异常：

org.geotools.data.DataSourceException
    at org.geotools.gce.geotiff.GeoTiffReader.<init>(GeoTiffReader.java:266)
    at org.geotools.gce.geotiff.GeoTiffFormat.getReader(GeoTiffFormat.java:307)
    at org.geotools.gce.geotiff.GeoTiffFormat.getReader(GeoTiffFormat.java:76)
    at ImageLab.displayLayers(ImageLab.java:74)
    at ImageLab.getLayersAndDisplay(ImageLab.java:36)
    at ImageLab.main(ImageLab.java:25)
Caused by: org.geotools.data.DataSourceException
    at org.geotools.gce.geotiff.GeoTiffReader.getHRInfo(GeoTiffReader.java:500)
    at org.geotools.gce.geotiff.GeoTiffReader.<init>(GeoTiffReader.java:253)
    ... 5 more
Caused by: java.lang.NullPointerException
    at org.geotools.referencing.cs.AbstractCS.axisUsingUnit(AbstractCS.java:503)
    at org.geotools.referencing.cs.DefaultCartesianCS.usingUnit(DefaultCartesianCS.java:291)
    at org.geotools.coverage.grid.io.imageio.geotiff.GeoTiffMetadata2CRSAdapter.createUserDefinedPCS(GeoTiffMetadata2CRSAdapter.java:750)
    at org.geotools.coverage.grid.io.imageio.geotiff.GeoTiffMetadata2CRSAdapter.createProjectedCoordinateReferenceSystem(GeoTiffMetadata2CRSAdapter.java:245)
    at org.geotools.coverage.grid.io.imageio.geotiff.GeoTiffMetadata2CRSAdapter.createCoordinateSystem(GeoTiffMetadata2CRSAdapter.java:181)
    at org.geotools.gce.geotiff.GeoTiffReader.getHRInfo(GeoTiffReader.java:330)
    ... 6 more

Exception in thread "main" java.lang.NullPointerException
    at ImageLab.displayLayers(ImageLab.java:75)
    at ImageLab.getLayersAndDisplay(ImageLab.java:36)
    at ImageLab.main(ImageLab.java:25)

对于投影数据，确实需要添加一些代码，因为 well.My 光栅文件位于 WyonmingNAD27 : 32056 投影系统中

Answer 1

我的代码中缺少以下依赖项。

<dependency>
<groupId>org.geotools</groupId>
<artifactId>gt-epsg-hsql</artifactId>
<version>${geotools.version}</version>
</dependency>

我添加了上面的依赖之后，CRS的问题就解决了

计算栅格中的波段数

Calculating number of bands in a raster

raster

geotools