Tổng quan nghiên cứu

Biến động đất trồng lúa tại vùng Đồng bằng sông Cửu Long đang diễn ra phức tạp dưới áp lực chuyển đổi cơ cấu cây trồng, đô thị hóa và tác động của biến đổi khí hậu. Theo số liệu thống kê ngành, tổng diện tích đất lúa toàn vùng năm 2010 đạt 4.120,18 nghìn ha, sau đó sụt giảm xuống 3.951,00 nghìn ha vào giai đoạn 2014 - 2015 và phục hồi nhẹ lên 4.030,75 nghìn ha vào cuối năm 2015. Trong bối cảnh đó, tỉnh Hậu Giang với diện tích tự nhiên khoảng 1.600 km2, nằm ở toạ độ từ 9°30'35'' đến 10°19'17'' vĩ độ Bắc, là một trong những trung tâm sản xuất lúa gạo trọng điểm nhưng cũng chịu nhiều biến động mạnh mẽ về cơ cấu sử dụng đất nông nghiệp.

Vấn đề đặt ra cho công tác quản lý tài nguyên là các phương pháp thống kê, kiểm kê truyền thống thường tốn kém chi phí, chu kỳ cập nhật dài và khó xác định chính xác vị trí không gian của các biến động. Đề tài tập trung giải quyết bài toán theo dõi, đánh giá định lượng sự biến động diện tích đất trồng lúa trên địa bàn tỉnh Hậu Giang giai đoạn 2010 – 2015 thông qua việc tích hợp tư liệu viễn thám đa thời gian và công nghệ thông tin địa lý.

Mục tiêu cụ thể của công trình bao gồm: xây dựng quy trình giải đoán chuỗi ảnh viễn thám phản ánh hiện trạng canh tác; thành lập bản đồ hiện trạng và bản đồ biến động đất trồng lúa hai thời điểm 2010 và 2015; định lượng mức độ tăng giảm diện tích tại 8 đơn vị hành chính cấp huyện; đồng thời kiểm chứng độ tin cậy của kết quả giải đoán so với số liệu kiểm kê thực tế của Sở Tài nguyên và Môi trường. Kết quả nghiên cứu xác nhận tổng diện tích đất trồng lúa của tỉnh năm 2015 đã giảm 4,18% so với năm 2010, cung cấp cơ sở dữ liệu không gian chuẩn xác phục vụ công tác quy hoạch sử dụng đất bền vững và đảm bảo an ninh lương thực địa phươ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 dựa trên sự tích hợp giữa lý thuyết phản xạ phổ điện từ của thảm thực vật trong viễn thám quang học và khung phân loại sử dụng đất nông nghiệp theo quy định của Luật Đất đai năm 2013 cùng các hướng dẫn đánh giá đất đai từ Tổ chức Lương thực và Nông nghiệp Liên Hợp Quốc.

Theo nguyên lý viễn thám, mỗi lớp phủ bề mặt có một đặc trưng phản xạ phổ riêng biệt tại các bước sóng khác nhau. Thảm thực vật khỏe mạnh hấp thụ mạnh bức xạ ánh sáng đỏ tại kênh phổ 0,620 – 0,670 μm phục vụ quá trình quang hợp của diệp lục và phản xạ rất mạnh ở dải sóng cận hồng ngoại tại bước sóng 0,841 – 0,876 μm do cấu trúc tế bào lá. Dựa trên đặc tính này, chỉ số khác biệt thực vật chuẩn hóa NDVI được tính toán theo công thức:

$$NDVI = \frac{NIR - RED}{NIR + RED}$$

Trong đó, NIR là giá trị bức xạ phổ cận hồng ngoại và RED là giá trị phổ dải đỏ. Giá trị NDVI biến thiên từ -1 đến +1, phản ánh trực tiếp mật độ sinh khối và chu kỳ sinh trưởng của cây trồng.

Bên cạnh đó, nghiên cứu áp dụng khung phân loại đất nông nghiệp theo quy chuẩn hiện hành, phân định rõ đất chuyên trồng lúa nước (canh tác từ 2 đến 3 vụ/năm) với đất trồng cây hàng năm khác, đất trồng cây lâu năm, đất nuôi trồng thủy sản và đất phi nông nghiệp. Mô hình chuỗi thời gian NDVI đa thời kỳ được sử dụng làm cơ sở để nhận diện đường cong sinh thái mùa vụ, phân tách vụ Đông Xuân, Hè Thu và Thu Đông dựa trên nhịp điệu sinh học của cây lúa nước.

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

Nguồn dữ liệu đầu vào sử dụng chuỗi ảnh vệ tinh MODIS cấp độ sản phẩm MOD09Q1 được cung cấp bởi Cơ quan Hàng không và Vũ trụ Hoa Kỳ với độ phân giải không gian 250 m và độ phân giải thời gian tổ hợp 8 ngày. Cỡ mẫu dữ liệu bao gồm 65 cảnh ảnh liên tục cho chuỗi thời gian 2010 – 2011 và 65 cảnh ảnh cho chuỗi thời gian 2015 – 2016, đảm bảo bao phủ trọn vẹn 3 vụ canh tác lúa trong năm và chu kỳ chuyển tiếp đầu năm sau. Dữ liệu thực địa gồm 40 điểm khảo sát không gian phân bố trên 8 huyện, thị xã, thành phố được thu thập bằng thiết bị định vị vệ tinh GPS cầm tay kết hợp phỏng vấn lịch thời vụ của người dân.

Phương pháp chọn mẫu áp dụng kỹ thuật chọn mẫu ngẫu nhiên phân tầng theo các tiểu vùng sinh thái nông nghiệp đặc trưng, đảm bảo tính đại diện cho vùng chuyên canh lúa thâm canh 3 vụ, vùng lúa - màu luân canh và vùng chuyển đổi sang cây ăn trái hoặc mía đường.

Quy trình phân tích dữ liệu được thực hiện tuần tự qua các bước:

  • Hiệu chỉnh hình học và nắn chuyển hệ tọa độ ảnh từ vĩ độ/kinh độ nguyên bản sang hệ quy chiếu UTM WGS-84 múi chiếu 48 Bắc nhằm đồng nhất với bản đồ địa chính địa phương.
  • Lọc nhiễu ảnh bằng thuật toán trung vị Median kích thước cửa sổ 3x3 nhằm loại bỏ các giá trị dị biệt do mây mù và tán xạ khí quyển mùa mưa.
  • Phân loại không kiểm soát bằng thuật toán ISODATA trên phần mềm ENVI để nhóm tự động các điểm ảnh thành 7 nhóm phổ ban đầu, sau đó hiệu chỉnh thành 6 đối tượng lớp phủ thực tế.
  • Phân loại có kiểm soát dựa trên việc khoanh các vùng quan tâm mẫu đại diện ROI, đối chiếu với chuỗi biến thiên NDVI dạng hình sin và số liệu khảo sát thực địa.
  • Chuyển đổi dữ liệu phân loại định dạng raster sang vector, tích hợp vào phần mềm MapInfo để biên tập bản đồ hiện trạng sử dụng đất lúa tỷ lệ 1:50.000 và tính toán diện tích không gian.

Lý do lựa chọn chuỗi ảnh MOD09Q1 kết hợp bộ công cụ ENVI và MapInfo là nhờ độ phân giải thời gian 8 ngày vượt trội, giúp theo dõi liên tục tiến độ mùa vụ ở quy mô vùng mà không bị gián đoạn bởi mây mù, đồng thời tiết kiệm 70% thời gian và ngân sách so với phương pháp đo đạc truyền thống.

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 không gian và giải đoán chuỗi ảnh viễn thám đa thời gian giai đoạn 2010 – 2015 trên địa bàn tỉnh Hậu Giang đã làm sáng tỏ 4 phát hiện quan trọng:

Thứ nhất, tổng diện tích đất trồng lúa toàn tỉnh có xu hướng suy giảm rõ rệt. Năm 2010, diện tích đất trồng lúa giải đoán được duy trì ở mức cao trên toàn bộ các vùng đồng bằng phù sa màu mỡ. Đến năm 2015, tổng diện tích đất trồng lúa đã giảm 4,18% so với năm 2010, phản ánh sự chuyển dịch đất đai sang các mục đích kinh tế khác.

Thứ hai, biến động diện tích đất lúa diễn ra không đồng đều giữa các địa phương, trong đó huyện Phụng Hiệp ghi nhận mức sụt giảm diện tích lớn nhất. Diện tích đất trồng lúa năm 2015 của huyện Phụng Hiệp giảm tới 6,03% so với năm 2010. Nguyên nhân chính là do người dân chủ động chuyển đổi các diện tích lúa năng suất bấp bênh sang vùng chuyên canh mía nguyên liệu với quy mô hơn 10.300 ha và cải tạo thành các vườn cây ăn trái có múi đặc sản.

Thứ ba, xuất hiện xu hướng tăng diện tích đất lúa cục bộ tại một số tiểu vùng sinh thái, điển hình là huyện Long Mỹ. Trái ngược với xu thế giảm chung của toàn tỉnh, diện tích đất trồng lúa tại huyện Long Mỹ năm 2015 đã tăng 9,61% so với năm 2010. Sự gia tăng này bắt nguồn từ các dự án hoàn thiện hệ thống đê bao thủy lợi khép kín, ngăn mặn từ triều cường biển Tây, cho phép chuyển đổi diện tích đất hoang hóa, đất một vụ lúa mùa kém hiệu quả sang thâm canh lúa 2 đến 3 vụ ổn định.

Thứ tư, đồ thị chuỗi thời gian NDVI đã bóc tách rõ nét nhịp điệu canh tác nông nghiệp. Vùng đất trồng lúa 3 vụ thể hiện rõ đường cong biến thiên NDVI dạng hình sin với 3 chu kỳ hoàn chỉnh trong năm. Giá trị NDVI đạt mức cực tiểu từ 0,15 đến 0,25 vào các giai đoạn làm đất, sạ cấy ngập nước hoặc ngay sau thu hoạch; sau đó tăng vọt lên mức đỉnh từ 0,75 đến 0,88 tại giai đoạn lúa đẻ nhánh làm đòng rộ; và hạ dần khi lúa chín vàng. Ngược lại, các vùng cây lâu năm và vườn cây ăn trái luôn duy trì chỉ số NDVI ổn định ở mức cao từ 0,60 đến 0,85 quanh năm mà không có chu kỳ dao động sâu.

Thảo luận kết quả

Sự biến động giảm diện tích đất lúa tại Hậu Giang trong giai đoạn 2010 – 2015 phản ánh đúng định hướng chuyển dịch cơ cấu nông nghiệp sang nâng cao giá trị gia tăng. Tỉnh đã hình thành các vùng chuyên canh cây có múi trên 10.000 ha, vùng khóm Cầu Đúc 1.500 ha và vùng nuôi trồng thủy sản tập trung 1.400 ha. Đồng thời, việc mở rộng hệ thống giao thông liên tỉnh như Quốc lộ 61, tuyến nối Vị Thanh – Cần Thơ cùng các cụm công nghiệp tại huyện Châu Thành và Châu Thành A đã thu hồi một phần diện tích đất lúa ven trục lộ.

Khi so sánh với các nghiên cứu ứng dụng viễn thám giám sát lũ và thảm thực vật tại vùng hạ lưu sông Mê Kông, kết quả nghiên cứu này tương đồng cao về mặt quy luật phản xạ phổ của cây lúa nước. Việc sử dụng sản phẩm tổ hợp 8 ngày MOD09Q1 đã giải quyết triệt để bài toán thiếu dữ liệu do mây mù nhiệt đới – vốn là rào cản lớn khi dùng ảnh vệ tinh Landsat chu kỳ 16 ngày.

Dữ liệu nghiên cứu được biểu diễn trực quan thông qua:

  • Bản đồ hiện trạng sử dụng đất lúa năm 2010 và 2015 thể hiện không gian phân bố cụ thể của từng loại hình canh tác lúa 2 vụ, lúa 3 vụ và lúa - màu.
  • Biểu đồ hình cột so sánh biến động diện tích đất lúa giữa hai mốc thời gian của từng huyện, thị xã, thành phố.
  • Bảng ma trận sai số phân loại đất đai với hệ số chính xác tổng thể đạt 86,4% và hệ số kiểm định Kappa đạt 0,82, khẳng định độ tin cậy học thuật cao.
  • Biểu đồ tương quan tuyến tính giữa diện tích giải đoán từ viễn thám và số liệu kiểm kê đất đai của Sở Tài nguyên và Môi trường tỉnh Hậu Giang cho thấy hệ số tương quan đạt mức xấp xỉ 0,92, chứng minh tính khả thi của việc ứng dụng công nghệ viễn thám trong kiểm kê đất đai định kỳ.

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

Dựa trên kết quả nghiên cứu biến động đất trồng lúa giai đoạn 2010 – 2015, luận văn đưa ra 4 nhóm giải pháp mang tính ứng dụng thực tiễn cao:

Thứ nhất, thiết lập hệ thống giám sát tự động biến động đất nông nghiệp bằng công nghệ viễn thám và GIS định kỳ 6 tháng một lần. Đơn vị chủ trì là Sở Tài nguyên và Môi trường phối hợp cùng Sở Nông nghiệp và Phát triển nông thôn tỉnh Hậu Giang, triển khai hoàn thiện trong giai đoạn 2026 – 2028. Mục tiêu là kiểm soát 100% diện tích đất lúa bị chuyển đổi tự phát sang nuôi thủy sản hoặc chuyển nhượng trái phép, giảm thiểu sai số kiểm kê từ phương pháp báo cáo hành chính thủ công.

Thứ hai, quy hoạch khoanh vùng và bảo vệ nghiêm ngặt 32.000 ha vùng nguyên liệu lúa chất lượng cao, trọng tâm tại hai cánh đồng mẫu lớn quy mô 300 - 500 ha thuộc huyện Vị Thủy và Châu Thành A. Ủy ban nhân dân tỉnh Hậu Giang cùng các phòng ban chuyên môn cấp huyện cần ban hành quy chế quản lý phân vùng sử dụng đất từ năm 2026, duy trì diện tích gieo trồng lúa ổn định ở mức trên 85.000 ha toàn tỉnh nhằm bảo đảm an ninh lương thực và đáp ứng nhu cầu gạo xuất khẩu.

Thứ ba, nâng cấp cơ sở dữ liệu không gian đất đai tích hợp trên nền tảng WebGIS dùng chung. Trung tâm Công nghệ thông tin Tài nguyên Môi trường cần số hóa và đồng bộ 100% các lớp dữ liệu bản đồ hiện trạng, bản đồ biến động đất lúa định dạng vector tỷ lệ 1:25.000 trước năm 2027. Việc này giúp các cơ quan quản lý tại 8 huyện, thị xã, thành phố dễ dàng truy vấn thuộc tính, theo dõi biến động thời gian thực và hỗ trợ cấp phép chuyển mục đích sử dụng đất chính xác.

Thứ tư, đẩy mạnh ứng dụng tư liệu vệ tinh đa phổ thế hệ mới kết hợp viễn thám Radar khẩu độ tổng hợp. Khuyến nghị các viện, trường và cơ quan chuyên môn tích hợp dữ liệu ảnh Sentinel-1 và Sentinel-2 độ phân giải 10 m nhằm khắc phục hoàn toàn hiện tượng che phủ mây trong mùa mưa bão, nâng độ chính xác phân loại diện tích đất lúa lên mức trên 92% trong các kỳ kiểm kê đất đai định kỳ 5 năm tiếp theo.

Đố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ị thực tiễn cho 4 nhóm đối tượng chính:

Nhóm thứ nhất là cán bộ quản lý 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, Văn phòng Đăng ký đất đai và Phòng Tài nguyên Môi trường cấp huyện. Luận văn cung cấp quy trình công nghệ chuẩn để theo dõi biến động lớp phủ đất đai, giúp rút ngắn thời gian điều tra thống kê hằng năm và nâng cao độ chính xác của hồ sơ địa chính.

Nhóm thứ hai là các chuyên gia, viện nghiên cứu và đơn vị tư vấn quy hoạch đất đai. Tài liệu cung cấp cơ sở dữ liệu chuỗi thời gian NDVI thực tế và phương pháp luận tích hợp ENVI - MapInfo để xây dựng các phương án phân vùng sử dụng đất nông nghiệp thích ứng với biến đổi khí hậu tại vùng Đồng bằng sông Cửu Long.

Nhóm thứ ba là nghiên cứu sinh, học viên cao học và sinh viên các chuyên ngành Quản lý đất đai, Địa tin học, Trắc địa - Bản đồ, Khoa học Môi trường và Nông nghiệp. Luận văn đóng vai trò là tài liệu tham khảo chi tiết về thuật toán xử lý ảnh số, kỹ thuật khoanh vùng mẫu ROI, phân loại ISODATA và đánh giá ma trận sai số.

Nhóm thứ tư là các tổ chức phát triển nông nghiệp, hiệp hội lương thực và doanh nghiệp đầu tư cánh đồng mẫu lớn. Doanh nghiệp có thể ứng dụng phương pháp giải đoán chuỗi ảnh viễn thám để ước tính tiến độ xuống giống, dự báo mùa vụ thu hoạch và quản lý vùng nguyên liệu lúa chất lượng cao với quy mô tập trung.

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

Tại sao ảnh vệ tinh MODIS MOD09Q1 lại được ưu tiên lựa chọn để theo dõi lúa nước thay vì ảnh Landsat hay SPOT?

Ảnh MOD09Q1 sở hữu độ phân giải thời gian tổ hợp 8 ngày vượt trội, cung cấp chuỗi dữ liệu liên tục 46 đến 65 cảnh ảnh mỗi năm. Điều này cho phép theo dõi sát sao từng giai đoạn sinh trưởng ngắn ngày của cây lúa nước. Trong khi đó, ảnh Landsat có chu kỳ lặp 16 ngày thường xuyên bị mây mù che phủ trong mùa mưa bão tại Đồng bằng sông Cửu Long, gây gián đoạn chuỗi dữ liệu theo dõi.

Chỉ số NDVI bóc tách chu kỳ sinh trưởng của các mùa vụ lúa như thế nào?

Chỉ số NDVI biến động theo chu kỳ sinh học của cây lúa: giá trị thấp khoảng 0,10 đến 0,25 khi ruộng ngập nước làm đất hoặc sau thu hoạch; tăng nhanh và đạt đỉnh cực đại từ 0,75 đến 0,90 khi lúa đẻ nhánh, làm đòng và trổ bông; sau đó hạ thấp khi lúa chín. Dựa vào số lượng đỉnh hình sin trong năm, nhà nghiên cứu phân biệt chính xác khu vực lúa 3 vụ, lúa 2 vụ hay luân canh lúa - màu.

Thuật toán lọc trung vị Median cửa sổ 3x3 xử lý sai số dữ liệu ảnh như thế nào?

Toán tử lọc Median 3x3 dịch chuyển qua từng điểm ảnh, lấy giá trị trung vị của 9 pixel lân cận để thay thế cho pixel trung tâm. Phương pháp này loại bỏ hiệu quả các điểm ảnh nhiễu do tán xạ khí quyển, đám mây nhỏ hoặc lỗi đường truyền vô tuyến mà vẫn bảo toàn nguyên vẹn các đường ranh giới và cấu trúc địa vật của thửa đất.

Độ chính xác của kết quả giải đoán đất lúa từ viễn thám so với số liệu kiểm kê thực tế đạt bao nhiêu?

Kết quả kiểm chứng ma trận sai số phân loại đạt độ chính xác tổng thể 86,4% cùng chỉ số Kappa đạt 0,82. Khi đối chiếu với số liệu kiểm kê đất đai năm 2010 và 2015 của Sở Tài nguyên và Môi trường tỉnh Hậu Giang, hệ số tương quan tuyến tính đạt trên 0,92, khẳng định tính chuẩn xác và độ tin cậy khoa học cao của phương pháp.

Quy trình tích hợp giữa phần mềm xử lý ảnh ENVI và GIS MapInfo được thực hiện ra sao?

Phần mềm ENVI đảm nhận toàn bộ khâu tiền xử lý, hiệu chỉnh hình học, tính toán chỉ số NDVI và phân loại có kiểm soát chuỗi ảnh đa phổ. Kết quả sau phân loại được xuất sang định dạng vector dạng tệp SHP hoặc MIF, sau đó nạp vào MapInfo để biên tập bản đồ chuyên đề, kết nối bảng dữ liệu thuộc tính và truy vấn diện tích từng thửa đất.

Kết luận

  • Luận văn đã xây dựng thành công quy trình ứng dụng chuỗi ảnh viễn thám MODIS MOD09Q1 kết hợp phần mềm ENVI 4.8 và MapInfo 10.5 để giám sát biến động đất trồng lúa đa thời gian tại tỉnh Hậu Giang.
  • Kết quả giải đoán định lượng xác nhận tổng diện tích đất trồng lúa toàn tỉnh giai đoạn 2010 – 2015 suy giảm 4,18%, phản ánh sự chuyển dịch sang đất chuyên canh cây ăn quả, mía đường và đất phi nông nghiệp.
  • Huyện Phụng Hiệp ghi nhận mức giảm diện tích đất lúa sâu nhất với tỷ lệ 6,03%, trong khi huyện Long Mỹ tăng trưởng 9,61% diện tích nhờ các công trình thủy lợi ngăn mặn khép kín.
  • Kiểm định ma trận sai số phân loại đạt độ chính xác 86,4% và hệ số tương quan đạt trên 0,92 so với số liệu kiểm kê của Sở Tài nguyên và Môi trường, khẳng định tính khả thi vượt trội của công nghệ viễn thám.
  • Đóng góp chính của đề tài là cung cấp bộ bản đồ hiện trạng không gian chuẩn xác và giải pháp công nghệ chi phí thấp, hỗ trợ đắc lực cho công tác quản lý đất đai địa phương.

Kế hoạch tiếp theo cần tập trung chuyển giao quy trình phân tích cho cán bộ địa chính cấp huyện và tích hợp nguồn dữ liệu vệ tinh độ phân giải cao Sentinel trước năm 2028. Các cơ quan quản lý và nhà nghiên cứu quan tâm có thể tiếp cận, ứng dụng toàn bộ khung dữ liệu này để phục vụ công tác quy hoạch sử dụng đất bền vững ngay từ hôm nay.