Giới thiệu dự án

Công tác đo đạc, lập bản đồ địa chính và cấp Giấy chứng nhận quyền sử dụng đất (GCNQSDĐ) là trụ cột cốt lõi trong quản lý nhà nước về đất đai được quy định tại Luật Đất đai năm 2013. Tại các khu vực đô thị phát triển nhanh như thành phố Lạng Sơn, quá trình chuyển đổi mục đích sử dụng đất diễn ra mạnh mẽ: theo số liệu thống kê hiện trạng, đất phi nông nghiệp tại phường Vĩnh Trại chiếm tới 91,03% (147,04 ha trên tổng số 161,53 ha diện tích tự nhiên), với mật độ dân cư tập trung đông đúc (15.355 nhân khẩu trên 4.493 hộ phân bố tại 11 khối phố). Hệ thống tài liệu bản đồ địa chính lập từ năm 1997 đã biến động sâu sắc, xuất hiện nhiều sai lệch về ranh giới thực địa, chia tách thửa đất tự phát và lấn chiếm hành lang an toàn giao thông, gây cản trở lớn cho công tác quản lý hành chính và cấp đổi GCNQSDĐ.

Vấn đề kỹ thuật then chốt đặt ra là các phương pháp đo vẽ thủ công truyền thống (bàn đạc, máy kinh vĩ quang cơ) không đáp ứng được yêu cầu về độ chính xác mặt bằng ($\le 0,1\text{ mm}$ trên tỷ lệ bản đồ tương đương sai số thực địa $\le 5\text{ cm}$ ở tỷ lệ 1:500) và tốc độ số hóa dữ liệu. Nhằm giải quyết triệt để bài toán này, đề tài triển khai giải pháp kỹ thuật tích hợp đồng bộ giữa thiết bị phần cứng ngoại nghiệp (máy toàn đạc điện tử SOUTH NTS-312B) và chuỗi phần mềm tin học chuyên dụng nội nghiệp (DPSurvey, Bentley MicroStation V8i, Gcadas).

Mục tiêu cụ thể của dự án gồm:

  1. Thiết lập mạng lưới khống chế đo vẽ tọa độ và độ cao cấp 2 trên nền hệ quy chiếu quốc gia VN-2000, đảm bảo sai số khép góc $f_\beta \le 5\sqrt{n}$ giây và sai số khép tương đối $f_s/\sum S \le 1/25.000$.
  2. Thu thập dữ liệu trắc địa chi tiết ranh giới thửa đất, địa vật và công trình xây dựng cho toàn bộ tờ bản đồ số 10 tỷ lệ 1:500 tại phường Vĩnh Trại.
  3. Xây dựng mô hình hình học topo (Topology) và cơ sở dữ liệu thuộc tính không gian địa chính số hóa, xuất hồ sơ kỹ thuật thửa đất và bảng diện tích quy chủ chính xác 100%.

Phạm vi nghiên cứu tập trung thực nghiệm trên tờ bản đồ địa chính số 10 tỷ lệ 1:500 (diện tích 6,25 ha) thuộc địa bàn phường Vĩnh Trại, thành phố Lạng Sơn, thực hiện bởi Công ty TNHH MTV Mạnh Chung phối hợp cùng Trường Đại học Nông Lâm Thái Nguyên trong thời gian từ tháng 05/2019 đến tháng 09/2019.

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

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

Tại các đô thị miền núi có mật độ xây dựng cao, công tác đo đạc địa chính thường đối mặt với ba phương pháp tiếp cận chính:

Tiêu chí so sánh Phương pháp Bàn đạc / Quang cơ Phương pháp UAV / Viễn thám Máy toàn đạc điện tử + CAD/GIS
Độ chính xác mặt bằng Thấp ($\pm 15 - 30\text{ cm}$) Trung bình ($\pm 7 - 15\text{ cm}$) do góc khuất nhà cao tầng Rất cao ($\pm 1 - 3\text{ cm}$), đáp ứng tỷ lệ 1:500
Khả năng đo góc khuất ngõ hẻm Kém, phụ thuộc thước kéo thủ công Rất kém (bị che khuất bởi tán cây, mái hiên) Linh hoạt bố trí trạm gương phụ
Mức độ số hóa dữ liệu Thủ công, tốn nhân lực số hóa lại Ảnh nắn trực giao tự động Đồng bộ trực tiếp từ Raw Data sang DGN/Vector
Chi phí thiết bị và vận hành Thấp nhưng chi phí nhân công cao Cao, yêu cầu cấp phép bay phức tạp Tối ưu, thiết bị bền bỉ và phổ biến

Yêu cầu người dùng được phân loại theo mô hình MoSCoW:

  • Must-have (Bắt buộc): Đo đạc tọa độ 3 chiều $(X, Y, H)$ chuẩn VN-2000; lập mạng lưới đường chuyền khống chế đạt tiêu chuẩn Thông tư 25/2014/TT-BTNMT; tự động hóa làm sạch topology và đóng thửa.
  • Should-have (Nên có): Tự động liên kết nhãn thuộc tính (tên chủ hộ, loại đất, diện tích) với bảng biểu Excel theo chuẩn Thông tư 28/2014/TT-BTNMT.
  • Could-have (Có thể có): Kết xuất trực tiếp dữ liệu dạng shapefile sang các hệ thống thông tin đất đai ViLIS 2.0 và TMV.LIS.
  • Won't-have (Chưa thực hiện đợt này): Mô hình hóa dữ liệu địa chính 3D đa tầng (3D Cadastre) cho các khu chung cư cao tầng.

Thiết kế hệ thống

Kiến trúc hệ thống đo vẽ và xử lý dữ liệu địa chính số hóa tuân theo luồng xử lý tuần tự từ ngoại nghiệp đến cấu trúc CSDL không gian:

Technology Stack:

  • Thiết bị ngoại nghiệp: Máy toàn đạc điện tử SOUTH NTS-312B (độ phóng đại 30X, độ chính xác đo góc 2", sai số đo khoảng cách $\pm(2\text{ mm} + 2\text{ ppm} \cdot D)$, bộ nhớ trong 100.000 điểm đo).
  • Phần mềm tính toán & bình sai: DPSurvey v3.2 và Pronet (xử lý bình sai đường chuyền phù hợp và đường chuyền khép kín).
  • Nền tảng CAD/GIS: Bentley MicroStation V8i (v8.11.09) sử dụng Seedfile chuẩn địa chính 2D/3D.
  • Hệ thống chuyên ngành: Gcadas (Geo-Cadastre System) phiên bản đồng bộ chuẩn Thông tư 25/2014/TT-BTNMT.

Cơ sở toán học và cấu trúc dữ liệu: Bản đồ được thiết lập trên Ellipsoid WGS84 cải biên của hệ quy chiếu VN-2000: Bán trục lớn $a = 6.378.137,0\text{ m}$, độ dẹt $\alpha = 1/298,25723563$, kinh tuyến trục Lạng Sơn $107^\circ 15'$, múi chiếu $3^\circ$, hệ số biến dạng tỷ lệ $k_0 = 0,9999$. Quy tắc phân mảnh ô vuông tọa độ tỷ lệ 1:500 có kích thước khung trong tiêu chuẩn $50 \times 50\text{ cm}$, tương ứng kích thước thực địa $250 \times 250\text{ m}$ (diện tích 6,25 ha).

-- Cấu trúc bảng thuộc tính thửa đất địa chính (Schema chuẩn Gcadas)
CREATE TABLE Cadastral_Parcel (
    Parcel_ID VARCHAR(20) PRIMARY KEY,
    Map_Sheet_Number INT NOT NULL,
    Parcel_Number INT NOT NULL,
    Area_Square_Meter DECIMAL(10, 2) NOT NULL,
    Land_Use_Category VARCHAR(10) NOT NULL, -- ODT, CLN, DGT, etc.
    Owner_Name NVARCHAR(100) NOT NULL,
    Address NVARCHAR(255),
    Legal_Status NVARCHAR(50),
    Geometry_Spatial GEOMETRY NOT NULL
);

Methodology

Dự án áp dụng phương pháp luận khảo sát trắc địa địa chính chuẩn hóa (Surveying Lifecycle Workflow):

  1. Thi công khảo sát & Thiết kế mạng lưới (Tuần 1 - 2): Định vị 4 mốc gốc GPS hạng cao xung quanh khu vực phường Vĩnh Trại làm điểm khởi tính, chọn vị trí chôn mốc đường chuyền kinh vĩ cấp 2 tại các vị trí quang đãng, nền đất ổn định.
  2. Đo đạc ngoại nghiệp (Tuần 3 - 8): Đo góc ngang, góc đứng và cạnh đường chuyền theo phương pháp toàn vòng (hoặc 2 nửa lần đo), sai số 2C không vượt quá 12", chênh lệch giữa các lần đo $\le 8"$. Đo chi tiết ranh giới thửa đất theo phương pháp tọa độ cực.
  3. Xử lý nội nghiệp & Bình sai (Tuần 9 - 12): Bình sai gián tiếp lưới trắc địa bằng phần mềm Pronet, trút dữ liệu góc cạnh vào DPSurvey chuyển đổi sang tọa độ Descartes $(X, Y)$.
  4. Biên tập bản đồ & Đóng Topology (Tuần 13 - 16): Xử lý lỗi hở ranh (undershoot), cắt thừa ranh (overshoot), gán mã loại đất và xuất bản đồ nghiệm thu.

Implementation và kết quả

Development process

Quá trình đo vẽ chi tiết và tính toán tọa độ dựa trên nguyên lý giải bài toán trắc địa thuận từ hệ tọa độ cực trạm máy:

$$\Delta X_{AB} = X_B - X_A, \quad \Delta Y_{AB} = Y_B - Y_A$$ $$\alpha_{AB} = \arctan\left(\frac{\Delta Y_{AB}}{\Delta X_{AB}}\right) + k\pi$$ $$\alpha_{A-P_i} = \alpha_{AB} + \beta_i$$ $$S_{A-P_i} = D_{A-P_i} \cdot \cos(v_i)$$ $$X_{P_i} = X_A + S_{A-P_i} \cdot \cos(\alpha_{A-P_i})$$ $$Y_{P_i} = Y_A + S_{A-P_i} \cdot \sin(\alpha_{A-P_i})$$ $$H_{P_i} = H_A + S_{A-P_i} \cdot \tan(v_i) + i_m - l_g$$

Trong đó: $A(X_A, Y_A, H_A)$ là trạm máy; $B$ là điểm định hướng ngắm; $P_i$ là điểm chi tiết ranh thửa; $D_{A-P_i}$ là khoảng cách nghiêng; $v_i$ là góc nghiêng; $i_m$ là chiều cao máy; $l_g$ là chiều cao gương.

# Thuật toán tính tọa độ điểm chi tiết từ trạm máy toàn đạc điện tử
import math

def calculate_polar_point(station_coord, azimuth_rad, beta_deg, slope_dist, vert_angle_deg, inst_h, prism_h):
    x_st, y_st, h_st = station_coord
    beta_rad = math.radians(beta_deg)
    v_rad = math.radians(vert_angle_deg)
    
    # Góc phương vị từ trạm máy tới điểm đo
    total_azimuth = azimuth_rad + beta_rad
    
    # Khoảng cách ngang và chênh cao
    horizontal_dist = slope_dist * math.cos(v_rad)
    delta_x = horizontal_dist * math.cos(total_azimuth)
    delta_y = horizontal_dist * math.sin(total_azimuth)
    
    # Tọa độ không gian 3 chiều của điểm chi tiết
    target_x = x_st + delta_x
    target_y = y_st + delta_y
    target_h = h_st + (horizontal_dist * math.tan(v_rad)) + inst_h - prism_h
    
    return round(target_x, 3), round(target_y, 3), round(target_h, 3)

# Ví dụ thực thi với trạm đo A01 ngắm định hướng A02
station_A01 = (2415620.450, 485120.320, 265.400)
azimuth_init = math.radians(45.250)
point_detail = calculate_polar_point(station_A01, azimuth_init, 32.150, 45.820, 1.250, 1.520, 1.600)
print(f"Toạ độ điểm chi tiết (X, Y, H): {point_detail}")

Testing và validation

Số liệu đo lưới đường chuyền kinh vĩ cấp 2 tại phường Vĩnh Trại được kiểm tra nghiêm ngặt theo các tiêu chí kỹ thuật quy phạm:

Thông số kiểm tra Tiêu chuẩn TT 25/2014 Kết quả thực tế đo đạc Đánh giá
Sai số trung phương vị trí điểm ($M_p$) $\le 5,0\text{ cm}$ $2,8\text{ cm}$ Đạt tiêu chuẩn
Sai số khép tương đối đường chuyền ($f_s/\sum S$) $\le 1/25.000$ $1/32.400$ Vượt chỉ tiêu
Sai số khép góc đường chuyền ($f_\beta$) $\le 5\sqrt{n}" = 15"$ ($n=9$) $7,5"$ Đạt tiêu chuẩn
Dao động 2C trong 1 lần đo $\le 12"$ $6" - 8"$ Đạt tiêu chuẩn
Sai số trung phương đo khoảng cách ($M_s$) $\le 1,2\text{ cm}$ $0,6\text{ cm}$ Đạt tiêu chuẩn

Trong quá trình xử lý topology trên Gcadas:

  • Tự động phát hiện và sửa chữa 142 lỗi hở đỉnh ranh giới thửa đất (Node tolerance $\le 0,02\text{ m}$).
  • Khắc phục 28 lỗi chồng lấn ranh giới giữa đất ở đô thị (ODT) và đất giao thông (DGT).
  • Hoàn thiện 100% việc tạo vùng tự động (Generate Polygons) và gán nhãn số thửa, diện tích cho toàn bộ tờ bản đồ số 10.

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

Dự án đã nghiệm thu thành công sản phẩm bản đồ địa chính với các chỉ số định lượng:

  • Khối lượng thực hiện: Hoàn thành trọn vẹn tờ bản đồ số 10 tỷ lệ 1:500, bao quát diện tích 6,25 ha với 186 thửa đất độc lập.
  • Cơ cấu phân loại đất trên tờ bản đồ số 10: Đất ở tại đô thị (ODT) chiếm 72,4% diện tích đo vẽ (148 thửa), đất công trình sự nghiệp và hành chính chiếm 14,2%, đất giao thông nội bộ chiếm 13,4%.
  • Chất lượng hình học: 100% ranh giới công trình kiên cố (nhà bê tông, nhà gạch) được thể hiện chính xác theo mép tường ngoài, sai số ranh giới đo lặp kiểm tra tại thực địa nhỏ hơn $\pm 3\text{ cm}$.

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

  1. Chuẩn hóa quy trình liên thông số (Digital Pipeline): Loại bỏ hoàn toàn khâu chuyển vẽ thủ công bằng tay vốn dễ sinh sai số thô bạo. Toàn bộ chuỗi từ máy toàn đạc NTS-312B trút qua cáp RS-232/USB sang DPSurvey, tự động lọc điểm và đẩy vào lớp Level chuẩn (Level 10 - ranh thửa, Level 15 - nhà, Level 20 - giao thông) trong MicroStation V8i.
  2. Hiệu năng xử lý không gian vượt trội: Giảm 60% thời gian xử lý nội nghiệp so với phương pháp vẽ thủ công trên AutoCAD; giảm 85% các lỗi nhập sai số liệu thuộc tính nhờ tính năng đồng bộ trực tiếp từ bảng điều tra ngoại nghiệp Excel vào Gcadas.
  3. So sánh với giải pháp trước đây:
    • So với phương pháp bàn đạc truyền thống: Tốc độ hoàn thiện tăng gấp 4,5 lần; độ chính xác vị trí điểm tăng 80%.
    • So với phương pháp sử dụng GPS-RTK đơn lẻ: Giải quyết triệt để hiện tượng mất tín hiệu vệ tinh GNSS trong các ngõ ngách sâu dưới 2 m và các nhà cao tầng san sát tại trung tâm thành phố Lạng Sơn.

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

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

  • Cấp đổi GCNQSDĐ đồng loạt: Dữ liệu bản đồ số tờ 10 cung cấp hồ sơ kỹ thuật thửa đất chính xác cho 148 hộ gia đình tại khối 4 và khối 5 phường Vĩnh Trại, loại trừ hoàn toàn tình trạng tranh chấp lấn chiếm ranh giới ngõ đi chung.
  • Phục vụ đền bù giải phóng mặt bằng: Phục vụ công tác mở rộng tuyến đường giao thông kết nối quốc lộ 4B và quy hoạch chỉnh trang đô thị quanh khu vực chợ Đông Kinh và Ga Lạng Sơn.

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

  • Hạ tầng phần cứng: Máy tính trạm (Workstation) CPU Core i5/i7 thế hệ 8 trở lên, RAM tối thiểu 8GB, ổ cứng SSD 256GB, màn hình độ phân giải Full HD để biên tập CAD. Máy in khổ lớn A0/A1 phục vụ in ấn bản đồ giấy nghiệm thu.
  • Môi trường phần mềm: Hệ điều hành Windows 7/10 64-bit, Bentley MicroStation V8i (SELECTseries 3), bộ cài Gcadas 2015 đã đăng ký bản quyền, phần mềm DPSurvey, phông chữ TCVN3 / Unicode địa chính.
Hướng dẫn thiết lập môi trường biên tập Gcadas trên MicroStation V8i:
1. Sao chép Seedfile chuẩn VN-2000 (Seed2D_LangSon_107_15.dgn) vào thư mục: C:\Bentley\Workspace\Standards\seed\
2. Cài đặt Driver cáp trút máy toàn đạc USB-to-Serial Prolific PL2303.
3. Khởi động MicroStation V8i -> Chọn Utilities -> MDX Applications -> Browse nạp tập tin Gcadas.ma.
4. Chọn menu Gcadas -> Kết nối CSDL Access/SQL -> Thiết lập đơn vị hành chính: Lạng Sơn / TP Lạng Sơn / Phường Vĩnh Trại.
5. Thiết lập tỷ lệ hiển thị 1:500 và tải bảng phân lớp chuẩn Thông tư 25/2014.

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

  • Hạn chế kỹ thuật: Phương pháp toàn đạc điện tử quang học phụ thuộc lớn vào điều kiện thời tiết (sương mù vùng cao, mưa lớn làm tán xạ tia hồng ngoại EDM) và đòi hỏi thông hướng ngắm quang học trực tiếp giữa trạm máy và gương phản xạ, dẫn đến việc phải tăng dày nhiều điểm trạm phụ trong các hẻm dân cư phức tạp.
  • Hướng phát triển công nghệ:
    1. Kết hợp công nghệ bay chụp không người lái (UAV LiDAR / Photogrammetry) sử dụng trạm tham chiếu hoạt động liên tục CORS 4G RTK để quét địa hình tổng thể, kết hợp máy toàn đạc đo bù các góc khuất tầng trệt.
    2. Nâng cấp hệ thống CSDL địa chính lên mô hình WebGIS điện toán đám mây theo kiến trúc Microservices, cho phép cán bộ địa chính phường tra cứu và cập nhật biến động đất đai thời gian thực trên thiết bị di động.

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

  • Cơ quan quản lý nhà nước (UBND Phường & Phòng TN&MT): Sở hữu bộ CSDL bản đồ số có tính pháp lý cao, giảm 70% thời gian thụ lý hồ sơ tranh chấp, cấp mới và chuyển nhượng đất đai.
  • Kỹ sư Trắc địa - Bản đồ & Doanh nghiệp đo đạc: Tiếp cận quy trình kỹ thuật chuẩn hóa, tối ưu hóa năng suất lao động và giảm thiểu sai sót kỹ thuật trong khâu biên tập topology.
  • Sinh viên & Giảng viên ngành Quản lý đất đai: Tài liệu tham khảo toàn diện kết hợp giữa cơ sở lý thuyết trắc địa công trình và kỹ năng vận hành các phần mềm nghiệp vụ chuyên ngành.
  • Người dân sử dụng đất: Ranh giới thửa đất được xác lập minh bạch, chính xác đến từng centimet, bảo vệ quyền và lợi ích hợp pháp về tài sản gắn liền với đất.

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

1. Cần những điều kiện phần cứng và phần mềm cụ thể nào để triển khai biên tập bản đồ địa chính?
Hệ thống yêu cầu máy toàn đạc điện tử có độ chính xác đo góc từ 2" đến 5" (như SOUTH NTS-312B, Leica TCR, Topcon GTS), máy tính Windows Core i5, RAM $\ge 8\text{ GB}$, phần mềm MicroStation V8i tích hợp module Gcadas và DPSurvey cùng file seed chuẩn hệ tọa độ VN-2000 kinh tuyến trục địa phương.

2. Làm thế nào để kiểm soát và triệt tiêu sai số khi đo vẽ chi tiết trong khu vực đô thị nhà cao tầng san sát?
Áp dụng phương pháp bố trí lưới đường chuyền khép kín hoặc phù hợp, đo lặp 2 nửa lần đo (thuận kính và đảo kính để khử sai số 2C và góc nghiêng $i$), hạn chế chiều dài ngắm chi tiết dưới $150\text{ m}$ ở tỷ lệ 1:500 và sử dụng gương sào ngắn kết hợp dọi tâm quang học chính xác.

3. Làm cách nào để tích hợp tệp dữ liệu bản đồ DGN từ MicroStation vào hệ thống cơ sở dữ liệu quốc gia?
Gcadas cung cấp sẵn bộ công cụ kết xuất dữ liệu không gian sang định dạng ESRI Shapefile (*.shp) hoặc XML theo chuẩn mô hình dữ liệu ViLIS 2.0 và TMV.LIS, bảo toàn toàn bộ bảng thuộc tính chủ sở hữu, số tờ, số thửa và mã mục đích sử dụng đất.

4. Khi xuất hiện lỗi hở ranh (undershoot) hoặc bắt dính nhầm đỉnh thửa trên Gcadas thì khắc phục như thế nào?
Sử dụng công cụ Tìm và Sửa lỗi tự động trên menu Topology của Gcadas, thiết lập bán kính tìm kiếm lân cận (Search Tolerance) trong khoảng $0,01 - 0,03\text{ m}$. Đối với các vị trí phức tạp, chuyển sang chế độ sửa thủ công bằng công cụ bắt điểm Keypoint Snap trong MicroStation V8i.

5. Chi phí đầu tư cho một tổ đội đo đạc địa chính số hóa hoàn chỉnh là bao nhiêu?
Chi phí đầu tư thiết bị toàn đạc điện tử và phụ kiện dao động từ 60 - 120 triệu VNĐ, bản quyền phần mềm và máy tính đồ họa khoảng 30 - 50 triệu VNĐ. Với năng suất đo vẽ đạt 1,5 - 2,5 ha/ngày đối với tỷ lệ 1:500, thời gian thu hồi vốn đầu tư (ROI) trung bình đạt từ 4 đến 6 tháng thực hiện dự án.

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ờ bản đồ số 10 tỉ lệ 1:500 tại phường Vĩnh Trại, thành phố Lạng Sơn, tỉnh Lạng Sơn" đã hoàn thành xuất sắc các mục tiêu nghiên cứu và thực tiễn đề ra. Dự án chứng minh tính ưu việt vượt trội của việc tích hợp thiết bị đo đạc điện tử hiện đại SOUTH NTS-312B cùng giải pháp phần mềm liên thông MicroStation V8i - Gcadas, tạo ra sản phẩm bản đồ địa chính số hóa đạt độ chính xác cao, đúng tiêu chuẩn Thông tư 25/2014/TT-BTNMT. Kết quả của khóa luận không chỉ trực tiếp phục vụ công tác quản lý đất đai và cấp GCNQSDĐ tại địa phương mà còn cung cấp tài liệu kỹ thuật giá trị cho các nghiên cứu ứng dụng công nghệ thông tin trong ngành Địa chính - Trắc địa.