การระบุพิกัด ให้ข้อมูลภาพถ่ายจากดาวเทียม GeoTiff (assign coordinate)

การระบุพิกัด ให้ข้อมูลภาพถ่ายจากดาวเทียม GeoTiff (assign coordinate)

ตั้งแต่ 1 พ.ย. 56 เป็นต้นมา ได้มีโอกาสย้ายตำแหน่งมาทำงานฝ่ายใหม่ ซึ่งได้ทำข้อมูล Raster ตลอด

ช่วงนี้กำลังปรับสีของข้อมูลภาพดาวเทียมอยู่ ที่นี้ปัญหามันก็เกิด เนื่องจากการปรับสี บางครั้งไม่ใช่ปรับใน Software Remote sensing แล้วจะจบ ถ้าใครเคยทำงานจะรู้ว่าบางพื้นที่ ปรับยากมาก พอปรับต้นไม้เห็นชัด พวกบ้าน สิ่งปลูกสร้าง ก็ขาวเวอร์ซะจน มองไม่เห็นรายละเอียดเลย ถ้าปรับให้สิ่งปลูกสร้างเห็นชัด ต้นไม้ก็มืดมองไม่เห็นอีก เลยต้องพึ่งบริการ Photoshop โดยปรกติแล้ววิธีการของผม ถ้าจะเอาภาพแบบ เทพๆ เลยนะ จะใช้วิธีการ ปรับใน Envi ก่อน ภาพแรก ก็ปรับให้ เห็น สิ่งปลูกสร้างชัดเจน แล้ว save ออกมาเป็น 1 ภาพ ต่อมาก็ปรับให้ต้นไม้ชัดอีก บ้านจะขาวสว่างใสก็ช่างมัน save ออกมาเป็นอีกหนึ่งภาพ โดยคลุมโทนสีให้เท่ากัน ถ้าจะออกเหลืองก็ออกทั้งสองภาพ จากนั้นนำมาเปิดใน Photoshop เปิดขึ้นมา แล้วเอาทั้งสองภาพมาซ้อนกัน ถ้าเอา layer ต้นไม้ไว้ชั้นบน layer สิ่งปลูกสร้าง ผมก็จะใช้เครื่องมือยางลบ eraser ลบสิ่งปลูกสร้างที่สว่างจ้าใน layer ต้นไม้ออก ก็จะเห็นสิ่งปลูกสร้างใน Layer สิ่งปลูกสร้างที่ปรับมาสวยงามแทน (ถ้าใครใช้ Photoshop เป็น คงนึกภาพออก) จากนั้นจึง Merge layer แล้ว save

ทีนี้มันมีปัญหาที่เกิดคือ ภาพต้นฉบับมันมี Header ที่เป็นระบบพิกัดติดอยู่ เช่น UTM Zone 47 N พอ save แล้วมันหายไป  ตัวอย่างheader

 

ข้อมูลที่บอก Coordinate System is หายไป  มันจะเหลือแค่  เครื่องหมาย  ‘  ‘ แทน

แต่ข้อมูลพิกัดขอบภาพ Upper Left Lower Right  4 มุม ยังอยู่ (ลืมบอกไปว่าภาพ ต้องมี file .tfw คู่กันด้วยนะครับ ถึงจะมีค่าพิกัด 4 มุมอยู่ )แต่มันจะไม่สามารถนำมาใช้กับ extension image Boundary ที่ใช้สำหรับทำ Index ที่ผมเคยแนะนำมาก่อนหน้านี้ ดังนั้นจึงต้องสรรหาวิธีมาจัดการ จริงๆแล้วใช้ Arcgis ก็สามารถทำได้โดยการเปิดภาพแล้วเลือก Properties จากนั้นกำหนดระบบพิกัดให้มันซะ แต่ถ้ามีภาพหลายร้อย file ก็คงนั่งคลิ๊กกันมือหงิก ก็เลยพยายามหาคำสั่งมาช่วยจัดการ ถ้าเก่งทางด้านเขียนโปรแกรมก็คงง่าย เขียน program ขึ้นมาตัวหนึ่งให้มันรัน Auto ไปเลย

แต่เราเป็นนักภูมิศาสตร์ แค่ใช้ Software เป็น ไม่เก่งขนาดเขียน Code เอง จึงต้องทำงานแบบ Semi -Auto

วิธีการที่เจออันดับแรกเลย คือ Gdalwarp แต่มันต้องสร้าง file ใหม่ขึ้นมาซึ่งไม่สะดวก ช้า เปลืองเนื้อที่เก็บไฟล์

คำสั่งคือ >  gdalwarp -s_srs epsg: 32647 -t_srs epsg:32647 input.tif output.itf

-s_srs epsg:32647 คือบอกมันว่าระบบพิกัดข้อมูลต้นฉบับของเราเป็น UTM Zone 47 N นะ

-t_srs epsg:32647 คือบอกมันว่าระบบพิกัดข้อมูลผลัพธ์ของเราเป็น UTM Zone 47 N นะ

แล้วก็ชื่อ inputfile outputfile

ต่อมาก็เจอวิธีการแบบไม่ต้อง สร้าง file ขึ้นมาใหม่เลย คือการ assign ระบบพิกัดเข้าไปในภาพเดิมเลย

คำสั่งคือ > Gdal_edit -a_srs epsg:32647 input.tif  สั้นๆ แค่นี้แหละ

ใช้คำสั่งวนลูป for %i in (*.tif) do gdal_edit -a_srs epsg:32647 %i

หมายเหตุ ควรจะตรวจสอบให้แน่ใจว่ามีระบบพิกัดอยู่แบบเดียว ใน folder ไม่ผสมกันนะ แล้วจะหาว่าไม่เตือน

Leave a Reply

Your email address will not be published. Required fields are marked *