GCP และงานสนาม

ทำไมต้องรู้

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

ตำแหน่งภาพ GCP และจุดตรวจ

แยกหลักฐานเรื่องตำแหน่งเป็นสามกลุ่ม แต่ละกลุ่มต้องมีที่มา หน่วย ระบบอ้างอิง และสถานะคุณภาพ

RTK (real-time kinematic) แก้ตำแหน่งระหว่างบินด้วยข้อมูลจากสถานีฐานหรือเครือข่าย ส่วน PPK (post-processed kinematic) บันทึกข้อมูลดิบแล้วคำนวณตำแหน่งภายหลัง สถานะ “fixed” บนจอบอกคุณภาพของการแก้สัญญาณ แต่ไม่ครอบคลุมความผิดพลาดจากพิกัดสถานีฐานผิด การตั้งกล้อง หรือการประมวลผลภาพ

วางและรังวัดจุด

เป้าต้องเห็นชัดและชี้จุดเดียวกันได้ในหลายภาพ เลือกพื้นมั่นคง ไม่ถูกใบไม้หรือรถบัง และไม่ใช้เงาเป็นจุดอ้างอิง กระจายจุดให้ครอบคลุมขอบและกลางพื้นที่ รวมทั้งที่สูงและต่ำ ไม่กระจุกอยู่ใกล้จุดขึ้นบินด้านเดียว คู่มือ OpenDroneMap แนะนำให้มี GCP อย่างน้อย 5 จุด (8–10 จุดสำหรับงานใหญ่) และแต่ละจุดควรเห็นในภาพอย่างน้อย 3 ภาพ แต่ไม่มีจำนวนสากลที่ใช้ได้กับทุกงาน จำนวนและการกระจายต้องมาจากขนาดพื้นที่ ภูมิประเทศ และเกณฑ์ที่ตกลงไว้

ระบบความสูง

GNSS ให้ ความสูงเหนือทรงรี ($h$) ซึ่งเป็นผิวทางคณิตศาสตร์ แต่งานสำรวจและวิศวกรรมใช้ ความสูงออร์โทเมตริก ($H$) ที่อ้างอิงระดับน้ำทะเลเฉลี่ย สองค่านี้ต่างกันด้วย ค่าความสูงยีออยด์ ($N$) ตามความสัมพันธ์ที่ NOAA NGS อธิบายไว้

$H=h-N$

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

points = [  # id, บทบาท, E, N, h เสาอากาศ (m), ความสูงเสา p (m), เครื่องหักเสาแล้วหรือยัง
    ("G1", "control", 500.00, 1000.00, 52.430, 2.000, False),
    ("G2", "control", 700.00, 1010.00, 53.105, 2.000, False),
    ("C1", "check", 540.00, 1030.00, 50.610, 2.000, True),
    ("C2", "check", 660.00, 1120.00, 52.880, 1.800, False),
]
GEOID_N = -30.000   # ค่าสมมติสำหรับฝึก ต้องอ่านจาก TGM2017 ตามตำแหน่งจริง

for pid, role, e, n, h_ant, pole, corrected in points:
    h_ground = h_ant if corrected else h_ant - pole
    ortho = h_ground - GEOID_N
    print(f"{pid} {role:<7} h_ground {h_ground:.3f} m  H {ortho:.3f} m")

controls = [p[0] for p in points if p[1] == "control"]
checks = [p[0] for p in points if p[1] == "check"]
print("used for adjustment:", controls, "| kept for checking:", checks)

G1 control h_ground 50.430 m  H 80.430 m
G2 control h_ground 51.105 m  H 81.105 m
C1 check   h_ground 50.610 m  H 80.610 m
C2 check   h_ground 51.080 m  H 81.080 m
used for adjustment: ['G1', 'G2'] | kept for checking: ['C1', 'C2']

ตรวจข้อมูลก่อนกลับ

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

import math

checks = {  # id: (ค่ารังวัด E, N, Z), (ค่าจากแบบจำลอง E, N, Z)
    "C1": ((500.00, 1000.00, 10.00), (500.03, 999.96, 10.08)),
    "C2": ((600.00, 1000.00, 10.00), (600.00, 1000.06, 9.98)),
    "C3": ((500.00, 1100.00, 11.00), (499.92, 1100.06, 11.04)),
}
for pid, (ref, model) in checks.items():
    de, dn, dz = (round(m - r, 3) for m, r in zip(model, ref))
    print(f"{pid}: dE {de:+.2f}  dN {dn:+.2f}  horizontal {math.hypot(de, dn):.2f} m  dZ {dz:+.2f} m")

C1: dE +0.03  dN -0.04  horizontal 0.05 m  dZ +0.08 m
C2: dE +0.00  dN +0.06  horizontal 0.06 m  dZ -0.02 m
C3: dE -0.08  dN +0.06  horizontal 0.10 m  dZ +0.04 m

ปฏิบัติการประจำโมดูล

ข้อผิดพลาดที่พบบ่อย

สรุป

• ตำแหน่งภาพช่วยการประมวลผล GCP ใช้ปรับแบบจำลอง และจุดตรวจกันไว้ตรวจผลอย่างอิสระ
• จุดต้องเห็นชัด กระจายทั่วพื้นที่ และมีทะเบียนที่ตรวจย้อนกลับได้
• หักความสูงเสาครั้งเดียว แล้วแปลงด้วย $H=h-N$ จากแบบจำลองยีออยด์ที่เหมาะสม เช่น TGM2017
• ตรวจข้อมูลก่อนออกจากพื้นที่ และตัดสินใจบินซ่อมจากหลักฐาน

แบบฝึกตรวจความเข้าใจ

1. ความสูงเสาอากาศ 45.620 m เสาสูง 1.800 m เครื่องยังไม่หักเสา ความสูงของจุดบนพื้นเท่าใด
2. $h$ = 43.820 m และ $N$ = −28.500 m (ค่าสมมติ) ความสูงออร์โทเมตริกเท่าใด
3. จุดตรวจมี ΔE = −0.06 m และ ΔN = +0.08 m คลาดแนวราบเท่าใด
4. ถ้านำจุดตรวจไปใช้ปรับแบบจำลอง จะเกิดอะไรกับการประเมินความถูกต้อง
5. สถานะ fixed ของ RTK รับรองความถูกต้องของแผนที่ได้หรือไม่ เพราะอะไร