Tag: coordinates

Update: โปรแกรมแปลงค่าพิกัดภูมิศาสตร์ Geographic Calculator (GeoCalc) บนเครื่องคิดเลข Casio fx-9860G II SD

Geographic Calculator

สืบเนื่องจากตอนก่อนหน้านี้ผมได้นำเสนอโปรแกรมแปลงพิกัด Geographic Calculator แบบไม่ได้ใช้ไลบรารีช่วยเรื่อง User Interface โปรแกรมมีลักษณะง่ายๆ เปิดมาเจอเมนูเลือกลักษณะที่จะคำนวณ จากนั้นโปรแกรมจะถามค่าพิกัดที่ต้องการแปลงแล้วคำนวณให้ ข้อดีคือใช้ง่าย ข้อเสียถ้าป้อนข้อมูลผิดพลาด จะย้อนกลับไม่ได้ ต้องเดินหน้าผิดไปจนจบ แล้วค่อยย้อนกลับมาอีกที

เปลี่ยนรูปแบบการติดต่อกับผู้ใช้ด้วยไลบรารี MyLib

ไลบรารี MyLib เป็นไลบรารีภาษาซีเล็กๆที่ผู้พัฒนาใช้นาม hayzel ได้เขียนไว้เพื่อใช้บนเครื่องคิดเลข Casio fx-9860G II SD ผมนำมาใช้และชอบ ทำให้มีแรงใจที่จะเขียนโปรแกรมบนเครื่องคิดเลขรุ่นเทพรุ่นนี้ได้มาหลายโปรแกรม และก็เหมือนเดิมครับว่าโปรแกรมที่ผมเขียนนั้นใช้งานได้ฟรี (Freely Usability) เพื่อใช้ในแวดวงงานสำรวจทำแผนที่ตลอดจนงานสำรวจเพื่อการก่อสร้างก็ตาม

ก็ขอตั้งชื่อโปรแกรมเล็กๆสำหรับแปลงพิกัดบนพื้นหลักฐาน WGS84 นี้ว่า “GeoCalc Extra” ก่อนอื่นสูตรที่ใช้ในการคำนวณผมใช้ไลบรารีชื่อ mgrs สามารถคำนวณแปลงพิกัดในระบบพิกัด UTM, geographic, UPS และ MGRS ได้ ขนาดไม่ใหญ่มากนัก สามารถคอมไพล์และบิวท์มาใส่เครื่องคิดเลขรุ่นนี้ได้ ทั้ง MyLib และ mgrs เป็นโปรแกรมเปิดโค้ด ฟรีทั้งคู่

ดาวน์โหลดโปรแกรม (Download)

ไปที่หน้าดาวน์โหลด (Download) มองหาโปรแกรมบนเครื่องคิดเลข Casio fx-9860G II SD  ชื่อโปรแกรม GeoCalc Extra จากนั้นทำการดาวน์โหลดมาจะได้ไฟล์ชื่อ “GEOCALC.G1A” แล้วทำการ copy ไฟล์ตัวนี้ไปยังเครื่องคิดเลขด้วยโปรแกรม Casio FA-124  หรือ copy ผ่านทางตัว SD Card ที่มากับเครื่องคิดเลข

เริ่มใช้โปรแกรม

ที่ “Main Menu” ของเครื่องคิดเลขเลื่อนไปหาไอคอนของโปรแกรมดังรูปด้านล่าง เมื่อกดคีย์ “EXE” แล้วจะเข้าเมนูของโปรแกรม ดังนี้

F1 – Set เลือกรายการคำนวณแปลงพิกัดระหว่าง UTM, Geographic หรือ MGRS

F2 – Src (Source) ป้อนค่าพิกัดที่ต้องการแปลงพิกัด

F3 – Calc คำนวณแปลงพิกัดพร้อมแสดงผลลัพธ์

F5 – Info แสดงเครดิตไลบรารีที่โปรแกรมใช้งาน

F6 – Exit ออกจากโปรแกรม

เลือกรายการคำนวณ (Menu)

ที่เมนูหลักกดคีย์ F1 – Set เพื่อเข้าไปเลือกรายการคำนวณ จะเห็น ระบบพิกัดเริ่มต้น (Source)   และระบบพิกัดปลายทาง (Target) ส่วนด้านล่าง MGRS Precision จะเป็นความละเอียดของระบบพิกัด MGRS (Military Grid Reference System) เลือกได้ 6 ระดับคือ 0, 2, 4, 6, 8, 10

ตัวอย่างที่ 1 แปลงค่าพิกัดจากค่าพิกัดภูมิศาสตร์ไปยังค่าพิกัดยูทีเอ็ม (Geographic to UTM)

ตั้งค่าระบบพิกัดเริ่มต้นและปลายทาง

ตั้งค่าระบบพิกัดต้นทางและปลายทางดังรูป จากนั้นกดคีย์ F6 – OK เพื่อออก

ป้อนค่าพิกัด

กลับมาที่เมนูหลักของโปรแกรมอีกครั้ง กดคีย์ F2 – Src เพื่อป้อนค่าพิกัดภูมิศาสตร์

ป้อนค่าพิกัดละติจูด Latitude 39°57’9.34803″N โดยการป้อน 39-57-9.34803N ค่าลองจิจูด 75°9’54.75490″W ป้อนค่า 75-9-54.75490W เสร็จแล้วกดคีย์ F6 – OK เพื่อออกไปคำนวณ

คำนวณแปลงพิกัด

กลับมาที่เมนูหลักของโปรแกรม กดคีย์ F3 – Calc เพื่อคำนวณจะได้ผลลัพธ์ดังรูปด้านล่าง โปรแกรม

จะแสดงค่าพิกัดเริ่มต้นให้และค่าพิกัดปลายทางคือยูทีเอ็มให้ พร้อมทั้งบอกโซนของยูทีเอ็มให้ กดคีย์ F6 – Done เพื่อออก

ตัวอย่างที่ 2 แปลงค่าพิกัดจากค่าพิกัดยูทีเอ็มไปยังค่าพิกัด MGRS (UTM to MGRS)

กำหนดค่าพิกัดยูทีเอ็ม (UTM) ดังนี้ N: 2642783.110, E: 232030.949 UTM Zone No: 46N กลับมาที่เมนูหลักของโปรแกรม กดคีย์ F1 – Set เพื่อเปลี่ยนรายการคำนวณ ตั้งค่าตามรูปด้านล่าง 

เสร็จแล้วกดคีย์ F6 – OK เพื่อออก

กลับมาเมนูหลักของโปรแกรมกดีย์ F2 – Src เพื่อออกป้อนค่าพิกัดยูทีเอ็มดังนี้ จากนั้นกดคีย์ F6 – OK เพื่อออก

กลับมาที่เมนูหลักของโปรแกรม กดคีย์ F3 – Calc เพื่อคำนวณ จะได้ผลลัพธ์ดังรูปด้านล่าง

ตัวอย่างที่ 3 แปลงค่าพิกัดจากค่าพิกัด MGRS ไปยังค่าพิกัดภูมิศาสตร์ (MGRS to Geographic)

กำหนดค่าพิกัด MGRS: 46QCK0907425049 ส่วนขั้นตอนจะขอรวบรัดแสดงด้วยรูปภาพ

จัดเก็บข้อมูลและเรียกมาใช้ภายหลัง

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

ก่อนจะใช้งานได้ต้องมีการเตรียมโฟลเดอร์บน SDCard ดังต่อไปนี้  คือดึง SDCard จากเครื่องคิดเลขมาเสียบบนคอมพิวเตอร์ แล้วทำการสร้างโฟลเดอร์ชื่อ “svdata” ดังรูป แต่ถ้ามีการสร้างมาแล้วก็ไม่จำเป็นต้องทำอะไร

จากนั้นนำ SDCard มาเสียบบนเครื่องคิดเลขอีกครั้ง เมื่อนำไปใช้งานได้สักพักถ้าเอามาเปิดอีกครั้งจะเห็นไฟล์หลายๆไฟล์ มีนามสกุลเป็น “CFG”  หมายถึง config ตัวอย่างถ้าใช้โปรแกรมคำนวณโค้งแปลงพิกัดภูมิศาสตร์นี้ไฟล์ที่จัดเก็บข้อมุลคือ “GEOCALC.CFG

ก็ตามที่สัญญาไว้ว่าจะไล่รื้อโปรแกรมเก่าๆ ที่ลงมาหลายๆตอนหน้านี้ด้วยระบบติดต่อผู้ใช้ตามไลบรารี mylib ที่ผมใช้อยู่ ติดตามกันตอนต่อไปครับ

ติดปีกเครื่องคิดเลขเทพ Casio fx 9860G II SD ด้วยโปรแกรมภาษาซีบน AddIn ตอนที่ 7 โปรแกรมคำนวณโค้งราบ (Horizontal Curve)

ติดปีกเครื่องคิดเลขเทพ Casio fx 9860G II SD ด้วยโปรแกรมภาษาซีบน AddIn ตอนที่ 7 โปรแกรมคำนวณโค้งราบ (Horizontal Curve)

โปรแกรมคำนวณโค้งราบ HCurve สำหรับเครื่องคิดเลข Casio fx-9860 G

ช่วงนี้อยู่ในชุดซีรี่ย์โปรแกรมเครื่องคิดเลข Casio fx-9860G ที่ใช้ไลบรารี MyLib ต่อไปขอนำเสนอโปรแกรมคำนวณโค้งราบ (Simple Horizontal Curve) ตัวโปรแกรมพัฒนาด้วยภาษาซี  ใช้เป็นโปรแกรม AddIn ดังรูปด้านล่าง

โปรแกรม HCurve คำนวณโค้งราบ (Simple Horizontal Curve)

องค์ประกอบของโค้งราบ (Elements of Horizontal Curve)

องค์ประกอบของโค้งราบ (Elements of Curve)

คำนิยาม (Abbrivations)

R – Radius คือรัศมีของโค้งราบ รัศมีจะตั้งฉากกับเส้นสัมผัสวงกลมเสมอ

PC –  Point of Curvature คือจุดเริ่มต้นโค้ง บางครั้งเรียกว่า BC (beginning of curve) หรือ TC (tangent to curve)

PI – Point of Intersection คือจุดตัดของ tangent 2 เส้น

PT – Point of Tangency คือจุดสิ้นสุดโค้ง บางครั้งเรียกว่า EC (end of curve) หรือ CT (curve to tangent)

POC – Point of Curve คือจุดบนโค้งในตำแหน่งใดก็ตาม

L – Length of Curve คือความยาวโค้งวัดตามโค้งจากจุด PC ไปจุด PT

T – Tangent Distance หรือ Tangent Length คือเส้นตรงที่สัมผัสโค้งวัดจากจุด PC ไปจุด PI หรือวัดจาก จุด PI ไปจุด PT

คุณสมบัติของโปรแกรม

โค้งราบ (Simple Horizontal Curve) ในงานสำรวจใช้ในงานสำรวจเพื่อการก่อสร้าง (Construction Survey) โดยที่ให้ตำแหน่ง (Setting out) งานก่อสร้างถนน ทางรถไฟ ผู้ใข้ป้อนข้อมูลองค์ประกอบของโค้งราบที่โปรแกรมต้องการจนครบ โปรแกรมสามารถคำนวณองค์ประกอบโค้งที่เหลือและสามารถคำนวณค่าพิกัดบนเส้น Center line หรือกระทั่งสามารถ offset ไปด้านซ้ายหรือด้านขวาก็ได้เช่นเดียวกัน สามารถคำนวณค่าพิกัดแบบทุกช่วงระยะ (interval) ให้ค่าพิกัดออกมาเป็นบัญชีรายการได้ โปรแกรมนี้ออกแบบและพัฒนามาเพื่อช่วยช่างสำรวจให้สามารถนำโปรแกรมไปตรวจสอบข้อมูลโค้งราบได้ด้วยตัวเอง

ในปัจจุบันการวางโค้งในสนาม ไม่ได้ถูกจำกัดเช่นสมัยแต่ก่อนเนื่องจากเครื่องมือเช่นอุปกรณ์กล้องสำรวจ Total Station ทันสมัยสามารถให้ตำแหน่งจากค่าพิกัดได้เลย ไม่เหมือนสมัยแต่ก่อนที่มีแต่กล้อง Theodolite และเทปวัดระยะ ที่การวางโค้งต้องอาศัยการตั้งกล้องที่จุด PC หรือ PI หรือแม้กระทั่ง PT แล้วทำการเปิดมุมและดึงเทปไปตามคอร์ด ปัจจุบันไม่ต้องทำอย่างนั้นแล้ว สามารถตั้งกล้องที่ไหนก็ได้ที่สะดวก

ก็เนื่องจากความทันสมัยของกล้อง Total Station นี้เอง ในอนาคตไม่ไกลการโอนข้อมูลจากกล้องไปที่ออฟฟิศหรือจากออฟฟิศไปที่กล้องจะเป็นเรื่องธรรมดามากผ่านทางออนไลน์ ซึ่งเรื่องนี้ข้อดีก็มีมากมาย แต่ข้อเสียที่เกิดขึ้นคือช่างสำรวจจะมีเวลาใช้สมองคิดเรื่อง geometry น้อยลงเพราะโปรแกรมบนกล้องทำให้หมด ทำให้ขาดการฝึกฝนทักษะในด้านนี้ ซึ่งผมคิดว่าอาจจะทำให้คุณภาพของบุคลากรด้านสำรวจของเราด้อยลงในอนาคต ดังนั้นการใช้เครื่องคิดเลขมาช่วยอาจจะช่วยฝึกฝนทักษะได้บ้างในจุดนี้

ดาวน์โหลดโปรแกรม (Download)

ไปที่หน้าดาวน์โหลด มองหาโปรแกรมคำนวณโค้งราบ HCurve  จะได้ไฟล์ชื่อ “HCURVEEX.G1A” เมื่อดาวน์โหลดมาแล้วโอนเข้าเครื่องคิดเลขผ่านทางโปรแกรม FA-124 หรือ SD Card

ใช้งานฟรี (Freely Usable)

ก็ยังเหมือนเดิมไม่ว่าจะโปรแกรมบนเครื่องคอมพิวเตอร์หรือบนเครื่องคิดเลข คือให้ใช้งานได้ฟรี

เริ่มต้นใช้งาน

ที่เครื่องคิดเลขกดคีย์ “Menu” เพื่อเข้าสู่หน้า AddIn หรือ  Main Menu ของเครื่องคิดเลข เลื่อนลูกศรลงไปด้านล่างๆจะเห็นไอคอนของโปรแกรมดังรูปด้านล่าง กดคีย์ “EXE” เพื่อเข้าสู่โปรแกรม

จะเห็นหน้าเมนูหลักของโปรแกรมดังนี้

ก็เหมือนโปรแกรมคำนวณโค้งดิ่ง (VCurve) ที่ผ่านมาคือมีเมนูย่อยเรียงจากซ้ายไปขวา สัมพันธ์กับการกดคีย์ F1-F6 เรียงกันไป

เมนูหลัก (Main Menu)

F1 (Know) – Known สำหรับป้อนจุดที่ทราบค่า station และค่าพิกัด

F2 (Angl) – Angle สำหรับป้อนที่เกี่ยวกับมุมเช่น back tangent azimuth และมุมเบี่ยงเบน (Deflection Angle) ตลอดจนถึงทิศทางของโค้ง (Curve Direction)

F3 (Elem) – Elements สำหรับป้อนองค์ประกอบของโค้งเช่นรัศมีหรือความยาวโค้าง

F4 (Info) – Information สำหรับคำนวณหาข้อมูลพื้นฐานของโค้งทั้งหมด

F5 (Calc) – Calculate สำหรับคำนวณหาค่าพิกัดโค้งได้หลายรูปแบบเช่นกำหนดสถานี ระยะ offset ตลอดจนคำนวณจากช่วงระยะทาง (interval) ที่กำหนดให้

F6 (Exit) – Exit ออกจากโปรแกรม

ตัวอย่างการคำนวณโค้งราบ (Example)

มาดูวิธีการใช้งานจากตัวอย่างจะเข้าใจง่ายที่สุด

เลือกและป้อนสถานีที่ทราบและค่าพิกัด (Known Station and Know Coordinates)

ที่เมนูหลักกดคีย์ F1 (Know) แก้ไขค่าตามโจทย์ตัวอย่างที่ 1 ดังนี้

ในที่นี้สถานีและค่าพิกัดกำหนดที่จุด PI ทั้งคู่ เมื่อป้อนค่าเสร็จแล้วกดคีย์ F6 (OK) เพื่อออก

เลือกป้อนมุมอะซิมัทและมุมเบี่ยงเบน(Known Tangent and Deflection Angle)

ที่เมนูหลักกดคีย์ F2 (Angl) เพื่อเลือกอะซิมัทของเส้นสัมผัสที่ทราบค่าและมุมเบี่ยงเบน ตลอดจนป้อนทิศทางของโค้งราบว่าเลี้ยวซ้ายหรือเลี้ยวขวา จากตัวอย่างเลือกและป้อนค่าดังนี้

เสร็จแล้วกดคีย์ F6 (OK) เพื่อจัดเก็บค่าและออก

ป้อนองค์ประกอบโค้งราบ (Elements of Curve)

ที่เมนูหลักกดคีย์ F3 (Elem)  เลือกว่าจะป้อนค่ารัศมีโค้ง (Radius) หรือว่าความยาวโค้ง (Length of Curve) ในที่นี้เลือกรัศมีโค้งและป้อนค่า 201.950 เมตร

คำนวณหาข้อมูลพื้นฐานโค้งราบ (Curve Information)

ที่เมนูหลักกดคีย์ F4 (Info)

แสดงองค์ประกอบของโค้งและข้อมูลพื้นฐาน ค่าพิกัดของจุด PC, PI และ PT ตลอดจนจุดศูนย์กลางของโค้ง เนื่องจากจอภาพมีขนาดเล็กดังนั้นใช้การกดคีย์ F1 (PgUp) หรือ F2 (PaDn) เพื่อเลื่อนดูหน้าก่อนหน้านี้หรือหน้าถัดไป กดคีย์ F6 (Done) เพื่อออก

การคำนวณโค้งราบ (Horizontal Curve Calculation)

ที่เมนูหลักกดคีย์ F5 (Calc)จะเห็นเมนูย่อยอีกเมนูคือเมนูสำหรับคำนวณโค้งราบ

จะมีเมนูดังนี้

F1 (Sta) – Station คำนวณหาค่าพิกัดเมื่อกำหนดสถานี

F2 (INT) – Interval คำนวณหาค่าพิกัดสถานีเมื่อกำหนดช่วงระยะทาง (Interval) ให้

F4 (Info) – Information คำนวณข้อมูลพื้นฐานโค้งราบ โดยที่ผลลัพธ์เหมือนกับเมนู Info บนเมนูหลัก

F5 (Plot) – Plot Curve วาดรูปร่างโค้งราบ

F6 (Done) ออกจากเมนูคำนวณโค้งราบ

คำนวณหาค่าพิกัดเมื่อกำหนดสถานี (Calculate Coordinates of Station)

ที่เมนูคำนวณโค้งราบกดคีย์ F1 (Sta) จะมีไดอะล็อกให้ป้อนสถานี ตัวอย่างนี้ต้องการทราบค่าพิกัดของสถานี 17+200 โดยที่ offset ไปด้านซ้าย 8 เมตร ป้อนข้อมูลดังรูป ถ้าไม่ต้องการคำนวณหรือเก็บข้อมูลที่ป้อนก็กดคีย์ F5 (Canc) เพื่อ Cancel ออกไป หรือต้องการเก็บค่าแต่ไม่คำนวณก็กดคีย์ F6 (OK) ออกไป ถ้าต้องการคำนวณก็กดคีย์ F1 (Calc) จะได้ผลลัพธ์ดังรูปถัดไป กดคีย์ “EXE” เพื่อออก

คำนวณหาค่าพิกัดสถานีแบบกำหนดช่วงระยะทาง (Interval Calculation)

ที่เมนูคำนวณโค้งกดคีย์ F2 (INT) ในที่นี้ต้องการคำนวณทุกๆระยะ 25 เมตร โดยคำนวณในแนว Center Line (ไม่มีการ offset ไปซ้ายหรือขวา) กดคีย์ F1 (Calc) เพื่อคำนวณ

จะได้ผลลัพธ์ เริ่มตั้งแต่ PC (17+151.314), Sta 17+175, Sta 17+200, Sta 17+225, Sta 17+250, Sta 7+275, Sta 17+300 และสุดท้ายที่ PT (17+313.794) กด F1 (PgUp) เพื่อเลื่อนไปหน้าก่อนหน้านี้ และกด F2 (PgDn) เพื่อไปดูหน้าถัดไป กด F6 (Done)

 

วาดรูปโค้งราบ (Plot Curve)

จากเมนูคำนวณโค้งราบ กดคีย์ F5 (Plot) จะมีเมนูย่อยลงไปอีกสำหรับย่อ F1 (Z-) ขยาย F2 (Z+) ดึงรูปไปด้านซ้าย F3 (Lt) ดึงรูปไปด้านขวา F4 (Rt) กดคีย์ F5 (>) เพื่อไปเมนูย่อยอีกเมนูด้านขวา เมนูด้านขวาจะมีดึงรูปลง F1 (Dn) หรือดึงรูปขึ้น F2 (Up) กดคีย์ F6 (Exit) เพื่อออกจากเมนู

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

จัดเก็บข้อมูลและเรียกมาใช้ภายหลัง

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

ก่อนจะใช้งานได้ต้องมีการเตรียมโฟลเดอร์บน SDCard ดังต่อไปนี้  คือดึง SDCard จากเครื่องคิดเลขมาเสียบบนคอมพิวเตอร์ แล้วทำการสร้างโฟลเดอร์ชื่อ “svdata” ดังรูป แต่ถ้ามีการสร้างมาแล้วก็ไม่จำเป็นต้องทำอะไร

จากนั้นนำ SDCard มาเสียบบนเครื่องคิดเลขอีกครั้ง เมื่อนำไปใช้งานได้สักพักถ้าเอามาเปิดอีกครั้งจะเห็นไฟล์หลายๆไฟล์ มีนามสกุลเป็น “CFG”  หมายถึง config ตัวอย่างถ้าใช้โปรแกรมคำนวณโค้งราบนี้ไฟล์ที่จัดเก็บข้อมุลคือ “HCVEX.CFG

Surveyor Pocket Tools – ทดสอบโปรแกรมการแปลงพิกัดบน State Plane Coordinate System (SPC)

Surveyor Pocket Tools – ทดสอบโปรแกรมการแปลงพิกัดบน State Plane Coordinate System (SPC)

ในฐานะช่างสำรวจในย่าน AEC บ้านเรา ส่วนใหญ่จะคุ้นเคยกับระบบพิกัดที่ส่วนใหญ่ใช้เส้นโครงแผนที่ Transverse Mercator กันส่วนใหญ่ แต่มาเลเซียนั้นต่างออกไปเนื่องจากมีพื้นที่ที่ยาวเฉียงๆ ทั้งสองพื้นที่คั่นด้วยทะเลจีนใต้ พื้นที่แรกอยู่บนเกาะบอร์เนียวอีกพื้นที่หนึ่งติดกับประเทศไทย ทางมาเลเซียใช้เส้นโครงแผนที่ Oblique Mercator ซึ่งเป็นเส้นโครงแผนที่ค่อนข้างซับซ้อนกว่าอันอื่น

เรามาลองไปทัศนศึกษาที่หรัฐอเมริกาดูกัน สหรัฐอเมริกาเป็นประเทศใหญ่ มีระบบพิกัดและเส้นโครงแผนที่ที่หลากหลายมาก ในฐานะช่างสำรวจพอจะเป็นความรู้ประดับบ่ากันไว้นิดๆหน่อยๆ ไม่ถือว่าเหลือบ่ากว่าแรงจนต้องแบกหาม มาทัศนาผ่านทางโปรแกรมแปลงพิกัด Transform Coordinate ที่อยู่ในชุด Surveyor Pocket Tools

State Plane Coordinate System (SPC)

  • คือกลุ่มระบบพิกัดของสหรัฐอเมริกาที่รวบรวมระบบพิกัดที่ใช้ในแต่ละพื้นที่หรือรัฐทั้งหมด 124 โซน โดยที่รัฐที่มีพื้นที่ติดกันจำนวน 110 โซน และอื่นๆที่เหลือได้แกเช่นอลาสก้า ฮาวาย เปอร์โตริโก้และหมู่เกาะยูเอสเวอร์จิน
  • State Plane Coordinate System (SPC) ปัจจุบันคือ North America Datum 1983 (NAD83) ใช้ทรงรี GRS80 ต่างจากของเดิมคือ NAD27 ที่ใช้ทรงรี Clark 1866

AcroRd32_2017-03-22_08-20-10.jpg

  • เส้นโครงแผนที่ (Map projection) ใช้อยู่ 3 ประเภทคือ
    1. Transverse Mercator (TM) สามารถรักษา scale factor คงที่ได้ในแนวแกนเหนือ-ใต้ จึงนิยมใช้สำหรับพื้นที่ที่ยาวจากเหนือไปใต้
    2. Lambert Conformal Conic (LCC) เนื่องจากสามารถรักษา scale factor ให้คงที่ได้ในแนวแกนตะวันออก-ตะวันตก จึงนิยมใช้กับพื้นที่มีความยาวจากตะวันออกไปตะวันตก
    3. Oblique Mercator (OM) นิยมใช้กับพื้นที่ที่ยาวเฉียงแบบทะแยง ใช้อยู่รัฐเดียวคืออลาสก้า

โปรแกรมแปลงพิกัด Transform Coordinate

  • ถ้าจะดาวน์โหลดก็ดูที่ช่องยาวๆด้านขวาของ blog ในขณะที่เขียนบทความนี้ เป็นรุ่น 0.66 build 501 ดาวน์โหลดแล้วก็ติดตั้ง เมื่อคลิกโปรแกรมและรันจะเห็นไดอะล๊อกหน้าแรกรวมโปรแกรมชุดของ Surveyor Pocket Tools ตามรูปด้านล่างที่ไฮไลต์ไว้คือ “Transform Coordinates” 

Surveyor Pocket Tools_2017-03-21_20-45-41

  • มาดูหน้าโปรแกรม ผม update ข้างในไปมาก ทั้งๆที่ตอนแรกตั้งใจจะเขียนให้ใช้งานได้เฉพาะเส้นโครงแผนที่ Transverse Mercator ไปๆมาๆ ก็ไปไกลเกินกว่าที่ตั้งใจไว้ตอนแรก หน้าตาก็เรียบง่ายแต่เพิ่มเครื่องมือจัดเก็บค่าพิกัดและเรียกใช้ค่าพิกัด รวมถึงเครื่องมือปักหมุดบน google maps & google earth

ทดสอบการแปลงพิกัดบนเส้นโครงแผนที่ Oblique Mercator (OM)

  • จะทดสอบการแปลงพิกัดจากค่าแลตติจูดและลองจิจูด จากทรงรี GRS80 ไปยังระบบพิกัดฉาก NAD83 สำหรับ Alaska zone 1 นี้ใช้ Oblique Mercator ซึ่งพารามิเตอร์ของเส้นโครงแผนที่ OM ตามตารางด้านล่าง
    •  Alazka zone 1 – NAD83
    • Scale factor at the projection’s center = 0.9999
    • Longitude of the projection’s center = 133º 40′ W
    • Latitude of the projection’s center = 57º 0′ N
    • Azimuth at the projection’s center = 323º07’48.3685″
    • Angle from Rectified to Skew Grid =  323º07’48.3685″
    • False Easting (meters) = 5000000
    • False Northing (meters) = -5000000

AcroRd32_2017-03-22_15-44-02

ข้อมูลทดสอบโปรแกรม

  • ข้อมูลทดสอบผมจะใช้หมุดของทางการ ที่สะดวกมากสามารถจะดูหมุดที่ไหนก็ได้ ข้อมูลแต่ละหมุดจะเรียกว่า data sheet เข้าไปดูได้ที่เว็บไซต์ของ NGS ตาม ลิ๊งค์ นี้ครับ
  • อันดับแรกเลือกรัฐก่อน เลือก “Alaska” ต่อจากนั้นเลือก county ให้เลือก “AK|103|KETCHIKAN GATEWAY BOROUGH” จากนั้นคลิกที่ปุ่ม “submit” แล้วจะเห็นหมุดถูกลิสต์ออกมามากพอสมควร ลองเรียงหมุดดูตาม latitude ผมคลิกเลือก “Latitude” แล้วคลิกที่ปุ่ม “Re-Sort-By” ตอนนี้จะเลือกหมุดที่มีหมายเลข PID = UV5754 เลือกแล้วคลิกที่ปุ่ม “Get data sheets” จะได้ data sheet ออกมา
  • ลองดู data sheet  ของหมุด  UV5754 ด้านล่าง

chrome_2017-03-22_08-52-04

  • ที่ผมลากสี่เหลี่ยมไว้ด้านบนคือค่าพิกัดในระบบ geographic ของ NAD83 จะใช้ค่านี้มาทดสอบ เลือกระบบพิกัดและป้อนข้อมูลเข้าไปดังรูป Latitude = 55° 6′ 45.20173″N Longitude = 131° 43′ 58.97516″W

Surveyor Pocket Tools_2017-03-22_09-02-26

  • ด้านขวาปลายทางเลือก Group = “Projected Coordinate System” เลือก Datum = “NAD83 (National Spatial Reference System 2007)” จากนั้นเลือก System = “NAD83(NSRS2007) / Alaska zone 1” พร้อมแล้วคลิกลูกศรขวาเพื่อทำคำนวณ

Surveyor Pocket Tools_2017-03-22_09-07-59

  • จะได้ค่า North = 366701.8435, East = 942069.8596, Grid scale factor = 0.9999085667,  Convergence = 1°36’31.76352″ ซึ่งผลลัพธ์ตรงกันกับ data sheet

ทดสอบการแปลงพิกัดบน Lambert Conformal Conic (LCC)

  • สำหรับเส้นโครงแผนที่ Lambert Conformal Conic (LCC) ส่วนใหญ่จะเป็นแบบรอยตัดสองรอย (secant) บนทรงรี มากกว่าจะเป็นแบบสัมผัส ดังนั้นจะมี Latitude of parallel อยู่ตรงสองรอยตัดด้านบนและด้านล่าง สำหรับข้อมูลทดสอบจะเลือกรัฐ “Oregon” โซนด้านเหนือ พารามิเตอร์สำหรับการแปลงพิกัดมีดังนี้
    • Oregon North Zone (Designation 3601)
      • Oregon State Plane North – NAD 1983
      • Lambert Conformal Conic Two Standard Parallel Projection (Secant)
      • Central Meridian: -120° 30′ W
      • Latitude of Origin: 43° 40′ N
      • Standard Parallel (South): 44° 20′ N
      • Standard Parallel (North): 46° N
      • False Northing: 0.000 m
      • False Easting: 2 500 000.000 m

ข้อมูลทดสอบ

  • ข้อมูลทดสอบจะเข้าไปที่เว็บไซต์ของ NGS แต่ครั้งนี้จะเข้าไปใช้ interactive map แล้วค้นหาชื่อหมุด ตาม ลิ๊งค์ นี้ ด้านซ้ายมือจะมีระบบค้นหาเลือกค้นหาด้วย “PID” ป้อนชื่อหมุด “AJ8179” เมื่อเจอแล้วคลิกที่หมุดในแผนที่แล้ว แล้วคลิก “data sheet” จะได้รายละเอียดมาดังรูป

chrome_2017-03-22_17-16-10

  • ป้อนค่าพิกัดฉากของหมุด AJ8179 เข้าไปดังรูป คลิกที่รูปลูกศรขวา เพื่อจะแปลงพิกัดจาก LCC ไป geographic แต่ในขณะเดียวกันโปรแกรมจะคำนวณค่า grid scale factor และ convergence มาให้ด้วย

Surveyor Pocket Tools_2017-03-22_17-15-44

  • ได้ผลลัพธ์ดังนี้ ตรวจดูกับ data sheet จะได้ค่า grid scale factor และ convergence ตรงกัน ส่วนค่าพิกัดค่าลองจิจูดตรงกัน แต่ค่าแลตติจูดต่างกันเล็กน้อยมากที่ทศนิยมที่ 5

Surveyor Pocket Tools_2017-03-22_17-29-22

ข้อจำกัดของโปรแกรม

  • ในตอนนี้ยังไม่สนับสนุนหน่วยฟุต
  • สนับสนุนการคำนวณ grid scale factor และ convergence ให้เฉพาะเส้นโครงแผนที่ Transverse Mercator, Lambert Conformal และ Oblique Mercator ส่วนเส้นโครงแผนที่อื่นๆคำนวณให้เฉพาะการแปลงค่าพิกัดเท่านั้น เนื่องจากระบบพิกัดในโลกนี้หลากหลายมากมาย ทำให้การทดสอบข้อมูลต้องทยอยทำไปเรื่อยๆ
  • ยังไม่ได้ทดสอบข้อมูลระบบพิกัดของยุโรป European Terrestrial Reference System (ETRS)
  • มาถึงตอนนี้คงพอหอมปากหอมคอ สังเกตว่าผมไม่ได้ทดสอบเส้นโครงแผนที่ TM เนื่องจากบ้านเราใช้กันอยู่และคุ้นเคยกันดีอยู่แล้ว พบกันตอนต่อไปครับ

Surveyor Pocket Tools – Update เพิ่มพื้นหลักฐาน Myanmar Datum 2000 (MMD 2000)

  • Surveyor Pocket Tools ก่อนหน้านี้ได้เพิ่มพื้นหลักฐาน (Datum) ของเพื่อนบ้าน AEC ของเราไปแล้วคือ สปป.ลาว (Lao National Datum 1997) กับเวียดนาม (VN-2000) มาตอนนี้มาอัพเดทของเมียนมา (Myanma Dtum 2000) ประเทศเพื่อนบ้านของเรา
  • รู้สึกว่าปี 2000 ที่ผ่านมาเข้าสหัสวรรษใหม่ จะมีระบบพิกัดของบางประเทศที่ประกาศใช้ในปีนั้น เลขสวยและขลัง 🙂 ถึงแม้ในบ้านเราจะพยายามปรับมาใช้พื้นหลักฐาน WGS84 กันมาก แต่เพื่อนบ้านเราก็ยังใช้ระบบพิกัดของประเทศตัวเองอยู่พอสมควร เนื่องจากแผนที่พื้นฐานหรือฐานข้อมูลยังอ้างอิงกับระบบนี้อยู่
  • ผมไปค้นหาระบบพิกัด Myanmar Datum 2000 ในเสิร์ชเอนจิ้น แต่ก็ต้องผิดหวังเพราะไม่มีในสารบบ ผมคิดว่าช่วงเมียนมาร์ ประกาศใช้ช่วงนั้นยังเป็นรัฐบาลทหารหรือกระไร ถึงไม่มีการตอบรับ ดังนั้นหมายเลขพื้นหลักฐาน หมายเลขเส้นโครงแผนที่ UTM จะไม่มีในฐานข้อมูลของ EPSG แม้แต่น้อย
  • Myanmar Datum 2000 จะพารามิเตอร์ในการแปลงพิกัดเพียงสามพารามิเตอร์คือ Translation DX,DY,DZ เท่านั้น

ข้อกำหนดของ Myanmar Datum 2000

  • ใช้ทรงรี Everest 1830 
    • Semi major axis (a) = 6377276.345 m
    • Semi minor axis(b) = 6356075.434 m
    • Flattening = 1 / 300.8017
  • ค่าพารามิเตอร์สำหรับการแปลงพิกัด (Transformation Parameters) จาก MMD2000 ไปยัง WGS84 :
    • DX = 246.632; DY = 784.833; DZ = 276.923
  • ใช้เส้นโครงแผนที่ UTM ซึ่งจะมีสองโซนคือ
    • UTM zone 46N มี CM = 93°; False Easting = 500000; False Northing = 0; K0 = 0.9996
    • UTM zone 47N มี CM = 99°; False Easting = 500000; False Northing = 0; K0 = 0.9996

เพิ่มพื้นหลักฐาน MMD2000 บน Surveyor Pocket Tools

  • สามารถใช้ได้บนโปรแกรมที่ Build 451 ขึ้นไป รูปด้านล่างเป็น “Transform Coordinate” แปลงพิกัดทีละจุด
Surveyor Pocket Tools support “Myanmar Datum 2000”
  • ถ้าต้องการแปลงค่าพิกัดแบบอ่านจากไฟล์ก็ใช้ “File Transform Coordinates” ก็รองรับและใช้ได้
Surveyor Pocket Tools support “Myanmar Datum 2000”
  • Proj.4 string  ที่ใช้สำหรับการแปลงพิกัดที่ส่งไปให้ไลบรารี pyproj เป็นดังนี้

+proj=utm +zone=47 +a=6377276.345 +b=6356075.41314023 +towgs84=246.632,784.833,276.923,0,0,0,0 +units=m +no_defs

  • พบกันตอนใหม่ครับผู้อ่านทุกท่าน

Surveyor Pocket Tools – Update เพิ่มพื้นหลักฐาน VN-2000 (Vietnam)

Vietnam 2000 Datum (VN-2000)

  • ผมมีโอกาสได้จับงานเอกสารงานในกรุงโฮจิมินท์ของเวียดนาม (แต่ยังไม่ได้ไปทำหรอกครับ) ซึ่งในงานนั้นมีงานสำรวจอ้างอิงถึงระบบพิกัดของเวียดนามคือ VN-2000 ที่งานนี้จะต้องใช้อ้างอิงระบบพิกัดนี้ ก็เลยไปค้นหาในอินเตอร์เน็ตก็ได้พบใบประกาศใช้ระบบพิกัดนี้ ของทางราชการของเวียดนามคือ Ministry of Natural Resources and Environment (MONRE) ผมดาวน์โหลดมาแต่กลับไปค้นหาลิ๊งค์ต้นฉบับไม่เจอ ถ้าต้องการดูก็ตามลิ๊งค์นี้ครับ
  • ผมเห็นบทความของของอาจารย์ในมหาวิทยาลัยของเวียดนามท่านนี้ได้พูดถึงระบบพิกัด VN-2000 ของเวียดนามกับ EGM2008 ว่าระบบพิกัดของเวียดนามที่ประกาศใช้ ตามที่ผมได้กล่าวไปแล้วข้างต้นจาก MONRE  นอกจากจะมีการเลื่อน (Translation) แกน DX,DY,DZ แล้วยังมีการหมุน (Rotation) รอบแกน RX,RY,RZ และยังมี scale ด้วย ใช้พารามิเตอร์ครบทั้งเจ็ดตัวในการ Transformation การแปลงพิกัดแบบเจ็ดพารามิเตอร์เพื่อดึงค่าความสูงจีออยด์ (Geoid separation) จาก EGM2008 มาใช้ก็ดูจะยุ่งยาก ผมจับใจความได้ว่า แค่แปลงพิกัด VN-2000 ไป WGS-84 ใช้พารามิเตอร์แค่การเลื่อน DX,DY,DZ สามพารามิเตอร์ก็น่าจะเพียงพอ ผมก็เห็นด้วยเพราะ EGM2008 แบบละเอียดสุดคือ ขนาดกริด 1 min x 1 min ประมาณ 1.8 กม. x 1.8 กม. การแปลงพิกัดใช้สามพารามิเตอร์ได้ค่าพิกัดที่คลาดเคลื่อนบวกลบ สัก 10-20 เมตร คงไม่ใช่ปัญหา
  • แต่ถ้าสำหรับงานรังวัด GNSS เนื่องจากเป็นงาน Precise Surveying ค่าพารามิเตอร์ต้องครบเจ็ดตัว ถึงจะได้ค่าพิกัดที่ละเอียดในระดับมิลลิเมตร

ข้อกำหนดของ VN-2000 Datum

  • ใช้ทรงรี (Ellipsoid) : WGS84
  • หมายเลข EPSG พื้นหลักฐาน VN-2000 : 6756
  • ค่าพารามิเตอร์สำหรับการแปลงพิกัด (Transformation Parameters) จาก VN-2000 ไปยัง WGS84 :
    • Translation Dx = –191.90441429; DY = –39.30318279; DZ = –111.45032835
    • Rotation RX = –0.00928836; RY = 0.01975479; RZ = –0.00427372 (หน่วยมุมสำหรับการหมุนนี้เป็นหน่วยฟิลิปดา (second)
    • สเกล (scale) ต้องระวังนิดหนึ่งเพราะหน่วยเป็น ppm (หนึ่งส่วนต่อหนึ่งล้าน) ถ้าเขียนไว้ว่า k = 1.000000252906278 แสดงว่าเขียนไว้เต็มๆ ไม่ใช่ ppm เขียนใหม่ได้ = 1 + 0.252906278 / 10 สังเกตชุดหลังหารด้วยด้วยหนึ่งล้าน ตัวเลขที่เอาไปป้อนในโปรแกรมแปลงพิกัดคือตัวเลข Δs = 0.252906278 (ppm)
  • ใช้เส้นโครงแผนที่ UTM ซึ่งจะมีสองโซนคือ
    • UTM zone 48N มี CM = 105°; False Easting = 500000; False Northing = 0; K0 = 0.9996 หมายเลขของ EPSG = 3405
    • UTM zone 49N มี CM = 111°; False Easting = 500000; False Northing = 0; K0 = 0.9996 หมายเลขของ EPSG = 3406

เพิ่มพื้นหลักฐาน VN-2000 บน Surveyor Pocket Tools

  • ที่เพิ่มพื้นหลักฐานของ VN-2000 ใช้ได้บนโปรแกรมที่ Build 451 เป็นต้นไป ถ้าผู้อ่านดาวน์โหลดโปรแกรมรุ่นก่อนหน้านี้ก็มาดาวน์โหลดกันใหม่

  • Proj.4 string ที่ผมใช้ในการแปลงพิกัดสำหรับโซน 48N เขียนได้ดังนี้

+proj=utm +zone=48 +ellps=WGS84 +towgs84=-191.90441429,-39.30318279,-111.45032835,-0.00928836,0.01975479,-0.00427372,0.252906278 +units=m +no_defs

  • ตอนนี้โปรแกรมในรุ่นนี้ผมเขียนสถานะการใช้ระบบพิกัดไว้ในไฟล์ xml เล็กๆ เพื่อที่เวลาปิดโปรแกรมไป เวลาเปิดมาใช้งานใหม่จะได้เอาระบบพิกัดที่เคยใช้มาแสดง ไม่ต้องเลือกใหม่ให้เสียเวลา
  • ในส่วน “File Transform Coordinate” ก็อัพเดทเหมือนกัน

  • พบกันตอนหน้าครับ

Surveyor Pocket Tools – เปิดตัวโปรแกรมแปลงไฟล์พิกัดข้ามพื้นหลักฐาน File Transform Coordinates

File Transform Coordinates

  • โปรแกรมแปลงพิกัดข้ามพื้นหลักฐาน Transform Coordinates ที่ผมนำเสนอมาก่อนนั้น แปลงพิกัดได้ทีละจุด อาจจะไม่สะดวกถ้าผู้อ่านมีจุดตั้งแต่ 5-10 จุดขึ้นไป ทางออกผมเลยเขียนโปรแกรมเพิ่มอีกตัวเข้ามา โดยอ่านไฟล์พิกัดที่ต้องการแปลง โดยที่ไฟล์นั้นจะเก็บไว้ในรูปแบบ CSV ที่มีตัวคั่นเป็นเครื่องหมายคอมม่า
  • ขอตั้งชื่อโปรแกรมเป็น File Transform Coordinates โดยใส่คำว่า File นำหน้าชื่อโปรแกรมตัวเดิม เพื่อให้สื่อความหมาย ว่าแปลงพิกัดจากไฟล์

ดาวน์โหลดและติดตั้ง Surveyor Pocket Tools

  • ตอนนี้รุ่นโปรแกรม Surveyor Pocket Tools ที่มีโปรแกรมย่อย File Transform Coordinates มาด้วยนั้น ต้อง build 447 ขึ้นไป ให้ดูด้านขวาที่ช่องดาวน์โหลด (Download) คลิกและก็ดาวน์โหลดมา unzip และทำการติดตั้งได้ง่ายๆ ไม่กี่คลิก
  • เมื่อติดตั้งแล้วจะได้ไอคอน Surveyor Pocket Tools ออกมา เรียกใช้งานก็ดับเบิ้ลคลิกมาได้จะได้ดังรูปด้านล่าง จะสังเกตว่าโปรแกรมในรุ่นนี้มีโปรแกรมย่อยคือ “File Transform Coordinates

โฟลเดอร์ข้อมูลตัวอย่างทดสอบ

  • ก่อนจะไปต่อผมขอเกริ่นถึงโฟลเดอร์สำหรับข้อมูลตัวอย่างที่จะนำมาทดสอบในโปรแกรม ผมทำไอคอน “Example folder” ไว้ดังรูป

  • เมื่อดับเบิ้ลคลิกเข้าไป ข้อมูลนี้เมื่อโปรแกรมถูกติดตั้งแล้วจะไปอยู่โฟลเดอร์ที่ซ่อนของวินโดส์ “C:\Users\ชื่อผู้ใช้\Appdata\Roaming\Surveyor Pocket Tools” และจะเห็นโฟลเดอร์ “example folder” ให้ดับเบิ้ลคลิกเข้าไปจะได้ดังรูปด้านล่าง

วิธีการใช้งานโปรแกรม File Transfer Coordinates ในเบื้องต้น

  • เปิดโปรแกรม  File Transfer Coordinates แล้วจะเห็นหน้าตาโปรแกรมดังรูปด้านล่าง

  • หน้าต่างไดอะล็อกจะประกอบไปด้วยพาเนลด้านซ้ายและพาเนลด้านขวา ด้านซ้ายสำหรับนำไฟล์ข้อมูลเข้า (input) และด้านขวาสำหรับแสดงผลลัพธ์ (output) ที่ได้จากการแปลงพิกัด ผมทำตารางข้อมูลสไตล์ลายม้าลายให้ดูสวยงามและอ่านง่าย ด้านซ้ายจัดเป็นสี navy blue  ด้านขวาสีโทนเขียวๆ
  • แต่ละด้านจะประกอบไปด้วย ช่องเลือกระบบพิกัดเลือกกรุ๊ป (Group) ก่อนว่าเป็นพิกัดบนเส้นโครงแผนที่หรือพิกัดภูมิศาสตร์บนทรงรี จากนั้นจะเลือกพื้นหลักฐาน (Datum) ตามมา และสุดท้ายจะระบบพิกัด (System)
  • ส่วนการอ่านไฟล์ข้อมูลค่าพิกัดที่ต้องการแปลงจากไฟล์ CSV จะอ่านมาก่อน แล้วค่อยมากำหนดว่าคอลัมน์ไหนเป็นค่า Northing/Latitude หรือว่า Easting/Longitude ทีหลัง
  • และรูปแบบของมุม ว่าเป็นดีกรี (Degree) หรือ รูปมุมแยกมีทศนิยมที่ฟิลิปดา (DD MM SS.SSSS) หรือมุมแยกแบบทศนิยมที่ลิปดา (DD MM.MMMM)
  • ล่างสุดเป็นตารางสำหรับแสดงไฟล์นำเข้า CSV

  • ส่วน button มีปุ่มไอคอนเปิดไฟล์ข้อมูล ด้านบนขวามือ มีปุ่มสำหรับคำนวณรูปลูกศรชี้ไปทางขวา และปุ่มไอคอนรูปโลโก้ excel สำหรับ export จัดเก็บผลลัพธ์ออกเป็นไฟล์ในรูปแบบ Microsoft Excel

เปิดไฟล์ข้อมูลตัวอย่างเพื่อทดสอบ

  • สถานที่ของโฟลเดอร์ที่เด็บข้อมูลตามที่ผมเกริ่นไปตอนต้นบทความ ที่โปรแกรม File Transform Coordinates คลิกทีไอคอนรูปโฟลเดอร์ด้านบนขวา จะเห็นไดอะล็อกเลือกไฟล์ ให้เลือกปลายทางมาที่โฟลเดอร์ดังกล่าว

  • เลือกไฟล์ “boundary2-utm47n-indian1975.csv” เมื่อเปิดแล้วจะเห็นไออะล็อก แสดงว่าในไฟล์มีอะไรอยู่บ้างดังรูปด้านล่าง จะเห็น header “Name,Northing,Easting” ซึ่งผมในฐานะคนเขียนโปรแกรม แนะนำว่าการมี Header นั้นมาจะใช้งานโปรแกรมได้สะดวกที่สุด แต่การไม่ใส่ก็ไม่ได้ผิดอะไรครับ แต่ต้องคลิกมากกว่าเดิมเมื่อเปิดไฟล์มาแล้ว

แปลงค่าพิกัดจากพื้นหลักฐาน Indian 1975 ไปยังพื้นหลักฐาน WGS84

  • ค่าพิกัดในไฟล์ชุดนี้เป็นค่าพิกัดฉากในระบบ “UTM Zone 47N” อยู่บนพื้นหลักฐาน “Indian 1975”  ต้องการแปลงไประบบพิกัดฉาก “UTM Zone 47N” ของพื้นหลักฐาน “WGS84

  • ขั้นตอนต่อไปก็มาดูข้อมูลเมื่อเขียนลงตารางม้าลายแล้วจะเป็นอย่างไร ด้านซ้ายจะเขียนลงตารางให้ทั้งหมด ถ้ามี header โปรแกรมจะจัดคอลัมน์ให้ แต่ใน Header อย่างเช่นค่าพิกัด “Y” ต้องมีข้อความแบบนี้ “Northing/North/Latitude” ส่วนตารางด้านขวาจะลอกตารางด้านซ้ายมาทั้งหมด แต่มีคอลัมน์ว่าง ที่จะรอการแปลงพิกัดมาใส่ ดูรูปด้านล่าง

ปรับระบบพิกัดต้นทางและปลายทาง

  • จะเห็นว่าตอนนี้ไฟล์ข้อมูลเป็นค่าพิกัดฉาก UTM ตามที่ผมกล่าวไปแล้ว แต่ในโปรแกรมยังเป็น “Geographic Coordinate System” แบบค่าพิกัดเป็น “Latitude/Longitude” อยู่ ดังนั้นเราจะมาเปลี่ยนระบบพิกัดเป็น
    • Group “Projected Coordinate System
    • System “UTM Zone 47N
    • Datum “Indian 1975

พิมพ์ค้นหาระบบพิกัดสะดวกรวดเร็ว

  • ก่อนหน้านี้เวลาค้นหาพื้นหลักฐาน Datum หรือว่า Projection เช่นพื้นหลักฐาน WGS84 ค่อนข้างจะใช้เวลาเพราะมีตัวย่อยๆเยอะมาก ต้องใช้เมาส์สกรอลล์ สลับกับการใช้คีย์บอร์ดกดตัวหน้าช่วย ตอนนี้ผมปรับปรุงเพิ่มตัวช่วย ขอแค่รู้คำก็ให้พิมพ์ลงไป ที่ด้านล่างตรง Datum ผมพิมพ์คำว่า Indian จะมีพื้นหลักฐานที่เกี่ยวข้องกันมาโผล่สามอย่างคือ “Indian 1960” , “Indian 1954”, “Indian 1975” ให้คลิกเลือกที่ “Indian 1975

  • เลือก Group ซึ่งจะมีแค่สองอย่าง ไม่ต้องพิมพ์ใช้เมาส์คลิกไปที่ “Projected Coordinate System” จากนั้นมาดูที่ System ซึ่ง Indian 1975 จะมี UTM แค่สองโซนคือ 47N กับ 48N ดังนั้นใช้เมาส์คลิกเลือก “Indian 1975 / UTM zone 47N

  • มาดูปลายทางด้นขวาบ้าง จะเปลี่ยนระบบพิกัดเช่นเดียวกัน ดังนี้
    • Group “Projected Coordinate System
    • System “UTM Zone 47N
    • Datum “WGS84
  • ที่ชอง Datum ผมพิมพ์ world จะเห็นมีสามรายการที่เกี่ยวข้องขึ้นมาดังรูป เราเลือก “World Geodetic System 1984

  • จะได้ผลลัพธ์ดังนี้

  • ต่อไปจะเลือกเส้นโครงแผนที่ ที่ช่อง System ผมพิมพ์ 47N จะเห็นมีรายการที่เกี่ยวข้องขึ้นมาสองรายการ รายการแรกที่เราต้องการคือ “WGS84 / UTM zone 47N” ส่วนรายการที่สองเป็นระบบพิกัดของเมียนมา ที่ใช้กับงานขุดเจาะแก๊สในอ่าวเมาะตะมะ ชื่อเต็มๆคือ “Moattama 92 / UTM zone 47N”

  • ต่อไปเลือกคอลัมน์ที่จะให้ผลลัพธ์ไปออกที่ตาราง จะเห็นคอลัมน์ว่างๆมารออยู่ 4 คอลัมน์ ส่วนชื่อเรียกเป็นหมายเลขคือ Col4, Col5, Col6 และ Col7 เรียงกันไป ตัวนี้โปรแกรมเลือกมาให้เป็นค่าปริยายครับ
  • ความหมายคือต้องการให้ค่า Northing ไปเขียนที่  คอลัมน์ที่ 4 (Col4) ต้องการให้ค่า Easting ไปเขียนที่ คอลัมน์ 5 (Col5) เป็นต้น แต่ถ้าผู้ใช้จะสลับก็ได้นะครับ

  • การตั้งรูปแบบมุมไม่มีครับ เพราะต้องการแปลงค่าพิกัดแค่ระบบพิกัดฉากเท่านั้น

คำนวณ Grid Scale Factor & Convergence

  • ผมได้เพิ่มรายการคำนวณ Grid Scale Factor (GSF) & Convergence สูตรคำนวณสองอย่างนี้ไม่มีไลบรารีตัวไหนทำให้ ต้องมาเขียนเอง เปรียบเทียบค่า Grid Scale Factor แล้วกับโปรแกรม Blue marble “Geographic Calculator” และ Trimble “Coordinates Calculator” ค่าใกล้เคียงกันมาก แตกต่างกันที่ทศนิยมที่ 9 
  • แต่กับ Convergence ทำให้ผมแปลกใจเนื่องจาก สองโปรแกรม Geographic Calculator & Coordinates Calculator และที่ผมคำนวณมา  ให้ค่าที่ต่างกันที่ทศยิมที่ 4 ทำไมเป็นอย่างนั้น ทำให้ผมไม่มั่นใจ ทั้งๆที่สูตรมีแค่บรรทัดเดียวเท่านั้น
  • เพิ่มเติมอีกนิดว่า Convergence เป็นมุมต่างระหว่างทิศเหนือจริงกับทิศเหนือในระบบพิกัดฉาก สมัยเมื่อ 20 กว่าปีที่แล้ว ที่ยังไม่มี GPS/GNSS ผมเคยไปรังวัดอะซิมัทภาคทิศจากดาวและดวงอาทิตย์ จำได้ว่าใช้มุม convergence ที่คำนวณด้วยมือสมัยนั้น แล้วเอาค่าอะซิมัทที่รังวัดมาได้มาลบออกด้วยมุม convergence จะได้อะซิมัทบนระบบพิกัดฉาก แล้วสมัยนี้เอาไปใช้อะไรกันบ้าง นึกไม่ออกครับ

คำนวณและแสดงผล

  • ตอนนี้ตั้งค่าทุกอย่างพร้อมแล้วจะทำการคำนวณ ก็คลิกที่ไอคอนรูปลูกศร คำนวณแปลงพิกัดจากซ้ายไปขวา ได้ผลลัพธ์ดังรูปด้านล่าง จะได้ค่าพิกัดฉาก UTM บน WGS84 พร้อมคำนวณ Grid Scale Factor และ Convergence

จัดเก็บไฟล์ผลลัพธ์ในรูปแบบ Excel

  • เมื่อคำนวณแล้ว สิ่งที่ต้องการก็คือไฟล์ที่เก็บข้อมูลค่าพิกัดที่ได้จากการแปลงพิกัดในรูปแบบ Excel ซึ่งเป็นโปรแกรมสามัญประจำเครืองคอม คลิกที่ไอคอนรูป Excel + ลูกศรสีเหลือง โปรแกรมจะแสดงไดอะล็อก เพื่อถามที่จัดเก็บไฟล์ ป้อนชื่อไฟล์
  • มาเปิดดูไฟล์ผลลัพธ์กัน ผมใช้ LibreOffice Calc เพื่อเปิดไฟล์ผลลัพธ์ดังกล่าว  (LibreOffice โปรแกรมเปิดซอร์สเดี๋ยวนี้เก่งกว่าแต่ก่อนมาก)

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

ตัวอย่างที่ 2

  • จะแปลงพิกัดจากระบบพิกัดฉาก UTM zone 48N ของ Lao National Datum 1997 ไปยังค่าพิกัดภูมิศาสตร์ บน WGS84
  • เปิดไฟล์ข้อมูลอีกครั้งที่โฟลเดอร์ไฟล์ตัวอย่างข้อมูล “example folder” เลือกไฟล์ “coordinates-lao 1997-utm 48n.csv

  • ในไฟล์นี้ผมใส่  header ไว้ด้วย

  • เลือกระบบพิกัดผมพิมพ์คำว่า lao โปรแกรมจะเลือกที่เกี่ยวข้องมาให้ดังรูป จะเห็นมีไอเท็ม “Lao 1993” และ “Lao National Datum 1997” เลือกอย่างหลังครับ

  • เลือกระบบพิกัดให้ได้ตามรูปด้านล่าง และรูปแบบของมุมเลือก DD MM SS.SSSS

  • คลิกที่ไอคอนรูปลูกศรเพื่อทำการคำนวณ จะได้ค่าพิกัดภูมิศาสตร์ บน WGS84

  • จัดเก็บเข้าไฟล์ excel แล้วเปิดดูที่ไฟล์ จะเห็นที่คอลัมน์ “Latitude”, “Longitude” ที่โปรแกรมคำนวณมาให้

  • ฟีเจอร์คำนวณ Scale Factor และ Convergence ได้เพิ่มไปในโปรแกรม Transform Coordinates ที่คำนวณจุดต่อจุด และการพิมพ์เพื่อเลือกระบบพิกัดก็ทำได้เช่นเดียวกัน

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