Quelle est la bonne façon de convertir une orthoimage USGS du format WorldImage en GeoTiff en mosaïque pour GeoServer?

10

L'USGS fournit des quadrangles d'ortho-images de différentes résolutions au format WorldImage (* .tif avec un fichier mondial * .tfw). Certains de ces quadrangles peuvent être utilisés par GeoServer (2.1.0 - 2.1.2) en créant un magasin WorldImage et en publiant la couche (certains lancent une exception lors de la tentative de publication de la couche). Lors de la publication d'une couche GeoTiff qui a été convertie à partir d'une couche WorldImage en cours de fonctionnement, j'obtiens l'exception suivante:

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

J'ai essayé d'utiliser gdal_translate de la manière suivante (remplacez 1234 par des valeurs réelles): gdal_translate -of GTiff -a_srs "EPSG:2239" -co "TILED=YES" 1234.tif ../geotiff/1234.tifLa conversion réussit et je peux créer avec succès un magasin GeoTiff dans GeoServer, mais je ne peux pas publier la couche en raison de l'erreur décrite ci-dessus. Cela m'amène à croire que j'essaie de publier des couches (dans les deux cas) qui ne sont pas correctement formatées. Ainsi, quelle est la bonne façon de convertir une orthoimage [de travail] USGS au format WorldImage en un GeoTiff en mosaïque pour GeoServer 2.1.2?

Mat
la source
Je pense que nous devons savoir quel code EPSG vous utilisez réellement car c'est probablement le bit important
Ian Turton
Modifié pour afficher un EPSG que j'ai essayé avec un ensemble de données 2007 extrait d'EarthExplorer.
Matt
Savez-vous quel est le SRS du tiff d'entrée?
Ian Turton
1
Il semble que vous ne soyez pas arrivé à la cause principale, car Parameter message can't be nullune erreur Apache est générée lorsqu'un FeedbackMessage n'a pas de chaîne de message (voir ici: bit.ly/zBUjy3 pour les détails sanglants). Je vous suggère de publier toute la trace de la pile afin que nous puissions être en mesure d'identifier le message que GeoServer essaie de nous dire.
MerseyViking
3
gdal_translate -a_srsne reprojete pas, il assigne (ou écrase) simplement la projection de l'image. Si l'image n'a aucune information de délimitation, elle supposera 1 pied US (pour EPSG: 2239) par pixel et placera l'image à l'origine du système de coordonnées. Si le cadre de délimitation est dans une unité autre que les pieds américains, votre image peut être n'importe où ou à n'importe quelle taille en fonction des unités et de l'origine de la projection dans laquelle elle a été numérisée.
MerseyViking

Réponses:

1

Pourquoi ne pas simplement télécharger dans geotiff? USGS avait cette option. La liste du site USGS pour ce format est The National Map et US Topo.

lewis
la source