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?
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.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.Réponses:
Si vous obtenez des géopdf de USGS et que vous souhaitez convertir en géotiff - et que vous avez une boîte Linux avec une capacité perl ou mod_perl et également des bibliothèques GDAL / OGR - vous pouvez les exécuter par lots avec le script geopdf2gtiff.pl comme indiqué ici: https: //github.com/godfreja/qtxastir/blob/master/scripts/geopdf2gtiff.pl
la source
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.
la source