Qual é a maneira correta de converter uma USGS Orthoimage do formato WorldImage em GeoTiff para GeoServer lado a lado?

10

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?

Matt
fonte
Acho que precisamos de saber qual o código EPSG você está realmente usando como é provavelmente a pouco importante
Ian Turton
Editado para mostrar um EPSG que tentei com um conjunto de dados de 2007 recuperado do EarthExplorer.
Matt
Você sabe qual é o SRS do tiff de entrada?
Ian Turton
1
Parece que você não chegou à causa raiz porque ocorreuParameter message can't be null um 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.
MerseyViking
3
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.
MerseyViking

Respostas:

1

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.

Lewis
fonte