Tổng quan nghiên cứu

Việc giám sát biến động lớp phủ mặt đất đóng vai trò then chốt trong quản lý tài nguyên thiên nhiên, quy hoạch đô thị và ứng phó biến đổi khí hậu. Tuy nhiên, tại các khu vực nhiệt đới gió mùa như Việt Nam, độ che phủ mây quanh năm lên tới hơn 50% ở phần lớn các thời điểm thu nhận ảnh, gây cản trở nghiêm trọng đến việc trích xuất thông tin từ ảnh viễn thám quang học. Đối với thành phố Hà Nội với tổng diện tích khoảng 3.300 km2, trong đó đất nông nghiệp chiếm tới 56,6% và diện tích trồng lúa nước chiếm 60,9% khu vực nông nghiệp, đặc tính biến động theo mùa vụ của cây trồng đòi hỏi phải có chuỗi quan sát liên tục để nhận dạng chính xác.

Luận văn thạc sĩ chuyên ngành Khoa học máy tính của tác giả Mẫn Đức Chúc, dưới sự hướng dẫn của Tiến sĩ Nguyễn Thị Nhật Thanh tại Trường Đại học Công nghệ – Đại học Quốc gia Hà Nội, đã tập trung giải quyết bài toán phân loại lớp phủ mặt đất tại khu vực có độ che phủ mây cao và biến động mùa vụ phức tạp. Đề tài thuộc khuôn khổ Chương trình Khoa học Công nghệ Vũ trụ Việt Nam với mã số tài trợ VT-UD/06/16-20. Mục tiêu nghiên cứu cụ thể là xây dựng quy trình ghép ảnh không mây tối ưu từ dữ liệu vệ tinh Landsat 8 Surface Reflectance giai đoạn 2013–2016 và ứng dụng mô hình học máy kết hợp (ensemble learning) để nâng cao độ chính xác phân loại cho năm 2016. Kết quả nghiên cứu đã được công bố trên tạp chí quốc tế uy tín International Journal of Remote Sensing, mang lại giải pháp công nghệ quan trọng giúp nâng độ chính xác toàn diện đạt 83,91% và hệ số Kappa đạt 0,79, đồng thời giảm mức sai số ước tính diện tích lúa thực tế xuống chỉ còn 4,79% so với số liệu thống kê chính thức.

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 ứng dụng nền tảng lý thuyết phản xạ phổ viễn thám quang học và các thuật toán học máy hiện đại trong bài toán phân loại đa lớp. Cơ sở lý thuyết tập trung vào các khái niệm chính gồm: phản xạ bề mặt (Surface Reflectance), chỉ số thực vật khác biệt chuẩn hóa (NDVI), chỉ số đất trống (BSI), và phương pháp tổng hợp điểm ảnh tốt nhất (Best-Available-Pixel - BAP). Sự biến thiên của các chỉ số NDVI và BSI theo thời gian phản ánh rõ chu kỳ sinh trưởng của thực vật và đặc điểm thổ nhưỡng, hỗ trợ phân tách 7 lớp phủ chính: lúa nước, cây hoa màu, cỏ và cây bụi, cây thân gỗ, đất trống, bề mặt không thấm và mặt nước.

Mô hình học máy áp dụng bao gồm 5 thuật toán cơ sở: Hồi quy Logistic (Logistic Regression), Máy vector hỗ trợ với hàm nhân RBF và nhân tuyến tính (SVM-RBF, SVM-Linear), Mạng nơ-ron truyền thẳng đa tầng (Multilayer Perceptron - MLP) và Thuật toán tăng cường độ dốc cực đại (eXtreme Gradient Boosting - XGBoost). Bên cạnh đó, lý thuyết học kết hợp (Ensemble Learning) theo cơ chế bỏ phiếu đa số (Majority Voting) được vận dụng để tích hợp thế mạnh phân loại của từng mô hình đơn lẻ, giúp tối ưu hóa ranh giới phân tách trên không gian đặc trưng phức tạp gồm 35 kênh phổ - thời gian kết hợp với 7 giá trị độ lệch chuẩn trung bình (MSD).

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

Nguồn dữ liệu đầu vào bao gồm 54 cảnh ảnh Landsat 8 Surface Reflectance (L8SR) độ phân giải không gian 30 m thu thập từ cổng thông tin USGS Earth Explorer trong giai đoạn 4 năm từ 2013 đến 2016. Dữ liệu sau đó được ghép thành 27 ảnh hoàn chỉnh bao phủ toàn bộ lãnh thổ Hà Nội. Thuật toán ghép ảnh BAP cải tiến chấm điểm từng điểm ảnh tiềm năng dựa trên 4 tiêu chí: khoảng cách thời gian năm (Year score mở rộng 3 năm hỗ trợ 2013-2015), ngày trong năm (DOY score theo phân phối Gauss), độ mờ khí quyển (Opacity score theo hàm Sigmoid dựa trên 4 mức aerosol) và khoảng cách tới mây hoặc bóng mây (bán kính 50 pixel). Quy trình này tạo ra chuỗi 5 ảnh tổng hợp gần như sạch mây hoàn toàn (độ phủ sạch mây trên 99%) tương ứng với các ngày thứ 137, 153, 169, 265 và 281 của năm 2016.

Cỡ mẫu nghiên cứu được thiết kế theo phương pháp lấy mẫu phân tầng ngẫu nhiên (Stratified Random Sampling) dựa trên bản đồ hiện trạng sử dụng đất 11 tầng của Sở Tài nguyên và Môi trường Hà Nội. Tổng số mẫu thu thập gồm 5.079 điểm huấn luyện và 2.748 điểm kiểm thử độc lập, được kiểm chứng và gán nhãn chính xác thông qua ảnh vệ tinh độ phân giải cao trên Google Earth kết hợp khảo sát thực địa. Kỹ thuật kiểm định chéo 10 nếp (10-fold cross validation) cùng tìm kiếm lưới (Grid Search) được áp dụng trên tập huấn luyện để tinh chỉnh siêu tham số tối ưu trước khi đánh giá mô hình trên tập kiểm thử độc lập.

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

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

Thứ nhất, phương pháp ghép ảnh BAP cải tiến đã loại bỏ hiệu quả tác động của mây trên ảnh vệ tinh quang học. Trước khi xử lý, tỷ lệ che phủ mây trung bình trên 5 ảnh mục tiêu là 20,54%, trong đó ngày 169 bị mây che phủ nặng nhất lên tới 73,81%. Sau khi áp dụng thuật toán BAP với các ảnh hỗ trợ (chủ yếu từ năm 2015 đóng góp 72,11% số điểm ảnh thay thế), toàn bộ 5 ảnh composite đều đạt độ phủ sạch mây trên 99%, tạo ra chuỗi dữ liệu đầu vào liền mạch.

Thứ hai, việc sử dụng chuỗi ảnh thời gian (Time-series Composites) mang lại hiệu quả vượt trội hoàn toàn so với phân loại trên từng ảnh đơn lẻ. Khi phân loại bằng thuật toán XGBoost, mô hình chuỗi thời gian đạt độ chính xác toàn diện cao hơn trung bình 10,03% và chỉ số Kappa cao hơn 0,13 so với các mô hình ảnh đơn lẻ. Các ảnh đơn lẻ có kết quả rất bấp bênh, độ chính xác dao động từ 68,68% đến 73,81% do không bao quát được chu kỳ sinh trưởng của các lớp phủ động.

Thứ ba, mô hình tổ hợp 5 thuật toán học máy (Ensemble) đạt hiệu năng tổng thể cao nhất với độ chính xác toàn diện OA đạt 83,96% và hệ số Kappa đạt 0,79, vượt qua tất cả các mô hình đơn lẻ như SVM-Linear (81,94%), MLP (82,88%), SVM-RBF (83,41%) và XGBoost (83,84%). Đặc biệt, mô hình Ensemble khắc phục triệt để điểm yếu của XGBoost và MLP tại các lớp phủ khó như đất trống (F1-score tăng từ 0,23 lên 0,46) và cỏ/cây bụi (F1-score đạt 0,49).

Thứ tư, khi đối chiếu với số liệu thống kê của Cục Thống kê Hà Nội về diện tích gieo cấy lúa vụ xuân hè (khoảng 99.454 ha), bản đồ phân loại từ mô hình Ensemble cho kết quả ước tính 104.218 ha, sai số tuyệt đối chỉ là 4.764 ha (tương đương sai số 4,79%). Trong khi đó, các kịch bản ảnh đơn lẻ có mức chênh lệch rất lớn, dao động từ 7.811 ha đến 16.785 ha với sai số từ 7,85% đến 16,88%.

Thảo luận kết quả

Sự vượt trội của mô hình chuỗi thời gian xuất phát từ khả năng nắm bắt trọn vẹn đặc trưng biến đổi phổ theo mùa của thực vật. Tại các thời điểm đơn lẻ như ngày 137, 153 và 265, lúa non trong giai đoạn xanh thường bị nhầm lẫn nghiêm trọng với hoa màu và cỏ cây bụi. Đến ngày 169, các cánh đồng lúa bị ngập nước chuẩn bị cấy lại phản xạ phổ tương tự mặt nước. Ở ngày 281, sau khi thu hoạch, gốc rạ và đất khô bị nhận diện nhầm thành đất trống hoặc bề mặt không thấm. Khi xếp chồng 35 kênh phổ của 5 thời điểm cùng 7 kênh độ lệch chuẩn trung bình, mô hình đã tách biệt hoàn hảo lớp lúa nước với điểm F1 đạt 0,90.

Mô hình kết hợp Ensemble đã phát huy cơ chế bù trừ khuyết điểm giữa các thuật toán. Trong khi XGBoost xử lý xuất sắc các lớp chiếm diện tích lớn như mặt nước và lúa nước nhưng lại suy giảm hiệu năng ở các lớp có ít mẫu như đất trống (chỉ có 75 mẫu huấn luyện), thì SVM-RBF với hàm phân tách phi tuyến trong không gian nhiều chiều đã bù đắp xuất sắc sự thiếu hụt này. Sự nhất quán giữa kết quả viễn thám và số liệu thống kê kinh tế - xã hội khẳng định độ tin cậy vượt bậc của phương pháp luận đề xuất.

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

Dựa trên kết quả thực nghiệm của luận văn, 4 nhóm giải pháp và khuyến nghị kỹ thuật được đề xuất nhằm nâng cao hiệu quả ứng dụng viễn thám trong quản lý lãnh thổ:

  1. Triển khai chuẩn hóa quy trình tiền xử lý và ghép ảnh tự động theo phương pháp BAP mở rộng cho toàn bộ hệ thống cơ sở dữ liệu ảnh vệ tinh quang học tại các cơ quan viễn thám quốc gia. Đích ngắm là đạt tỷ lệ làm sạch mây trên 95% trên các chuỗi ảnh Landsat và Sentinel-2 phục vụ quan trắc biến động lớp phủ mặt đất thường niên, hoàn thành tích hợp hệ thống trong vòng 12 tháng do Trung tâm Vũ trụ Việt Nam và Cục Viễn thám Quốc gia chủ trì.
  2. Ứng dụng mô hình học máy kết hợp đa thuật toán (Ensemble Learning) vào quy trình kiểm kê hiện trạng sử dụng đất định kỳ 5 năm một lần tại Sở Tài nguyên và Môi trường các tỉnh thành phố phía Bắc. Giải pháp này giúp duy trì độ chính xác toàn diện trên 85%, giảm thiểu 50% chi phí và nhân lực khảo sát thực địa thủ công giai đoạn 2026–2030.
  3. Tích hợp dữ liệu lịch mùa vụ nông nghiệp địa phương vào thuật toán lựa chọn ảnh thành phần nhằm khắc phục hiện tượng loang lổ phổ do chênh lệch pha sinh trưởng của cây trồng khi ghép ảnh đa năm. Chi cục Trồng trọt và Bảo vệ Thực vật phối hợp các đơn vị công nghệ thông tin xây dựng lịch số hóa chi tiết theo từng tiểu vùng sinh thái trong vòng 6 tháng.
  4. Nghiên cứu mở rộng kết hợp dữ liệu ảnh quang học với dữ liệu ảnh vi sóng (Radar SAR Sentinel-1) không phụ thuộc thời tiết để bổ sung thông tin cấu trúc bề mặt, hướng tới mục tiêu giảm sai số phân loại đất trống và bề mặt nhân tạo xuống dưới 5% trong vòng 2 năm tới.

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

Nội dung và phương pháp trong luận văn đem lại giá trị học thuật và ứng dụng thực tiễn cho nhiều nhóm độc giả:

  • Học viên cao học và nghiên cứu sinh ngành Khoa học máy tính, Hệ thống thông tin địa lý (GIS) và Viễn thám: Cung cấp tài liệu tham khảo chuẩn mực về quy trình xử lý dữ liệu ảnh vệ tinh đa thời gian, trích xuất đặc trưng không gian - phổ và kỹ thuật lập trình học máy chuyên sâu với thư viện Scikit-learn, Numpy, GDAL.
  • Chuyên viên phân tích dữ liệu không gian tại các viện nghiên cứu và trung tâm viễn thám: Nắm vững thuật toán ghép ảnh BAP cải tiến để xử lý triệt để bài toán mây che phủ dày đặc tại các vùng nhiệt đới.
  • Cán bộ quản lý quy hoạch tại Sở Tài nguyên và Môi trường, Sở Nông nghiệp và Phát triển Nông thôn: Ứng dụng bản đồ hiện trạng lớp phủ và phương pháp ước tính diện tích mùa vụ có độ chính xác cao vào công tác kiểm kê đất đai và giám sát an ninh lương thực.
  • Kỹ sư phát triển phần mềm và giải pháp trí tuệ nhân tạo (AI): Tận dụng mô hình Ensemble kết hợp 5 thuật toán cơ sở để xây dựng các sản phẩm tự động hóa phân tích ảnh vệ tinh thương mại phục vụ nông nghiệp thông minh và quản lý đô thị.

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

Tại sao việc phân loại lớp phủ bằng ảnh viễn thám quang học tại miền Bắc Việt Nam lại gặp nhiều khó khăn?
Khu vực miền Bắc và đặc biệt là Hà Nội có đặc trưng khí hậu nhiệt đới gió mùa với tần suất che phủ mây quanh năm rất lớn, tới 15 trên 27 ảnh vệ tinh có độ che mây trên 50%. Tình trạng thiếu ảnh sạch mây liên tục khiến việc quan sát các giai đoạn phát triển sinh học của cây trồng bị gián đoạn, gây nhầm lẫn phổ nghiêm trọng giữa lúa, hoa màu và đất trống.

Phương pháp ghép ảnh Best-Available-Pixel (BAP) cải tiến trong nghiên cứu có ưu điểm gì nổi bật?
Thuật toán BAP cải tiến tính toán điểm số tổng hợp dựa trên 4 chỉ tiêu chặt chẽ: khoảng cách năm hỗ trợ (mở rộng 3 năm), độ lệch ngày trong năm, độ mờ khí quyển aerosol và khoảng cách an toàn tới mây/bóng mây bán kính 50 pixel. Nhờ đó, phương pháp tạo ra chuỗi 5 ảnh composite năm 2016 đạt độ sạch mây trên 99%, đảm bảo tính đồng nhất về phổ và hình học cho khu vực nghiên cứu.

Tại sao phân loại trên chuỗi ảnh thời gian (Time-series) lại chính xác hơn nhiều so với ảnh đơn lẻ?
Chuỗi ảnh thời gian ghi lại đầy đủ các pha biến đổi phản xạ phổ theo chu kỳ canh tác thực tế. Ví dụ, lúa nước trải qua các giai đoạn làm đất ngập nước, đẻ nhánh xanh tốt và thu hoạch để lộ đất trống. Việc tích hợp 35 kênh phổ cùng 7 kênh độ lệch chuẩn trung bình giúp mô hình phân biệt rõ ràng lúa nước với các lớp phủ thực vật khác, nâng độ chính xác toàn diện tăng 10,03% so với ảnh đơn.

Mô hình học kết hợp (Ensemble Learning) giải quyết bài toán phân loại như thế nào?
Mô hình kết hợp 5 thuật toán (XGBoost, Hồi quy Logistic, SVM-RBF, SVM-Linear, MLP) theo cơ chế bỏ phiếu đa số. Kỹ thuật này phát huy thế mạnh nhận diện các lớp diện tích lớn của XGBoost, đồng thời sử dụng khả năng phân tách không gian phi tuyến của SVM để cải thiện vượt bậc độ chính xác tại các lớp có ít mẫu như đất trống (F1 đạt 0,46 so với 0,23 của XGBoost đơn lẻ).

Kết quả phân loại diện tích lúa trong luận văn có độ tin cậy thực tế ra sao?
Kết quả phân loại lớp lúa nước của mô hình Ensemble đạt chỉ số F1 là 0,90. Khi đối chiếu trên toàn địa bàn Hà Nội, diện tích lúa ước tính từ bản đồ phân loại đạt 104.218 ha, sai số chênh lệch chỉ 4.764 ha (tương đương 4,79%) so với số liệu thống kê chính thức 99.454 ha của Cục Thống kê Hà Nội, khẳng định độ chính xác xuất sắc và tính ứng dụng thực tế cao.

Kết luận

  • Luận văn đã giải quyết triệt để thách thức mây che phủ dày đặc trong viễn thám quang học bằng thuật toán ghép ảnh BAP cải tiến, tạo ra chuỗi ảnh Landsat 8 năm 2016 đạt độ sạch mây trên 99%.
  • Chứng minh tính ưu việt vượt trội của việc kết hợp chuỗi dữ liệu đa thời gian (35 kênh phổ và 7 kênh MSD), giúp nâng độ chính xác toàn diện lên hơn 10,03% so với phương pháp phân loại trên ảnh đơn lẻ truyền thống.
  • Xây dựng thành công mô hình học máy kết hợp Ensemble từ 5 thuật toán cơ sở, đạt độ chính xác toàn diện OA 83,91% và hệ số Kappa 0,79, vượt qua tất cả các mô hình phân loại đơn lẻ.
  • Cung cấp bản đồ lớp phủ mặt đất độ phân giải 30 m cho thành phố Hà Nội năm 2016 với sai số ước tính diện tích lúa nước chỉ 4,79% so với số liệu thống kê nhà nước.
  • Mở ra hướng nghiên cứu mới về việc tích hợp thông tin lịch mùa vụ nông nghiệp và kết hợp dữ liệu vệ tinh Radar SAR đa nguồn để tối ưu hóa việc phân tách các lớp phủ phức tạp.

Quý độc giả, các nhà nghiên cứu và cơ quan quản lý quan tâm đến công nghệ xử lý ảnh viễn thám và mô hình học máy ứng dụng có thể khai thác toàn văn tài liệu và mã nguồn thuật toán để triển khai hiệu quả vào các dự án giám sát tài nguyên môi trường trong giai đoạn tới.