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.tif
Konwersja 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?
Parameter message can't be null
błą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ć.gdal_translate -a_srs
nie 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.Odpowiedzi:
Jeśli otrzymujesz geopdf z USGS i chcesz przekonwertować na geotiff - i masz Linux-a z funkcją perl lub mod_perl, a także biblioteki GDAL / OGR - możesz wsadowo uruchomić je za pomocą skryptu geopdf2gtiff.pl, jak pokazano tutaj: https: //github.com/godfreja/qtxastir/blob/master/scripts/geopdf2gtiff.pl
źródło
Dlaczego nie pobrać w geotiff? USGS miała tę opcję. Witryna USGS dla tego formatu to The National Map i US Topo.
źródło