วิธีที่ถูกต้องในการแปลง USGS Orthoimage จากรูปแบบ WorldImage เป็น GeoTiff แบบเรียงต่อกันสำหรับ GeoServer คืออะไร?


10

USGS ให้ความละเอียด orthoimage quadrangles ที่แตกต่างกันในรูปแบบ WorldImage (* .tif พร้อมกับไฟล์ * .tfw world world ประกอบ) quadrangles เหล่านี้บางส่วนสามารถใช้ได้โดย 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 ได้สำเร็จ แต่ไม่สามารถเผยแพร่เลเยอร์ได้เนื่องจากข้อผิดพลาดที่อธิบายไว้ข้างต้น สิ่งนี้ทำให้ฉันเชื่อว่าฉันกำลังพยายามเผยแพร่เลเยอร์ (ทั้งสองกรณี) ที่จัดรูปแบบไม่ถูกต้อง ดังนั้นวิธีที่ถูกต้องในการแปลง [ทำงาน] USO orthoimage ในรูปแบบ WorldImage เป็น GeoTiff แบบเรียงต่อกันสำหรับ GeoServer 2.1.2 คืออะไร?


ฉันคิดว่าเราจำเป็นต้องรู้รหัส epsg ที่คุณใช้จริง ๆ ซึ่งอาจเป็นบิตสำคัญ
Ian Turton

แก้ไขเพื่อแสดง EPSG ที่ฉันได้ลองกับชุดข้อมูล 2007 ที่ดึงมาจาก EarthExplorer
Matt

คุณรู้หรือไม่ว่า SRS ของอินพุต tiff คืออะไร
Ian Turton

1
ดูเหมือนว่าคุณไม่ได้เป็นต้นเหตุเพราะParameter message can't be nullเป็นข้อผิดพลาด Apache ที่เกิดขึ้นเมื่อ FeedbackMessage ไม่มีสตริงข้อความ (ดูที่นี่: bit.ly/zBUjy3สำหรับรายละเอียดเต็มไปด้วยเลือด) ฉันแนะนำให้คุณโพสต์การติดตามสแต็กทั้งหมดเพื่อให้เราสามารถระบุได้ว่าข้อความใดที่ GeoServer พยายามบอกเรา
MerseyViking

3
gdal_translate -a_srsไม่ได้ปฏิเสธมันเพียง แต่กำหนด (หรือเขียนทับ) การฉายภาพ หากภาพไม่มีข้อมูลที่มีขอบเขตมันจะถือว่า 1 US foot (สำหรับ 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

ทำไมไม่เพียงดาวน์โหลดใน geotiff? USGS เคยมีตัวเลือกนั้น เว็บไซต์ USGS สำหรับรูปแบบนี้คือ The National Map และ US Topo

โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.