Tổng quan nghiên cứu

Khu vực Đồng bằng sông Cửu Long sở hữu hơn 128.537 ha rừng ngập mặn, chiếm trên 60% tổng diện tích rừng ngập mặn của cả nước, đóng vai trò huyết mạch trong việc bảo vệ bờ biển, chống xói lở và bảo tồn tài nguyên sinh học ven biển. Riêng tại tỉnh Cà Mau, diện tích rừng đạt gần 80.000 ha, tập trung lớn nhất tại huyện Ngọc Hiển với 36.613 ha và huyện Năm Căn. Tuy nhiên, hệ sinh thái ven biển này đang chịu sức ép nghiêm trọng từ các hoạt động chuyển đổi mục đích sử dụng đất sang nuôi trồng thủy sản tự phát, sạt lở bờ biển phía Đông với cường độ hơn 20 mét mỗi năm, xen kẽ quá trình bồi tụ 50 đến 80 mét mỗi năm tại bãi bồi Mũi Cà Mau.

Phương pháp điều tra thực địa truyền thống tốn kém nhiều nhân lực và chi phí, trong khi viễn thám quang học thường xuyên gặp trở ngại lớn do mây mù che phủ liên tục trong mùa mưa tại vùng nhiệt đới. Trước thách thức này, nghiên cứu đặt mục tiêu ứng dụng công nghệ viễn thám radar khẩu độ tổng hợp (SAR) để xác định, phân loại đối tượng rừng ngập mặn tại Mũi Cà Mau. Nghiên cứu thực hiện khảo sát định lượng đặc trưng tán xạ ngược của dữ liệu radar băng C (ENVISAT ASAR APP) và băng L (ALOS PALSAR) phân cực kép, đồng thời kết hợp ảnh radar với ảnh quang học SPOT 5 để đánh giá khả năng nâng cao độ chính xác phân loại. Phạm vi nghiên cứu được giới hạn tại khu vực Mũi Cà Mau với nguồn tư liệu đồng bộ giai đoạn 2010 - 2014, mang lại giải pháp công nghệ quan trắc độc lập với điều kiện thời tiết, phục vụ đắc lực công tác quản lý và phát triển rừng bền vững.

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

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng lý thuyết tương tác sóng điện từ siêu cao tần của hệ thống radar khẩu độ tổng hợp (SAR) và lý thuyết phản xạ phổ quang học đa phổ. Nguyên lý viễn thám radar dựa trên cơ chế phát xung vi sóng và thu nhận năng lượng tán xạ ngược từ bề mặt đất, được định lượng bằng hệ số tán xạ ngược (sigma nought) tính theo đơn vị decibel (dB). Mức độ tương tác của sóng radar phụ thuộc chặt chẽ vào bước sóng và cấu hình phân cực (HH, HV, VV):

  • Bước sóng băng L (~23,6 cm) có khả năng xuyên thấu sâu qua tầng tán lá cây rừng ngập mặn, phản xạ mạnh mẽ từ thân cây và tương tác tán xạ kép (double-bounce) với mặt nước hoặc bùn ngập bên dưới.
  • Bước sóng băng C (~5,6 cm) chủ yếu tương tác với các cấu trúc nhỏ ở bề mặt tán lá phía trên (tán xạ khối bề mặt).

Bên cạnh đó, nghiên cứu tích hợp chỉ số thực vật chuẩn hóa NDVI và phân tích thành phần chính PCA từ kênh phổ quang học SPOT 5. Kỹ thuật PCA giúp cô đọng hơn 80% lượng thông tin phổ vào thành phần chính đầu tiên (PC1), giảm thiểu hiện tượng tương quan giữa các kênh phổ và triệt tiêu nhiễu quang phổ trong phân loại thực phủ.

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

Nghiên cứu sử dụng hệ thống dữ liệu viễn thám đa nguồn, đa thời gian bao gồm:

  • 3 cảnh ảnh ALOS PALSAR phân cực kép HH/HV (kênh L, độ phân giải 12,7 m thu nhận vào tháng 8, 10 và 11 năm 2010).
  • 2 cảnh ảnh ENVISAT ASAR APP phân cực kép HH/VV và HH/HV (kênh C, độ phân giải 12,01 m và 12,19 m các năm 2010, 2012).
  • 1 cảnh ảnh quang học SPOT 5 đa phổ 4 kênh (độ phân giải 10 m thu nhận ngày 03/01/2011).

Về cỡ mẫu và phương pháp chọn mẫu, nghiên cứu thiết lập 246 điểm kiểm chứng độc lập phân bố theo mạng lưới hệ thống dạng mắt lưới 1 km x 1 km, tương ứng với 984 phần tử ảnh (kích thước mỗi điểm mẫu thực tế khoảng 30 m x 30 m, tương đương 4 pixel để loại trừ sai số định vị GPS cầm tay 10 m). Mẫu huấn luyện được thu thập thông qua đợt khảo sát thực địa từ ngày 30/11/2011 đến 06/12/2011 kết hợp bản đồ kiểm kê hiện trạng rừng năm 2010.

Quy trình xử lý ảnh bao gồm định chuẩn bức xạ tuyệt đối bằng bộ công cụ BEST và NEST, hiệu chỉnh hình học với sai số trung phương RMS dưới 0,75 pixel, khử nhiễu đốm bằng bộ lọc đa thời gian 9x9 kết hợp lọc không gian Gamma 3x3. Nghiên cứu triển khai 4 thuật toán phân loại có giám sát: Khoảng cách tối thiểu (MiD), Khoảng cách Mahalanobis (MaD), Gần đúng nhất (MLC) và Máy véc-tơ hỗ trợ (Support Vector Machine - SVM). Thuật toán SVM được lựa chọn làm trọng tâm nhờ khả năng tối ưu hóa siêu phẳng phân tách trên không gian đa chiều, hạn chế tối đa lỗi phân loại nhầm khi dữ liệu có độ phân tán lớn và phân bố không chuẩn. Hậu xử lý sử dụng phép lọc Majority Analysis cửa sổ 3x3 để loại bỏ các điểm ảnh dị biệt.

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

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

Quá trình phân tích định lượng dữ liệu viễn thám tại Mũi Cà Mau đã mang lại 4 phát hiện khoa học nổi bật:

  1. Khả năng phân cấp mật độ vượt trội của sóng radar băng L: Dữ liệu ALOS PALSAR (băng L) phân cực kép HH/HV đạt độ chính xác toàn cục trên 75% khi phân loại rừng ngập mặn theo 3 cấp độ che phủ (dày từ 70% trở lên, trung bình từ 50% đến 70%, thưa từ 30% đến 50%). Ngược lại, ảnh ENVISAT ASAR (băng C) đạt độ chính xác dưới 58%, gặp rất nhiều khó khăn trong việc phân biệt các trạng thái rừng do hiện tượng bão hòa tín hiệu tán xạ trên tầng tán.

  2. Nhận diện chính xác rừng mắm tiên phong trên bãi bồi: Dữ liệu radar băng L phân cực kép thể hiện khả năng nhận diện xuất sắc các quần xã mắm trắng (Avicennia alba) mới phát triển tại các bãi bồi thường xuyên ngập triều. Tín hiệu tán xạ kép đặc trưng giữa thân cây mắm non và nền bùn ướt giúp tách biệt rõ ràng đối tượng này với nền mặt nước biển và bãi triều trơ trụi.

  3. Hiệu quả nâng cao độ chính xác khi kết hợp ảnh quang học và radar: Đối với bài toán phân loại chi tiết 6 lớp đối tượng, việc tích hợp ALOS PALSAR với các kênh NDVI, PC1 và SWIR của ảnh SPOT 5 giúp tăng độ chính xác toàn cục lên trên 82% và hệ số Kappa tăng hơn 0,15 so với việc chỉ sử dụng ảnh radar đơn lẻ. Tuy nhiên, đối với bài toán 4 lớp đối tượng tổng quát, việc kết hợp chỉ cải thiện độ chính xác khoảng 2% đến 3% so với ảnh radar đã qua xử lý lọc nhiễu.

  4. Ưu thế vượt bậc của thuật toán SVM: Trong tất cả các kịch bản thử nghiệm phân loại 4 lớp và 6 lớp, thuật toán SVM luôn đạt kết quả cao nhất, vượt trội hơn thuật toán Maximum Likelihood (MLC) từ 4% đến 7% về độ chính xác toàn cục, đồng thời giảm thiểu tối đa hiện tượng nhầm lẫn giữa rừng thưa và mô hình mặt nước nuôi tôm xen kẽ rừng.

Thảo luận kết quả

Khả năng phân loại của radar phụ thuộc trực tiếp vào cơ chế tán xạ vi sóng theo cấu trúc hình học của thực vật. Sóng băng L với bước sóng dài 23,6 cm dễ dàng xuyên qua tầng lá để tương tác với các nhánh cây lớn và hệ thống rễ chống phức tạp của đước đôi (Rhizophora apiculata), tạo ra tín hiệu tán xạ khối mạnh mẽ ở kênh phân cực chéo HV. Ngược lại, bước sóng băng C chỉ dài 5,6 cm bị tán xạ tán lá chi phối, khiến độ chênh lệch hệ số tán xạ giữa rừng dày và rừng trung bình chỉ dao động từ 0,5 dB đến 1,2 dB, dẫn đến chồng lấn phổ nghiêm trọng trên biểu đồ tán xạ phân cực.

Khi biểu diễn trên ma trận sai số và không gian phân cực HH-HV, các đối tượng mặt nước sông rạch và đất dân cư thể hiện giá trị tách biệt rõ rệt. Tuy nhiên, việc phân biệt chi tiết tới cấp loài giữa đước và mắm trưởng thành trên ảnh radar đơn thời điểm còn gặp nhiều thách thức do độ gồ ghề bề mặt tương đồng. Sự bổ sung kênh hồng ngoại sóng ngắn SWIR và chỉ số NDVI từ ảnh SPOT 5 cung cấp thêm thông tin độ ẩm và lượng diệp lục, giúp bóc tách hoàn hảo ranh giới thực phủ phức tạp. Kết quả này hoàn toàn nhất quán với các công trình nghiên cứu quốc tế tại khu vực Nam Mỹ và Đông Nam Á, khẳng định tiềm năng to lớn của viễn thám radar trong giám sát sinh thái đất ngập nước.

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

Nhằm nâng cao hiệu quả quản lý, giám sát và bảo vệ tài nguyên rừng ngập mặn tại khu vực Đồng bằng sông Cửu Long, nghiên cứu đưa ra 4 khuyến nghị hành động cụ thể:

  1. Xây dựng quy trình tự động hóa tích hợp dữ liệu viễn thám radar và quang học: Chuẩn hóa quy trình tiền xử lý định chuẩn sigma chuẩn hóa và nắn chỉnh hình học tự động, rút ngắn thời gian thành lập bản đồ hiện trạng rừng xuống dưới 48 giờ sau khi nhận ảnh vệ tinh, bảo đảm độ chính xác bản đồ đạt trên 85%. Thời gian triển khai trong giai đoạn 2026 - 2028, chủ thể thực hiện là Cục Viễn thám Quốc gia phối hợp cùng Viện Địa lý Tài nguyên TP.HCM.

  2. Thiết lập hệ thống trạm quan trắc viễn thám định kỳ giám sát xói lở và bồi tụ: Ứng dụng ảnh radar đa thời gian băng L và băng C để theo dõi diễn biến 254 km đường bờ biển Cà Mau theo chu kỳ 3 tháng/lần, kịp thời cảnh báo sạt lở tại dải bờ biển phía Đông và khoanh vẽ diện tích rừng mắm mới bồi tụ tại Mũi Cà Mau. Thời gian triển khai từ quý 1 năm 2027 do Sở Tài nguyên và Môi trường tỉnh Cà Mau chủ trì.

  3. Nâng cấp mô hình phân loại ứng dụng máy học và học sâu: Nghiên cứu phát triển các mô hình phân loại tiên tiến dựa trên thuật toán SVM kết hợp mạng nơ-ron tích chập (CNN) trên dữ liệu radar phân cực toàn phần (Full-Polarimetry), đặt mục tiêu nâng độ chính xác phân biệt cấp loài cây lên trên 80% trên toàn bộ 80.000 ha rừng ngập mặn Cà Mau. Thời gian thực hiện giai đoạn 2027 - 2029 do các trường đại học và viện nghiên cứu chuyên ngành đảm trách.

  4. Ứng dụng bản đồ phân loại viễn thám vào kiểm soát quy hoạch lâm - ngư kết hợp: Sử dụng kết quả phân loại radar để kiểm tra, giám sát tỷ lệ che phủ rừng tối thiểu 60% trong các mô hình nuôi tôm sinh thái theo đúng quy định nhà nước tại các huyện Ngọc Hiển và Năm Căn. Thời gian thực hiện định kỳ hàng năm từ năm 2026 dưới sự giám sát của Chi cục Kiểm lâm và Ban Quản lý Vườn Quốc gia Mũi Cà Mau.

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

Nội dung và kết quả nghiên cứu của luận văn mang lại giá trị ứng dụng thiết thực cho 4 nhóm đối tượng chính:

  1. Cơ quan quản lý nhà nước về tài nguyên và kiểm lâm: Cán bộ Chi cục Kiểm lâm, Sở Tài nguyên và Môi trường có thể ứng dụng trực tiếp quy trình phân loại ảnh radar để cập nhật dữ liệu kiểm kê rừng thường niên mà không bị gián đoạn bởi thời tiết mùa mưa, phục vụ công tác quy hoạch bảo vệ rừng và phòng chống sạt lở ven biển.

  2. Kỹ sư và chuyên gia phân tích viễn thám - GIS: Các chuyên viên xử lý ảnh vệ tinh có thể khai thác tài liệu chi tiết về kỹ thuật định chuẩn bức xạ radar trên phần mềm NEST và BEST, phương pháp lọc nhiễu đa thời gian Spk9 kết hợp bộ lọc không gian Gamma 3x3, cùng quy trình cấu hình tham số tối ưu cho thuật toán SVM.

  3. Giảng viên, nghiên cứu sinh và học viên cao học: Các nhà nghiên cứu thuộc chuyên ngành Bản đồ, Viễn thám, Hệ thông tin Địa lý, Lâm nghiệp và Khoa học Môi trường có thể sử dụng luận văn làm tài liệu tham khảo chuẩn mực về phương pháp luận kết hợp ảnh quang học - radar và thiết kế mạng lưới kiểm chứng thực địa 246 điểm mẫu hệ thống.

  4. Các tổ chức bảo tồn và phát triển bền vững: Các ban quản lý khu dự trữ sinh quyển, tổ chức phi chính phủ trong và ngoài nước có thể ứng dụng các số liệu phân cấp mật độ rừng để xây dựng báo cáo đánh giá dịch vụ hệ sinh thái, ước tính trữ lượng sinh khối và lập dự án chi trả dịch vụ môi trường rừng ngập mặn.

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

  1. Tại sao ảnh radar băng L lại phân loại rừng ngập mặn hiệu quả hơn ảnh radar băng C? Ảnh radar băng L có bước sóng dài khoảng 23,6 cm, cho phép sóng vi ba xuyên thấu qua vòm tán lá rậm rạp để tương tác trực tiếp với các cành lớn, thân cây và hệ rễ chống ngập nước, tạo ra độ tương phản tán xạ rõ nét theo mật độ sinh khối. Trong khi đó, sóng băng C có bước sóng ngắn 5,6 cm bị tán xạ ngay tại bề mặt tán lá, nhanh chóng đạt ngưỡng bão hòa nên khó phân biệt cấu trúc rừng rậm rạp.

  2. Việc kết hợp dữ liệu ảnh radar và ảnh quang học mang lại lợi ích gì trong bài toán phân loại rừng? Sự kết hợp này phát huy tối đa ưu thế tương hỗ giữa hai nguồn dữ liệu. Ảnh radar ALOS PALSAR cung cấp thông tin về cấu trúc hình học và độ gồ ghề của thân cây, còn ảnh quang học SPOT 5 cung cấp thông tin phản xạ phổ diệp lục (NDVI) và độ ẩm (SWIR). Sự kết hợp này giúp nâng độ chính xác toàn cục cho bài toán phân loại chi tiết 6 lớp đối tượng lên trên 82%, tăng hơn 7% so với việc chỉ xử lý trên một nguồn ảnh riêng lẻ.

  3. Kỹ thuật lọc nhiễu đốm đa thời gian đóng vai trò như thế nào trong xử lý ảnh radar? Nhiễu đốm là hiện tượng giao thoa sóng ngẫu nhiên làm suy giảm độ tương phản của ảnh SAR. Bằng việc áp dụng bộ lọc đa thời gian Spk9 trên chuỗi 3 thời điểm ảnh ALOS PALSAR cùng quỹ đạo, giá trị độ lệch chuẩn của các khóa giải đoán giảm rõ rệt, giúp làm mịn vùng đồng nhất, giữ nguyên ranh giới đường bờ và cải thiện độ chính xác phân loại của thuật toán SVM thêm khoảng 5%.

  4. Dữ liệu radar có khả năng phân loại chi tiết đến từng loài cây rừng ngập mặn hay không? Dữ liệu radar phân cực kép đơn thời điểm chưa thể phân biệt hoàn toàn giữa các cá thể đước đôi và mắm trưởng thành do cấu trúc tán có độ nhám tương đương nhau. Tuy nhiên, ảnh radar băng L lại có khả năng nhận diện rất đặc hiệu đối với loài mắm trắng tiên phong mọc trên các bãi bồi ven biển nhờ cơ chế tán xạ kép nổi bật giữa thân cây non và mặt bùn lỏng thường xuyên ngập triều.

  5. Thuật toán SVM có ưu điểm gì vượt trội so với thuật toán Maximum Likelihood (MLC) truyền thống? Thuật toán phân loại xác suất MLC phụ thuộc chặt chẽ vào giả định dữ liệu mẫu phải tuân theo phân bố chuẩn đa chiều. Trong khi đó, thuật toán SVM hoạt động trên nguyên lý tối ưu hóa siêu phẳng phân tách trong không gian đặc trưng nhiều chiều, không đòi hỏi dữ liệu tuân theo phân bố chuẩn. Nhờ vậy, SVM kiểm soát nhiễu đốm tốt hơn và nâng độ chính xác phân loại toàn cục cao hơn MLC từ 4% đến 7%.

Kết luận

  1. Khẳng định tính hiệu quả và độ tin cậy cao của công nghệ viễn thám radar khẩu độ tổng hợp SAR băng L (ALOS PALSAR) trong việc xác định và phân cấp mật độ che phủ rừng ngập mặn độc lập với điều kiện thời tiết tại Mũi Cà Mau.
  2. Chứng minh giải pháp lọc nhiễu đa thời gian Spk9 kết hợp lọc không gian Gamma 3x3 giúp triệt tiêu nhiễu đốm hiệu quả, duy trì độ phân giải biên cạnh và bảo đảm tính đồng nhất cho các lớp khóa giải đoán.
  3. Thiết lập thành công mô hình kết hợp đa nguồn giữa ảnh radar băng L và ảnh quang học đa phổ SPOT 5, nâng độ chính xác phân loại 6 lớp đối tượng lên trên 82% với giải thuật học máy SVM.
  4. Cung cấp bộ dữ liệu 246 điểm kiểm chứng thực địa chuẩn hóa, tạo cơ sở khoa học vững chắc cho công tác quan trắc diễn biến bãi bồi và đường bờ biển tại tỉnh Cà Mau.
  5. Định hướng lộ trình nghiên cứu tiếp theo trong giai đoạn 2026 - 2030 hướng tới khai thác dữ liệu SAR phân cực toàn phần (Full-Polarimetry) kết hợp kỹ thuật giao thoa InSAR và mô hình học sâu (Deep Learning) để ước tính trữ lượng carbon sinh khối rừng ngập mặn.

Các cơ quan quản lý lâm nghiệp và địa phương cần đẩy mạnh ứng dụng công nghệ viễn thám radar tích hợp vào hệ thống giám sát tài nguyên định kỳ, góp phần hiện đại hóa công tác quản lý và bảo tồn bền vững hệ sinh thái rừng ngập mặn ven biển Việt Nam.