Apa cara yang benar untuk mengonversi Orthoimage USGS dari format WorldImage ke ubin GeoTiff untuk GeoServer?

10

USGS menyediakan resolusi ortoimage quadrangles berbeda dalam format WorldImage (* .tif dengan file dunia * .tfw yang menyertai). Beberapa quadrangles ini dapat digunakan oleh GeoServer (2.1.0 - 2.1.2) dengan membuat toko WorldImage dan menerbitkan layer (beberapa dari mereka mengeluarkan pengecualian ketika mencoba untuk mempublikasikan layer). Saat menerbitkan lapisan GeoTiff yang telah dikonversi dari lapisan WorldImage yang saat ini berfungsi, saya mendapatkan pengecualian berikut:

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

Saya mencoba menggunakan gdal_translate dengan cara berikut (gantikan 1234 dengan nilai real): gdal_translate -of GTiff -a_srs "EPSG:2239" -co "TILED=YES" 1234.tif ../geotiff/1234.tifKonversi berhasil, dan saya berhasil membuat toko GeoTiff di GeoServer, tetapi tidak bisa menerbitkan lapisan karena kesalahan yang dijelaskan di atas. Ini membuat saya percaya bahwa saya mencoba untuk mempublikasikan lapisan (dalam kedua kasus) yang tidak diformat dengan benar. Jadi, apa cara yang benar untuk mengonversi orthoimage USGS [yang berfungsi] dalam format WorldImage ke GeoTiff ubin untuk GeoServer 2.1.2?

Mat
sumber
Saya pikir kita perlu tahu kode epsg mana yang sebenarnya Anda gunakan karena ini mungkin yang paling penting
Ian Turton
Diedit untuk menunjukkan EPSG yang saya coba dengan kumpulan data 2007 yang diambil dari EarthExplorer.
Matt
Apakah Anda tahu apa SRS masukan tiff?
Ian Turton
1
Sepertinya Anda tidak sampai ke akar penyebab karena Parameter message can't be nullmerupakan kesalahan Apache yang dilemparkan ketika FeedbackMessage tidak memiliki string pesan (lihat di sini: bit.ly/zBUjy3 untuk detail berdarah). Saya sarankan Anda memposting seluruh jejak tumpukan sehingga kami mungkin dapat mengidentifikasi pesan apa yang coba disampaikan GeoServer kepada kami.
MerseyViking
3
gdal_translate -a_srstidak memproyeksi ulang, itu hanya memberikan (atau menimpa) proyeksi gambar. Jika gambar tidak memiliki informasi pembatas, ia akan mengasumsikan 1 kaki AS (untuk EPSG: 2239) per piksel dan menempatkan gambar pada asal sistem koordinat. Jika kotak pembatas berada dalam unit selain dari kaki AS, maka gambar Anda dapat berada di mana saja atau pada ukuran apa pun tergantung pada unit dan asal dari proyeksi apa pun yang didigitalkan.
MerseyViking

Jawaban:

1

Mengapa tidak mengunduh saja di geotiff? USGS dulu memiliki opsi itu. Situs USGS daftar untuk format ini adalah The National Map dan US Topo.

Lewis
sumber