Hi,
As part
of our agriculture-based project we need to store and process
Landsat 8 imagery for large areas in various countries in
Europe.
We are
interested in using GeoMesa as solution since it seems to
provide many of the features we need to store, visualize and
process (using the additional Geoserver plugins) the large
amount raster images we have to consider in the project .
The
natural earth example can be visualized without problem,
however when it comes to visualizing the ingested Landsat 8
image via a WMS request I am getting the following:
Error:
Error rendering coverage on the fast path
org.locationtech.geomesa.shade.google.common.util.concurrent.UncheckedExecutionException:
java.lang.IllegalArgumentException: duplicate key: 0.002526
java.lang.IllegalArgumentException:
duplicate key: 0.002526
duplicate
key: 0.002526
</ServiceException></ServiceExceptionReport>
Somethings
regarding the steps I followed for ingesting the Landsat
rasters:
- We
need to store the 11 (B1 to B11) bands that compose a
Landsat 8 scene, therefore we created one layer per band
(AFAIK multiband rasters are not supported yet by Geomesa)
- The
Landsat 8 scene is reprojected to EPSG:4326 (from
EPSG:32630): gdalwarp -t_srs EPSG:4326
LC82020322015275LGN00_B3.TIF
LC82020322015275LGN00_B3_4326.TIF
-
Image pyramid creation: gdal_retile.py
-v -pyramidOnly -r lanczos -levels 4 -ps 675 675 -targetDir
./ LC82020322015275LGN00_B3_4326.TIF
- Rasters ingestion: geomesa ingestraster -u
geomesa -p geomesa -t LANDSAT8_B3 -f "./1/" -fmt geotiff
Same procedure for pyramid
levels 2, 3 and 4
Do you have any idea about
what may be causing the duplicate key problem?
I am attaching also the
output from gdalinfo command applied to the LC82020322015275LGN00_B3.TIF image in case the
info might be of help.
Best regards