Tổng quan nghiên cứu

Nghiên cứu về biến động sử dụng đất đóng vai trò then chốt trong công tác quản trị tài nguyên và quy hoạch không gian bền vững. Tại huyện Thạch Thành, tỉnh Thanh Hóa, tổng diện tích tự nhiên đạt 55.921,72 ha, trong đó đất nông nghiệp chiếm ưu thế với 46.200,88 ha (chiếm 82,62%), đất phi nông nghiệp chiếm 8.221,35 ha (chiếm 14,70%) và đất chưa sử dụng là 1.499,49 ha (chiếm 2,68%). Trước áp lực từ quá trình công nghiệp hóa, sự hình thành các cụm công nghiệp Vân Du, Thạch Quảng cùng mạng lưới hạ tầng huyết mạch như Quốc lộ 45 dài 20,5 km và đường Hồ Chí Minh dài 13,3 km, cơ cấu lớp phủ mặt đất tại địa phương đã diễn ra nhiều biến động phức tạp. Các phương pháp thống kê và kiểm kê đất đai truyền thống thường có chu kỳ dài từ 5 đến 10 năm, tốn kém kinh phí và khó xác định trực quan vị trí không gian biến động. Đề tài luận văn thạc sĩ chuyên ngành Quản lý đất đai đã ứng dụng tích hợp công nghệ Hệ thống thông tin địa lý (GIS) và viễn thám vệ tinh đa thời gian nhằm đánh giá chi tiết biến động sử dụng đất giai đoạn 2014 - 2017 trên địa bàn 28 xã, thị trấn của huyện Thạch Thành. Kết quả nghiên cứu xây dựng thành công bộ dữ liệu số với độ chính xác phân loại tổng thể đạt 86,00% năm 2014 và 91,00% năm 2017, chỉ số Kappa tương ứng là 0,83 và 0,89. Công trình mang giá trị thực tiễn cao, giúp rút ngắn hơn 60% thời gian giám sát thực địa, tạo lập cơ sở khoa học vững chắc phục vụ điều chỉnh quy hoạch sử dụng đất cấp huyện giai đoạn tiếp theo.

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 hai trụ cột lý thuyết cốt lõi là lý thuyết tương tác bức xạ điện từ trong viễn thám đa phổ và lý thuyết phân tích không gian địa lý trong GIS. Theo nguyên lý phản xạ phổ, mỗi đối tượng bề mặt trái đất như thực vật, mặt nước, đất trống và công trình xây dựng đều sở hữu đặc trưng phản xạ riêng biệt tại từng dải sóng điện từ. Sự kết hợp giữa kênh phổ cận hồng ngoại (NIR) và dải sóng ánh sáng đỏ (Red) qua chỉ số thực vật NDVI cho phép nhận diện chính xác trạng thái sinh trưởng của thảm phủ thực vật. Bên cạnh đó, mô hình phân tích không gian GIS ứng dụng nguyên lý đại số bản đồ raster và vector, hỗ trợ tích hợp đa lớp thông tin địa lý và thiết lập ma trận chuyển dịch trạng thái sử dụng đất (LUCC) giữa các thời kỳ. Ba khái niệm chuyên ngành nền tảng được chuẩn hóa gồm: Phân loại ảnh có kiểm định (Supervised Classification), Ma trận nhầm lẫn (Confusion Matrix) và Chỉ số Kappa (Kappa Coefficient) đánh giá độ chính xác phân loại ảnh số.

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

Nguồn dữ liệu sơ cấp được thu thập từ 100 điểm khảo sát thực địa bằng thiết bị GPS cầm tay phân bố trên 28 xã, thị trấn để xây dựng tệp mẫu huấn luyện và đối soát kết quả giải đoán. Phương pháp chọn mẫu áp dụng là chọn mẫu ngẫu nhiên phân tầng đại diện cho 6 loại hình sử dụng đất chính trên địa bàn. Nguồn dữ liệu thứ cấp gồm ảnh vệ tinh Landsat 8 OLI/TIRS (mã cảnh 127-046) thu nhận ngày 19/09/2014 và ngày 08/09/2017 với độ phân giải không gian 30 m, kết hợp số liệu kiểm kê đất đai năm 2014 và thống kê đất đai năm 2017 của địa phương. Quy trình xử lý ảnh số được thực hiện trên phần mềm ENVI: nắn chỉnh hình học về hệ tọa độ WGS84 với sai số trung phương dưới 0,5 pixel, tăng cường chất lượng ảnh và phân loại có kiểm định theo thuật toán Xác suất cực đại (Maximum Likelihood). Thuật toán này được lựa chọn vì tính toán xác suất dựa trên phân bố chuẩn Gauss của vector độ sáng, mang lại độ tin cậy quang phổ tối ưu nhất cho ảnh Landsat. Sau đó, chức năng phân tích không gian của phần mềm ArcGIS được sử dụng để chồng xếp bản đồ hai thời kỳ 2014 và 2017 nhằm chiết xuất ma trận biến động và biên tập bản đồ biến động đất đai tỷ lệ 1:25.000.

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 loại và kiểm định ảnh vệ tinh Landsat 8 đã chiết xuất thành công hiện trạng đất đai huyện Thạch Thành thành 6 nhóm sử dụng chính với độ tin cậy khoa học rất cao. Năm 2014, độ chính xác phân loại tổng thể đạt 86,00% cùng chỉ số Kappa là 0,83; đến năm 2017, độ chính xác đạt 91,00% cùng chỉ số Kappa là 0,89. Độ chính xác năm 2017 cao hơn nhờ thời gian bay chụp vệ tinh gần sát thời điểm điều tra thực địa 100 điểm mẫu.

Giai đoạn 2014 - 2017 ghi nhận sự chuyển dịch diện tích cụ thể giữa các loại hình đất:

  • Đất nông nghiệp ngắn ngày: Diện tích đất trồng cây hàng năm khác giảm mạnh nhất với 54,36 ha; diện tích đất trồng lúa nước giảm 7,84 ha do chuyển đổi sang cây lâu năm và xây dựng công trình hạ tầng.
  • Đất trồng cây lâu năm và lâm nghiệp: Đất trồng cây lâu năm tăng 7,59 ha nhờ định hướng phát triển cây ăn quả và cây công nghiệp có giá trị cao; diện tích đất lâm nghiệp giảm nhẹ 15,76 ha do hoạt động khai thác rừng sản xuất và chu kỳ luân canh.
  • Đất mặt nước và đất xây dựng: Diện tích đất mặt nước chuyên dùng và nuôi trồng thủy sản tăng 13,92 ha; đất xây dựng sau khi chuẩn hóa hiện trạng phân loại ghi nhận sự sắp xếp lại với mức giảm 56,45 ha trên bề mặt giải đoán theo hiện trạng khai thác thực tế.

Thảo luận kết quả

Nguyên nhân căn bản của sự biến động bắt nguồn từ quá trình chuyển dịch cơ cấu kinh tế của huyện Thạch Thành từ thuần nông sang công nghiệp, dịch vụ và thương mại dọc tuyến Quốc lộ 45 dài 20,5 km. Việc tăng diện tích cây lâu năm và mở rộng mặt nước phản ánh rõ chủ trương tái cơ cấu ngành nông nghiệp, nâng cao giá trị gia tăng trên từng hecta canh tác.

Về phương thức trực quan hóa, toàn bộ luồng luân chuyển không gian được mô tả rõ nét qua bảng ma trận biến động hai chiều và biểu đồ hình cột phân bố diện tích trong môi trường GIS. Bảng ma trận phản ánh chi tiết số lượng hecta chuyển từ đất lúa sang đất phi nông nghiệp và từ đất cây hàng năm sang đất cây lâu năm. So với nhiều công trình nghiên cứu tại các tỉnh trung du miền núi phía Bắc đạt chỉ số Kappa bình quân từ 0,78 đến 0,85, kết quả của luận văn đạt Kappa 0,89 đã khẳng định tính ưu việt của thuật toán Maximum Likelihood. Giải pháp tích hợp này giúp tiết kiệm hơn 50% kinh phí và nhân lực so với phương pháp đo vẽ thủ công truyền thống.

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

Dựa trên kết quả phân tích biến động không gian, bốn nhóm giải pháp hành động cụ thể được khuyến nghị thực hiện:

  1. Xây dựng và chuẩn hóa cơ sở dữ liệu GIS địa chính động: Ủy ban nhân dân huyện Thạch Thành cùng Phòng Tài nguyên và Môi trường cần triển khai hệ thống quản trị dữ liệu đất đai số cho toàn bộ 28 xã, thị trấn trong giai đoạn 2019 - 2021, đảm bảo số hóa 100% hồ sơ địa chính và giảm tỷ lệ sai lệch thông tin diện tích xuống dưới 2%.
  2. Ứng dụng ảnh viễn thám độ phân giải cao giám sát tài nguyên rừng: Hạt Kiểm lâm huyện Thạch Thành phối hợp cùng các Ban Quản lý rừng ứng dụng ảnh vệ tinh Sentinel-2 có độ phân giải 10 m giai đoạn 2020 - 2025 nhằm giám sát liên tục toàn bộ diện tích đất lâm nghiệp, phát hiện kịp thời biến động thảm phủ rừng và kiểm soát diện tích suy giảm dưới mức 15 ha mỗi năm.
  3. Siết chặt quy hoạch chuyển đổi đất lúa và đất cây hàng năm: Ủy ban nhân dân các xã và thị trấn Kim Tân, Vân Du cần tăng cường kiểm tra việc thực thi quy hoạch sử dụng đất tầm nhìn 2021 - 2030, duy trì ổn định tối thiểu 5.500 ha đất trồng lúa chất lượng cao, hạn chế tình trạng chuyển đổi tự phát đất nông nghiệp sang đất phi nông nghiệp.
  4. Nâng cao năng lực công nghệ địa không gian cho cán bộ địa phương: Sở Tài nguyên và Môi trường tỉnh Thanh Hóa tổ chức định kỳ 2 khóa tập huấn mỗi năm về xử lý ảnh vệ tinh, vận hành phần mềm ArcGIS và thiết bị GPS cho 100% cán bộ địa chính cấp xã, nâng tỷ lệ giải quyết hồ sơ biến động đất đai đúng hạn lên trên 95%.

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

Luận văn mang lại giá trị tham khảo chuyên môn sâu sắc cho 4 nhóm đối tượng trọng tâm:

  1. Cơ quan quản lý Nhà nước về đất đai và tài nguyên môi trường: Lãnh đạo và cán bộ Phòng Tài nguyên và Môi trường huyện Thạch Thành, Chi nhánh Văn phòng Đăng ký đất đai có thể ứng dụng trực tiếp quy trình xử lý ảnh Landsat và ma trận GIS vào các kỳ kiểm kê đất đai 5 năm một lần, giúp tiết kiệm hơn 40% ngân sách và thời gian điều tra.
  2. Các tổ chức tư vấn quy hoạch và phát triển đô thị - nông thôn: Các viện nghiên cứu, kỹ sư quy hoạch có thể sử dụng bộ bản đồ biến động tỷ lệ 1:25.000 để mô hình hóa xu hướng phát triển không gian tại 28 xã, gắn kết hài hòa với quy hoạch vùng đệm rừng quốc gia Cúc Phương trong giai đoạn 2021 - 2030.
  3. Giảng viên, nhà nghiên cứu và sinh viên ngành Quản lý đất đai, GIS, Viễn thám: Đề tài cung cấp case study chuẩn mực từ quy trình tiền xử lý ảnh trên phần mềm ENVI, thu thập 100 mẫu kiểm định thực địa bằng GPS đến thuật toán tính hệ số Kappa và biên tập bản đồ chuyên đề.
  4. Các tổ chức bảo tồn tài nguyên và phát triển nông lâm nghiệp: Áp dụng phương pháp phân tích lớp phủ viễn thám để đánh giá tác động sinh thái, bảo vệ môi trường bền vững trên hơn 25.000 ha rừng và đồi núi tại khu vực phía Tây tỉnh Thanh Hóa.

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

  1. Vì sao luận văn lựa chọn tư liệu ảnh vệ tinh Landsat 8 để đánh giá biến động đất đai? Trả lời: Tư liệu Landsat 8 cung cấp ảnh đa phổ hoàn toàn miễn phí với độ phân giải không gian 30 m và chu kỳ chụp lặp 16 ngày, rất phù hợp với phạm vi quy mô huyện có diện tích 55.921,72 ha. Dữ liệu gồm 11 kênh phổ giúp phân tách rõ nét 6 loại hình sử dụng đất chính mà không phát sinh chi phí mua ảnh thương mại đắt tiền.

  2. Thuật toán Xác suất cực đại (Maximum Likelihood) có ưu thế gì trong nghiên cứu này? Trả lời: Thuật toán tính toán xác suất phân bố chuẩn Gauss của từng pixel thuộc vào các lớp phủ mẫu huấn luyện. Phương pháp này mang lại độ tin cậy vượt trội khi bản đồ năm 2014 đạt độ chính xác 86,00% (Kappa = 0,83) và năm 2017 đạt 91,00% (Kappa = 0,89), hạn chế tối đa hiện tượng phân loại nhầm giữa các đối tượng phổ gần nhau.

  3. Cỡ mẫu 100 điểm kiểm tra thực địa có bảo đảm tính đại diện khoa học cho toàn huyện? Trả lời: Cỡ mẫu 100 điểm được thiết kế theo phương pháp chọn mẫu ngẫu nhiên phân tầng trải rộng trên 28 xã, thị trấn. Việc kết hợp định vị tọa độ bằng GPS cầm tay trên 6 nhóm lớp phủ giúp đối soát trung thực giữa dữ liệu giải đoán và thực tế, đáp ứng đầy đủ tiêu chuẩn kiểm định của ngành trắc địa bản đồ.

  4. Nguyên nhân nào dẫn đến việc diện tích đất trồng cây hàng năm khác giảm 54,36 ha? Trả lời: Diện tích đất cây hàng năm giảm 54,36 ha chủ yếu do người dân chuyển đổi cơ cấu cây trồng sang đất trồng cây lâu năm có giá trị kinh tế cao hơn (tăng 7,59 ha) như cây ăn quả, mía tím, đồng thời mở rộng diện tích nuôi trồng thủy sản (tăng 13,92 ha) và chuyển dịch một phần sang phát triển hạ tầng cụm công nghiệp.

  5. Quy trình tích hợp GIS và viễn thám trong đề tài có khả năng chuyển giao cho địa phương khác không? Trả lời: Quy trình hoàn toàn có thể nhân rộng. Các bước chuẩn hóa từ xử lý ảnh số trên ENVI đến phân tích ma trận chồng lớp trên ArcGIS áp dụng hiệu quả cho tất cả các huyện miền núi có quy mô diện tích từ 40.000 ha đến 100.000 ha nhằm hiện đại hóa công tác giám sát tài nguyên đất đai.

Kết luận

  • Luận văn đã hoàn thành toàn diện việc xây dựng cơ sở dữ liệu số và đánh giá chính xác biến động sử dụng đất trên diện tích 55.921,72 ha của huyện Thạch Thành giai đoạn 2014 - 2017.
  • Khẳng định tính ưu việt của công nghệ viễn thám Landsat 8 và GIS với độ chính xác phân loại đạt 86,00% năm 2014 (Kappa = 0,83) và 91,00% năm 2017 (Kappa = 0,89).
  • Lượng hóa rõ nét xu hướng chuyển dịch với mức giảm 54,36 ha đất cây hàng năm khác, giảm 7,84 ha đất lúa nước, đồng thời tăng 7,59 ha đất cây lâu năm và 13,92 ha đất mặt nước.
  • Cung cấp luận cứ khoa học thực tiễn và đề xuất quy trình ứng dụng GIS đồng bộ cho 28 xã, thị trấn phục vụ lập quy hoạch sử dụng đất giai đoạn 2021 - 2030.
  • Khuyến nghị các cơ quan chuyên môn tiếp tục cập nhật tư liệu viễn thám Sentinel-2 chu kỳ 12 tháng một lần trong giai đoạn 2020 - 2025 để chủ động theo dõi diễn biến tài nguyên đất.

Hãy chủ động ứng dụng công nghệ viễn thám và GIS để hiện đại hóa công tác quản trị tài nguyên đất đai tại địa phương của bạn ngay hôm nay!