Giới thiệu dự án

Tài nguyên rừng đóng vai trò cốt lõi trong việc duy trì cân bằng sinh thái, bảo vệ nguồn nước, chống xói mòn đất và giảm thiểu tác động tiêu cực của biến đổi khí hậu toàn cầu. Tại vùng Đông Bắc Việt Nam, đặc biệt là tỉnh Cao Bằng với địa hình đồi núi hiểm trở và hệ sinh thái karst (núi đá vôi) đặc trưng, công tác quản trị và giám sát tài nguyên rừng đặt ra những thách thức lớn đối với ngành lâm nghiệp. Huyện Bảo Lạc có diện tích tự nhiên $92.063,68\text{ ha}$ gồm 1 thị trấn và 16 xã, địa hình bị chia cắt sâu sắc bởi hệ thống sông suối và các dãy núi cao, gây khó khăn nghiêm trọng cho hoạt động tuần tra, kiểm kê thực địa bằng phương pháp thủ công truyền thống.

+-----------------------------------------------------------------------------------+
|                            BỐI CẢNH VÀ VẤN ĐỀ NGHIÊN CỨU                         |
+-----------------------------------------------------------------------------------+
|  1. Địa bàn nghiên cứu: Huyện Bảo Lạc, tỉnh Cao Bằng (Tổng diện tích 92.063,68 ha)|
|  2. Phương pháp cũ: Đo vẽ dã ngoại thủ công, thống kê giấy phân tán, độ trễ cao  |
|  3. Vấn đề cốt lõi: Sai số không gian lớn, thiếu CSDL số chuẩn hóa theo VN-2000   |
|  4. Giải pháp ứng dụng: Tích hợp công nghệ 3S (GIS - Viễn thám SPOT 4 - GPS)     |
+-----------------------------------------------------------------------------------+

Phương pháp kiểm kê lâm nghiệp truyền thống phụ thuộc vào việc đo đạc ngoại nghiệp dã ngoại bằng la bàn, thước dây và bản đồ giấy tĩnh. Cách tiếp cận này bộc lộ nhiều điểm nghẽn kỹ thuật:

  • Thời gian xử lý kéo dài và chi phí nhân công cao: Không thể cập nhật biến động theo chu kỳ ngắn hạn.
  • Độ chính xác không gian thấp: Khó bóc tách chính xác ranh giới các trạng thái rừng ở vùng hiểm trở, vùng núi đá vôi không có đường tiếp cận.
  • Thiếu tính đồng bộ dữ liệu: Không tích hợp được lớp dữ liệu không gian (spatial data) với cơ sở dữ liệu thuộc tính (attribute database) theo quy chuẩn quốc gia.

Đề tài "Đánh giá tình hình diễn biến diện tích và trạng thái rừng bằng phần mềm GIS tại huyện Bảo Lạc, tỉnh Cao Bằng giai đoạn 2012–2014" được triển khai nhằm giải quyết các tồn tại trên thông qua ứng dụng công nghệ tích hợp 3S (Remote Sensing - Viễn thám, Geographic Information System - Hệ thông tin địa lý, Global Positioning System - Hệ thống định vị toàn cầu).

Mục tiêu nghiên cứu cụ thể

  1. Chuẩn hóa công nghệ giải đoán ảnh viễn thám: Khai thác tư liệu ảnh vệ tinh SPOT 4 (độ phân giải không gian $20\text{ m}$) năm 2014 kết hợp hệ thống định vị GPS ngoại nghiệp để phân loại chính xác các trạng thái rừng và kiểu sử dụng đất.
  2. Thành lập cơ sở dữ liệu hiện trạng rừng 2014: Xác định chi tiết diện tích các loại đất lâm nghiệp, loại rừng theo 17 đơn vị hành chính cấp xã, theo chức năng quy hoạch 3 loại rừng (đặc dụng, phòng hộ, sản xuất) và theo chủ thể quản lý.
  3. Phân tích không gian đánh giá biến động giai đoạn 2012–2014: Ứng dụng kỹ thuật chồng xếp bản đồ (spatial overlay/intersect) giữa hai kỳ dữ liệu 2012 và 2014 để định lượng diện tích tăng/giảm của từng trạng thái rừng.
  4. Biên tập hệ thống bản đồ thành quả: Xây dựng bản đồ hiện trạng và bản đồ biến động tài nguyên rừng cấp huyện (tỷ lệ 1:50.000) và cấp xã (tỷ lệ 1:10.000) trên nền Hệ quy chiếu quốc gia VN-2000.

Phạm vi và giới hạn đề tài: Nghiên cứu bao phủ $92.063,68\text{ ha}$ diện tích tự nhiên của huyện Bảo Lạc trong khung thời gian từ 2012 đến 2014. Đề tài tập trung vào việc giải đoán phổ ảnh vệ tinh quang học SPOT 4 kết hợp thuật toán phân loại có giám định (Supervised Classification), chuyển đổi mô hình dữ liệu Raster sang Vector để tính toán không gian trên phần mềm ArcView GIS 3.2a, ArcGIS 10.1 và biên tập chuyên đề trên MapInfo Professional 10.0.

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

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

Trước khi công nghệ 3S được áp dụng rộng rãi, công tác theo dõi diễn biến rừng tại Cao Bằng thực hiện theo Quyết định số 78/2002/QĐ-BNN-KL chủ yếu dựa trên báo cáo hành chính từ cấp kiểm lâm địa bàn tổng hợp lên. Khi có biến động do cháy rừng, khai thác hoặc chuyển mục đích làm nương rẫy, việc cập nhật ranh giới lô rừng thường ước lượng bằng mắt hoặc đo vẽ sơ sài, dẫn đến hiện tượng sai lệch số liệu thực tế so với sổ sách quản lý.

Tiêu chí so sánh Đo vẽ thủ công & Bản đồ giấy Phương pháp Viễn thám & GIS (Đề tài ứng dụng) Đo đạc viễn thám UAV / LiDAR cao cấp
Độ phân giải / Chi tiết Thấp, sai số biên lớn ($> 50\text{ m}$) Trung bình - Cao (Ảnh SPOT $20\text{ m}$, bản đồ 1:10.000) Rất cao ($0,05 - 0,5\text{ m}$)
Thời gian cập nhật $6 - 12\text{ tháng}$ cho cấp huyện $5 - 10\text{ ngày}$ xử lý số liệu toàn huyện $1 - 3\text{ ngày}$ theo từng khu vực quét
Chi phí triển khai Tốn kém nhân lực đi thực địa dã ngoại Tiết kiệm $60 - 75%$ chi phí vận hành Chi phí thiết bị và bay quét rất cao
Khả năng phân tích không gian Rất hạn chế, không truy vấn tự động Tự động chồng xếp, tính toán diện tích bằng thuật toán Rất cao, phân tích mô hình 3D tán cây
Tính khả thi thực tế Lạc hậu, dễ gây sai lệch thống kê Tối ưu cho cơ quan quản lý cấp Tỉnh / Huyện Chưa khả thi trên diện rộng ($>90.000\text{ ha}$)

Phân tích yêu cầu hệ thống theo mô hình MoSCoW:

  • Must-have (Bắt buộc phải có): Chuẩn hóa toàn bộ dữ liệu trên Hệ quy chiếu quốc gia VN-2000 (kinh tuyến trục phù hợp với Cao Bằng theo Thông tư 973/2001/TT-TCĐC); thuật toán phân loại ảnh viễn thám đạt độ tin cậy $>80%$; tính toán chính xác số liệu biến động diện tích theo từng lô và từng xã.
  • Should-have (Nên có): CSDL thuộc tính liên kết tự động giữa mã trạng thái rừng ($1100, 1150, 1210, 1220, 1230, 1240, 2000, 2020, 2050, 3000$) với bảng tra cứu trữ lượng và đường kính, chiều cao bình quân ($D_{bq}, H_{bq}$).
  • Could-have (Có thể mở rộng): Tích hợp giao diện truy vấn không gian trực quan cho cán bộ Hạt Kiểm lâm không chuyên sâu về GIS.
  • Won't-have (Chưa thực hiện kỳ này): Nhận dạng tự động loại cây rừng ở quy mô dưới tán (sub-canopy) bằng mô hình Deep Learning.

Thiết kế hệ thống

Kiến trúc luồng xử lý dữ liệu lâm nghiệp được thiết kế theo chu trình 4 tầng: Thu thập $\rightarrow$ Tiền xử lý $\rightarrow$ Phân loại số & Chuyển đổi không gian $\rightarrow$ Chồng xếp & Biên tập bản đồ chuyên đề.

Technology Stack & Công cụ sử dụng:

  • ERDAS IMAGINE 9.2: Module xử lý ảnh số, nắn chỉnh tọa độ hình học (Geometric Correction), tổ hợp dải phổ (Layer Stack), trích xuất chữ ký phổ (Signature Editor), phân loại có kiểm định (Supervised Classification).
  • ArcView GIS 3.2a / ArcGIS 10.1 (Spatial Analyst): Xử lý không gian, chuyển đổi raster sang vector polygon, thực hiện phép giao cắt hình học (Geoprocessing Intersect Wizard).
  • MapInfo Professional 10.0: Thiết lập quan hệ bảng thuộc tính, biên tập bản đồ thành quả theo tỷ lệ 1:10.000 và 1:50.000, tạo lập các lớp nhãn (thematic mapping).
  • Thiết bị dã ngoại: GPS Garmin 60CSx, thước đo cao Blumleiss, la bàn địa chất, thước dây xác định ô tiêu chuẩn ($500\text{ m}^2 - 1.000\text{ m}^2$).

Thiết kế cấu trúc cơ sở dữ liệu không gian (Forest Inventory Geodatabase Schema):

-- Cấu trúc bảng thuộc tính dữ liệu biến động tài nguyên rừng
CREATE TABLE Forest_Dynamics_BaoLac (
    FID INT PRIMARY KEY,
    Polygon_Geometry GEOMETRY,
    Xa_ID VARCHAR(50) NOT NULL,
    TieuKhu INT NOT NULL,
    Khoanh INT NOT NULL,
    Lo INT NOT NULL,
    Status_2012_Code VARCHAR(10),      -- Mã trạng thái rừng 2012 (VD: 1100, 1210, 2020)
    Status_2014_Code VARCHAR(10),      -- Mã trạng thái rừng 2014 (VD: 1100, 1220, 2000)
    Forest_Function VARCHAR(20),       -- Chức năng: 'Phong ho', 'San xuat', 'Dac dung'
    Management_Owner VARCHAR(50),      -- Chủ quản lý: 'Ho gia dinh', 'Cong dong', 'To chuc'
    Area_Ha DECIMAL(10, 4),            -- Diện tích tính theo Hecta
    D_bq DECIMAL(5, 2),                -- Đường kính bình quân (cm)
    H_bq DECIMAL(5, 2),                -- Chiều cao bình quân (m)
    Volume_M3_Ha DECIMAL(8, 2),        -- Trữ lượng bình quân (m3/ha)
    Change_Type VARCHAR(50)            -- Phân loại biến động: 'Rừng phục hồi', 'Mất rừng do nương rẫy',...
);

Methodology

Phương pháp luận thực hiện tích hợp quy trình điều tra viễn thám chuẩn của Viện Điều tra Quy hoạch Rừng (VĐTQHR) kết hợp quy chuẩn quản lý kiểm lâm:

  1. Phương pháp kế thừa dữ liệu: Thu thập bộ bản đồ quy hoạch 3 loại rừng năm 2014, số liệu kiểm kê rừng 2010 và bản đồ nền địa hình số hóa tỷ lệ 1:25.000 của Chi cục Kiểm lâm tỉnh Cao Bằng.
  2. Phương pháp điều tra ô tiêu chuẩn (OTC) dã ngoại: Thiết lập các ô tiêu chuẩn đại diện trên 17 xã nhằm thu thập mẫu phổ cho từng trạng thái rừng. Mỗi kiểu thảm thực vật lấy tối thiểu 30 mẫu ngoại nghiệp, đo đếm các chỉ tiêu lâm học cốt lõi ($D_{bq}, H_{bq}$, độ tàn che, tổ thành loài cây mục đích).
  3. Phương pháp phân loại ảnh số có giám định (Supervised Classification): Áp dụng nguyên lý hàm mật độ xác suất cực đại (Maximum Likelihood) nhằm gán từng pixel ảnh vào nhóm trạng thái tương ứng dựa trên vector khoảng cách Mahalanobis trong không gian phổ đa chiều.

Implementation và kết quả

Development process

Quá trình triển khai kỹ thuật được thực hiện tuần tự qua các giai đoạn xử lý số hóa nghiêm ngặt:

Bước 1: Tiền xử lý ảnh viễn thám SPOT 4 trong ERDAS IMAGINE

  • Nhập ảnh gốc định dạng nhị phân thông qua giao diện Import/Generic Binary.
  • Tiến hành ghép dải phổ (Interpreter/Utilities/Layer Stack) nhằm tạo ảnh tổng hợp màu giả (False Color Composite) phục vụ nhận dạng thảm thực vật: Kênh Hồng ngoại gần (NIR), Đỏ (Red) và Xanh lục (Green).
  • Nắn chỉnh tọa độ hình học (Geometric Correction) dựa trên các điểm khống chế mặt đất (Ground Control Points - GCPs) trích xuất từ lớp bản đồ mạng lưới thủy văn (sông, suối) tỷ lệ 1:25.000. Hệ tọa độ chuẩn: UTM WGS 84 North, Zone 48, sau đó chuyển đổi về VN-2000. Thuật toán nội suy mẫu điểm lân cận gần nhất (Nearest Neighbor) được áp dụng nhằm bảo toàn giá trị số phổ (Digital Number - DN) nguyên bản.
  • Cắt ảnh theo đường ranh giới hành chính huyện Bảo Lạc thông qua công cụ Inquire Box và Sub-set nhằm tối ưu hóa dung lượng bộ nhớ và tăng tốc độ xử lý tính toán.

Bước 2: Trích xuất mẫu phổ và phân loại có giám định

  • Sử dụng công cụ Signature Editor kết hợp vùng quan tâm AOI (Area of Interest) để khoanh vẽ mẫu cho từng trạng thái rừng:
    • Rừng tự nhiên gỗ dày/trung bình (LRG): Màu xanh lá cây đậm.
    • Rừng thưa/phục hồi/rừng non (LRN): Màu xanh lá cây nhạt.
    • Đất trống có cây tái sinh (Ic, Ib - T): Màu nâu nhạt.
    • Đất trống không có cây (Ia) / Đất nông nghiệp / Khu dân cư (N): Màu vàng chanh.
    • Vùng bóng núi, mây che phủ (K): Màu đen/trắng sáng.
    • Mặt nước sông suối, hồ (S): Màu xanh nước biển sẫm.
  • Khởi chạy thuật toán Maximum Likelihood với tham số phân loại chuẩn. Sau khi hoàn thành, ma trận điểm ảnh được xử lý qua bộ lọc không gian Neighborhood Majority Filter 3x3 liên tục từ 4 đến 6 lần để khử nhiễu đốm hạt (Salt-and-Pepper noise), gộp các pixel đơn lẻ vào mảng liền kề có diện tích lớn hơn.
+-----------------------------------------------------------------------------------+
|               KỸ THUẬT PHÂN LOẠI ẢNH VÀ CHUYỂN ĐỔI VECTOR GIS                     |
+-----------------------------------------------------------------------------------+
|  1. Raster Classification: Maximum Likelihood -> Majority Filter 3x3 (4-6 lần)   |
|  2. Format Conversion: ERDAS (*.img) -> ESRI Grid (*.grd) -> Shapefile (*.shp)   |
|  3. Spatial Overlay: Intersect (Layer 2012 .shp, Layer 2014 .shp)                |
|  4. Attribute Computing: Area = Area(Shape)/10000 (Ha), SQL Code Mapping          |
+-----------------------------------------------------------------------------------+

Bước 3: Chuyển đổi định dạng và phân tích không gian trong GIS

Xuất tệp sau phân loại từ ERDAS sang định dạng ESRI Grid, sau đó nạp vào ArcView GIS 3.2a / ArcGIS với module Spatial Analyst để chuyển đổi tự động từ Raster Grid sang Vector Shapefile (Convert to Shapefile).

Thực hiện thuật toán giao cắt không gian giữa hai kỳ bản đồ thông qua script tự động hóa xử lý:

# Script Python (ArcPy) tự động hóa tính toán biến động diện tích rừng
import arcpy

# Thiết lập môi trường làm việc
arcpy.env.workspace = "C:/Forestry_Data/BaoLac_GIS"
arcpy.env.overwriteOutput = True

input_feature_2012 = "Forest_Status_2012.shp"
input_feature_2014 = "Forest_Status_2014.shp"
output_intersect = "Forest_Dynamics_2012_2014.shp"

# 1. Thực hiện phép giao cắt không gian giữa 2 kỳ bản đồ
arcpy.analysis.Intersect([input_feature_2012, input_feature_2014], output_intersect, "ALL")

# 2. Thêm trường tính toán diện tích chuẩn Hecta (ha)
arcpy.management.AddField(output_intersect, "Area_Ha", "DOUBLE", 10, 4)
arcpy.management.AddField(output_intersect, "Change_Code", "TEXT", field_length=20)

# 3. Tính toán diện tích từ hình học Topology thực tế
arcpy.management.CalculateField(
    output_intersect, 
    "Area_Ha", 
    "!shape.area@HECTARES!", 
    "PYTHON_9.3"
)

# 4. Xác định mã biến động trạng thái lô rừng
code_block = """
def classify_change(s12, s14):
    if s12 == s14:
        return "Unchanged"
    elif s12.startswith('20') and s14.startswith('11'):
        return "Rung_Phuc_Hoi"
    elif s12.startswith('11') and s14.startswith('20'):
        return "Rung_Suy_Thoai"
    elif s14.startswith('12'):
        return "Rung_Trong_Moi"
    else:
        return "Bien_Dong_Khac"
"""

arcpy.management.CalculateField(
    output_intersect,
    "Change_Code",
    "classify_change(!STATUS_12!, !STATUS_14!)",
    "PYTHON_9.3",
    code_block
)
print("Xử lý giao cắt không gian và tính toán biến động rừng hoàn tất thành công.")

Testing và validation

Công tác kiểm chứng độ chính xác phân loại được tiến hành đối chiếu giữa kết quả giải đoán ảnh SPOT 4 và 120 điểm kiểm tra thực địa bằng định vị vệ tinh GPS trên địa bàn các xã Cốc Pàng, Huy Giáp, Hưng Thịnh, Cô Ba và Xuân Trường.

Trạng thái rừng giải đoán Số mẫu thực địa kiểm chứng Số mẫu phân loại chính xác Độ chính xác người dùng (User's Accuracy) Độ chính xác nhà sản xuất (Producer's Accuracy)
Rừng tự nhiên gỗ (1110) 40 35 $87,5%$ $89,7%$
Rừng núi đá vôi (1150) 25 22 $88,0%$ $84,6%$
Rừng trồng các loại (1200) 20 16 $80,0%$ $76,2%$
Đất trống, trảng cỏ (2020) 20 17 $85,0%$ $89,5%$
Đất có cây tái sinh (2030) 15 12 $80,0%$ $75,0%$
Tổng thể ma trận nhầm lẫn 120 102 Độ chính xác tổng thể: 85,0% Chỉ số Kappa: 0,81

Chỉ số Kappa đạt $0,81$ chứng minh phương pháp giải đoán ảnh viễn thám kết hợp phân loại giám định Maximum Likelihood và lọc nhiễu không gian đáp ứng đầy đủ tiêu chuẩn kỹ thuật điều tra quy hoạch lâm nghiệp cấp quốc gia.

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

1. Hiện trạng tài nguyên rừng và đất lâm nghiệp huyện Bảo Lạc năm 2014

Tổng diện tích tự nhiên toàn huyện đạt $92.063,68\text{ ha}$, cơ cấu sử dụng đất và hiện trạng rừng năm 2014 được xác định cụ thể:

  • Đất có rừng: Đạt $45.121,23\text{ ha}$, chiếm tỷ lệ $49,01%$ diện tích tự nhiên.
    • Rừng tự nhiên: Chiếm ưu thế tuyệt đối với $43.661,01\text{ ha}$ ($47,42%$ diện tích tự nhiên). Trong đó, rừng gỗ núi đất chiếm $35.337,62\text{ ha}$ ($38,38%$, trữ lượng bình quân $135\text{ m}^3\text{/ha}$, $H_{bq}=11\text{ m}$, $D_{bq}=26\text{ cm}$, độ tàn che $0,4 - 0,6$); Rừng trên núi đá vôi chiếm $8.223,39\text{ ha}$ ($8,93%$, trữ lượng bình quân $32\text{ m}^3\text{/ha}$, $H_{bq}=5,5\text{ m}$, $D_{bq}=12\text{ cm}$).
    • Rừng trồng: Tổng diện tích $1.460,22\text{ ha}$ ($1,59%$ diện tích tự nhiên), gồm rừng trồng gỗ có trữ lượng ($80,67\text{ ha}$, trữ lượng $50\text{ m}^3\text{/ha}$ gồm Thông, Sa mộc, Mỡ, Keo); rừng trồng chưa có trữ lượng ($82,99\text{ ha}$); rừng tre luồng/trúc ($1.199,53\text{ ha}$); rừng cây đặc sản Hồi, Quế, Dẻ ($97,03\text{ ha}$).
  • Đất chưa có rừng: $37.605,22\text{ ha}$ ($40,84%$ diện tích tự nhiên). Trong đó:
    • Đất trống không có gỗ tái sinh ($I_a, I_b$): $24.629,47\text{ ha}$ ($65,49%$ diện tích đất chưa có rừng).
    • Đất có cây gỗ tái sinh ($I_c$): $6.425,11\text{ ha}$ ($17,08%$ diện tích đất chưa có rừng, mật độ tái sinh đạt $500 - 1.000\text{ cây/ha}$, tỷ lệ cây mục đích chiếm $55 - 65%$).
    • Núi đá không có rừng: $6.521,93\text{ ha}$ ($17,34%$).
    • Nương rẫy trong lâm nghiệp: $1,71\text{ ha}$; đất khác trong lâm nghiệp: $27\text{ ha}$.
  • Đất ngoài lâm nghiệp (Nông nghiệp, đất ở thổ cư, sông suối): $9.337,23\text{ ha}$ ($10,15%$).
+-----------------------------------------------------------------------------------+
|            CƠ CẤU HIỆN TRẠNG TÀI NGUYÊN RỪNG HUYỆN BẢO LẠC NĂM 2014               |
+-----------------------------------------------------------------------------------+
|  [Đất có rừng: 45.121,23 ha (49,01%)]                                             |
|    |-- Rừng tự nhiên: 43.661,01 ha (47,42%) -> Rừng gỗ (35.337 ha) + Núi đá (8.223 ha)|
|    \-- Rừng trồng:    1.460,22 ha ( 1,59%) -> Tre nứa (1.199 ha) + Gỗ/Đặc sản      |
|  [Đất chưa có rừng: 37.605,22 ha (40,84%)]                                        |
|    |-- Đất trống không cây tái sinh (Ia, Ib): 24.629,47 ha                        |
|    |-- Đất có cây gỗ tái sinh mục đích (Ic):   6.425,11 ha                        |
|    \-- Núi đá trọc không có rừng:              6.521,93 ha                        |
|  [Đất khác ngoài lâm nghiệp: 9.337,23 ha (10,15%)]                                |
+-----------------------------------------------------------------------------------+

2. Phân bố diện tích theo chức năng 3 loại rừng và chủ thể quản lý

  • Theo chức năng quy hoạch: Toàn huyện không có diện tích quy hoạch rừng đặc dụng. Tổng diện tích lâm phận phòng hộ đạt $41.830,25\text{ ha}$ (trong đó đất có rừng phòng hộ là $25.270,33\text{ ha}$). Tổng diện tích lâm phận sản xuất đạt $40.862,61\text{ ha}$ (trong đó đất có rừng sản xuất là $19.848,31\text{ ha}$).
  • Theo chủ quản lý: Hộ gia đình cá nhân quản lý diện tích lớn nhất với $42.865,25\text{ ha}$ ($23.711,10\text{ ha}$ đất có rừng và $16.700,78\text{ ha}$ đất chưa có rừng); Cộng đồng dân cư thôn bản quản lý $26.127,45\text{ ha}$ ($11.288,05\text{ ha}$ rừng tự nhiên); Các tổ chức khác quản lý $15,31\text{ ha}$.

3. Phân tích biến động diện tích rừng giai đoạn 2012–2014

Kết quả phân tích không gian giữa bản đồ 2012 và 2014 cho thấy tổng diện tích đất có rừng toàn huyện tăng $120,13\text{ ha}$ (bình quân tăng $40,04\text{ ha/năm}$).

Phân loại đất / loại rừng Diện tích 2012 (ha) Diện tích 2014 (ha) Biến động (ha) Tỷ lệ biến động (%) Tốc độ BQ (ha/năm)
A. Rừng tự nhiên 43.450,12 43.661,01 +210,89 +0,49% +70,30
- Rừng gỗ tự nhiên 35.293,50 35.337,62 +44,12 +0,13% +14,71
- Rừng trên núi đá vôi 8.156,62 8.323,39 +166,77 +2,04% +55,59
B. Rừng trồng 1.530,94 1.460,22 -70,72 -4,62% -23,57
- Rừng trồng có trữ lượng 170,12 80,67 -89,45 -52,58% -29,82
- Rừng trồng chưa trữ lượng 120,43 82,99 -37,44 -31,09% -12,48
- Rừng tre luồng, trúc sào 1.152,32 1.199,53 +47,21 +4,10% +15,74
- Rừng cây đặc sản 108,11 97,03 -11,08 -10,25% -3,69
C. Đất chưa có rừng 37.590,87 37.605,22 +14,35 +0,04% +4,78
- Đất không cây tái sinh ($I_a, I_b$) 24.650,66 24.629,47 -21,19 -0,09% -7,06
D. Đất khác ngoài lâm nghiệp 9.471,71 9.337,23 -134,48 -1,42% -44,83
Tổng diện tích tự nhiên 92.063,68 92.063,68 0,00 0,00% 0,00 
+-----------------------------------------------------------------------------------+
|               BIẾN ĐỘNG RỪNG THEO ĐƠN VỊ HÀNH CHÍNH (2012 - 2014)                 |
+-----------------------------------------------------------------------------------+
|  Top tăng rừng tự nhiên:                                                          |
|    1. Xã Hưng Thịnh : +271,08 ha (BQ +90,40 ha/năm) -> Khoanh nuôi tái sinh tốt   |
|    2. Xã Cô Ba      : +181,44 ha (BQ +60,50 ha/năm)                              |
|    3. Xã Hưng Đạo   : +113,54 ha (BQ +37,80 ha/năm)                              |
|  Top giảm rừng tự nhiên:                                                          |
|    1. Xã Cốc Pàng   : -213,99 ha (BQ -71,30 ha/năm) -> Cháy rừng, áp lực nương rẫy|
|    2. Xã Phan Thanh : -128,38 ha (BQ -42,80 ha/năm)                              |
|    3. Xã Hồng An    : -75,57 ha  (BQ -25,20 ha/năm)                              |
+-----------------------------------------------------------------------------------+
  • Phân tích nguyên nhân biến động:
    • Tích cực: Diện tích rừng tự nhiên toàn huyện tính theo tổng các xã tăng ròng $524,33\text{ ha}$ nhờ hiệu quả vốn đầu tư từ Dự án Trồng mới 5 triệu ha rừng (Dự án 661), Chương trình 135, Dự án Xóa bỏ cây thuốc phiện (DA06/CP) và Dự án Phát triển Nông thôn Cao Bằng - Bắc Kạn do Liên minh Châu Âu (EU) tài trợ. Ý thức bảo vệ rừng cộng đồng thôn bản và việc giao đất giao rừng cho hộ gia đình được đẩy mạnh.
    • Tiêu cực: Rừng trồng gỗ giảm mạnh $89,45\text{ ha}$ do người dân khai thác trắng khi đến tuổi thành thục kinh tế nhưng chưa kịp trồng lại rừng sau khai thác. Tại các xã vùng sâu biên giới như Cốc Pàng (giảm $213,99\text{ ha}$ rừng tự nhiên), Phan Thanh (giảm $128,38\text{ ha}$), áp lực xâm lấn đốt rừng làm nương rẫy và nguy cơ cháy rừng trong mùa hanh khô vẫn diễn biến phức tạp.

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

Nghiên cứu mang lại những cải tiến mang tính bước ngoặt cho công tác quản lý tài nguyên lâm nghiệp tại địa phương:

+-----------------------------------------------------------------------------------+
|                            CÁC ĐÓNG GÓP KỸ THUẬT NỔI BẬT                          |
+-----------------------------------------------------------------------------------+
|  1. Chuẩn hóa quy trình 3S khép kín: Tiền xử lý -> Phân loại -> Vector -> GIS     |
|  2. Tối ưu thời gian: Rút ngắn 75% chu kỳ cập nhật CSDL so với phương pháp cũ    |
|  3. Xây dựng bản đồ số 100% trên nền Hệ quy chiếu quốc gia VN-2000                 |
|  4. Định vị chính xác các điểm nóng suy thoái rừng cấp lô/khoảnh/tiểu khu        |
+-----------------------------------------------------------------------------------+
  1. Chuẩn hóa quy trình tích hợp 3S khép kín: Khắc phục triệt để tình trạng phân mảnh dữ liệu giữa đo đạc dã ngoại và tổng hợp số liệu. Quy trình liên hoàn từ xử lý ảnh số SPOT 4 trong ERDAS đến phân tích không gian trong ArcView và xuất bản đồ chuyên đề trong MapInfo tạo nên một chuẩn kỹ thuật đồng nhất.
  2. Nâng cao hiệu suất và độ chính xác không gian: Thay vì dựa vào ước tính thủ công, việc ứng dụng thuật toán Maximum Likelihood kết hợp bộ lọc không gian Majority Filter 3x3 và nắn chỉnh hình học theo mạng thủy văn giúp nâng độ chính xác tổng thể lên $85,0%$, giảm thiểu sai số phân loại biên lô rừng trên địa hình núi đá phức tạp.
  3. Thành lập hệ thống CSDL không gian đồng bộ: Toàn bộ dữ liệu của $92.063,68\text{ ha}$ đất đai thuộc 17 xã/thị trấn được chuyển đổi hoàn chỉnh sang Hệ quy chiếu quốc gia VN-2000, cung cấp hệ thống bản đồ thành quả tỷ lệ 1:50.000 (cấp huyện) và 1:10.000 (cấp xã) phục vụ trực tiếp cho các kỳ quy hoạch sử dụng đất và kiểm kê rừng tiếp theo.

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

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

Hệ thống cơ sở dữ liệu và bản đồ biến động rừng được triển khai ứng dụng trực tiếp tại các cơ quan quản lý:

  • Hạt Kiểm lâm huyện Bảo Lạc: Sử dụng bản đồ biến động 1:10.000 để xác định chính xác tọa độ các vùng rừng bị suy giảm nghiêm trọng tại xã Cốc Pàng (tiểu khu 3) và xã Phan Thanh (tiểu khu 48) nhằm khoanh vùng tuần tra trọng điểm, ngăn chặn kịp thời các hành vi phá rừng làm nương rẫy.
  • Ủy ban nhân dân huyện Bảo Lạc: Căn cứ vào số liệu $24.629,47\text{ ha}$ đất trống đồi trọc ($I_a, I_b$) và $6.425,11\text{ ha}$ đất có cây tái sinh ($I_c$) để phân bổ nguồn vốn ngân sách trồng rừng kinh tế, khoanh nuôi xúc tiến tái sinh tự nhiên trong kế hoạch phát triển kinh tế - xã hội hàng năm.
  • Chi cục Kiểm lâm tỉnh Cao Bằng: Tích hợp dữ liệu số huyện Bảo Lạc vào hệ thống thông tin quản lý rừng (VFRS) toàn tỉnh, phục vụ công tác chi trả dịch vụ môi trường rừng (PES) chính xác cho các chủ rừng là hộ gia đình và cộng đồng dân cư thôn bản.
+-----------------------------------------------------------------------------------+
|                        LỘ TRÌNH VÀ YÊU CẦU TRIỂN KHAI HỆ THỐNG                     |
+-----------------------------------------------------------------------------------+
|  Phần cứng: Máy trạm Workstation (Core i5/i7, 8GB+ RAM, GPU 2GB+, GPS Garmin)     |
|  Phần mềm : ERDAS IMAGINE 9.2+, ArcGIS Desktop 10.x, MapInfo Professional 10.0+   |
|  Lộ trình : Giai đoạn 1 (Nhập liệu/Nắn ảnh) -> Giai đoạn 2 (Phân loại/Thực địa)   |
|             -> Giai đoạn 3 (Chồng xếp không gian) -> Giai đoạn 4 (Bàn giao CSDL)  |
+-----------------------------------------------------------------------------------+

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

Hạn chế kỹ thuật hiện tại

  • Đặc trưng phổ ảnh quang học: Tư liệu ảnh SPOT 4 có độ phân giải không gian $20\text{ m}$, do đó chưa bóc tách triệt để các diện tích rừng trồng non phân tán, trồng xen dưới tán rừng tự nhiên hoặc các mô hình vườn đồi nhỏ lẻ ($<0,1\text{ ha}$).
  • Nhiễu khí tượng và địa hình: Khí hậu cận nhiệt đới ẩm vùng cao thường xuyên xuất hiện mây bao phủ và hiện tượng bóng núi đá vôi vào mùa đông, đòi hỏi phải loại trừ thủ công bằng mặt nạ mây và kiểm chứng thực địa bổ sung.
  • Độ trễ cập nhật tư liệu ảnh: Chu kỳ thu thập ảnh vệ tinh thương mại định kỳ phụ thuộc vào kinh phí cấp phát theo từng giai đoạn dự án.

Hướng nâng cấp và mở rộng

  1. Tích hợp tư liệu viễn thám đa nguồn miễn phí độ phân giải cao: Chuyển dịch sang khai thác ảnh Sentinel-2 ($10\text{ m}$) và Landsat 8/9 Operational Land Imager (OLI), kết hợp dữ liệu siêu cao từ PlanetScope ($3\text{ m}$) để nâng cao độ phân giải không gian và rút ngắn chu kỳ quan trắc xuống $5 - 10\text{ ngày/lần}$.
  2. Ứng dụng Trí tuệ nhân tạo và Điện toán đám mây: Xây dựng mô hình phân loại tự động sử dụng thuật toán Random Forest, Support Vector Machine (SVM) hoặc mạng nơ-ron tích chập (Convolutional Neural Networks - CNN) trên nền tảng Google Earth Engine (GEE), cho phép xử lý dữ liệu lớn (Big Geospatial Data) theo thời gian thực.
  3. Kết hợp thiết bị bay không người lái (UAV/Drone): Trang bị UAV tích hợp cảm biến đa phổ (Multispectral) phục vụ đo đạc kiểm chứng nhanh các điểm nóng cháy rừng và kiểm kê trữ lượng lâm phần ở cấp độ từng cây.

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

+-----------------------------------------------------------------------------------+
|                        GIÁ TRỊ MANG LẠI CHO CÁC NHÓM ĐỐI TƯỢNG                    |
+-----------------------------------------------------------------------------------+
|  [Sinh viên & Học viên cao học] -> Giáo trình thực hành chuẩn về xử lý ảnh 3S    |
|  [Kỹ sư GIS & Lâm nghiệp]      -> Bộ script mẫu, quy trình nắn chỉnh & vector hóa |
|  [Cơ quan Quản lý Nhà nước]   -> CSDL chuẩn hóa VN-2000, tối ưu hóa ngân sách    |
|  [Nhà nghiên cứu Sinh thái]    -> Chuỗi số liệu định lượng về thảm phủ karst      |
+-----------------------------------------------------------------------------------+
  • Sinh viên và Học viên ngành Lâm nghiệp/Quản lý Đất đai/GIS: Tiếp cận tài liệu tham khảo hoàn chỉnh về phương pháp luận kết hợp giữa lâm học cổ điển và công nghệ địa tin học hiện đại; nắm vững quy trình xử lý dữ liệu viễn thám thực tế.
  • Kỹ sư GIS và Cán bộ kỹ thuật lâm nghiệp: Ứng dụng trực tiếp quy trình nắn chỉnh ảnh, ma trận phân loại Maximum Likelihood và bộ code tự động hóa giao cắt không gian để triển khai cho các địa bàn miền núi khác.
  • Cơ quan quản lý Nhà nước (Kiểm lâm, Phòng Nông nghiệp & PTNT, UBND các cấp): Sở hữu công cụ số hóa trực quan, chính xác để chỉ đạo sản xuất lâm nghiệp, lập quy hoạch 3 loại rừng và quản lý chi trả dịch vụ môi trường rừng minh bạch.
  • Các nhà nghiên cứu môi trường và biến đổi khí hậu: Khai thác bộ số liệu thực nghiệm tin cậy về sinh khối, trữ lượng rừng gỗ và đặc tính phục hồi của hệ sinh thái rừng núi đá vôi vùng Đông Bắc.

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 để vận hành hệ thống CSDL này là gì?

Hệ thống vận hành mượt mà trên máy tính để bàn hoặc laptop trang bị vi xử lý Intel Core i5/i7 (tốc độ $\ge 2.5\text{ GHz}$), RAM tối thiểu $8\text{ GB}$ (khuyến nghị $16\text{ GB}$ cho xử lý ảnh số lớn), ổ cứng SSD dung lượng trống tối thiểu $50\text{ GB}$ và card đồ họa rời. Về phần mềm, cần cài đặt ERDAS IMAGINE (phiên bản 9.2 hoặc mới hơn), ArcGIS Desktop 10.x / ArcView GIS 3.2a và MapInfo Professional 10.0+ trên hệ điều hành Windows.

2. Làm thế nào để loại bỏ sai số do bóng núi đá vôi và mây che phủ trên ảnh viễn thám?

Trong quy trình tiền xử lý tại ERDAS, vùng mây sáng và bóng núi được khoanh vẽ thành các lớp đối tượng riêng biệt (Signature class K1, K2). Kỹ thuật viên kết hợp sử dụng mô hình số hóa độ cao (DEM) để phân tích góc chiếu sáng mặt trời và tiến hành đối chiếu với dữ liệu kiểm chứng ngoại nghiệp bằng máy định vị GPS Garmin tại thực địa để gán lại mã trạng thái rừng chính xác trước khi xuất kết quả vector.

3. Quy trình chuyển đổi dữ liệu tọa độ từ WGS-84 sang VN-2000 được thực hiện ra sao?

Dữ liệu ảnh viễn thám gốc sau khi nắn chỉnh hình học trên phép chiếu UTM WGS-84 Zone 48N sẽ được chuyển đổi sang Hệ quy chiếu quốc gia VN-2000 bằng phần mềm chuyển đổi tọa độ chuẩn của Bộ Tài nguyên & Môi trường ban hành (theo Thông tư 973/2001/TT-TCĐC), với kinh tuyến trục và độ lệch giả lập phù hợp cho tỉnh Cao Bằng nhằm đảm bảo tính thống nhất pháp lý khi in ấn bản đồ thành quả.

4. Tại sao diện tích rừng trồng có trữ lượng lại có xu hướng giảm trong giai đoạn 2012–2014?

Diện tích rừng trồng có trữ lượng giảm $89,45\text{ ha}$ (từ $170,12\text{ ha}$ năm 2012 xuống $80,67\text{ ha}$ năm 2014) chủ yếu do chu kỳ khai thác gỗ rừng trồng (Keo, Mỡ, Thông) của các hộ gia đình đến tuổi thành thục kinh tế. Sau khi khai thác trắng, các diện tích này chưa kịp khép tán hoặc đang trong giai đoạn chuẩn bị đất tái trồng rừng, do đó trên ảnh viễn thám được xếp vào nhóm trạng thái đất chưa có rừng ($I_a, I_b$) hoặc rừng trồng chưa có trữ lượng.

5. Khả năng tích hợp CSDL này vào Hệ thống theo dõi diễn biến rừng Quốc gia (VFRS) như thế nào?

Toàn bộ lớp dữ liệu không gian định dạng Shapefile (.shp) và bảng thuộc tính tương thích hoàn toàn với cấu trúc cơ sở dữ liệu địa lý (Geodatabase) của Cục Kiểm lâm. Dữ liệu có thể nạp trực tiếp vào phần mềm VFRS thông qua các công cụ nhập liệu chuẩn, giúp đồng bộ hóa dữ liệu từ cấp huyện lên cấp tỉnh và Trung ương mà không cần thực hiện tái cấu trúc bảng.

Kết luận

Đề tài "Đánh giá tình hình diễn biến diện tích và trạng thái rừng bằng phần mềm GIS tại huyện Bảo Lạc, tỉnh Cao Bằng giai đoạn 2012–2014" đã giải quyết trọn vẹn bài toán định lượng và không gian hóa tài nguyên lâm nghiệp vùng cao Đông Bắc. Thông qua việc làm chủ công nghệ viễn thám SPOT 4, thuật toán phân loại có giám định Maximum Likelihood và kỹ thuật phân tích không gian GIS, nghiên cứu đã xây dựng thành công bộ cơ sở dữ liệu số toàn diện cho $92.063,68\text{ ha}$ diện tích tự nhiên toàn huyện.

Kết quả khẳng định tính ưu việt vượt trội của công nghệ 3S trong việc rút ngắn thời gian điều tra, giảm chi phí nhân lực và nâng cao độ chính xác tổng thể lên $85,0%$. Dữ liệu chỉ rõ diện tích đất có rừng năm 2014 đạt $45.121,23\text{ ha}$ ($49,01%$), ghi nhận mức tăng ròng $120,13\text{ ha}$ so với năm 2012, đồng thời xác lập vị trí các điểm nóng suy thoái rừng tại Cốc Pàng và Phan Thanh. Đây là nguồn luận cứ khoa học và công cụ thực thi vô cùng giá trị, phục vụ trực tiếp cho các nhà quản lý, kỹ sư lâm nghiệp và chính quyền địa phương trong chiến lược phát triển bền vững tài nguyên rừng và ứng phó với biến đổi khí hậu trong tương lai.