Giới thiệu dự án

Công tác quản lý nhà nước về đất đai đóng vai trò then chốt trong phát triển kinh tế - xã hội. Tại Việt Nam, quá trình tích tụ ruộng đất, dồn điền đổi thửa và cấp đổi Giấy chứng nhận quyền sử dụng đất (GCNQSDĐ) theo Luật Đất đai đặt ra yêu cầu cấp bách về việc hiện đại hóa công tác trắc địa bản đồ. Theo thống kê từ Bộ Tài nguyên và Môi trường (TN&MT), việc ứng dụng công nghệ đo đạc số hóa giúp giảm thiểu tới 40% sai sót pháp lý và tăng tốc độ xử lý thủ tục hành chính đất đai hơn 60% so với phương pháp thủ công.

Tại xã Mai Trung, huyện Hiệp Hòa, tỉnh Bắc Giang (diện tích tự nhiên 1.050,14 ha, bao gồm 670,21 ha đất nông nghiệp và 379,8 ha đất phi nông nghiệp), việc phân mảnh thửa đất sau nhiều giai đoạn canh tác cùng với sự biến động ranh giới thực địa gây khó khăn lớn cho công tác quản lý. Phương pháp đo đạc truyền thống bằng máy kinh vĩ hoặc toàn đạc điện tử đơn thuần bộc lộ nhiều điểm nghẽn: phụ thuộc góc nhìn thông hướng, tốn nhiều nhân lực dẫn lưới khống chế và thời gian đo kéo dài.

Đề tài "Ứng dụng công nghệ GNSS-RTK trong công tác đo đạc bản đồ địa chính phục vụ công tác dồn điền đổi thửa, kê khai đăng ký, cấp giấy chứng nhận quyền sử dụng đất tại xã Mai Trung, huyện Hiệp Hòa, tỉnh Bắc Giang" được thực hiện nhằm giải quyết dứt điểm các bài toán trên.

Mục tiêu dự án

  1. Mục tiêu tổng quát: Xây dựng cơ sở dữ liệu không gian và thành lập hệ thống bản đồ địa chính chính quy tỷ lệ 1:1000 phục vụ đăng ký, kê khai, cấp GCNQSDĐ tại xã Mai Trung, huyện Hiệp Hòa, tỉnh Bắc Giang.
  2. Mục tiêu cụ thể:
    • Ứng dụng thiết bị GNSS-RTK (Real-Time Kinematic) đa tần số thu nhận và xử lý tọa độ thực địa cho 39 điểm khống chế đo vẽ và 1.500 điểm đo chi tiết.
    • Ứng dụng nền tảng MicroStation V8i kết hợp phần mềm biên tập địa chính chuyên dụng gCadas 2015 để tự động hóa xử lý dữ liệu đồ họa, tạo vùng, đánh số thửa và xuất hồ sơ địa chính.
    • Đảm bảo độ chính xác vị trí điểm ranh giới thửa đất đạt chuẩn sai số mặt phẳng $\le 15\text{ cm}$ đối với khu dân cư và $\le 22,5\text{ cm}$ đối với đất nông nghiệp theo quy định tại Thông tư 25/2014/TT-BTNMT.

Phương pháp tiếp cận và giới hạn phạm vi

  • Giải pháp: Sử dụng trạm máy tĩnh (Base Station) đặt tại điểm khống chế địa chính quốc gia và máy động (Rover Station) thu sóng vi sai qua sóng vô tuyến UHF công suất 25W, kết hợp quy trình xử lý topo tự động trên phần mềm gCadas.
  • Phạm vi không gian: Địa bàn xã Mai Trung, huyện Hiệp Hòa, tỉnh Bắc Giang.
  • Phạm vi thời gian: Số liệu đo đạc thực tế được chuẩn hóa và biên tập hoàn thiện tờ bản đồ địa chính số DC28 tỷ lệ 1:1000.
  • Giới hạn kỹ thuật: Công nghệ GNSS-RTK phụ thuộc vào độ thông thoáng vệ tinh (chỉ số PDOP < 3); các góc khuất tán cây rậm rạp hoặc khe tường hẹp được phối hợp đo bù bằng thước thép và máy toàn đạc điện tử.

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

Phân tích hiện trạng và so sánh giải pháp

Trước khi ứng dụng GNSS-RTK, công tác thành lập bản đồ địa chính nông thôn chủ yếu dựa vào máy toàn đạc điện tử (Total Station) hoặc ảnh hàng không (Photogrammetry). Bảng so sánh dưới đây làm rõ tính ưu việt của giải pháp lựa chọn:

Tiêu chí đánh giá Đo truyền thống (Toàn đạc điện tử) Ảnh hàng không (Photogrammetry) Ứng dụng công nghệ GNSS-RTK
Yêu cầu thông hướng Bắt buộc phải thông hướng giữa các trạm máy Cần thông thoáng từ trên cao, dễ bị che khuất Không cần thông hướng giữa các điểm đo chi tiết
Nhân lực vận hành 3 - 4 người/tổ (người đứng máy, người gương) Đội bay, tổ bay, đội đo khống chế ngoại nghiệp 1 - 2 người/tổ (cầm sào Rover gắn sổ tay điện tử)
Tốc độ thu thập 150 - 250 điểm/ngày Nhanh toàn diện nhưng mất thời gian giải đoán ảnh 600 - 900 điểm/ngày (tức thời 1 - 3 giây/điểm)
Độ chính xác mặt phẳng $\pm (2\text{mm} + 2\text{ppm})$ $\pm 10\text{ cm} \div 30\text{ cm}$ (phụ thuộc GSD và GCP) $\pm (10\text{mm} + 1\text{ppm})$ (đạt chuẩn Thông tư 25)
Chi phí triển khai lưới Rất cao (phải dẫn đường chuyền cấp 1, cấp 2) Cao (chi phí bay chụp và bình sai lưới khống chế ảnh) Tối ưu (đặt trực tiếp Base lên mốc gốc, không cần lập lưới đo vẽ dày đặc)

Yêu cầu người dùng theo mô hình MoSCoW

  • Must have (Bắt buộc): Chuyển đổi chính xác tọa độ từ WGS-84 sang VN-2000 (kinh tuyến trục $107^\circ 00'$, múi chiếu $3^\circ$). Tự động làm sạch dữ liệu hình học (lỗi dangles, undershoots, duplicate lines). Xuất đúng biểu mẫu: Trích lục thửa đất, Bản mô tả ranh giới mốc giới, Đơn đăng ký cấp GCNQSDĐ.
  • Should have (Nên có): Tự động phát hiện lỗi trùng đè ranh thửa (Overlaps), kiểm tra sai số khép cạnh, gán thông tin thuộc tính hàng loạt từ file Excel/Access.
  • Could have (Có thể có): Liên kết dữ liệu không gian trực tiếp với hệ quản trị cơ sở dữ liệu địa chính tập trung VBDLIS/VILIS.
  • Won't have (Chưa thực hiện đợt này): Quét 3D laser mặt đất (TLS) và tích hợp hệ thống WebGIS 3D thời gian thực.

Thiết kế hệ thống và công nghệ sử dụng

Hệ thống xử lý được thiết kế theo quy trình khép kín từ khâu thu nhận sóng vệ tinh đến xuất hồ sơ pháp lý:

                                                    [gCadas 2015 Module]
  • Technology Stack:
    • Phần cứng thu nhận: Máy GNSS RTK Kolida K5+ đa vệ tinh (GPS L1/L2/L5, GLONASS, BeiDou, Galileo), công suất phát Radio trong 25W, khoảng cách phát vi sai lên đến 12 km.
    • Nền tảng đồ họa CAD: Bentley MicroStation V8i (SELECTseries 3) phiên bản 08.11.09.
    • Phần mềm nghiệp vụ: gCadas 2015 (tích hợp các phân hệ eMap, eCadas, eData) phát triển bởi Công ty TNHH Trắc địa và Công nghệ Toàn Việt, tuân thủ Thông tư 28/2014/TT-BTNMT và Thông tư 25/2014/TT-BTNMT.
    • Cơ sở dữ liệu thuộc tính: Microsoft Access (.mdb) liên kết qua chuẩn kết nối ODBC.

Cơ sở dữ liệu và sơ đồ quan hệ thuộc tính

Implementation và kết quả

Quy trình triển khai kỹ thuật và thuật toán cốt lõi

1. Thuật toán chuyển đổi tọa độ Helmert 7 tham số (WGS-84 sang VN-2000)

Tọa độ thu nhận từ vệ tinh ở dạng trắc địa $(B, L, H){\text{WGS-84}}$ được chuyển đổi sang không gian trực chuẩn ba chiều $(X, Y, Z){\text{WGS-84}}$, sau đó áp dụng phép biến đổi 7 tham số của Cục Đo đạc và Bản đồ Việt Nam:

$$\begin{bmatrix} X \ Y \ Z \end{bmatrix}{\text{VN-2000}} = \begin{bmatrix} \Delta X \ \Delta Y \ \Delta Z \end{bmatrix} + (1 + m) \begin{bmatrix} 1 & -\omega_Z & \omega_Y \ \omega_Z & 1 & -\omega_X \ -\omega_Y & \omega_X & 1 \end{bmatrix} \begin{bmatrix} X \ Y \ Z \end{bmatrix}{\text{WGS-84}}$$

Sau đó, tọa độ được chiếu lên mặt phẳng kinh tuyến trục $107^\circ 00'$ (tỉnh Bắc Giang), múi chiếu $3^\circ$ với hệ số co giãn kinh tuyến trục $k_0 = 0,9999$.

2. Cấu trúc chuẩn hóa Layer/Level trong MicroStation V8i

Để đảm bảo bản đồ tuân thủ tiêu chuẩn kỹ thuật số liệu địa chính, hệ thống layer được quy định nghiêm ngặt:

Level 10: Ranh giới thửa đất (Color: 0 - White/Black, Weight: 0, Style: 0)
Level 15: Ranh giới nhà và công trình chính (Color: 2 - Green, Weight: 0, Style: 2)
Level 23: Tim đường và ranh giới giao thông (Color: 3 - Red, Weight: 0, Style: 0)
Level 32: Hệ thống mương, thủy lợi nội đồng (Color: 1 - Blue, Weight: 0, Style: 0)
Level 41: Nhãn thông tin thửa đất (Số tờ, số thửa, diện tích, loại đất)
Level 62: Khung bản đồ, lưới km và bảng chú dẫn

3. Xử lý Topology tự động và xây dựng vùng thửa đất

Dữ liệu đo ngoại nghiệp định dạng ASCII STT, X, Y, H, Code được import tự động qua script xử lý topology:

# Script mô phỏng thuật toán xử lý Topology và đóng vùng thửa đất (Clean & Build Algorithm)
class CadastralTopologyEngine:
    def __init__(self, tolerance=0.02):
        self.tolerance = tolerance # Dung sai bắt điểm 2cm theo Thông tư 25
        self.nodes = []
        self.edges = []
        self.polygons = []

    def snap_nodes(self, point_a, point_b):
        """Khử sai số undershoot/overshoot bằng thuật toán Snapping."""
        distance = ((point_a['x'] - point_b['x'])**2 + (point_a['y'] - point_b['y'])**2)**0.5
        if 0 < distance <= self.tolerance:
            point_b['x'] = point_a['x']
            point_b['y'] = point_a['y']
            return True
        return False

    def build_centroids(self, parcel_lines):
        """Xác định tâm thửa (Centroid) từ danh sách các chuỗi đường khép kín."""
        for loop in parcel_lines.find_closed_loops():
            area = loop.calculate_shoelace_area()
            centroid_x, centroid_y = loop.calculate_barycenter()
            self.polygons.append({
                "centroid": (centroid_x, centroid_y),
                "area_m2": round(area, 2),
                "status": "VALID_PARCEL"
            })
        return self.polygons

Kết quả kiểm định độ chính xác thực tế

Đề tài đã kiểm tra đối soát 39 điểm khống chế đo vẽ và 1.500 điểm chi tiết tại tờ bản đồ DC28 xã Mai Trung. Độ lệch vị trí mặt phẳng được đánh giá dựa trên công thức sai số trung phương (RMS Error):

$$M_p = \pm \sqrt{\frac{\sum (\Delta X_i^2 + \Delta Y_i^2)}{n}}$$

Nhóm đối tượng kiểm tra Số lượng điểm kiểm tra Sai số vị trí cho phép (TT 25/2014) Sai số thực tế đạt được ($M_p$) Đánh giá chất lượng
Điểm khống chế đo vẽ 39 điểm $\le 10,0\text{ cm}$ $\pm 1,42\text{ cm}$ Đạt yêu cầu 100%
Ranh thửa đất thổ cư 420 điểm $\le 15,0\text{ cm}$ $\pm 4,18\text{ cm}$ Xuất sắc
Ranh thửa đất nông nghiệp 1.080 điểm $\le 22,5\text{ cm}$ $\pm 6,35\text{ cm}$ Xuất sắc
Sai số tương hỗ cạnh ngắn (<5m) 150 cạnh $\le 4,0\text{ cm}$ $\pm 1,85\text{ cm}$ Đạt chuẩn tuyệt đối

Đổi mới và đóng góp

  1. Rút ngắn chu kỳ sản xuất bản đồ: Việc ứng dụng trạm GNSS-RTK phát sóng UHF 25W bán kính 12 km cho phép một tổ đo 2 người hoàn thành từ 600 - 900 điểm chi tiết/ngày, tăng 300% năng suất so với đo toàn đạc điện tử thông thường.
  2. Tự động hóa 100% khâu biên tập hồ sơ pháp lý: Thay vì vẽ và cắt trích lục thủ công trên MicroStation SE, module gCadas 2015 cho phép xuất hàng loạt:
    • Sơ đồ hình thể thửa đất (SĐHT).
    • Phiếu xác nhận kết quả đo đạc hiện trạng.
    • Bản mô tả ranh giới, mốc giới thửa đất theo mẫu quy chuẩn của Bộ TN&MT.
    • Hồ sơ đăng ký cấp GCNQSDĐ với độ chính xác diện tích tới $0,01\text{ m}^2$.
  3. Chuẩn hóa quy trình dồn điền đổi thửa: Đưa ra giải pháp đo định vị ranh giới mới ngay sau khi cày ải thực địa, giúp giải quyết tranh chấp ranh giới giữa các hộ gia đình với độ tin cậy tuyệt đối.
TIẾT KIỆM THỜI GIAN VÀ NHÂN LỰC:

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

Kịch bản triển khai tại địa phương

  • Bước 1: Khảo sát và đặt Base: Đặt trạm tĩnh Base tại điểm địa chính cơ sở 053410 tại trụ sở UBND hoặc vị trí cao thoáng, nhập thông số tọa độ VN-2000 gốc, kết nối máy phát Radio ngoài công suất 25W.
  • Bước 2: Đo chi tiết ranh thửa (Rover): Nhân viên đo đạc cầm sào Rover di chuyển theo ranh giới dồn điền đổi thửa của từng thôn (Trung Hòa, Trung Hưng, Bách Nhẫn). Khi sổ tay báo trạng thái Fix ($HRMS \le 0,015\text{m}$, $VRMS \le 0,025\text{m}$), bấm lưu tọa độ và gán Code (RANH, DUONG, MUONG, NHA).
  • Bước 3: Biên tập nội nghiệp: Cuối ngày trích xuất dữ liệu dạng text, import vào gCadas trên MicroStation V8i để kiểm tra hình học, sửa lỗi nút topo và gán thông tin chủ sử dụng đất.

Hiệu quả kinh tế và thời gian hoàn vốn (ROI)

  • Chi phí thiết bị: Đầu tư 01 bộ máy GNSS-RTK Kolida K5+ (Base + Rover) khoảng 120 - 150 triệu VNĐ.
  • Hiệu quả thu hồi: Giảm chi phí thuê nhân công dẫn lưới (tiết kiệm 45 triệu VNĐ/xã), rút ngắn 15 ngày công thực địa cho một đơn vị xã quy mô 1.000 ha. Thời gian hoàn vốn đầu tư thiết bị ước tính dưới 4 dự án đo vẽ cấp xã.

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

  • Hạn chế kỹ thuật:
    • Tín hiệu vệ tinh bị gián đoạn (mất trạng thái Fix sang Float) tại các khu vực dưới tán cây rậm rạp hoặc các ngõ hẹp có nhà cao tầng san sát.
    • Sóng Radio nội bộ UHF có thể bị cản trở bởi địa hình đồi gò phức tạp, đòi hỏi phải nâng cao cột ăng-ten Base hoặc dùng trạm lặp (Repeater).
  • Hướng phát triển:
    • Chuyển đổi phương thức truyền số gia hiệu chỉnh từ Radio UHF sang mạng trạm định vị vệ tinh quốc gia VNGEONET (NTRIP/VRS) qua sóng 4G/5G, loại bỏ nhu cầu thiết lập trạm Base tĩnh riêng biệt.
    • Tích hợp công nghệ bay quét không người lái UAV-RTK kết hợp cảm biến LiDAR để quét hiện trạng địa hình diện rộng, sau đó dùng GNSS-RTK đo kiểm tra và bổ sung các điểm ẩn.

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

  • Sinh viên ngành Quản lý đất đai / Trắc địa: Nắm vững quy trình khép kín từ lý thuyết trắc địa cao cấp (Helmert, VN-2000) đến thao tác phần mềm chuyên ngành (MicroStation, gCadas).
  • Kỹ sư khảo sát & Đơn vị tư vấn đo đạc: Tiếp cận giải pháp kỹ thuật tối ưu hóa nhân công, loại bỏ sai số luỹ tiến trong đo dẫn lưới khống chế.
  • Cán bộ địa chính cấp xã/huyện: Nhận bàn giao bộ cơ sở dữ liệu số hóa đồng bộ, sạch lỗi topo, dễ dàng cập nhật biến động đất đai và cấp đổi GCNQSDĐ nhanh chóng.
  • Người dân địa phương: Nhận GCNQSDĐ với sơ đồ ranh giới chính xác tuyệt đối, tránh tranh chấp ranh giới sau công tác dồn điền đổi thửa.

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 trơn tru hệ thống gCadas và MicroStation V8i?

Hệ thống yêu cầu máy trạm tính toán chạy hệ điều hành Windows 7/10/11 (64-bit), CPU tối thiểu Intel Core i5 thế hệ 8 trở lên, RAM 8GB (khuyến nghị 16GB để xử lý mượt các file DGN nặng trên 50MB), ổ cứng SSD lưu trữ dữ liệu và cài đặt MicroStation V8i SS3 cùng bộ cài gCadas 2015 có khóa cứng bản quyền.

2. Khi gặp hiện tượng mất trạng thái Fix ở khu vực ngõ hẹp thì xử lý ra sao?

Cần chuyển sang chế độ đo lệch tâm (Offset Measurement): Sử dụng thước thép kết hợp đo 2 điểm phụ ngoài khoảng thoáng có trạng thái Fix để xác định tọa độ giao hội điểm góc nhà/ranh giới khuất, hoặc phối hợp máy toàn đạc điện tử đặt trạm trên 2 điểm RTK đã được Fix chuẩn.

3. Tọa độ xuất từ sổ tay RTK có cần bình sai lại trong phần mềm chuyên dụng không?

Không cần bình sai mạng lưới đo vẽ chi tiết. Do máy tính trạm Base liên tục phát tín hiệu cải chính tọa độ thời gian thực đến Rover, dữ liệu tọa độ lưu trên Rover đã là tọa độ độ chính xác cao trong hệ VN-2000 (kinh tuyến trục địa phương), có thể nạp trực tiếp vào CAD.

4. gCadas 2015 xử lý như thế nào đối với các thửa đất có nhiều đỉnh phức tạp?

Phần mềm sử dụng thuật toán tính diện tích giải tích Gauss thông qua tọa độ các đỉnh ranh giới khép kín $S = \frac{1}{2} \left| \sum (X_i Y_{i+1} - X_{i+1} Y_i) \right|$, đảm bảo tính diện tích chính xác tuyệt đối đến hàng phần trăm mét vuông bất kể hình dạng thửa đất đa giác lồi hay lõm.

5. Dữ liệu biên tập từ MicroStation/gCadas có tương thích với hệ thống VBDLIS quốc gia không?

Hoàn toàn tương thích. Bản đồ số dạng .DGN và cơ sở dữ liệu thuộc tính .MDB có thể được chuyển đổi trực tiếp sang định dạng chuẩn GeoDatabase (.GDB) hoặc Shapefile (.SHP) thông qua các công cụ ETL của GIS để tích hợp đồng bộ vào hệ thống cơ sở dữ liệu đất đai quốc gia VBDLIS.

Kết luận

Đề tài tốt nghiệp của tác giả Hoàng Văn Đội đã chứng minh tính khoa học, hiệu quả kinh tế và giá trị thực tiễn vượt trội của việc ứng dụng công nghệ GNSS-RTK (Kolida K5+) kết hợp phần mềm MicroStation V8i và gCadas 2015 trong công tác thành lập bản đồ địa chính tỷ lệ 1:1000 tại xã Mai Trung, huyện Hiệp Hòa, tỉnh Bắc Giang.

Dự án không chỉ giải quyết triệt để bài toán kỹ thuật đo vẽ hiện trường với độ chính xác cao ($M_p \le 6,35\text{ cm}$, vượt trội so với hạn mức cho phép của Thông tư 25/2014/TT-BTNMT), mà còn chuẩn hóa toàn bộ quy trình biên tập, đóng vùng, đánh số thửa và xuất hồ sơ cấp GCNQSDĐ tự động. Đây là mô hình chuẩn mực có khả năng nhân rộng cao cho toàn bộ các địa phương đang triển khai công tác dồn điền đổi thửa và hiện đại hóa hệ thống hồ sơ địa chính trên cả nước.