Lớp phân loại Dữ liệu kiểm chứng Tổng hàng (điểm) ĐCX sản xuất (%) (1) (2) (3) (4) (5) (6) Đất trồng lúa (1) 26 1 0 0 2 1 30 86,67 Đất trồng cây HNK (2) 2 27 1 0 0 0 30 90,00 Đất rừng-cây lâu năm (3) 1 0 28 0 1 0 30 93,33 Đất xây dựng (4) 0 0 0 27 0 3 30 90,00 Đất mặt nước (5) 0 0 1 2 27 0 30 96,67 Đất chưa sử dụng (6) 1 2 1 0 0 26 30 86,67 Tổng cột (điểm) 30 30 31 29 30 30 180 ĐCX sử dụng (%) 86,67 90,00 90,32 93,10 90,00 86,67 ĐCX tổng thể 89,44% Chỉ số Kappa: 0,88
Trong bảng 4.8, đối với đất trồng lúa, trong 30 điểm kiểm chứng thì cũng có 2 điểm bị phân loại nhầm sang lớp mặt nước, do tại thời điểm chụp ảnh vệ tinh, đất trồng lúa nước có độ phản xạ phổ gần giống với đất mặt nước. Đối với đất chưa sử dụng, trong 30 điểm kiểm chứng thì cũng có 3 điểm bị phân loại nhầm sang lớp đất xây dựng. Đối với đất xây dựng có 4/30 bị nhầm lẫn sang đất chưa sử dụng do độ phản xạ phổ của hai lớp này gần với nhau.
Độ chính xác tổng thể ảnh phân loại năm 2018 là 89,44% tương ứng với chỉ số Kappa bằng 0,88, đạt yêu cầu.
Trước khi chuyển kết quả phân loại sang dạng bản đồ chúng ta cần phải thực hiện một số thao tác cần thiết như khái quát hóa các lớp thông tin, véc tơ hóa ảnh sau khi phân loại …
- Ghép nhóm các đối tượng: là việc ghép các lớp phân loại có cùng tính chất giống nhau thành một nhóm. Với một loại đối tượng có thể có nhiều giá trị độ xám có độ chênh lệch nhau lớn, còn khi chọn vùng mẫu chỉ chọn những vùng đặc trưng dẫn đến việc cùng một đối tượng nhưng ở các lớp khác nhau vì vậy phải gộp chúng lại về cùng một đối tượng để thể hiện trên bản đồ.
- Phân tích theo đa số và thiểu số: Đây là phương pháp phân tích để gộp các pixel lẻ tẻ hoặc phân loại lẫn trong các lớp vào chính lớp chứa nó. Phép lọc trong ENVI là phép lọc Majority/Minority Analysis. Sau khi nhập kích thước của sổ lọc thì giá trị của pixel trung tâm sẽ được thay thế bởi giá trị của pixel chiếm đa số (cách phân tích theo đa số) hoặc bằng giá trị của pixel chiếm thiểu số (cách phân tích theo thiểu số). Kết quả sẽ cho ta một ảnh mới. Kết quả của phép lọc không nâng cao được độ chính xác của kết quả phân loại mà chỉ cải thiện về mặt hình thức.
- Chuyển kết quả phân loại sang dạng véc tơ: Để biên tập và thành lập bản đồ sử dụng đất, từ kết quả phân loại ảnh ta phải chuuyển kết quả phân loại sang dạng véctơ. Hai ảnh vệ tinh năm 2010 và 2018 sau khi phân loại được chuyển sang dạng véctơ và lưu ở định dạng file 2010.evf và 2018.evf. Sau đó chuyển sang dạng shape file để nhập vào phần mềm Arc GIS.
4.3.3. Thành lập bản đồ biến động sử dụng đất
4.3.3.1. Đối soát thực địa
Thực chất kết quả phân loại ảnh thu được hiện trạng lớp phủ tại thời điểm phân loại, do vậy phải tiến hành đối soát thực địa để thành lập bản đồ sử dụng đất.
Tiến hành in kết quả phân loại và đối soát với thực địa để chỉnh lý kết quả phân loại đúng với hiện trạng sử dụng đất. Việc đối soát thực địa thực hiện với ảnh năm 2018, ảnh năm 2010 đối soát với bản đồ hiện trạng sử dụng đất.
Từ kết quả chỉnh lý tiến hành biên tập, thành lập bản đồ sử dụng đất.
4.3.3.2. Biên tập bản đồ sử dụng đất
Ảnh phân loại sau khi được chuyển sang dạng véctơ được nhập vào phần mềm ArcGis để biên tập thành lập bản đồ sử dụng đất. Tiến hành đổi màu, biên tập các yếu tố để thành lập bản đồ sử dụng đất.
Sau khi ảnh được phân loại, đánh giá độ chính xác vẫn cần phải kiểm tra, đối soát các lớp bản đồ với thực địa, sau đó chỉnh lý bản đồ theo đúng hiện trạng sử dụng đất. Kết quả thu được hai bản đồ sử dụng đất khu vực nghiên cứu năm 2010 và 2018 (thể hiện ở phụ lục 3 và phụ lục 4).
Hình 4.8 và 4.9 là sơ đồ sử dụng đất khu vực nghiên cứu năm 2010 và năm 2018 được thu nhỏ.
Để thống kê diện tích từng loại đất trên bản đồ ta sử dụng công cụ tính diện tích của các đối tượng dạng vùng trong phần mềm ArcGis.
Để đánh giá độ chính xác của bản đồ sử dụng đất giải đoán từ ảnh vệ tinh ta so sánh diện tích trên bản đồ giải đoán so với số liệu thống kê của khu vực nghiên cứu.
Kết quả thống kê diện tích giải đoán ảnh năm 2010, năm 2018 và so sánh chênh lệch giữa diện tích thống kê và diện tích giải đoán các loại đất trên bản đồ hiện trạng sử dụng đất năm 2010, năm 2018 được thể hiện trong bảng 4.9 và bảng 4.10.