Giới thiệu dự án

Công tác quản lý nhà nước về đất đai đòi hỏi hệ thống tư liệu không gian và thuộc tính đạt độ chính xác cao, đồng bộ và có giá trị pháp lý tuyệt đối. Theo số liệu thống kê từ Bộ Tài nguyên và Môi trường (Bộ TN&MT), việc hoàn thiện cơ sở dữ liệu địa chính số giúp giảm hơn 60% thời gian xử lý thủ tục hành chính về đất đai và hạn chế tối đa tình trạng tranh chấp ranh giới. Tại địa bàn xã Tứ Quận, huyện Yên Sơn, tỉnh Tuyên Quang — khu vực có địa hình trung du miền núi phức tạp với sự đan xen giữa đất ở, đất vườn đồi và đất trồng cây lâu năm — việc sử dụng các phương pháp đo đạc thủ công truyền thống bộc lộ nhiều hạn chế về thời gian, nhân lực và sai số lũy tích.

Vấn đề cốt lõi đặt ra là sự thiếu hụt bản đồ địa chính chính quy tỷ lệ lớn (1:1000) được số hóa hoàn toàn theo Hệ quy chiếu và Hệ tọa độ quốc gia VN-2000, gây khó khăn cho công tác đăng ký đất đai, cấp Giấy chứng nhận quyền sử dụng đất và tích hợp cơ sở dữ liệu vào hệ sinh thái quản lý địa chính số.

Mục tiêu dự án

  1. Thiết lập mạng lưới khống chế đo vẽ: Xây dựng và bình sai lưới đường chuyền kinh vĩ cấp 1, cấp 2 gồm 68 điểm khống chế (kết hợp 3 điểm địa chính cấp cao làm điểm khởi tính) đảm bảo sai số vị trí điểm sau bình sai $M_p \le 5\text{ cm}$.
  2. Đo vẽ chi tiết thực địa: Ứng dụng phương pháp tọa độ cực thông qua máy toàn đạc điện tử Topcon GTS-236 để thu thập toàn bộ các yếu tố ranh giới thửa đất, địa vật và công trình nhân tạo.
  3. Số hóa và biên tập tờ bản đồ địa chính số 11 tỷ lệ 1:1000: Sử dụng tích hợp bộ giải pháp phần mềm MicroStation V8i, FAMIS và Gcadas để làm sạch dữ liệu hình học (topology), tạo tâm thửa, đánh số thửa và gán nhãn thuộc tính theo Thông tư 25/2014/TT-BTNMT.

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

  • Phương pháp: Kết hợp đo ngắm quang - điện tử ngoại nghiệp độ chính xác cao với công nghệ xử lý đồ họa vector, phân lớp dữ liệu tự động nội nghiệp.
  • Phạm vi: Tờ bản đồ địa chính số 11, tỷ lệ 1:1000 (kích thước khung chuẩn $50 \times 50\text{ cm}$, tương ứng diện tích thực địa $25\text{ ha}$) tại xã Tứ Quận, huyện Yên Sơn, tỉnh Tuyên Quang; kinh tuyến trục quy định $106^\circ00'$, múi chiếu $3^\circ$.

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

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

Tiêu chí kỹ thuật Phương pháp Kinh vĩ quang cơ & Thước thép Phương pháp Bay chụp UAV / Ảnh hàng không Giải pháp Toàn đạc điện tử + CAD/GIS
Độ chính xác mặt bằng Thấp ($M_p > 15\text{ cm}$ do sai số đọc số và đo dài) Trung bình ($5 - 10\text{ cm}$, phụ thuộc điểm GCP) Cao ($M_p \le 5\text{ cm}$, đáp ứng tỷ lệ 1:1000)
Khả năng quan sát địa vật Hạn chế, thao tác thủ công nặng nhọc Bị che khuất mạnh bởi tán cây rậm rạp Linh hoạt, đo trực tiếp góc/cạnh tại góc thửa
Tốc độ số hóa dữ liệu Rất chậm (nhập liệu tay, dễ sai sót) Nhanh nhưng khâu lọc đám mây điểm phức tạp Tự động hóa qua file ASCII/SDR vào CAD
Chi phí triển khai Thấp về thiết bị, cao về nhân công Cao về thiết bị và bản quyền phần mềm xử lý Cân bằng, hiệu quả kinh tế và độ tin cậy cao

Yêu cầu người dùng theo khung MoSCoW

  • Must have: Lưới đường chuyền khống chế đạt sai số khép góc $f_\beta \le 2m_\beta\sqrt{n}$; bản đồ thể hiện đủ 3 yếu tố (số thửa, diện tích, loại đất); phân lớp chuẩn level 1-63.
  • Should have: Tự động sửa lỗi tiếp biên và đóng vùng thửa đất qua thuật toán MRFClean; gán nhãn động từ cơ sở dữ liệu thuộc tính.
  • Could have: Xuất dữ liệu đồng bộ sang định dạng ViLIS, eLIS phục vụ tích hợp CSDL cấp huyện.
  • Won't have (trong phạm vi này): Xây dựng mô hình thông tin công trình địa chính 3D (3D Cadastre).

Thiết kế hệ thống và Cơ sở toán học

Cơ sở toán học của bản đồ được xác lập trên Hệ quy chiếu và Hệ tọa độ quốc gia VN-2000 với các tham số chuẩn:

  • Ellipsoid quy chiếu: WGS-84 toàn cầu hiệu chỉnh cho Việt Nam ($a = 6.378.137\text{ m}$; độ dẹt $\alpha = 1/298,25723563$).
  • Phép chiếu: Lưới chiếu hình trụ ngang đồng góc UTM, múi chiếu $3^\circ$, hệ số co dãn chiều dài tại kinh tuyến trục $k_0 = 0,9999$.
  • Kinh tuyến trục tỉnh Tuyên Quang: $106^\circ00'$.

Technology Stack và Thông số phiên bản

  • Phần cứng thu thập dữ liệu: Máy toàn đạc điện tử Topcon GTS-236 (độ chính xác đo góc $2"$, đo cạnh $\pm(2\text{mm} + 2\text{ppm}\cdot D)$); Hệ thống định vị GNSS GPS South 6 tần số.
  • Môi trường đồ họa nền: Bentley MicroStation V8i (v08.11.09.357), định dạng tệp chuẩn .dgn.
  • Phần mềm chuyên ngành địa chính: FAMIS (Field Work and Cadastral Mapping Integrated Software) Ver 2006; Module tiện ích Gcadas 2014.
  • Hệ quản trị CSDL thuộc tính: Microsoft Access / dBase liên kết qua ODBC.

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

Dự án áp dụng quy trình kiểm soát chất lượng 2 tầng (Ngoại nghiệp $\leftrightarrow$ Nội nghiệp) theo mô hình tuần hoàn khép kín:

Implementation và kết quả

Quá trình xử lý và Thuật toán cốt lõi

1. Thuật toán xác định tọa độ phẳng trạm đo chi tiết (Phương pháp tọa độ cực)

Số liệu đo từ Topcon GTS-236 bao gồm: góc bằng $\beta$, góc đứng $V$ (hoặc góc thiên đỉnh $Z$), và khoảng cách nghiêng $D$. Bộ xử lý tính toán chuyển đổi cạnh ngang và tọa độ điểm chi tiết $P(X_P, Y_P, H_P)$ từ trạm máy $A(X_A, Y_A, H_A)$ ngắm định hướng về $B(X_B, Y_B)$:

$$\begin{cases} S_{AP} = D_{AP} \cdot \cos V_P = D_{AP} \cdot \sin Z_P \ \alpha_{AP} = \alpha_{AB} + \beta_P \ \Delta X_{AP} = S_{AP} \cdot \cos \alpha_{AP} \ \Delta Y_{AP} = S_{AP} \cdot \sin \alpha_{AP} \ X_P = X_A + \Delta X_{AP} \ Y_P = Y_A + \Delta Y_{AP} \ H_P = H_A + S_{AP} \cdot \cot Z_P + i_m - l_g \end{cases}$$

Trong đó: $\alpha_{AB}$ là góc phương vị định hướng; $i_m$ là chiều cao máy; $l_g$ là chiều cao gương sào.

import math

def calculate_polar_point(station_x, station_y, station_h, azimuth_ab, beta, slope_dist, zenith_angle, inst_h, prism_h):
    """
    Tính tọa độ điểm chi tiết theo phương pháp tọa độ cực
    """
    # Chuyển đổi độ sang radian
    rad_zenith = math.radians(zenith_angle)
    rad_azimuth = math.radians(azimuth_ab + beta)
    
    # Tính khoảng cách ngang và chênh cao
    horizontal_dist = slope_dist * math.sin(rad_zenith)
    delta_x = horizontal_dist * math.cos(rad_azimuth)
    delta_y = horizontal_dist * math.sin(rad_azimuth)
    delta_h = (horizontal_dist / math.tan(rad_zenith)) + inst_h - prism_h
    
    point_x = station_x + delta_x
    point_y = station_y + delta_y
    point_h = station_h + delta_h
    
    return round(point_x, 3), round(point_y, 3), round(point_h, 3)

# Ví dụ thực nghiệm trạm đo A01 ngắm định hướng A02
x_p, y_p, h_p = calculate_polar_point(
    station_x=2415620.125, station_y=564310.850, station_h=45.20,
    azimuth_ab=48.5211, beta=15.2340, slope_dist=124.56,
    zenith_angle=88.7520, inst_h=1.45, prism_h=1.60
)

2. Xử lý Topology và Đóng vùng thửa đất trong FAMIS/MicroStation

Dữ liệu đo sau khi trút từ thiết bị ở định dạng ASCII/Raw Data được chuẩn hóa thành tệp .dat với cấu trúc:

Point_ID, X_Coordinate, Y_Coordinate, Elevation, Feature_Code
101, 2415654.321, 564345.112, 45.32, RANH_THUA
102, 2415680.145, 564360.874, 45.40, RANH_THUA
103, 2415672.890, 564395.230, 45.85, NHA_G

Quy trình tự động hóa Topology:

  1. MRFClean (Node Snapping & Undershoot/Overshoot Trimming): Xử lý triệt để sai số hở đỉnh ($< 0,02\text{ m}$) và bắt dính các đường biên thửa chồng lấn.
  2. Centroid Creation: Khởi tạo tâm thửa tự động, xác định quan hệ lân cận không gian giữa các đối tượng thửa đất.
  3. Attribute Association: Tự động gán nhãn gồm 3 tham số cơ bản: Số thứ tự thửa, Diện tích ($m^2$), và Mã loại đất (ONT, CLN, HNK, LUC...).

Testing và Validation

Đánh giá độ chính xác lưới khống chế đường chuyền

Chỉ tiêu kỹ thuật kiểm định Quy chuẩn TT 25/2014/TT-BTNMT Kết quả thực tế tại khu đo Đánh giá
Sai số khép góc ($f_\beta$) $\le \pm 5" \sqrt{n} = \pm 18,7"$ $+8,5"$ Đạt chuẩn ($< 50%$ giới hạn)
Sai số khép tương đối ($f_s / [S]$) $\le 1/25.000$ (Cấp 1); $\le 1/10.000$ (Cấp 2) $1/32.400$ Vượt yêu cầu độ chính xác
Sai số vị trí điểm sau bình sai ($M_p$) $\le 5,0\text{ cm}$ $2,1\text{ cm} - 3,4\text{ cm}$ Đảm bảo độ tin cậy tuyệt đối
Sai số trung phương đo cạnh ($M_s$) $\le 0,05\text{ m}$ $0,018\text{ m}$ Đạt tiêu chuẩn kỹ thuật
Hạn sai đo góc (Dao động 2C) $\le 12"$ $\le 6"$ Ổn định quang học cao

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

  • Sản phẩm số hóa: Hoàn thiện $100%$ tờ bản đồ địa chính số 11 tỷ lệ 1:1000 xã Tứ Quận, huyện Yên Sơn với 184 thửa đất khép vùng Topology chuẩn, không xảy ra lỗi trùng đè hoặc sót thửa.
  • Bộ hồ sơ kỹ thuật: Xuất tự động toàn bộ trích đo địa chính thửa đất, bảng tổng hợp diện tích theo chủ sử dụng và sổ mục kê đất đai đồng bộ với cơ sở dữ liệu địa chính cấp xã.

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

Cải tiến công nghệ và Hiệu quả định lượng

Tiêu chí so sánh Quy trình truyền thống (Kinh vĩ + CAD thuần) Mô hình FAMIS + MicroStation V8i Cải tiến tích hợp (Topcon + FAMIS + Gcadas)
Thời gian đo ngoại nghiệp 100% thời gian cơ sở Giảm 25% nhờ trút dữ liệu Giảm 45% nhờ tự động hóa mã trạm đo
Khâu sửa lỗi hình học Thủ công từng nét vẽ (rất lâu) Bán tự động (dễ sót lỗi) Tự động 98% thông qua MRFClean
Độ chính xác tọa độ điểm Phụ thuộc người vẽ ($\pm 0,2\text{ mm}$ bản đồ) Chính xác số hóa ($\pm 0,001\text{ m}$) Tuyệt đối theo trị đo số ($\pm 0,001\text{ m}$)
Khả năng kết xuất CSDL Không hỗ trợ xuất tự động Xuất bảng biểu cơ bản Đồng bộ hóa 100% chuẩn ViLIS/eLIS
  • Đóng góp chuyên ngành: Chuẩn hóa quy trình liên kết liền mạch từ thiết bị đo dã ngoại thế hệ mới vào hệ sinh thái phần mềm MicroStation, thiết lập tài liệu kỹ thuật mẫu cho các đề tài quy hoạch và quản lý đất đai khu vực miền núi phía Bắc.

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

Kịch bản ứng dụng trong quản lý nhà nước

  • Cấp mới và cấp đổi Giấy chứng nhận QSDĐ: Cung cấp đường ranh giới và diện tích thửa đất chính xác tuyệt đối, giảm thiểu tranh chấp đất đai giáp ranh giữa các hộ gia đình.
  • Thu hồi, bồi thường giải phóng mặt bằng: Trích xuất dữ liệu ranh giới địa chính phục vụ tính toán áp giá đền bù cho các dự án giao thông, thủy lợi trên địa bàn huyện Yên Sơn.
  • Tích hợp WebGIS và Cơ sở dữ liệu quốc gia: Tệp dữ liệu vector .dgn và bảng thuộc tính sẵn sàng chuyển đổi sang định dạng không gian GeoDatabase (.gdb, .shp) phục vụ kết nối CSDL đất đai quốc gia.

Yêu cầu triển khai hệ thống

  • Yêu cầu phần cứng: Workstation Intel Core i5 thế hệ 8 trở lên, RAM tối thiểu 8GB, màn hình độ phân giải Full HD để thao tác bản đồ vector.
  • Yêu cầu phần mềm: MicroStation V8i (SelectSeries), bộ gõ TCVN3 cho hệ font địa chính (VnTime, VnArial), phần mềm chuyên ngành FAMIS và Gcadas 2014.

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

Hạn chế kỹ thuật

  • Tầm nhìn quang học bị giới hạn: Tại các khu vực đồi chè và vườn cây ăn quả rậm rạp, máy toàn đạc điện tử đòi hỏi phải phát quang thông hướng ngắm, làm tăng thời gian đo ngoại nghiệp.
  • Ràng buộc tương thích nền tảng: Phần mềm FAMIS được phát triển trên kiến trúc đồ họa 32-bit cũ, đôi khi gặp xung đột font chữ trên các hệ điều hành 64-bit hiện đại.

Hướng phát triển

  1. Ứng dụng công nghệ bay quét LiDAR / UAV RTK: Kết hợp máy bay không người lái để thu thập mô hình số độ cao (DEM) và bình đồ ảnh trực chính (Orthophoto) tại các vùng địa hình đồi núi hiểm trở, giảm 80% công tác phát quang thực địa.
  2. Nâng cấp sang giải pháp OpenGIS / QGIS Cadastre: Tái cấu trúc cơ sở dữ liệu sang định dạng không gian mở (PostGIS, GeoServer) hỗ trợ chia sẻ dữ liệu trực tuyến theo chuẩn OGC.

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

  • Sinh viên & Học viên chuyên ngành Quản lý đất đai / Trắc địa: Nắm vững quy trình công nghệ kết hợp giữa đo đạc trắc địa truyền thống với các phần mềm chuyên ngành hiện đại.
  • Kỹ sư Trắc địa & Cán bộ địa chính cấp xã: Sở hữu quy trình thực hành chi tiết về xử lý số liệu đo toàn đạc, bình sai lưới và biên tập bản đồ đạt chuẩn kỹ thuật của Bộ TN&MT.
  • Cơ quan Quản lý Nhà nước (UBND xã Tứ Quận, Phòng TN&MT huyện Yên Sơn): Tiếp nhận bộ bản đồ địa chính tờ số 11 tỷ lệ 1:1000 số hóa đồng bộ, có tính pháp lý cao để phục vụ công tác quản lý chuyên môn.
  • Chủ sử dụng đất tại địa phương: Được xác lập quyền sử dụng đất rõ ràng trên bản đồ chính quy, hạn chế tối đa rủi ro tranh chấp.

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

1. Yêu cầu kỹ thuật tối thiểu để vận hành phần mềm FAMIS và Gcadas trên máy tính hiện nay là gì?

Hệ thống yêu cầu máy tính chạy hệ điều hành Windows 7/10 (hỗ trợ tương thích ứng dụng 32-bit), cài đặt bộ font TCVN3 (.VnTime), thư viện Visual C++ Redistributable và môi trường CAD nền tảng Bentley MicroStation V8i hoặc SE.

2. Giới hạn sai số cho phép khi đo vẽ bản đồ địa chính tỷ lệ 1:1000 vùng nông thôn là bao nhiêu?

Theo Thông tư 25/2014/TT-BTNMT, sai số trung phương vị trí điểm ranh giới thửa đất biểu thị trên bản đồ địa chính tỷ lệ 1:1000 so với điểm khống chế đo vẽ gần nhất không được vượt quá $0,1\text{ m}$ ($10\text{ cm}$) đối với khu vực đất nông thôn đồng bằng và $0,15\text{ m}$ ($15\text{ cm}$) đối với vùng đồi núi.

3. Dữ liệu bản đồ trên MicroStation V8i có tích hợp trực tiếp vào hệ thống ViLIS/eLIS được không?

Có. Module Gcadas tích hợp sẵn chức năng trích xuất thuộc tính không gian và phi không gian sang định dạng cơ sở dữ liệu trung gian (.xml, .mdb, .shp), giúp nạp trực tiếp vào hệ thống thông tin đất đai ViLIS và eLIS mà không làm biến dạng dữ liệu.

4. Cần thực hiện kiểm định máy toàn đạc điện tử Topcon GTS-236 theo chu kỳ như thế nào?

Máy toàn đạc điện tử phải được kiểm định sai số góc dọi tâm quang học, sai số 2C, góc bù tự động và hằng số gương ($K$) định kỳ 6 tháng/lần tại các phòng kiểm định trắc địa đạt chuẩn trước khi đưa vào thực địa.

5. Chi phí đầu tư giải pháp công nghệ toàn đạc + phần mềm chuyên ngành so với bay chụp UAV?

Giải pháp toàn đạc điện tử kết hợp phần mềm MicroStation/Gcadas có chi phí thiết bị ban đầu thấp hơn 40-50% so với tổ hợp UAV tích hợp LiDAR/RTK và phần mềm xử lý ảnh chuyên sâu (Pix4D/Agisoft), rất phù hợp cho các đơn vị đo đạc vừa và nhỏ triển khai theo từng tờ bản đồ độc lập.

Kết luận

Đề tài tốt nghiệp "Ứng dụng công nghệ tin học và máy toàn đạc điện tử thành lập bản đồ địa chính tờ số 11 tỷ lệ 1:1000 xã Tứ Quận, huyện Yên Sơn, tỉnh Tuyên Quang" đã giải quyết triệt để bài toán hiện đại hóa công tác biên tập bản đồ địa chính tại địa phương. Việc tích hợp chuẩn xác giữa máy toàn đạc điện tử Topcon GTS-236, hệ thống định vị GPS South cùng bộ phần mềm MicroStation V8i, FAMIS và Gcadas 2014 đảm bảo toàn bộ quy trình đo đạc, bình sai và xây dựng Topology đạt chuẩn theo Thông tư 25/2014/TT-BTNMT. Kết quả của dự án là nguồn tư liệu pháp lý và kỹ thuật quan trọng, đóng góp thiết thực vào công cuộc số hóa cơ sở dữ liệu đất đai và chuyển đổi số ngành quản lý tài nguyên môi trường.