วิธีเรียก Service จาก ArcGIS Server โดยใช้โปรแกรม QGIS

_______หลายๆ คนเคยอยากที่จะต้องการเรียก service จาก arcgis server หรือ arcgis onine บนโปรแกรม QGIS ซึ่งเราก็สามารถทำได้ 2 วิธีโดยใช้ Python Console และอีกแบบก็จะใช้ GDAL/OGR ครับ
วิธีที่ 1 Python Console
1. In the QGIS Python console ไปที่ Plugins >> Python console
2. จากนั้นก็พิมพ์ qgis.utils.iface.addRasterLayer("http://server.arcgisonline.com/ArcGIS/rest/services/World_Street_Map/MapServer?f=json&pretty=true","esriLayer")

วิธีที่ 2 GDAL/OGR (gdal_translate)
1.ก็พิมพ์ลงใน terminal => gdal_translate "http://server.arcgisonline.com/ArcGIS/rest/services/World_Street_Map/MapServer?f=json&pretty=true" esriStreetLayer.xml -of WMS
2.จากนั้นจะได้ xml มาก็ให้ลากใส่ใน QGIS ได้เลยครับ

Reference : http://hub.qgis.org/wiki/quantum-gis/Arcgis_rest

การ 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)

การตั้งค่าสำหรับเปิดข้อมูลภาพจากดาวเทียมเพื่อให้พร้อมใช้งานในโปรแกรม QGIS

_____หลายท่านคงเคยประสบปัญหาเหมือนๆกันคือว่า เวลาที่เปิดข้อมูลภาพเข้าไปในโปรแกรม QGIS แล้วเรายังต้องมา set ค่าและทำการ enhancement ให้กับข้อมูลภาพอีก ซึ่งเปิดทุกคร้งก็จะต้องทำทุกครั้ง แน่นอนมันก็ต้องเหนื่อยแน่แหละ ยังงั้นเราลองมาดูวิธีที่ทำให้ตอนเราเปิดข้อมูลภาพมาแล้ว สีของภาพชัดเจนพร้อมใช้งานเลย

1. ทำการ Set ค่าดังภาพ โดยไปที่ Setting > Options

2. จากนั้นลองเปิดข้อมูลภาพเข้ามาดูครับ

ขอเชิญเข้าร่วมกลุ่ม QGIS (Quantum GIS) Thailand Community [Facebook]

_____เมื่อหลายสัปดาห์ที่แล้วผมได้สร้างกลุ่มเปิดขึ้นมาเพื่อ แชร์ความรู้เทคนิคใหม่ๆ ทางด้าน gis และรวมถึงท่านที่ติดปัญหาเกี่ยวกับการใช้งานทางด้าน open source gis ก็ช่วยๆกันตอบน่ะครับ อีกอย่างช่วยๆ กันเผยแพร่น่ะครับ กลุ่มเราจะได้กลุ่มใหญ่ๆที่มั่นคงครับ ยังไงก็ช่วยกันติดตามและแผยแพร่ครับ :)

Click เข้าร่วมกลุ่ม

วิธีที่ 3 สำหรับการแก้ปัญหา Encoding ของ Quantum GIS 1.8.0-2 โดยการแก้ไขไฟล์ qgis.bat

_______เมื่ออาทิตย์ที่แล้วผมได้ทดลองการแก้ไขปัญหา Encoding ภาษาไทยที่ทุกคนพบปัญหาบ่อยๆ ครับ จังหวะพอดีพี่ไหลส่งข้อความมาหาว่า "เฮ้ยชัย พี่ได้วิธีแก้ไขปัญหา Encoding อีกวิธีหนึ่งแล้วน่ะ" ก็คือการเข้าไปเพิ่มเข้าไปบรรทัดหนึ่งในไฟล์ qgis.bat (C:\Program Files\Quantum GIS Lisboa\bin\) ซึ่งมีขั้นตอนดังนี้ครับ

1. ภาพที่เรายังไม่ได้ทำการแก้ไขปัญหา encoding ครับ



2. เข้าไปแก้ไขไฟล์ที่ C:\Program Files\Quantum GIS Lisboa\bin\qgis.bat ดังภาพ



3. เปิดไฟล์ qgis.bat ด้วย Notepad++ แล้วเพิ่ม Script เข้าไปดังภาพ
โดยระบุ Script ลงไปดังนี้
SET SHAPE_ENCODING=ISO-8859-11



4. จากนั้นก็ลองเปิดชั้นข้อมูล แล้วเปิดตารางมา รับรองว่าได้ชัวร์ครับ



5. ถึงเราจะต้อง save ไฟล์หรือทำอะไรกับข้อมูล เมื่อเปิดมาอีกครังก็ยังเป็นภาษาไทยเหมือนเดิมครับ

ดูค่าพิกัดของข้อมูล 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 อะไรบ้างด้วย แจ่ม

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

NUGIS Developer Meeting in Doi Angkhang Chiang Mai, January 2013

_______พวกเราได้มีโอกาสพักผ่อนและหาแรงบัลดาลใจในการพัฒนางานที่รับผิดชอบ และได้พูดคุยกันในงานที่เราจะต้องช่วยกันพัฒนาต่อไป

Expression based labeling in Quantum GIS

_______มีหลายคนน่าจะเคยคิดเหมือนกันกับผมว่า ตอนเวลาเราจะแสดง label ของข้อมูลโดยอยากจะแสดงทั้ง 2 fileds ให้เป็น 2 บรรทัด และอาจจะใส่ข้อความบ้างนิดหน่อย ที่นี่มีคำตอบครับ "prov_name" || ' \nเสียชีวิต: ' || "sum" || ' ราย'

Reference: http://planet.qgis.org/planet/tag/labeling/

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

Review Quantum GIS for Windows 8 Enterprise

_______พอทราบข่าวว่า windows 8 ที่มอออกมาให้ใช้ได้แล้ว (ที่จริงทราบตั้งแต่อาทิตย์ที่แล้ว) ผมก็เลยรีบทำการติดตั้งใหม่ซ่ะเลย โดย windows 8 ที่ผม Review นั้นส่วนใหญ่ผมขอในส่วนของ Open Source GIS แล้วกันน่ะครับ ว่ามันมีความเสถียรภาพหรือไม่อย่างไรครับ เรามาดูกันแต่ละรูปเลยดีกว่าครับ

1. ภาพนี้เป็นภพแรกที่เข้ามาใน windows 8 ครับ เมื่อเรา login เข้ามาก็จะพบกับ apps ของ windows 8 ดูเหมือนจะทำมาเพื่อระบบสัผัสโดยเฉพาะเลยน่ะเนี่ย

2. ภาพนี้ก็เป็นการตั้งค่าต่างของระบบครับ ยกตัวอย่างเช่น Personalize ครับ

3. จะเป็น Apps ที่เกี่ยวกับการตรวจสอบสภาพอากาศครับ ที่จริงมีเยอะมากแต่ยกตัวอย่างบางส่วนให้ดูครับ

4. เรามาเริ่ม review ในส่วนโปรแกรมของทางด้าน Open Source GIS กันครับ

5. เปิดโปรแกรม QGIS ใน windows 8 ครับซึ่งเมื่อเข้ามาแล้วลักษณะก็เป็นเหมือนกับ windows 7 นั่นแหละครับ แต่อาจจะดูสวยกว่านิดหน่อยครับ โดยชื่อโปรแกรมจะอยู่ที่ตรงการของ title bar ประมาณนั้นครับ

6. เมื่อเปิดตาราง attribute ก็จะเป็นภาษาไทยด้วยครับ

7. คำนวณ Route ทั้งจังหวัดก็ OK ครับ ไม่มีหน่วงเลย

8. ทำ density ทั้งจังหวัดก็โอเคครับ ไม่มีแฮงค์(มีแต่คนครับ) แต่อาจจะใช้เวลานิดหน่อยครับ

9. เปิดไฟล์ข้อมูลน้ำท่วมทั้งประเทศก็สบายหายห่วงเลย แจ่มมากครับ

สรุปแล้วสำหรับ windows 8 สำหรับโปรแกรมทางด้าน Open Source GIS โดยรวมแล้วผมถือว่าโอเคเลยครับ และที่สำคัญ Theme สำหรับ windows 8 ก็สวยมากด้วยครับ เวลาประมวลผลข้อมูลทางด้าน GIS ก็ถือว่ามีความเสถียรครับ ซึ่งผมยังไม่เจอปัญหาแต่อย่างใดครับ ถ้าท่านไหนที่สนใจ windows 8 อยู่และทำงานทางด้านนี้ก็แนะนำเลยครับ แต่ Driver ของ Note Book บางยี่ห้ออย่างเช่น Dell ของผมเองก็ยังไม่ Driver มาเยครับ ซึ่งภาพที่ผม capture มานั้นเป็นการติดตั้งกับเครื่อง PC ที่ทำงานน่ะครับ :)

Kernel Density Estimation with Quantum GIS [SAGA Toolbox]

_______ใน QGIS 1.7.3 ยังไม่มี Tools ที่เป็นแบบ Kernel Density แต่มาในเวอร์ชั่น 1.8.0 ก็เริ่มมี Tools ที่เป็นแบบ Density เช่น Heatmap ครับ ซึ่งวันนี้ผมก็ได้มีเทคนิควิธีการใหม่ๆ ในการ import tools จาก lib ที่หลากหลายมาแนะนำครับ ยกตัวอย่างเช่น SAGA Algorithms เป็นต้น

Data Example by gistnu

Select Feature Transparent with Quantum GIS

_______เราคงเคยทราบกันดีอยู่แล้วว่าการ select feature โดยใช้ QGIS ลักษณะ feature ที่ถูกเลือกจะเป็นสีเหลือง แต่คราวนี้ผมจะขอแนะนำลูกเล่นใหม่สำหรับการ select feature แบบว่าเมื่อเลือกแล้ว feature ก็จะเป็นลักษณะโปร่งใส (transparent)มาดูวิธีทำกันเลยดีกว่าครับ :)
1. ภาพการ select feature แบบปกติทั่วไป

2. เริ่มการตั้งค่าโดยไปที่ Setting >> Options แล้วเลือกแถบเมนู General

3. จากนั้นให้คลิกเลือกปุ่มสีในส่วนของ Selection Color

4. จะปรากฏหน้าต่าง Selection Color ให้เปลี่ยน Alpha Channel เป็น 0

5. เพียงเท่านี้ก็ select feature ก็จะเป็นโปร่งใสแล้วครับ

ลองไปทำดูกันน่ะครับ \m/