การเรียกดูชื่อฟิวด์และประเภทของข้อมูล Shapefile ด้วย Python GDAL/OGR

การเปิดข้อมูล ESRI Personal GeoDatabase (*.mdb files) ผ่านโปรแกรม QGIS

_______หลายคนเคยเพิ่มชั้นข้อมูลแบบ ESRI Shapefile เข้าไปในโปรแกรม QGIS หรือรู้ไม่ว่าโปรแกรม QGIS เองนั่นสามารถเปิดข้อมูล GIS ได้หลากหลายรูปแบบนามสกุล เช่น GeoJSON, GML, KML หรืออื่นๆ ที่ไฟล์นั่นมีการรองรับของ OGR ซึ่งวันนี้ผมก็จะขอยกตัวอย่างการเปิดข้อมูลแบบ ESRI Personal GeoDatabase กันดูครับว่ามีวิธีการอย่างไร


1.เพิ่มชั้นข้อมูลแบบ Vector Layers เลือกแบบ File

2. จากนั้นก็เข้าไปยังไฟล์ ESRI Personal GeoDatabase ที่เราเก็บไว้ แล้วเลือก File type เป็นแบบ ESRI Personal GeoDatabase ดังภาพ

3. เมื่อเพิ่มเข้ามาแล้วก็จะมีหน้าต่างปรากฏขึ้นมาให้เราเลือกข้อมูลที่เราสนใจ แล้วแต่ว่าเราจะเลือกหมดก็ได้ขึ้นอยู่ที่เราเอง ซึ่งไฟล์ mdb เองนั้นสามารถเก็บประเภอของข้อมูลได้หลายประเภท

4. ผลลัพธ์ที่ได้ ดังภาพ

http://www.gdal.org/drv_pgeo.html

การแปลงข้อมูล 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

ดูค่าพิกัดของข้อมูล Shapefile ด้วย ogrinfo

________หลายครั้งที่ผมเคยประสบปัญหาเกี่ยวกับ ค่าพิกัดหรือ projection ของข้อมูล shapefile ถ้าเกิดเราไม่ทราบว่าข้อมูลมันเป็นค่าพิกัดอะไร ซึ่งจะทำให้อาจเกิดความผิดพลาดไปกันใหญ่ เตลิดเปิดเปิงไปไกลเลยทีเดียว ยกตัวอย่างเช่น เรามีข้อมูล shapefile อยู่ชุดนึง แต่ดันไม่มีไฟล์ที่เป็นนามสกุล .prj ซะงั้น ซึ่งเมื่อมันอาจจะทำให้เราไม่ทราบค่าพิกัดข้อมูลชุดนั้น(หรือ shapefile ไม่ได้ define นั้นเอง)ซึ่งผมจะลองเปรียบเทียบ QGIS และ ArcGIS ให้ดูว่าถ้าข้อมูล shapefile ไม่มี .prj จะเป็นอย่างไรครับ แต่ ogrinfo ช่วยเราได้ครับ ซึ่งเราจำเป็นต้องติดตั้ง lib ของ GDAL/OGR ก่อนหรือไม่ก็ติดตั้ง FWTools ครับ ตั้งอย่างเช่น ogrinfo -al -so your.shp เพียงแค่นี้ก็ได้แล้วครับ

เปิดข้อมูล shapefile ใน arcgis

ไม่มีข้อมูลค่าพิกัดของ shapefile

ใน QGIS ก็จะเด้งถามเหมือนกันครับ

นี่แหละของจริงครับ ถามได้ บอกได้

บอกได้ถึงว่ามี Field อะไรบ้างด้วย แจ่ม

เอาไปลองดูกันเด้อครับ พี่น้องไทย...

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