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.tif
la 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?
Parameter message can't be null
produce 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.gdal_translate -a_srs
no 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.Respuestas:
Si obtiene geopdf de USGS y desea convertir a geotiff, y tiene un cuadro de Linux con capacidad perl o mod_perl y también bibliotecas GDAL / OGR, puede ejecutarlos por lotes con el script geopdf2gtiff.pl como se muestra aquí: https: //github.com/godfreja/qtxastir/blob/master/scripts/geopdf2gtiff.pl
fuente
¿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.
fuente