Jaki jest prawidłowy sposób konwersji USGS Orthoimage z formatu WorldImage na kafelkowy GeoTiff dla GeoServer?


10

USGS zapewnia różne rozdzielczości czworokątów ortoobrazowych w formacie WorldImage (* .tif z dołączonym plikiem światowym * .tfw). Niektóre z tych czworokątów mogą być używane przez GeoServer (2.1.0 - 2.1.2) poprzez utworzenie magazynu WorldImage i opublikowanie warstwy (niektóre z nich rzucają wyjątek podczas próby opublikowania warstwy). Podczas publikowania warstwy GeoTiff przekonwertowanej z obecnie działającej warstwy WorldImage otrzymuję następujący wyjątek:

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

Próbowałem użyć gdal_translate w następujący sposób ( zamień 1234 na wartości rzeczywiste): gdal_translate -of GTiff -a_srs "EPSG:2239" -co "TILED=YES" 1234.tif ../geotiff/1234.tifKonwersja się powiodła i mogę z powodzeniem utworzyć magazyn GeoTiff w GeoServer, ale nie mogę opublikować warstwy z powodu błędu opisanego powyżej. To prowadzi mnie do przekonania, że ​​próbuję opublikować warstwy (w obu przypadkach), które nie są odpowiednio sformatowane. Zatem, jaki jest prawidłowy sposób konwersji [działającego] ortoobrazu USGS w formacie WorldImage na kafelkowy GeoTiff dla GeoServer 2.1.2?


Myślę, że musimy wiedzieć, którego kodu epsg tak naprawdę używasz, ponieważ jest to prawdopodobnie najważniejszy element
Ian Turton

Edytowane, aby pokazać plik EPSG, którego próbowałem z zestawem danych 2007 pobranym z EarthExplorer.
Matt

Czy wiesz, czym jest SRS wejściowy Tiff?
Ian Turton

1
Wygląda na to, że nie dotarłeś do głównej przyczyny, ponieważ Parameter message can't be nullbłąd Apache jest generowany, gdy komunikat FeedbackMessage nie ma ciągu komunikatu (szczegóły: krw.ly/zBUjy3 ). Sugeruję opublikowanie śledzenia całego stosu, abyśmy mogli określić, jaki komunikat GeoServer próbuje nam przekazać.
MerseyViking,

3
gdal_translate -a_srsnie wykonuje ponownej projekcji, jedynie przypisuje (lub zastępuje) projekcję obrazu. Jeśli obraz nie zawiera informacji o obrzeżach, przyjmuje 1 stopę amerykańską (dla EPSG: 2239) na piksel i umieszcza obraz na początku układu współrzędnych. Jeśli obwiednia znajduje się w jednostce innej niż stopy amerykańskie, wówczas obraz może być w dowolnym miejscu lub w dowolnym rozmiarze, w zależności od jednostek i pochodzenia dowolnej projekcji, w której został on przetworzony.
MerseyViking

Odpowiedzi:



1

Dlaczego nie pobrać w geotiff? USGS miała tę opcję. Witryna USGS dla tego formatu to The National Map i US Topo.

Korzystając z naszej strony potwierdzasz, że przeczytałeś(-aś) i rozumiesz nasze zasady używania plików cookie i zasady ochrony prywatności.
Licensed under cc by-sa 3.0 with attribution required.