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 là nền tảng cốt lõi trong quản lý nhà nước về tài nguyên và môi trường. Theo định hướng hiện đại hóa hệ thống hồ sơ địa chính quốc gia và Thông tư số 25/2014/TT-BTNMT của Bộ Tài nguyên và Môi trường, việc ứng dụng công nghệ định vị vệ tinh GNSS kết hợp các phần mềm chuyên ngành xử lý không gian (MicroStation V8i, gCADAS Pro, DPSURVEY) là giải pháp bắt buộc nhằm nâng cao tính chuẩn xác, tinh giảm nhân lực và số hóa 100% dữ liệu đất đai.

Thực tế tại khu vực miền núi như xã Hà Hiệu, huyện Ba Bể, tỉnh Bắc Kạn (diện tích tự nhiên $39{,}98\text{ km}^2$, địa hình chia cắt phức tạp bởi các dãy núi đá vôi và đồi feralit độ dốc cao), tài liệu bản đồ địa chính cũ thành lập từ năm 2006 theo các dự án 3PAD và CSSP đã bộc lộ nhiều sai lệch, biến động ranh giới do quá trình chia tách, chuyển đổi mục đích sử dụng đất nông - lâm nghiệp của 12 thôn bản.

Vấn đề thực tiễn và điểm nghẽn kỹ thuật (Problem Statement)

  • Địa hình hiểm trở và độ che phủ thực vật lớn: Gây suy giảm tín hiệu vệ tinh GNSS, tăng chỉ số PDOP (Position Dilution of Precision) và cản trở thông hướng khi đo chuyền quang học truyền thống.
  • Biến động ranh giới phức tạp: Tình trạng xen cài giữa đất thổ cư, đất lúa nước và đất lâm nghiệp (trồng cây Mỡ) đòi hỏi độ chính xác tiếp biên cao.
  • Rời rạc giữa số liệu ngoại nghiệp và biên tập nội nghiệp: Dữ liệu đo đạc thực địa định dạng thô từ sổ tay điện tử dễ phát sinh lỗi hở vùng (undershoot), bắt vượt (overshoot) khi thành lập mô hình không gian topology.

Mục tiêu dự án cụ thể

  1. Đánh giá toàn diện điều kiện tự nhiên, kinh tế - xã hội và hiện trạng quản lý đất đai tại xã Hà Hiệu, huyện Ba Bể.
  2. Thiết kế và thành lập mạng lưới khống chế đo vẽ gồm 04 điểm địa chính cơ sở và 20 điểm đo vẽ GNSS đạt tiêu chuẩn kỹ thuật cấp 1, cấp 2.
  3. Thực hiện đo đạc chi tiết toàn bộ các góc thửa, ranh giới sử dụng đất thuộc tờ bản đồ số 54.
  4. Xử lý bình sai mạng lưới bằng phần mềm DPSURVEY, biên tập bản đồ địa chính tỷ lệ 1/1000 và chuẩn hóa quan hệ không gian topology bằng gCADAS Pro chạy trên nền Bentley MicroStation V8i.
  5. Lập hồ sơ kỹ thuật thửa đất, thống kê diện tích đất đai và xuất bản bản đồ số phục vụ cấp Giấy chứng nhận quyền sử dụng đất.

Giải pháp kỹ thuật và kết quả kỳ vọng

  • Giải pháp: Kết hợp công nghệ định vị vệ tinh GNSS tĩnh (xây dựng lưới khống chế) và GNSS RTK kết hợp toàn đạc điện tử (đo chi tiết), tự động hóa xử lý nội nghiệp bằng chuỗi công cụ DPSURVEY $\rightarrow$ MicroStation V8i $\rightarrow$ gCADAS Pro.
  • Chỉ tiêu định lượng: Sai số trung phương vị trí điểm sau bình sai $M_p \le 0{,}05\text{ m}$; sai số tương đối cạnh $\le 1/50.000$; sai số khép góc $f_\beta \le 2m_\beta\sqrt{n}$; 100% thửa đất khép vùng kín không lỗi topology; thời gian xử lý dữ liệu giảm trên 60% so với phương pháp thủ công.
  • Phạm vi áp dụng: Mảnh bản đồ địa chính tờ số 54 tỷ lệ 1/1000, xã Hà Hiệu, huyện Ba Bể, tỉnh Bắc Kạn (quy mô khung chuẩn $50\text{ cm} \times 50\text{ cm}$, tương đương diện tích thực địa $25\text{ ha}$).

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

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

Ma trận ưu tiên yêu cầu kỹ thuật (MoSCoW)

  • Must-have (Bắt buộc): Chuyển đổi chính xác hệ tọa độ quốc gia VN-2000 (kinh tuyến trục $106^\circ00'$, múi chiếu $3^\circ$); sai số vị trí điểm khống chế $\le 5\text{ cm}$; 100% thửa đất có thuộc tính: Số tờ, số thửa, diện tích, loại đất, tên chủ sở hữu.
  • Should-have (Cần thiết): Xử lý lọc điểm tự động từ tệp .csv/.crv; gán nhãn quy chủ tự động từ cơ sở dữ liệu bảng tính.
  • Could-have (Có thể mở rộng): Tích hợp thông tin không gian lên nền WebGIS quản lý đất đai cấp xã.
  • Won't-have (Không thực hiện): Đo đạc các đối tượng ngầm (đường dây, ống dẫn nước sâu) ngoài phạm vi quản lý địa chính bề mặt.

Thiết kế hệ thống

graph TD
    A["Vệ tinh GNSS (GPS/GLONASS)"] --> B["Trạm Base / Rover RTK"]
    B --> C["Sổ tay điện tử (.CSV / .CRV)"]
    C --> D["DPSURVEY v3.1: Bình sai lưới"]
    D --> E["MicroStation V8i: Nền tảng Đồ họa"]
    E --> F["gCADAS Pro v2.0: Xử lý Chuyên ngành"]
    F --> G["Tạo Topology & Sửa lỗi Hở/Trùng"]
    F --> H["Liên kết CSDL Thuộc tính (MS Access/Oracle)"]
    G --> I["Bản đồ Địa chính Tờ số 54 (DGN)"]
    H --> I
    I --> J["Hồ sơ kỹ thuật thửa đất & Sổ mục kê"]

Công nghệ và thông số kỹ thuật (Technology Stack)

  • Hệ thống phần cứng: Máy thu GNSS 2 tần số (Trimble/South Galaxy), Máy toàn đạc điện tử Topcon GTS-235W, sổ tay điện tử Field Controller.
  • Nền tảng phần mềm:
    • Bentley MicroStation V8i (SELECTseries 3/4): Hệ thống CAD hạt nhân quản lý các lớp đồ họa (Levels từ 1 đến 63).
    • gCADAS Pro v2.0: Module chuyên ngành địa chính hỗ trợ quản lý cấu trúc thửa đất, tự động tạo vùng topology, gán số thửa và biên tập bản đồ.
    • DPSURVEY v3.1: Module tính toán trắc địa, bình sai lưới đường chuyền mặt bằng và lưới độ cao.
    • Microsoft Excel / Text Editor: Chuẩn hóa dữ liệu trung gian dạng file bảng phân cách (.txt, .csv, .crv).

Cơ sở toán học và mô hình dữ liệu (Database Design)

  • Hệ quy chiếu: VN-2000, Elipxoid WGS-84 quy ước quốc tế hiệu chỉnh:
    • Bán trục lớn: $a = 6.378.137{,}0\text{ m}$
    • Độ dẹt: $\alpha = 1 / 298{,}25723563$
    • Vận tốc góc tự quay: $\omega = 7292115{,}0 \times 10^{-11}\text{ rad/s}$
    • Hằng số trọng trường: $GM = 3986005 \times 10^8\text{ m}^3/\text{s}^2$
    • Điểm gốc tọa độ: $N_{00}$ đặt tại Viện Công nghệ Địa chính, Hoàng Quốc Việt, Hà Nội.
    • Kinh tuyến trục tỉnh Bắc Kạn: $106^\circ00'$, múi chiếu $3^\circ$ (hệ số biến dạng chiều dài $k_0 = 0{,}9999$).
Trường dữ liệu Kiểu dữ liệu Kích thước Mô tả
ID_THUA BigInteger 8 bytes Khóa chính định danh duy nhất thửa đất
SO_TO Integer 4 bytes Số hiệu tờ bản đồ địa chính (Tờ số 54)
SO_THUA Integer 4 bytes Số thứ tự thửa đất trong tờ bản đồ
DIEN_TICH Double Precision 8 bytes Diện tích tính toán hình học ($m^2$, làm tròn 1 chữ số thập phân)
LOAI_DAT Varchar 10 Mã mục đích sử dụng đất (ONT, CLN, LUC, RSX, HNK, DGT)
TEN_CHU NVarchar 150 Tên chủ hộ hoặc tổ chức sử dụng đất
DIA_CHI NVarchar 255 Thôn, bản thuộc xã Hà Hiệu
GEOMETRY Polygon Variable Cấu trúc không gian vùng ranh khép kín

Implementation và kết quả

Quy trình triển khai kỹ thuật (Development Process)

Giai đoạn 1: Ngoại nghiệp (Khảo sát, chôn mốc, đo lưới tĩnh GNSS & đo chi tiết RTK)
Giai đoạn 2: Xử lý dữ liệu thô (Trút sổ tay -> Chuyển đổi tọa độ -> Bình sai DPSURVEY)
Giai đoạn 3: Biên tập không gian (Trút điểm MicroStation -> Nối ranh -> Tạo Topology gCADAS)
Giai đoạn 4: Hoàn thiện hồ sơ (Gán thuộc tính quy chủ -> Xuất bản bản đồ tờ số 54)

Thuật toán tính tọa độ cực và bình sai lưới

Khi sử dụng máy toàn đạc kết hợp trạm máy tĩnh định vị từ điểm gốc $A_1(X_{A1}, Y_{A1})$ định hướng tới điểm kiểm tra $P$, tọa độ phẳng $(X_P, Y_P)$ được xác định theo mô hình:

$$\Delta X_{A1-P} = S \cdot \cos \alpha_{A1-P}$$

$$\Delta Y_{A1-P} = S \cdot \sin \alpha_{A1-P}$$

$$X_P = X_{A1} + \Delta X_{A1-P}, \quad Y_P = Y_{A1} + \Delta Y_{A1-P}$$

Trong đó $S$ là khoảng cách ngang thực tế và $\alpha_{A1-P}$ là góc phương vị cạnh đo.

Sai số khép góc $f_\beta$ của tuyến đường chuyền $n$ đỉnh được kiểm tra theo chuẩn:

$$f_\beta = \left| \sum \beta_{thuc} - \sum \beta_{ly_thuyet} \right| \le 2m_\beta \sqrt{n}$$

(với $m_\beta = \pm 5''$ là sai số trung phương đo góc thiết bị).

Xử lý dữ liệu chuyển đổi định dạng

Dữ liệu thô từ máy RTK được trích xuất qua cổng USB dưới dạng chuỗi byte, lưu trữ vào tệp cấu trúc .crv/.csv, sau đó chuẩn hóa thành cấu trúc không gian qua script tiền xử lý:

# Script chuan hoa du lieu toa do tho tu may do GNSS sang dinh dang trut MicroStation
import csv

def parse_survey_data(input_file, output_dgn_txt):
    """
    Chuyen doi file so lieu do tu GNSS (.crv/.csv) sang tap diem chuan:
    Point_ID, Northing (X), Easting (Y), Height (Z), Code
    """
    with open(input_file, mode='r', encoding='utf-8') as infile, \
         open(output_dgn_txt, mode='w', encoding='utf-8') as outfile:
        reader = csv.reader(infile)
        for row in reader:
            if not row or len(row) < 4:
                continue
            # Trich xuat 4 cot co ban: PointName, Northing, Easting, Height
            point_id = row[0].strip()
            northing = float(row[1].strip()) # Toa do X
            easting  = float(row[2].strip()) # Toa do Y
            height   = float(row[3].strip()) # Do cao H
            code     = row[4].strip() if len(row) > 4 else "RanhDat"
            
            # Ghi ra dinh dang chuan dau vao cho modul trut diem gCADAS
            outfile.write(f"{point_id} {northing:.4f} {easting:.4f} {height:.4f} {code}\n")

if __name__ == "__main__":
    parse_survey_data("06-05-2022.crv", "SoLieuTrutDiem_To54.txt")

Quy trình biên tập bản đồ trên gCADAS Pro & MicroStation V8i

  1. Thiết lập môi trường làm việc: Khởi tạo tệp .dgn theo chuẩn hạt nhân CAD chuẩn thiết kế địa chính; gán kinh tuyến trục $106^\circ00'$, chọn đơn vị hành chính: Tỉnh Bắc Kạn $\rightarrow$ Huyện Ba Bể $\rightarrow$ Xã Hà Hiệu.
  2. Trút điểm chi tiết: Gọi lệnh nhập dữ liệu tọa độ từ tệp văn bản đã xử lý, hệ thống tự động giải bài toán phân lớp đồ họa: Điểm mốc (Level 1), ranh giới thửa (Level 10), nhà ở và công trình xây dựng (Level 21), thủy hệ (Level 32), giao thông (Level 43).
  3. Thiết lập Topology và đóng vùng thửa đất: Sử dụng thuật toán xây dựng đồ thị phẳng định hướng để tự động tính toán nút giao, cắt lọc các đoạn thừa (overshoot) và kéo dài điểm hở (undershoot) với dung sai ngưỡng $\epsilon \le 0{,}002\text{ m}$.
  4. Gán nhãn quy chủ tự động: Kết nối bảng chỉ dẫn diện tích với cơ sở dữ liệu địa chính (.dbf/.mdb), tự động xuất nhãn thửa gồm 4 thông tin cốt lõi: Số hiệu thửa / Diện tích / Mã loại đất / Tên chủ sở hữu.

Kiểm tra, đánh giá và chuẩn hóa dữ liệu (Testing & Validation)

                            DUNG SAI ĐO ĐẠC ĐỊA CHÍNH TỜ 54

Mọi sai số khép cạnh, khép góc đều nằm sâu bên trong giới hạn cho phép của Thông tư 25/2014/TT-BTNMT, bảo đảm tính pháp lý tuyệt đối để nghiệm thu kỹ thuật cấp Sở Tài nguyên và Môi trường.

Kết quả đạt được

Tờ bản đồ địa chính số 54 tỷ lệ 1/1000 xã Hà Hiệu được hoàn thành toàn diện với đầy đủ các lớp đối tượng địa lý và thông tin thuộc tính thửa đất.

Mục đích sử dụng đất Mã đất Số lượng thửa Tổng diện tích ($m^2$) Tỷ lệ cơ cấu (%)
Đất ở tại nông thôn ONT 48 18.240,5 7,30%
Đất trồng lúa nước LUC 112 84.650,2 33,86%
Đất trồng cây lâu năm CLN 35 42.180,0 16,87%
Đất sản xuất lâm nghiệp (Rừng Mỡ) RSX 26 91.540,8 36,62%
Đất giao thông DGT 8 9.320,5 3,73%
Đất sông ngòi, kênh rạch SON 4 4.068,0 1,62%
Tổng cộng 233 250.000,0 100,00%

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

  1. Khắc phục triệt để vùng lõm tín hiệu vệ tinh: Đề xuất quy trình đan kết linh hoạt: Sử dụng GNSS đo tĩnh để lập mạng lưới khống chế cấp cao tại các đỉnh đồi quang đãng $\rightarrow$ Phát triển đường chuyền toàn đạc điện tử hoặc GNSS RTK đo bù tại các khe tụ thủy, tán rừng rậm.
  2. Tự động hóa luồng dữ liệu (Data Pipeline): Loại bỏ hoàn toàn khâu vẽ ranh thủ công; chuyển đổi trực tiếp định dạng từ máy đo sang tệp số liệu chuẩn hóa, tăng năng suất biên tập bản đồ lên 300% (từ trung bình 12 ngày/tờ xuống còn 4 ngày/tờ).
  3. Kiểm soát tính toàn vẹn không gian (Spatial Integrity): Ứng dụng engine topology của gCADAS Pro giúp phát hiện và tự động xử lý các vi hạt thửa (sliver polygons), loại bỏ sai số diện tích do trùng đè ranh giới giữa các hộ dân.
QUY TRÌNH TRUYỀN THỐNG:

QUY TRÌNH ĐỀ TÀI:

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

Tình huống ứng dụng thực địa

  • Cấp mới và đổi Giấy chứng nhận quyền sử dụng đất: Bản đồ tờ 54 cung cấp số liệu kích thước hình học, tọa độ đỉnh thửa chính xác đến từng centimet, loại bỏ hoàn toàn tranh chấp ranh giới giữa các chủ hộ.
  • Phục vụ bồi thường, giải phóng mặt bằng: Phục vụ công tác quy hoạch mở rộng nâng cấp Quốc lộ 279 và Tỉnh lộ 212 đoạn đi qua xã Hà Hiệu.
  • Thành lập hồ sơ địa chính số: Đồng bộ trực tiếp vào hệ thống cơ sở dữ liệu VILIS/VBDLIS của Chi nhánh Văn phòng Đăng ký Đất đai huyện Ba Bể.

Phân tích hiệu quả kinh tế (ROI & Chi phí)

  • Tiết kiệm nhân công: Tổ đo giảm từ 5 kỹ thuật viên xuống còn 2 kỹ thuật viên cho một ca máy.
  • Rút ngắn thời gian khảo sát ngoại nghiệp: Giảm 65% số ngày bám tuyến thực địa tại vùng núi.
  • Thời gian hoàn vốn thiết bị (ROI): Tối ưu hóa chu kỳ luân chuyển máy móc, đạt điểm hoàn vốn đầu tư công nghệ trong vòng 1,2 năm triển khai các gói thầu đo đạc địa chính.

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

Hạn chế kỹ thuật

  • Tín hiệu GNSS RTK vẫn bị suy giảm cục bộ (mất Fix chuyển sang chế độ Float) tại những vị trí hẻm núi sâu có vách đá vôi dựng đứng hoặc dưới tán cây lâm nghiệp rậm rạp.
  • Công tác xác minh quy chủ đất đai vẫn phụ thuộc vào sự có mặt và hợp tác của người dân địa phương tại thời điểm đo vẽ thực địa.

Hướng phát triển

  • Tích hợp công nghệ bay chụp ảnh hàng không bằng máy bay không người lái (UAV/Drone trang bị cảm biến LiDAR) để quét bình đồ địa chính khu vực đồi núi phức tạp.
  • Phát triển module WebGIS ứng dụng mã nguồn mở (GeoServer, MapLibre GL) phục vụ công khai quy hoạch và tra cứu thông tin thửa đất trực tuyến cho người dân.

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

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

1. Yêu cầu cấu hình phần cứng và phần mềm để xử lý tệp bản đồ DGN trên MicroStation V8i và gCADAS Pro?

Hệ thống khuyến nghị bao gồm: CPU Intel Core i5 thế hệ 8 trở lên (hoặc tương đương), RAM tối thiểu 8GB (khuyến nghị 16GB), ổ cứng SSD dung lượng trống tối thiểu 20GB, Card đồ họa rời tối thiểu 2GB VRAM, hệ điều hành Windows 10/11 64-bit có cài đặt sẵn nền tảng Microsoft .NET Framework 3.5/4.8 và khóa cứng bản quyền gCADAS Pro.

2. Khi đo vẽ chi tiết bằng GNSS RTK ở khu vực đồi núi bị mất sóng Fix thì xử lý như thế nào?

Khi bộ thu chuyển sang trạng thái Float (độ chính xác giảm xuống hàng mét), kỹ thuật viên chuyển sang phương pháp đo bù: Sử dụng máy toàn đạc điện tử đặt trạm tại điểm đã Fix tọa độ gần nhất, hoặc thiết lập cọc phụ thông hướng để đo giao hội/tọa độ cực đến điểm ranh giới thửa đất bị che khuất.

3. Tại sao phải chuyển đổi số liệu qua tệp trung gian đuôi .crv/.csv trước khi trút vào gCADAS?

Sổ tay điều khiển GNSS lưu trữ nhiều thông tin phụ trợ (PDOP, số lượng vệ tinh, thời gian đo, trạng thái giải mã). Việc đưa qua định dạng .crv/.csv giúp trích lọc chính xác 4 trường dữ liệu chuẩn [Point_ID, X, Y, H] tương thích hoàn toàn với cấu trúc bảng thuộc tính đồ họa của MicroStation V8i.

4. gCADAS Pro xử lý các lỗi hình học không gian (Topology Errors) bằng cơ chế nào?

Phần mềm sử dụng thuật toán quét không gian (spatial sweep algorithm) với tham số dung sai tự định nghĩa (Tolerance $\approx 0{,}002\text{ m}$). Hệ thống tự động phát hiện các đoạn đường giao nhau không có đỉnh nút, điểm bắt hở ranh (undershoot) và tự động kéo dài nối khép kín để đóng vùng diện tích (Polygon creation).

5. Dự án có thể mở rộng cho toàn bộ các xã khác thuộc huyện Ba Bể được không?

Hoàn toàn khả thi. Do toàn bộ mạng lưới đã được thiết kế thống nhất trên Ellipsoid WGS-84, lưới chiếu UTM/Gauss-Kruger kinh tuyến trục $106^\circ00'$ tỉnh Bắc Kạn, mô hình này có thể nhân rộng đồng bộ cho toàn bộ 11 xã, thị trấn trên địa bàn huyện nhằm xây dựng hệ thống cơ sở dữ liệu đất đai đa mục tiêu đồng nhất.

Kết luận

Đề tài tốt nghiệp "Ứng dụng công nghệ GPS, tin học và các phần mềm chuyên ngành trong công tác thành lập mảnh bản đồ địa chính tờ số 54 tỷ lệ 1/1000 xã Hà Hiệu, huyện Ba Bể, tỉnh Bắc Kạn" đã chứng minh tính ưu việt vượt trội của việc số hóa quy trình quản lý đất đai. Việc kết hợp chặt chẽ giữa thiết bị định vị vệ tinh hiện đại (GNSS) và hệ thống phần mềm chuyên dụng (DPSURVEY, MicroStation V8i, gCADAS Pro) không chỉ giải quyết triệt để bài toán đo đạc trên địa hình đồi núi phức tạp mà còn đảm bảo sản phẩm bản đồ đạt độ chính xác tuyệt đối theo Thông tư 25/2014/TT-BTNMT.

Kết quả thành lập tờ bản đồ số 54 với 233 thửa đất khép kín không lỗi topology là minh chứng thực tế cho khả năng tự động hóa và đồng bộ dữ liệu địa chính. Đây là tài liệu khoa học và thực tiễn quan trọng, đóng góp trực tiếp vào công tác quản lý tài nguyên, quy hoạch sử dụng đất và thúc đẩy kinh tế - xã hội tại địa phương phát triển bền vững.