ตัวเลขพิกัดต้องมาพร้อมความหมาย
พิกัด 100.5, 13.7 อาจเป็นลองจิจูด/ละติจูด แต่ถ้าโปรแกรมอ่านเป็นเมตร จุดจะอยู่ผิดตำแหน่งทันที Coordinate Reference System หรือ CRS บอกวิธีเชื่อมตัวเลขกับตำแหน่งบนโลก รวมถึงระบบอ้างอิง การฉายแผนที่ และหน่วยที่เกี่ยวข้อง
ก่อนรวม orthomosaic, เส้นขอบเขต และจุดสำรวจ ให้บันทึก CRS ต้นทาง ขอบเขตพิกัด หน่วย และที่มาของแต่ละไฟล์ แยก CRS ของชั้นข้อมูลออกจาก CRS ที่โครงการใช้แสดงผล QGIS สามารถแสดงชั้นข้อมูลต่าง CRS ร่วมกันได้ แต่ต้องรู้ CRS ต้นทางที่ถูกต้อง ตาม QGIS Working with Projections
Geographic กับ projected ใช้คนละหน่วย
| ประเภท | รูปแบบที่พบ | ใช้อย่างไรในงานนี้ |
|---|---|---|
| Geographic CRS | ละติจูด/ลองจิจูด หน่วยองศา เช่น WGS 84 EPSG:4326 | แลกเปลี่ยนตำแหน่งและนำเข้าข้อมูล GNSS โดยตรวจลำดับแกน |
| Projected CRS | Easting/Northing มักเป็นเมตร เช่น UTM | วิเคราะห์งานเฉพาะพื้นที่โดยเลือกเขตและ datum ให้เหมาะสม |
| Vertical reference | ความสูงเหนือ ellipsoid, geoid หรือระบบระดับเฉพาะ | ระบุแยกจากพิกัดแนวราบก่อนเทียบ DEM กับหมุดระดับ |
UTM แบ่งเป็นเขตและซีกโลก ไม่ใช่รหัสเดียวสำหรับทุกตำแหน่ง เลือกจากพื้นที่ใช้งานของ CRS และข้อกำหนดงาน ไม่เดาจากตัวเลขที่ดูคล้ายกัน ส่วน datum ใช้กำหนดกรอบอ้างอิงของตำแหน่ง ระบบที่มีชื่อ projection เหมือนกันแต่ datum ต่างกันยังอาจให้พิกัดต่างกัน
สำหรับผู้เริ่มต้นให้อ่านต่อจาก พื้นฐานพิกัดการนำทาง แล้วบันทึกชื่อ CRS เต็มกับรหัส EPSG คู่กัน การอ่านไฟล์บางชนิดใช้ลำดับแกนต่างกัน จึงตรวจว่าคอลัมน์ใดคือ X/longitude และคอลัมน์ใดคือ Y/latitude ก่อนนำเข้า
Assign CRS กับ reproject แก้คนละปัญหา
Assign/Set Layer CRS ใช้เมื่อรู้ว่าตัวเลขในไฟล์อยู่ใน CRS ใด แต่ป้ายกำกับหายหรือผิด มันเปลี่ยนวิธีตีความตัวเลข ไม่ได้คำนวณตัวเลขใหม่ ส่วน reproject ใช้แปลงจาก CRS ที่ถูกต้องไปยัง CRS เป้าหมายและเขียนค่าพิกัดใหม่ การเปลี่ยน project CRS เพื่อดูแผนที่ก็ไม่เท่ากับส่งออกไฟล์ใหม่
ตัวอย่าง: CSV เก็บ lon=100.5000, lat=13.7000 และผู้จัดทำระบุ WGS 84 ให้ตั้ง X เป็น lon, Y เป็น lat และ CRS ต้นทาง EPSG:4326 ถ้าต้องวิเคราะห์ใน WGS 84 / UTM zone 47N ให้ส่งออกสำเนาเป็น EPSG:32647 อย่ากำหนด EPSG:32647 ให้ตัวเลของศาโดยตรง เพราะจะถูกอ่านเป็นระยะเมตร
ก่อนแปลง datum ให้ตรวจ transformation ที่ QGIS เลือก พื้นที่ใช้งาน และ grid ที่ต้องใช้ หากข้อมูลต้นทางไม่มีที่มาพอจะระบุ CRS อย่าไล่เปลี่ยนรหัสจนดูซ้อนกับแผนที่ฐานแล้วถือว่ายืนยันแล้ว ดู QGIS: Datum Transformations
ความสูง GNSS กับระดับน้ำทะเลไม่ใช่ค่าเดียวกัน
ความสูง h อ้างอิงผิว ellipsoid และความสูง H แบบ orthometric อ้างอิง geoid ส่วน N เป็นระยะต่างระหว่างสองผิว โดยความสัมพันธ์หลักคือ H = h − N ต้องใช้แบบจำลองและกรอบอ้างอิงที่เข้ากัน ข้อมูลที่แสดงคำว่า altitude หรือ elevation อย่างเดียวจึงยังตีความไม่ครบ NOAA NGS: Height Relationships
ตัวอย่างสมมติสำหรับฝึกเครื่องหมาย: h = 52.40 m, N = −30.20 m จึงได้ H = 52.40 − (−30.20) = 82.60 m ค่า N นี้สร้างเพื่อฝึกคำนวณ ไม่ใช่ค่า geoid ของจุดในแบบฝึก QGIS และห้ามนำไปแก้ข้อมูลพื้นที่จริง
การเปลี่ยน CRS แนวราบสองมิติไม่ได้ยืนยันว่าแปลงความสูงแล้ว งานจริงต้องทราบ vertical datum, หน่วย, แบบจำลอง geoid, วิธีแปลง และคุณภาพหมุดอ้างอิง หากไม่มีหลักฐาน ให้เก็บความสูงเดิมและระบุว่ายังไม่เทียบระดับข้ามชุด
แบบฝึกบนคอมพิวเตอร์: นำเข้าพิกัดแล้วส่งออกสำเนา
- สร้างไฟล์
crs-practice.csvด้วยข้อความด้านล่าง จุดทั้งหมดเป็นข้อมูลฝึก ไม่ใช่หมุดสำรวจ - ใน QGIS ใช้ Add Delimited Text Layer เลือก X=
lon, Y=latและ EPSG:4326 ตรวจว่าแถวทั้งสามถูกอ่านเป็นจุด - เปิด Layer Properties บันทึก CRS, extent และจำนวนจุด ตั้ง project CRS เป็น EPSG:32647 เพื่อดูการแสดงผลในระบบเป้าหมาย
- คลิกขวาชั้นข้อมูล เลือก Export → Save Features As บันทึก GeoPackage ใหม่และกำหนด EPSG:32647 เมนูอาจต่างเล็กน้อยตามรุ่น
- ตรวจ CRS และ extent ของชั้นที่ส่งออก ค่าพิกัดควรเปลี่ยนเป็น Easting/Northing หน่วยเมตร ขณะที่ตำแหน่งยังตรงกับต้นทางเมื่อแสดงในโครงการเดียวกัน
- บันทึกโครงการ
.qgzและโน้ตว่าขั้นใดแค่เปลี่ยนการแสดงผล ขั้นใดสร้างข้อมูลใหม่ แล้วคำนวณตัวอย่างความสูงแยกในชีต
id,lon,lat
A,100.5000,13.7000
B,100.5010,13.7000
C,100.5000,13.7010
ผลที่ควรได้: CSV ต้นฉบับยังเป็นองศา, GeoPackage มีสามจุดใน EPSG:32647 และทั้งสองชั้นซ้อนตำแหน่งกันได้เมื่อ CRS ถูกต้อง ค่า N สมมติใช้คำนวณแยกได้ 82.60 เมตร การซ้อนกันของสำเนาที่แปลงจากข้อมูลเดียวกันตรวจ workflow การแปลงเท่านั้น ไม่ได้พิสูจน์ความถูกต้องภาคพื้น
ตรวจหาสาเหตุก่อนเลื่อนข้อมูลให้ตรง
ถ้าชั้นข้อมูลอยู่คนละทวีป ให้ตรวจลำดับแกน หน่วย และ CRS ก่อน หากอยู่ใกล้กันแต่เหลื่อมทั้งชุด ให้ตรวจ datum และ transformation ถ้าระดับสูงต่างกันเกือบคงที่ ให้ตรวจ vertical reference แต่รูปแบบการเหลื่อมเพียงอย่างเดียวยังวินิจฉัยสาเหตุไม่ได้
อย่าขยับ raster ให้ตรงกับภาพฐานเพื่อซ่อนปัญหา เพราะแผนที่ฐานเองอาจมีอายุหรือความคลาดเคลื่อนต่างกัน เก็บข้อมูลดิบและวิธีแก้ แล้วใช้จุดอ้างอิงที่ตรวจสอบได้เมื่อต้องยืนยันตำแหน่งจริง
บรรณานุกรมและเรียนต่อ
- QGIS contributors. Working with Projections, QGIS 3.44 User Guide. คู่มือ
- Milbert, D. G., & Smith, D. A. Converting GPS Height into NAVD88 Elevation with the GEOID96 Geoid Height Model. NOAA National Geodetic Survey. หลักการความสูง
ตรวจเอกสาร 8 กันยายน 2569 อ้าง NOAA เพื่ออธิบายความสัมพันธ์ของชนิดความสูง ไม่ได้แนะนำให้ใช้แบบจำลอง geoid ของสหรัฐฯ กับพื้นที่ไทย แบบฝึก QGIS เป็นขั้นตอนให้ผู้เรียนทำ ยังไม่มีผลการรันแนบ เรียนต่อ RTK, PPK, GCP และจุดตรวจสอบ