Tổng quan nghiên cứu

Lúa gạo là nguồn lương thực thiết yếu bảo đảm an ninh lương thực toàn cầu, trong đó khoảng 90% diện tích canh tác (tương đương 140 triệu hecta) tập trung tại khu vực châu Á. Tại Việt Nam, Đồng bằng sông Hồng là một trong những vựa lúa trọng điểm với quy mô diện tích tự nhiên khoảng 15.000 km2 trải dài trên 11 tỉnh thành phố. Hàng năm, khu vực này đóng góp sản lượng lương thực ổn định trên 6,3 triệu tấn với hai vụ lúa chính gồm vụ Đông Xuân và vụ Hè Thu, trong đó diện tích lúa vụ Đông Xuân chiếm từ 80% đến 90% tổng diện tích gieo cấy lúa tại các tỉnh miền Bắc.

Tuy nhiên, đặc thù khí hậu nhiệt đới gió mùa tại Đồng bằng sông Hồng thường xuyên xuất hiện mưa bão, lũ lụt và độ che phủ mây dày đặc quanh năm. Thực trạng này gây cản trở nghiêm trọng cho các phương pháp kiểm kê nông nghiệp truyền thống cũng như các nghiên cứu sử dụng ảnh vệ tinh quang học đơn thuần. Việc thiếu hụt dữ liệu không gian cập nhật kịp thời làm giảm hiệu lực ra quyết định chỉ đạo sản xuất của các cơ quan quản lý nhà nước và hạn chế năng lực ứng phó thiên tai của các hộ nông dân.

Nhằm khắc phục những khoảng trống kỹ thuật này, nghiên cứu đặt mục tiêu xây dựng quy trình ứng dụng học máy hiện đại trên nền tảng viễn thám đa nguồn, kết hợp giữa ảnh quang học Landsat 8 và ảnh radar khẩu độ tổng hợp Sentinel-1A. Phạm vi nghiên cứu bao quát toàn bộ 11 tỉnh thành Đồng bằng sông Hồng trong giai đoạn từ năm 2015 đến năm 2018. Kết quả nghiên cứu có ý nghĩa thực tiễn to lớn khi thiết lập phương pháp giám sát lúa liên tục với độ chính xác tổng quát trên 91% và thu hẹp sai số diện tích so với số liệu thống kê xuống mức 7,06%, phục vụ đắc lực cho công tác quản trị nông nghiệp số và giảm nhẹ rủi ro thiên tai.

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ừ của công nghệ viễn thám quang học và viễn thám radar khẩu độ tổng hợp. Viễn thám quang học Landsat 8 sử dụng hai bộ cảm biến OLI và TIRS với 11 kênh phổ, độ phân giải không gian 30 mét và chu kỳ chụp lặp 16 ngày. Viễn thám radar Sentinel-1A hoạt động ở băng tần C với chu kỳ chụp lặp từ 6 đến 12 ngày, cung cấp dữ liệu phân cực kép VV và VH với độ phân giải không gian 5x20 mét, sở hữu khả năng xuyên mây vượt trội giúp quan trắc bề mặt Trái Đất liên tục trong mọi điều kiện thời tiết.

Về mặt hiện tượng học thực vật, chu kỳ sinh trưởng của cây lúa từ giai đoạn cấy ngập nước, đẻ nhánh, làm đòng, trổ bông đến khi chín được phản ánh rõ nét qua các chỉ số thực vật chuyên biệt. Các chỉ số này bao gồm chỉ số khác biệt thực vật chuẩn hóa NDVI, chỉ số thực vật tăng cường EVI và chỉ số nước bề mặt LSWI.

Về mặt thuật toán, nghiên cứu ứng dụng hai mô hình học máy giám sát chủ đạo:

  • Thuật toán máy vectơ hỗ trợ SVM: Xác định siêu phẳng phân tách tối ưu trong không gian nhiều chiều thông qua hàm nhân phi tuyến kernel, đặc biệt phù hợp cho bài toán phân lớp nhị phân lúa và phi lúa.
  • Thuật toán cây quyết định tăng cường XGBoost: Tối ưu hóa hàm mục tiêu có thành phần chính quy hóa, hỗ trợ xử lý xuất sắc các tập dữ liệu có độ thưa lớn do điểm ảnh bị mây che phủ.

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

Nghiên cứu thiết lập quy trình tiền xử lý và phân tích dữ liệu tự động trên nền tảng điện toán đám mây Google Earth Engine. Cỡ mẫu và phương pháp chọn mẫu được thiết kế nghiêm ngặt theo từng loại ảnh:

  • Đối với ảnh Landsat 8: Áp dụng phương pháp lấy mẫu phân tầng với tổng cộng 1.277 điểm mẫu huấn luyện (gồm 530 điểm lúa và 747 điểm phi lúa) cùng 378 điểm mẫu kiểm tra độc lập (gồm 108 điểm lúa và 270 điểm phi lúa), được gán nhãn chính xác dựa trên ảnh vệ tinh độ phân giải cao Google Earth.
  • Đối với ảnh Sentinel-1A: Triển khai phương pháp lấy mẫu tích hợp đa cấp (phân tầng cấp tỉnh, chọn cụm cấp huyện và xã, chọn ngẫu nhiên cấp điểm ảnh). Tập dữ liệu huấn luyện đạt quy mô 37.680 điểm (gồm 23.794 điểm lúa và 13.886 điểm phi lúa), tập kiểm tra đạt 1.398 điểm (gồm 549 điểm lúa và 849 điểm phi lúa).

Lý do lựa chọn phương pháp phân tích: Thuật toán SVM được lựa chọn cho dữ liệu radar nhờ khả năng phân tách ranh giới phi tuyến mạnh mẽ trên chuỗi thời gian phân cực kép. Thuật toán XGBoost được áp dụng cho ảnh Landsat 8 nhằm giải quyết bài toán dữ liệu khuyết thiếu khi ghép ảnh trung bình 12 tháng với 84 đặc trưng phổ.

Quy trình tiền xử lý ảnh Sentinel-1A bao gồm bốn bước: cập nhật quỹ đạo chính xác, loại bỏ nhiễu nhiệt, hiệu chỉnh bức xạ và hiệu chỉnh địa hình bằng mô hình số độ cao SRTM 30 mét. Chuỗi thời gian nghiên cứu tập trung vào vụ Đông Xuân từ ngày 1 tháng 2 đến ngày 30 tháng 5 qua các năm 2015, 2016 và 2018.

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

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

Thứ nhất, mô hình học máy XGBoost áp dụng cho ảnh quang học Landsat 8 đạt độ chính xác tổng quát là 90,5% vào năm 2015 và 89,4% vào năm 2016. Độ hồi tưởng đối với lớp lúa đạt mức cao từ 91% đến 94%. Tuy nhiên, mức độ sai khác diện tích ước tính từ bản đồ phân lớp so với số liệu thống kê của Tổng cục Thống kê có sự dao động lớn, từ 7,06% (chênh lệch 38.261 hecta năm 2015) tăng lên 15,42% (chênh lệch 80.554 hecta năm 2016) do ảnh hưởng của lượng mây che phủ dày đặc trong năm 2016.

Thứ hai, mô hình học máy SVM kết hợp tập đặc trưng đa thời gian phân cực kép VV và VH từ Sentinel-1A mang lại hiệu năng cao nhất với độ chính xác tổng quát đạt 91,2%, hệ số Kappa đạt 0,81 và chỉ số F1-Score đạt 0,89. Độ chính xác của lớp lúa đạt 91% và độ hồi tưởng đạt 84%, vượt trội hoàn toàn so với mô hình chỉ sử dụng kênh đơn phân cực VV (độ chính xác tổng quát 82,4%) hoặc kênh đơn phân cực VH (độ chính xác tổng quát 88,1%).

Thứ ba, việc sử dụng chỉ số tỷ lệ phân cực chuẩn hóa (VV-VH)/(VV+VH) mang lại độ chính xác 83% cho lớp phi lúa, hỗ trợ phân biệt hiệu quả giữa các cánh đồng lúa nước với bề mặt đầm lầy, thảm thực vật tự nhiên và vùng nuôi trồng thủy sản ven biển tại Quảng Ninh và Nam Định.

Thứ tư, phương pháp phân lớp lúa theo từng tháng cho phép bóc tách chi tiết thời điểm gieo cấy lúa sớm (tháng 2), lúa chính vụ (tháng 3, tháng 4) và lúa thu hoạch muộn (tháng 5), tạo tiền đề cho việc xây dựng bản đồ canh tác lúa động theo chu kỳ sinh trưởng.

Thảo luận kết quả

Hiệu năng vượt trội của ảnh radar Sentinel-1A được giải thích bởi cơ chế tán xạ sóng radar thay đổi nhịp nhàng theo hình thái học cây lúa. Trong giai đoạn cấy ngập nước, mặt ruộng phẳng tạo ra hiện tượng phản xạ gương làm giảm thiểu hệ số tán xạ ngược. Khi cây lúa phát triển qua giai đoạn đẻ nhánh và ôm đòng, sinh khối thân lá tăng vọt dẫn đến hiện tượng tán xạ khối mạnh mẽ, làm gia tăng đột biến giá trị tán xạ ngược ở kênh phân cực chéo VH.

Về phương diện trực quan hóa, toàn bộ động thái phát triển mùa vụ được mô hình hóa chi tiết qua đồ thị đường biểu diễn giá trị tán xạ ngược (đơn vị decibel) theo thời gian từ tháng 2 đến tháng 5. Dữ liệu phân loại được tổng hợp trực quan thành bản đồ chuyên đề độ phân giải 10x10 mét phủ kín 11 tỉnh thành và bảng ma trận nhầm lẫn chuẩn hóa, minh chứng rõ ràng tỷ lệ phân loại đúng và giảm thiểu tối đa các điểm nhận diện sai giữa lúa và hoa màu ngắn ngày.

So với các nghiên cứu trước đây vốn chủ yếu dựa vào việc phân ngưỡng chỉ số EVI trên ảnh MODIS độ phân giải thô 250-500 mét tại Đồng bằng sông Cửu Long, phương pháp đề xuất trong luận văn đã giải quyết triệt để bài toán mây mù tại miền Bắc, nâng cao độ phân giải không gian lên 10 mét và thiết lập khung giám sát định kỳ tự động hóa hoàn toàn.

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

Nhằm chuyển giao và phát huy tối đa giá trị khoa học của đề tài vào thực tiễn quản lý sản xuất nông nghiệp, các giải pháp cụ thể được đề xuất như sau:

  • Tự động hóa hệ thống tiền xử lý dữ liệu viễn thám: Xây dựng đường ống tự động thu nạp, lọc nhiễu nhiệt, hiệu chỉnh địa hình và trích xuất đặc trưng chuỗi thời gian ảnh Sentinel-1A trên nền tảng Google Earth Engine. Mục tiêu rút ngắn 85% thời gian xử lý dữ liệu, hoàn thành xuất bản đồ sơ bộ trong vòng 24 giờ sau khi vệ tinh chụp ảnh, do Trung tâm Công nghệ Giám sát hiện trường FIMO chủ trì thực hiện từ quý 1 năm 2026.
  • Ứng dụng mô hình SVM đa phân cực trong dự báo mùa vụ định kỳ: Triển khai công cụ giám sát lúa hàng tháng phục vụ Cục Trồng trọt và Chi cục Bảo vệ thực vật tại 11 tỉnh Đồng bằng sông Hồng. Mục tiêu kiểm soát sai số ước tính diện tích gieo trồng dưới mức 5% so với số liệu thực địa trước quý 3 năm 2026.
  • Xây dựng hệ sinh thái cảnh báo và đánh giá ngập úng khẩn cấp: Tích hợp bản đồ lúa viễn thám 10 mét với dữ liệu lượng mưa và mực nước thực thời để định lượng diện tích lúa bị úng ngập sau bão lũ. Ban Chỉ đạo Quốc gia về Phòng chống thiên tai phối hợp cùng các viện nghiên cứu đưa hệ thống vào vận hành thường trực trước mùa mưa bão năm 2026, hỗ trợ ban hành phương án cứu trợ nông dân trong vòng 48 giờ.
  • Chuẩn hóa cơ sở dữ liệu mẫu thực địa quốc gia: Bộ Nông nghiệp và Phát triển Nông thôn phối hợp với các trường đại học triển khai thu thập tối thiểu 100.000 điểm mẫu mặt đất chuẩn hóa trên phạm vi toàn quốc, cập nhật thường niên vào kho dữ liệu lớn phục vụ huấn luyện các thuật toán học máy và học sâu giai đoạn 2026 - 2028.

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

Nội dung và kết quả thực nghiệm của luận văn mang lại giá trị học thuật và ứng dụng chuyên sâu cho bốn nhóm đối tượng:

  • Học viên cao học, nghiên cứu sinh và nhà khoa học chuyên ngành Hệ thống Thông tin, Địa tin học (GIS), Khoa học Dữ liệu và Viễn thám: Tài liệu cung cấp cơ sở lý thuyết hoàn chỉnh về thuật toán SVM, XGBoost, kỹ thuật xử lý ảnh radar chuỗi thời gian và phương thức lập trình trên Google Earth Engine.
  • Cán bộ quản lý tại Bộ Nông nghiệp và Phát triển Nông thôn, Tổng cục Thống kê: Vận dụng quy trình lập bản đồ viễn thám để số hóa công tác thống kê diện tích gieo cấy, theo dõi tiến độ sản xuất lúa và tối ưu hóa chuỗi cung ứng vật tư nông nghiệp.
  • Chuyên viên tại các cơ quan Khí tượng Thủy văn và Phòng chống Thiên tai: Ứng dụng bản đồ phân bố lúa để đánh giá rủi ro ngập lụt, khoanh vùng thiên tai và xây dựng kịch bản điều tiết hồ chứa, bảo vệ mùa màng.
  • Doanh nghiệp công nghệ nông nghiệp và tổ chức tài chính - bảo hiểm nông nghiệp: Khai thác dữ liệu phân lớp lúa độ phân giải cao để thẩm định diện tích bảo hiểm mùa màng, dự báo sản lượng thu mua và quản lý rủi ro tín dụng nông thôn.

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

Ảnh vệ tinh radar Sentinel-1A có ưu điểm vượt trội gì so với ảnh quang học Landsat 8 tại khu vực Đồng bằng sông Hồng? Ảnh radar Sentinel-1A sử dụng sóng vi ba băng tần C có khả năng xuyên qua mây mù, mưa phùn và hoạt động liên tục cả ngày lẫn đêm. Trong điều kiện Đồng bằng sông Hồng thường xuyên bị mây che phủ với tỷ lệ khuyết ảnh quang học rất cao, Sentinel-1A bảo đảm nguồn dữ liệu ổn định với chu kỳ lặp từ 6 đến 12 ngày, khắc phục triệt để khoảng trống dữ liệu 16 ngày của Landsat 8.

Sự kết hợp giữa hai kênh phân cực VV và VH mang lại lợi ích gì cho bài toán phân lớp lúa? Kênh phân cực kép kết hợp giữa VV và VH giúp mô hình đạt độ chính xác tổng quát cao nhất là 91,2%. Trong đó, kênh phân cực chéo VH cực kỳ nhạy cảm với cấu trúc hình thái và tán xạ khối của thân lá cây lúa, còn kênh cùng phân cực VV phản ánh rõ nét tương tác giữa mặt nước ngập và gốc lúa, giúp tách biệt hoàn toàn lúa với các thảm phủ thực vật khác.

Thuật toán học máy XGBoost phát huy vai trò như thế nào khi phân loại ảnh quang học? XGBoost là thuật toán tăng cường độ dốc tối ưu hóa việc xây dựng các cây quyết định tuần tự. Khi xử lý ảnh quang học Landsat 8 tổng hợp 12 tháng với 84 đặc trưng phổ, XGBoost có khả năng xử lý hiệu quả các ma trận dữ liệu có nhiều điểm ảnh bị khuyết do lọc mây, giúp duy trì độ chính xác tổng quát trên 89,4%.

Phương pháp chọn mẫu trong luận văn được thiết kế thế nào để tránh thiên lệch dữ liệu? Đối với dữ liệu Sentinel-1A, tác giả áp dụng quy trình chọn mẫu tích hợp đa cấp gồm phân tầng theo 11 tỉnh thành, chọn cụm tại các huyện/xã trọng điểm và chọn mẫu ngẫu nhiên trên bản đồ. Tập huấn luyện với 37.680 điểm và tập kiểm tra với 1.398 điểm được phân tách độc lập hoàn toàn, bảo đảm tính đại diện khách quan cho toàn bộ tiểu vùng sinh thái.

Nguyên nhân chính dẫn đến sai khác giữa diện tích ước tính từ viễn thám và số liệu của Tổng cục Thống kê là gì? Sự chênh lệch từ 7,06% đến 15,42% bắt nguồn từ đặc thù ruộng lúa manh mún, bờ vùng bờ thửa nhỏ hẹp tại miền Bắc tạo ra hiện tượng điểm ảnh hỗn hợp (mixed pixel). Bên cạnh đó, tình trạng xen canh hoa màu ngắn ngày cùng độ trễ thời gian trong phương pháp báo cáo thống kê hành chính thủ công tại các địa phương cũng tạo ra sai số đối sánh.

Kết luận

  • Luận văn đã chứng minh tính khả thi và hiệu quả vượt bậc của việc ứng dụng học máy kết hợp dữ liệu viễn thám radar Sentinel-1A trong việc phân lớp và giám sát lúa liên tục tại Đồng bằng sông Hồng.
  • Mô hình máy vectơ hỗ trợ SVM trên tập dữ liệu đa thời gian phân cực kép VV và VH đạt độ chính xác tổng quát 91,2%, hệ số Kappa 0,81 và độ chính xác lớp lúa đạt 91%.
  • Quy trình tự động hóa ghép ảnh và trích xuất đặc trưng trên nền tảng Google Earth Engine đã giải quyết triệt để rào cản mây mù của ảnh vệ tinh quang học truyền thống.
  • Đóng góp tập dữ liệu chuẩn hóa quy mô lớn với hơn 39.000 điểm mẫu thực địa phủ kín 11 tỉnh thành, cung cấp tài liệu tham khảo khoa học có giá trị cao cho ngành Hệ thống thông tin và Địa tin học.
  • Đặt nền móng vững chắc cho các nghiên cứu tiếp theo mở rộng sang kiến trúc học sâu Deep Learning và mạng nơ-ron tích chập đa thời gian giai đoạn 2026 - 2030.

Các cơ quan quản lý nông nghiệp, viện nghiên cứu và doanh nghiệp công nghệ hãy chủ động tiếp cận và ứng dụng giải pháp viễn thám học máy này nhằm nâng cao hiệu quả điều hành sản xuất và thúc đẩy tiến trình chuyển đổi số nông nghiệp quốc gia.