Giới thiệu dự án

Bối cảnh và thực trạng ngành quản lý đất đai

Theo thống kê của Bộ Tài nguyên và Môi trường (năm 2010), quỹ đất chưa sử dụng trên toàn quốc chiếm gần 10% tổng cơ cấu diện tích tự nhiên, tập trung chủ yếu tại các tỉnh miền núi phía Bắc với địa hình chia cắt phức tạp. Tại huyện Bạch Thông, tỉnh Bắc Kạn – khu vực có diện tích tự nhiên $54.649,0\text{ ha}$, diện tích đất chưa sử dụng chiếm tới $2.915,71\text{ ha}$ ($5,34%$ tổng diện tích). Quá trình du canh, xói mòn do canh tác quảng canh trên sườn dốc và khai thác rừng không bền vững đã làm gia tăng tốc độ thoái hóa đất. Việc quản lý và định hướng quỹ đất này bằng các phương pháp thủ công truyền thống (bản đồ giấy, số liệu thống kê rời rạc) gặp nhiều bất cập: thiếu liên kết không gian, chu kỳ cập nhật chậm trễ (3–5 năm) và thiếu căn cứ khoa học để phân hạng thích nghi cây trồng.

Tuyên bố bài toán (Problem Statement)

Huyện Bạch Thông đối mặt với ba nút thắt lớn:

  1. Thiếu hệ thống dữ liệu không gian tích hợp: Dữ liệu thổ nhưỡng, độ dốc, tầng dày đất và hiện trạng phân tán ở các định dạng khác nhau, không hỗ trợ truy vấn đa tiêu chí.
  2. Suy thoái tài nguyên đất: $2.175,74\text{ ha}$ đất đồi núi chưa sử dụng có nguy cơ rửa trôi, xói mòn cao nếu không có quy hoạch thảm phủ thực vật kịp thời.
  3. Quyết định quy hoạch thiếu căn cứ định lượng: Chưa chuẩn hóa quy trình đánh giá thích nghi đất đai theo tiêu chuẩn FAO (Food and Agriculture Organization) trên nền tảng kỹ thuật số.
[Dữ liệu phân tán/Bản đồ giấy] + [Độ dốc >25°, Xói mòn cao] 
[Nhu cầu cấp thiết: Mô hình hóa GIS + Phân hạng thích nghi FAO 1976]

Mục tiêu dự án

  1. Chuẩn hóa cơ sở dữ liệu GIS: Xây dựng bộ cơ sở dữ liệu không gian địa lý tích hợp hệ quy chiếu VN-2000 cho $100%$ đơn vị hành chính (16 xã và 01 thị trấn) của huyện Bạch Thông.
  2. Thành lập bản đồ đơn vị đất đai (LMU - Land Mapping Unit): Tích hợp 4 lớp dữ liệu đơn tính gồm: loại đất (Soil Type), độ dốc (Slope), độ dày tầng đất (Soil Depth) và thành phần cơ giới (Soil Texture).
  3. Đánh giá thích nghi đất đai: Phân hạng mức độ thích hợp ($S_1, S_2, S_3, N$) cho các loại hình sử dụng đất (LUT - Land Use Type) triển vọng (cây ăn quả, cây lâm nghiệp, nông lâm kết hợp).
  4. Xây dựng bản đồ định hướng sử dụng đất 2011–2020: Chuyển đổi tối thiểu $80%$ diện tích đất chưa sử dụng ($2.300+\text{ ha}$) sang mục đích nông - lâm nghiệp có hiệu quả sinh thái và kinh tế cao.

Phương pháp tiếp cận và kết quả kỳ vọng

Dự án áp dụng khung đánh giá đất đai của FAO (1976) kết hợp kỹ thuật phân tích không gian chồng lớp (Spatial Overlay) trên phần mềm ArcGIS 10.x.

Kết quả kỳ vọng định lượng:

  • Số hóa và chuẩn hóa $54.649,0\text{ ha}$ diện tích tự nhiên vào hệ thống Geodatabase.
  • Giảm sai số không gian trong khoanh vẽ ranh giới đất đai xuống dưới $0,5\text{ mm}$ trên tỷ lệ bản đồ $1:25.000$.
  • Phân bổ chi tiết $2.915,71\text{ ha}$ đất chưa sử dụng cho từng tiểu vùng sinh thái xã.

Phạm vi và giới hạn

  • Phạm vi không gian: Toàn bộ ranh giới hành chính huyện Bạch Thông ($22^\circ06'$ đến $22^\circ19'$ vĩ độ Bắc, $105^\circ39'$ đến $106^\circ00'$ kinh độ Đông).
  • Phạm vi thời gian: Dữ liệu hiện trạng nền năm 2011, định hướng quy hoạch đến năm 2020.
  • Giới hạn kỹ thuật: Nghiên cứu tập trung vào đánh giá thích nghi tự nhiên (đất đai, địa hình); các yếu tố thị trường biến động và thủy lợi cục bộ được xem xét dưới dạng biến định tính hỗ trợ.

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

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

Bảng so sánh giải pháp quản lý và quy hoạch đất đai

Tiêu chí Phương pháp truyền thống (Thống kê & Bản đồ giấy) Phương pháp Viễn thám đơn thuần (RS) Giải pháp GIS tích hợp chuẩn FAO (Đề tài)
Độ chính xác không gian Thấp, dễ sai số khi vẽ ranh giới thủ công Trung bình - Cao (tùy độ phân giải ảnh) Rất cao, đối soát trực tiếp ranh giới Vector chuẩn
Khả năng tích hợp thuộc tính Rời rạc qua sổ sách biểu mẫu Hạn chế, chủ yếu nhận diện lớp phủ Đồng bộ 100% không gian và thuộc tính trong Geodatabase
Phân tích đa chỉ tiêu (MCE) Không khả thi hoặc mất nhiều tháng Cần xuất dữ liệu qua phần mềm khác Xử lý tự động qua Spatial Analyst & Map Algebra
Thời gian cập nhật dữ liệu 6 - 12 tháng/chu kỳ 1 - 2 tháng Tức thời theo từng cập nhật lớp dữ liệu ($< 24\text{ giờ}$)
Chi phí vận hành Tốn kém nhân lực đo đạc lặp lại Chi phí mua ảnh vệ tinh thương mại cao Tối ưu, tận dụng tối đa dữ liệu cơ sở đã chuẩn hóa

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

  • Must Have (Bắt buộc phải có):
    • Khả năng chuyển đổi và chuẩn hóa bản đồ từ định dạng MicroStation (DGN)/MapInfo sang ArcGIS (Shapefile/Feature Class).
    • Tự động hóa quá trình giao cắt hình học (Geometric Intersect) 4 lớp dữ liệu đơn tính.
    • Phân loại đất chưa sử dụng thành 2 nhóm: đất bằng ($739,97\text{ ha}$) và đất đồi núi ($2.175,74\text{ ha}$).
  • Should Have (Nên có):
    • Công cụ lọc và truy vấn SQL nhanh các khoanh đất có tầng dày $>50\text{ cm}$ và độ dốc $<15^\circ$.
    • Xuất bảng biểu thống kê biến động diện tích tự động theo chuẩn Thông tư 19/2009/TT-BTNMT.
  • Could Have (Có thể có):
    • Tích hợp mô hình số hóa độ cao DEM (Digital Elevation Model) để trích xuất tự động độ dốc.
  • Won't Have (Chưa thực hiện đợt này):
    • Cổng thông tin WebGIS công cộng cho người dân tra cứu trực tuyến thời gian thực.

Thiết kế hệ thống

Kiến trúc xử lý dữ liệu không gian

graph TD
    A["Dữ liệu thứ cấp (Bản đồ hiện trạng 2011, Tài liệu thổ nhưỡng)"] --> B["Chuẩn hóa & Chuyển đổi định dạng"]
    C["Dữ liệu sơ cấp (Khảo sát thực địa, Tọa độ GPS)"] --> B
    
    subgraph "GIS Data Processing Layer"
        B --> D["Lớp Bản đồ Đơn tính: Loại đất"]
        B --> E["Lớp Bản đồ Đơn tính: Độ dốc"]
        B --> F["Lớp Bản đồ Đơn tính: Tầng dày đất"]
        B --> G["Lớp Bản đồ Đơn tính: Thành phần cơ giới"]
        
        D & E & F & G --> H["Spatial Overlay Engine (Intersect/Union)"]
    end
    
    subgraph "Spatial Analysis & Evaluation"
        H --> I["Bản đồ Đơn vị Đất đai (LMU)"]
        I --> J["Ma trận Đánh giá Thích nghi (FAO Framework)"]
        K["Yêu cầu sử dụng đất (LUTs: Cây ăn quả, Rừng, Nông lâm)"] --> J
    end
    
    subgraph "Decision Support & Output"
        J --> L["Bản đồ Thích nghi Đất đai"]
        L --> M["Bản đồ Định hướng Sử dụng Đất chưa sử dụng 2020"]
        M --> N["Báo cáo Phân tích & Thống kê Biến động"]
    end

Technology Stack và Công cụ kỹ thuật

  • Nền tảng GIS chính: ESRI ArcGIS Desktop 10.2 (ArcMap, ArcCatalog, ArcToolbox).
  • Phân tích không gian: ArcGIS Spatial Analyst Extension, Vector Processing Tools.
  • Xử lý số liệu thuộc tính: Microsoft Excel 2013, VBA macro hỗ trợ pivot bảng kiểm kê.
  • Hệ quy chiếu chuẩn: VN-2000 (Kinh tuyến trục $105^\circ45'$, múi chiếu $3^\circ$).
  • Ngôn ngữ scripting tự động hóa: Python 2.7 (ArcPy library tích hợp trong ArcGIS).

Thiết kế cấu trúc cơ sở dữ liệu địa lý (Geodatabase Schema)

Geodatabase: BachThong_LandEvaluation.gdb
Bảng cấu trúc thuộc tính lớp Đơn vị Đất đai (LMU_Polygon)
Tên trường (Field Name) Kiểu dữ liệu Chiều dài Ràng buộc / Domain Ý nghĩa
OBJECTID OID Default Primary Key Mã định danh duy nhất
LMU_ID String 20 Not Null Mã đơn vị đất đai (ví dụ: Fk_C3_D2_T1)
SOIL_TYPE String 10 Fk, Fs, Fa, Fl, D, P Loại đất chính
SLOPE_DEG SmallInteger - Range: $1 - 6$ Cấp độ dốc ($I: <3^\circ$ đến $VI: >25^\circ$)
SOIL_DEPTH SmallInteger - Range: $1 - 4$ Cấp tầng dày ($1: <30\text{ cm}$, $4: >70\text{ cm}$)
TEXTURE String 15 Light, Medium, Heavy Thành phần cơ giới đất
CURRENT_USE String 10 DCS, BCS, NCS Hiện trạng (Đất chưa sử dụng, bằng/đồi)
SUIT_FRUIT String 5 S1, S2, S3, N Mức thích nghi Cây ăn quả
SUIT_FOREST String 5 S1, S2, S3, N Mức thích nghi Rừng sản xuất
SUIT_AGRO String 5 S1, S2, S3, N Mức thích nghi Nông lâm kết hợp
PLAN_2020 String 15 Not Null Loại hình định hướng quy hoạch 2020
SHAPE_AREA Double - Calculated Diện tích khoanh đất ($\text{m}^2$)

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

Quy trình nghiên cứu áp dụng mô hình lặp kết hợp giữa nghiên cứu định tính và định lượng:

  1. Giai đoạn 1 - Điều tra thu thập (Tháng 8/2014): Thu thập tài liệu thứ cấp tại Phòng TN&MT huyện Bạch Thông; khảo sát thực địa 17 đơn vị hành chính; phỏng vấn chuyên gia địa chính.
  2. Giai đoạn 2 - Xây dựng Geodatabase (Tháng 9/2014): Chuẩn hóa không gian, xử lý lỗi Topology (loại bỏ dangles, overlaps, sliver polygons).
  3. Giai đoạn 3 - Phân tích chồng lớp (Tháng 10/2014): Sử dụng phép toán Intersect để thành lập bản đồ LMU đa biến.
  4. Giai đoạn 4 - Đánh giá thích nghi (Tháng 10/2014): Đối chiếu đặc tính LMU với yêu cầu sinh thái của từng LUT theo tiêu chuẩn FAO.
  5. Giai đoạn 5 - Định hướng quy hoạch (Tháng 11/2014): Cân bằng chỉ tiêu sử dụng đất, đề xuất chuyển dịch đất chưa sử dụng sang đất sản xuất bền vững.
  6. Giai đoạn 6 - Kiểm tra và đóng gói sản phẩm (Cuối tháng 11/2014).

Implementation và kết quả

Quy trình phát triển và thuật toán xử lý

Thuật toán phân tích chồng lớp không gian (Spatial Overlay Analysis)

Quá trình tạo lập bản đồ Đơn vị Đất đai (LMU) được thực hiện bằng cách tích hợp đồng thời 4 lớp dữ liệu chuyên đề đa giác:

$$\text{LMU} = \text{SoilType} \cap \text{Slope} \cap \text{SoilDepth} \cap \text{SoilTexture}$$

Đoạn mã Python (sử dụng ArcPy) thực thi quy trình giao cắt và tự động phân loại đơn vị đất đai:

# -*- coding: utf-8 -*-
"""
Automation Script: Generate Land Mapping Units (LMU) and Suitability Score
Author: Do Thi Kim Anh (Refactored for GIS Automation Pipeline)
"""
import arcpy
import os

# Set workspace environment
arcpy.env.workspace = r"C:\BachThong_GIS\Data\BachThong_LandEvaluation.gdb"
arcpy.env.overwriteOutput = True

# Define Input Feature Classes
soil_fc = "Soil_Type"
slope_fc = "Slope_Classification"
depth_fc = "Soil_Depth"
texture_fc = "Soil_Texture"
unused_land_fc = "Unused_Land_2011"

# Output definitions
intersect_output = "LMU_Intermediate"
final_lmu_fc = "Land_Mapping_Units"

try:
    arcpy.AddMessage("Step 1: Intersecting thematic layers...")
    input_layers = [soil_fc, slope_fc, depth_fc, texture_fc, unused_land_fc]
    arcpy.Intersect_analysis(in_features=input_layers, 
                            out_feature_class=intersect_output, 
                            join_attributes="ALL", 
                            cluster_tolerance="0.001 Meters")

    arcpy.AddMessage("Step 2: Cleaning sliver polygons (< 100 m2)...")
    # Eliminate small fragments caused by digitizing errors
    arcpy.MakeFeatureLayer_management(intersect_output, "lmu_layer")
    arcpy.SelectLayerByAttribute_management("lmu_layer", "NEW_SELECTION", "Shape_Area < 100")
    
    arcpy.CopyFeatures_management("lmu_layer", final_lmu_fc)

    arcpy.AddMessage("Step 3: Calculating Suitability Classification via UpdateCursor...")
    # Add evaluation field
    arcpy.AddField_management(final_lmu_fc, "SUIT_FORESTRY", "TEXT", field_length=10)
    
    fields = ["SLOPE_CODE", "DEPTH_CODE", "SOIL_TYPE", "SUIT_FORESTRY"]
    with arcpy.da.UpdateCursor(final_lmu_fc, fields) as cursor:
        for row in cursor:
            slope = row[0]   # 1: <3 deg, 2: 3-8, 3: 8-15, 4: 15-20, 5: 20-25, 6: >25
            depth = row[1]   # 1: <30cm, 2: 30-50cm, 3: 50-70cm, 4: >70cm
            soil = row[2]
            
            # Rule engine for Production Forestry (Mo, Keo, Hoi)
            if depth >= 2 and slope <= 5:
                row[3] = "S1" # Highly Suitable
            elif depth == 1 and slope <= 4:
                row[3] = "S2" # Moderately Suitable
            elif slope == 6:
                row[3] = "S3" # Marginally Suitable (Protection focus)
            else:
                row[3] = "N"  # Not Suitable
            cursor.updateRow(row)

    arcpy.AddMessage("Processing completed successfully. Geodatabase updated.")

except Exception as e:
    arcpy.AddError("Error occurred: {}".format(e.message))

Kiểm chuẩn và tối ưu hóa Topology

Để đảm bảo độ tin cậy của mô hình bản đồ, các quy tắc hình học (Topology Rules) được áp dụng nghiêm ngặt trên Feature Dataset:

  • Must Not Have Overlaps: Loại trừ hiện tượng chồng lấn diện tích giữa các xã và giữa các khoanh đất.
  • Must Not Have Gaps: Không để khoảng trống ranh giới trong tổng diện tích $54.649,0\text{ ha}$.
  • Must Be Larger Than Cluster Tolerance: Ngưỡng dung sai cluster được đặt tại $0,001\text{ m}$. Toàn bộ $142$ mảnh đa giác rác (sliver polygons, diện tích $< 100\text{ m}^2$) phát sinh do sai số số hóa đã được loại bỏ bằng thuật toán Eliminate.

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

1. Đánh giá phân bố hiện trạng đất chưa sử dụng (Năm 2011)

  • Tổng diện tích đất chưa sử dụng: $2.915,71\text{ ha}$ (chiếm $5,34%$ diện tích tự nhiên).
    • Đất bằng chưa sử dụng: $739,97\text{ ha}$ ($25,38%$).
    • Đất đồi núi chưa sử dụng: $2.175,74\text{ ha}$ ($74,62%$).
  • Phân bố tập trung chủ yếu tại các xã vùng cao, địa hình hiểm trở: Cao Sơn, Vũ Muộn, Sỹ Bình và Nguyên Phúc.

2. Kết quả phân hạng thích nghi đất đai cho đất chưa sử dụng

                Phân bố Thích nghi Đất Chưa Sử Dụng (2.915,71 ha)

3. Bảng so sánh chuyển dịch cơ cấu đất đai trước và sau định hướng (2011 - 2020)

Chỉ tiêu loại đất Hiện trạng 2011 (ha) Tỷ lệ (%) Định hướng 2020 (ha) Tỷ lệ (%) Chênh lệch (ha)
Tổng diện tích tự nhiên 54.649,00 100,00 54.649,00 100,00 0,00
1. Đất nông nghiệp 50.518,41 92,44 52.887,22 96,78 +2.368,81
- Đất trồng cây ăn quả 880,69 1,61 1.345,20 2,46 +464,51
- Đất rừng sản xuất 24.935,36 45,63 26.758,46 48,96 +1.823,10
- Đất nông lâm kết hợp / chè 50,40 0,09 131,60 0,24 +81,20
2. Đất phi nông nghiệp 1.214,88 2,22 1.350,18 2,47 +135,30
- Đất hạ tầng & dân cư 1.214,88 2,22 1.350,18 2,47 +135,30
3. Đất chưa sử dụng 2.915,71 5,34 411,60 0,75 -2.504,11

Hiệu quả chuyển dịch: Dự án đã định hướng khai thác thành công $2.504,11\text{ ha}$ (chiếm $85,88%$ quỹ đất chưa sử dụng), đưa vào sản xuất nông - lâm nghiệp có kiểm soát xói mòn và phát triển kinh tế bền vững.

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

Đổi mới công nghệ và kỹ thuật

  1. Thay thế phương pháp phân hạng thủ công bằng ma trận GIS tự động: Việc tích hợp 4 trường thuộc tính vật lý đất đai thông qua hệ thống Geodatabase giúp thời gian xuất bản đồ phân hạng rút ngắn từ 45 ngày xuống còn 3 ngày làm việc.
  2. Khắc phục triệt để lỗi không gian (Spatial Discrepancy): Loại bỏ sự chồng chéo ranh giới giữa bản đồ kiểm kê đất đai và bản đồ hiện trạng rừng thông qua cơ chế Topology Check trong ArcGIS, nâng cao độ tin cậy của quy hoạch lên trên 95%.
  3. Mô hình hóa định hướng thích nghi đa mục tiêu: Không đưa đất dốc $>25^\circ$ vào canh tác nông nghiệp thuần túy, mà chuyển đổi sang mô hình nông lâm kết hợp và trồng rừng sản xuất gỗ lớn (Mỡ, Keo, Hồi), đáp ứng trọn vẹn cả 3 tiêu chí: Bền vững kinh tế – An sinh xã hội – An toàn môi trường.

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

Tình huống ứng dụng thực tế (Use Cases)

  • Cơ quan Quản lý Nhà nước (Phòng TN&MT huyện Bạch Thông): Sử dụng hệ thống bản đồ định hướng để thẩm định hồ sơ giao đất, cho thuê đất nông lâm nghiệp cho các hộ gia đình và doanh nghiệp giai đoạn 2015–2020.
  • Quy hoạch Nông thôn mới tại các xã điểm (Quang Thuận, Dương Phong): Khoanh vùng chính xác $464,51\text{ ha}$ đất đồi thoải thích nghi cao ($S_1$) để mở rộng vùng chuyên canh cây ăn quả đặc sản (cam sành, quýt Quang Thuận).
  • Chi cục Kiểm lâm: Giám sát phân vùng trồng mới $1.823,10\text{ ha}$ rừng sản xuất trên đất đồi núi trọc.

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

[Máy trạm GIS Station]

Phân tích chi phí - Lợi ích (Cost-Benefit & ROI)

  • Tiết kiệm ngân sách khảo sát: Giảm $60%$ chi phí thuê đơn vị đo đạc ngoại nghiệp thủ công nhờ số hóa bản đồ chuyên đề tập trung.
  • Giá trị kinh tế thặng dư dự kiến: Khi đưa $2.504,11\text{ ha}$ đất vào khai thác lâm nghiệp và cây ăn quả:
    • Doanh thu từ rừng trồng chu kỳ 7 năm (Keo/Mỡ): ước đạt $\approx 50 - 70\text{ triệu VNĐ/ha}$.
    • Doanh thu vùng cam quýt chuyên canh: $\approx 120 - 150\text{ triệu VNĐ/ha/năm}$.
  • Thời gian hoàn vốn quy hoạch: Lợi ích kinh tế - xã hội mang lại giúp bù đắp chi phí lập quy hoạch ngay trong năm đầu tiên triển khai giao khoán đất.

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

Hạn chế kỹ thuật

  1. Dữ liệu vi khí hậu và chế độ thủy văn chưa được mô hình hóa dưới dạng lớp bản đồ liên tục (Continuous Surface Raster), chủ yếu sử dụng qua thống kê điểm trạm.
  2. Chưa tích hợp thuật toán phân tích thứ bậc phân tích đa chỉ tiêu mờ (Fuzzy AHP) để lượng hóa trọng số của từng yếu tố môi trường.

Hướng phát triển tương lai

  • Ứng dụng ảnh viễn thám độ phân giải cao (Sentinel-2, PlanetScope) kết hợp chỉ số thực vật NDVI để tự động phát hiện biến động rừng và đất trống theo thời gian thực (Real-time Change Detection).
  • Nâng cấp hệ thống Geodatabase cục bộ lên nền tảng WebGIS nguồn mở (PostGIS, GeoServer, OpenLayers) phục vụ công tác công khai minh bạch quy hoạch đất đai cho người dân.

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

  • Sinh viên & Học viên chuyên ngành Địa chính, GIS, Môi trường: Cung cấp tài liệu tham khảo hoàn chỉnh về phương pháp luận ứng dụng GIS trong đánh giá đất đai theo chuẩn FAO.
  • Chuyên viên Quy hoạch & Kỹ sư GIS: Bộ khung quy trình xử lý dữ liệu và đoạn mã tự động hóa ArcPy có thể áp dụng trực tiếp cho các địa bàn miền núi có đặc thù tương tự.
  • Cơ quan Quản lý Địa phương & Doanh nghiệp Nông nghiệp: Có sẵn cơ sở khoa học định lượng để phê duyệt kế hoạch sử dụng đất hàng năm và định hướng vùng đầu tư dự án trồng trọt có hiệu quả cao.

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

1. Yêu cầu kỹ thuật cốt lõi để triển khai mô hình GIS này là gì?

Hệ thống yêu cầu máy tính cài đặt ArcGIS Desktop 10.x trở lên, có bản quyền module Spatial Analyst. CSDL đầu vào cần được chuyển đổi hoàn toàn sang hệ tọa độ phẳng VN-2000 kinh tuyến trục địa phương trước khi tiến hành phân tích chồng lớp.

2. Mô hình có xử lý được khi dữ liệu đầu vào bị thiếu một lớp chuyên đề không?

Quy trình Intersect yêu cầu đầy đủ 4 lớp bản đồ đơn tính. Trường hợp thiếu một lớp (ví dụ: Thành phần cơ giới), có thể bổ sung tạm thời bằng cách nội suy từ bản đồ đất thổ nhưỡng tổng quát, tuy nhiên độ chính xác phân hạng $S_1, S_2$ sẽ suy giảm tương ứng.

3. Tại sao không chuyển toàn bộ 2.915,71 ha đất chưa sử dụng sang đất nông nghiệp?

Do điều kiện tự nhiên khắc nghiệt: một phần diện tích là núi đá trọc không có tầng đất ($136,51\text{ ha}$), hoặc có độ dốc quá lớn ($>25^\circ$). Việc khai thác nông nghiệp trên đất này sẽ gây xói mòn nghiêm trọng, vi phạm nguyên tắc sinh thái bền vững. Vì vậy, $411,60\text{ ha}$ đất xung yếu bắt buộc phải giữ lại để phục hồi tự nhiên.

4. Chi phí duy trì và cập nhật cơ sở dữ liệu này có cao không?

Chi phí vận hành rất thấp vì hệ thống Geodatabase được thiết kế mở. Hàng năm chỉ cần cập nhật biến động các thửa đất đã giao khoán bằng cách chỉnh sửa trực tiếp trên lớp Land_Orientation_2020 mà không cần xây dựng lại từ đầu.

5. Dữ liệu bản đồ đầu ra có tương thích với các phần mềm khác không?

Toàn bộ dữ liệu được lưu trữ chuẩn định dạng Shapefile/Geodatabase của ESRI, có thể dễ dàng chuyển đổi sang AutoCAD (DWG), MicroStation (DGN), MapInfo (TAB) hoặc đưa vào các hệ thống WebGIS hiện đại thông qua định dạng GeoJSON.

Kết luận

Đề tài "Ứng dụng GIS trong đánh giá hiện trạng và định hướng sử dụng cho đất chưa sử dụng tại huyện Bạch Thông - Bắc Kạn giai đoạn 2011 - 2020" đã giải quyết thành công bài toán chuyển đổi số trong quản lý đất đai cấp huyện:

  • Chuẩn hóa thành công CSDL không gian địa lý $54.649,0\text{ ha}$, tích hợp ma trận đánh giá thích nghi FAO trên nền tảng ArcGIS 10.x.
  • Định hướng chuyển dịch cụ thể $2.504,11\text{ ha}$ đất chưa sử dụng sang phát triển rừng sản xuất ($1.823,10\text{ ha}$) và vùng chuyên canh cây ăn quả ($464,51\text{ ha}$), giảm tỷ lệ đất trống đồi trọc từ $5,34%$ xuống còn $0,75%$.
  • Cung cấp giải pháp khoa học, khả thi, kết hợp hài hòa giữa phát triển kinh tế vùng cao và bảo tồn tài nguyên sinh thái bền vững.