การ Capture ข้อมูลภาพจาก Google Map มีค่าพิกัดติดมาด้วย โดยใช้โปรแกรม Quantum GIS

_______การ Capture ข้อมูลภาพจาก google earth หรือ google maps นั้นสามารถทำได้หลากหลายวิธีมากเลยครับ แต่บางวิธีมันก็ยากและจะต้องนำภาพที่ capture มาต่อ(mosaic) ก่อนจึงจะนำมาตรึงค่าพิกัด (georeferencer) ซึ่งหลายท่านน่าจะคิดออกน่ะครับว่า ยกตัวอย่างเช่น เราจะ save ภาพจาก google earth มาจะเป็นแค่ภาพธรรมดามันไม่มีค่าพิกัดและทำให้เราต้องเสียเวลาเป็นอันมากในการทำตรึงค่าพิกัด (georeferencer) และวันนี้ผมมีวิธีการที่ง่ายๆ โดยใช้เจ้า QGIS นี่แหละครับ ซึ่งเรา save ภาพก็จะมีค่าพิกัดติดมาด้วยทันทีครับ จึงทำให้ไม่ต้องเสียเวลามาตรึงค่าพิกัดครับ โดยมีวิธีการดังนี้

1. เปิดโปรแกรม QGIS แล้วไปที่เมนู Setting >> Project Properties... แล้วทำการเลือก Projection เป็น Google Mercator (EPSG:900913)

2. จากนั้นให้เราใช้ Plugin "Openlayers Plugin" (ถ้าใครที่ยังไม่ได้ติดตั้ง Plugin ก็ให้ทำการติดตั้งด้วยน่ะครับ)

3. ถ้าท่านใดที่มี Block ในพื้นของตนเองอยู่แล้วก็ให้เอามาซ้อนทันได้เลยครับ เพื่อที่จะได้ Capture ได้ตรงตามที่ต้องการครับ

4. จากนั้นก็จะเริ่ม Capture ข้อมูลภาพ โดยไปที่เมนู File >> Save as Image... (สามารถเลือกประเภทของไฟล์ภาพหลายนามสกุลครับ ซึ่งผมเลือกเป็นแบบ *.tif)

5. ทำดั้งข้อที่ 4 และเลื่อนพื้นที่ให้ครบตามพื้นที่ที่ต้องการ ดังภาพ

6. ภาพทั้งหมดที่ Capture โดยเปิดขึ้นมาครับดังภาพครับ

7.ถ้าท่านไหนที่ต้องการจะ mosaic image หรือจะ Warp Projection ให้เป็น Projection ที่ต้องการก็แนะนำให้ใช้ GDAL เลยครับ โดยเรียกใช้จาก lib ของ QGIS ก็ได้ หรือจะใช้ FWTools ก็ดีครับ

Download Data Example (Powered by Google Drive)

การแปลงข้อมูล ESRI Shapefile (*.shp) ไปเป็นข้อมูลแบบ GeoJSON โดยใช้ ogr2ogr

_______การแปลงข้อมูลไปเป็น GeoJSON นั่นสามารถทำได้หลากหลายวิธีมากครับ ไม่ว่าจะใช้โปรแกรมแปลง เช่น QGIS หรือสามารถที่จะใช้บนเว็บออนไลน์ก็ได้เช่นกัน แต่วันนี้ผมขอแนะนำการแปลงข้อมูลโดยใช้ ogr2ogr ซึ่งเป็นหนึ่งในไลบราลี่ของ GDAL/OGR(Geospatial Data Abstraction Library)  ซึ่งท่านสามรถใช้ผ่านโปรแกรม OSGeo4W Shell หรือ FWTools ได้เลยครับ

1. ตัวอย่างข้อมูลแบบ ESRI Shapefile

2. ทดสอบว่า ogr2ogr ใช้งานได้หรือไม่

3. จากนั้นเราต้องพิมพ์คำสั่งลงไปดังด้านล่างนี้
>>> $ ogr2ogr -f GeoJSON outputfile.json inputfile.shp
ท่านไหนที่ต้องการแปลงระบบพิกัดด้วยก็ระบุ -t_srs แล้วตามด้วย EPSG Code เช่น -t_srs EPSG:4326 ครับ

4. ผลที่ได้จากการแปลง ESRI Shapefile เป็น GeoJSON ครับ

สามารถศึกษาเพิ่มเติมได้ที่ http://www.gdal.org/ogr2ogr.html

การแปลงข้อมูลภาพไฟล์นามสกุล MrSID ให้ไปเป็น GeoTIFF โดยใช้ FWTools และ QGIS

_____อย่างที่ทราบกันดีครับว่า ส่วนใหญ่แล้วข้อมูลภาพถ่ายทางอากาศนั้นจะเป็นไฟล์แบบ MrSID ซึ่งเป็นไฟล์ที่ผ่านการบีบอักข้อมูลให้เหลือขนาดไฟล์อยู่แค่นิดเดียว เพื่อสะดวกต่อการจัดเก็บและส่งให้กันได้ง่าย แต่ในขณะทีบางโปรแกรมนั้นไม่อาจสามารถเปิดไฟล์ประเภทนี้ได้ (แต่ QGIS เปิดได้น่ะครับ) ดังนั้นเราจึงต้องแปลงกลับไปให้อยู่ในรูปแบบไฟล์ที่สามารถใช้ได้หลากหลายโปรแกรม อย่างเช่น GeoTiff เป็นต้น แต่การแปลงกลับมานี่คือจะต้องดูพื้นที่ในเครื่องเราสักหน่อยน่ะครับ ถ้าแปลงไฟล์หลายๆไฟล์นี้ก็กินพื้นที่ในเครื่องเราพอสมควรครับ ซึ่งโปรแกรมที่ผมจะแนะนำในวันนี้มีอยู่ 2 ตัวคือ FWTools และ QGIS ใชัได้ทั้งสองแบบเลย มาเริ่มกันดีกว่าครับ

FWTools
1. เปิด FWTols แล้วใช้คำสั่ง gdal_translate -of GTiff image.sid image.tif

2. ผลลัพธ์จาก openev

QGIS (Quantum GIS)
1. Raster >> Conversion >> Translate (Convert format)

2. ระบุพารามิเตอร์ดังภาพ

การแปลงข้อมูลภาพไฟล์นามสกุล ECW ให้ไปเป็น GeoTIFF โดยใช้ FWTools เพื่อให้สามารถเปิดกับโปรแกรม QGIS ได้

_______ผมได้ลองใช้ QGIS 1.8.0 เปิดข้อมูลภาพถ่ายทางอากาศที่มีรูปแบบไฟล์เป็น *.ecw แล้วเปิดไม่ได้ ซึ่งไฟล์ *.ecw ประเภทนี้มีชื่อเต็มว่า "Enhanced Compression Wavelet" หรือ "ERDAS Compress Wavelets" เป็นการบีบอัดไฟล์อีกแบบหนึ่งซึ่งถูกพัฒนาโดย ERDAS หลายๆ ท่านอาจคงรู้จักดีในโปรแกรม ERDAS IMAGINE แต่ปัจจุบันนี้บริษัท Intergraph ได้ซื้อไปแล้ว เรามาเข้าประเด็นกันดีกว่าครับฮ่าๆ เมื่อก่อนนั้นตอน QGIS 1.7.4 ยังสามารถเปิดได้อยู่ (ถ้าใครสนใจที่จะ Downgrade ลงไปน่ะครับ) ซึ่งผมเข้าใจว่าน่าจะเป็นที่การ Config GDAL ในเวอร์ชั่น 1.9 จึงทำให้เกิดปัญหานี้ แต่เรามีวิธีแก้ไขปัญหานี้โดยการใช้ FWTools โดยการแปลงไฟล์ให้ไปเป็นไฟล์อื่นแทนเช่น GeoTIFF ครับ โดยการเรียกใช้ gdal_translate Example: gdal_translate -of GTiff image.ecw image.tif (แต่เมื่อเราแปลงไฟล์เสร็จแล้วขนาดไฟล์จะมีขนาดใหญ่ขึ้นน่ะครับ)

1. ใช้ QGIS 1.8.0 เปิดไฟล์ *ecw แล้วเปิดไม่ได้

2. จงแปลงไฟล์ไฟล์เป็น GeoTIFF ด้วยสุดยอดโปรแกรมเทพด้วย FWTools ซ่ะ

3. ผลลัพธ์ที่ได้จากการแปลงไฟล์

4. ผลลัพธ์ที่ได้จากการแปลงไฟล์

5. ข้อพึงระวังเมื่อเราแปลงไฟล์ได้แล้ว ไฟล์ใหม่ที่ได้มาจะมีขนาดใหญ่ขึ้นดังภาพ สาเหตุเพราะว่า *ecw นั้นเป็นไฟล์ที่บีบอัดภาพมาให้ขนาดเล็กลงครับ

Powered by gistnu

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

Conversions from PostGIS to KML formats [ogr2ogr]

_______ห่างหายกันไปนานครับช่วงนี้ ผมพยายามที่จะแปลงประเภทของข้อมูลจากที่อยู่ใน PostgreSQL/PostGIS ให้เป็น KML ที่สามารถเปิดในโปรแกรม Google Earth เรามาลองกันเลยครับซึ่งขั้นแรกเราต้องมีโปรแกรม FWTools โดยใช้ ชุดคำสั่ง ogr2ogr ครับ Example: ogr2ogr -f "KML" c:\workspace\chingchai.kml PG:"host=localhost user=postgres dbname=postgis20 password=****" -sql "select * from prov" จากนั้นเราก็จะได้ไฟล์ kml มาก็จะสามารเปิดกับ Google Earth ได้ครับ

Powered by NUGIS

สร้างเส้นชั้นความสูงจาก DEM โดยใช้ FWTools [gdal_contour]

_____เมื่อก่อนนั้น...เคยคิดว่าการทำเส้นชั้นความสูงเป็นเรื่องที่ยากมาก แต่เดี๋ยวนี้มันเป็นอะไรที่ง่ายและสะดวกมากเลยครับ ซึ่งวันนี้ผมก็ขอแนะนำโปรแกรมอีกตัวหนึ่งที่มีชื่อว่า "FWTools" ซึ่งเป็น Open Source GIS โดยการรวมเอาเครื่องมือและชุดคำสั่ง ไลบรารี่ต่างๆ ได้แก่ penEV, GDAL, MapServer, PROJ.4 และ OGDI เพื่อเป็นเครื่องมือสำหรับจัดการข้อมูลเชิงพื้นที่ครับ และผมก็ขอแนะนำการสร้างเส้นชั้นความสูง โดยการเรียกใช้ไลบรารี่ของ GDAL ที่ชื่อว่า gdal_contour ครับ โดยตัวอย่างตามนี้เลยครับ gdal_contour -a elev dem.tif contour.shp -i 10.0

1. สร้างเส้นชั้นความสูงจาก DEM ช่วงชั้นความสูงเท่ากับ 20 เมตร และกำหนดให้สร้าง Attribute มีชื่อว่า ELEV


2.ภาพหลังจากการประมวลผลเสร็จแล้ว


3.จงเปิดมันขึ้นมาด้วยโปรแกรม OpenEV ดังภาพ



5. และที่สำคัญจงเปิดเส้นชั้นความสูงด้วย Quantum GIS ด้วยน่ะครับ


ศึกษาข้อมูลเพิ่มเติมได้ที่
http://fwtools.maptools.org/
http://www.gdal.org/gdal_utilities.html

แสดงรายละเอียดข้อมูลภาพด้วย gdalinfo

____เนื่องจากข้อมูลภาพหรือข้อมูล Raster บางที่เราจำเป็นที่จะต้องทราบรายละเอียดของข้อมูล เช่น Coordinate System, Pixel Size, Corner Coordinate และรายละเอียดต่างๆของข้อมูลภาพ โดยปกติแล้วผมก็จะใช้โปรแกรม FWTools ซึ่งจะมี Library GDAL อยู่โดยใช้คำสั่งตามนี้ gdalinfo image.tif หรือไม่ก็ใช้ QGIS ก็ได้เหมือนกันครับโดยไปที่เมนู Raster > Miscellaneous > Information ซึ่งท่านที่สนใจสามารถลองนำไปใช้ได้ครับ





Reference: http://www.gdal.org/