Как правильно преобразовать USGS Orthoimage из формата WorldImage в мозаичный GeoTiff для GeoServer?

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?

— Matt
источник

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

— Ian Turton

Отредактирован, чтобы показать EPSG, который я пробовал с набором данных 2007 года, полученным из EarthExplorer.

— Мэтт

Вы знаете, что такое SRS входного TIFF?

— Ян Тертон

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

— MerseyViking

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

— MerseyViking

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

— cm1
источник

Почему бы просто не скачать в Geotiff? USGS имел такую возможность. Список USGS на сайте для этого формата: Национальная карта и Топо США.

— Льюис
источник