การทำให้ QGIS ตีความพิกัดเป็น lat-long แทนที่จะเป็น lat-long?


13

ฉันกำลังทำงานกับไฟล์ TRMM 3B31 เกี่ยวกับการตกตะกอนของ NASA เมื่อฉันโหลดมันลงบน QGIS เลือกการฉาย WGS84 พร้อมกับไฟล์รูปร่างอื่น ๆ แผนที่จะแสดงในแนวตั้งเช่น QGIS ตีความเป็นละติจูดว่าอะไรคือเส้นแวงและในทางกลับกัน ฉันได้ลองใช้ตัวเลือก gdal_translate -a_ullr แล้วเลือกพิกัดที่ดี แต่มันไม่ได้แก้ปัญหา: ภาพถูกบีบอัดและยืดออกเพื่อตอบสนองความต้องการของฉัน แต่นี่ไม่ใช่สิ่งที่ฉันต้องการ

เมื่อคุณโหลดไฟล์. csv คุณสามารถเลือกว่าจะต้องเป็นละติจูดและลองจิจูดอะไรมีวิธีการทำเช่นเดียวกันสำหรับ rasters หรือไม่

เนื่องจากฉันมีแรสเตอร์ 193 มันจะดีมากถ้าคุณสามารถแนะนำคำสั่งที่สามารถทำซ้ำได้

ป้อนคำอธิบายรูปภาพที่นี่

ภาพแสดงแผนที่ของประเทศอินเดียและแถบแรกของไฟล์ 3B31 เช่นการตกตะกอนของพื้นผิว เส้นสีดำสองเส้นคือเส้นขนาน 40 และ -40 ฉันต้องการให้พวกมันเป็นแนวนอนและด้วยเหตุนี้ข้อมูลฝนจึงคาบเกี่ยวกับอินเดีย


ภาพทดสอบจะมีประโยชน์มาก
user30184

@ user30184 gis.stackexchange.com/questions/89261/trmm-data-downloadให้การเชื่อมโยงมากมาย
AndreJ

คำตอบ:


7

ฉันลองไฟล์ 3b31 hdf จากmiradorและพบการวางแนวเดียวกันกับที่คุณพบ ดังนั้นฉันจึงสร้าง GCP เพื่อให้ทวนเข็มนาฬิกาหมุน:

gdal_translate -of VRT HDF4_SDS:UNKNOWN:"3b31.hdf":0 0.vrt
gdal_translate -of VRT -gcp 0 0 -180 -40 -gcp 160 0 -180 40 -gcp 0 720 180 -40 -gcp 160 720 180 40 0.vrt 1.vrt
gdalwarp -r bilinear -t_srs EPSG:4326 1.vrt 0.tif

ด้วยผลลัพธ์นี้:

ป้อนคำอธิบายรูปภาพที่นี่

การหมุนเวียนของฉันควรเป็นไปตามแท็กข้อมูลเมตา Origin = SOUTHWEST


นี่คือสิ่งที่ฉันกำลังมองหาถึงแม้ว่าฉันจะไม่สามารถใช้คำสั่งที่คุณเขียนในวิธีที่เหมาะสม คุณช่วยอธิบายการใช้งานในรายละเอียดเพิ่มเติมได้ไหม?
เปาโลซานตินี

วางไว้ในไฟล์ข้อความชื่อrotate.batในโฟลเดอร์เดียวกับไฟล์ HDF ตรวจสอบชื่อไฟล์ต้นฉบับเริ่มเชลล์ OSGEO4W นำทางไปยังโฟลเดอร์ที่มี cd ... และเรียกใช้แบตช์ การนำทาง DOS แบบเก่าที่ดีนั้นมีประโยชน์เสมอ ;-)
AndreJ

3

ฉันดาวน์โหลด 3B42RT.2016031806.7.7day.tif และ. tfw แต่ฉันไม่เห็นอะไรแปลก ๆ ในการอ้างอิงทางภูมิศาสตร์

ป้อนคำอธิบายรูปภาพที่นี่

รายงาน Gdalinfo ของภาพนี้คือ

gdalinfo 3B42RT.2016031806.7.7day.tif
Driver: GTiff/GeoTIFF
Files: 3B42RT.2016031806.7.7day.tif
       3B42RT.2016031806.7.7day.tfw
Size is 1440, 480
Coordinate System is `'
Origin = (-180.000000000000000,60.000000000000000)
Pixel Size = (0.250000000000000,-0.250000000000000)
Metadata:
  TIFFTAG_DATETIME=2016:03:18 10:25:21
  TIFFTAG_DOCUMENTNAME=./image/3B42RT.2016031806.7.7day.tif
  TIFFTAG_IMAGEDESCRIPTION=IDL TIFF file
  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
  TIFFTAG_SOFTWARE=IDL 7.1.1, ITT Visual Information Solutions
  TIFFTAG_XRESOLUTION=100
  TIFFTAG_YRESOLUTION=100
Image Structure Metadata:
  COMPRESSION=LZW
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (-180.0000000,  60.0000000)
Lower Left  (-180.0000000, -60.0000000)
Upper Right ( 180.0000000,  60.0000000)
Lower Right ( 180.0000000, -60.0000000)
Center      (   0.0000000,   0.0000000)
Band 1 Block=1440x2 Type=UInt16, ColorInterp=Gray

ดูเหมือนว่าไฟล์ 3b42 tif มีการอ้างอิงทางภูมิศาสตร์ที่ถูกต้องในขณะที่ไฟล์ 3b31 hdf ยังไม่มี
AndreJ

มันดูเหมือนว่าค่อนข้าง HDF มีความซับซ้อนและไม่ได้รับการสนับสนุนที่ดีเลิศที่ GDAL trac.osgeo.org/gdal/wiki/HDF
user30184

บางส่วนที่คุณพบใน Mirador นั้นมีการอ้างอิงทางภูมิศาสตร์การเลือกคุณสมบัติ NetCDF ในหน้าต่างการปรับแต่งอื่น ๆ นั้นไม่ใช่
เปาโลซานตินี

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