Giới thiệu dự án

Công tác đo đạc, lập bản đồ địa chính và xây dựng cơ sở dữ liệu đất đai đóng vai trò then chốt trong quản lý nhà nước về tài nguyên, kinh tế và quy hoạch đô thị. Theo thống kê của ngành Quản lý Đất đai, hệ thống hồ sơ địa chính số giúp giảm hơn 70% thời gian xử lý thủ tục hành chính liên quan đến cấp Giấy chứng nhận quyền sử dụng đất (GCNQSDĐ), giải quyết tranh chấp và kiểm kê quỹ đất.

Tại phường Tam Thanh – trung tâm kinh tế, chính trị và dịch vụ của thành phố Lạng Sơn với diện tích tự nhiên 234,11 ha và mật độ dân số lên tới 9.400,18 người/km² (tỷ lệ đất phi nông nghiệp chiếm 80,24%), quá trình đô thị hóa diễn ra nhanh chóng đã dẫn đến biến động ranh giới thửa đất trên 80% so với bộ bản đồ đo vẽ từ năm 2003. Bản đồ giấy cũ bị xuống cấp cơ học, sai lệch ranh giới thực địa và không đáp ứng được yêu cầu chuẩn hóa cơ sở dữ liệu đất đai theo quy định mới.

[Máy Toàn đạc Leica / GPS RTK ComNav T300] 
   [MicroStation V8i & Gcadas 2017] 
[Tạo Topology & Sửa lỗi]    [Gán nhãn, số thửa, loại đất]
 [Tờ Bản đồ Địa chính số 11 (1:500) & Hồ sơ địa chính chuẩn]

Mục tiêu dự án

  1. Khảo sát tổng thể điều kiện tự nhiên, kinh tế - xã hội và đánh giá hiện trạng quản lý sử dụng đất tại phường Tam Thanh, thành phố Lạng Sơn.
  2. Xây dựng, đo đạc và bình sai lưới khống chế đo vẽ gồm 46 điểm đường chuyền kinh vĩ dựa trên 4 điểm địa chính gốc cấp cao (LS27, LS28, LS31, LS32) trong hệ quy chiếu quốc gia VN-2000.
  3. Ứng dụng tích hợp công nghệ toàn đạc điện tử và các giải pháp tin học trắc địa (DPSurvey v2.7, MicroStation V8i, Gcadas 2017) để đo vẽ chi tiết và biên tập số hóa hoàn chỉnh tờ bản đồ địa chính số 11 tỷ lệ 1:500.
  4. Đánh giá tính chính xác, phân tích ưu nhược điểm của quy trình công nghệ số và đề xuất giải pháp tối ưu hóa tiến độ đo vẽ địa chính đô thị.

Giải pháp kỹ thuật và phạm vi

  • Phương pháp tiếp cận: Ứng dụng quy trình đo đạc kỹ thuật số tích hợp: sử dụng công nghệ định vị vệ tinh GPS RTK kết hợp máy toàn đạc điện tử Leica để đo chi tiết theo phương pháp tọa độ cực; xử lý bình sai lưới bằng phần mềm DPSurvey 2.7; tự động hóa biên tập không gian và thuộc tính thửa đất trên nền tảng MicroStation V8i và tiện ích Gcadas 2017.
  • Phạm vi thực hiện: Tờ bản đồ địa chính số 11, tỷ lệ 1:500, diện tích 6,25 ha (kích thước khung trong tiêu chuẩn 50 x 50 cm tương ứng ô vuông thực địa 0,25 x 0,25 km) tại phường Tam Thanh, thành phố Lạng Sơn.
  • Yêu cầu kỹ thuật đạt được: Đảm bảo sai số vị trí điểm sau bình sai $m_p \le 0,05\text{ m}$, sai số khép góc $f_\beta \le 2m_\beta\sqrt{n}$, tuân thủ nghiêm ngặt quy định tại Thông tư số 25/2014/TT-BTNMT của Bộ Tài nguyên và Môi trường.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Trước khi triển khai, việc đánh giá so sánh các phương pháp thành lập bản đồ địa chính là yếu tố quyết định chất lượng và tính khả thi của dự án:

Tiêu chí so sánh Đo vẽ bàn đạc / Kinh vĩ quang cơ Đo vẽ ảnh hàng không / UAV phối hợp Toàn đạc điện tử kết hợp CAD/GIS (Dự án chọn)
Độ chính xác mặt bằng Thấp ($> 0,15\text{ m}$), dễ tích lũy sai số đọc Trung bình - Cao ($0,05 - 0,10\text{ m}$), phụ thuộc GCP Rất cao ($< 0,03 - 0,05\text{ m}$), kiểm soát trực tiếp từng điểm góc
Khả năng đo khu vực đô thị nén Kém, tốc độ đo chậm Hạn chế do nhà cao tầng, ban công che khuất ranh đất Tối ưu, linh hoạt đặt trạm máy và gương phụ luồn lách
Mức độ tự động hóa dữ liệu Thủ công 100%, ghi sổ giấy Xử lý ảnh tự động nhưng cần biên vẽ bổ sung lớn Trút số liệu số (.GSI, .txt), tự động tính tọa độ và phun điểm
Chi phí & tính khả thi Nhân công cao, thời gian kéo dài Chi phí bay chụp cao, khó xin cấp phép bay đô thị Tận dụng thiết bị sẵn có, hiệu quả kinh tế cao nhất

Phân loại yêu cầu hệ thống theo mô hình MoSCoW

  • Must have (Bắt buộc): Hệ tọa độ phẳng VN-2000, kinh tuyến trục 107°15', múi chiếu 3° ($k = 0,9999$); bình sai chặt chẽ lưới đường chuyền; tuân thủ quy chuẩn lớp đối tượng (Level 1-63) trong môi trường MicroStation V8i.
  • Should have (Nên có): Tự động phát hiện và xử lý lỗi hình học thửa đất (bắt quá mút, hở đỉnh, đè ranh); tự động đánh số thửa và xuất bảng nhãn thửa đất.
  • Could have (Có thể có): Khả năng chuyển đổi dữ liệu thuộc tính trực tiếp sang hệ thống thông tin đất đai đa mục tiêu (VILG/eData).
  • Won't have (Chưa thực hiện): Khảo sát mô hình 3D địa chính (3D Cadastre) cho các công trình ngầm và không gian cao tầng.

Thiết kế hệ thống

Hệ thống công nghệ được thiết kế theo chuỗi xử lý dữ liệu khép kín:

[Ngoại nghiệp: Leica Total Station / ComNav T300]
     [Tệp dữ liệu thô: .GSI / .txt]
      [Phần mềm DPSurvey Version 2.7]
   (Bình sai lưới, Module TDDC tính XP, YP, H)
        [Tệp dữ liệu chuẩn: *.TXT]
   [Nền tảng CAD: MicroStation V8i (Seed2D.dgn)]
     [Phần mềm chuyên dụng: Gcadas 2017]
  (eMap -> Tạo vùng -> Gán nhãn -> Đóng khung BĐĐC)
   [Bộ sản phẩm: Tờ BĐĐC số 11 & CSDL Thuộc tính]

Ngăn xếp công nghệ (Technology Stack)

  • Thiết bị đo: Máy toàn đạc điện tử Leica (độ phóng đại ống kính 30x, đo cạnh điện tử EDM sai số $2\text{ mm} + 2\text{ ppm}$); Máy định vị vệ tinh GPS RTK ComNav T300.
  • Phần mềm trắc địa: DPSurvey phiên bản 2.7 (xử lý số liệu trắc địa, trút số liệu GSI, tính chuyển tọa độ, bình sai lưới phụ thuộc và tự do).
  • Nền tảng thiết kế & CSDL: Bentley MicroStation V8i (SelectSeries 3/4), Gcadas 2017 (gồm các phân hệ eMap biên tập bản đồ, eCadas lập hồ sơ địa chính, eData quản trị dữ liệu thuộc tính).
  • Tiện ích chuyển đổi: Total Commander quản lý tệp và chuẩn hóa phần mở rộng dữ liệu trắc địa.

Phương pháp luận (Methodology)

Dự án áp dụng quy trình kiểm soát chất lượng 4 bước dựa trên quy chuẩn kỹ thuật quốc gia:

  1. Khảo sát - Thiết kế kỹ thuật: Thu thập bản đồ địa giới hành chính 364 (tỷ lệ 1:10.000), bản đồ địa hình 1994 (1:25.000) và bản đồ hiện trạng sử dụng đất 2014 (1:2.000). Xác định 4 điểm mốc gốc địa chính gồm LS27, LS28, LS31, LS32.
  2. Triển khai ngoại nghiệp: Thiết kế lưới đường chuyền kinh vĩ cấp 1 và cấp 2 gồm 46 điểm khống chế đo vẽ. Sử dụng máy GPS RTK và toàn đạc Leica đo ngắm 2 lượt (đi và về), mỗi lượt 2 nửa lần đo.
  3. Nội nghiệp & Tự động hóa: Bình sai gián tiếp bằng DPSurvey, chuyển tọa độ vào seedfile DGN, dùng Gcadas tự động xây dựng vùng không gian (Topology), lọc lỗi tự động và gán thông tin thuộc tính (mục đích sử dụng, chủ sở hữu, diện tích).
  4. Kiểm tra, nghiệm thu (QA/QC): Thực hiện quy trình 3 cấp kiểm tra theo Thông tư 49/2016/TT-BTNMT, bảo đảm không có sai sót hình học hay nhầm lẫn loại đất.

Triển khai thực hiện và kết quả

Quy trình đo vẽ và thuật toán xử lý

1. Thuật toán tính tọa độ cực các điểm chi tiết ngoại nghiệp

Khi đặt máy tại trạm khống chế đã biết tọa độ $A(X_A, Y_A)$, ngắm định hướng về trạm $B(X_B, Y_B)$ để đưa bàn độ bằng về $0^\circ 00' 00''$, tọa độ của điểm chi tiết $P$ được xác định theo công thức:

$$\alpha_{A-P} = \alpha_{A-B} + \beta$$

$$X_P = X_A + S \cdot \cos\alpha_{A-P}$$

$$Y_P = Y_A + S \cdot \sin\alpha_{A-P}$$

Trong đó: $S$ là khoảng cách ngang từ trạm máy $A$ đến điểm gương $P$; $\beta$ là góc kẹp bằng đo được; $\alpha_{A-P}$ là góc phương vị cạnh đo $A-P$.

2. Xử lý tệp dữ liệu thô (.GSI) từ máy toàn đạc Leica

Tệp dữ liệu ngoại nghiệp dạng GSI-16 được phân tích cú pháp để trích xuất chỉ số khối phục vụ tính toán tọa độ:

import math

def parse_gsi_line(line: str) -> dict:
    """
    Phân tích dòng dữ liệu GSI từ máy toàn đạc Leica.
    Word 11: Point ID; Word 21: Horizontal Angle; Word 22: Vertical Angle; Word 31: Slope Distance
    """
    parts = line.strip().split()
    data = {}
    for part in parts:
        word_id = part[:2]
        if word_id == "11":
            data['point_id'] = part[7:].lstrip('0')
        elif word_id == "21":  # Góc bằng (gon hoặc dms)
            raw_angle = float(part[7:]) / 100000.0
            data['hz_angle'] = raw_angle
        elif word_id == "31":  # Khoảng cách nghiêng (mm)
            data['slope_dist'] = float(part[7:]) / 1000.0
    return data

def compute_polar_coords(x_station: float, y_station: float, 
                         azimuth_station: float, hz_deg: float, dist_m: float) -> tuple:
    """Tính tọa độ mặt bằng điểm chi tiết theo phương pháp tọa độ cực."""
    total_azimuth_rad = math.radians((azimuth_station + hz_deg) % 360)
    x_point = x_station + dist_m * math.cos(total_azimuth_rad)
    y_point = y_station + dist_m * math.sin(total_azimuth_rad)
    return round(x_point, 4), round(y_point, 4)

Ví dụ chuyển đổi bản ghi GSI sang tọa độ phẳng:
Trạm máy KV01 (X = 2417318.041, Y = 527192.115)
Ngắm điểm định hướng LS27 -> Đưa bàn độ về 0°00'00"
Đo điểm chi tiết 101: Góc bằng = 45°12'18", Cạnh ngang = 35.420 m
Kết quả tính toán: X_101 = 2417343.018 m, Y_101 = 527217.251 m

Thử nghiệm và đánh giá độ chính xác

Hệ thống lưới gồm 50 điểm (4 điểm địa chính gốc và 46 điểm kinh vĩ) được bình sai chặt chẽ bằng phần mềm DPSurvey 2.7.

Kết quả bình sai lưới khống chế đo vẽ (Hệ tọa độ VN-2000, Kinh tuyến trục 107°15')

STT Tên điểm mốc Tọa độ X (m) Tọa độ Y (m) Độ cao h (m) Sai số $m_x$ (m) Sai số $m_y$ (m) Sai số vị trí $m_p$ (m)
1 LS27 (Gốc) 2418313.470 527632.180 268.450 - - -
2 LS28 (Gốc) 2418012.115 526985.430 255.120 - - -
3 KV01 2417318.041 527192.115 252.314 0.008 0.007 0.011
4 KV14 2417552.890 527340.612 258.902 0.012 0.011 0.016
5 KV26 (Max) 2417890.120 527510.334 261.445 0.021 0.019 0.028

Các chỉ số kiểm định sai số hình học

  • Số lượng tam giác khép hình: 46 tam giác.
  • Sai số khép tương đối tam giác nhỏ nhất: $1/748.835$ (Tam giác $LS27 - KV04 - KV01$).
  • Sai số khép tương đối tam giác lớn nhất: $1/11.672$ (Tam giác $KV46 - KV45 - KV44$, vẫn vượt xa chỉ tiêu cho phép $1/2.500$).
  • Tỷ số Ratio kiểm định: Lớn nhất đạt $72,670$ ($KV22 - KV21$), nhỏ nhất đạt $2,880$ ($KV05 - LS28$).
  • Chiều dài cạnh trung bình: $S_{\text{tb}} = 489,66\text{ m}$.
            [KV21]
           /      \
      Ratio: 72.67 \
         /          \
  Tam giác khép hình đạt chuẩn 1/748.835

Kết quả biên tập tờ bản đồ số 11 tỷ lệ 1:500

Sau khi hoàn tất khâu đo đạc ngoại nghiệp và bình sai, dữ liệu được nạp vào MicroStation V8i và Gcadas 2017:

  • Xử lý hình học tự động: Thuật toán của Gcadas quét sạch 100% lỗi đồ họa (các lỗi cắt quá mút, lỗi góc hở $d < 0,001\text{ m}$, lỗi đè vạch phân lô).
  • Tạo tâm thửa và gán nhãn: Toàn bộ các thửa đất trong diện tích 6,25 ha được tạo vùng (Polygonize) tự động, gán thông tin số hiệu thửa, diện tích (tính toán chính xác tới $0,1\text{ m}^2$) và mã loại đất theo Thông tư 28/2014/TT-BTNMT (ODT, TSC, DSN, CLN,...).
  • Xuất bản đồ và hồ sơ: Hoàn thiện tờ bản đồ số 11 với đầy đủ khung tiêu chuẩn, lưới tọa độ kilômét, bảng chú dẫn, phân lớp đối tượng giao thông, thủy hệ, công trình kiên cố và ranh giới quy hoạch.

Đổi mới và đóng góp kỹ thuật

So sánh với các giải pháp hiện hành

Chỉ số đánh giá Phương pháp cổ điển (Kinh vĩ + AutoCAD) Phương pháp đo toàn đạc độc lập Chuỗi tích hợp (GPS RTK + Leica + DPSurvey + Gcadas)
Thời gian nội nghiệp 12 - 15 ngày/tờ 1:500 5 - 7 ngày/tờ 1:500 1,5 - 2 ngày/tờ 1:500 (Giảm > 65%)
Tỷ lệ sót lỗi topology Cao ($5 - 8%$ do vẽ tay) Trung bình ($2 - 3%$) $0%$ (Gcadas kiểm tra tự động)
Tự động hóa thuộc tính Thủ công từng nhãn Thủ công bằng AutoLISP Tự động liên kết cơ sở dữ liệu bảng mã
Khả năng cập nhật biến động Khó khăn, phải vẽ lại Cập nhật cục bộ Chỉnh lý tức thời trên các phân lớp DGN

Các cải tiến kỹ thuật nổi bật

  1. Liên thông dòng dữ liệu (Data Pipeline): Thiết lập quy trình chuyển đổi tệp GSI sang TXT chuẩn hóa thông qua DPSurvey, hạn chế tuyệt đối sai sót do nhập liệu thủ công.
  2. Nâng cao độ chính xác đo vẽ đô thị: Kết hợp GPS RTK bắt điểm tọa độ mốc gốc với toàn đạc điện tử đo chi tiết ngóc ngách, loại bỏ ảnh hưởng của hiện tượng đa đường dẫn tín hiệu (multipath) trong không gian phố núi nén chặt.
  3. Chuẩn hóa cấu trúc lớp dữ liệu: Thiết lập hệ thống Level/Color/LineStyle trong MicroStation V8i chuẩn theo Thông tư 25/2014/TT-BTNMT, tạo tiền đề đồng bộ dữ liệu vào hệ sinh thái VILG cấp tỉnh.

Ứng dụng thực tế và triển khai

Kịch bản ứng dụng thực tế

  • Cấp đổi GCNQSDĐ: Cung cấp hồ sơ kỹ thuật thửa đất chính xác, minh bạch ranh giới để giải quyết tình trạng lấn chiếm hành lang đường bộ, đất công trên địa bàn phường Tam Thanh.
  • Quy hoạch và phát triển hạ tầng: Làm nền tảng phân tích không gian cho các dự án mở rộng lòng đường, ngầm hóa hệ thống điện, thoát nước thải dọc tuyến sông Kỳ Cùng.
  • Xử lý tranh chấp đất đai: Cung cấp bằng chứng đo đạc kỹ thuật số có giá trị pháp lý cao, xác định ranh giới sai lệch với độ chuẩn xác từng centimet.

Lộ trình triển khai và tính toán kinh tế

   Tuần 1-2                 Tuần 3-6                Tuần 7-10               Tuần 11-12
(Thu thập tài liệu,     (Đo 46 điểm kinh vĩ,    (Bình sai, sửa topology, (Kiểm tra 3 cấp theo
 4 mốc địa chính)        đo chi tiết thửa đất)   tạo nhãn, đóng khung)    Thông tư 49/2016)
  • Hiệu quả đầu tư (ROI): Giảm 40% chi phí nhân công ngoại nghiệp nhờ tăng tốc độ đo đạc toàn đạc điện tử kết hợp GPS RTK; giảm 60% thời gian xử lý dữ liệu nội nghiệp nhờ tự động hóa bằng Gcadas. Bộ dữ liệu có vòng đời khai thác hiệu quả từ 5 – 10 năm.

Hạn chế và hướng phát triển

Hạn chế kỹ thuật

  • Vùng che khuất tự nhiên: Địa hình phường Tam Thanh có núi đá đan xen khu dân cư, một số điểm góc thửa bị cây cối và nhà tầng che khuất tầm ngắm quang học, buộc phải tăng dày các trạm máy phụ.
  • Tính năng 3D: Quy trình hiện tại chỉ xử lý bản đồ 2D ($X, Y$), chưa tích hợp mô hình độ cao số (DEM) hoặc bản đồ địa chính 3D cho các khu nhà phức hợp cao tầng.

Hướng phát triển

  • Ứng dụng thiết bị quét 3D Laser Scanner mặt đất (TLS) hoặc UAV trang bị cảm biến LiDAR quét xuyên tán cây cho các khu vực địa hình đồi núi phức tạp.
  • Nâng cấp cơ sở dữ liệu lên nền tảng WebGIS và tích hợp kiến trúc điện toán đám mây để người dân có thể tra cứu thông tin quy hoạch trực tuyến qua thiết bị di động.

Đối tượng hưởng lợi

  • Sinh viên & Giảng viên ngành Quản lý đất đai/Trắc địa: Nguồn tài liệu thực tế về tích hợp phần cứng máy đo và phần mềm đồ họa trắc địa chuyên dụng.
  • Đo đạc viên & Kỹ sư GIS: Bộ khung quy trình xử lý lỗi hình học và bình sai tự động áp dụng được ngay cho các dự án đo vẽ đô thị tỷ lệ lớn (1:500, 1:200).
  • Chính quyền địa phương (UBND phường Tam Thanh & TP Lạng Sơn): Sở hữu bộ bản đồ số chuẩn hóa, quản lý hiệu quả 234,11 ha đất đai, ngăn chặn lãng phí tài nguyên và thất thoát thuế đất.
  • Người dân địa phương: Giảm thiểu tranh chấp giáp ranh, rút ngắn thời gian làm thủ tục cấp phép xây dựng và sang tên chuyển nhượng quyền sử dụng đất.

Câu hỏi thường gặp

1. Cần cấu hình phần cứng và phần mềm như thế nào để vận hành quy trình này?

Máy tính cấu hình tối thiểu: CPU Intel Core i5 thế hệ 8 trở lên, RAM 8GB, Card đồ họa rời hoặc tích hợp hỗ trợ DirectX 11, cài đặt hệ điều hành Windows 10/11 64-bit; phần mềm gồm MicroStation V8i (SS3), Gcadas 2017 và DPSurvey 2.7.

2. Độ chính xác của lưới đường chuyền có đáp ứng yêu cầu đo vẽ tỷ lệ 1:500 ở đô thị không?

Hoàn toàn đáp ứng. Theo Bảng 4.6, sai số trung phương vị trí điểm lớn nhất là $m_p = 0,028\text{ m}$ (tại điểm KV26), nhỏ hơn nhiều so với hạn mức cho phép $0,05\text{ m}$ theo Thông tư 25/2014/TT-BTNMT.

3. Quy trình xử lý khi gặp sai số khép góc đường chuyền vượt hạn sai cho phép?

Nếu sai số khép góc $f_\beta > 2m_\beta\sqrt{n}$, tiến hành kiểm tra lại góc đo từng trạm máy trong tệp dữ liệu thô, xác định trạm có chênh lệch giữa 2 nửa lần đo $> 20''$ để tổ chức đo ngắm lại ngay ngoài thực địa trước khi đưa vào bình sai tổng thể.

4. Dữ liệu biên tập từ Gcadas có thể xuất sang định dạng GIS khác như ArcGIS hay QGIS không?

Có. Từ tệp *.dgn của MicroStation/Gcadas, có thể chuyển đổi trực tiếp sang định dạng Shapefile (*.shp) hoặc Geodatabase thông qua công cụ Data Interoperability / FME để tích hợp vào các hệ thống thông tin địa lý GIS chuyên sâu.

5. Chi phí triển khai đo vẽ số hóa cho một phường đô thị quy mô tương tự là bao nhiêu?

Chi phí phụ thuộc vào diện tích tự nhiên, mật độ thửa và tỷ lệ bản đồ. Với khu vực đô thị như phường Tam Thanh (khoảng 60 - 62 tờ tỷ lệ 1:500), định mức dự toán được áp dụng theo Thông tư của Bộ TN&MT về định mức kinh tế - kỹ thuật đo đạc lập bản đồ địa chính, tiết kiệm khoảng 30 - 35% so với phương pháp đo vẽ thủ công truyền thống.

Kết luận

Đề tài đã hoàn thành xuất sắc các mục tiêu nghiên cứu và ứng dụng thực tiễn trong việc thành lập tờ bản đồ địa chính số 11 tỷ lệ 1:500 tại phường Tam Thanh, thành phố Lạng Sơn. Việc kết hợp đồng bộ giữa máy toàn đạc điện tử Leica, công nghệ định vị GPS RTK cùng hệ sinh thái phần mềm DPSurvey 2.7, MicroStation V8i và Gcadas 2017 đã chứng minh tính ưu việt vượt trội về độ chính xác (sai số vị trí mốc $m_p < 0,03\text{ m}$), tự động hóa xử lý dữ liệu và loại bỏ hoàn toàn lỗi cấu trúc hình học thửa đất.

Kết quả của đề tài không chỉ là tài liệu khoa học có giá trị tham khảo chuyên sâu cho sinh viên, kỹ sư ngành Quản lý Đất đai và Trắc địa - Bản đồ, mà còn là sản phẩm kỹ thuật số sẵn sàng tích hợp vào cơ sở dữ liệu đất đai quốc gia, phục vụ thiết thực cho tiến trình chuyển đổi số và hiện đại hóa công tác quản trị đô thị tại Việt Nam.