¿Cuál es la forma correcta de convertir una Ortoimagen de USGS del formato WorldImage a GeoTiff en mosaico para GeoServer?

10

USGS proporciona cuadrángulos de ortoimagen de diferente resolución en el formato WorldImage (* .tif con un archivo mundial * .tfw adjunto). GeoServer puede utilizar algunos de estos cuadrángulos (2.1.0 - 2.1.2) creando una tienda WorldImage y publicando la capa (algunos de ellos arrojan una excepción al intentar publicar la capa). Al publicar una capa GeoTiff que se ha convertido de una capa WorldImage que funciona actualmente, obtengo la siguiente excepción:

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

Intenté usar gdal_translate de la siguiente manera (sustituya 1234 con valores reales): gdal_translate -of GTiff -a_srs "EPSG:2239" -co "TILED=YES" 1234.tif ../geotiff/1234.tifla conversión se realiza correctamente y puedo crear con éxito una tienda GeoTiff en GeoServer, pero no puedo publicar la capa debido al error descrito anteriormente. Esto me lleva a creer que estoy tratando de publicar capas (en ambos casos) que no están formateadas correctamente. Por lo tanto, ¿cuál es la forma correcta de convertir una ortoimagen USGS [en funcionamiento] en formato WorldImage en un GeoTiff en mosaico para GeoServer 2.1.2?

Mate
fuente
Creo que necesitamos saber qué código epsg está usando realmente, ya que probablemente sea el bit importante
Ian Turton
Editado para mostrar un EPSG que he probado con un conjunto de datos de 2007 recuperado de EarthExplorer.
Matt
¿Sabes cuál es el SRS del tiff de entrada?
Ian Turton
1
Parece que no ha llegado a la causa raíz porque se Parameter message can't be nullproduce un error de Apache cuando un FeedbackMessage no tiene una cadena de mensaje (consulte aquí: bit.ly/zBUjy3 para obtener los detalles sangrientos). Le sugiero que publique todo el seguimiento de la pila para que podamos identificar qué mensaje GeoServer está tratando de decirnos.
MerseyViking
3
gdal_translate -a_srsno reproyecta, simplemente asigna (o sobrescribe) la proyección de la imagen. Si la imagen no tiene información delimitadora, asumirá 1 pie estadounidense (para EPSG: 2239) por píxel y colocará la imagen en el origen del sistema de coordenadas. Si el cuadro de límite se encuentra en una unidad distinta de los pies de los Estados Unidos, entonces su imagen podría estar en cualquier lugar o en cualquier tamaño dependiendo de las unidades y el origen de cualquier proyección que se digitalizó en.
MerseyViking

Respuestas:

1

¿Por qué no simplemente descargar en geotiff? USGS solía tener esa opción. La lista de USGS del sitio para este formato son The National Map y US Topo.

Luis
fuente