Sử dụng ảnh viễn thám xác định ngưỡng chỉ số cho phép phát hiện sớm cháy rừng tại huyện tĩnh gia, tỉnh thanh hóa

LỜI CẢM ƠN Trƣớc tiên, em chân thành cảm ơn Ban giám hiệu Trƣờng Đại học Lâm Nghiệp Hà Nội, q thầy ngồi ngành Khoa học mơi trƣờng tận tình dạy, truyền đạt kiến thức quý báu cho em suốt bốn năm học Em xin chân thành cảm ơn thầy PGS.TS Nguyễn Hải Hịa, giảng viên Bộ mơn Kỹ thuật Mơi trƣờng - Khoa Quản lý tài nguyên Rừng Môi trƣờng, Trƣờng Đại học Lâm Nghiệp Hà Nội, ngƣời hƣớng dẫn, góp ý để em hồn thành Khóa luận Em xin chân thành cảm ơn Ban Quản lý Rừng Phịng Hộ huyện Tĩnh Gia, tỉnh Thanh hóa hỗ trợ, giúp đỡ để em hoàn thành Khóa luận Tơi xin cảm ơn giúp đỡ tập thể lớp K60 – QTNV ngày tháng ngồi dƣới giảng đƣờng đại học Tuy có nhiều cố gắng nhƣng thời gian trình độ chuyên mơn cịn hạn chế nên chắn Khóa luận cịn nhiều thiếu sót hạn chế định Kính mong nhận đƣợc ý kiến đóng góp thầy cô bạn Hà Nội, ngày10 tháng năm 2019 Sinh viên thực Nguyễn Trung Đức i MỤC LỤC LỜI CẢM ƠN i MỤC LỤC ii DANH MỤC CHỮ VIÊT TẮT v DANH MỤC CÁC BẢNG vi DANH MỤC CÁC HÌNH vii ĐẶT VẤN ĐỀ Chƣơng I TỔNG QUAN VẤN ĐỀ NGHIÊN CỨU 1.1 Viễn thám hệ thống thông tin địa lý (GIS) 1.1.1 Khái quát chung hệ thống thông tin địa lý 1.2 Khái quát chung viễn thám 1.2.1 Khái niệm viễn thám 1.2.2 Nguyên lý hoạt động viễn thám 1.3 Ứng dụng GIS viễn thám nghiên cứu giới Việt Nam 11 1.4 Tính cấp thiết vấn đề nghiên cứu 13 Chƣơng II MỤC TIÊU, NỘI DUNG VÀ PHƢƠNG PHÁP NGHIÊN CỨU 15 2.1 Mục tiêu nghiên cứu 15 2.1.1 Mục tiêu chung 15 2.1.2 Mục tiêu cụ thể 15 2.2 Phạm vi nghiên cứu 15 2.2.1 Phạm vi không gian: Nghiên cứu phạm vi huyện Tĩnh Gia Thanh Hóa 15 2.2.2 Phạm vi thời gian: Nghiên cứu giai đoạn từ năm 2010 - 2016 15 2.2.3 Phạm vi nội dung: Khảo sát, xác định thay đổi diện tích rừng, số đánh giá trƣớc sau cháy rừng theo không gian thời gian khu vực nghiên cứu 15 2.2.4 Phạm vi phƣơng pháp: 15 2.3 Nội dung nghiên cứu 15 2.3.1 Nghiên cứu trạng rừng hoạt động phòng chữa cháy rừng huyện Tĩnh Gia, tỉnh Thanh Hóa 15 ii 2.3.2 Nghiên cứu lựa chọn số viễn thám cho phép phát sớm cháy rừng khu vực nghiên cứu 16 2.3.4 Nghiên cứu đề xuất giải pháp nâng cao hiệu cơng tác phịng chữa cháy rừng khu vực nghiên cứu 16 2.4 Phƣơng pháp nghiên cứu 16 2.4.1 Kế thừa số liệu thứ cấp 16 2.4.2 Phƣơng pháp thu thập ngoại nghiệp xử lý số liệu 17 2.4.3 Nguyên nhân gây nên tình trạng cháy rừng cơng tác quản lý môi trƣờng khu vực nghiên cứu 23 2.4.4 Đề xuất giải pháp nâng cao hiệu quản lí chất lƣợng mơi trƣờng khơng khí khu vực nghiên cứu 23 Chƣơng III ĐIỀU KIỆN TỰ NHIÊN, KINH TẾ XÃ HỘI KHU VỰC NGHIÊN CỨU 25 3.1 Điều kiện tự nhiên xã hội huyện Tĩnh Gia 25 3.1.1 Vị trí địa lý 25 3.1.2 Khí hậu 26 3.1.3 Địa hình 27 3.1.4 Thủy văn 28 3.1.5.Tài nguyên đất 28 3.1.6 Nhiệt độ 29 3.2 Đặc điểm kinh tế - xã hội 29 3.2.1 Đặc điểm phát triển kinh tế 29 3.2.2 Đặc điểm phát triển xã hội 31 Chƣơng IV KẾT QUẢ NGHIÊN CỨU VÀ THẢO LUẬN 32 4.1 Hiện trạng rừng hoạt động phòng chữa cháy rừng khu vực huyện Tĩnh Gia, tỉnh Thanh Hóa 32 4.1.1 Hiện trạng rừng huyện Tĩnh Gia 32 4.1.2 Hoạt động ban quản lý rừng cơng tác phịng chữa cháy rừng địa bàn 33 4.1.3 Xây dựng đồ đánh giá độ xác trạng rừng huyện Tĩnh Gia 34 4.2 Lựa chọn xác định ngƣỡng số viễn thám cho phép phát sớm cháy rừng khu vực nghiên cứu 41 4.2.1 Lựa chọn số viễn thám cho phép phát sớm cháy rừng 41 iii 4.2.2 Xác định ngƣỡng số viễn thám cho phép phát sớm cháy rừng khu vực nghiên cứu 41 4.2.3 Đánh giá độ xác biến động diện tích cháy rừng từ ảnh đánh giá với số liệc thực địa 42 4.3 Xây dựng đồ ngƣỡng cháy với số viễn thám cho phép phát sớm cháy rừng khu vực nghiên cứu 44 4.4 Giải pháp phòng chống cháy rừng khu vực nghiên cứu 46 4.4.1 Giải pháp kỹ thuật 46 4.4.2 Giải pháp thể chế sách 47 4.4.3 Giải pháp kinh tế - xã hội 48 Phần V KẾT LUẬN, TỒN TẠI VÀ KIẾN NGHỊ 50 5.1 Kết luận 50 5.3 Kiến nghị 51 TÀI LIỆU THAM KHẢO PHỤ LỤC iv DANH MỤC CHỮ VIÊT TẮT GIS Hệ Thống Thông Tin Địa Lý (Geographic Information System) PCCCR Phòng Chống Chữa Cháy Rừng RS Viễn Thám(Remote Sensing) UBND Ủy Ban Nhân Dân VLC Vật Liệu Cháy v DANH MỤC CÁC BẢNG Bảng 1.1 Thông số Landsat TM Landsat 10 Bảng 2.1: Dữ liệu ảnh Landsat GIS sử dụng nghiên cứu 17 Bảng 3.1 Nhiệt độ qua tháng huyện Tĩnh Gia 29 Bảng 4.1 Diện tích rừng đất lâm nghiệp huyện Tĩnh Gia 32 Bảng 4.2 Thống kê diện tích loại hình sử dụng đất huyện Tĩnh Gia 36 Bảng 4.3 Đánh giá độ xác đồ trạng rừng huyện Tĩnh Gia tháng năm 2010 38 Bảng 4.4 Đánh giá độ xác đồ trạng rừnghuyện Tĩnh Gia tháng 10 năm 2010 38 Bảng 4.5 Đánh giá độ xác đồ trạng rừnghuyện Tĩnh Gia tháng năm 2015 39 Bảng 4.6 Đánh giá độ xác đồ trạng rừnghuyện Tĩnh Gia tháng năm 2015 39 Bảng 4.7 Đánh giá độ xác đồ trạng rừnghuyện Tĩnh Gia tháng năm 2016 39 Bảng 4.8 Đánh giá độ xác đồ trạng rừnghuyện Tĩnh Gia tháng 12 năm 2016 39 Bảng 4.9: Giá trị thống kê biến động số vị trí cháy rừng Tĩnh Gia giai đoạn 2010 - 2016 41 Bảng 4.10: Kết kiểm chứng biến động diện tích suy thối rừng từ ảnh với điều tra thực địa 42 vi DANH MỤC CÁC HÌNH Hình 1.1 Các thành phần GIS Hình 1.2 Sơ đồ nguyên lý thu nhận hình ảnh viễn thám Hình 3.1 Bản đồ ranh giới huyện Tĩnh Gia- tỉnh Thanh Hóa 26 Hình 4.1 Hiện trạng sử dụng đất năm 2010 (Landsat ETM+/2010) 35 Hình 4.2 Hiện trạng sử dụng đất năm 2015 (Landsat OLI/TILS/2015) 36 Hình 4.3 Hiện trạng sử dụng đất năm 2016 (Landsat OLI/TILS/2016) 36 Hình 4.4 Ngƣỡng số cho phép phát sớm cháy rừng năm 2010 45 Hình 4.5 Ngƣỡng số cho phép phát sớm cháy rừng năm 2015 45 Hình 4.6 Ngƣỡng số cho phép phát sớm cháy rừng năm 2016 46 vii ĐẶT VẤN ĐỀ Rừng hệ sinh thái có nhiều chức sinh học nhƣ mang lại giá trị kinh tế to lớn Có phân bố rộng khắp bề mặt đất có tổng hợp nhiều loài sinh vật, rừng đƣợc xem nơi chứa đựng lớn trái đất Là hệ sinh thái có tính đa dạng sinh học cao so sánh với hệ sinh thái lại Khu vực rừng mƣa nhiệt đới với việc bao phủ 7% diện tích mặt nhƣng lại chiếm tới 50 ÷ 70% tổng số loài sinh vật toàn giới Rừng chứa đựng nguồn đa dạng sinh học lớn giữ vai trị khơng thể thay đổi việc ổn định tính đa dạng sinh học Trái đất Việc phát yếu tố ảnh hƣởng tới rừng đƣợc coi trình quan trọng việc đánh giá tài nguyên đặc biệt lĩnh vực có liên quan đến vịng tuần hồn cácbon biến đổi khí hậu Cháy rừng ln đƣợc coi tác nhân ảnh hƣởng đến xếp đặc điểm sinh thái rừng, yếu tố hình thái cấu trúc nhƣ thuộc tính rừng, bao gồm: đặc điểm tính đa dạng sinh học nhƣ vịng tuần hồn sinh hóa lồi sinh vật Từ cho thấy việc xác định đƣợc khu vực mức độ tác động cháy rừng vô quan trọng Điều giúp hiểu thay đổi môi trƣờng khí hậu diễn khu vực Đồng thời mang lại nhƣng nguồn liệu hữu ích cho việc xác định xu hƣớng mặt không gian thời gian quản lý tài nguyên rừng Huyện Tĩnh Gia huyện ven biển thuộc tỉnh Thanh Hóa, nằm miền Trung Việt Nam Khu vực nơi tập trung lớn diện tích rừng tỉnh với nhiều loài thực vật quý hiếm, nhƣng khu vực chịu ảnh hƣởng lớn cháy rừng Cùng với việc tần suất xuất nhiều nhƣ hậu mà mang lại ngày khốc liệt ảnh hƣởng không nhỏ không với ban quản lý rừng với ngƣời dân sinh sống khu vực ảnh hƣởng mang lại Theo số liệu báo cáo tính năm 2010 có tới 12 vụ cháy diễn với diện tích rừng bị cháy khoảng 390,61ha hay năm 2015 114ha Những số liệu cháy rừng thu thập đƣợc chủ yếu thông qua việc đánh giá trực tiếp khu vực cháy Việc thực đánh giá theo phƣơng pháp không làm nhiều công sức, tiền bạc nhƣ có tính liên tục cập nhật mà cịn có hiệu chƣa thực cao Từ việc làm để đánh giá xác đƣợc liệu liên quan đến cháy rừng nhƣ: diện tích cháy ,vùng cháy cách xác có đƣợc cập nhật liên tục trở nên vô cần thiết quan trọng công tác quản lý rừng Ngày nay, việc ứng dụng viễn thám quản lý tài nguyên trở nên phổ biến đƣợc sử dụng rộng khắp nhiều lĩnh vực đặc biệt quản lý tài nguyên rừng ƣu điểm trội Việc ứng dụng liệu viễn thám phát sớm cháy rừng không giúp cho việc quản lý cháy rừng đơn giản mà đem lại hiệu cao so với phƣơng pháp thông thƣờng Xuất phát từ nhu cầu thực tiễn trên, đề tài: “Sử Dụng Ảnh Viễn Thám Xác Định Ngƣỡng Chỉ Số Cho Phép Phát Hiện Sớm Cháy Rừng Tại Huyện Tĩnh Gia, Tỉnh Thanh Hóa” đƣợc thực Chƣơng I TỔNG QUAN VẤN ĐỀ NGHIÊN CỨU 1.1 Viễn thám hệ thống thông tin địa lý (GIS) 1.1.1 Khái quát chung hệ thống thông tin địa lý Hệ thống thông tin địa lý (Geographic Information System ) bắt đầu đƣợc sử dụng rộng rãi nƣớc phát triểm thập niên qua, đâu dạng ứng dụng công nghệ tin học nhằm mô tả giới thực mà lồi sống, tìm hiểu, khai thác Với tính ƣu việt, kĩ thuật GIS ngày đƣợc ứng dụng nhiều lĩnh vực nghiên cứu quản lí, đặc biệt quản lí quy hoạch sử dụng, khai thác nguồn tài nguyên cách bền vững hợp lí Khái niệm hệ thống thơng tin địa lý: Hệ thống thông tin địa lý (GIS) ngành khoa học đƣợc ứng dụng rộng rãi nên có nhiều định nghĩa khác GIS Theo Ducke (1979) định nghĩa, GIS hệ thống thông tin, sở liệu bao gồm quan sát đặc trƣng phân bố không gian, hoạt động kiện xác định khoảng không nhƣ điểm, đƣờng, vùng Sự đời hệ thống thông tin địa lý: Lĩnh vực hệ thống thông tin địa lý (GIS) bắt đầu vào năm 1960 máy tính khái niệm ban đầu địa lý định lƣợng tính tốn xuất Công việc GIS ban đầu bao gồm nghiên cứu quan trọng cộng đồng học thuật Sau đó, Trung tâm Dữ liệu Phân tích Địa lý Quốc gia Michael Goodchild dẫn đầu thức nghiên cứu chủ đề khoa học thông tin địa lý chủ chốt nhƣ phân tích khơng gian hình dung Những nỗ lực thúc đẩy cách mạng định lƣợng giới khoa học địa lý đặt móng cho GIS Tác phẩm tiên phong Roger Tomlinson để khởi xƣớng, lên kế hoạch phát triển hệ thống thông tin địa lý Canada dẫn tới GIS giới vào năm 1963 Chính phủ Canada ủy quyền cho Tomlinson tạo kiểm kê quản lý nguồn tài ngun thiên nhiên Ơng hình dung việc sử dụng máy vi tính để kết hợp liệu tài nguyên thiên nhiên từ Thực tốt cơng tác giao đất, khốn rừng, có sách ƣu tiên gia đình sống gần rừng, cạnh rừng nhận đƣợc đất khoán lâu dài, giải tốt vấn đề tranh chấp đất đai giao khốn, có chế độ đãi ngộ hợp lý với cán công chức lao động hộ gia đình giữ cho kinh tế phát triển vƣờn, rừng - Thiết lập mơ hình quản lý cháy rừng sở cộng đồng Rất nhiều hoạt động sử dụng lửa ngƣời dân xung quanh rừng gây khơng thể kiểm sốt đƣợc Vì để làm tốt cơng tác quản lý phịng chống chá rừng cần đƣợc thực kiểm soát nguồn lửa xuất rừng dựa vào cộng đồng dân cƣ sinh sống địa phƣơng có rừng Thực liên kết thành viên cộng đồng , gắn trách nhiệm bảo vệ rừng để không bị cháy rừng với lợi ích mà cộng đồng đạt đƣợc Từ thúc đẩy ý thức ngƣời dân nhƣ nâng cao đƣợc hiệu công tác quản lý cháy rừng 49 Phần V KẾT LUẬN, TỒN TẠI VÀ KIẾN NGHỊ 5.1 Kết luận Từ kết nghiên cứu, đề tài có số kết luận nhƣ sau: - Huyện Tĩnh Gia có tổng diện tích tự nhiên 4.828,67ha, diện tích rừng đất lâm nghiệp 17.723,39ha chiếm 38,6% diện tích huyện Tài nguyên rừng huyện Tĩnh Gia chủ yếu trạng thái rừng Thơng, Keo, Bạch đàn rừng tự nhiên - Tình hình cháy rừng địa bàn huyện Tĩnh Gia phức tạp Diện tích rừng trồng địa bàn nghiên cứu chủ yếu rừng thông nhựa nên khả bắt lửa xảy cháy cao Các vụ cháy thƣờng xảy xã Nguyên Bình, Trúc Lâm, Phú Sơn, Tân Dân vào tháng 5,6 năm - Việc ứng dụng công nghệ viễn thám GIS xây dựng đồ trạng rừng nhƣ đồ giám sát phân tích tài nguyên rừng trở nên phổ biến Việt Nam Đề tài xây dựng đƣợc sở liệu đồ trạng rừng cho đối tƣợng thời điểm trƣớc sau cháy năm 2010,2015, 2016; với độ tin cậy 92,7% Kết nghiên cứu lựa chọn số viễn thám để xác định ngƣỡng cháy rừng khu vực nghiên cứu NDVI, MIRBI, EVI, SAVI, NBR với độ tin cậy từ 68, ÷ 95,65% Trong số NDVI, NBR đƣợc dùng để xây dựng đồ thể phát sớm cháy rừng cho khu vực nghiên cứu Đối với số NDVI ngƣỡng giá trị phát cháy rừng dao động 0,014 ÷ 0,369; số MIRBI 0,023 ÷ 0,512; số EVI 0,037÷ 0,585; số NBR 0,028 ÷ 0,776 số SAVI 0,031 ÷ 0,238 Kết kiểm chứng cho thấy khả quan việc sử dụng số viễn thám để phát sớm cháy rừng cho khu vực huyện Tĩnh Gia, tỉnh Thanh Hóa Nguyên nhân chủ yếu vấn đề cháy rừng khu vực chủ yếu ý thức bất cẩn ngƣời dân với việc dùng lửa hoạt động canh tác mâu thuẫn cá nhân xuất phát từ lợi ích rừng đất lâm nghiệp đốt rừng để phá hoại tài sản Bên cạnh số nguyên nhân khác q trình biến 50 đổi khí hậu xảy khiến nhiệt độ khơng khí khu vực cao gây tình trạng cháy rừng Biện pháp phổ biến đƣợc ứng dụng hiệu thực xử lý vật liệu cháy dƣới tán rừng biện pháp đốt cháy trƣớc có điều khiển, xây dựng tƣờng cản lửa Trong thời gian tới cần xây dựng hệ thống tự động dự báo sớm cháy rừng khu vực nghiên cứu làm tăng tính xác nhƣ giảm thiểu đƣợc thiệt hại có cháy rừng xảy khu vực 5.2 Tồn Đề tài thực huyện Tĩnh Gia, tỉnh Thanh Hóa thu đƣợc kết định song bên cạnh cịn tồn khó khăn định Quá trình khảo sát thực địa diễn thời gian ngắn nên chƣa có nhiều thời gian để kiểm chứng Ảnh viễn thám sử dụng báo cáo nghiên cứu chất lƣợng kém, độ phân giải chƣa cao, sai số xuất trình thành lập đồ Bên cạnh thời tiết diễn biến thất thƣờng nên ảnh hƣởng phần đến độ xác kết Phân bố điểm khảo sát cịn chƣa ảnh hƣởng đến q trình nghiên cứu 5.3 Kiến nghị Cần có hệ thống dự báo sớm cháy rừng khu vực nghiên cứu để đối chiếu lại với kết xây dựng đồ phân cấp chất lƣợng khơng khí để tăng độ tin cậy kết nghiên cứu Cần có tƣ liệu viễn thám với độ xác cao nhằm giảm sai số trình thành lập đồ 51 TÀI LIỆU THAM KHẢO TÀI LIỆU TIẾNG VIỆT: Đỗ Anh Tuấn (2001), Nghiên cứu số nguyên tắc giải pháp quản lý khu Bảo tồn thiên nhiên Pù Mát [10] Nguyễn Hải Hòa, Nguyễn Hữu An (2016), Ứng dụng ảnh viễn thám Landsat GIS xây dựng đồ sinh khối trữ lƣợng cácbon rừng trồng Keo Lai (acacia hybrid) huyện Yên Lập, tỉnh Phú Thọ, T/C Khoa học Công nghệ Lâm nghiệp TÀI LIỆU TIẾNG ANH: [1] Evaluation of Spectral Indices for Assessing Fire Severity in Australian Temperate Forests [2] Oliver Fernando Gomez, 1999: Change Detection of Vegetation Using Landsat Imagery [3] Evaluation of Spectral Indices for Assessing Fire Severity in Australian Temperate Forests [5] Assessment of different spectral indices in the rednear-infrared spectral domain for burned land discrimination [6] Burned Area Detection and Burn Severity Assessment of a Heathland Fire in Belgium Using Airborne Imaging Spectroscopy (APEX) [7] Rubén Ramo, A data approach for global burned area mapping [8] Aitor Bastarrika Mapping Burned Area from Landsat TM/ETM+ data with two – phase algorith [9] Burn Severity estimation from remotelt sensed data: Performance of simulation verus empirical models [11] Burned Area Detection and Burn Severity Assessment of a Heathland Fire in Belgium Using Airborne Imaging Spectroscopy (APEX) [12] Shudong Wang, Muhammad Hasan Ali Baig,Suhong Li, Estimating the area burned by agricultural fires from Landsat Data using the Vegetation Difference Index [13] Lennert Schepers, Birgen Haest Burned Area Detection and Burn Severity Assessment of a Heathland Fire in Belgium Using Airborne Imaging Spectroscopy (APEX) [14] Evaluation of Spectral Indices for Assessing Fire Severity in Australian Temperate Forests [15] Zhuongtin Wu, Vegetation Burn Serverity Mapping using Landsat8 and Worldview - [17] Justin Epting, Evaluation of remotely sensed indicies for assessing burn serverity in interior Alaska using Landsat TM and ETM+ PHỤ LỤC Phụ lục 01: Giá trị số viễn thám trƣớc sau thời điểm cháy STT Kinh độ Vĩ độ Trƣớc cháy Sau cháy NDVI EVI SAVI MIRBI NBR NDVI EVI SAVI MIRBI NBR 105,7649 19,5112 0,450 0,405 0,283 1,409 0,189 0,340 0,140 0,199 1,303 0,415 105,7643 19,5109 0,418 0,297 0,253 1,476 0,299 0,261 0,123 0,147 1,404 0,300 105,7646 19,5109 0,443 0,394 0,276 1,454 0,185 0,332 0,137 0,193 1,348 0,430 105,7641 19,5107 0,404 0,326 0,241 1,454 0,257 0,279 0,117 0,158 1,362 0,312 105,7643 19,5107 0,439 0,337 0,272 1,479 0,274 0,288 0,135 0,162 1,368 0,334 105,7646 19,5106 0,446 0,426 0,278 1,436 0,156 0,350 0,138 0,206 1,332 0,449 105,7635 19,5104 0,423 0,359 0,257 1,439 0,269 0,301 0,126 0,175 1,376 0,330 105,7638 19,5104 0,395 0,352 0,233 1,455 0,194 0,299 0,113 0,171 1,417 0,372 105,7641 19,5104 0,419 0,342 0,254 1,455 0,222 0,292 0,124 0,167 1,371 0,370 10 105,7643 19,5104 0,445 0,414 0,278 1,414 0,159 0,341 0,138 0,200 1,310 0,441 11 105,7629 19,5101 0,400 0,337 0,237 1,533 0,200 0,288 0,115 0,161 1,469 0,418 12 105,7632 19,5101 0,399 0,281 0,237 1,536 0,293 0,246 0,115 0,136 1,437 0,309 12 105,7635 19,5101 0,403 0,302 0,240 1,526 0,321 0,261 0,117 0,148 1,431 0,290 14 105,7638 19,5101 0,408 0,332 0,244 1,500 0,233 0,284 0,119 0,160 1,432 0,381 15 105,7641 19,5101 0,417 0,334 0,252 1,487 0,239 0,286 0,123 0,161 1,429 0,388 16 105,7643 19,5101 0,425 0,336 0,260 1,493 0,207 0,288 0,127 0,162 1,393 0,410 17 105,7629 19,5098 0,411 0,342 0,246 1,525 0,212 0,291 0,120 0,165 1,442 0,414 18 105,7632 19,5098 0,427 0,309 0,261 1,510 0,281 0,269 0,128 0,151 1,398 0,332 19 105,7635 19,5098 0,420 0,310 0,255 1,524 0,262 0,268 0,125 0,150 1,383 0,351 20 105,7638 19,5098 0,424 0,340 0,258 1,517 0,227 0,290 0,127 0,163 1,397 0,391 21 105,7641 19,5098 0,434 0,315 0,267 1,502 0,254 0,273 0,132 0,151 1,354 0,368 22 105,7643 19,5098 0,423 0,315 0,257 1,509 0,239 0,274 0,126 0,152 1,381 0,376 23 105,7652 19,5098 0,421 0,366 0,255 1,443 0,225 0,311 0,125 0,178 1,417 0,408 24 105,7632 19,5096 0,430 0,305 0,263 1,539 0,248 0,267 0,130 0,148 1,377 0,373 25 105,7635 19,5096 0,422 0,308 0,256 1,541 0,238 0,267 0,125 0,147 1,395 0,379 26 105,7638 19,5096 0,445 0,316 0,277 1,520 0,253 0,273 0,137 0,151 1,377 0,376 27 105,7640 19,5096 0,442 0,341 0,275 1,508 0,227 0,293 0,136 0,164 1,347 0,398 28 105,7643 19,5096 0,444 0,331 0,276 1,505 0,246 0,286 0,137 0,160 1,365 0,380 29 105,7646 19,5096 0,429 0,324 0,263 1,482 0,262 0,282 0,129 0,159 1,376 0,353 30 105,7649 19,5096 0,423 0,347 0,257 1,471 0,257 0,298 0,126 0,170 1,407 0,362 31 105,7652 19,5096 0,424 0,298 0,257 1,498 0,333 0,262 0,126 0,146 1,426 0,300 32 105,7692 19,5095 0,432 0,475 0,265 1,477 0,112 0,395 0,130 0,234 1,406 0,529 33 105,7643 19,5093 0,444 0,300 0,276 1,547 0,246 0,265 0,137 0,145 1,359 0,374 34 105,7646 19,5093 0,433 0,281 0,265 1,548 0,284 0,248 0,131 0,136 1,383 0,324 35 105,7649 19,5093 0,423 0,286 0,256 1,545 0,283 0,253 0,125 0,138 1,405 0,336 36 105,7652 19,5093 0,442 0,274 0,272 1,554 0,354 0,244 0,135 0,133 1,415 0,306 37 105,7655 19,5093 0,438 0,331 0,269 1,533 0,248 0,287 0,133 0,159 1,423 0,419 38 105,7692 19,5093 0,412 0,323 0,247 1,489 0,225 0,291 0,121 0,160 1,429 0,398 39 105,7695 19,5093 0,410 0,330 0,247 1,471 0,226 0,296 0,120 0,164 1,439 0,400 40 105,7698 19,5093 0,441 0,342 0,275 1,475 0,234 0,305 0,136 0,170 1,377 0,387 41 105,7701 19,5093 0,490 0,355 0,323 1,480 0,266 0,313 0,164 0,177 1,270 0,357 42 105,7703 19,5093 0,465 0,357 0,305 1,468 0,223 0,315 0,153 0,182 1,231 0,336 43 105,7706 19,5093 0,441 0,473 0,282 1,448 0,081 0,391 0,140 0,239 1,261 0,463 44 105,7646 19,5090 0,416 0,314 0,250 1,530 0,213 0,275 0,122 0,152 1,390 0,406 45 105,7649 19,5090 0,418 0,295 0,252 1,536 0,244 0,260 0,123 0,142 1,408 0,386 46 105,7652 19,5090 0,428 0,284 0,260 1,534 0,283 0,251 0,128 0,137 1,396 0,357 47 105,7655 19,5090 0,420 0,292 0,253 1,538 0,243 0,258 0,124 0,140 1,412 0,389 48 105,7692 19,5090 0,406 0,331 0,244 1,477 0,223 0,296 0,119 0,165 1,437 0,394 49 105,7695 19,5090 0,411 0,311 0,248 1,463 0,259 0,282 0,121 0,156 1,437 0,367 50 105,7698 19,5090 0,421 0,284 0,256 1,467 0,299 0,262 0,125 0,143 1,429 0,329 51 105,7701 19,5090 0,439 0,285 0,272 1,467 0,305 0,263 0,135 0,143 1,385 0,318 52 105,7703 19,5090 0,476 0,314 0,309 1,466 0,283 0,284 0,156 0,158 1,299 0,335 53 105,7706 19,5090 0,436 0,398 0,281 1,436 0,136 0,343 0,139 0,205 1,239 0,362 54 105,7709 19,5090 0,489 0,579 0,328 1,398 0,022 0,451 0,167 0,288 1,149 0,548 55 105,7646 19,5088 0,425 0,361 0,258 1,505 0,168 0,309 0,127 0,174 1,403 0,467 56 105,7649 19,5088 0,415 0,334 0,250 1,495 0,208 0,289 0,122 0,162 1,420 0,421 57 105,7652 19,5087 0,423 0,298 0,257 1,469 0,293 0,262 0,126 0,147 1,393 0,330 58 105,7655 19,5087 0,420 0,292 0,255 1,468 0,304 0,256 0,125 0,143 1,370 0,307 59 105,7692 19,5087 0,428 0,282 0,262 1,507 0,279 0,259 0,129 0,139 1,400 0,343 60 105,7695 19,5087 0,434 0,290 0,267 1,511 0,267 0,265 0,132 0,143 1,418 0,360 61 105,7698 19,5087 0,376 0,319 0,218 1,495 0,157 0,287 0,105 0,158 1,458 0,386 62 105,7701 19,5087 0,375 0,332 0,218 1,493 0,145 0,298 0,105 0,165 1,443 0,402 63 105,7703 19,5087 0,430 0,367 0,265 1,456 0,183 0,322 0,131 0,184 1,369 0,418 64 105,7706 19,5087 0,438 0,335 0,274 1,467 0,230 0,297 0,135 0,167 1,320 0,354 65 105,7709 19,5087 0,446 0,419 0,293 1,406 0,115 0,353 0,146 0,218 1,151 0,384 66 105,7649 19,5085 0,417 0,378 0,254 1,481 0,151 0,321 0,124 0,185 1,394 0,458 67 105,7652 19,5085 0,428 0,328 0,264 1,457 0,258 0,285 0,130 0,164 1,356 0,350 68 105,7655 19,5085 0,420 0,291 0,256 1,496 0,290 0,255 0,126 0,143 1,364 0,310 69 105,7692 19,5085 0,365 0,328 0,208 1,495 0,087 0,293 0,099 0,161 1,473 0,404 70 105,7695 19,5085 0,310 0,317 0,168 1,520 0,027 0,284 0,079 0,157 1,505 0,370 72 105,7701 19,5085 0,333 0,333 0,186 1,467 0,112 0,296 0,088 0,169 1,474 0,332 73 105,7703 19,5085 0,414 0,388 0,251 1,439 0,173 0,336 0,123 0,195 1,392 0,425 74 105,7706 19,5085 0,433 0,342 0,268 1,460 0,226 0,305 0,132 0,171 1,386 0,384 75 105,7709 19,5085 0,440 0,348 0,280 1,469 0,211 0,306 0,139 0,175 1,282 0,335 76 105,7643 19,5082 0,408 0,360 0,245 1,465 0,193 0,307 0,120 0,177 1,409 0,417 77 105,7646 19,5082 0,428 0,300 0,263 1,490 0,269 0,265 0,129 0,148 1,377 0,349 78 105,7649 19,5082 0,446 0,369 0,279 1,478 0,172 0,314 0,138 0,180 1,347 0,451 81 105,7695 19,5082 0,251 0,319 0,129 1,491 0,011 0,283 0,059 0,160 1,556 0,302 82 105,7698 19,5082 0,334 0,287 0,186 1,492 0,191 0,260 0,088 0,146 1,504 0,249 83 105,7701 19,5082 0,393 0,286 0,233 1,479 0,262 0,261 0,112 0,144 1,430 0,297 84 105,7703 19,5082 0,426 0,357 0,260 1,450 0,224 0,314 0,128 0,178 1,390 0,404 85 105,7706 19,5082 0,421 0,396 0,256 1,444 0,163 0,342 0,125 0,197 1,408 0,453 86 105,7643 19,5079 0,421 0,287 0,255 1,514 0,259 0,255 0,125 0,140 1,403 0,364 87 105,7646 19,5079 0,417 0,269 0,251 1,528 0,278 0,240 0,123 0,130 1,416 0,348 88 105,7649 19,5079 0,422 0,329 0,257 1,502 0,192 0,284 0,126 0,159 1,400 0,430 90 105,7695 19,5079 0,392 0,283 0,230 1,511 0,270 0,256 0,111 0,141 1,463 0,264 91 105,7698 19,5079 0,408 0,268 0,244 1,536 0,332 0,246 0,119 0,133 1,422 0,267 92 105,7700 19,5079 0,419 0,268 0,254 1,515 0,286 0,248 0,124 0,133 1,401 0,322 93 105,7703 19,5079 0,425 0,362 0,259 1,503 0,173 0,318 0,127 0,179 1,391 0,448 94 105,7646 19,5077 0,402 0,364 0,239 1,513 0,130 0,311 0,116 0,176 1,439 0,476 95 105,7649 19,5077 0,404 0,392 0,242 1,512 0,085 0,329 0,117 0,188 1,444 0,512 96 105,7692 19,5076 0,365 0,309 0,209 1,498 0,177 0,276 0,100 0,151 1,484 0,362 97 105,7695 19,5076 0,395 0,252 0,232 1,558 0,305 0,234 0,112 0,124 1,454 0,303 98 105,7698 19,5076 0,402 0,248 0,239 1,558 0,271 0,232 0,116 0,122 1,442 0,346 99 105,7700 19,5076 0,409 0,300 0,245 1,546 0,202 0,272 0,119 0,146 1,446 0,421 100 105,7692 19,5074 0,400 0,276 0,235 1,562 0,204 0,252 0,114 0,134 1,472 0,407 101 105,7695 19,5074 0,394 0,275 0,231 1,563 0,208 0,252 0,112 0,134 1,458 0,413 102 105,7698 19,5074 0,404 0,326 0,241 1,544 0,171 0,291 0,117 0,159 1,446 0,453 104 105,7295 19,3839 0,427 0,508 0,271 1,206 0,036 0,283 0,133 0,221 1,241 0,455 110 105,7312 19,3839 0,400 0,530 0,244 1,442 0,002 0,367 0,118 0,239 1,368 0,564 111 105,7315 19,3839 0,433 0,548 0,272 1,454 0,018 0,378 0,134 0,247 1,318 0,582 112 105,7295 19,3836 0,430 0,486 0,271 1,119 0,149 0,251 0,133 0,207 1,281 0,394 116 105,7306 19,3836 0,410 0,420 0,250 1,438 0,157 0,294 0,122 0,188 1,381 0,443 117 105,7309 19,3836 0,418 0,425 0,258 1,445 0,161 0,302 0,126 0,192 1,371 0,452 118 105,7312 19,3836 0,428 0,517 0,268 1,410 0,079 0,357 0,132 0,234 1,363 0,535 119 105,7315 19,3836 0,449 0,614 0,287 1,367 0,021 0,412 0,142 0,279 1,273 0,583 120 105,7298 19,3834 0,397 0,478 0,242 1,133 0,100 0,276 0,117 0,214 1,376 0,432 121 105,7300 19,3834 0,393 0,505 0,238 1,291 0,042 0,313 0,115 0,222 1,365 0,475 122 105,7303 19,3834 0,410 0,462 0,251 1,369 0,112 0,300 0,122 0,203 1,375 0,424 123 105,7306 19,3834 0,423 0,406 0,262 1,398 0,241 0,280 0,128 0,183 1,368 0,369 124 105,7309 19,3834 0,458 0,396 0,296 1,426 0,304 0,280 0,148 0,183 1,314 0,326 125 105,7312 19,3834 0,460 0,440 0,297 1,422 0,228 0,311 0,148 0,201 1,278 0,398 126 105,7315 19,3834 0,452 0,576 0,290 1,372 0,063 0,393 0,144 0,264 1,279 0,552 127 105,7295 19,3831 0,464 0,459 0,302 1,091 0,266 0,248 0,151 0,205 1,257 0,357 128 105,7297 19,3831 0,437 0,466 0,275 1,250 0,215 0,280 0,136 0,204 1,348 0,407 129 105,7300 19,3831 0,436 0,455 0,274 1,331 0,212 0,286 0,135 0,200 1,330 0,396 130 105,7303 19,3831 0,453 0,444 0,291 1,364 0,233 0,287 0,145 0,197 1,283 0,371 131 105,7306 19,3831 0,454 0,420 0,293 1,381 0,274 0,280 0,146 0,188 1,255 0,329 132 105,7309 19,3831 0,472 0,413 0,310 1,413 0,320 0,284 0,156 0,189 1,230 0,281 133 105,7312 19,3831 0,477 0,430 0,315 1,389 0,276 0,300 0,159 0,198 1,206 0,331 134 105,7315 19,3831 0,459 0,494 0,297 1,385 0,159 0,345 0,148 0,225 1,297 0,461 135 105,7297 19,3828 0,467 0,494 0,305 1,217 0,193 0,288 0,153 0,218 1,266 0,422 136 105,7300 19,3828 0,443 0,468 0,281 1,292 0,221 0,287 0,139 0,207 1,321 0,391 137 105,7303 19,3828 0,454 0,449 0,292 1,332 0,238 0,290 0,145 0,201 1,278 0,365 138 105,7306 19,3828 0,468 0,442 0,307 1,340 0,282 0,294 0,154 0,201 1,229 0,323 139 105,7309 19,3828 0,482 0,420 0,321 1,379 0,338 0,287 0,163 0,193 1,228 0,280 140 105,7312 19,3828 0,463 0,409 0,302 1,406 0,296 0,286 0,151 0,189 1,216 0,292 141 105,7315 19,3828 0,455 0,439 0,293 1,393 0,216 0,313 0,146 0,203 1,267 0,391 142 105,7317 19,3828 0,444 0,570 0,283 1,378 0,061 0,390 0,140 0,260 1,287 0,549 143 105,7297 19,3826 0,465 0,520 0,303 1,087 0,157 0,293 0,151 0,234 1,241 0,437 144 105,7300 19,3826 0,460 0,494 0,298 1,166 0,186 0,287 0,149 0,224 1,235 0,418 145 105,7303 19,3826 0,491 0,459 0,332 1,275 0,235 0,292 0,169 0,209 1,105 0,344 146 105,7306 19,3826 0,477 0,469 0,315 1,308 0,263 0,309 0,159 0,213 1,180 0,330 147 105,7309 19,3826 0,461 0,435 0,299 1,326 0,296 0,298 0,149 0,202 1,247 0,300 148 105,7312 19,3826 0,456 0,421 0,295 1,319 0,251 0,295 0,147 0,196 1,236 0,318 149 105,7315 19,3826 0,459 0,475 0,299 1,307 0,174 0,333 0,149 0,221 1,249 0,408 150 105,7317 19,3826 0,461 0,612 0,300 1,299 0,052 0,412 0,150 0,283 1,254 0,548 151 105,7300 19,3823 0,480 0,523 0,320 1,106 0,143 0,294 0,162 0,234 1,145 0,436 152 105,7303 19,3823 0,506 0,483 0,353 1,306 0,235 0,294 0,182 0,217 1,031 0,351 153 105,7306 19,3823 0,491 0,483 0,335 1,290 0,212 0,311 0,170 0,218 1,100 0,363 154 105,7309 19,3823 0,462 0,454 0,302 1,294 0,220 0,306 0,151 0,209 1,159 0,330 155 105,7312 19,3823 0,444 0,420 0,284 1,310 0,251 0,295 0,140 0,197 1,271 0,317 156 105,7303 19,3820 0,494 0,532 0,337 1,170 0,126 0,309 0,172 0,240 1,076 0,445 157 105,7306 19,3820 0,475 0,493 0,316 1,245 0,148 0,306 0,159 0,224 1,132 0,411 158 105,7309 19,3820 0,452 0,503 0,291 1,270 0,161 0,324 0,145 0,228 1,242 0,409 159 105,7312 19,3820 0,438 0,464 0,278 1,288 0,203 0,312 0,137 0,217 1,312 0,385 160 105,7315 19,3820 0,443 0,475 0,283 1,288 0,171 0,329 0,140 0,221 1,261 0,397 161 105,7306 19,3817 0,471 0,522 0,312 1,192 0,084 0,321 0,157 0,240 1,107 0,445 162 105,7309 19,3817 0,481 0,540 0,321 1,233 0,102 0,344 0,162 0,247 1,141 0,444 163 105,7312 19,3817 0,448 0,548 0,288 1,259 0,102 0,349 0,143 0,249 1,234 0,462 164 105,7315 19,3817 0,439 0,572 0,279 1,297 0,078 0,383 0,138 0,265 1,281 0,501 165 105,7309 19,3815 0,466 0,561 0,307 1,242 0,096 0,358 0,153 0,253 1,197 0,464 166 105,7312 19,3815 0,460 0,632 0,299 1,235 0,042 0,402 0,149 0,287 1,228 0,528 Phụ lục 2: Giá trị biến động số viễn thám đánh giá khu vực cháy Thời điểm trƣớc cháy Trung bình ( Mean) Độ lệch chuẩn ( SD) Nhỏ ( Min) Lớn ( Max) Thời điểm sau cháy NDVI SAVI NBR MIRBI EVI 0,457 0,218 0,620 1,569 0,368 0,042 0,031 0,064 0,088 0,193 0,117 0,204 0,586 0,327 0,772 NDVI SAVI NBR MIRBI EVI 0,284 0,139 0,234 1,488 0,233 0,051 0,073 0,038 0,148 0,105 0,064 1,246 0,222 0,087 0,061 0,051 1,044 0,103 1,815 0,568 0,562 0,315 0,701 1,801 0,523