Geoserver with Image mosaicking plugin

_______สำหรับผู้ที่สนใจในงานด้าน Internet GIS หลายคนคงเคยคิดว่าจะทำอย่างไรดี จะให้ข้อมูลที่เป็น basemap อย่างเช่น ข้อมูลภาพดาวเทียม ภาพถ่ายทางอากาศ ทั่วทั้งพื้นที่ที่เราต้องการ เพราะว่าไฟล์ขนาดใหญ่มากจึงทำให้มีปัญหาในการต่อภาพ (mosaic) ดังนั้นผมจึงแนะนำอีกวิธีหนึ่งที่ง่ายในการทำโดยใช้ "Image mosaicking plugin" ของ GeoServer โดยการทำดังต่อไปนี้

1. เตรียมข้อมูลภาพที่ต้องการจะแสดงผล และจะเปิดข้อมูลนั้นด้วย QGIS :)

2. Copy Data ไปยัง dir ของ Geoserver

3. ทำการสร้าง Layers >> โดยใช้ ImageMosaic

4. ทำการ Publish

5. Click Layer Preview

Download Data Example : Click
Reference:http://docs.geoserver.org/latest/en/user/tutorials/image_mosaic_plugin/imagemosaic.html

การ 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 1.9 และ QGIS 2.0 สามารถเปิดข้อมูลภาพ ไฟล์นามสกุล ECW ได้

_______ก่อนอื่นเลยก็ต้องแนะนำประเภทของข้อมูลภาพนี้ก่อนครับคือ ECW (Enhanced Compression Wavelet) ซึ่งพัฒนามาจากบริษัท ERDAS แต่ตอนนี้เป็นของ Intergraph แล้วครับ เท่าที่ผมได้ใช้โปรแกรม QGIS มาน่ะครับ ตั้งแต่เวอร์ชั่น 1.8 หรือ 1.9 ขึ้นไป จะไม่สามารถเปิดข้อมูลภาพ ECW นี้ได้ แต่วันนี้ผมมีวิธีการที่ทำให้สามารถเปิดได้แล้วครับ ดังนี้

0. แสดงภาพที่ Error ที่เกิดขึ้นครับ

1. ให้ดาวน์โหลดไฟล์ gdal-19-1600-ecw จากลิงค์นี้ Click  หรือเข้าไปโหลดที่ http://www.gisinternals.com/sdk/

2. เมื่อดาวน์โหลดแล้วก็ดับเบิลคลิ๊กไฟล์นั้นเลย แล้วคลิกปุ่ม Next

 
3. คลิกปุ่ม accept แล้วคลิกปุ่ม Next

4. คลิกที่ปุ่ม Complete แล้วคลิกปุ่ม Next
 

5. คลิกปุ่ม Install เริ่มการติดตั้งโปรแกรม

  

6. คลิกปุ่ม Finish สิ้นสุดการติดตั้ง

7. เมื่อติดตั้งเสร็จแล้วให้ไปที่โฟลเดอร์ C:\Program Files\GDAL\ ก็จะเห็นไฟล์ libecwj2.dll ให้เรา copy

8. แล้วนำมาวางไว้ที่ C:\Program Files\QGIS Dufour\bin ดังภาพ

9. จากนั้นให้ไปที่ C:\Program Files\GDAL\gdalplugins ก็จะเห็นไฟล์ gdal_ECW_JP2ECW.dll ให้เรา copy

10. แล้วนำมาวางไว้ที่ C:\Program Files\QGIS Dufour\bin\gdalplugins\1.9 ดังภาพ

11. ปิด-เปิดโปรแกรมใหม่ ก็จะสามารถเปิดไฟล์ ECW ได้แล้วครับ

 

 Reference: http://www.faunalia.com

การแปลงข้อมูลภาพไฟล์นามสกุล 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

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

สร้างเส้นชั้นความสูงจาก 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/