Qual è il modo corretto per convertire un USGS Orthoimage dal formato WorldImage in GeoTiff piastrellato per GeoServer?


10

USGS fornisce quadrangoli di ortoimmagini con risoluzione diversa nel formato WorldImage (* .tif con un file mondiale * .tfw di accompagnamento). Alcuni di questi quadrangoli possono essere utilizzati da GeoServer (2.1.0 - 2.1.2) creando un archivio WorldImage e pubblicando il layer (alcuni generano un'eccezione quando si tenta di pubblicare il layer). Quando pubblico un livello GeoTiff che è stato convertito da un livello WorldImage attualmente funzionante, ottengo la seguente eccezione:

2012-01-10 08:38:33,228 WARN [referencing.factory] - Axis elements found in a wkt definition, the force longitude first axis order hint might not be respected:
PROJCS["WGS84 / Simple Mercator", GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS_1984", 6378137.0, 298.257223563]], PRIMEM["Greenwich", 0.0], UNIT["degree", 0.017453292519943295]], PROJECTION["Mercator_1SP_Google"], PARAMETER["latitude_of_origin", 0.0], PARAMETER["central_meridian", 0.0], PARAMETER["scale_factor", 1.0], PARAMETER["false_easting", 0.0], PARAMETER["false_northing", 0.0], UNIT["m", 1.0], AXIS["x", EAST], AXIS["y", NORTH], AUTHORITY["EPSG","54004"]]
2012-01-10 08:38:39,674 INFO [geoserver.web] - Error saving layer
java.lang.NullPointerException
    at org.geoserver.catalog.impl.CoverageStoreInfoImpl.getFormat(CoverageStoreInfoImpl.java:42)
...

Ho provato a utilizzare gdal_translate nel modo seguente (sostituire 1234 con valori reali): gdal_translate -of GTiff -a_srs "EPSG:2239" -co "TILED=YES" 1234.tif ../geotiff/1234.tifla conversione ha esito positivo e posso creare correttamente un archivio GeoTiff in GeoServer, ma non riesco a pubblicare il livello a causa dell'errore sopra descritto. Questo mi porta a credere che sto provando a pubblicare livelli (in entrambi i casi) che non sono formattati correttamente. Quindi, qual è il modo corretto di convertire un ortogaggio USGS [funzionante] in formato WorldImage in un GeoTiff piastrellato per GeoServer 2.1.2?


Penso che dobbiamo sapere quale codice epsg stai effettivamente usando in quanto è probabilmente la parte importante
Ian Turton

Modificato per mostrare un EPSG che ho provato con un set di dati del 2007 recuperato da EarthExplorer.
Matt,

Sai qual è l'ingresso SRS di TIFF?
Ian Turton

1
Sembra che tu non sia arrivato alla causa principale perché Parameter message can't be nullè stato generato un errore Apache quando un FeedbackMessage non ha una stringa di messaggio (vedi qui: bit.ly/zBUjy3 per i dettagli cruenti). Ti suggerisco di pubblicare l'intera traccia dello stack in modo che possiamo essere in grado di identificare quale messaggio GeoServer sta cercando di dirci.
MerseyViking,

3
gdal_translate -a_srsnon riproiettare, semplicemente assegna (o sovrascrive) la proiezione dell'immagine. Se l'immagine non ha informazioni di delimitazione, assumerà 1 piede USA (per EPSG: 2239) per pixel e posizionerà l'immagine all'origine del sistema di coordinate. Se il rettangolo di selezione si trova in un'unità diversa dai piedi statunitensi, l'immagine potrebbe essere ovunque o di qualsiasi dimensione a seconda delle unità e dell'origine di qualsiasi proiezione in cui è stata digitalizzata.
MerseyViking

Risposte:



1

Perché non scaricare semplicemente in geotiff? USGS aveva quell'opzione. L'elenco USGS del sito per questo formato è The National Map e US Topo.

Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.