Khai thác dữ liệu vệ tinh phục vụ việc theo dõi biến động thực phủ

Đại Học Quốc Gia Tp.Hồ Chí Minh TRƯỜNG ĐẠI HỌC BÁCH KHOA PHẠM THỊ MAI THY KHAI THÁC DỮ LIỆU VỆ TINH PHỤC VỤ VIỆC THEO DÕI BIẾN ĐỘNG THỰC PHỦ Chuyên ngành: Địa Tin Học Mã số ngành: 2.16.00 LUẬN VĂN THẠC SĨ TP.HỒ CHÍ MINH, tháng 11 năm 2006 MỤC LỤC TÓM TẮT LUẬN VĂN i DANH SÁCH BẢNG BIỂU v DANH SÁCH CÁC HÌNH VẼ vi VIẾT TẮT viii PHẦN I TỔNG QUAN Chương 1.MỞ ĐẦU I LÝ DO CHỌN ĐỀ TÀI II MỤC TIÊU ĐỀ TÀI III Ý NGHĨA KHOA HỌC IV CƠ SỞ TÀI LIỆU Chương 2.TỔNG QUAN VỀ VIỄN THÁM 2.1.TỔNG QUAN VỀ VIỄN THÁM 2.1.1 Các sở viễn thám 2.1.2 Lịch sử viễn thám 2.2 CÁC THỰC TRẠNG VÀ ỨNG DỤNG VIỄN THÁM TẠI VIỆT NAM 10 Chương 3.CƠ SỞ KHOA HỌC TRONG NGHIÊN CỨU BIẾN ĐỘNG THỰC PHỦ 12 3.1 Tổng quan 12 3.2 Thuật toán NDVI (Normalized Difference Vegetation Index) 13 3.3 Thuật toán TVI (Transform Vegetation Index) 14 3.4 Thuật toán RVI (Ratio vegetion index) 14 3.5 Thuật toán EVI (The Enhanced Vegetation Index): 15 3.6 Năm thuật tốn phân tích biến động thực phủ Zhan 16 PHẦN II ỨNG DỤNG DỮ LIỆU MODIS VÀO NGHIÊN CỨU BIẾN ĐỘNG THỰC PHỦ 18 Chương TỔNG QUAN KHU VỰC ĐBSCL VÀ ĐẶC TÍNH DỮ LIỆU MODIS 19 4.1 TỔNG QUAN ĐỒNG BẰNG SÔNG CỬU LONG 19 4.1.1 Vị trí địa lý 19 4.1.2 Điều kiện thủy văn 20 4.1.3 Điều kiện khí hậu 21 Trang ii 4.1.4 Hiện trạng sử dụng đất 21 4.2 TẦM QUAN TRỌNG CỦA NGHIÊN CỨU THỰC PHỦ 22 4.2.1 Khái niệm thực phủ 22 4.2.2 Tầm quan trọng việc lập đồ thực phủ 22 4.3 ĐẶC TÍNH CỦA VỆ TINH TERRA VÀ AQUA MANG BỘ CẢM BIẾN MODIS 23 4.3.1.Vệ tinh TERRA mang cảm biến MODIS 23 4.3.1.1 Đặc điểm vệ tinh TERRA 23 4.3.1.2.Đặc điểm cảm biến MODIS: .25 4.3.1.3 Đặc điểm trạm thu ảnh MODIS Viện Vật Lý Điện Tử 28 4.3.2 Vệ tinh AQUA mang cảm biến MODIS 29 4.3.2.1 Đặc điểm vệ tinh AQUA .29 4.3.2.2 Đặc điểm cảm biến MODIS .32 4.3.3.Ứng dụng ảnh vệ tinh MODIS 32 4.3.3.1 Ứng dụng giới 32 4.3.3.2 Ứng dụng Việt Nam 33 Chương 5.CÁC BƯỚC TIỀN XỬ LÝ 36 5.1 GIẢI PHÁP BIẾN ĐỔI DỮ LIỆU 37 5.1.1 Định chuẩn MODIS: 37 5.1.2 Hiệu chỉnh méo bow-tie: 37 5.1.3 Hiệu chỉnh khí 39 5.1.4 Nắn chỉnh hình học 40 5.2 DỮ LIỆU KIỂM TRA ĐÁNH GIÁ ĐỘ CHÍNH XÁC 41 5.2.1 Nguồn liệu kiểm tra 41 5.2.2 Phương pháp lấy mẫu kiểm tra 42 5.2.2.1 Lấy mẫu kiểm tra từ liệu ảnh vệ tinh Landsat ETM+ 42 5.2.2.2 Lấy mẫu kiểm tra từ mẫu thực địa 43 5.2.3 Công thức đánh giá độ xác 45 Chương 6.PHÂN TÍCH BIẾN ĐỘNG 46 6.1 LOẠI BỎ MÂY 46 Trang iii 6.2 NĂM PHƯƠNG PHÁP PHÂN TÍCH BIẾN ĐỘNG THỰC PHỦ ZHAN 48 6.2.1 Phương pháp phân chia không gian đỏ -hồng ngoại gần (The Red-NIR space partitioning method) 48 6.2.1.1 Phương pháp thực 48 6.2.1.2 Kết thực 52 6.2.1.3 Đánh giá độ xác 55 6.2.2 Phương pháp vector thay đổi không gian đỏ-hồng ngoại gần (The redNIR space change vector method) 56 6.2.2.1 Phương pháp thực 56 6.2.2.2 Kết thực 60 6.2.2.3 Đánh giá độ xác 62 6.2.3 Phương pháp đặt ngưỡng không gian delta (The modified delta space thresholding method) 63 6.2.3.1 Phương pháp thực 63 6.2.3.2 Kết thực 65 6.2.3.3 Đánh giá độ xác 67 6.2.4 Phương pháp cấu trúc thay đổi (Changes in spatial textures) 68 6.2.4.1 Phương pháp thực 68 6.2.4.2 Kết thực 72 6.2.5 Phương pháp nhận biết yếu tố tuyến tính bị thay đổi (Changes in linear features) 74 6.3 ĐÁNH GIÁ 75 Chương 7.KẾT LUẬN VÀ KIẾN NGHỊ 78 7.1 KẾT LUẬN 78 7.2 KIẾN NGHỊ 79 TÀI LIỆU THAM KHẢO 80 TÀI LIỆU TIẾNG VIỆT 80 TÀI LIỆU NƯỚC NGOÀI 80 PHỤ LỤC 83 Trang iv DANH SÁCH BẢNG BIỂU Bảng 4.1 Bảng thông số phổ ảnh vệ tinh MODIS 26 Bảng 5.1 Số điểm kiểm tra để đánh giá độ xác .43 Bảng 6.1 Phân tích biến động thời điểm 46 Bảng 6.2 Bảng tra LUT cho phương pháp phân chia không gian đỏ - hồng ngoại gần .50 Bảng 6.3 Số liệu biến động thực phủ sử dụng phương pháp phân chia không gian đỏ-hồng ngoại gần 54 Bảng 6.4 Độ xác loại biến động phương pháp phân chia không gian đỏ - hồng ngoại gần 55 Bảng 6.5 Bảng tra cho phương pháp vector biến động 58 Bảng 6.6 Số liệu biến động thực phủ sử dụng phương pháp vector thay đổi không gian đỏ-hồng ngoại gần .62 Bảng 6.7 Các sai số phân tích biến động phương pháp phân tích vector biến động khơng gian đỏ-hồng ngoại gần .63 Bảng 6.8: Bảng tra cho phương pháp ngưỡng không gian delta .65 Bảng 6.9 Số liệu biến động thực phủ sử dụng phương pháp ngưỡng không gian delta .67 Bảng 6.10 Sai số biến động phương pháp ngưỡng không gian delta 68 Bảng 6.11 Sai số loại biến động phương pháp phân tích biến động 76 Bảng 6.12 Bảng sai số toàn cục phương pháp 77 Trang v DANH SÁCH CÁC HÌNH VẼ Hình 1.1 Khu vực nghiên cứu-ĐBSCL Hình 2.1 Nguyên lý thu nhận kỹ thuật viễn thám Hình 2.2 Các bước sóng sử dụng viễn thám Hình 2.3 Các loại lượng sensor thu nhận Hình 2.4 Phổ phản xạ đối tượng Hình 2.5 Các pixel giá trị số Hình 3.1 Minh họa giá trị phản xạ hấp thụ thực vật tốt sâu bệnh 14 Hình 4.1 Vệ tinh TERRA 24 Hình 4.2 Quỹ đạo bay vệ tinh TERRA 25 Hình 4.3 Dữ liệu ảnh vệ tinh thu nhận vào ngày 16/11/2002, dãi sọc trắng khơng có liệu phủ vào ngày hôm sau 26 Hình 4.4 Trạm đo ảnh MODIS Viện Vật Lý Điện Tử 28 Hình 4.5 Vệ tinh AQUA cảm biến 31 Hình 4.6 Quỹ đạo bay vệ tinh AQUA 31 Hình 4.7 Màu đỏ khu vực giải đoán khu vực bị cháy rừng có nguy cháy .34 Hình 4.8 Ứng dụng MODIS quan sát diễn biến ngập lụt 35 Hình 5.1 Ảnh hưởng góc qt gây méo bow-tie 38 Hình 5.2 Ảnh vệ tinh MODIS bị méo bow-tie .38 Hình 5.3 Ảnh vệ tinh MODIS trước sau hiệu chỉnh Bow-tie .39 Hình 5.4 Phương pháp nắn chỉnh hình học theo phương pháp Peusodo Affine 40 Hình 5.5 Các điểm tic để nắn chỉnh ảnh có vị trí hệ tọa độ 41 Hình 5.6 Hai ảnh vệ tinh Landsat dùng để lấy điểm kiểm tra 42 Hình 5.7 Kiểu lấy mẫu hệ thống khơng thẳng hàng 43 Hình 5.8 Các điểm kiểm tra để đánh giá độ xác 44 Hình 6.1 Ảnh MODIS trước sau loại bỏ mây 47 Hình 6.2 Mối liên hệ đối tượng với không gian đỏ - hồng ngoại gần, 48 Hình 6.3 Phương pháp chồng lớp tính biến động cơng nghệ GIS 52 Trang vi Hình 6.4 Kết biến động sử dụng phương pháp phân chia không gian đỏ-hồng ngoại gần 53 Hình 6.5 Vector biến động khơng gian đỏ-hồng ngoại gần 57 Hình 6.6 Các bước thực gõ lệnh Envi 60 Hình 6.7 Kết phân tích biến động phương pháp vector thay đổi không gian đỏ-hồng ngoại gần .61 Hình 6.8 Khơng gian delta .64 Hình 6.9 Kết biến động sử dụng phương pháp đặt ngưỡng cho không gian delta 66 Hình 6.10 Cách thể CoV phương pháp cấu trúc thay đổi .70 Hình 6.11 Kết mẫu thực theo phương pháp cấu trúc thay đổi 71 Hình 6.12 Kết sử dụng phương pháp cấu trúc bị thay đổi 73 Hình 6.13 Phân tích yếu tố tuyến tính 74 Hình 6.14 Kết phương pháp Zhan 76 DANH SÁCH CÁC SƠ ĐỒ Sơ đồ 5.1 Quy trình thực trình nghiên cứu 36 Sơ đồ 6.1: Quy trình thực phương pháp phân chia không gian đỏ-hồng ngoại gần .49 Sơ đồ 6.2 Phân loại theo định 51 Sơ đồ 6.3: Quy trình thực phương pháp vector biến động .58 Sơ đồ 6.4 Các bước thực phương pháp cấu trúc thay đổi .69 DANH SÁCH CÁC PHỤ LỤC Phụ lục Các đồ thị phổ mục đích giải đốn mắt mục đích để chọn điểm kiểm tra .83 Phụ lục Các hình ảnh lấy mẫu kiểm tra thực địa 86 Phụ lục Source hàm tính tốn biến động áp dụng phương pháp vector thay đổi 89 Phụ lục Source tính tốn cho phương pháp đặt ngưỡng không gian delta 90 Phụ lục Source tính tốn cho phương pháp cấu trúc thay đổi 91 Trang vii VIẾT TẮT ĐBSCL: Đồng sông Cửu Long EVI: (The Enhanced Vegetation Index) số thực vật nâng cao LUT: Look-Up Table : Bảng Tra Cứu MCST (MODIS Characterization Support Team) Nhóm cung cấp thơng số thuộc tính ảnh vệ tinh MODIS MODIS: (Moderate Resolution Imaging Spectroradiometer) quang phổ kế đo ảnh độ phân giải trung bình MSS: (MultiSpectral Scanner) Bộ quét đa phổ NASA: (National Aeronautics and Space Administration) Cơ quan Hàng không Vũ trụ Hoa Kỳ NDVI: (Normalized Difference Vegetation Index): Chỉ số thực vật sai khác chuẩn RVI: (Ratio Vegetation Index) : Chỉ số tỷ số thực vật TVI: (Transform Vegetation Index) Chỉ số thực vật biến đổi Trang viii PHẦN I TỔNG QUAN HV: Phạm Thị Mai Thy Trang Chương MỞ ĐẦU I LÝ DO CHỌN ĐỀ TÀI Đồng sơng Cửu Long (Hình 1.1) thuộc hệ thống sông MêKong, với Lào, Campuchia , mối quan tâm nhà nghiên cứu Đây vùng rộng lớn qua 13 tỉnh, nơi sản xuất lúa gạo lớn Việt Nam có diện tích trồng rau quả, ăn trái rộng lớn Do vùng nhạy cảm với biến động môi trường tự nhiên (cháy rừng, lũ lụt ) người tác động (phá rừng, thị hóa) Việc quan sát bình diện rộng biến động thực phủ khu vực đồng sông Cửu Long vấn đề cần thiết cho nhà quản lý Hiện nay, ứng dụng ảnh vệ tinh với kỹ thuật phân tích ảnh số hiệu cho phân tích biến động, giúp cho nhà quản lý định nhanh chóng, xác khách quan Trước đây, người hay sử dụng ảnh vệ tinh Landsat, Spot, Radarsat việc quan sát biến động chủ yếu cục Mặc dù ảnh có độ phân giải cao 10m-30m giá thành mắc (1 thời điểm ảnh cho khu vực ĐBSCL tốn hết 3ảnh * 600USD = 1800 USD) trường hợp nghiên cứu vùng rộng lớn Ngoài ảnh vệ tinh có giới hạn tần số thu nhận ảnh: Landsat: 16 ngày, Spot: 26 ngày Đến năm 2001, ảnh MODIS bắt đầu đưa vào sử dụng khắc phục nhược điểm Ảnh MODIS nhanh chóng trở thành cơng cụ hữu ích quan sát biến động có quy mơ tồn cầu Ảnh vệ tinh MODIS thu nhận với tần suất phủ toàn ngày/ảnh cung cấp liệu kịp thời MODIS có 36 kênh phổ khoảng dải phổ rộng: 0.459 đến 14.385 µm thích hợp phân tích nhiều đối tượng: đất, thực phủ, mây, nước, nhiệt độ nước biển, nhiệt độ thềm lục địa, hàm lượng chlrophyl Đặc biệt giá thành ảnh MODIS rẻ (36USD – vùng rộng 1000x1000km) Trong luận văn sử dụng ảnh vệ tinh MODIS có kênh phổ: HV: Phạm Thị Mai Thy Trang Đa số loại biến động sử dụng phương pháp phân chia khơng gian đỏhồng ngoại gần có độ xác cao Ngoài kết biến động nước thành đất nông nghiệp/cây bụi phương pháp vector thay đổi khơng gian đỏ - hồng ngoại gần có độ xác cao so với phương pháp khác Do phân tích biến động thực phủ nên kết hợp phương pháp mang lại kết cao Phương pháp biến đổi cấu trúc áp dụng theo công thức công bố Zhan không kết mong muốn Qua tìm hiểu điều kiện Việt Nam điều kiện địa hình sử dụng đất bị chia nhỏ, khơng đồng (trung bình diện tích đất canh tác hộ dân khoảng 1ha) nên kết tốt Bảng 6.12 Bảng sai số toàn cục phương pháp Các phương pháp phân tích biến động Độ xác tồn cục Phân chia không gian đỏ-hồng ngoại gần 78.8 % Vector thay đổi không gian đỏ-hồng ngoại gần 73.7 % Đặt ngưỡng không gian delta 72.8 % Cấu trúc bị thay đổi Các yếu tố tuyến tính bị thay đổi Không xét Không biến động Phương pháp sử dụng cho ảnh vệ tinh MODIS cho độ xác 7080% việc sử dụng phân tích nhằm mục đích nhận biết tín hiệu ban đầu biến động để cảnh báo sớm Nếu cần độ xác phải cần liệu ảnh vệ tinh có độ xác cao để định xác HV: Phạm Thị Mai Thy Trang 77 Chương KẾT LUẬN VÀ KIẾN NGHỊ 7.1 KẾT LUẬN Qua trình nghiên cứu thực luận văn, tác giả tìm hiểu tính hữu ích ảnh vệ tinh MODIS phân tích biến động thực phủ với phạm vi rộng lớn Luận văn nêu lên phương pháp nghiên cứu biến động thực phủ thấy mạnh phương pháp Zhan Phương pháp truyền thống phân tích biến động thực phủ phân loại thời điểm chồng ghép, phân tích biến động thực phủ Phương pháp Zhan vừa kết hợp với phương pháp truyền thống vừa mở chân trời cách phân tích biến động thực phủ, mang lại kết nhiều tính thuyết phục Bên cạnh luận văn đề xuất số phương pháp phân tích biến động thực phủ Zhan phù hợp với điều kiện ĐBSCL (Bảng 6.2, Bảng 6.5, Bảng 6.8 ) Luận văn cịn đề quy trình xử lý ảnh vệ tinh MODIS nhà quản lý dễ dàng ứng dụng phân tích biến động thực phủ (Sơ đồ 5.1, Sơ đồ 6.1, Sơ đồ 6.3) Ngoài luận văn lập đồ phân tích biến động thực phủ khu vực ĐBSCL (Hình 6.4, Hình 6.7, Hình 6.9, Hình 6.12) Độ xác phương pháp tính tốn qua cho thấy phương pháp “Phương pháp phân chia không gian đỏ - hồng ngoại gần” có độ xác cao nhất, đáng tin cậy Phương pháp “Phương pháp cấu trúc thay đổi” không thành công trường hợp điều kiện Việt Nam Cuối luận văn xây dựng số module phần mềm ENVI ngơn ngữ IDL để hỗ trợ việc tính tốn nhanh chóng thực phủ, sử dụng thuật tốn nêu HV: Phạm Thị Mai Thy Trang 78 7.2 KIẾN NGHỊ Do giới hạn điều kiện thu thập nên liệu vệ tinh MODIS có tương đối thời gian khoảng cách hai ảnh ngắn nên khơng thấy nhiều biến động Nếu có điều kiện thu thập liệu vệ tinh MODIS nhiều thời điểm kết luận văn mang lại nhiều biến động Các thông số bảng tra kinh nghiệm tác giả việc kết hợp với ảnh vệ tinh Landsat độ phân giải 15m thành lập nên Do thơng số hồn thiện thêm cho xác kết hợp phân tích thêm với ảnh vệ tinh tốt độ phân giải, độ tương phản kênh ví dụ Spot, Aster Các điểm kiểm tra độ xác từ thực địa ít, khơng trải đều, chủ yếu điểm lấy ảnh vệ tinh có độ phân giải cao Các điểm thực địa tập trung chủ yếu huyện Chợ Mới tỉnh An Giang, ven biển tỉnh Trà Vinh điểm dọc theo tuyến đường tỉnh Long An, Tiền Giang, Đồng Tháp Nếu tác giả có điều kiện thực địa tồn ĐBSCL luận văn mang tính thuyết phục Hiện luận văn tìm thuật tốn thích hợp cho phân tích biến động thực phủ quy trình xử lý rút ngắn lại cho người xử lý chương trình tự động phần mềm ENVI Tuy nhiên người sử dụng phải biết chút phần mềm ENVI Do hướng nghiên cứu tác giả lập chương trình quản lý chạy độc lập để người dùng khơng chun sử dụng dễ dàng Do liệu vệ tinh MODIS thu nhận phạm vi rộng lớn nên không tránh khỏi ảnh hưởng mây Kỹ thuật xử lý mây hướng nghiên cứu sau tác giả, cách sử dụng nhiều ảnh thời gian gần áp dụng thuật toán loại bỏ mây HV: Phạm Thị Mai Thy Trang 79 TÀI LIỆU THAM KHẢO TÀI LIỆU TIẾNG VIỆT • Bộ kế hoạch đầu tư, http://www.mpi.gov.vn/Bộ Tài Nguyên Môi Trường, http://www.rsc.gov.vn/ • Hồng Viết Giao, Nguyễn Hùng Sơn, Hà Vĩnh Tân, Doãn Hà Phong, Nguyễn Quang Huy, Nguyễn Tiến Dũng, 2004 Báo cáo: Xây dựng sở liệu vệ tinh TERRA-MODIS, phục vụ phòng chống ngập lụt theo dõi môi trường Đồng Bằng Sông Cửu Long Viện Vật Lý Điện Tử • Lê Văn Trung, 2002 báo cáo ‘Nghiên cứu hệ thống tích hợp viễn thám GIS xây dựng cập nhật liệu không gian (phục vụ thành lập sở liệu cho SAGOGIS) Sở Khoa Học Công Nghệ Môi Trường Tp.Hồ Chí Minh • Nguyễn Ngọc Thạch, 1997 Viễn Thám Trong Nghiên Cứu Tài Nguyên Môi Trường NXB Khoa Học & Kỹ Thuật Hà Nội • Tổng cục thống kê Việt Nam, http://www.gso.gov.vn/ TÀI LIỆU NƯỚC NGỒI • Alan Strahler, Doug Muchoney, Jordan Borak, Mark Friedl, Sucharita Gopal, Eric Lambin and Aaron Moody, MODIS Land Cover and Land cover Change – Algorithm Theoretical Basis Document • Alfredo Heute, Chris Justice, Wim Van Leeuwen, 1999 MODIS Vegetation Index, Algorithm Theoretical Basis Document • Bo_Cai Gao, Ping Yang, Wei Han, Rong-Rong and Warren J.Wiscombe, Aug 2002 An Algorithm Using Visible and 1.38-µm Channels to Retrieve Cirrus Cloud Reflectances From Aircraft and Satellite Data IEEE Transactions on GeoScience and Remote Sensing, Vol 40, No.8 • Canada Centre of Remote Sensing, 1998, Fundamentals of Remote Sensing http://www.ccrs.nrcan.gc.ca/ HV: Phạm Thị Mai Thy Trang 80 • Christopher O.Justice, Eric Vermote, John R.G.Townshend, Ruth Defrices, David P.Roy, Dorothy K.Hall, Vincent V.Salomonson, Jeffrey L.Privette, George Riggs, Alan Strahler, July 1998 The Moderate Resolution Imaging Spectroradiometer (MODIS): Land remote sensing for Global Change Research IEEE Transactions on Geoscience and Remote Sensing, Vol.36, No.4 • Deering, D W., Rouse, J W., Haas, R H., and Schell, J A., 1975 Measuring "for age production" of grazing Units from Landsat MSS data Proceedings of the 10th International Symposium on Remote Sensing of Environment, II, pp 1169-1178 • Japan Association of Remote Sensing (JARS), 1999 Remote Sensing Notes • John A.Richards, 1994 Remote Sensing Digital Image Analysis Springer Verlag • John Weier David Herring, Measuring Vegetation, the Earth Observatory, http://earthobservatory.nasa.gov/Library/MeasuringVegetation/ • M.Hansen, R.DeFries, C.Dimiceli, C.Huang, R.Sohlberg, X.Zhan, J.Townshend, 1998 Red and Infrared space partitioning for detecting land cover change IEEE Transactions on Geoscience and Remote Sensing • MODIS Website, http://modis.gsfc.nasa.gov/ • NOAA Coastal Services Center, http://www.csc.noaa.gov/crs/lca/ • Roger M Mc Coy, 2005 Field Methods The Guilford Press • Rouse, J W., R H Haas, J A Schell, and D W Deering, 1973 Monitoring vegetation systems in the Great Plains with ERTS Third ERTS Symposium, NASA SP-351 I, 309-317 • Russell G.Congalton, Kass Green Assessing the Accuracy of Remotely Sensed Data: Principles and Practices Lewis Publishers • W.Paul.Menzen, Suzanne W.Seemann, Jun Li, Liam E.Gumley, 2002 MODIS Asmospheric Profile Retrieval- Algorithm Theoretical Basis Document HV: Phạm Thị Mai Thy Trang 81 • X.Zhan, C.Huang, J.Townshend, R.DeFries, M.Hansen, C.Dimiceli, R.Sohlberg, J.Hewson-Scardelletti, A.Tompkins, 1998 Land cover change detection with change vector in the red and near-infrared reflectance space IEEE Transactions on Geoscience and Remote Sensing • X.Zhan, R.A Sohlberg, 2002 Detection of land cover changes using MODIS 250m data Remote Sensing for Environment 83 (2002) page 336- 360 • X.Zhan, R.DeFries, J.Townshend, C.Dimiceli, M.Hansen, C.Huang, R.Sohlberg, 2000 The 250m global land cover change product from the Moderate Resolution Imaging Spectroradiometer of the NASA’s earth observing system International Journal of Remote Sensing – 1433-1460 No & 7, Vol.21 HV: Phạm Thị Mai Thy Trang 82 PHỤ LỤC Phụ lục Các đồ thị phổ mục đích giải đốn mắt mục đích để chọn điểm kiểm tra 80 Thực vật 70 Đất trống 60 Giá trị phản xạ (%) Nước sông Nước sông đục 50 40 30 20 10 0.4 0.6 0.8 1.0 1.2 1.4 1.6 Bước sóng (μm) 1.8 2.0 2.2 2.4 2.6 đường thẳng đứng bước sóng kênh 2, 3, (0.56, 0.66, 0.83 μm) ảnh vệ tinh Landsat ETM+ Mẫu đồ thị phổ nông nghiệp HV: Phạm Thị Mai Thy Trang 83 Mẫu đồ thị phổ đất trống Mẫu đồ thị phổ rừng HV: Phạm Thị Mai Thy Trang 84 Mẫu đồ thị phổ nước HV: Phạm Thị Mai Thy Trang 85 Phụ lục Các hình ảnh lấy mẫu kiểm tra thực địa Ruộng lúa bị ngập nước vào thời điểm tháng 11 đất trồng lúa vào tháng năm HV: Phạm Thị Mai Thy Trang 86 Đất trồng lúa vào tháng 11 trở thành đất trống (đã gặt) vào tháng năm HV: Phạm Thị Mai Thy Trang 87 Đất ni tơm q trình bị bỏ trống HV: Phạm Thị Mai Thy Trang 88 Phụ lục Source hàm tính tốn biến động áp dụng phương pháp vector thay đổi HV: Phạm Thị Mai Thy Trang 89 Phụ lục Source tính tốn cho phương pháp đặt ngưỡng không gian delta HV: Phạm Thị Mai Thy Trang 90 Phụ lục Source tính tốn cho phương pháp cấu trúc thay đổi HV: Phạm Thị Mai Thy Trang 91 ... với độ phân giải 250m tốt cho biến động thực phủ Với yêu cầu thực tiễn đề tài nghiên cứu có tên: ? ?Khai thác liệu ảnh vệ tinh MODIS việc theo dõi biến động thực phủ? ?? Hình 1.1 Khu vực nghiên cứu-ĐBSCL... đồ thực phủ biến động khu vực nghiên cứu độ xác phương pháp thực III Ý NGHĨA KHOA HỌC Việc nghiên cứu ứng dụng ảnh vệ tinh MODIS vấn đề lớn mang tính tồn cầu việc quan sát biến động thực phủ. .. Nhận định thực phủ xác định rõ ranh giới từ hoạt động kiểm tra (khám phá biến động thực phủ) thể Những biến động thực phủ phép đo trực tiếp số lượng môi trường sống bị hay tăng lên Dữ liệu mơ tả