Create color relief with GDAL

______เมื่อสร้าง Hillshade แล้วก็มาสร้าง color relief กันดูบ้างครับ ซึ่งวิธีการไม่ได้ยากอะไรเลยครับ ตัวอย่าง gdaldem color-relief DEM.tif elevation.txt colordem.tif

Create Hillshade by gdaldem

_____สร้าง Hillshade ง่ายๆ โดยใช้ GDAL-based DEM utilities ตัวอย่าง gdaldem hillshade dem.tif shade.tif

Enhance image with FWTools

_______ผมได้ลองใช้ FWTools ในการ enhance ข้อมูลภาพดาวเทียม โดยใช้ชุดคำสั่ง gdalenhance ซึ่งผมได้วิธีทำแบบนี้มาจากสื่อการสอนของพี่โอ๊ต sawarinO ครับ แล้วผมก็ลองเอามาทำดูบ้างก็น่าสนใชดีครับ ซึ่งวิธีการก็ไม่ยากเลยครับ แต่ปัจจุบันนี้มันก็มี software ที่เราไม่ต้องพิมพ์ให้มันยากขนาดนั้นครับ แต่ผมชอบแบบ hardcore มากกว่าครับเพราะว่ามันรู้สึกสนุกดีครับในการเขียน command line อย่าง FWTools ซึ่งก็ต้องขอบคุณบล๊อกของพี่โอ๊ตมากเลยครับที่ได้แนะนำการใช้ชุดคำสั่ง gdal/ogr ครับ

ภาพก่อนที่จะทำการ enhance

ชุดคำสั่งการใช้ gdalenhance

ภาพที่ได้หลังจากการ enhance

ภาพที่ได้หลังจากการ enhance

Reference: http://www.gdal.org & http://sawarino.blogspot.com

Clip Data ESRI Shapefiles using FWTools [ogr2ogr]

______สำหรับคนที่ชอบความ Hardcore แบบผมครับ ที่จริงการตัดข้อมูลไม่จำเป็นต้องทำแบบผมก็ได้ เราก็แค่ใช้ฟังก์ชั่นที่เค้ามีให้เลยก็ได้ครับ แต่ถ้าจะลองทำแบบผมก็ตามนี้เลยครับ ogr2ogr -clipsrc plk_prov.shp newflu52.shp flu52.shp

Reference: http://www.gdal.org/ogr2ogr.html
Data Example: NUGIS

Conversion from MapInfo Tab file to ESRI Shapefile using FWTools [ogr2ogr]

_______ช่วงนี้ยังคงอยู่กับ gdal/ogr ครับ และยังทำงานอยู่กับ ogr2ogr เพื่อที่แปลงไฟล์จาก MapInfo ให้ไปเป็น ESRI Shapefile เคยมีน้องที่ทำงานอยู่บริษัทเอกชน หลายคนถามบ่อยๆ ว่าผมอยากจะเอาไฟล์จากโปรแกรม MapInfo ไปเปิดใน ArcGIS หรือ QGIS ได้อย่างไรครับ ซึ่งผมก็พยายามที่จะให้น้องเค้าทำตามนี้เลยครับโดยใช้ FWTools แล้วให้ใช้ชุดคำสั่ง ogr2ogr ครับ ยกตัวอย่างเช่น ogr2ogr -f "ESRI Shapefile" plkprovince.shp Province_phitsanulok.TAB หรืออีกตัวอย่าง ogr2ogr -f "ESRI Shapefile" new_mexico.shp Mexico.TAB ถ้าต้องการระบุค่าพิกัดใหม่ก็ตามนี้เลย ogr2ogr -f "ESRI Shapefile" -t_srs "EPSG: 4269" Mapinfo2Shp\LincsBoundary.shp Mapinfo2Shp\Lincolnshire.tab

Reference: http://www.mapwindow.org/phorum/read.php?4,21913
Data Example: http://gw.vtrenz.net/?NYEWEN5YSB