Giới thiệu dự án

Quy hoạch không gian nông - lâm nghiệp tại các địa bàn miền núi đòi hỏi độ chính xác cao về đặc tính địa hình nhằm giảm thiểu nguy cơ xói mòn, sạt lở và thoái hóa tài nguyên đất. Theo thống kê của Bộ Nông nghiệp và Phát triển Nông thôn, địa hình đồi núi chiếm tới hơn 74% diện tích tự nhiên Việt Nam; trong đó, phần lớn diện tích đất dốc đối mặt với áp lực rửa trôi nghiêm trọng vào mùa mưa và hạn hán cục bộ trong mùa khô. Huyện Hòa An, tỉnh Cao Bằng sở hữu tổng diện tích tự nhiên $60.710,33\text{ ha}$ với $21$ đơn vị hành chính (20 xã và 1 thị trấn). Đây là địa bàn có cấu trúc địa hình dạng lòng máng chia cắt phức tạp: địa hình đồi núi đất chiếm $63%$, núi đá vôi dốc đứng chiếm $20%$, và vùng bồn địa thung lũng chỉ chiếm $17%$.

Thực trạng canh tác truyền thống tại Hòa An còn mang tính tự phát, thiếu các công cụ định lượng khoa học về độ dốc, dẫn đến việc phân bố cây trồng chưa tối ưu và năng suất cây trồng trên đất dốc chưa tương xứng với tiềm năng. Phương pháp xây dựng bản đồ độ dốc thủ công từ bản đồ địa hình giấy trước đây bộc lộ nhiều hạn chế: thời gian thành lập kéo dài, sai số nội suy lớn và không thể tích hợp linh hoạt vào các hệ thống quản trị dữ liệu số.

Đề tài “Ứng dụng phần mềm Vertical Mapper trong xây dựng bản đồ độ dốc phục vụ công tác định hướng phát triển nông lâm nghiệp huyện Hòa An, tỉnh Cao Bằng” được thực hiện nhằm chuyển đổi toàn diện quy trình thành lập bản đồ địa hình sang mô hình số hóa không gian 3D.

+-----------------------------------------------------------------------------+
|                            MỤC TIÊU DỰ ÁN                                   |
+-----------------------------------------------------------------------------+
| 1. Xử lý & chuẩn hóa mô hình số độ cao DEM (ASTER GDEM 30m) trên VN-2000.   |
| 2. Ứng dụng Vertical Mapper 3.0 phân tích không gian, nội suy 5 cấp độ dốc. |
| 3. Xây dựng CSDL GIS không gian - thuộc tính cho 21 xã/thị trấn huyện Hòa An.|
| 4. Định hướng cơ cấu cây trồng, lâm nghiệp phòng hộ theo từng đai dốc.      |
+-----------------------------------------------------------------------------+

Dự án giới hạn phạm vi nghiên cứu trên toàn bộ ranh giới hành chính huyện Hòa An ($60.710,33\text{ ha}$), tích hợp cơ sở dữ liệu hiện trạng sử dụng đất $55.115,01\text{ ha}$ đất nông - lâm nghiệp (chiếm $90,78%$ diện tích tự nhiên) nhằm cung cấp công cụ ra quyết định chuẩn xác cho công tác quy hoạch phát triển bền vữ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 công nghệ nội suy số, các phương pháp truyền thống và công cụ phân tích không gian hiện hữu được đánh giá toàn diện:

Tiêu chí so sánh Phương pháp thủ công (Bản đồ giấy) Hệ thống ArcGIS / ArcView 3.x Giải pháp MapInfo + Vertical Mapper 3.0
Nguồn dữ liệu vào Bản đồ địa hình giấy, đo lưới thủ công DEM raster, Shapefile vector ASTER GDEM, DEM raster, MapInfo TAB
Thời gian xử lý 20 - 30 ngày/huyện 5 - 7 ngày/huyện 2 - 3 ngày/huyện
Độ chính xác nội suy Thấp, phụ thuộc chủ quan Rất cao, thuật toán phức tạp Cao, thuật toán tối ưu hóa lưới Grid
Yêu cầu phần cứng Không yêu cầu Cấu hình máy trạm cao Cấu hình phổ thông, xử lý mượt mà
Tích hợp CSDL Địa chính Khó đồng bộ hóa Tốt, nhưng chi phí bản quyền lớn Rất tốt, tương thích hoàn toàn CSDL địa chính

Yêu cầu kỹ thuật của hệ thống được phân loại theo mô hình MoSCoW:

  • Must-have: Chuyển đổi chuẩn xác hệ quy chiếu từ WGS84 sang VN-2000 (kinh tuyến trục $105^\circ45'$, múi chiếu $3^\circ$); nội suy và phân loại tự động 5 cấp độ dốc ($0^\circ-8^\circ$, $8^\circ-15^\circ$, $15^\circ-20^\circ$, $20^\circ-25^\circ$, $>25^\circ$); truy vấn diện tích độ dốc theo từng xã.
  • Should-have: Mô hình hóa địa hình 3D (3D Surface Drape); tạo bản đồ chuyên đề đa lớp (độ dốc, thủy hệ, ranh giới giao thông, địa danh).
  • Could-have: Tự động hóa trích xuất báo cáo thống kê diện tích sang định dạng Microsoft Excel .xls.
  • Won't-have: Mô phỏng dòng chảy lũ quét thời gian thực (Real-time dynamic hydro-modeling) trong giai đoạn hiện tại.

Thiết kế hệ thống

Kiến trúc xử lý dữ liệu được thiết kế theo luồng dữ liệu 4 tầng chuyên biệt:

[ASTER GDEM 30m / Bản đồ 364 / DGN] 
[Chuẩn hóa Hệ tọa độ VN-2000 & Cắt ranh giới huyện Hòa An]
[Phân tích Spatial Analysis & Nội suy Grid Surface (.GRD)]
[Vector hóa Polygon, Gán thuộc tính 5 Cấp dốc & CSDL Địa chính]
[Bản đồ chuyên đề độ dốc + Báo cáo thống kê diện tích 21 Xã]

Technology Stack ứng dụng trong dự án:

  • Hệ điều hành: Microsoft Windows 7 / Windows XP Professional SP3.
  • GIS Platform: MapInfo Professional v9.5 (Quản trị CSDL không gian và thuộc tính dạng .TAB, .DAT, .MAP, .ID, .IND).
  • Surface Modeling Module: Vertical Mapper v3.0 (Nội suy Grid, Spatial Analysis, phân tích độ dốc).
  • Data Conversion Tool: Global Mapper v14.0 (Truy xuất ASTER GDEM, chuyển đổi raster sang định dạng Vertical Mapper .GRD).
  • CAD/Data Integration: MicroStation SE (Xử lý bản đồ số địa giới 364, chuyển đổi file .DGN sang .MIF/.MID).
  • Báo cáo & Thống kê: Microsoft Excel 2007.

Cấu trúc Cơ sở Dữ liệu Bản đồ Độ dốc được chuẩn hóa:

Tên lớp dữ liệu Định dạng Kiểu đối tượng Mục đích quản lý
HoaAn-DEM2_Slope .TAB / .GRD Grid / Region Phân cấp và hiển thị dải màu cho 5 cấp độ dốc
HoaAn_Songsuoi .TAB Polyline Quản lý thủy hệ: Sông Bằng, Sông Hiến và các suối chính
Ranhgioi_giaothong .TAB Polyline / Line Hệ thống giao thông (Quốc lộ, Tỉnh lộ 203) và địa giới xã
HoaAn_Text .TAB Text / Point Ghi chú địa danh xã, thị trấn, đỉnh núi cao
Khung_Bandoc .TAB Line / Text Quản lý khung bản đồ, tỷ lệ thước, lưới kinh vĩ độ VN-2000
Ghi_chu .TAB Text Chú giải bản đồ, thông tin nguồn dữ liệu và pháp lý

Methodology

Phương pháp luận của đề tài áp dụng mô hình tuần tự dạng thác nước kỹ thuật (Technical Waterfall Pipeline) gồm 4 giai đoạn rõ ràng:

[Thu thập dữ liệu] ➔ [Xử lý & Nội suy Grid] ➔ [Xây dựng CSDL & Biên tập] ➔ [Đánh giá & Định hướng]
  (20/01 - 04/2013)      (05/2013 - 09/2013)        (10/2013 - 01/2014)          (02/2014 - 30/04/2014)

Quy trình quản lý chất lượng (QA) áp dụng nguyên tắc kiểm tra đóng khớp ranh giới hình học: Tổng diện tích các vùng đa giác độ dốc nội suy từ bản đồ không được lệch quá $0,5%$ so với số liệu kiểm kê đất đai hành chính năm 2012 của huyện ($60.710,33\text{ ha}$).

Implementation và kết quả

Development process

Quá trình triển khai kỹ thuật tập trung vào giải thuật nội suy bề mặt số và tính toán góc dốc (Slope Calculation). Mô hình độ cao số DEM dạng Raster sau khi trích xuất từ ASTER GDEM có kích thước ô lưới cell size $30\text{m} \times 30\text{m}$.

Thuật toán tính toán độ dốc trong Vertical Mapper sử dụng ma trận láng giềng $3 \times 3$ ô lưới (Moving Window Grid):

+----+----+----+
| z1 | z2 | z3 |
+----+----+----+
| z4 | z0 | z5 |
+----+----+----+
| z6 | z7 | z8 |
+----+----+----+

Độ thay đổi độ cao theo trục $X$ ($p$) và trục $Y$ ($q$) tại điểm trung tâm $z_0$ được tính bằng đạo hàm bậc một:

$$p = \frac{\partial z}{\partial x} = \frac{(z_3 + 2z_5 + z_8) - (z_1 + 2z_4 + z_6)}{8 \times \Delta x}$$

$$q = \frac{\partial z}{\partial y} = \frac{(z_6 + 2z_7 + z_8) - (z_1 + 2z_2 + z_3)}{8 \times \Delta y}$$

Góc dốc $\theta$ (đơn vị độ) được xác định bằng công thức:

$$\text{Slope} (\theta) = \arctan\left(\sqrt{p^2 + q^2}\right) \times \frac{180}{\pi}$$

Sau khi tính toán tự động giá trị góc dốc liên tục ($0^\circ \le \theta \le 74^\circ$), module thực hiện tái phân loại (Reclassify) và gán mã cấp độ dốc theo cấu trúc logic:

' Pseudo-code tái phân loại cấp độ dốc trên MapInfo / Vertical Mapper
If Slope_Value >= 0 And Slope_Value < 8 Then
    CapDoDoc = 1  ' Cấp I: Dốc nhẹ (< 8 độ)
ElseIf Slope_Value >= 8 And Slope_Value < 15 Then
    CapDoDoc = 2  ' Cấp II: Dốc vừa (8 - 15 độ)
ElseIf Slope_Value >= 15 And Slope_Value < 20 Then
    CapDoDoc = 3  ' Cấp III: Dốc hơi mạnh (15 - 20 độ)
ElseIf Slope_Value >= 20 And Slope_Value < 25 Then
    CapDoDoc = 4  ' Cấp IV: Dốc mạnh (20 - 25 độ)
ElseIf Slope_Value >= 25 Then
    CapDoDoc = 5  ' Cấp V: Dốc rất mạnh (> 25 độ)
End If

Trích xuất thống kê diện tích tự động theo đơn vị hành chính xã được thực hiện qua câu truy vấn SQL không gian (Spatial SQL Query) trong MapInfo Professional:

Select Xa_Table.Ten_Xa, 
       HoaAn_Slope.CapDoDoc, 
       Sum(Area(Overlap(Xa_Table.obj, HoaAn_Slope.obj), "hectare")) As DienTich_Ha
From Xa_Table, HoaAn_Slope
Where Xa_Table.obj Intersects HoaAn_Slope.obj
Group By Xa_Table.Ten_Xa, HoaAn_Slope.CapDoDoc
Into ThongKe_DoDoc_Xa

Testing và validation

Kết quả nội suy từ mô hình DEM và bản đồ độ dốc được kiểm định thông qua việc đối chiếu với số liệu tọa độ thực địa đo bằng máy định vị GPS cầm tay Garmin tại các vị trí mốc cao độ đặc trưng thuộc xã Bạch Đằng, xã Bế Triều và Thị trấn Nước Hai.

Bảng đối chuẩn tổng hợp diện tích giữa dữ liệu không gian bản đồ và dữ liệu kiểm kê thống kê địa chính:

Hạng mục kiểm tra Số liệu Thống kê Địa chính (2012) Số liệu Nội suy Bản đồ GIS Độ lệch (Delta) Sai số tương đối (%)
Tổng diện tích tự nhiên $60.710,33\text{ ha}$ $60.708,85\text{ ha}$ $-1,48\text{ ha}$ $0,0024%$ (Đạt chuẩn)
Diện tích xã Bạch Đằng $6.006,74\text{ ha}$ $6.006,51\text{ ha}$ $-0,23\text{ ha}$ $0,0038%$
Diện tích xã Dân Chủ $5.571,54\text{ ha}$ $5.571,62\text{ ha}$ $+0,08\text{ ha}$ $0,0014%$
Diện tích TT Nước Hai $122,61\text{ ha}$ $122,60\text{ ha}$ $-0,01\text{ ha}$ $0,0081%$

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

Hệ thống đã xây dựng thành công bộ bản đồ phân cấp độ dốc tỷ lệ 1:25.000 cho huyện Hòa An với 5 dải cấp độ dốc phân bố cụ thể:

  • Cấp I ($0^\circ - 8^\circ$ - Dốc nhẹ): Chiếm khoảng $17,2%$ diện tích ($10.442,18\text{ ha}$), phân bố tập trung tại dải thung lũng bồn địa ven sông Bằng và sông Hiến thuộc 8 xã trung tâm (Bế Triều, Hoàng Tung, TT Nước Hai, Hồng Việt...).
  • Cấp II ($8^\circ - 15^\circ$ - Dốc vừa): Chiếm $12,8%$ diện tích ($7.770,92\text{ ha}$), phân bố ở các vùng chân đồi tiếp giáp đồng bằng.
  • Cấp III ($15^\circ - 20^\circ$ - Dốc hơi mạnh): Chiếm $18,4%$ diện tích ($11.170,70\text{ ha}$), dạng đồi bát úp lượn sóng.
  • Cấp IV ($20^\circ - 25^\circ$ - Dốc mạnh): Chiếm $19,6%$ diện tích ($11.899,23\text{ ha}$), địa hình đồi núi đất phân cắt.
  • Cấp V ($> 25^\circ$ - Dốc rất mạnh, dải giá trị $25^\circ - 74^\circ$): Chiếm $32,0%$ diện tích ($19.425,82\text{ ha}$), tập trung ở các dãy núi đá vôi phía Đông, Đông Bắc (Bạch Đằng, Hà Trì, Ngũ Lão) và các đai rừng phòng hộ dốc đứng.

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

  1. Tự động hóa hoàn toàn quy trình xây dựng bản đồ địa hình chuyên đề: Thay thế quy trình đo lưới thủ công bằng việc kết hợp giải thuật nội suy Grid Surface của Vertical Mapper 3.0 với dữ liệu viễn thám độ cao ASTER GDEM, rút ngắn thời gian xử lý dữ liệu từ 30 ngày xuống còn 3 ngày làm việc (giảm $90%$ thời gian thực hiện).
  2. Độ chính xác và độ mịn không gian cao: Khác phục triệt để nhược điểm đứt gãy đường bình độ của phương pháp vector truyền thống; cung cấp khả năng mô phỏng bề mặt 3D giúp các nhà quản lý quan sát trực quan hướng dốc và độ dốc thực tế.
  3. Tích hợp đa chiều cơ sở dữ liệu: Liên kết thành công cơ sở dữ liệu không gian địa hình với CSDL thuộc tính kinh tế - xã hội của 21 đơn vị hành chính cấp xã, tạo nền tảng phục vụ quản lý chuyên sâu thay vì các bản đồ giấy đơn lớp tĩnh.

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

Căn cứ trên kết quả phân tích 5 cấp độ dốc, đề tài đề xuất ma trận định hướng sử dụng đất nông - lâm nghiệp bền vững cho huyện Hòa An:

+-----------------------------------------------------------------------------+
|               KHUNG ĐỊNH HƯỚNG QUY HOẠCH NÔNG - LÂM NGHIỆP                  |
+-----------------------------------------------------------------------------+
| • Cấp I (0° - 8°): Vùng thung lũng bồn địa                                  |
|   -> Chuyên canh lúa nước chất lượng cao (2 lúa - 1 màu), vùng trồng thuốc   |
|      lá xuất khẩu, rau củ ôn đới thương phẩm.                               |
| • Cấp II (8° - 15°): Vùng sườn đồi thoải                                    |
|   -> Canh tác ngô lai, đậu tương, lạc; phát triển cây ăn quả có múi kết     |
|      hợp các biện pháp làm đất theo đường đồng mức chống xói mòn.           |
| • Cấp III (15° - 20°): Vùng đồi lượn sóng                                   |
|   -> Phát triển mô hình Nông - Lâm kết hợp (Agroforestry), trồng chè, cây    |
|      công nghiệp ngắn ngày, xen canh cây họ đậu cải tạo đất.                |
| • Cấp IV (20° - 25°): Vùng đồi núi cao                                      |
|   -> Trồng rừng sản xuất (keo, thông, bạch đàn), phát triển đồng cỏ chăn    |
|      nuôi đại gia súc (trâu, bò thịt) tập trung.                            |
| • Cấp V (> 25°): Vùng núi đá vôi & đồi dốc đứng                              |
|   -> Khoanh nuôi tái sinh tự nhiên, bảo vệ tuyệt đối 45.401,56 ha rừng      |
|      phòng hộ đầu nguồn, nghiêm cấm nương rẫy du canh.                      |
+-----------------------------------------------------------------------------+

Đánh giá hiệu quả kinh tế - xã hội (ROI): Việc ứng dụng bản đồ số độ dốc giúp tiết kiệm khoảng $65%$ chi phí ngân sách dành cho công tác đo đạc, điều tra thực địa trong các kỳ quy hoạch sử dụng đất cấp huyện, đồng thời hạn chế tối đa tổn thất tài nguyên do quy hoạch sai vị trí cây trồng.

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

  • Hạn chế:
    • Độ phân giải không gian của ASTER GDEM ($30\text{m} \times 30\text{m}$) ở một số vùng địa hình karst đá vôi phức tạp chưa bóc tách được các khe rãnh cục bộ có bề rộng dưới $10\text{m}$.
    • Hệ thống chưa tích hợp cơ chế tự động cập nhật biến động thảm phủ thực vật theo thời gian thực (Real-time temporal monitoring).
  • Hướng phát triển:
    • Nâng cấp nguồn dữ liệu đầu vào sử dụng ảnh LiDAR hoặc dữ liệu UAV độ phân giải cao ($< 1\text{m}$) cho các khu vực trọng điểm đô thị hóa và sạt lở nguy hiểm.
    • Phát triển hệ thống WebGIS mã nguồn mở (GeoServer + Leaflet/OpenLayers) giúp cán bộ phòng Tài nguyên & Môi trường và người dân 21 xã có thể tra cứu trực tiếp thông tin độ dốc qua trình duyệt web.

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

  • Sinh viên & Học viên ngành Địa chính, Quản lý Đất đai, GIS: Nắm bắt quy trình thực hành chuẩn về xử lý mô hình DEM, làm chủ kỹ thuật nội suy Grid trên Vertical Mapper và MapInfo Professional.
  • Kỹ sư GIS & Đo đạc Bản đồ: Tiếp cận phương pháp chuyển đổi linh hoạt giữa các định dạng dữ liệu viễn thám raster (.GRD) và vector địa chính (.TAB, .DGN), xử lý chuẩn hóa hệ quy chiếu VN-2000.
  • Cơ quan quản lý Nhà nước (UBND huyện Hòa An, Sở TN&MT Cao Bằng): Sở hữu công cụ số hóa chuẩn xác làm căn cứ lập quy hoạch sử dụng đất giai đoạn 2015 - 2020 và xây dựng các chương trình mục tiêu quốc gia về Nông thôn mới.
  • Nhà nghiên cứu Sinh thái & Nông nghiệp: Có tập số liệu thực chứng về phân bố không gian địa hình phục vụ các mô hình dự báo xói mòn đất (RUSLE), phân tích thủy văn lưu vực sông Bằng.

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

1. Yêu cầu cấu hình hệ thống tối thiểu để vận hành giải pháp là gì?

Máy tính cá nhân trang bị CPU Intel Core 2 Duo hoặc Core i3 trở lên, RAM tối thiểu 2GB, ổ cứng trống 10GB và chạy hệ điều hành Windows XP/7/8/10. Phần mềm cần cài đặt: MapInfo Professional (v9.5 trở lên), module tích hợp Vertical Mapper 3.0 và Global Mapper v14.

2. Làm thế nào để xử lý triệt để sai số dịch chuyển hệ quy chiếu giữa WGS84 và VN-2000?

Dữ liệu ASTER GDEM gốc ở hệ tọa độ WGS84 toàn cầu (Elipsoid WGS84). Khi đưa vào Global Mapper, cần thiết lập thông số chuyển đổi 7 tham số (Datum Transformation) về Hệ quy chiếu Quốc gia VN-2000 với kinh tuyến trục địa phương tỉnh Cao Bằng là $105^\circ45'$, múi chiếu $3^\circ$ (hoặc $6^\circ$), sai số dịch chuyển vị trí sau chuyển đổi đạt mức $< 0,05\text{m}$.

3. Ưu điểm nổi bật của Vertical Mapper so với công cụ Spatial Analyst trên ArcGIS là gì?

Vertical Mapper được thiết kế gọn nhẹ, tương thích trực tiếp và liền mạch với cấu trúc bảng .TAB của MapInfo - nền tảng GIS tiêu chuẩn được ngành Địa chính và Tài nguyên Môi trường Việt Nam sử dụng rộng rãi, giúp việc chồng ghép bản đồ hiện trạng sử dụng đất với lưới độ dốc diễn ra nhanh chóng mà không cần chuyển đổi định dạng trung gian phức tạp.

4. Cơ sở dữ liệu bản đồ độ dốc này có thể cập nhật định kỳ như thế nào?

Khi có dữ liệu viễn thám độ cao mới hoặc bản đồ địa hình tỷ lệ lớn hơn (1:10.000, 1:5.000), người quản trị chỉ cần chạy lại luồng xử lý Grid Interpolation trong Vertical Mapper và thực hiện cập nhật lại bảng liên kết thuộc tính HoaAn-DEM2_Slope trong MapInfo mà không cần thay đổi cấu trúc CSDL cốt lõi.

5. Chi phí đầu tư triển khai phương pháp này so với phương pháp truyền thống ra sao?

Giải pháp tận dụng nguồn ảnh viễn thám độ cao ASTER GDEM miễn phí toàn cầu kết hợp phần mềm chuyên dụng, giúp giảm trên $60%$ chi phí nhân công và vật tư so với đo đạc thủ công ngoài thực địa, mang lại tỷ suất hoàn vốn đầu tư (ROI) rất cao cho các dự án quy hoạch cấp huyện.

Kết luận

Đề tài đã hoàn thành xuất sắc mục tiêu ứng dụng phần mềm Vertical Mapper 3.0 kết hợp MapInfo Professional v9.5 để xây dựng bản đồ độ dốc số hóa tỷ lệ 1:25.000 và hệ thống CSDL GIS địa hình cho huyện Hòa An, tỉnh Cao Bằng. Kết quả phân cấp 5 cấp độ dốc cung cấp luận cứ khoa học thực tiễn xác đáng cho việc tổ chức lại không gian sản xuất nông - lâm nghiệp, bảo vệ $45.401,56\text{ ha}$ đất rừng phòng hộ và phát huy tối đa tiềm năng canh tác vùng thung lũng bồn địa. Đây là mô hình kỹ thuật chuẩn xác, hiệu quả cao, sẵn sàng chuyển giao và nhân rộng cho các địa bàn miền núi có điều kiện tự nhiên tương đồng trên cả nước.