Tổng quan nghiên cứu

Vườn Quốc gia Ba Bể tại tỉnh Bắc Kạn là một trong những khu dự trữ sinh quyển trọng yếu của Việt Nam, sở hữu hệ sinh thái tự nhiên phong phú với khoảng 1.281 loài thực vật thuộc 162 họ và 672 chi. Nhiều loài trong số đó được ghi nhận trong Sách Đỏ Việt Nam và thế giới, đòi hỏi một cơ chế giám sát nghiêm ngặt và khoa học. Tuy nhiên, công tác điều tra và thành lập bản đồ hiện trạng rừng tại đây bằng các phương pháp đo vẽ trắc địa truyền thống gặp rất nhiều rào cản do địa hình núi đá vôi hiểm trở, độ dốc lớn và diện tích trải rộng qua nhiều xã vùng đệm. Hạn chế này dẫn đến việc thiếu hụt cơ sở dữ liệu không gian đồng bộ và làm chậm trễ các quyết định quy hoạch bảo tồn lâm nghiệp.

Nhằm giải quyết triệt để vấn đề trên, luận văn tập trung nghiên cứu ứng dụng công nghệ Hệ thống thông tin địa lý (GIS) và Viễn thám (Remote Sensing) để xây dựng bản đồ trạng thái rừng chuẩn xác tại khu vực Vườn Quốc gia Ba Bể. Mục tiêu cụ thể của nghiên cứu là thu thập, tiền xử lý và giải đoán ảnh vệ tinh SPOT-5 chụp vào tháng 4 năm 2009 với độ phân giải cao 2.5m, từ đó thiết lập cơ sở dữ liệu không gian kết hợp dữ liệu thuộc tính về trạng thái tài nguyên rừng. Nghiên cứu giới hạn phạm vi không gian tại vùng lõi vườn quốc gia và các xã vùng đệm liên quan thuộc huyện Ba Bể. Về mặt ý nghĩa thực tiễn, việc ứng dụng giải pháp công nghệ này giúp giảm hơn 60% chi phí và thời gian điều tra thực địa, đồng thời nâng cao độ chính xác trong công tác quản lý tài nguyên lên trên 85%, cung cấp luận cứ khoa học vững chắc phục vụ chiến lược quy hoạch lâm nghiệp bền vững tại địa phương.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Luận văn được xây dựng trên nền tảng lý thuyết phân loại rừng theo Thông tư số 34/2009/TT-BNNPTNT của Bộ Nông nghiệp & Phát triển nông thôn và tiêu chuẩn phân loại trạng thái rừng của Loeschau năm 1966. Khung phân loại này chia thảm thực vật thành 4 nhóm chính: Nhóm I là đất chưa có rừng (Ia, Ib, Ic), Nhóm II là rừng phục hồi (IIa, IIb), Nhóm III là rừng thứ sinh phân theo cấp trữ lượng từ 50 m3/ha đến trên 300 m3/ha (IIIa1, IIIa2, IIIa3, IIIb, IIIc), và Nhóm IV là rừng nguyên sinh giàu trữ lượng.

Bên cạnh đó, nghiên cứu vận dụng lý thuyết bức xạ điện từ và viễn thám quang học để đánh giá phản xạ phổ của thảm thực vật thông qua chỉ số thực vật NDVI (Normalized Difference Vegetation Index) với dải giá trị chuẩn từ -1 đến +1. Trong mô hình GIS, dữ liệu được tổ chức dưới hai cấu trúc cốt lõi: dữ liệu không gian Raster (mạng lưới ô vuông pixel biểu thị mức độ phản xạ phổ) và dữ liệu Vector (các đối tượng hình học Point, Line, Polygon). Cơ chế liên kết dữ liệu không gian với bảng thuộc tính phi không gian thông qua cấu trúc dữ liệu Shapefile đóng vai trò trung tâm trong phân tích địa lý.

Phương pháp nghiên cứu

Nguồn dữ liệu chính của nghiên cứu là cảnh ảnh vệ tinh SPOT-5 mã hiệu F-48-44-A với dung lượng 371 MB, độ phân giải không gian 2.5m được chụp vào tháng 4 năm 2009, kết hợp với các lớp dữ liệu thủy văn và bản đồ hiện trạng lâm nghiệp thời kỳ trước. Cỡ mẫu nghiên cứu bao gồm 15 điểm khống chế mặt đất (Ground Control Points - GCP) được lựa chọn phân bố đồng đều trên toàn bộ cảnh ảnh phục vụ nắn chỉnh hình học, cùng với 100 điểm kiểm chứng ngoại nghiệp phân tầng ngẫu nhiên theo các dạng thực phủ trên thực địa để đánh giá độ chính xác phân loại.

Về phương pháp phân tích, nghiên cứu ứng dụng phép chiếu UTM trên hệ quy chiếu WGS 84 để nắn chỉnh không gian hình học, khống chế sai số toàn phương RMS ở mức dưới 1 pixel (nhỏ hơn 2.5m). Quá trình giải đoán ảnh sử dụng phương pháp phân loại có chọn mẫu với thuật toán khoảng cách nhỏ nhất tới số trung bình (Minimum Distance to Mean) trên phần mềm ENVI 4.5. Lý do lựa chọn thuật toán này là khả năng tính toán nhanh khoảng cách Euclidean từ mỗi pixel tới vector trung bình của từng vùng mẫu quan tâm (ROI), tối ưu hóa độ phân tách giữa các trạng thái rừng có phổ tương đồng trên địa hình đồi núi phức tạp. Kết quả giải đoán sau đó được vector hóa và biên tập trên hệ thống ArcGIS 9.2. Timeline nghiên cứu được thực hiện đồng bộ từ khâu thu thập ảnh năm 2009 đến hoàn thiện biên tập bản đồ và thử nghiệm hệ thống web bản đồ vào năm 2011.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Nghiên cứu đã mang lại 4 phát hiện quan trọng về cấu trúc hiện trạng tài nguyên rừng tại khu vực Vườn Quốc gia Ba Bể:

Thứ nhất, kết quả phân loại đã bóc tách rõ rệt 5 trạng thái rừng và thảm phủ chủ đạo, trong đó rừng tự nhiên trên núi đá vôi (thuộc nhóm IIb và III) chiếm tỷ lệ áp đảo với hơn 65% tổng diện tích tự nhiên của vùng lõi vườn quốc gia.

Thứ hai, việc nắn chỉnh ảnh với 15 điểm GCP và lọc tăng cường độ tương phản trên phần mềm ENVI 4.5 giúp nâng cao khả năng phân biệt giữa rừng phục hồi sau nương rẫy (trạng thái IIa) và các trảng cây bụi đất trống (trạng thái Ib), đạt độ chính xác giải đoán tổng thể trên 85% khi đối chiếu với 100 điểm kiểm định dã ngoại.

Thứ ba, quy trình chuyển đổi từ cấu trúc Raster sang Vector dạng Shapefile đã chuẩn hóa thành công 5 lớp dữ liệu chuyên đề gồm: lớp đất rừng, mạng lưới giao thông, hệ thống thủy văn sông hồ, mô hình số độ cao DEM và lớp ranh giới dân cư - đất nông nghiệp.

Thứ tư, nghiên cứu phát hiện sự suy thoái cục bộ của diện tích rừng tự nhiên tại khu vực ranh giới giáp ranh giữa vùng đệm và vùng lõi thuộc địa bàn 8 xã liên quan, nơi có hoạt động canh tác nông nghiệp chiếm khoảng 15% đến 20% quỹ đất ven hồ.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự phân bố manh mún của một số dải rừng phục hồi là do áp lực dân sinh và tập quán canh tác nương rẫy kéo dài của cộng đồng địa phương trước khi vườn quốc gia được khoanh vùng bảo vệ nghiêm ngặt. Đặc thù địa hình núi đá dốc đứng tạo ra hiện tượng bóng che cục bộ trên ảnh viễn thám, nhưng việc kết hợp chỉ số thực vật NDVI và phương pháp phân lớp ROI đã bù trừ hiệu quả sai số quang học này.

Trong thực tế, dữ liệu trạng thái rừng có thể được trực quan hóa tối ưu thông qua biểu đồ tròn mô tả cơ cấu tỷ lệ diện tích các trạng thái rừng và bảng ma trận nhầm lẫn (Confusion Matrix) chi tiết hóa sai số giữa kết quả giải đoán và thực địa. So với các nghiên cứu trước đây chỉ sử dụng ảnh Landsat độ phân giải 30m với sai số biên vẽ lên tới 10% đến 15%, việc ứng dụng ảnh SPOT-5 độ phân giải 2.5m trong luận văn này đã tăng độ chi tiết không gian gấp 12 lần, cho phép phân định chính xác các khoảnh rừng có diện tích nhỏ dưới 0.1 ha.

Đề xuất và khuyến nghị

Dựa trên kết quả nghiên cứu, 4 giải pháp hành động cụ thể được đề xuất nhằm nâng cao năng lực quản lý rừng:

Một là, chuẩn hóa và cập nhật định kỳ cơ sở dữ liệu ảnh viễn thám SPOT với chu kỳ 12 tháng một lần, nhằm duy trì mục tiêu kiểm soát 100% các biến động thảm phủ rừng tại vùng lõi và vùng đệm. Chủ thể thực hiện là Ban Quản lý Vườn Quốc gia Ba Bể phối hợp cùng Chi cục Kiểm lâm tỉnh Bắc Kạn, triển khai từ quý I năm 2012.

Hai là, triển khai tích hợp bản đồ số trạng thái rừng lên nền tảng Web GIS ứng dụng ngôn ngữ JavaScript và PHP, nhằm cắt giảm 40% thời gian tra cứu và tổng hợp báo cáo chuyên đề cho cán bộ kiểm lâm. Chủ thể thực hiện là Trung tâm Công nghệ Thông tin trực thuộc Sở Tài nguyên và Môi trường, hoàn thành trong lộ trình 6 tháng.

Ba là, thiết lập 50 trạm giám sát cố định sử dụng thiết bị định vị vệ tinh GPS tại các khu vực ranh giới xung yếu giữa đất rừng đặc dụng và đất nông nghiệp của 8 xã vùng đệm, phấn đấu kéo giảm tỷ lệ xâm lấn rừng tự nhiên xuống dưới 1.5% mỗi năm. Chủ thể thực hiện là Hạt Kiểm lâm Ba Bể, thực hiện liên tục giai đoạn 2012 - 2015.

Bốn là, tổ chức các khóa đào tạo nâng cao kỹ năng xử lý phần mềm ArcGIS 9.2 và ENVI 4.5 cho 100% cán bộ kỹ thuật và kiểm lâm địa bàn trong thời gian 3 tháng, đảm bảo 100% nhân sự làm chủ công nghệ giải đoán ảnh số và quản trị cơ sở dữ liệu không gian. Chủ thể thực hiện là Ủy ban nhân dân huyện Ba Bể phối hợp với các cơ sở đào tạo chuyên ngành.

Đối tượng nên tham khảo luận văn

Luận văn là tài liệu tham khảo có giá trị chuyên môn cao cho 4 nhóm đối tượng cụ thể:

Nhóm thứ nhất là ban quản lý các khu bảo tồn thiên nhiên và vườn quốc gia trên toàn quốc, có thể tiếp cận quy trình ứng dụng ảnh viễn thám độ phân giải cao 2.5m để giám sát thảm phủ và phân định ranh giới bảo tồn 100% diện tích rừng đặc dụng.

Nhóm thứ hai là các kỹ sư GIS, chuyên viên đo đạc viễn thám và chuyên gia tin học địa lý, có thể tham khảo phương pháp xử lý ảnh số trên ENVI 4.5, quy trình nắn chỉnh 15 điểm GCP đạt sai số RMS dưới 1 pixel và giải pháp vector hóa dữ liệu raster sang shapefile.

Nhóm thứ ba là học viên cao học, nghiên cứu sinh và giảng viên chuyên ngành Hệ thống thông tin, Lâm nghiệp, Quản lý tài nguyên môi trường, sử dụng luận văn làm tài liệu nghiên cứu điển hình về tích hợp công nghệ không gian và xây dựng web bản đồ trực quan bằng ngôn ngữ JavaScript.

Nhóm thứ tư là các cơ quan quản lý nhà nước và chính quyền địa phương cấp huyện, xã, có thể khai thác bộ số liệu thống kê hiện trạng diện tích rừng tại 8 xã vùng đệm để xây dựng quy hoạch kinh tế - xã hội bền vững trong tầm nhìn 5 đến 10 năm tới.

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

Ưu điểm vượt trội của ảnh vệ tinh SPOT-5 độ phân giải 2.5m so với các nguồn ảnh vệ tinh khác là gì? Ảnh vệ tinh SPOT-5 với độ phân giải không gian 2.5m và dung lượng 371 MB cung cấp độ chi tiết quang học cao gấp 12 lần so với ảnh Landsat 30m truyền thống. Điều này cho phép nhận diện chính xác các lô rừng nhỏ dưới 0.1 ha và phân tách rõ nét các trạng thái thực phủ phức tạp trên địa hình núi đá vôi dốc đứng.

Quy trình nắn chỉnh không gian ảnh vệ tinh trong nghiên cứu được kiểm soát độ chính xác như thế nào? Nghiên cứu sử dụng 15 điểm khống chế mặt đất (GCP) trải đều trên toàn bộ cảnh ảnh F-48-44-A theo hệ quy chiếu UTM WGS 84. Sai số toàn phương trung bình RMS được kiểm soát nghiêm ngặt ở mức dưới 1 pixel (tương đương dưới 2.5m thực tế), đảm bảo độ tin cậy hình học cao nhất cho bản đồ.

Thuật toán Minimum Distance to Mean hoạt động ra sao trong việc giải đoán trạng thái rừng? Thuật toán tính toán khoảng cách Euclidean từ mỗi điểm ảnh chưa xác định tới vector trung bình của các vùng mẫu ROI đã được khảo sát dã ngoại. Phương pháp này tối ưu hóa tốc độ xử lý trên phần mềm ENVI 4.5 và phân loại thành công 100% các pixel thành 5 nhóm trạng thái rừng chính.

Luận văn căn cứ vào những quy chuẩn pháp lý và kỹ thuật nào để phân loại trạng thái rừng? Luận văn kết hợp khung pháp lý của Thông tư số 34/2009/TT-BNNPTNT ngày 10/6/2009 của Bộ Nông nghiệp & Phát triển nông thôn và hệ thống phân loại Loeschau năm 1966. Tiêu chuẩn này phân cấp khoa học các trạng thái từ đất trống Ia, Ib đến rừng thứ sinh IIIa, IIIb với trữ lượng từ 50 đến trên 300 m3/ha.

Hệ thống Web GIS được xây dựng bằng công nghệ nào và mang lại lợi ích gì? Hệ thống được phát triển bằng ngôn ngữ lập trình PHP kết hợp JavaScript, cho phép xuất bản trực tuyến 5 lớp bản đồ chuyên đề động. Người dùng có thể dễ dàng tra cứu, truy vấn thông tin không gian và thuộc tính của hơn 1.281 loài thực vật mà không cần cài đặt các phần mềm GIS chuyên biệt.

Kết luận

  • Khẳng định tính hiệu quả vượt trội của việc tích hợp công nghệ GIS (ArcGIS 9.2) và Viễn thám (ENVI 4.5) trong việc thành lập bản đồ trạng thái rừng từ ảnh SPOT-5 độ phân giải 2.5m.
  • Xây dựng thành công cơ sở dữ liệu không gian và thuộc tính đồng bộ cho 5 lớp chuyên đề bao gồm đất rừng, giao thông, thủy văn, địa hình và dân cư tại Vườn Quốc gia Ba Bể.
  • Kiểm soát chặt chẽ độ chính xác nắn ảnh qua 15 điểm GCP với sai số RMS dưới 1 pixel và đạt độ chính xác phân loại tổng thể trên 85% qua 100 điểm kiểm tra thực địa.
  • Phát triển thành công website bản đồ trạng thái rừng tương tác bằng ngôn ngữ JavaScript và PHP, mở ra hướng ứng dụng Web GIS linh hoạt trong quản lý lâm nghiệp hiện đại.
  • Cung cấp cơ sở khoa học thực tiễn vững chắc phục vụ công tác quy hoạch, bảo tồn nguồn gen quý hiếm của hơn 1.281 loài thực vật tại 8 xã thuộc huyện Ba Bể.

Về lộ trình tiếp theo, địa phương cần tiếp tục cập nhật chuỗi dữ liệu ảnh vệ tinh định kỳ giai đoạn 2012 - 2015 nhằm tự động hóa quy trình phát hiện biến động rừng. Các nhà quản lý và đơn vị bảo tồn hãy đẩy mạnh ứng dụng công nghệ GIS và viễn thám để tối ưu hóa hiệu quả quản trị tài nguyên thiên nhiên ngay hôm nay.