Converting a black border to transparent in geotiff

27 03 2009

ข้อมูลภาพดาวเทียม หรือภาพถ่ายทางอากาศต่างๆ เมื่อทำการปรับแก้ให้ถูกต้องตามระบบพิกัดแล้ว หากต้องการนำขึ้น Google Earth มักจะเกิดปัญหาไม่พึงประสงค์คือจะมีขอบภาพสีดำเกิดขึ้นเนื่องจากภาพไม่อยู่ในแนวสี่เหลี่ยมพอดี ขอบภาพสีดำนี้โดยมากแล้้วค่า digital value นั้นจะมีค่าเป็นศูนย์ (บางโปรแกรมให้เป็น 999) บางจุดในภาพที่เป็นข้อมูลภาพแม้ว่าเราจะเห็นเป็นสีดำ แต่ก็ยังพอมีค่า digital value อยู่บ้างเพียงแต่ตามนุษย์ไม่สามารถมองเห็นความแตกต่างในระดับนั้นได้ เมื่อเกิดขอบภาพสีดำขึ้นในขณะที่นำไปซ้อนทับบน Google Earth ทำให้ขอบที่เป็นสีดำบดบังรายละเอียดด้านหลัง และไม่สวยงาม วันนี้จะมาแก้ปัญหานี้กัน
การแก้ปัญหานี้ ทำได้โดยเราจะทำการแปลงส่วนขอบภาพที่เป็นสีดำนั้นให้โปร่งแสงซะ เมื่อแสดงใน Google Earth ก็จะได้หมดปัญหาเรื่องบดบัง และความสวยงาม

  • ข้อมูลภาพ: Radar image [DAT_01.001]
  • โปรแกรมอรรถประโยชน์สำหรับนัก GIS:  FWTools 2.2.6

ขั้นตอนพอสรุปได้เหลือ 3 ขั้นตอนในการ นำภาพข้อมูล Radar image ขึ้นสู่ Google Earth

  1. ทำการ enhace ภาพข้อมูล Radar
  2. แปลงข้อมูลขอบสีดำให้โปร่งแสง
  3. จัดเตรียมขึ้น Google Earth

มาดูวิธีโดยละเอียดกัน

  1. ทำการ enhance ภาพข้อมูล Radar : [Image Enhancement]
    ใช้ gdalenhance ในการจัดการโดยมีคำสั่งดังนี้
    >gdalenhance -ot byte -equalize DAT_01.001 DAT_01_EQ.oo1
    ขั้นตอนนี้จะได้ผลลัพธ์เป็นไฟล์ DAT_01_EQ.001 

    Equalization

    Equalization

  2. แปลงข้อมูลขอบสีดำให้โปร่งแสง : [Black border transparent]
    ใช้ gdalwarp ในการจัดการโดยมีคำสั่งดังนี้
    >gdalwarp -t_srs epsg:4326 -r bilinear -srcnodata 0 -dstalpha -of GTiff DAT_01_EQ.001 radar_map.tif
    คำอธิบาย: ใช้ gdalwarp แปลงข้อมูลในอยู่ในระบบพิกัด ภูมิศาสตร์ สร้างภาพใหม่ด้วยวิธี bilinear กำหนดค่าที่ไม่มีข้อมูล(ขอบภาพสีดำ)เท่ากับ ศูนย์ และให้ข้อมูลนี้ไปอยู่ในชั้นข้อมูลอัลฟา ผลลัพธ์เป็นไฟล์ geotiff กระทำกับไฟล์ DAT_01_EQ.001 กำหนดผลลัพธ์ชื่อ radar_map.tif 

    transparent black border data

    transparent black border data

  3. จัดเตรียมขึ้น Google Earth : [Prepare for Google Earth]
    ใช้ gdal2tiles ในการจัดการโดยมีคำสั่งดังนี้
    >gdal2tiles -title “RADAR” radar_map.tif radar
    คำอธิบาย: gdal2tiles จะทำ tilling เพื่อแบ่งข้อมูลเป็นส่วนย่อยๆเพื่อให้รวดเร็วในการเปิดบน Google Earth หรือ Google Maps กระทำกับไฟล์ radar_map.tif และให้ผลลัพธ์ อยู่ในโฟล์ดเดอร์ชื่อ radar 

    Gdal to tiles

    Gdal to tiles

เมื่อผ่านทั้ง 3 กระบวนการแล้ว ผลลัพธ์ที่ได้จะอยู่ในโฟล์ดเดอร์ radar จะมีไฟล์ที่ชื่อ doc.kml เมื่อเรียกไฟล์นี้ใช้งาน จะทำการเปิด Google Earth และแสดงภาพโดยที่ขอบสีดำจะถูกทำให้โปร่งแสงสามารถเห็นข้อมูลที่อยู่ด้านหลังได้


Actions

Information

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s




%d bloggers like this: