Dữ liệu NDVI được tính tốn riêng cho từng tổ hợp 8-ngày đã được ghép lại sử dụng chức năng Layer Stacking của ENVI để tạo thành ảnh đa kênh thời gian (46 kênh). Sản phẩm chỉ số thực vật được sử dụng để tạo bản đồ và theo dõi biến động lớp phủ theo chỉ số thực vật. Đồ thị chuỗi thời gian theo các chỉ số thực vật tại các khu vực có lớp phủ khác nhau được xây dựng với chức năng Z-profile của ENVI.
3.2.3. Tạo ảnh tổ hợp NDVI theo tháng
Mặc dù bộ dữ liệu MODIS đã được xử lý tới mức 3 tức là đã làm giảm bớt hoặc loại bỏ được mây nhưng vẫn còn hiện tượng mức độ mây che phủ tương đối lớn. Điều này gây khó khăn cho việc theo dõi biến động cũng như phân loại lớp phủ. Do vậy, học viên tiếp tục tiến hành thêm bước lọc mây và tạo ảnh tổ hợp theo tháng.
Ngưỡng giá trị phản xạ đã được lựa chọn dựa trên sự tham khảo các tài liệu trong nước và quốc tế và dựa trên sự phân tích thực tế bộ dữ liệu ảnh MODIS thu thập được. Theo đó, giá trị phản xạ > 0,2 đối với kênh 1 (Red) và > 0,5 đối với kênh
2 (Near infared) được coi là mây và thuật toán lọc mây đã được tiến hành sử dụng cơng cụ tính tốn Band Math của phần mềm ENVI. Việc lọc mây được thực hiện trên ảnh NDVI. Các pixel được coi là mây sẽ được gán giá trị đặc biệt (No data) là - 999 để phân biệt với các đối tượng khác, đồng thời phục vụ cho việc tổ hợp ảnh 16 - ngày và tổ hợp ảnh theo tháng sau này. Các pixel được giải đốn khơng phải là mây sẽ được trả về giá trị NDVI thực tế của nó.
Thuật tốn dùng để lọc và gán giá trị - 999 cho mây sử dụng công cụ Band Math trên phần mềm ENVI có dạng như sau:
((b1 GT 0.2) AND (b2 GT 0.5))*(-999) + b3*((b1 LE 0.2) OR (b2 LE 0.5))
Trong đó: b1: là giá trị của pixel ở band 1 trên ảnh gốc 8-ngày; b2: là giá trị của pixel ở band 2 trên ảnh gốc 8-ngày;
b3: là giá trị tính tốn NDVI trên ảnh tổ hợp 8-ngày của pixel đó; Thuật tốn trên được hiểu là nếu giá trị của pixel trên kênh 1 lớn hơn 0,2 và trên kênh 2 lớn hơn 0,5 thì gán giá trị cho pixel đó là - 999, nếu khơng thì trả về giá trị NDVI thực của pixel đó.