O USGS fornece diferentes quadrângulos de ortoimagem de resolução no formato WorldImage (* .tif com um arquivo mundial * .tfw). Alguns desses quadrângulos podem ser usados pelo GeoServer (2.1.0 - 2.1.2) criando uma loja WorldImage e publicando a camada (alguns deles lançam uma exceção ao tentar publicar a camada). Ao publicar uma camada GeoTiff que foi convertida de uma camada WorldImage atualmente em funcionamento, recebi a seguinte exceção:
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)
...
Tentei usar gdal_translate da seguinte maneira (substitua 1234 por valores reais): gdal_translate -of GTiff -a_srs "EPSG:2239" -co "TILED=YES" 1234.tif ../geotiff/1234.tifA conversão foi bem-sucedida e posso criar com êxito um repositório GeoTiff no GeoServer, mas não posso publicar a camada devido ao erro descrito acima. Isso me leva a acreditar que estou tentando publicar camadas (nos dois casos) que não estão formatadas corretamente. Portanto, qual é a maneira correta de converter uma ortoimagem [em funcionamento] do USGS no formato WorldImage em um GeoTiff lado a lado para GeoServer 2.1.2?
Parameter message can't be nullum erro do Apache quando um FeedbackMessage não possui uma sequência de mensagens (veja aqui: bit.ly/zBUjy3 para obter detalhes sórdidos). Sugiro que você publique todo o rastreamento da pilha, para que possamos identificar qual mensagem o GeoServer está tentando nos dizer.gdal_translate -a_srsnão reprojeta, apenas atribui (ou sobrescreve) a projeção da imagem. Se a imagem não tiver informações de contorno, ela assumirá 1 pé americano (para EPSG: 2239) por pixel e colocará a imagem na origem do sistema de coordenadas. Se a caixa delimitadora está em uma unidade diferente pés dos EUA, em seguida, sua imagem pode ser qualquer lugar ou em qualquer tamanho, dependendo das unidades e origem de qualquer projeção foi digitalizado no.Respostas:
Se você está obtendo geopdf do USGS e deseja converter em geotiff - e possui uma caixa Linux com capacidade perl ou mod_perl e também bibliotecas GDAL / OGR - você pode executá-las em lote com o script geopdf2gtiff.pl, como mostrado aqui: https: //github.com/godfreja/qtxastir/blob/master/scripts/geopdf2gtiff.pl
fonte
Por que não fazer o download no geotiff? O USGS costumava ter essa opção. A lista USGS do site para este formato é The National Map e US Topo.
fonte