Який правильний спосіб перетворити USGS Orthoimage з формату WorldImage в кахельний GeoTiff для GeoServer?


10

USGS надає ортографічні чотирикутники різної роздільної здатності у форматі WorldImage (* .tif із супровідним світовим файлом * .tfw). Деякі з цих чотирикутників може використовуватися GeoServer (2.1.0 - 2.1.2) шляхом створення магазину WorldImage та публікації шару (деякі з них кидають виняток при спробі опублікувати шар). При публікації шару GeoTiff, який був перетворений з поточно працюючого шару WorldImage, я отримую таке виняток:

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

Я спробував використовувати gdal_translate наступним чином (замінити 1234 реальними значеннями): gdal_translate -of GTiff -a_srs "EPSG:2239" -co "TILED=YES" 1234.tif ../geotiff/1234.tifПеретворення вдалося, і я можу успішно створити сховище GeoTiff в GeoServer, але не можу опублікувати шар через помилку, описану вище. Це призводить мене до думки, що я намагаюся публікувати шари (в обох випадках), які не форматовані належним чином. Таким чином, який правильний спосіб перетворити [робочий] орфографічний малюнок USGS у форматі WorldImage в кахельний GeoTiff для GeoServer 2.1.2?


Я думаю, нам потрібно знати, який код epsg ви використовуєте, оскільки це, мабуть, важливий біт
Ян Тертон

Відредагований для показу EPSG, який я пробував із набором даних 2007 року, отриманим із EarthExplorer.
Метт

Чи знаєте ви, що таке SRS вхідного Tiff?
Ян Тертон

1
Схоже, ви не потрапили до першопричини, оскільки Parameter message can't be nullпомилка Apache викинута, коли у FeedbackMessage немає рядка повідомлення (див. Тут: bit.ly/zBUjy3 для деталей Gory ). Я пропоную вам опублікувати весь слід стека, щоб ми могли виявити, яке повідомлення намагається нам сказати GeoServer.
MerseyViking

3
gdal_translate -a_srsне перепроектує, він просто призначає (або перезаписує) проекцію зображення. Якщо зображення не має обмежувальної інформації, воно передбачає 1 фут США (для EPSG: 2239) на піксель і розмістить зображення на початку координатної системи. Якщо обмежувальний ящик знаходиться в одиниці, окрім ніг США, то ваше зображення може бути в будь-якому місці або будь-якого розміру, залежно від одиниць та походження будь-якої проекції, в якій воно було оцифровано.
MerseyViking

Відповіді:


2

Якщо ви отримуєте geopdf з USGS і хочете конвертувати в geotiff - і у вас є вікно Linux з можливістю perl або mod_perl, а також бібліотеки GDAL / OGR - ви можете пакетно запустити їх за допомогою сценарію geopdf2gtiff.pl, як показано тут: https: //github.com/godfreja/qtxastir/blob/master/scripts/geopdf2gtiff.pl


1

Чому б просто не завантажити в геотиф? У цьому варіанті були USGS. Список USGS на сайті для цього формату - Національна карта та US Topo.

Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.