Tổng quan nghiên cứu

Huyện Thạch Thành là một huyện miền núi nằm ở phía Tây Bắc tỉnh Thanh Hóa, sở hữu tổng diện tích tự nhiên 55.921,72 ha và được phân chia thành 28 đơn vị hành chính cấp xã, bao gồm 26 xã và 02 thị trấn. Trong giai đoạn 2014 - 2017, khu vực này chứng kiến tốc độ đô thị hóa nhanh chóng và sự chuyển dịch mạnh mẽ của cơ cấu kinh tế từ nông nghiệp sang công nghiệp - dịch vụ. Sự xuất hiện của các trục giao thông huyết mạch như Quốc lộ 45 dài 20,5 km, đường Hồ Chí Minh dài 13,3 km cùng việc mở rộng cụm công nghiệp Vân Du và Thạch Quảng đã tạo áp lực lớn lên tài nguyên đất đai, dẫn đến sự biến đổi liên tục của các loại hình lớp phủ bề mặt.

Vấn đề cốt lõi đặt ra cho công tác quản lý nhà nước tại địa phương là các phương pháp thống kê, kiểm kê đất đai truyền thống thường tốn nhiều thời gian, kinh phí và thiếu tính trực quan về mặt không gian. Đề tài tập trung giải quyết bài toán ứng dụng công nghệ Hệ thống Thông tin Địa lý (GIS) và Viễn thám nhằm xây dựng bản đồ hiện trạng và định lượng chính xác biến động sử dụng đất trên địa bàn toàn huyện trong giai đoạn 2014 - 2017.

Mục tiêu cụ thể của công trình là phân loại và đánh giá sự dịch chuyển giữa 46.200,88 ha đất nông nghiệp, 8.221,35 ha đất phi nông nghiệp và 1.509 ha đất chưa sử dụng. Về mặt ý nghĩa, nghiên cứu cung cấp một giải pháp công nghệ có độ tin cậy cao với độ chính xác phân loại ảnh đạt từ 86,00% đến 91,00%, giúp cắt giảm khoảng 40% thời gian thực hiện điều tra ngoại nghiệp. Đây là cơ sở khoa học và thực tiễn vững chắc phục vụ việc điều chỉnh quy hoạch, kế hoạch sử dụng đất kỳ cuối 2016 - 2020 huyện Thạch Thành theo đúng định hướng của Luật Đất đai năm 2013.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Luận văn xây dựng trên nền tảng của hai khung lý thuyết khoa học hiện đại: Lý thuyết bức xạ điện từ trong viễn thám và Lý thuyết mô hình hóa không gian địa lý của GIS. Về mặt cấu trúc phân tích, đề tài áp dụng mô hình phân loại ảnh đa phổ có kiểm định kết hợp với mô hình so sánh biến động sau phân loại thông qua ma trận chuyển dịch hai chiều.

Hệ thống lý thuyết tích hợp 5 khái niệm chuyên ngành then chốt:

  • Viễn thám (Remote Sensing): Khoa học thu nhận và phân tích thông tin về các đối tượng trên bề mặt Trái Đất từ khoảng cách xa thông qua việc đo lường năng lượng bức xạ và phản xạ sóng điện từ.
  • Hệ thống thông tin địa lý (GIS): Tổ hợp công cụ phần cứng, phần mềm, dữ liệu không gian và con người nhằm thu thập, lưu trữ, truy vấn, phân tích và hiển thị các thông tin gắn liền với tọa độ địa lý.
  • Biến động sử dụng đất (LUCC): Quá trình biến đổi trạng thái tự nhiên và công năng của các loại đất dưới tác động tương hỗ giữa các quá trình tự nhiên và hoạt động kinh tế - xã hội của con người.
  • Chỉ số thực vật khác biệt chuẩn hóa (NDVI): Chỉ số phổ biểu thị mối tương quan giữa kênh cận hồng ngoại và kênh đỏ nhằm định lượng mật độ và sức khỏe của thảm thực vật.
  • Thuật toán xác suất cực đại (Maximum Likelihood Classification): Phương pháp phân loại thống kê có kiểm định dựa trên hàm mật độ xác suất phân bố chuẩn Gauss để gán từng pixel vào lớp phủ có xác suất thuộc về cao nhất.

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

Nghiên cứu kết hợp đồng bộ giữa nguồn dữ liệu thứ cấp và sơ cấp để đảm bảo tính chuẩn xác của kết quả:

  • Nguồn dữ liệu thứ cấp: Sử dụng hai cảnh ảnh vệ tinh Landsat 8 mã hiệu 127-046 chụp vào ngày 19/9/2014 và ngày 08/9/2017 với độ phân giải không gian 30 m; hồ sơ kiểm kê đất đai năm 2014, thống kê đất đai năm 2017 và bản đồ địa chính của 28 xã, thị trấn thuộc huyện Thạch Thành.
  • Nguồn dữ liệu sơ cấp và cỡ mẫu: Tác giả thiết lập cỡ mẫu gồm 100 điểm khảo sát thực địa trên toàn huyện bằng thiết bị định vị vệ tinh GPS cầm tay.
  • Phương pháp chọn mẫu: Áp dụng phương pháp chọn mẫu phân tầng ngẫu nhiên theo từng loại hình sử dụng đất chính như đất lúa, đất cây hàng năm khác, đất cây lâu năm, đất rừng, đất nuôi trồng thủy sản và đất xây dựng. Phương pháp này đảm bảo tính đại diện cao trên một vùng địa hình phức tạp có độ cao dao động mạnh từ 15 m đến 825 m.

Lý do lựa chọn phương pháp phân tích: Tác giả sử dụng phần mềm ENVI để tiền xử lý ảnh gồm nắn chỉnh hình học về hệ quy chiếu WGS84, tăng cường chất lượng ảnh và phân loại bằng thuật toán Maximum Likelihood. Thuật toán này được lựa chọn vì có khả năng phân tách phổ tối ưu giữa các lớp phủ thực vật phức tạp tại vùng miền núi. Sau đó, phần mềm ArcGIS được ứng dụng để thực hiện phép chồng xếp không gian giữa hai bản đồ số 2014 và 2017 nhằm trích xuất ma trận biến động. Toàn bộ timeline nghiên cứu được tiến hành bài bản từ khâu thu thập ảnh vệ tinh, xử lý số liệu, kiểm chứng 100 điểm thực địa đến khi hoàn thiện tập bản đồ chuyên đề tỷ lệ 1:25.000 trong giai đoạn 2017 - 2018.

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

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

Kết quả xử lý ảnh viễn thám đa thời gian và phân tích GIS đã xác định rõ sự chuyển dịch cơ cấu của các loại hình sử dụng đất tại huyện Thạch Thành trong giai đoạn 2014 - 2017:

  • Thứ nhất, diện tích đất trồng cây nông nghiệp ngắn ngày có xu hướng giảm. Trong đó, đất trồng lúa nước giảm 7,84 ha (tương đương mức giảm khoảng 0,11%) và đất trồng cây hàng năm khác ghi nhận mức giảm đáng kể 54,36 ha so với thời điểm năm 2014.
  • Thứ hai, đất trồng cây lâu năm có xu hướng mở rộng với mức tăng 7,59 ha, phản ánh sự chuyển đổi diện tích đất đồi thấp kém hiệu quả sang trồng các loại cây ăn quả và cây công nghiệp. Ngược lại, diện tích đất lâm nghiệp giảm 15,76 ha do các hoạt động khai thác rừng trồng và chuyển mục đích sử dụng đất.
  • Thứ ba, diện tích đất mặt nước chuyên dùng và nuôi trồng thủy sản tăng thêm 13,92 ha, tập trung chủ yếu tại các xã ven lưu vực sông Bưởi nhằm mở rộng các hồ chứa thủy lợi và cơ sở hạ tầng nuôi cá nước ngọt.
  • Thứ tư, đất xây dựng và phi nông nghiệp diễn ra sự biến động tái cơ cấu với 56,45 ha đất chuyển dịch công năng để xây dựng hệ thống hạ tầng giao thông và các công trình phúc lợi công cộng.

Thảo luận kết quả

Nguyên nhân chính dẫn đến các biến động nêu trên bắt nguồn từ tốc độ chuyển dịch cơ cấu kinh tế địa phương. Việc triển khai các dự án nâng cấp Quốc lộ 45, tỉnh lộ 516, tỉnh lộ 522 cùng việc phát triển các cơ sở chế biến nông - lâm sản đã thúc đẩy quá trình chuyển đổi đất nông nghiệp sang đất phi nông nghiệp. Khi so sánh với các nghiên cứu biến động đất đai tại các địa bàn trung du miền núi lân cận, xu hướng giảm dần đất cây hàng năm và tăng diện tích cây lâu năm, đất hạ tầng tại Thạch Thành hoàn toàn phù hợp với quy luật phát triển kinh tế nông thôn mới.

Dữ liệu nghiên cứu được biểu diễn trực quan thông qua bảng ma trận biến động hai chiều giữa các cặp loại đất và biểu đồ hình cột so sánh biến động cơ cấu diện tích. Kết quả đánh giá độ chính xác phân loại bằng 100 điểm kiểm chứng thực tế cho thấy:

  • Bản đồ hiện trạng năm 2014 đạt độ chính xác tổng thể 86,00% với hệ số Kappa k = 0,83.
  • Bản đồ hiện trạng năm 2017 đạt độ chính xác tổng thể 91,00% với hệ số Kappa k = 0,89.

Độ chính xác của ảnh năm 2017 cao hơn năm 2014 do thời điểm chụp ảnh gần sát với thời gian khảo sát thực địa. Khi đối soát diện tích giải đoán với số liệu thống kê đất đai chính thức của huyện Thạch Thành, tỷ lệ chênh lệch đều nằm trong giới hạn cho phép dưới 5%, khẳng định phương pháp tích hợp viễn thám và GIS hoàn toàn đủ độ tin cậy để thay thế các phương pháp trắc địa truyền thống trong giám sát diện rộng.

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

Nhằm tối ưu hóa công tác quản trị tài nguyên đất đai và bảo đảm sử dụng đất bền vững, luận văn đưa ra 4 nhóm giải pháp trọng tâm:

  1. Xây dựng và hoàn thiện cơ sở dữ liệu địa chính số hóa đa mục tiêu bằng phần mềm GIS cho toàn bộ 28 xã, thị trấn, đặt mục tiêu đạt 100% hồ sơ đất đai được quản lý trên nền tảng không gian số trước năm 2025 dưới sự điều hành của Phòng Tài nguyên và Môi trường huyện Thạch Thành.
  2. Thiết lập quy chế kiểm soát nghiêm ngặt việc chuyển mục đích sử dụng đất trồng lúa nước và đất lâm nghiệp sang đất phi nông nghiệp, khống chế tỷ lệ suy giảm diện tích đất lúa dưới 0,5%/năm trong giai đoạn 2020 - 2030, do Ủy ban nhân dân huyện Thạch Thành phối hợp với Sở Tài nguyên và Môi trường tỉnh Thanh Hóa thực hiện.
  3. Ứng dụng công nghệ viễn thám giám sát định kỳ 6 tháng một lần đối với 1.509 ha đất chưa sử dụng và các khu vực có nguy cơ sạt lở ven sông Bưởi nhằm phát hiện kịp thời các hành vi lấn chiếm đất công, giao nhiệm vụ cho Chi nhánh Văn phòng Đăng ký đất đai triển khai từ năm 2021.
  4. Nâng cao tỷ lệ cấp mới và cấp đổi Giấy chứng nhận quyền sử dụng đất đạt trên 98% cho các hộ sản xuất nông - lâm nghiệp, đồng thời tổ chức 12 đợt tập huấn chuyên sâu về kỹ năng phân tích ảnh vệ tinh và vận hành ArcGIS cho đội ngũ cán bộ địa chính cấp xã đến hết năm 2023.

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

Công trình nghiên cứu mang lại giá trị ứng dụng thiết thực cho 4 nhóm đối tượng chính:

  1. Cơ quan quản lý nhà nước về tài nguyên và môi trường cấp tỉnh và huyện (Sở Tài nguyên và Môi trường Thanh Hóa, UBND huyện Thạch Thành): Ứng dụng quy trình tích hợp viễn thám để cập nhật nhanh bản đồ hiện trạng cho 55.921,72 ha diện tích tự nhiên, phục vụ đắc lực cho công tác lập quy hoạch sử dụng đất định kỳ 5 năm và 10 năm.
  2. Cán bộ địa chính và chuyên viên quy hoạch đô thị cấp cơ sở: Sử dụng các phương pháp phân tích chồng lớp không gian của GIS để theo dõi diễn biến chuyển đổi 56,45 ha đất xây dựng và quản lý chặt chẽ quỹ đất công cộng tại 28 xã, thị trấn.
  3. Giảng viên, học viên cao học và sinh viên các ngành Quản lý đất đai, Hệ thống thông tin địa lý, Bản đồ học và Khoa học môi trường: Khai thác luận văn như một tài liệu tham khảo chuẩn mực về phương pháp giải đoán ảnh Landsat 8, kỹ thuật kiểm tra 100 điểm GPS ngoại nghiệp và thuật toán tính hệ số Kappa đạt trên 0,80.
  4. Các đơn vị tư vấn đo đạc bản đồ và tổ chức bảo vệ môi trường: Tiếp cận phương pháp tích hợp ảnh vệ tinh đa thời gian để đánh giá diễn biến suy giảm rừng, phân tích tác động của xói mòn và đề xuất các giải pháp bảo vệ lưu vực sông Bưởi một cách bền vững.

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

  1. Thuật toán Maximum Likelihood mang lại lợi thế gì khi phân loại ảnh viễn thám huyện Thạch Thành? Thuật toán phân loại xác suất cực đại tính toán hàm phân bố xác suất Gauss cho từng pixel dựa trên mẫu tập huấn, giúp phân biệt rõ nét các lớp thực vật có phổ gần nhau. Nhờ đó, bản đồ năm 2017 đạt độ chính xác tới 91,00% và chỉ số Kappa k = 0,89 trên nền địa hình chia cắt phức tạp.

  2. Vì sao diện tích đất trồng cây hàng năm khác tại Thạch Thành lại giảm 54,36 ha trong giai đoạn 2014 - 2017? Sự suy giảm 54,36 ha đất cây hàng năm xuất phát từ quá trình chuyển dịch cơ cấu nông nghiệp sang trồng cây ăn quả lâu năm có hiệu quả kinh tế cao hơn (tăng 7,59 ha), đồng thời một phần diện tích được thu hồi để phát triển giao thông và mở rộng hạ tầng khu công nghiệp Vân Du.

  3. Kỹ thuật đánh giá độ chính xác của bản đồ giải đoán trong luận văn được thực hiện như thế nào? Tác giả đã tiến hành định vị và kiểm tra 100 điểm mẫu ngoài thực địa bằng thiết bị GPS cầm tay trên khắp 28 đơn vị hành chính. Kết quả đối soát ma trận sai số xác nhận độ chính xác bản đồ năm 2014 đạt 86,00% (k = 0,83) và bản đồ năm 2017 đạt 91,00% (k = 0,89), đáp ứng đầy đủ tiêu chuẩn khoa học.

  4. Ảnh vệ tinh Landsat 8 có đáp ứng tốt yêu cầu thành lập bản đồ địa hình cấp huyện không? Với độ phân giải không gian 30 m cùng 11 kênh phổ đa dạng, ảnh Landsat 8 hoàn toàn đáp ứng độ chính xác để thành lập bản đồ hiện trạng và biến động tỷ lệ 1:25.000 cho địa bàn rộng 55.921,72 ha, giúp giảm trên 40% chi phí và nhân lực so với phương pháp đo trắc địa truyền thống.

  5. Việc ứng dụng công nghệ GIS mang lại giá trị gia tăng nào cho dữ liệu viễn thám? Sau khi viễn thám phân loại thành công các lớp phủ, GIS cung cấp các công cụ phân tích không gian mạnh mẽ như chồng lớp bản đồ và truy vấn thuộc tính để tạo lập ma trận biến động hai chiều. Nhờ đó, người quản lý xác định được chính xác vị trí và diện tích chuyển đổi của từng loại đất như 13,92 ha mặt nước hay 15,76 ha đất lâm nghiệp.

Kết luận

  • Luận văn đã xây dựng thành công quy trình công nghệ tích hợp GIS và viễn thám Landsat 8 để thành lập bản đồ hiện trạng và bản đồ biến động sử dụng đất huyện Thạch Thành giai đoạn 2014 - 2017 ở tỷ lệ 1:25.000.
  • Kết quả kiểm chứng thực địa trên 100 điểm mẫu đạt độ chính xác toàn thể từ 86,00% đến 91,00%, hệ số Kappa đạt 0,83 - 0,89 khẳng định tính ưu việt và độ tin cậy của thuật toán Maximum Likelihood.
  • Nghiên cứu đã định lượng rõ nét xu hướng giảm 54,36 ha đất cây hàng năm, giảm 7,84 ha đất lúa nước, tăng 7,59 ha đất cây lâu năm và tăng 13,92 ha đất mặt nước do tác động của chuyển dịch cơ cấu kinh tế.
  • Đóng góp giải pháp khoa học thực tiễn giúp các cơ quan quản lý rút ngắn thời gian cập nhật dữ liệu, nâng cao hiệu lực thi hành Luật Đất đai năm 2013 và phục vụ công tác quy hoạch giai đoạn 2020 - 2025.
  • Các cơ quan quản lý đất đai và độc giả quan tâm hãy ứng dụng ngay mô hình tích hợp GIS và ảnh vệ tinh đa thời gian này để tối ưu hóa công tác giám sát tài nguyên đất đai tại địa phương.