Wie kann ein USGS-Orthobild vom WorldImage-Format in gekacheltes GeoTiff für GeoServer konvertiert werden?

10

USGS bietet Orthoimage-Vierecke mit unterschiedlicher Auflösung im WorldImage-Format (* .tif mit einer zugehörigen * .tfw-Weltdatei). Einige dieser Vierecke können von GeoServer (2.1.0 - 2.1.2) verwendet werden, indem ein WorldImage-Speicher erstellt und die Ebene veröffentlicht wird (einige von ihnen lösen beim Versuch, die Ebene zu veröffentlichen, eine Ausnahme aus). Beim Veröffentlichen einer GeoTiff-Ebene, die von einer derzeit funktionierenden WorldImage-Ebene konvertiert wurde, tritt die folgende Ausnahme auf:

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)
...

Ich habe versucht, gdal_translate folgendermaßen zu verwenden (1234 durch echte Werte ersetzen): gdal_translate -of GTiff -a_srs "EPSG:2239" -co "TILED=YES" 1234.tif ../geotiff/1234.tifDie Konvertierung ist erfolgreich, und ich kann erfolgreich einen GeoTiff-Speicher in GeoServer erstellen, den Layer jedoch aufgrund des oben beschriebenen Fehlers nicht veröffentlichen. Dies lässt mich glauben, dass ich versuche, Ebenen (in beiden Fällen) zu veröffentlichen, die nicht richtig formatiert sind. Was ist also der richtige Weg, um ein [funktionierendes] USGS-Orthobild im WorldImage-Format in ein gekacheltes GeoTiff für GeoServer 2.1.2 zu konvertieren?

Matt
quelle
Ich denke, wir müssen wissen, welchen Epsg-Code Sie tatsächlich verwenden, da dies wahrscheinlich das wichtige Bit ist
Ian Turton
Bearbeitet, um ein EPSG anzuzeigen, das ich mit einem 2007-Datensatz aus EarthExplorer ausprobiert habe.
Matt
Wissen Sie, was der SRS des Eingangstiffs ist?
Ian Turton
1
Es sieht so aus, als ob Sie nicht zur Hauptursache gelangen müssen, da Parameter message can't be nullein Apache-Fehler ausgelöst wird, wenn eine FeedbackMessage keine Nachrichtenzeichenfolge enthält (siehe hier: bit.ly/zBUjy3 für die wichtigsten Details). Ich schlage vor, dass Sie den gesamten Stack-Trace veröffentlichen, damit wir möglicherweise feststellen können, welche Nachricht GeoServer uns mitteilen möchte.
MerseyViking
3
gdal_translate -a_srswird nicht neu projiziert, sondern weist lediglich die Projektion des Bildes zu (oder überschreibt sie). Wenn das Bild keine Begrenzungsinformationen enthält, wird 1 US-Fuß (für EPSG: 2239) pro Pixel angenommen und das Bild am Ursprung des Koordinatensystems platziert. Befindet sich der Begrenzungsrahmen in einer anderen Einheit als US-Fuß, kann Ihr Bild je nach Einheit und Ursprung der Projektion, in der es digitalisiert wurde, überall oder in beliebiger Größe sein.
MerseyViking

Antworten:

1

Warum nicht einfach in Geotiff herunterladen? USGS hatte früher diese Option. Die USGS-Liste der Site für dieses Format lautet The National Map und US Topo.

Lewis
quelle