PostGis-Convert Lat,Lng (4326) to Geom UTM wgs84 zone 47 (32647)

UPDATE place SET the_geom = ST_Transform(ST_SetSRID(ST_MakePoint(cast(lng as double precision), cast(lat as double precision)), 4326),32647);

HospitalOS GIS

ระบบสารสนเทศภูมิศาสตร์เพื่อบันทึกพิกัดหลังคาเรือนสำหรับหน่วยงานที่ใช้โปรแกรมของ HospitalOS (HospitalOS GIS)

รายละเอียด

- โปรแกรม GIS for HospitalOS พัฒนาต่อจาก อ.ไกรวุฒิ แก้วชาลุน เวอร์ชั่น original http://www.kriwoot.com/gisjhcis/ โดยสามารถบันทึกพิกัดหลังคาเรือน และสืบค้น จากชื่อ-สกุล รหัส ปชช. สืบค้นประชากรตามกลุ่มอายุ และการหาหลังคาเรือนจากรัศมีที่กำหนดได้ เท่านั้น เป็นโปรแกรมส่วนขยายเพิ่มทำให้สามารถทำการบันทึกพิกัดหลังคาเรือน สำหรับหน่วยงานที่ใช้โปรแกรม ของ HospitalOS

- ตารางที่เกี่ยวข้องกับการทำงานคือ t_health_family (ข้อมูลบุคคล) และ t_health_home (ข้อมูลหลังคาเรือน) โดยจะต้องเพิ่มข้อมูล lat และ lng ในตาราง t_health_home (รายละเอียดตามคู่มือ)

- ทดสอบการทำงาน Web Server ใช้ Mapserver for Window (MS4W) หรือจะใช้ Appsev หรืออื่นๆ ก็ได้ ฐานข้อมูลที่ทดสอบระบบเป็น PostgreSQL 9.0 Web Browser Google Chrome รุ่น 22.0.1229.94

- Username และ Password ใช้ gis เหมือนกัน

- มีปัญหาการติดตั้งโปรแกรมหรือมีข้อสงสัย mail มาที่ ajarn.veerapong@gmail.com


คู่มือการใช้งาน || Source Code โปรแกรม

HosxP GIS

 ระบบสารสนเทศภูมิศาสตร์เพื่อบันทึกพิกัดหลังคาเรือนสำหรับหน่วยงานที่ใช้โปรแกรมของ HosxP (HosxP GIS)


- โปรแกรม GIS for HOSxP พัฒนาต่อจาก อ.ไกรวุฒิ แก้วชาลุน เวอร์ชั่น original http://www.kriwoot.com/gisjhcis/ โดยสามารถบันทึกพิกัดหลังคาเรือน และสืบค้น จากชื่อ-สกุล รหัส ปชช. สืบค้นประชากรตามกลุ่มอายุ และการหาหลังคาเรือนจากรัศมีที่กำหนดได้ เท่านั้น เป็นโปรแกรมส่วนขยายเพิ่มทำให้สามารถทำการบันทึกพิกัดหลังคาเรือน สำหรับหน่วยงานที่ใช้โปรแกรม ของ HOSxP

- ตารางที่เกี่ยวข้องกับการทำงานคือ person(ข้อมูลบุคคล) และ house(ข้อมูลหลังคาเรือน) โดยจะใช้ข้อมูล latitude และ longitude ในตาราง house (รายละเอียดตามคู่มือ)

- ทดสอบการทำงาน Web Server ใช้ Mapserver for Window (MS4W) หรือจะใช้ Appsev หรืออื่นๆ ก็ได้ ฐานข้อมูลที่ทดสอบระบบเป็น Mysql (character-set = utf8) และ Web Browser Google Chrome รุ่น 22.0.1229.94

- Username และ Password ใช้ gis เหมือนกัน

- มีปัญหาการติดตั้งโปรแกรมหรือมีข้อสงสัย mail มาที่ ajarn.veerapong@gmail.com


คู่มือการใช้งาน || Source Code โปรแกรม

การติดตั้งโปรแกรม ArcView 3.3 บน Window XP, Vista, และ 7 แบบ 32 และ 64 bit

Download ไฟล์ทั้งหมดจาก  www.ajarnveerapong.com/download/arcview32&64bit.zip
ทำ การแตก zip แล้วทำตามคู่มือ หรือคลิกขวาเลือกเปิด  install_arcview_32Bit.bat สำหรับ 32 bit หรือ  install_arcview_64Bit.bat สำหรับ 64 bit ครับ

FWTools การ merge ภาพหลายๆ ไฟล์ด้วย gdal_merge

ที่มา:http://sawarino.blogspot.com/2009/07/gdalmerge.html
            การใช้งาน gdal_merge อย่างง่าย โดยจะขอทดลองกับข้อมูล GDEM โดยการใช้งาน gdal_merge นั้นจะแตกกต่างกับคำสั่ง gdal อื่นๆ คือ ต้องระบุผลลัพธ์ (Output file) ก่อนแล้วจึงระบุไฟล์ต้นฉบับ การใช้งานเบื้องต้นเป็นไปตาม syntax ดังนี้ครับ
#gdal_merge [option] [output file] [input files] ซึ่งตัวอย่างการใช้งานดังนี้ครับ
#gdal_merge -of GTiff -o output.tif abc1.tif abc2.tif abc3.tif abc4.tif
จากตัวอย่างดังกล่าว จะเห็นว่าการระบุไฟล์ต้นฉบับที่มากกว่า 1 ไฟล์นั้นค่อนข้างยุ่งยาก เราสามารถใช้เครื่องหมาย * เข้ามาช่วยในคำสั่งได้ดังนี้ครับ
#gdal_merge -of GTiff -o sample.tif *.tif

Resampling ภาพด้วย gdalwarp

ที่มา:http://sawarino.blogspot.com/2011/08/resampling-gdalwarp.html

การ resampling ข้อมูลภาพเพื่อให้มีขนาดหรือจำนวนจุดภาพเท่าที่ต้องการ ตัวอย่างเช่น ภาพเดิมมีขนาด 300x300 จุดภาพ แต่ต้องการเพิ่มให้เป็น 1200x1200 จุดภาพ เรามาดูวิธีการกันเลยนะครับ
คำสั่งหรือโปรแกรมประยุกต์ที่ใช้คือ gdalwarp ส่วน option หรือ ทางเลือกที่จะใช้คือ -ts width height ดังนี้ครับ

>>>gdalwarp -ts width height inputfile outputfile

ตัวอย่างการใช้งานมีดังนี้ครับ

>>>gdalwarp -ts 1200 1200 abc300.tif abc1200.tif

ปล. หากต้องการระบุวิธี resampling ให้เพิ่มทางเลือก -r xxxx นะครับ เช่น -r bilinear, cubic, cubicspline หรือ lanczos. แต่ถ้าไม่ระบุโปรแกรมจะใช้วิธี nearest เป็นค่าเริ่มต้นครับ

การแปลงไฟล์แผนที่ GIS จาก MapInfo เป็น Shapefile

           การแปลงไฟล์แผนที่ GIS จาก MapInfo เป็น shapefile จริงแล้วจะต้องใช้เครื่องมือที่มีอยู่ในโปรแกรม MapInfo หรือโปรแกรมเสริที่จะรับไฟล์ดังกล่าวมาแสดง อย่างเช่นใน ArcView เป็นต้น แต่ก็ยังมีโปรแกรมเฉพาะเล็กๆที่ไม่ต้องพึ่งโปรแกรมดังกล่าว โดยสามารถแปลงได้ด้วยอย่างดายและสะดวกได้แก่ ogr2gui 
          Download จาก http://www.ajarnveerapong.com/download/ogr2gui.zip

การแปลงค่าพิกัดจากองศา ลิปดา ฟิลิปดา เป็น แบบองศาทศนิยม

 ที่มา :
http://vconsultgis.blogspot.com/2009/06/blog-post.html

การแปลงค่าพิกัดจาก องศา ลิปดา ฟิลิปดา (hddd?ddd’ddd.ddd’’) และ องศา ลิปดา ฟิลิปดา (hddd?ddd.ddd’) 

เป็น แบบองศาทศนิยม (hddd.ddd?)

ตัวอย่างที่ 1 องศา ลิปดา ฟิลิปดา (hddd?ddd’ddd.ddd’’)

Lat  13? 45’ 53.7’’ N

Long  100? 32' 17.9’’ E

ขั้นตอนที่ 1 หาค่าลิปดา นำค่าลิปดาหารด้วย 60 จะได้

Lat  13? 45’ 53.7’’ N

ลิปดา = 45/60

ลิปดา = 0.75

Long  100? 32' 17.9’’ E

ลิปดา = 32/60

ลิปดา = 0.533

ขั้นตอนที่ 2  หาค่าฟิลิปดา นำค่าฟิลิปดาหารด้วย 3600 จะได้

Lat  13? 45’ 53.7’’ N

          ฟิลิปดา = 53.7/3600

          ฟิลิปดา = 0.014916

Long  100? 32' 17.9’’ E

ฟิลิปดา = 17.9/3600

          ฟิลิปดา = 0.004972

ขั้นตอนที่ 3 นำค่าลิปดาจากขั้นตอนที่ 1 +  ค่าฟิลิปดาจากขั้นตอนที่ 2

Lat  13? 45’ 53.7’’ N

ลิปดา = 45/60

ลิปดา = 0.75

Lat  13? 45’ 53.7’’ N

          ฟิลิปดา = 53.7/3600

          ฟิลิปดา = 0.014916

Lat 13? 45’ 53.7’’ N

= 13? + (0.75+0.014916)

                   = 13? + 0.764916

ดังนั้น       Lat  = 13.764916? N

Long  100? 32' 17.9’’ E

ลิปดา = 32/60

ลิปดา = 0.533

Long  100? 32' 17.9’’ E

ฟิลิปดา = 17.9/3600

          ฟิลิปดา = 0.004972

Long  100? 32' 17.9’’ E

= 100? + (0.533+0.004972)

                   = 100? + 0.537972

ดังนั้น       Long  = 100.537972? E

สรุป

Lat  13? 45’ 53.7’’ N               =          Lat  13.764916? N

Long  100? 32' 17.9’’ E           =          Long  100.537972? E

มีอะไรใน worldfile

ไฟล์ GeoTiff เป็นไฟล์ ที่รู้จักกันดี ในหมู่นักทำแผนที่ และ ภูมิสารสนเทศน์
เป็นไฟล์ ที่เก็บ ภาพ ได้ หลายชั้นข้อมูล พร้อมทั้ง มีข้อมูลเชิงตำแหน่ง ด้วย
ทั่วไป ก็มี ด้วยกัน สองแบบ

1. GeoTiff แบบที่มี worldfile แบบนี้ จะมีด้วย กัน สอง ไฟล์ ไปไหนไปด้วยกัน ไม่งั้นทำงานไม่ได้ ไฟล์แรก เก็บภาพ นามสกุล .tif อีกไฟล์ เป็น ไฟล์ ที่เก็บข้อมูลเชิงตำแหน่งของภาพ นามสกุล .wtf เรียกว่า world file
2. GeoTiff แบบที่ มี world file อยู่บนหัวไฟล์ tiff เลย แบบนี้ มีไฟล์เดียว เวลาใช้งาน ก็มาอ่านค่า จากส่วนหัวนี้ก่อน ไปเป็นข้อมูลเชิงตำแหน่ง

ทีนี้ ลองมาดูว่า worldfile เก็บค่าอะไรบ้าง
worldfile มีอยู่ด้วยกัน 6 บรรทัด เช่น (สามารถเปิดอ่านด้วย text editor ทั่วไปได้)
20.17541308822119
0.00000000000000
0.00000000000000
-20.17541308822119
424178.11472601280548
4313415.90726399607956
ความหมายของ ทั้ง 6 บรรทัด มีดังนี้

line 1: x-dimension of a pixel in map units
line 2: rotation parameter
line 3: rotation parameter
line 4: NEGATIVE of y-dimension of a pixel in map units
line 5: x-coordinate of center of upper left pixel
line 6: y-coordinate of center of upper left pixel
ในการใช้งาน
ตัวอย่างเช่น โปรแกรม Arcview เรียกข้อมูลภาพ ขึ้นมา นั่นหมายถึง กำลังนำภาพ เข้าสู่ระบบพิกัดโลก ที่เรียกว่า การฉายแผนที่ ซึ่งต้องอาศัย พารามิเตอร์จำนวน 6 ตัว ในสมการ affine transformation

x1 = Ax + By + C
y1 = Dx + Ey + F
โดยที่
x1 = ค่าพิกัด แกน x ของภาพ บน แผนที่
y1 = ค่าพิกัด แกน y ของภาพ บน แผนที่
x = จำนวน column ของจุดภาพ ในภาพ
y = จำนวน row ของจุดภาพ ในภาพ
A = มาตราส่วนแกน x
B, D = เทอมการหมุน
C, F = แปลงค่าพิกัดจากแผนที่ จากแบบ จากจุดภาพบนซ้าย (การนับจำนวนจุดภาพ ของจอภาพ นับแบบ แกนx จากซ้ายมาขวา แกนy จากบนมาล่าง)
E = ค่าลบ มาตราส่วนแกน y

ดังนั้น ค่าต่างๆ ใน worldfile สามารถแปลงมาเป็น พารามิเตอร์ทั้ง 6 ได้ ดังนี้
The transformation parameters are stored in the world file in this order:
20.17541308822119 – A
0.00000000000000 – D
0.00000000000000 – B
-20.17541308822119 – E
424178.11472601280548 – C
4313415.90726399607956 – F
เมื่อ อ่าน ข้อมูลใน worldfile เข้าไปในโปรแกรม ก็จะสามารถ ทราบได้ ว่า ภาพ นั้น อยู่ตรงส่วนไหนของแผนที่ หรืออยู่ตรงส่วนใดของโลกนั่นเอง

การแปลง Shapefile เป็น KML ด้วย ogr2ogr

 shp2kml.bat
set GDAL_DATA=C:\Program Files\FWTools2.4.7\data
C:\"Program Files"\FWTools2.4.7\bin\ogr2ogr -f "KML" -s_srs "epsg:32647" -t_srs "epsg:4326" C:\nksw_index\index.kml C:\nksw_index\index.shp

limitation of shapefile

  เมื่อวานมีน้องคนหนึ่งโทรมาถาว่าทำไม่ dump ข้อมูลจาก database มาเป็น shapefile แล้วไม่สำเร็จสักที นั่งดูกันอยู่สักพักก็ถึงบางอ้อ เพราะว่า shapefile ก็มีข้อจำกัดที่ผมคิดว่าหลายๆคนยังไม่รู้ ข้อจำกัดที่ว่าได้แก่

         1. shapefile 1 ไฟล์ประกอบด้วยไฟล์ย่อยหลักๆคือ .shp, .shx, .dbf และ prj (srs info) โดยหนึ่ง shapefile จะบรรจุข้อมูลได้ไม่เกิน 2 GB
          2. shapefile 1 ไฟล์บรรจุข้อมูลได้ทั้งหมด 65535 records เกินไม่ได้
          3. ใน 1 recordเก็บข้อมูลได้มากสุด 1000 bytes หรือ ไม่เกิน 32 field
          4. ชื่อ Field limited แ ค่ 8 ตัวอักษร
          5. เก็บตัวอักษรใน field ได้มากสุด 254 byte
         6. ไม่ support Unicode

อ่านรายละเอียดเพิ่มเติมที่ http://www.esri.com/library/whitepapers/pdfs/shapefile.pdf

Create KML From Mysql SHAPE

update  province_geo set kml=astext(SHAPE)

SELECT REPLACE( FieldName ,'คำที่ค้นหา','คำที่เอามาแทน' ) AS NewData FROM TableName

update  amphoe_geo set kml=astext(SHAPE)

update  amphoe_geo set kml=replace(kml,')','');

update  amphoe_geo set kml=replace(kml,'(','');

update  amphoe_geo set kml=replace(kml,'MULTIPOLYGON','');

update  amphoe_geo set kml=replace(kml,'POLYGON','');

update  amphoe_geo set kml=replace(kml,',','|');

update  amphoe_geo set kml=replace(kml,' ',',');

update  amphoe_geo set kml=replace(kml,'|',' ');

Shape File to MySQL With ogr2ogr

ogr2ogr -f MySQL MySQL:my_database,user=root,password=mypassword TM_WORLD_BORDERS_SIMPL-0.1.shp -nln my_table -nlt MULTIPOLYGON -update -overwrite -lco GEOMETRY_NAME=my_polygon_column

ระบบสารสนเทศภูมิศาสตร์รายงานสถานการณ์โรคไม่ติดต่อและโรคเรื้อรัง (NCD GIS)

คู่มือการใช้งาน

ผลงานที่ผ่านมา

• สำนักงานสาธารณสุขจังหวัดร้อยเอ็ด

• สำนักงานสาธารณสุขจังหวัดตาก

รายละเอียดเพิ่มเติม

• Mobile : 086-645-1435
• Mail : ajarn.veerapong@gmail.com

ระบบสารสนเทศภูมิศาสตร์รายงานสถานการณ์โรคที่ต้องเฝ้าระวังทางระบาดวิทยา (R506 GIS)

คู่มือการใช้งาน

ผลงานที่ผ่านมา
- สำนักงานป้องกันควบคุมโรคที่ 10 เชียงใหม่
- สำนักงานป้องกันควบคุมโรคที่ 7 อุบลราชธานี
- สำนักงานป้องกันควบคุมโรคที่ 6 ขอนแก่น
- สำนักงานสาธารณสุขจังหวัดร้อยเอ็ด
- สำนักงานสาธารณสุขจังหวัดตาก

รายละเอียดเพิ่มเติม

• Mobile : 086-645-1435
• Mail : ajarn.veerapong@gmail.com

บริการรับทำแผนที่ เส้นชั้นความสูง (Contour) มาตราส่วน 1 : 4,000

เส้นชั้นความสูง (Contour Line) คือ เส้นที่แสดงลักษณะความสูงต่ำของพื้นที่ เป็นเส้นจินตนาการของระดับที่คงที่บนพื้นดิน เป็นเส้นที่ได้จากการลากเส้นคงที่ผ่านจุดต่าง ๆ บนพื้นดินที่มีค่าระดับเท่ากัน เส้นชั้นความสูงที่มีค่าเป็นบวก คือเส้นที่แสดงค่าความสูงเหนือระดับน้ำทะเลปานกลาง ส่วนเส้นชั้นความสูงที่มีค่าเป็นลบเป็นเส้นชั้นความสูงที่แสดงค่าความสูงใต้ ระดับน้ำทะเลปานกลาง ส่วนโครงการจัดทำแผนที่เพื่อการบริหารทรัพยากรธรรมชาติและทรัพย์สิน ของกระทรวงเกษตรและสหกรณ์ที่กลุ่มบริการแผนที่และภาพถ่ายออร์โธสี ให้บริการเส้นชั้นความสูงในปัจจุบัน มีรายละเอียดดังนี้
 เส้นชั้นความสูง มีคุณสมบัติดังนี้

1. จัดเก็บในรูปแบบของ Vector Files (.shp)
2. เส้นชั้นความสูงสร้างจากแบบจำลองระดับสูงแบบจุด
3. วัดจุดระดับความสูงที่เป็นตัวแทนภูมิประเทศทุกๆ ระยะ 4 เมตร

ตัวอย่างเส้นชั้นความสูง (Contour

รายละเอียดเพิ่มเติม
Mobile : 086-645-1435
Mail : ajarn.veerapong@gmail.com Website: www.ajarnveerapong.com

GIS กับ ระบบฐานข้อมูล

 ที่มา http://www.resgat.net/modules.php?name=News&file=article&sid=32
 
อันที่จริง ระบบ GIS เป็นฐานข้อมูลหรือ Database ชนิดหนึ่ง แต่เป็นฐานข้อมูลที่เกิดขึ้นมาชนิดที่เป็นฐานข้อมูลที่รวมข้อมูลประเภทข้อความและข้อมูลภาพหรือ graphic เข้าไปด้วย ซึ่งต่างจากระบบฐานข้อมูลทั่วไปเพราะในสมัยนั้นคำว่า ฐานข้อมูลหมายถึงข้อมูลที่เป็นตัวอักษรหรือตัวเลขแต่เพียงอย่างเดียว

ฉะนั้นในห้วงแรกเมื่อระบบ GIS เกิดขึ้นมาจึงต้องมีการเขียน Function เพื่อจัดการด้าน Graphic พร้อมๆกับการจัดการระบบระบบฐานข้อมูล ทำให้ผู้ที่ได้ใช้ระบบ GIS ในสมัยแรกต้องเรียนรู้ คำสั่งด้านการจัดการ GIS พร้อมๆกับคำสั่งด้านการจัดการฐานข้อมูลไปด้วย

ในปัจจุบันแนวความคิดที่จะนำระบบ GIS เข้าสู่ฐานข้อมูลอย่างเต็มตัวได้ปลุกกระแสให้ ผู้พัฒนาระบบฐานข้อมูลดังๆ อย่าง Oracle สนใจที่จะนำข้อมูล GIS เข้าสู่ในระบบฐานข้อมูลของตน จึงก่อให้เกิด ฐานข้อมูลประเภทใหม่ที่สามารถจัดเก็บข้อมูลประเภท Graphic ได้หรือที่เรียกว่า object-relational database ขึ้นมา

เมื่อข้อมูล GIS เข้าสู่ระบบฐานข้อมูลที่เป็น object-relational database สิ่งที่ได้ตามมาก็คือ ความปลอดภัย การจัดการอย่างมีระบบและการทำงานแบบหลายผู้ใช้ ตามแบบอย่างระบบฐานข้อมูลทั่วไป

แต่สิ่งที่สำคัญที่สุดของการนำเอาข้อมูล GIS เข้าสู่ระบบฐานข้อมูล object-relational database ก็คือเราสามารถจัดการกับข้อมูลโดยใช้คำสั่งที่เรียกว่า SQL หรือ Structure query language ซึ่งเป็นภาษาที่มีประสิทธิภาพมากและเป็นภาษาในยุคที่ 4 ได้

นั่นคือเราสามารถจะใช้คำสั่ง SQL ในการสืบค้น การวิเคราะห์ การจัดการข้อมูล ของข้อมูล GIS ซึ่งก่อให้เกิดประโยชน์มากมาย ในแง่ความเป็นมาตรฐานของการจัดการข้อมูล

แล้วเราจะลองใช้ ระบบฐานข้อมูล GIS ที่เป็นฐานข้อมูลจริงๆได้อย่างไรล่ะครับ นี่คงเป็นคำถามของผู้ที่อ่านมาถึงตรงนี้ ก็ขอตอบเลยว่า เราสามารถใช้ freeware ที่ชื่อว่า PostGIS ครับโดยที่ PostGIS จะทำงานร่วมกับระบบฐานข้อมูลที่เป็น object-relational database ของ PostgreSQL (ซึ่งก็เป็นระบบฐานข้อมูลประเภท OpenSource เช่นเดียวกับ MySQL ซึ่งเป็นระบบฐานข้อมูล OpenSourceยอดนิยมใน Internetครับ) โดยที่ PostGIS จะทำให้เราสามารถใช้คำสั่ง SQL จัดการกับข้อมูล GIS บน PostgreSQL ได้ครับ

การนำ text จาก cad แปลงไปเป็น shape ไฟล์

ที่มา http://thaitopo.editboard.com/t31-topic

ใช้โปรแกรม AutoCAD Map ที่ชุดคำสั่ง Map>Tools>Export

ที่แท็บ Selection เลือก Text  

ที่แท็บ Data เลือก Select Attributes

ขยาย + Properties แล้วติ๊กเลือกที่ STRING (ถ้าต้องการค่า Geometry อื่นๆ ก็ติ๊กเลือกได้ โดยโปรแกรมจะสร้าง Field ใน Database ให้อัติโนมัติ)

เสร็จแล้ว คลิก Ok...ตัวโปรแกรม จะสร้าง .shp ไฟล์ ให้โดยมีตัวหนังสือ อยู่ใน Field ที่ชื่อ TEXTSTRING

การสร้างชั้นข้อมูล ArcView จาก text file

การเปลี่ยนโซน และค่าพิกัดจาก Geographic สู่ระบบ UTM WGS84 Zone 47 ด้วย FWTOOS

การเปลี่ยนโซน และค่าพิกัดจาก Geographic สู่ระบบ UTM WGS84 Zone 47

set GDAL_DATA=C:\Program Files\FWTools2.4.7\data
del spot_utm_wgs84z478.*
C:\"Program Files"\FWTools2.4.7\bin\ogr2ogr -f "ESRI Shapefile" -s_srs epsg:4326 -t_srs epsg:32647 spot_utm_wgs84z48.shp spot.shp