Luận văn thạc sĩ HUS nghiên cứu ứng dụng công nghệ viễn thám và GIS đánh giá tác động của nhiệt độ, độ ẩm đến lớp phủ thực vật thông qua chỉ số thực vật (NDVI) khu vực tây nguyên

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - PHẠM VĂN MẠNH NGHIÊN CỨU ỨNG DỤNG CÔNG NGHỆ VIỄN THÁM VÀ GIS ĐÁNH GIÁ TÁC ĐỘNG CỦA NHIỆT ĐỘ, ĐỘ ẨM ĐẾN LỚP PHỦ THỰC VẬT THÔNG QUA CHỈ SỐ THỰC VẬT (NDVI) KHU VỰC TÂY NGUYÊN LUẬN VĂN THẠC SĨ KHOA HỌC Hà Nội – 2013 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - PHẠM VĂN MẠNH NGHIÊN CỨU ỨNG DỤNG CÔNG NGHỆ VIỄN THÁM VÀ GIS ĐÁNH GIÁ TÁC ĐỘNG CỦA NHIỆT ĐỘ, ĐỘ ẨM ĐẾN LỚP PHỦ THỰC VẬT THÔNG QUA CHỈ SỐ THỰC VẬT (NDVI) KHU VỰC TÂY NGUYÊN Chuyên ngành: Bản đồ viễn thám Hệ thông tin địa lý Mã số: 60440214 LUẬN VĂN THẠC SĨ KHOA HỌC NGƯỜI HƯỚNG DẪN KHOA HỌC: PGS.TS PHẠM VĂN CỰ Hà Nội – 2013 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com LỜI CẢM ƠN Trước tiên, xin gửi lời cảm ơn tới thầy hướng dẫn luận văn tơi, phó giáo sư - tiến sĩ Phạm Văn Cự, thầy tạo điều kiện, tận tình hướng dẫn giúp đỡ tơi hồn thành tốt luận văn Sự hiểu biết sâu sắc khoa học, kinh nghiệm thầy tiền đề giúp đạt thành tựu kinh nghiệm quý báu Tôi xin gửi lời cảm ơn đến đồng nghiệp, cấp trên, thủ trưởng quan Cục Viễn thám quốc gia, Đài Viễn thám Trung ương tận tình giúp đỡ tạo nhiều điều kiện thuận lợi cho tơi q trình nghiên cứu Cuối tơi xin cảm ơn bạn bè gia đình bên tôi, cổ vũ động viên lúc khó khăn để vượt qua hồn thành tốt luận văn Tơi xin chân thành cảm ơn! LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com MỤC LỤC MỞ ĐẦU CHƯƠNG TỔNG QUAN VỀ PHƯƠNG PHÁP SỬ DỤNG DỮ LIỆU MODIS TRONG ĐÁNH GIÁ MỐI TƯƠNG QUAN GIỮA NHIỆT ĐỘ BỀ MẶT, ĐỘ ẨM KHƠNG KHÍ VÀ LỚP PHỦ THỰC VẬT 1.1 Tình hình nghiên cứu giới khu vực 1.2 Tình hình nghiên cứu nước 10 CHƯƠNG SỬ DỤNG DỮ LIỆU ẢNH MODIS TRONG ĐÁNH GIÁ NHIỆT ĐỘ BỀ MẶT, ĐỘ ẨM KHƠNG KHÍ VÀ CHỈ SỐ THỰC VẬT (NDVI) CỦA LỚP PHỦ THỰC VẬT 17 2.1 Đặc tính phổ ảnh MODIS 17 2.1.1 Ảnh vệ tinh MODIS 17 2.1.2 Thiết bị thu nhận MODIS 19 2.2 Phương pháp tính nhiệt độ bề mặt (LST) từ liệu MODIS 22 2.2.1 Cơ sở khoa học chiết xuất nhiệt độ bề mặt từ tư liệu ảnh viễn thám 22 2.2.2 Phương pháp chiết xuất nhiệt độ bề mặt từ tư liệu viễn thám 24 2.2.3 Phương pháp chiết xuất nhiệt độ bề mặt từ ảnh MODIS 27 2.3 Phương pháp tính độ ẩm khơng khí từ liệu MODIS 28 2.3.1 Các loại độ ẩm khơng khí 28 2.3.2 Phương pháp tính tốn độ ẩm khơng khí từ tư liệu viễn thám 30 2.4 Chiết xuất NDVI Lớp phủ thực vật từ liệu MODIS 31 2.4.1 Khái niệm số thực vật 31 2.4.2 Cơ sở vật lý số thực vật 32 2.4.3 Cơng thức tính số NDVI 33 2.4.4 Phương pháp xử lý liệu NDVI 34 ii LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com CHƯƠNG TÍCH HỢP VIỄN THÁM VÀ GIS PHÂN TÍCH MỐI TƯƠNG QUAN GIỮA NHIỆT ĐỘ BỀ MẶT, ĐỘ ẨM KHƠNG KHÍ VÀ LỚP PHỦ THỰC VẬT NĂM 2012 .36 3.1 Vùng nghiên cứu 36 3.1.1 Vị trí địa lý 36 3.1.2 Khí hậu 37 3.1.3 Dân cư 38 3.1.4 Kinh tế, tài nguyên, xã hội môi trường 39 3.2 Các liệu sử dụng 40 3.2.1 Dữ liệu ảnh MODIS 40 3.2.2 Dữ liệu khác 42 3.3 Xây dựng tổ hợp ảnh nhiệt độ trung bình theo năm 42 3.3.1 Tạo ảnh nhiệt độ bề mặt từ ảnh MODIS – MOD11A2 42 3.3.2 Tạo ảnh tổ hợp theo tháng, ảnh tổ hợp theo mùa ảnh tổ hợp theo năm 43 3.4 Xây dựng tổ hợp ảnh độ ẩm trung bình theo năm 46 3.4.1 Tảo ảnh số độ ẩm khơng khí từ ảnh MODIS - MOD07 46 3.4.2 Các ảnh tổ hợp số độ ẩm khơng khí theo tháng, theo mùa ảnh số độ ẩm theo năm 47 3.5 Xây dựng tổ hợp NDVI trung bình theo năm 50 3.5.1 Tạo ảnh số thực vật từ ảnh MODIS - MOD13A2 50 3.5.2 Các ảnh tổ hợp số thực vật NDVI theo tháng, theo mùa theo năm 51 3.6 Xây dựng mối tương quan nhiệt độ, độ ẩm lớp phủ thực vật 54 3.6.1 Phân tích khơng gian đánh giá mối tương quan nhiệt độ, độ ẩm đến lớp phủ thực vật thông qua số thực vật (NDVI) 54 3.6.2 Phân tích hồi quy 56 3.6.3 Sơ đồ khối phân tích, đánh giá tác động nhiệt độ bề mặt, độ ẩm khơng khí đến lớp phủ thực vật thơng qua số NDVI 58 iii LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 3.6.4 Mối tương quan nhiệt độ bề mặt (LST), độ ẩm khơng khí số thực vật (NDVI) 59 3.6.5 Thành lập đồ lớp phủ thực vật với tương quan LST-NDVI RHNDVI 65 3.7 Phân tích mối tương quan nhiệt độ, độ ẩm lớp phủ thực vật 70 KẾT LUẬN .71 TÀI LIỆU THAM KHẢO 72 iv LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com DANH MỤC CHỮ VIẾT TẮT GIS Geographic Information System (Hệ thống thông tin địa lý) MODIS Moderate Resolution Imaging Spectroradiometer (Đầu đo ảnh viễn thám độ phân giải trung bình gắn vệ tinh TERRA AQUA) NDVI Normalized Difference Vegetation Index (chỉ số thực vật) LST Land surface temperature (Nhiệt độ bề mặt) RH Relative humidity (Độ ẩm tương đối) HDF-EOS Hierarchical Data Format Earth Observing System (Định dạng ảnh MODIS) v LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com MỤC LỤC BẢNG BIỂU Bảng 2.1: Các thông số kỹ thuật ảnh MODIS 20 Bảng 2.2: Các ứng dụng kênh ảnh MODIS 21 Bảng 3.1: Các kênh phổ MODIS 41 Bảng 3.2: Sai số nhiệt độ bề mặt (LST) số liệu trạm quan trắc số liệu tính tốn ảnh MODIS 60 Bảng 3.3: Sai số độ ẩm khơng khí tương đối (RH) số liệu trạm quan trắc số liệu tính tốn ảnh MODIS 63 MỤC LỤC HÌNH Hình 1.1: Sơ đồ khơng gian nhiệt độ bề mặt – số thực vật quan hệ ý niệm với bay hơi, thoát nước phần trăm lớp phủ thực vật Hình 1.2: Chỉ số TVDI pixel ảnh [ Ts, NDVI] xác định tỷ lệ đường A = (Ts – Tsmin) B = (Tsmax – Tsmin) Hình 1.3: Biểu đồ quan hệ T – NDVI ảnh Landsat 2002 10 Hình 1.4: Biểu đồ quan hệ T – NDVI ảnh Aster 2006 10 Hình 1.5: Giá trị rìa khơ sử dụng tính tốn số TVDI – giá trị Ts cực đại xác định cho khoảng NDVI nhỏ rìa khơ xác định hồi quy tuyến tính 12 Hình 1.6: Sự phân bố không gian theo số TVDI cho vùng nghiên cứu, hình chữ nhật khu vực lấy mẫu rừng Đắc Lắk vòng tròn khu vực khơ hạn Bình Thuận 13 Hình 1.7: Đồ thị phân tán tương quan số TVDI từ liệu MODIS với số độ ẩm từ liệu Landsat ETM + kết hợp với mơ hình SVAT 14 Hình 1.8: Sự thay đổi theo thời gian số TVDI khu vực rừng Đắk Lắk khu vực khô hạn Bình Thuận qua mùa khơ 2000-2001; 2001-2002 2002-2003 14 Hình 2.1: Đồ thị đất (Soil line) 33 vi LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Hình 3.1: Vị trí hành Tây Nguyên đồ Việt Nam 36 Hình 3.2: Giao diện phần mềm ENVI 42 Hình 3.3: Tool chọn sản phẩm LST ảnh modis MOD11A2 43 Hình 3.4: Các thuật tốn tính tổ hợp ảnh MODIS - MOD11A2 43 Hình 3.5: Các ảnh tổ hợp nhiệt độ bề mặt năm 2012 46 Hình 3.6: Tool chọn sản phẩm để tính độ ẩm tương đối ảnh MODIS - MOD07 47 Hình 3.7: Các thuật tốn tính tổ hợp ảnh MODIS - MOD07 47 Hình 3.8: Các ảnh tổ hợp độ ẩm khơng khí tương đối năm 2012 50 Hình 3.9: Các thuật tốn tính tổ hợp ảnh MODIS – MOD13A2 51 Hình 3.10: Các ảnh tổ hợp số thực vật NDVI năm 2012 53 Hình 3.11: Các mơ hình vector raster 55 Hình 3.12: Vị trí trạm quan trắc khí tượng khu vực Tây Nguyên 59 Hình 3.13: Biểu đồ quan hệ LST - NDVI mùa khô (a); mùa mưa (b) 61 Hình 3.14: Biểu đồ quan hệ LST - NDVI năm 2012 61 Hình 3.15: Biểu đồ quan hệ RH - NDVI mùa khô (c); mùa mưa (d) 64 Hình 3.16: Biểu đồ quan hệ RH - NDVI năm 2012 64 Hình 3.17: Biểu đồ quan hệ (e) LST-NDVI (f) RH-NDVI đối tượng lớp phủ thực vật vào mùa khô 66 Hình 3.18: Biểu đồ quan hệ (g) LST-NDVI (h) RH-NDVI đối tượng lớp phủ thực vật vào mùa mưa 66 Hình 3.19: Biểu đồ quan hệ (i) LST-NDVI (k) RH-NDVI đối tượng lớp phủ thực vật năm 2012 67 Hình 3.20: Bản đồ lớp phủ thực vật Tây Nguyên mùa khơ; mùa mưa năm 2012 68 Hình 3.21: Bản đồ lớp phủ thực vật khu vực Tây Nguyên năm 2012 69 vii LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com MỞ ĐẦU Tính cấp thiết đề tài Rừng tài nguyên quan trọng Tây Nguyên, trình phát triển, rừng tài nguyên bị đe dọa tàn phá nhiều Trong năm gần đây, nỗ lực bảo vệ rừng Tây Nguyên có nhiều kết đáng ghi nhận, từ chủ trương tới hành động cụ thể Cùng với áp lực phá rừng lấy đất trồng công nghiệp, nỗ lực phát triển rừng đem lại thay đổi đáng kể diện tích chất lượng rừng Tây Nguyên Tuy vậy, đơn vị quản lý rừng chưa có cơng cụ hữu hiệu hỗ trợ cho công tác quản lý, giám sát chất lượng rừng, phản ánh xác biến động rừng khu vực quản lý, đặc biệt khu vực trọng điểm rừng đầu nguồn Với đặc thù địa hình, khí hậu sử dụng đất, rừng Tây Ngun có vai trị quan trọng môi trường phát triển kinh tế xã hội Với mơi trường rừng góp phần quan trọng việc điều hịa khí hậu, đảm bảo cân sinh thái, làm giảm nhẹ sức tàn phá khốc liệt thiên tai; Ngăn chặn xói mịn đất, bảo vệ mực nước ngầm; Là nơi lưu giữ nguồn gien động thực vật quý hiến Với phát triển kinh tế xã hội rừng nguồn cung cấp gỗ phục vụ cho nhu cầu đời sống sản xuất; Cung cấp mặt hàng lâm sản có giá trị xuất góp phần thu ngoại tệ phục vụ qua trình cơng nghiệp hóa, đại hóa đất nước; Cung cấp dược liệu quý có tác dụng chữa bệnh, cao sức khỏe cho người Ngoài rừng Tây Nguyên đứng trước thách thức lớn địi hỏi phải có nỗ lực giám sát, quản lý thường xuyên Trong điều kiện Tây Nguyên, sử dụng lực lượng có chưa thể đáp ứng yêu cầu nên hệ thống có khả cung cấp kịp thời, liên tục thông tin giám sát quản lý rừng cần thiết LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com (i) (k) Hình 3.19: Biểu đồ quan hệ (i) LST-NDVI (k) RH-NDVI đối tượng lớp phủ thực vật năm 2012 Từ biểu đồ quan hệ ta nhân thấy vùng đất ngập nước thường có nhiệt độ ấm vùng nước (ao, hồ) lớp phủ thực vật bên chúng có NDVI cao nước Đất cỏ, bụi có nhiệt độ bề mặt cao đất ngập nước Đất nơng nghiệp có NDVI cao Cây bụi, cỏ phụ thuộc vào tình trạng thủy lợi, giai đoạn tang trưởng mà có nhiệt độ, độ ẩm thay đổi Đối với rứng hỗn giao, kim, tre nứa rừng trồng có nhiệt độ bề mặt trung bình số NDVI, độ ẩm trung bình Rừng tự nhiên có nhiệt độ thấp, độ ẩm cao số NDVI cao 67 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Hình 3.20: Bản đồ lớp phủ thực vật Tây Nguyên mùa khô; mùa mưa năm 2012 68 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Hình 3.21: Bản đồ lớp phủ thực vật khu vực Tây Nguyên năm 2012 69 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 3.7 Phân tích mối tương quan nhiệt độ, độ ẩm lớp phủ thực vật Nhiệt độ bề mặt gián tiếp đến thực vật thơng qua q trình hóa học, sinh học lí học đất Nếu nhiệt độ đất cao thúc đẩy hoạt động sinh vật đất, đặc biệt hoạt động phân giải chất hữu vi sinh vật để cung cấp chất dĩnh dưỡng cho thực vật, mặt khác nhiệt độ đất ảnh hưởng đến khả trao đổi hấp thụ lý – hóa đất, qua cung cấp chất dinh dưỡng cách thuận lợi Đối với lớp phủ thực vật, nhiệt độ có ảnh hưởng đến hình thái, chức sinh lý khả sinh sản nhiệt độ thấp có ảnh hưởng đến hình thái Với hỗ trợ thông tin nhiệt độ bề mặt, độ ẩm khơng khí thực qua việc xem xét biểu đồ phân tán điểm phổ quan hệ LST - NDVI biểu đồ phân tán điểm phổ quan hệ RH - NDVI, lớp có đặc điểm tương tự gộp vào nhóm Phương pháp đặc biệt thích hợp cho trường hợp khơng biết rõ trước khu vực nghiên cứu thiếu quan trắc thực địa song hành vào thời điểm chụp ảnh Từ hình 3.17, 3.18, 3.19, 3.20 3.21 ta thấy Vào mùa khô nhiệt độ bề mặt cao, độ ẩm thấp, ngày khơ kéo dài nhiều tháng, hàm lượng nước bị mất, dẫn đến rụng khô, lúc trạng lớp phủ rừng thay đổi rõ rệt Vào mùa mưa, lượng mưa cao độ ẩm cao, sinh trưởng nhanh, hàm lượng nước ln cao, khơng có biểu tượng cho thấy bị thiếu nước cây, giai đoạn trạng thái lớp phủ rừng tốt Với liệu MODIS đa thời gian, ta phát hiện, theo dõi trạng thái cây, đặc biệt theo dõi hàm lượng nước thực vật thời điểm khác nhau, phục vụ cho quản lý giám sát thực vật phạm vi lớn Từ việc thành lập đồ số thực vật kết trên, ta áp dụng vào công tác quản lý, giám sát tài nguyên đất, rừng, nước Các ứng dụng đồ lớp phủ thực vật công tác quản lý xanh đô thị, giám sát phát triển lúa giám sát khô hạn… 70 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com KẾT LUẬN Nghiên cứu xây dựng hồn chỉnh bước tính tốn nhiệt độ bề mặt, độ ẩm khơng khí số thực vật từ ảnh MODIS phần mềm viễn thám GIS đơn giản Kết số liệu nhiệt bề mặt, độ ẩm khơng khí tính tốn bước đầu cho thấy có độ tin cậy cao thể mối tương quan chặt với liệu NDVI Đây nghiên cứu góp phần xác định mối tương quan số thực vật với yếu tố khí hậu như: nhiệt độ bề mặt (LST), độ ẩm khơng khí (RH) khu vực có số NDVI thấp ảnh hưởng đến sản xuất nông nghiệp, lâm nghiệp liên quan đến biến đổi khí hậu Ảnh vệ tinh MODIS cho thấy có ưu việt đáng kể sử dụng công tác giám sát rừng Tây Ngun Tuy có độ phân giải khơng gian mức trung bình ảnh vệ tinh lại có ưu lớn số lượng kênh phổ, tần suất chụp ảnh lặp lại độ trùm phủ lớn Kết hợp số nhiệt độ bề mặt, độ ẩm khơng khí số thực vật với cách hợp lý, hiệu cho phép xây dựng đồ lớp phủ thực vật biến động lớp phủ rừng cho toàn khu vực Tây Nguyên Hơn nữa, thông số môi trường rừng số thực vật, tăng trưởng thực vật, nhiệt độ bề mặt, độ ẩm đất, lượng bốc tính tốn cập nhật liên tục cho tồn khu vực nguồn thông tin bổ trợ cho nghiên cứu, giám sát rừng khác dự báo cháy, phát sâu bệnh, hay phát sớm rừng, suy thối rừng Ứng dụng cơng nghệ viễn thám theo dõi lớp phủ rừng nghiên cứu giới từ nhiều năm Với xuất ngày nhiều hệ thống vệ tinh quan sát Trái đất, việc kết hợp nhiều loại ảnh vệ tinh, đa độ phân giải nghiên cứu theo dõi lớp phủ rừng tiến hành Đặc biệt, kết hợp tư liệu ảnh vệ tinh MODIS với loại ảnh vệ tinh khác có độ phân giải cao Landsat, SPOT, hay VNREDsat Việt Nam giúp cải thiện đáng kể chất lượng phân loại rừng cung cấp thông tin cho việc theo dõi, giám sát, đánh giá biến động lớp phủ rừng quy mô lớn 71 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com TÀI LIỆU THAM KHẢO Tiếng Việt Trần Hùng (2007), “Sử dụng tư liệu MODIS theo dõi độ ẩm đất / thực vật bề mặt: Thử nghiệm với số mức độ khô hạn nhiệt độ - Thực vật (TVDI)”, Tạp chí Viễn thám Địa tin học, Số – 4/2007, tr 38-45 Huỳnh Thị Thu Hương, Trương Chí Quang, Trần Thanh Dân, (2012), “Ứng dụng ảnh MODIS theo dõi thay đổi nhiệt độ bề mặt đất tình hình khơ hạn vùng Đồng Bằng Sơng Cửu Long”, Tạp chí Khoa học-Trường Đại học Cần Thơ, 24a, 49-59 Dương Văn Khảm, Chu Minh Thu, “Ứng dụng ảnh vệ tinh Terra-Aqua (MODIS) việc tính tốn độ ẩm khơng khí độ phân giải cao”,Viện Khí tượng Thủy văn Nguyễn Xuân Lâm (2013), “Nghiên cứu ứng dụng công nghệ viễn thám để giám sát cảnh báo tác động biến đổi khí hậu nhằm chủ động phòng tránh giảm thiểu thiệt hại tai biến thiên nhiên”, Chương trình Khoa học Cơng nghệ cấp nhà nước biến đổi khí hậu, KHCN-BĐKH/11-15, Cục Viễn thám quốc gia Trần Công Minh (2007), “Khí hậu khí tượng đại cương”, Nhà xuất Đại học quốc gia Hà Nội Nguyễn Ngọc Thạch (2005), “Cơ sở viễn thám”, Trường Đại học khoa học Tự nhiênĐại học Quốc gia Hà Nội, Nhà xuất Nông nghiệp, Hà Nội Trần Thị Vân (2006), “Ứng dụng viễn thám nhiệt khảo sát đặc trưng nhiệt độ bề mặt đô thị với phân bố kiểu thảm phủ TP.Hồ Chí Minh”, Viện Mơi trường Tài nguyên, ĐHQG Hồ Chí Minh Trần Thị Vân, Nguyễn Hằng Hải (2011), “Quan hệ nhiệt số thực vật phân loại lớp phủ phục vụ đánh giá biến động đất đô thị”, Kỷ yếu Hội nghị Khoa học Công nghệ lần thứ 12, Hồ Chí Minh 72 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Tiếng Anh Bo-Cai Gao and Kaufman, Y J, (2003), “Water vapor retrievals using Moderate Resolution Imaging Spectroradiometer (MODIS) near-infrared channels Journal of Geophysical Research” 10 Brian D Wardlow (2007), “Analysis of Time-Series MODIS 250 m Vegetation Index Data for Crop Classification in the U.S.Central Great Plains”, University of Nebraska – Lincoln 11 Han.Y, Y Wang, Y Zhao (2010), “Estimating soil moisture conditions of the Greater Chanbai Moutain by Land Surface Temperature and NDVI”, IEE Transactions on Geoscience and Remote Sensing, v 48, n 6, p 2509-2515 12 Lambin E F, Ehrlich D (1996) The surface temperature-vegetation index space for land cover and land-cover change analysis International Journal of Remote Sensing, 17: 463–487 13 Ramon Solano, Kamel Didan (2010), “MODIS Vegetation Index User’s Guide (MOD13 Series)”, The University of Arizona 14 Sandholt I., Rasmussen K & Andersen J (2002) “A simple interpretation of the surface temperature/vegetation index space for assessment of surface moisture status”, Remote Sensing of Environment, Vol 79, pp 213-224 15 Wan Z, P Wang and L.X (2004), “Using MODIS Land surface temperature and Normalized Diference Vegetation index products for monitoring dought in the southern Great Plains”, USA, International Journal of remote sensing, v.25, 61-72 16 Zhengming Wan (2009), “MODIS Land Surface Temperature Products Users' Guide (MOD11 Series)”, ICESS, University of California, Santa Barbara 17 Website: https://wist.echo.nasa.gov 18 Website: http://earthexplorer.usgs.gov 19 http://vi.wikipedia.org/wiki/TâyNguyên 73 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com PHỤ LỤC PHỤ LỤC ẢNH CHỤP THỰC ĐỊA 74 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com PHỤ LỤC CÁC YẾU TỐ NHIỆT ĐỘ BỀ MẶT, ĐỘ ẨM KHƠNG KHÍ VÀ CHỈ SỐ THỰC VẬT NDVI ĐƯA VÀO QUÁ TRÌNH HỒI QUY TỪNG BƯỚC Độ ẩm khơng khí Tháng/Năm 10 11 12 Đơn vị : % Đà Lạt 2012 83 82 81 86 87 87 87 87 90 85 85 82 85.26737 Liên Khương 2012 76 74 75 78 81 83 85 84 89 83 83 78 80.93171 Bảo Lộc 2012 81 80 79 83 89 91 90 89 91 85 85 78 85.0567 Mdrăk 2012 87 84 81 79 78 71 74 72 85 84 87 86 80.76606 EakMat 2012 82 78 78 82 82 84 84 85 87 84 84 82 82.65539 Buôn Ma Thuật 2012 80 76 76 81 82 84 87 87 89 84 82 80 82.33333 EaHLeo 2012 83 79 79 82 83 91 90 89 88 84 84 82 84.40177 Buôn Hồ 2012 89 83 81 80 82 87 87 86 90 88 89 88 85.76053 Lăk 2012 79 75 76 79 78 80 83 81 86 81 81 78 79.77105 Pleiku 2012 78 77 76 80 83 90 90 90 88 81 80 76 82.33289 An Khê 2012 88 84 82 81 78 78 78 80 87 85 87 85 82.61816 YaLy 2012 76 77 78 84 85 89 90 90 91 83 82 74 83.20665 AyunPa 2012 80 76 72 76 76 76 78 78 85 81 81 78 78.16451 Đăk Tô 2012 73 71 74 80 83 90 90 92 90 85 85 78 82.76222 Kon Tum 2012 66 65 68 72 79 86 85 85 84 77 73 70 75.86328 Đăk Mil 2012 85 78 78 81 83 83 85 84 89 86 84 82 83.10239 Đăk Nông 2012 79 77 77 84 85 87 88 88 90 84 85 79 83.512 75 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Nhiệt độ bề mặt Tháng/Năm 10 11 12 Đơn vị: 0C Đà Lạt 2012 20.0 20.3 22.0 22.5 22.4 21.7 21.2 21.5 20.9 21.0 20.8 20.8 21.3 Liên Khương 2012 25.0 27.9 28.5 29.3 27.4 26.8 25.9 26.8 24.5 25.3 27.0 26.6 26.8 Bảo Lộc 2012 24.7 26.7 27.8 27.7 28.0 26.2 26.5 26.3 25.6 27.1 25.6 27.4 26.6 Mdrăk 2012 22.9 25.9 29.5 33.2 33.3 31.8 30.6 31.2 28.5 28.2 27.6 25.9 29.1 EakMat 2012 26.6 30.2 31.5 31.3 32.0 28.9 27.9 28.8 28.2 28.4 28.9 27.9 29.2 Buôn Ma Thuật 2012 26.6 30.1 30.4 31.0 31.3 28.4 27.8 29.0 27.4 28.6 29.6 28.5 29.1 EaHLeo 2012 24.5 28.2 29.6 28.9 29.5 26.0 25.4 25.7 25.2 27.0 27.6 25.8 26.9 Buôn Hồ 2012 25.1 28.8 29.4 32.5 31.1 27.6 27.2 27.2 27.0 27.7 29.0 27.2 28.3 Lăk 2012 27.7 31.0 32.3 32.7 31.8 29.8 29.0 29.2 28.9 29.6 30.4 28.4 30.1 Pleiku 2012 25.0 27.7 30.5 30.2 29.7 25.6 26.1 25.0 25.7 26.4 26.8 26.9 27.1 An Khê 2012 22.2 26.1 29.3 31.1 32.0 29.5 30.0 29.7 28.0 27.4 27.2 25.1 28.1 YaLy 2012 28.7 32.6 33.3 33.7 32.5 28.7 29.6 29.2 29.6 31.2 32.3 30.5 31.0 AyunPa 2012 27.1 30.8 33.4 35.5 35.2 32.4 32.4 32.3 30.4 31.0 31.1 29.2 31.7 Đăk Tô 2012 26.4 29.7 29.2 30.6 30.8 26.3 27.0 26.5 27.5 27.8 28.1 28.0 28.2 Kon Tum 2012 27.1 30.5 29.5 31.9 30.8 27.7 28.3 27.5 27.7 29.2 30.1 28.9 29.1 Đăk Mil 2012 24.3 28.0 27.5 27.1 28.4 27.9 26.1 27.1 25.6 26.5 28.7 26.7 27.0 Đăk Nông 2012 26.8 29.6 30.1 30.9 31.0 27.6 27.2 27.3 26.9 28.7 28.4 28.4 28.6 76 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Tọa độ X Tọa độ Y NDVI LST RH Tọa độ X Tọa độ Y NDVI LST RH 783029.3 1700859.76 0.819 23.28 82 782623.66 1700132.09 0.854 23.38 85 784093.01 1685258.65 0.799 20.91 80 785023.02 1687792.53 0.805 22.53 792602.71 1682776.65 0.842 20.56 83 794730.14 1672848.67 0.651 27.44 807849.26 813876.96 Tọa độ X Tọa độ Y LST RH 0.84 21.93 83 83 785324.6 1676018.51 0.734 22.10 74 793249.39 1684022.11 0.727 21.62 77 805354.99 1675690.14 0.853 21.59 86 72 786736.85 1675452.97 0.73 21.96 78 814549.26 1677003.61 0.846 22.89 83 1672494.1 0.852 20.45 83 794963.22 1676481.27 0.793 24.18 81 783026.03 1662555.46 0.826 21.27 82 1661502.4 0.849 17.65 80 805931.72 1676481.27 0.89 22.30 83 796489.08 1663212.2 0.673 27.02 68 0.51 31.62 65 816557.45 1675795.74 0.831 22.83 80 807653.56 1663212.2 0.786 19.02 79 799339.56 1648383.28 0.853 21.98 82 778167.71 1657629.17 0.769 22.18 77 817504.57 1660913.63 0.832 17.01 85 809622.11 1645192.14 0.622 29.12 72 790164.5 1656600.87 0.788 24.52 78 788279.9 1644166.91 0.624 28.19 64 826641.51 1645192.14 0.655 26.08 74 803189.59 1656600.87 0.762 24.84 75 800101.11 1643510.18 0.796 22.44 80 839406.06 1645546.71 0.865 21.75 82 800185.38 1637476.42 0.724 20.48 73 831952.71 1641539.98 0.824 21.67 79 Mùa 854652.61 1645546.71 0.858 21.7 80 Mùa 810284.39 1635386.97 0.694 26.74 70 838191.69 1646137.11 0.84 20.52 77 khô 781611.02 1620726.76 0.49 32.08 64 mưa 828392.97 1634690.48 0.787 21.31 76 Năm 849027.8 1647122.21 0.876 22.77 81 799339.56 1621081.33 0.492 31.21 64 836750.78 1634342.24 0.801 20.15 79 859863.91 1632017.33 0.822 22.32 80 816713.53 1621790.47 0.479 31.42 64 846849.8 1635038.73 0.855 21.79 81 844759.03 79 833023.79 81 855555.84 1636779.94 0.85 20.99 82 833266.18 1628405.29 0.79 23.83 80 20.9 83 872967.94 1632601.03 0.824 22.65 79 820131.5 1626763.46 0.62 27.06 71 867417.16 1620017.61 0.762 22.86 71 808891.43 1609268.82 0.61 26.42 72 805026.62 1625778.36 0.661 28.62 77 790829.86 1600516.22 0.853 22.73 80 775460.2 1609617.07 0.707 22.91 72 794518.88 1625121.62 0.633 29.49 73 802176.13 1600516.22 0.598 29.58 70 831875.39 1609268.82 0.767 21.72 75 782369.3 1624464.89 0.637 28.40 73 815649.82 1599452.51 0.643 29.69 72 847894.52 1611358.28 0.835 22.35 84 774816.86 1612315.31 0.808 23.69 77 828059.8 1599097.94 0.603 31.37 70 862520.68 1612403 0.761 21.14 75 797145.81 1610345.11 0.636 28.99 76 844370.05 1599807.08 0.431 30.74 63 876798.6 1610661.79 0.854 21.14 84 809623.76 1610345.11 0.642 28.28 74 879118 1598743.36 0.845 22.53 77 876450.36 1599866.29 0.859 21.94 86 824400.27 1608703.27 0.638 29.45 74 769555.61 1582433.11 0.652 29.47 76 867396.07 1582105.96 0.815 20.73 82 839833.52 1608046.54 0.805 23.81 82 789766.15 1642355.58 1621435.9 0.853 22.89 847206.62 1620726.76 0.874 788608.27 1689809.93 NDVI TB 1629390.4 0.858 21.26 77 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 797566.71 1580660.25 0.663 28.27 78 872271.45 1558425.51 0.812 21.92 82 855923.51 1607718.17 0.565 27.91 72 822032.09 1577469.12 0.647 30.00 72 836750.78 1558773.75 0.658 27.75 66 868401.45 1607718.17 0.866 22.09 83 836214.93 1577114.54 0.58 32.28 71 809587.91 1569569.25 0.661 29.07 68 880222.66 1607061.44 0.829 22.04 84 856425.46 1577823.69 0.871 22.78 81 765012.95 1563300.89 0.683 27.98 69 873655.32 1593926.76 0.793 23.87 79 877699.71 1576759.97 0.789 25.10 75 760834.04 1552157.15 0.641 28.86 67 865774.52 1590314.72 0.858 21.74 84 896491.97 1559031.43 0.504 31.23 66 770236.57 1549022.97 0.675 27.42 68 835564.75 1591956.56 0.621 29.94 75 0.83 22.30 78 793220.54 1550415.94 0.739 23.93 76 825385.37 1591628.19 0.572 30.13 68 845788.34 1544848.6 0.637 30.95 71 808543.18 1547281.76 0.695 28.44 70 814549.26 1591956.56 0.51 29.65 65 829832.65 1544848.6 0.667 29.73 73 826651.76 1544844.07 0.517 30.54 65 804041.52 1591628.19 0.537 30.17 65 803948.98 1543784.89 0.636 31.82 72 842670.89 1543102.86 0.698 27.53 73 778100.53 1590643.09 0.677 28.28 77 777001.6 1541657.46 0.662 31.80 73 859734.75 1544147.59 0.818 21.58 83 764965.84 1571597.8 0.561 29.64 71 780547.31 1523928.92 0.549 33.19 64 881673.99 1540316.93 0.675 29.06 69 775145.22 1570612.7 0.759 23.92 79 797212.13 1523928.92 0.665 31.56 76 898041.36 1538227.47 0.697 28.54 70 786309.7 1569627.6 0.553 29.99 67 824159.52 1523928.92 0.669 30.14 76 897344.87 1525342.52 0.729 22.67 74 795175.61 1567985.77 0.736 22.18 74 843306.34 1523928.92 0.619 31.47 75 882022.23 0.68 27.84 72 807981.92 1566672.3 0.636 29.81 67 861389.45 1522156.07 0.751 25.27 75 869137.28 1521163.62 0.669 27.51 72 820788.24 1566672.3 0.575 29.95 65 876636 1518964.93 0.629 28.37 72 853814.63 1519770.65 0.698 27.08 78 843117.19 1567000.66 0.644 28.92 76 1516837.5 0.575 30.04 63 837795.51 1517332.96 0.672 28.46 78 851983.1 1566343.93 0.828 23.60 79 907483.66 1513646.37 0.776 24.82 81 816900.99 1512457.57 0.641 29.07 75 867416.35 1564045.36 0.859 22.80 84 788702.43 1488471.84 0.319 34.60 60 796354.72 1515939.99 0.669 28.12 77 893357.35 1560433.32 0.554 29.42 69 806076.41 1489180.98 0.481 33.36 61 788693.4 1496090.2 0.841 23.05 85 903536.72 1552552.52 0.596 29.46 72 821677.52 1488471.84 0.357 35.61 64 806453.73 1491911.3 0.621 32.95 75 885148.17 1550253.95 0.575 30.42 72 836569.5 1489180.98 0.407 34.58 61 818642.2 1490170.09 0.609 31.65 74 869058.19 1547298.64 0.856 23.53 82 1487762.7 0.637 28.88 71 824910.55 1489821.85 0.522 31.12 72 858878.81 1545656.81 0.831 22.85 81 871672.01 1487053.55 0.537 30.82 67 838491.99 1489821.85 0.649 30.18 78 843117.19 1543686.61 0.587 30.50 76 898973.96 1478898.42 0.541 31.72 67 859038.26 67 825057.01 1547298.64 0.267 31.92 66 859616.6 1545203.17 896137.4 850397.76 1521860.1 1496090.2 0.5 30.09 78 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 912802.22 1475707.29 0.812 24.08 80 875405.63 1495741.96 0.579 30.15 67 811922.33 1546641.91 0.638 29.17 76 787284.15 1462588.17 0.46 32.74 64 898389.6 1501662.07 0.712 22.62 73 797802.55 1543686.61 0.637 29.96 79 799694.13 1461169.88 0.36 36.72 63 870182 1491563.06 0.542 28.33 65 774488.49 1545985.18 0.502 28.57 71 36.1 69 836054.3 1488428.88 0.581 29.04 65 780070.73 1533835.6 0.687 28.00 76 828768.94 1456560.46 0.547 32.44 78 818293.96 1489821.85 0.594 31.55 66 807325.19 1530880.29 0.609 29.34 74 842951.77 1456205.89 0.62 30.93 77 794613.51 1482857.01 0.697 25.89 73 853624.94 1528253.36 0.589 30.14 69 856780.03 1456560.46 0.645 30.64 75 1461962.5 0.502 28.96 63 861177.38 1533835.6 0.812 23.68 79 873444.86 1455851.32 0.703 25.49 78 853466.39 1423655.88 0.543 28.48 65 874640.43 1531865.39 0.624 29.17 76 902519.67 1450887.33 0.481 32.29 62 826651.76 1418780.5 0.595 27.94 67 883506.33 1531208.66 0.561 30.14 69 893300.83 1437768.21 0.788 25.65 80 810284.39 1415994.56 0.602 29.41 68 898939.58 1528910.09 0.881 23.56 84 780547.31 1420748.81 0.643 28.78 72 786603.94 1400323.68 0.847 23.45 85 906820.39 1522014.38 0.785 22.94 76 809622.11 1418621.38 0.42 33.32 66 802623.07 1392662.35 0.535 29.51 73 881207.77 1514461.94 0.625 30.15 74 827350.65 1416848.53 0.698 30.97 77 857645.29 1389179.94 0.581 28.73 73 865446.15 1512820.11 0.567 29.94 72 857843.75 1416493.96 0.682 28.91 75 874012.66 1386742.24 0.61 30.84 71 841146.99 1510521.54 0.608 30.75 73 885145.7 1413302.82 0.602 28.89 78 888638.82 1381866.86 0.613 29.36 72 824071.91 1509208.07 0.691 28.98 78 924148.49 1410111.68 0.832 23.50 83 902220.26 1365499.49 0.849 20.81 83 807653.56 1506581.13 0.665 29.15 77 884436.56 1384937.15 0.635 29.2 76 890380.03 1355052.23 0.845 20.68 81 794847.24 1503297.46 0.451 32.41 64 867417.16 1381746.02 0.859 22.27 84 824910.55 78 788936.64 1494759.92 0.784 24.91 78 830187.22 1381036.88 0.478 29.21 65 807150.22 1358882.89 0.78 22.75 80 797474.18 1483923.81 0.607 28.79 74 812104.11 1381391.45 0.579 29.92 72 797051.2 1357489.92 0.798 21.41 80 811593.96 1484908.91 0.452 32.04 63 792957.28 1378909.45 0.696 28.75 77 745859.64 1358534.65 0.758 22.14 78 825057.01 1484908.91 0.599 31.06 70 778774.45 1360471.77 28.8 75 750386.79 1334854.2 0.73 23.71 74 840490.26 1484908.91 0.531 31.21 70 796857.56 1357989.77 0.854 22.54 83 783469.77 1331720.02 0.68 26.99 72 858222.08 1483595.44 0.674 28.21 76 822741.23 1355507.78 0.627 28.85 76 808194.94 1327889.36 0.802 23.48 80 871356.76 1482610.34 0.537 30.59 74 877345.14 1345579.79 0.845 19.3 83 825607.04 1325103.42 0.829 23.86 81 883834.7 1481953.61 0.752 25.63 75 895428.25 1347352.65 0.831 21.0 83 841626.17 1320924.52 0.649 24.35 68 901238.15 1481953.61 0.685 28.87 74 816004.39 1457269.6 0.359 0.62 803667.8 1356445.2 0.745 24.54 79 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com MODIS Data Products Table MCD12C1 Combined Land Cover Type CMG 5600m Yearly MCD12Q1 Combined Land Cover Type Tile 500m Yearly MCD12Q2 Combined Land Cover Dynamics Tile 500m Yearly MCD15A2 Combined Leaf Area Index - FPAR Tile 1000m day MCD15A3 Combined Leaf Area Index - FPAR Tile 1000m day MCD43A1 Combined BRDF-Albedo Model Parameters Tile 500m 16 day MCD43A2 Combined BRDF-Albedo Quality Tile 500m 16 day MCD43A3 Combined Albedo Tile 500m 16 day MCD43A4 Combined Nadir BRDF-Adjusted Reflectance Tile 500m 16 day MCD43B1 Combined BRDF-Albedo Model Parameters Tile 1000m 16 day MCD43B2 Combined BRDF-Albedo Quality Tile 1000m 16 day MCD43B3 Combined Albedo Tile 1000m 16 day MCD43B4 Combined Nadir BRDF-Adjusted Reflectance Tile 1000m 16 day MCD43C1 Combined BRDF-Albedo Model Parameters CMG 5600m 16 day MCD43C2 Combined BRDF-Albedo Snow-free Quality CMG 5600m 16 day MCD43C3 Combined Albedo CMG 5600m 16 day MCD43C4 Combined Nadir BRDF-Adjusted Reflectance CMG 5600m 16 day MCD45A1 Combined Burned Area Tile 500m Monthly MOD09A1 MOD09CMG Terra Terra Surface Reflectance Bands 1–7 Surface Reflectance Bands 1–7 Tile CMG 500m 5600m day Daily MOD09GA Terra Surface Reflectance Bands 1–7 Tile 500/1000m Daily MOD09GQ Terra Surface Reflectance Bands 1–2 Tile 250m Daily MOD09Q1 Terra Surface Reflectance Bands 1–2 Tile 250m day MOD11A1 Terra Land Surface Temperature & Emissivity Tile 1000m Daily MOD11A2 Terra Land Surface Temperature & Emissivity Tile 1000m day MOD11B1 Terra Land Surface Temperature & Emissivity Tile 5600m Daily MOD11C1 Terra Land Surface Temperature & Emissivity CMG 5600m Daily MOD11C2 Terra Land Surface Temperature & Emissivity CMG 5600m day MOD11C3 Terra Land Surface Temperature & Emissivity CMG 5600m Monthly MOD11_L2 Terra Land Surface Temperature & Emissivity Swath 1000m MOD13A1 Terra Vegetation Indices Tile 500m 16 day MOD13A2 Terra Vegetation Indices Tile 1000m 16 day MOD13A3 Terra Vegetation Indices Tile 1000m Monthly MOD13C1 Terra Vegetation Indices CMG 5600m 16 day MOD13C2 Terra Vegetation Indices CMG 5600m Monthly MOD13Q1 Terra Vegetation Indices Tile 250m 16 day MOD14 Terra Thermal Anomalies & Fire Swath 1000m MOD14A1 Terra Thermal Anomalies & Fire Tile 1000m Daily 80 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com MOD14A2 Terra Thermal Anomalies & Fire Tile 1000m day MOD15A2 Terra Leaf Area Index - FPAR Tile 1000m day MOD17A2 Terra Gross Primary Productivity Tile 1000m day MOD17A3 Terra Net Primary Productivity Tile 1000m Yearly MOD44B Terra Vegetation Continuous Fields Tile 250m Yearly MOD44W Terra Land Water Mask Derived Tile 250m None MYD09A1 Aqua Surface Reflectance Bands 1–7 Tile 500m day MYD09CMG MYD09GA Aqua Aqua Surface Reflectance Bands 1–7 Surface Reflectance Bands 1–7 CMG Tile 5600m 500/1000m Daily Daily MYD09GQ Aqua Surface Reflectance Bands 1–2 Tile 250m Daily MYD09Q1 Aqua Surface Reflectance Bands 1–2 Tile 250m day MYD11A1 Aqua Land Surface Temperature & Emissivity Tile 1000m Daily MYD11A2 Aqua Land Surface Temperature & Emissivity Tile 1000m day MYD11B1 Aqua Land Surface Temperature & Emissivity Tile 5600m Daily MYD11C1 Aqua Land Surface Temperature & Emissivity CMG 5600m Daily MYD11C2 Aqua Land Surface Temperature & Emissivity CMG 5600m day MYD11C3 Aqua Land Surface Temperature & Emissivity CMG 5600m Monthly MYD11_L2 Aqua Land Surface Temperature & Emissivity Swath 1000m MYD13A1 Aqua Vegetation Indices Tile 500m 16 day MYD13A2 Aqua Vegetation Indices Tile 1000m 16 day MYD13A3 Aqua Vegetation Indices Tile 1000m Monthly MYD13C1 Aqua Vegetation Indices CMG 5600m 16 day MYD13C2 Aqua Vegetation Indices CMG 5600m Monthly MYD13Q1 Aqua Vegetation Indices Tile 250m 16 day MYD14 Aqua Thermal Anomalies & Fire Swath 1000m MYD14A1 Aqua Thermal Anomalies & Fire Tile 1000m Daily MYD14A2 Aqua Thermal Anomalies & Fire Tile 1000m day MYD15A2 Aqua Leaf Area Index - FPAR Tile 1000m day MYD17A2 MOD07 Aqua Terra, Aqua Gross Primary Productivity Atmospheric Profiles Tile Tile 1000m 1000m day day 81 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com ... PHẠM VĂN MẠNH NGHIÊN CỨU ỨNG DỤNG CÔNG NGHỆ VIỄN THÁM VÀ GIS ĐÁNH GIÁ TÁC ĐỘNG CỦA NHIỆT ĐỘ, ĐỘ ẨM ĐẾN LỚP PHỦ THỰC VẬT THÔNG QUA CHỈ SỐ THỰC VẬT (NDVI) KHU VỰC TÂY NGUYÊN Chuyên ngành: Bản đồ viễn. .. luận văn thạc sĩ: ? ?Nghiên cứu ứng dụng công nghệ viễn thám GIS đánh giá tác động nhiệt độ, độ ẩm đến lớp phủ thực vật thông qua số thực vật (NDVI) khu vực Tây Nguyên? ?? Mục tiêu nhiệm vụ nghiên cứu. .. cơng trình nghiên cứu tiêu biểu tác giả tiêu biểu giới nước vấn đề nghiên cứu đánh giá, phân tích tác động nhiệt độ, độ ẩm đến lớp phủ thực vật thông qua sổ thực vật dựa công nghệ viễn thám xem