Nghiên Cứu Ứng Dụng Kết Hợp Ảnh Radar và Ảnh Quang Học Để Chiết Tách Lớp Phủ Mặt Đất

Tóm tắt nghiên cứu

Báo cáo nghiên cứu khoa học cấp Bộ "Nghiên cứu khả năng ứng dụng kết hợp ảnh vệ tinh radar và ảnh quang học để chiết tách một số đối tượng lớp phủ mặt đất" do ThS. Chu Hải Tùng làm chủ nhiệm cùng tập thể cán bộ Trung tâm Viễn thám Quốc gia (Bộ Tài nguyên và Môi trường) thực hiện, giải quyết bài toán cốt lõi: Làm thế nào để tích hợp tối ưu ưu thế xuyên mây, nhạy cảm cấu trúc của sóng vi ba (SAR) với đặc trưng phản xạ phổ trực quan của ảnh quang học nhằm nâng cao độ chính xác trong phân loại lớp phủ bề mặt tại Việt Nam?

Về phương pháp, đề tài xây dựng quy trình công nghệ đa tầng bao gồm: nắn chỉnh hình học bằng mô hình vật lý kết hợp mô hình số độ cao (DEM); triệt tiêu nhiễu đốm (speckle noise) bằng phin lọc thích nghi; và áp dụng các kỹ thuật tổ hợp không gian màu (IHS, PCA, NDVI-SAR). Nghiên cứu được triển khai thực nghiệm trên nguồn dữ liệu đa thời gian (SPOT 2/3/5, Landsat TM/ETM+, ERS-1/2, ENVISAT/ASAR) tại hai khu vực trọng điểm phía Bắc (Hải Phòng - Quảng Ninh - Hải Dương) và phía Nam (TP. Hồ Chí Minh - Đồng Nai - Bà Rịa Vũng Tàu - Long An - Tiền Giang).

Kết quả chỉ ra rằng giải pháp hợp nhất dữ liệu radar – quang học giúp phân biệt rõ nét các đối tượng có cùng phổ phản xạ quang học nhưng khác nhau về độ gồ ghề và độ ẩm như rừng ngập mặn, đầm nuôi thủy sản, vùng chuyên canh lúa và khu dân cư nông thôn, nâng cao độ chính xác phân loại tổng thể thêm 5–10% so với sử dụng ảnh đơn lẻ. Đây là bước đệm công nghệ quan trọng chuẩn bị cho việc khai thác Trạm thu ảnh vệ tinh quốc gia.

Bối cảnh và tầm quan trọng

Nhu cầu giám sát lớp phủ mặt đất trong phát triển bền vững

Lớp phủ mặt đất (land cover) là trạng thái vật chất sinh thái và nhân tạo của bề mặt Trái Đất, bao gồm thảm thực vật, thổ nhưỡng, công trình xây dựng và mặt nước. Tại Việt Nam, quá trình công nghiệp hóa, đô thị hóa nhanh chóng cùng các áp lực phát triển kinh tế đã kéo theo sự suy giảm tài nguyên rừng, thu hẹp đất ngập nước tự nhiên, chuyển đổi ồ ạt đất nông nghiệp sang nuôi trồng thủy sản và nguy cơ ô nhiễm môi trường trên diện rộng. Việc theo dõi kịp thời, chính xác diễn biến lớp phủ mặt đất trở thành yêu cầu sống còn cho công tác quy hoạch không gian lãnh thổ và ứng phó biến đổi khí hậu.

Khoảng trống nghiên cứu và bài toán thời tiết nhiệt đới

Phương pháp viễn thám quang học truyền thống (Landsat, SPOT, ASTER) là công cụ quen thuộc trong thành lập bản đồ lớp phủ nhờ dải phổ nhìn thấy và cận hồng ngoại trực quan. Tuy nhiên, Việt Nam nằm trong vành đai khí hậu nhiệt đới gió mùa ẩm với lượng mây mù bao phủ quanh năm rất lớn, tạo nên các "khoảng trống dữ liệu" nghiêm trọng vào mùa mưa – thời điểm diễn ra nhiều biến động sinh thái và thiên tai lũ lụt nhất. Bên cạnh đó, ảnh quang học chỉ phản ánh tính chất hấp thụ bề mặt trong thị tần mà thiếu hoàn toàn thông tin về cấu trúc hình học, độ nhám và độ ẩm của vật thể.

Ngược lại, hệ thống cảm biến vi ba chủ động Synthetic Aperture Radar (SAR) hoạt động ở bước sóng dài (băng C, L, X), có khả năng xuyên qua mây mù, hoạt động bất kể ngày đêm và rất nhạy cảm với độ ẩm, độ gồ ghề bề mặt. Dẫu vậy, ảnh radar tồn tại các rào cản kỹ thuật lớn: biến dạng hình học phức tạp (co ngắn, chồng đè, bóng địa hình) và mức độ nhiễu đốm (speckle noise) cao. Các nghiên cứu trước đây ở Việt Nam chỉ dừng lại ở thử nghiệm đơn lẻ một loại tư liệu hoặc một phương pháp trộn ảnh cục bộ. Đề tài cấp Bộ này là công trình đầu tiên thiết lập một phương pháp luận đồng bộ, nghiên cứu toàn diện từ bản chất tán xạ, giải pháp tiền xử lý đến quy trình tích hợp quang học – radar chuyên biệt cho điều kiện lãnh thổ Việt Nam.

Phương pháp luận và Tiếp cận kỹ thuật

Thiết kế nghiên cứu và Nguồn tư liệu

Đề tài xây dựng một quy trình công nghệ khép kín từ dữ liệu gốc đến sản phẩm bản đồ chuyên đề tỷ lệ 1:100.000. Dữ liệu đầu vào kết hợp giữa các cảm biến quang học đa phổ và cảm biến radar không gian tại hai mốc thời gian (1996 và 2004):

  • Tư liệu ảnh quang học: SPOT 2, 3 XS (độ phân giải 20m), SPOT 5 XS (độ phân giải 10m), Landsat 5 TM / Landsat 7 ETM+ (độ phân giải 30m đa phổ, 15m toàn sắc).
  • Tư liệu ảnh Radar: ERS-1, ERS-2 SAR (băng C: 5.67 cm, phân cực VV, góc tới $23^\circ$), ENVISAT/ASAR (băng C, chế độ chuẩn IM, độ phân giải 30m), tham chiếu thêm RADARSAT-1 và ALOS-PALSAR (băng L: 23.6 cm).
   Ảnh Radar (ERS/ASAR)                   Ảnh Quang học (SPOT/Landsat)

Xử lý biến dạng hình học bằng Mô hình vật lý và DEM

Do đặc thù chụp quét nghiêng từ một phía, ảnh radar chịu ảnh hưởng nghiêm trọng của địa hình, tạo ra sai số dịch chuyển điểm $D$ theo công thức: $$D = \Delta h \cdot \tan(90^\circ - \eta)$$ (Trong đó $\Delta h$ là chênh cao địa hình, $\eta$ là góc nghiêng tia radar).

Với ảnh ERS có góc chụp trung tâm $23^\circ$, chênh cao $\Delta h = 100\text{ m}$ sẽ làm lệch vị trí đối tượng trên ảnh tới $235.6\text{ m}$. Đề tài đã áp dụng thuật toán nắn chỉnh trực giao bằng mô hình vật lý kết hợp mô hình số độ cao (DEM) nhằm loại trừ triệt để hiện tượng co ngắn phía trước (foreshortening), chồng đè (layover) và dịch chuyển tọa độ, đảm bảo chồng khít hình học tuyệt đối giữa ảnh radar và ảnh quang học.

Xử lý nhiễu đốm (Speckle Noise)

Nhiễu đốm phát sinh do sự giao thoa không cùng pha của các sóng phản xạ vi ba phân tán. Nghiên cứu áp dụng các phin lọc thích nghi theo vùng như Lee Filter và Sigma Filter. Kết quả phân tích tương quan giữa độ lệch chuẩn và giá trị trung bình cục bộ ($Z = X \cdot V$) cho thấy phin lọc Lee giữ lại được cấu trúc đường biên sắc nét của sông suối, ranh giới thửa đất mà vẫn làm mịn hiệu quả các vùng đồng nhất.

Chiến lược tổ hợp dữ liệu và Chiết tách thông tin

Đề tài thử nghiệm và chuẩn hóa các mô hình tích hợp không gian phổ:

  1. Biến đổi IHS (Intensity - Hue - Saturation): Thay thế kênh cường độ sáng $I$ bằng dữ liệu tán xạ ngược radar, giữ lại màu sắc $H$ và độ bão hòa $S$ từ ảnh quang học.
  2. Phân tích thành phần chính (PCA): Trộn các kênh phổ quang học với ảnh SAR đa thời gian để cô lập biến động.
  3. Tổ hợp màu giả đa nguồn: Xây dựng tổ hợp màu chuyên dụng (Red: Chỉ số thực vật NDVI; Green: Trung bình tán xạ radar $\text{TB_SAR}$; Blue: Kênh quang học SPOT 3/Landsat).

Phân loại đối tượng được tiến hành bằng thuật toán phân loại có kiểm định (Supervised Classification - Maximum Likelihood) kết hợp với giải đoán trực quan trên máy vi tính và đối soát thực địa bằng GPS.

Phát hiện chính của nghiên cứu

1. Phân biệt chính xác giữa Rừng ngập mặn và Vùng nuôi trồng thủy sản

Tại khu vực ven biển Hải Phòng - Quảng Ninh và Cần Giờ (TP.HCM), trên ảnh quang học đơn lẻ, các đầm nuôi thủy sản có thực vật trôi nổi hoặc bờ bao rậm rạp thường có giá trị phổ phản xạ rất giống với rừng ngập mặn thưa, gây nhầm lẫn lớn. Khi kết hợp với ảnh radar:

  • Rừng ngập mặn có cấu trúc tán lá phức tạp, thân cây ngập nước tạo nên hiệu ứng tán xạ góc kép (double-bounce scattering) và tán xạ khối cực mạnh $\rightarrow$ cho tông màu rất sáng trên ảnh SAR.
  • Mặt nước đầm tôm tạo ra hiện tượng phản xạ gương (specular reflection), sóng vi ba phản xạ ra xa khỏi ăng-ten thu $\rightarrow$ cho tông màu tối hoàn toàn. Sự tương phản này giúp phân lập ranh giới đầm tôm và rừng ngập mặn với độ chính xác đạt trên 90%.

2. Tách biệt Đất chuyên lúa nước và Đất hoa màu ngắn ngày

Ở vùng Đồng bằng sông Hồng và đồng bằng ven biển Nam Bộ:

  • Đất chuyên lúa có chu kỳ ngập nước nhân tạo tạo ra sự biến thiên tán xạ radar rất mạnh theo mùa vụ (tối sẫm khi cấy nước, tăng dần độ sáng khi lúa làm đòng và phát triển sinh khối tán).
  • Đất trồng rau màu, cây họ đậu có độ gồ ghề bề mặt ổn định và không có nền nước $\rightarrow$ tán xạ ngược ổn định. Sử dụng tổ hợp ảnh radar đa thời gian kết hợp kênh phổ quang học đã loại bỏ hoàn toàn sự nhầm lẫn giữa lúa nước và cây ngắn ngày.

3. Tăng cường nhận dạng Khu dân cư nông thôn

Khu dân cư nông thôn Việt Nam có đặc trưng nhà ở xen lẫn vườn cây ăn trái rậm rạp, khiến ảnh quang học rất dễ phân loại nhầm thành đất rừng hoặc vườn cây lâu năm. Tuy nhiên, các công trình xây dựng, tường nhà và mái ngói kim loại phản xạ theo cơ chế tán xạ góc vuông kim loại/bê tông cực mạnh trên ảnh radar (cho các điểm sáng chói corner reflector). Khi trộn màu IHS, các điểm dân cư nông thôn hiện rõ cấu trúc vân ảnh (texture) đặc thù, phân biệt rõ rệt với cây cối bao quanh.

4. Nâng cao độ chính xác phân loại tổng thể

Đánh giá định lượng thông qua chỉ số phân tách lớp (Transformed Divergence) cho thấy hầu hết các cặp lớp phủ đều đạt giá trị $> 1.9$ (mức phân tách rất tốt). Độ chính xác phân loại tổng quát trên các tập dữ liệu kết hợp đạt từ 74.6% đến 82.0%, vượt trội đáng kể so với mức 69.3% khi chỉ khai thác kênh phổ quang học đơn lẻ.

Đóng góp khoa học và Giá trị thực tiễn

Đóng góp về mặt lý thuyết và học thuật

  • Hệ thống hóa toàn diện cơ sở vật lý vi sóng ứng dụng cho viễn thám tài nguyên tại Việt Nam, phân tích chi tiết cơ chế tương tác vi ba (tán xạ bề mặt, tán xạ khối, tán xạ góc) đối với các hệ sinh thái nhiệt đới đặc thù.
  • Xây dựng hệ thống phân loại lớp phủ mặt đất 2 cấp gồm 17 lớp đối tượng, tương thích với tiêu chuẩn quốc tế (FAO LCC, CORINE) và phù hợp với khả năng phân giải của dữ liệu vệ tinh đa nguồn.

Đổi mới về mặt phương pháp luận

  • Hoàn thiện quy trình công nghệ kết hợp ảnh radar và quang học đạt chuẩn quy phạm, tích hợp thành công mô hình số độ cao (DEM) vào xử lý hình học SAR và tối ưu hóa các tham số lọc nhiễu thích nghi Lee/Sigma trước khi phân loại.
  • Chuẩn hóa bộ công thức tổ hợp dải phổ lai ghép (RGB kết hợp NDVI, SAR trung bình và quang học) phục vụ giải đoán bán tự động trên phần mềm máy tính.

Ý nghĩa thực tiễn và chính sách

  • Trực tiếp ứng dụng thành lập thành công 4 bộ bản đồ hiện trạng lớp phủ tỷ lệ 1:100.000 tại hai vùng kinh tế trọng điểm phía Bắc (Hải Phòng) và phía Nam (TP. Hồ Chí Minh) qua 2 thời kỳ (1996 và 2004), phục vụ giám sát biến động đô thị hóa và diện tích rừng.
  • Tạo tiền đề kỹ thuật, nguồn nhân lực và chuẩn bị công nghệ tiếp nhận, khai thác hiệu quả Trạm thu ảnh vệ tinh quốc gia (dự án ODA Pháp thuộc Bộ Tài nguyên và Môi trường), khai thác tối đa nguồn ảnh quang học (SPOT) và radar (ASAR).

Đối tượng hưởng lợi và Phạm vi ứng dụng

  • Cơ quan quản lý Tài nguyên & Môi trường: Bộ TN&MT, Sở TN&MT các tỉnh/thành có công cụ giám sát biến động lớp phủ đất định kỳ, kiểm soát chặt chẽ quá trình đô thị hóa và chuyển đổi mục đích sử dụng đất trái phép.
  • Ngành Nông nghiệp & Phát triển Nông thôn: Giám sát cơ cấu mùa vụ lúa nước, đánh giá biến động suy thoái rừng ngập mặn và quy hoạch hợp lý diện tích nuôi trồng thủy sản ven biển.
  • Cộng đồng Khoa học Địa tin học (GIS & RS): Cung cấp tài liệu tham khảo chuẩn mực về thuật toán xử lý dữ liệu vi ba và phương pháp tích hợp ảnh đa cảm biến.
  • Ban chỉ đạo Phòng chống Thiên tai: Ứng dụng khả năng quan sát xuyên mây để khoanh vùng ngập úng thời gian thực trong các đợt mưa bão, triều cường lớn.

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

1. Phát hiện quan trọng nhất của đề tài là gì?

Điểm đột phá nhất là việc chứng minh và định lượng hóa được khả năng của tổ hợp dữ liệu radar – quang học trong việc giải quyết triệt để sự nhầm lẫn giữa các đối tượng có cùng phản xạ phổ quang học nhưng khác nhau về hình thái cấu trúc vật lý: đặc biệt là cặp đối tượng Rừng ngập mặn – Vùng nuôi trồng thủy sản và Đất lúa – Đất rau màu/Làng quê.

2. Phương pháp nắn chỉnh hình học ảnh radar trong đề tài có điểm gì đặc biệt?

Thay vì dùng hàm đa thức thông thường (chỉ hiệu quả ở đồng bằng), đề tài sử dụng mô hình nắn ảnh vật lý kết hợp mô hình số độ cao (DEM). Phương pháp này mô hình hóa chính xác quỹ đạo vệ tinh và triệt tiêu sai số dịch chuyển do góc chụp nghiêng gây ra ở vùng có địa hình đồi núi, giúp ảnh radar chồng khít tuyệt đối lên lưới bản đồ tọa độ quốc gia.

3. Kết quả nghiên cứu có thể khái quát hóa (generalize) cho toàn quốc không?

Hoàn toàn có thể. Hai khu vực thử nghiệm (Hải Phòng - Quảng Ninh ở miền Bắc và TP.HCM - Đồng Nai - Tiền Giang ở miền Nam) đại diện cho đầy đủ các dạng địa hình phức tạp, hệ sinh thái nhạy cảm (rừng tự nhiên, rừng ngập mặn, bãi bồi ven biển) và các hoạt động kinh tế - nhân sinh đặc trưng của Việt Nam. Quy trình công nghệ được thiết kế chuẩn hóa để sẵn sàng áp dụng cho bất kỳ địa phương nào.

4. Bước phát triển tiếp theo của hướng nghiên cứu này là gì?

Hướng nghiên cứu tiếp theo sẽ tập trung vào khai thác dữ liệu Radar giao thoa (InSAR/DInSAR) để quan trắc lún sụt đô thị, sử dụng ảnh SAR đa phân cực toàn phần (Polarimetric SAR - PolSAR) và tích hợp trí tuệ nhân tạo (Deep Learning/CNNs) để tự động hóa hoàn toàn quá trình giải đoán lớp phủ.

5. Lợi ích kinh tế trực tiếp của giải pháp kết hợp này là gì?

Giúp tiết kiệm hàng tỷ đồng chi phí điều tra thực địa ngoại nghiệp, rút ngắn thời gian thành lập bản đồ hiện trạng từ nhiều tháng xuống vài tuần, và quan trọng nhất là loại bỏ rủi ro thiếu dữ liệu do thời tiết xấu, đảm bảo tính cập nhật liên tục cho các quyết định chính sách cấp bách.

Kết luận

Báo cáo nghiên cứu khoa học của ThS. Chu Hải Tùng và nhóm tác giả Trung tâm Viễn thám Quốc gia đã khẳng định tính tất yếu và hiệu quả vượt trội của phương pháp tích hợp ảnh vệ tinh radar và quang học trong giám sát tài nguyên môi trường. Công trình đã chuẩn hóa thành công quy trình công nghệ từ tiền xử lý hình học phức tạp, giải trừ nhiễu vi ba đến các thuật toán trộn màu và chiết tách chuyên đề tỷ lệ 1:100.000. Trong kỷ nguyên phát triển mạnh mẽ của công nghệ quan sát Trái Đất, việc làm chủ quy trình tích hợp cảm biến đa nguồn là chìa khóa mở đường cho việc xây dựng cơ sở dữ liệu không gian số chính xác, phục vụ đắc lực cho mục tiêu phát triển kinh tế - xã hội bền vững và chủ động ứng phó biến đổi khí hậu tại Việt Nam.