Read Raster Data Using GDAL

วันนี้จะมาแนะนำวิธีการอ่านข้อมูล raster ที่มีข้อมูลเชิงตำแหน่ง เช่น GeoTiff ,Tiff เป็นต้น โดยเครื่องมือที่ชื่อว่า gdallocationinfo โดยเราสามารถค้นหา ข้อมูลของ Raster จากตำแหน่งทางภูมิศาสตร์ได้เลย

ติดตั้ง GDAL

การติดตั้ง GDAL สามารถทำการติดตั้งได้ด้วยตัว python package manager ทั่วไป เช่น PIP Conda Mini Conda เป็นต้น

conda install -c conda-forge gdal
pip install gdal 

ตรวจสอบการติดตั้งด้วยคำสั่ง

gdal-config --version

>>> 3.5.3 เป็นต้น

เมื่อทำการติดตั้งสำเร็จการใช้งาน GDAL สามารถเรียกใช้งานได้ทั้งผ่าน command line และ บน python script โดยรันผ่าน OS หรือ subprocess ก็ได้เช่นกัน

ตัวอย่างการใช้งาน

#command line
gdallocationinfo -valonly input_raster.tif 50.3,25.5

#python
subprocess.run(['gdallocationinfo', '-valonly', raster_path_dsm, '-wgs84', str(longitude), str(latitude)], stdout=subprocess.PIPE)

#result 
pixel value

คำสั่ง gdallocationinfo ยังมีตัวเลือกอีกมากมายให้ปรับใช้สามารถค้นหาข้อมูลเพิ่มเติมได้ที่

https://gdal.org/programs/gdallocationinfo.html