27 15,37499 107,81626 Loại khác Loại khác ==
5.2.2.5. Xây dựng ảnh phân loại NDVI và TRR
+ Xây dựng ảnh phân loại NDVI
Để phân loại lớp phủ thực vật theo cấu trúc hình thái, chúng tôi đã sử dụng loại chỉ số này để phân loại trạng thái rừng trên t− liệu ảnh vệ tinh LANDSAT TM đ−ợc chụp qua các thời kỳ 1989 – 2001.
87
ở đây, chỉ số thực vật sử dụng trong máy tính đ−ợc tính theo công thức:
NDVI = [(NIR – Red)/(NIR + Red)]*127 +127. (5.5) Kết quả, thu đ−ợc ảnh chỉ số thực vật ở cả hai thời điểm 1989- 2001 và cho phép hiện trên màn hình nh− một kênh ảnh đen trắng.
Năm 1989 Năm 2001
Hình 5.13: ảnh chỉ số thực vật NDVI năm 1989, 2001
Muốn phân loại các đối t−ợng và thể hiện chúng bằng màu sắc, ta phải phân ng−ỡng NDVI cho các đối t−ợng đó. Mỗi ng−ỡng đặt cho nó một màu thích hợp và đ−ợc gọi là việc hiện màu giả.
Để chọn ng−ỡng thích hợp chúng tôi đã sử dụng ảnh phân loại có giám định năm 2001 chồng xếp với ảnh NDVI-2001. Trên mỗi trạng thái lớp phủ và lấy ngẫu nhiên 15 vị trí, xác định giá trị NDVI cho từng điểm, lấy giá trị trung bình làm chỉ số NDVI cho từng trạng thái. Dựa vào chỉ số trung bình
88
NDVI của các lớp phủ chúng tôi chia ng−ỡng cho các đối t−ợng t−ơng ứng (Bảng 5.7) sẽ tạo đ−ợc ảnh phân loại chỉ số thực vật của năm 2001. Ng−ỡng NDVI sẽ đ−ợc sử dụng để phân loại cho ảnh 1989
+ Xây dựng ảnh phân loại chỉ số tổng năng l−ợng phản xạ TRRI
Nh− đã giới thiệu ở mục (4.2.2.2), chỉ số tổng năng l−ợng phản xạ TRRI đ−ợc tính theo công thức (4.2). Vì vậy, để xây dựng ảnh chỉ số tổng năng l−ợng phản xạ TRRI ta phải nhập công thức trên vào máy cho cả 2 ảnh 1989_2001 , kết quả nhận đ−ợc 2 ảnh t−ơng ứng và cũng cho phép hiện bằng hình ảnh nh− một kênh ảnh đen trắng.
Năm 1989 Năm 2001
Hình 5.14:ảnh chỉ số tổng năng l−ợng phản xạ TRRI năm 1989,
2001
Tiến hành chồng xếp ảnh TRRI 2001 với ảnh phân loại có giám định năm 2001 để tìm ng−ỡng thích hợp cho các loại thảm rừng (Bảng 5.7).
89
Từ các ng−ỡng này có thể xây dựng ảnh phân loại chỉ số tổng năng l−ợng phản xạ TRRI cho thời điểm 1989.
Bảng 5.7: Ng−ỡng phân loại các loại thảm rừng
Ng−ỡng NDVI TRRI Loại thảm rừng
1 0-105 0-11 Các loại đất khác(Đất trống, thổ c−, n−ơng rẫy…)
2 105-137 11-19 Cỏ
3 137-190 19-27 Cây bụi và cây gỗ rải rác
4 190-202 27-35 Tre nứa
5 202-221 35-46 Rừng th−a 6 221-234 46-58 Rừng trung bình
7 234-255 58-87 Rừng kín
Để có ảnh phân loại năm 1989, chúng tôi xây dựng ảnh phân loại năm 1989 từ giá trị ng−ỡng phân loại tổ hợp NDVI và TRRI năm 2001.
Chúng tôi đã sử dụng ch−ơng trình ERDAT 8.5 và phần mềm ARCVIEW 3.2a để phân loại tự động, tính toán các chỉ số giá trị trên ảnh NDVI và TRRI. Ch−ơng trình xây dựng dựa trên nguyên lý là: dựa trên tệp ảnh đầu vào NDVI và TRRI ( trên ảnh 2001), ch−ơng trình sẽ tính cặp giá trị NDVI và TRRI t−ơng ứng trên ảnh 1989 sau đó so sánh với bảng ng−ỡng phân loại các đối t−ợng để xác định điểm ảnh thuộc vào đối t−ợng nào trong bảng phân loại.
Kết quả phân loại tại các thời điểm 1989 và 2001 đ−ợc thể hiện trên hình 5.15
90
Năm 1989 Năm 2001
Hình 5.15: ảnh phân loại năm 1989 và 2001
Vì t− liệu năm 1989 không có mẫu giám định, nên chúng tôi đã sử dụng ph−ơng pháp phân loại tổ hợp dựa vào ảnh phân loại chỉ số thực vật NDVI và ảnh phân loại chỉ số tổng năng l−ợng phản xạ TRRI. Nh− đã trình bày ở trên các chỉ số NDVI và TRRI đ−ợc tính từ giá trị các kênh ảnh năm 1989 và 2001 mà các kênh ảnh năm 1989 và 2001 đã đ−ợc chuẩn hoá. Trên cơ sở mức độ sai lệch cho phép của ảnh phân loại năm 2001 so với thực tế đã kiểm chứng thì ta cũng có kết quả ảnh phân loại năm 1989 với độ sai lệch cho phép.