khảo sát tình hình lũ lụt ở đồng bằng sông cửu long giai đoạn năm 2013 bằng ảnh modis

60 797 0
khảo sát tình hình lũ lụt ở đồng bằng sông cửu long giai đoạn năm 2013 bằng ảnh modis

Đang tải... (xem toàn văn)

Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống

Thông tin tài liệu

TRƢỜNG ĐẠI HỌC CẦN THƠ KHOA MÔI TRƢỜNG & TÀI NGUYÊN THIÊN NHIÊN BỘ MÔN TÀI NGUYÊN ĐẤT ĐAI - DANH XÀ RÂY KHẢO SÁT TÌNH HÌNH LŨ LỤT Ở ĐỒNG BẰNG SÔNG CỬU LONG GIAI ĐOẠN NĂM 2013 BẰNG ẢNH MODIS LUẬN VĂN KỸ SƢ NGÀNH QUẢN LÝ ĐẤT ĐAI Cần Thơ, 2014 TRƢỜNG ĐẠI HỌC CẦN THƠ KHOA MÔI TRƢỜNG & TÀI NGUYÊN THIÊN NHIÊN BỘ MÔN TÀI NGUYÊN ĐẤT ĐAI - DANH XÀ RÂY KHẢO SÁT TÌNH HÌNH LŨ LỤT Ở ĐỒNG BẰNG SÔNG CỬU LONG GIAI ĐOẠN NĂM 2013 BẰNG ẢNH MODIS LUẬN VĂN KỸ SƢ NGÀNH QUẢN LÝ ĐẤT ĐAI Mã ngành: 52850103 CÁN BỘ HƢỚNG DẪN SINH VIÊN THỰC HIỆN Ths: HUỲNH THỊ THU HƢƠNG DANH XÀ RÂY MSSV: 4115077 Lớp Quản Lý Đất Đai K37A2 Cần Thơ, 2014 TRƢỜNG ĐẠI HỌC CẦN THƠ KHOA MÔI TRƢỜNG VÀ TÀI NGUYÊN THIÊN NHIÊN BỘ MÔN TÀI NGUYÊN ĐẤT ĐAI NHẬN XÉT CỦA BỘ MÔN Chứng nhận chấp thuận báo cáo luận văn tốt nghiệp chuyên ngành Quản lý đất đai với đề tài: “Giải đoán tình hình lũ đồng sông Cửu Long ảnh MODIS giai đoạn năm 2013 ” Do sinh viên: Danh Xà Rây MSSV: 4115077 Lớp: Quản lý đất đai K37A2 thuộc Bộ môn Tài nguyên đất đai, Khoa Môi Trƣờng Tài Nguyên Thiên Nhiên - Trƣờng Đại Học Cần Thơ từ tháng 07/2014 đến tháng 12/2014 Nhận xét: ………………………………………………………………………………………… ………………………………………………………………………………………… ………………………………………………………………………………………… ………………………………………………………………………………………… ………………………………………………………………………………………… ………………………………………………………………………………………… ………………………………………………………………………………………… Cần Thơ, ngày ……tháng………năm 2014 Trƣởng môn i TRƢỜNG ĐẠI HỌC CẦN THƠ KHOA MÔI TRƢỜNG VÀ TÀI NGUYÊN THIÊN NHIÊN BỘ MÔN TÀI NGUYÊN ĐẤT ĐAI NHẬN XÉT CỦA CÁN BỘ HƢỚNG DẪN Chứng nhận chấp thuận báo cáo luận văn tốt nghiệp chuyên ngành Quản lý đất đai với đề tài: “Giải đoán tình hình lũ đồng sông Cửu Long ảnh MODIS giai đoạn năm 2013” Do sinh viên: Danh Xà Rây MSSV: 4115077 Lớp: Quản lý đất đai K37A2, thuộc Bộ môn Tài nguyên đất đai, Khoa Môi trƣờng Tài Nguyên Thiên Nhiên - Trƣờng Đại Học Cần Thơ thực từ tháng 07/2014 đến tháng 012/2014 Nhận xét: …………………………………………………………………………………… … ………………………………………………………………………………………… ………………………………………………………………………………………… ………………………………………………………………………………………… ………………………………………………………………………………………… ………………………………………………………………………………………… ………………………………………………………………………………………… ………………………………………………………………………………………… Cần Thơ, ngày…… tháng… năm 2014 Cán hƣớng dẫn Huỳnh Thị Thu Hƣơng ii TRƢỜNG ĐẠI HỌC CẦN THƠ KHOA MÔI TRƢỜNG VÀ TÀI NGUYÊN THIÊN NHIÊN BỘ MÔN TÀI NGUYÊN ĐẤT ĐAI NHẬN XÉT CỦA HỘI ĐỒNG BÁO CÁO Hội đồng chấm luận văn tốt nghiệp chấp thuận báo cáo luận văn tốt nghiệp chuyên ngành Quản lý đất đai với đề tài: “Giải đoán tình hình lũ đồng sông Cửu Long ảnh MODIS giai đoạn năm 2013” Do sinh viên: Danh Xà Rây MSSV: 4115077 Lớp: Quản lý đất đai K37A2, thuộc Bộ môn Tài nguyên đất đai, Khoa Môi trƣờng Tài Nguyên Thiên Nhiên - Trƣờng Đại Học Cần Thơ thực báo cáo trƣớc Hội đồng Ngày……….tháng………năm 2014 Báo cáo luận văn tốt nghiệp đƣợc Hội đồng đánh giá mức:…………………… Ý kiến hội đồng: ………………………………………………………………………………………… ………………………………………………………………………………………… ………………………………………………………………………………………… ………………………………………………………………………………………… ………………………………………………………………………………………… ………………………………………………………………………………………… Cần Thơ, ngày……tháng……năm 2014 Chủ tịch Hội đồng iii LỊCH SỬ CÁ NHÂN Họ tên: Danh Xà Rây Ngày sinh: 11/09/1993 Nơi sinh: huyện Long Mỹ, tỉnh Hậu Giang Quê quán: Ấp 4, xã Xà Phiên, huyện Long Mỹ, tỉnh Hậu Giang Họ tên cha: Danh Hiêl Họ tên mẹ: Thị Thời Tốt nghiệp Trung học phổ thông năm 2011, trƣờng phổ thông Dân Tộc Nội Trú Hậu Giang Vào học Trƣờng Đại học Cần Thơ năm 2011, ngành học Quản lý đất đai Tốt nghiệp Kỹ sƣ chuyên ngành Quản lý đất đai năm 2014 iv LỜI CAM ĐOAN Tôi cam đoan đề tài thực hiện, số liệu thu thập phân tích đề tài trung thực, không trùng với đề tài nghiên cứu khoa học Cần Thơ, ngày……tháng…….năm 2014 Danh Xà Rây v LỜI CẢM TẠ Trong suốt thời gian học tập, em vinh dự nhận đƣợc dạy dỗ rèn luyện kiến thức trƣờng Đại học Cần Thơ, quan tâm dạy dỗ tận tình thầy cô Bộ môn Tài nguyên đất đai, Khoa Môi trƣờng Tài nguyên thiên nhiên Đây tảng vô quan trọng, hành trang tri thức giúp em vững tin trình công tác tƣơng lai Em xin chân thành gửi lời cám ơn đến: Quý thầy cô Bộ môn Tài nguyên đất đai, khoa Môi trƣờng Tài nguyên thiên nhiên, trƣờng Đại học Cần Thơ trực tiếp giảng dạy hƣớng dẫn em suốt thời gian vừa qua Đặc biệt em xin gửi lời cám ơn sâu sắc đến cô Phan Kiều Diễm cô Nguyễn Thị Song Bình cố vấn học tập quan tâm giúp đỡ em thời gian vừa qua Em xin chân thành cám ơn cô Huỳnh Thị Thu Hƣơng trực tiếp hƣớng dẫn em thực đề tài, cô nhiệt tình hƣớng dẫn, đóng góp ý kiến quý báo tạo điều kiện tốt để em hoàn thành đề tài luận văn tốt nghiệp Cảm ơn ba mẹ dạy dỗ chăm sóc con, ba mẹ niềm động viên vững để vững tin bƣớc vào đời Một lần xin cám ơn bạn lớp Quản lý đất đai K37 động viên giúp đỡ trình thực luận văn tốt nghiệp Danh Xà Rây vi TÓM LƢỢC Lũ tƣợng tự nhiên xảy hàng năm Những năm gần ảnh hƣởng Biến đổi khí hậu với ấm lên toàn cầu làm cho lũ ngày trở nên phức tạp gây nên thiệt hại nghiêm trọng đến tài sản tính mạng ngƣời, sụt giảm sản lƣợng nông nghiệp diễn tiến lũ tăng nhanh Với mục tiêu: - Theo dõi trình diễn tiến lũ theo thời gian đồng sông Cửu Long năm 2013 - Xây dựng đồ trạng ngập lũ đồ phân vùng chế độ ngập lũ đồng sông Cửu Long năm 2013 Đề tài nghiên cứu“Giải đoán tình hình lũ lụt đồng sông Cửu Long giai đoạn năm 2013 ảnh MODIS” đƣợc thực với phƣơng pháp nghiên cứu sử dụng ảnh vệ tinh MODIS việc tính toán số EVI, LSWI DVEL Cùng việc thực trình phân loại kết hợp số EVI, LSWI, DVEL xác định vùng ngập lũ không ngập Qua phân loại không kiểm soát chuỗi ảnh DVEL nghiên cứu tiến hành đƣa phân tích, đánh giá cụ thể không gian nhƣ thời gian ngập lũ ĐBSCL năm 2013 Kết đạt đƣợc thông qua việc theo dõi trình diễn tiến lũ theo thời gian qua thành lập đồ ngập lũ giai đoạn đồ phân vùng chế độ ngập lũ ĐBSCL năm 2013 Kết Nghiên cứu cho thấy mối tƣơng quan ảnh theo không gian, thời gian với thay đổi lũ có độ xác tƣơng đối cao dựa so sánh tƣơng quan diện tích ngập mực nƣớc thực đo trạm thủy văn Tân Châu Châu Đốc cho thấy sử dụng ảnh vệ tinh đa thời gian MODIS - MOD09A1 để theo dõi diễn tiến lũ khu vực Đồng Bằng Sông Cửu Long vii MỤC LỤC NHẬN XÉT CỦA BỘ MÔN I NHẬN XÉT CỦA CÁN BỘ HƢỚNG DẪN II NHẬN XÉT CỦA HỘI ĐỒNG BÁO CÁO III LỊCH SỬ CÁ NHÂN IV LỜI CAM ĐOAN V LỜI CẢM TẠ VI TÓM LƢỢC VII MỤC LỤC VIII DANH MỤC TỪ VIẾT TẮT X DANH SÁCH BẢNG XI DANH SÁCH HÌNH XII MỞ ĐẦU CHƢƠNG I: LƢỢC KHẢO TÀI LIỆU 1.1 Đặc điểm vùng nghiên cứu 1.1.1 Vị trí địa lý 1.1.2 Địa hình, địa vật .2 1.1.3 Khí hậu thủy văn .2 1.1.4 Đất đai 1.1.5 Hiện trạng sử dụng đất đồng sông Cửu Long 1.2 Hiện trạng ngập lũ ĐBSCL giai đoạn 2000 đến 2012 1.2.1 Khái niệm lũ 1.2.2 iện trạng ngập lũ đồng Sông Cửu Long 1.2.3 Nguyên nhân gây lũ đồng sông Cửu Long .6 1.3 Giới thiệu tổng quan viễn thám .7 1.3.1 Khái niệm 1.3.2 Phân loại viễn thám .7 1.3.3 Nguyên lý hoạt động viii Trên sở kết phân loại trạng ngập lũ, đề tài phân tách đƣợc nhóm đối tƣợng là: 1) Ngập 2) Không ngập Hình 3.9 Sự thay đổi giá trị DVEL phân vùng chế độ ngập lũ ĐBSCL - Vùng 1: ngập quanh năm hay thủy sản chủ yếu tỉnh ven biển nhƣ: Kiên Giang (Kiên Lƣơng, An Biên), Cà Mau, Bạc liêu, Sóc Trăng (Vĩnh Châu), Trà Vinh (Trà Ôn, Duyên Hải, Cầu Ngang), Bến Tre (Ba Tri, Bình Đại) với diện tích khoảng 669.477 - Vùng 2: ngập từ 12/7 - 17/11 vùng có thời gian ngập lũ dài nhất, thời điểm bắt đầu ngập năm với diện tích khoảng 444.794 chủ yếu tỉnh đầu nguồn nhƣ Long An (Mộc hóa, Tân Thạnh), Đồng Tháp (Hông Ngự, Tam Nông, Thanh Bình), An Giang (An Phú, Châu Phú), Kiên Giang (Kiên Lƣơng, Vĩnh Thuận) phần nhỏ tỉnh quen biển - Vùng 3: ngập từ 5/8 - 17/11 vùng có thời gian ngập dài chiếm diện tích khoảng 424.981 chủ yếu tập chung tỉnh nhƣ Long An (Tân Hƣng, Vĩnh Hƣng), Đông Tháp (Lấp Vò, Tam Nông), An Giang (Tịnh Biên, Châu Phú), Cà Mau (Trần Văn Thời) nằm rải rác số tỉnh khác Kiên Giang, Cần Thơ, Sóc Trăng - Vùng 4: ngập từ 6/9 - 17/11 vùng ngập tƣơng đối trễ chiếm diện tích lớn khoảng 778.652 tập trung tỉnh Long An (Thủ Thừa, Tân Hƣng), Tiền Giang (Cái Bè, Cai Lậy), Cần Thơ (Vĩnh Thạnh, Thốt Nốt, Cờ Đỏ, Ô Môn), Hậu Giang (Vị Thanh, Vị Thủy) - Vùng 5: ngập từ 22/9 - 17/11 vùng có thời gian ngập trễ nhƣng chiếm diện tích cao khoảng 636.913 Qua hình 3.9 cho thấy diễn tiến lũ tỉnh đầu nguồn nhƣ Long An, Đồng Tháp, An Giang Sau lan rộng vùng lân cận nhƣ 32 Tiền Giang, Kiên Giang (Kiên Lƣơng, Hòn Đất, Giồng Riềng, An Biên), Vĩnh Long (Bình Tân, Tam Bình, Vũng Liêm), Hậu Giang (Vị Thanh, Vị Thủy, Long Mỹ), Sóc Trăng (Mỹ Xuyên, Mỹ Tú, Thạnh Trị) - Vùng 6: không ngập từ 4/7 - 17/11 chủ yếu tập trung rìa ranh giới tỉnh với Qua (hình 3.9) cho thấy vùng tập trung chủ yếu tỉnh Đồng Tháp (Tháp Mƣời), An Giang (Tân Phú, Thoại Sơn, Long Xuyên), Sóc Trăng (Kế Sách, Thạnh Trị), Bạc Liêu (Phƣớc Long) - Vùng 7: không ngập quanh năm Do vùng không ngập nên có nhiều kiểu sử dụng đất khác nhƣ canh tác lúa vụ, lâu năm, chuyên màu…, chiếm diện tích lớn khoảng 487.017 chủ yếu tỉnh An Giang (Tri Tôn), Trà Vinh, Bến Tre, Đồng Tháp Kết phân loại trạng phân loại vùng ngập lũ không kiểm soát có đánh giá tƣơng đối xác trạng diễn biến lũ khu vực ĐBSCL có phân biệt rõ ràng điểm ảnh liên quan đến nƣớc, nhằm cung cấp nhìn xác thay đổi không gian ngập lũ theo thời gian năm 2013 Diện tích ngập lũ thay đổi theo thi gian ĐBSCL đƣợc trình bày bảng 3.4 Với số liệu diện tích cụ thể giúp việc xác định giai đoạn lũ trở nên xác dễ dàng Bảng Bảng Diện tích phân vùng ngập lũ ĐBSCL năm 2013 STT Đối tƣợng phân loại Diện tích (ha) Thủy sản hay ngập quanh năm 669.477,47 Ngập từ 12/7 – 17/11 444.794,24 Ngập từ 5/8 – 17/11 424.981,24 Ngập từ 6/9 – 17/11 778.652,87 Ngập từ 22/9 – 17/11 636.913,75 Không ngập từ 4/7 – 17/11 573.975,83 Không ngập hoàn toàn 487.017,60 Tổng 401.581,30 33 Lũ biến đổi ngày không gian ngập lũ thay đổi liên tục theo thời gian, việc theo dõi lũ thông qua số liệu mực nƣớc đơn xác định thay đổi diện tích mực mặt nƣớc theo thời gian Không xác định đƣơc thời điểm ngập lũ đạt mức cực đại, xác định đƣợc cụ thể thời điểm đỉnh lũ đạt mức cực đại dựa vào (bảng 3.3) Tuy nhiên, với việc sử dụng liệu ảnh vệ tinh MODIS cho thấy thay đổi diện tích mặt nƣớc không gian lẫn thời gian, cung cấp nhìn rõ ràng trực quan phân loại đối tƣợng ngập nƣớc đƣợc thể hoàn toàn ảnh Hình 3.10 Diện tích phân vùng chế độ ngập lũ ĐBSCL năm 2013 Hình 3.9 thể không gian ngập cực đại thời điểm định năm 2013 ĐBSCL Thời gian ngập cực đại tỉnh dao động khoảng từ đầu tháng đến tháng 11 Thời điểm bắt đầu ngập lũ vào ngày 12/7 với diện tích ngập khoảng 444.794 nhƣng tăng lên nhanh chóng vào ngày 6/9 với diện tích khoảng 778.652 ha, đến khoảng tháng 11 lũ bắt đầu hạ xuống Theo Nguyễn Thị Bích Liên (2012), thời gian ngập lũ cuối tháng 7, đạt cực đại vào tháng 10 diện tích mặt nƣớc giảm đáng kể vào tháng 12 năm So với năm 2012 thời gian bắt đầu ngập sớm hơn, thời điểm ngập cực đại xảy chậm vào cuối tháng 10 nhƣng kết thúc lũ tƣơng đối giống với năm 2012 Nhƣ vậy, dựa vào đồ trạng sử dụng đất ĐBSCL năm 2013 (Nguyễn Thị Kim Yến, 2013) kết hợp (chồng lấp, gom nhóm phần mềm ENVI) đồ phân vùng chế độ ngập lũ ĐBSCL năm 2013 Kết thành lập đồ phân vùng ngập lũ ĐBSCL năm 2013 (hình 3.11) 34 Hình 3.11 Bản đồ phân vùng ngập lũ ĐBSCL năm 2013 - Vùng 1: ngập quanh năm hay thủy sản chiếm diện tích 669.477 chủ yếu phân bố vùng ven biển nhƣ Bến Tre (Bình Đại, Thạnh Phú), Trà Vinh (Duyên Hải), Sóc Trăng (Vĩnh Châu), Kiên Giang (Kiên Lƣơng, An Minh) phần lớn thành phố Cà Mau vùng có kiểu sử dụng đất vụ lúa mùa kết hợp luân canh nuôi tôm - Vùng 2: ngập lũ chiếm phần lớn diện tích khoảng 2285.023 ha, tập trung nhiều kiểu sử dụng đất khác nhƣ canh tác lúa vụ, lâu năm, chuyên màu…, phân bố tập chung chủ yếu tỉnh Long An, Đồng Tháp, Kiên Giang, Cần Thơ, Hậu Giang, Bến Tre, Sóc Trăng - Vùng 3: không ngập hoàn toàn tập trung chủ yếu tỉnh nhƣ An Giang (Tân Phú, Tịnh Biên, Châu Đốc), Long An (Châu Thành, Bến Lức, Cần Giuộc), Tiền Giang (Cái Bè, Cai 35 Lậy, Châu Thành, Mỹ Tho, Gò Công), Cần Thơ (Phong Điền), Hậu Giang (Châu Thành, Phụng hiệp), Sóc Trăng (Thạnh Trị, Cù Lao Dung), Bến tre Do vùng không ngập nên có nhiều kiểu sử dụng đất khác nhƣ canh tác lúa vụ, lâu năm, chuyên màu, rừng với diện tích khoảng 1210.888 3.2 Kết đánh giá độ tƣơng quan với d liệu quan trắc Để đánh giá đƣợc khả sử dụng ảnh vệ tinh MODIS việc giám sát lũ ĐBSCL năm 2013 Thì việc kiểm tra so sánh lại kết sau phân loại trở nên cần thiết, nghiên cứu tiến hành so sánh dựa số liệu mực nƣớc thủy văn ghi nhận đƣợc trạm quan trắc với diện tích ngập đƣợc giải đoán từ ảnh vệ tinh MODIS Số liệu mực nƣớc thủy văn trạm đo Tân Châu Châu Đốc đƣợc ghi nhận lại ngày phản ánh xác mức độ lũ năm Nghiên cứu tiến hành so sánh mối tƣơng quan diện tích ngập lũ ĐBSCL đƣợc giải đoán từ ảnh MOD09A1 với số liệu mực nƣớc thủy văn ghi nhận trạm quan trắc Tân Châu Châu Đốc (Phụ lục 3, 4) Số liệu quan trắc sử dụng so sánh nghiên cứu đƣợc thu thập khoảng thời gian bắt đầu mùa lũ từ đầu tháng 6/2013, cụ thể từ ngày 1/6/2013 - 31/10/2013 Kết tƣơng quan diện tích ngập lũ ĐBSCL với số liệu mực nƣớc thủy văn đƣợc thể hình (hình 3.12, hình 3.13) nhƣ sau: y = 359.12x + 893.9 R² = 0.6799 3000 2500 Diện tích ngập Đơn vị (NGhìn ha) Diện tích ngập Đơn vị (NGhìn ha) 3000 2000 1500 1000 500 0 y = 298x + 892.54 R² = 0.5821 2500 2000 1500 1000 500 Mực nƣớc thực đo trạm Châu Đốc (Đv: m) Mực nƣớc thực đo trạm Tân Châu (Đv: m) Hình 3.12 Tƣơng quan gi a diện tích ngập không che mây ĐBSC với số liệu quan trắc trạm thủy văn Châu Đốc Tân Châu 36 y = 420.27x + 831.42 R² = 0.7503 3000 Diện tích ngập Đơn vị (NGhìn ha) 3000 Diện tích ngập Đơn vị (NGhìn ha) y = 362.94x + 791.76 R² = 0.6957 2500 2000 1500 1000 500 2500 2000 1500 1000 500 0 Mực nƣớc thực đo trạm Châu Đốc (Đv: m) 4 Mực nƣớc thực đo trạm Tân Châu (Đv: m) Hình 3.13 Tƣơng quan gi a diện tích ngập có che mây ĐBSCL với số liệu quan trắc trạm thủy văn Châu Đốc Tân Châu Qua số liệu tƣơng quan diện tích ngập lũ có che mây diện tích ngập lũ không che mây trạm thủy văn Châu Đốc Tân Châu kết phân loại lũ che mây tƣơng quan với liệu phân loại lũ không che mây Kết có tƣơng quan thuận diện tích ngập lũ ĐBSCL với số liệu trạm quan trắc Tân Châu Châu Đốc với hệ số tƣơng quan không cao nhƣng chấp nhận đƣợc lần lƣợt 0,75 0,69 Số liệu mực nƣớc thủy văn diện tích ngập ĐBSCL tỷ lệ thuận với nhau, điều chứng tỏ mực nƣớc sông cao diện tích ngập lũ có xu hƣớng tăng theo ngƣợc lại mực nƣớc sông giảm lũ bắt đầu có tƣợng giảm xuống Với kết giải đoán, tính toán phân loại dựa kết hợp số EVI, LSWI DVEL việc xác định giám sát diễn biến lũ trở nên xác, cung cấp nhìn tổng quan trạng ngập lũ ĐBSCL với thay đổi không gian thời gian ngập năm 2013 Với kết so sánh, cho thấy việc tiến hành giám sát lũ khu vực ĐBSCL năm 2013 việc sử dụng ảnh vệ tinh MODIS hoàn toàn có sở tính khách quan 37 CHƢƠNG KẾT LUẬN VÀ KIẾN NGHỊ 4.1 Kết luận Qua kết phân loại từ ảnh MODIS với độ phân giải không gian 250m đề tài xác định đƣợc: - Diễn tiến lũ xảy ĐBSCL năm 2013, qua phát thời điểm mà bắt đầu ngập lũ cụ thể ngày 12/07/2013 với diện tích ngập lũ khoảng 1.359 nghìn đỉnh lũ đạt mức cực đại 2.767 nghìn vào ngày 01/11/2013 dần hạ xuống vào khoảng cuối tháng 11 tháng 12 - Thành lập đồ phân vùng ngập lũ ĐBSCL năm 2013 giải đoán từ ảnh MODIS đạt mức độ chi tiết cấp vùng, khả phát phân loại đƣợc đối tƣợng khác đó, có đối tƣợng ngập quanh năm hay vùng nuôi thủy sản chiếm 16,67% diện tích, đối tƣợng không ngập hoàn toàn chiếm 12,13% diện tích, đối tƣợng không ngập từ 4/7 - 17/11 chiếm 14,29% diện tích lại đối tƣợng ngập từ 12/7 - 17/11 chiếm 56,91% diện tích vùng Nghiên cứu cho thấy mối tƣơng quan ảnh theo không gian, thời gian với thay đổi lũ năm 2013 Kết nghiên cứu có độ tƣơng quan cao dựa so sánh tƣơng quan diện tích ngập mực nƣớc thực đo trạm thủy văn Tân Châu Châu Đốc với R2 Tân Châu Châu Đốc lần lƣợt 0,75 0,69 4.2 Kiến nghị Kết nghiên cứu đƣa nhìn tổng quan trạng ngập lũ nhƣ thay đổi không gian thời gian ngập lũ vào năm 2013 khu vực ĐBSCL Ảnh viễn thám khó xác định đƣợc độ sâu ngập, để đánh giá xác cƣờng độ ngập lũ cần có liệu hỗ trợ độ sâu ngập - Thời gian quan sát đối tƣợng nghiên cứu nên dài để quan sát đầy đủ diễn tiến ngập năm Qua đƣa kết đánh giá ảnh hƣởng lũ đến trình sản xuất nông nghiệp nhƣ dự báo tình hình lũ năm Nghiên cứu tìm hiểu phƣơng pháp loại bỏ ảnh hƣởng mây giúp cho kết giải đoán có độ xác 38 TÀI LIỆU THAM KHẢO Tài liệu nƣớc Hoàng Viết Giao Trần Bạch Giang ( 2004 ), Ứng dụng ảnh MODIS để theo dõi tàn ph rừng Miền Nam Brazil Huỳnh Thị Thu Hƣơng ctv, 2012 Ứng dụng ảnh MODIS theo dõi thay đổi nhiệt độ ề mặt đất tình hình hô hạn đồng ằng sông Cửu Long [24a 49-59] Lê Anh Tuấn, 2004 Phòng hống thiên tai Đại Học Cần Thơ Lê Anh Tuấn, 2008 Thủy văn môi trường Đại Học Cần Thơ Lê Sâm, 1996 Thủy văn ông trình NXB Khoa học kỹ thuật Hà Nội Lê Văn Trung ctv, 2006 Thự hành viễn th m NXB Đại Học Quốc Gia TP Hồ Chí Minh Nguyễn Đức Ngữ, 2009 Tr i đất ần húng ta! ãy liên ết hống lại iến đổi hí hậu Tổng Cụ môi trường Bộ Tài nguyên Môi trƣờng, Hà Nội Nguyễn Khắc Thời, 2011 Gi o trình iễn Th m Trƣờng Đại Học Nông Nghiệp Hà Nội Nguyễn Thị Kim Yến, 2013 Khảo s t ấu anh t Cửu Long Trƣờng Đại học Cần Thơ nông nghiệp đồng ằng sông Nguyễn Ngọc Thạch, 2005 Cơ sở viễn th m NXB Nông Nghiệp Hà Nội Nguyễn Thị Bích Liên, 2012 Ứng dụng ảnh vệ tinh MODIS gi m s t lũ đồng ằng sông Cửu Long Trƣờng Đại học Nông Lâm thành Phố Hồ Chí Minh Phan Thanh Nhàn, Võ Quang Minh, 2011 Theo dõi diễn tiến lũ lưu vự sông Mê ông làm sở dự o lũ Đồng Bằng Sông Cửu Long sử dụng ảnh viễn th m MODIS Kỷ yếu hội thảo Ứng dụng GIS toàn quốc 2011 Trần Nhƣ Hối, 2005 Xây dựng sở liệu mực nƣớc lũ vùng ngập lũ đồng sông Cửu Long nhằm đề xuất giải pháp khoa học xây dựng hệ thống đê bao Viện khoa học Thủy lợi miền Nam Trần Thanh Cảnh, 2000 Đồng Bằng Sông Cửu Long đón hào ỷ 21 NXB Văn Nghệ, thành phố Hồ Chí Minh Trần Tiễn Khanh, 2001 Nguyên nhân lũ lụt lớn Đồng ằng sông Cửu Long [Truy cập ngày 22/02/2013] Địa truy cập: http://www.vnbaolut.com/lulut_uni.htm 39 Viện Khoa học thủy lợi miền Nam, 2013 Tình hình lũ lụt đồng sông Cửu Long cuối tháng – đầu tháng 10/2103 Địa truy cập: http://www.siwrr.org.vn/tv3_files/ Dubaolu_thang9-10.2013 Võ Quang Minh, 2010 Gi o trình ỹ thuật viễn th m NXB Đại Học Cần Thơ Võ Quốc Thành, 2013 Quản lý giảm nhẹ t động lũ đồng ằng sông Cửu Long [Truy cập ngày 07/06/2013] Địa truy cập: http://luanvan.net.vn/luan-van/tieuluan-quan-ly-va-giam-nhe-tac-dong-cua-lu-dong-bang-song-cuu-long-36733 Vũ Hữu Long ctv, 2011 Sử dụng tư liệu ảnh vệ tinh MODIS nghiên ứu mùa vụ ây trồng, lập ảng đồ trạng iến động lớp phủ vùng đồng ằng Sông ồng giai đoạn 2008 – 2010 Hội thảo ứng dụng GIS toàn quốc Tài liệu nƣớc Islam S A., and Bala K S., Haque A., 2009 Flood inundation map of Bangladesh using MODIS surface reflectance data, 2nd International Conference on Water and Flood Management ICWFM 80: 245 - 265 Brochure, 2009 MODIS moderate – resolution imaging spectroradiometer Sakamoto T., et al., 2007 Detecting temporal changes in the extent of annual flooding within the Cambodia and the Vietnamese Mekong Delta from MODIS time-series imageo, Remote Sensing of Enviromen 109 (3): 295 - 313 Sakamoto T., et al., 2009 REVIEW Detection of Yearly Change in Farming Systems in the Vietnamese Mekong Delta from MODIS Time-Series Imagery Remote Sensing of Enviroment, JARQ 43 (3): 173 - 185 Zhuang R., Costelloe J., and Ryu D., 2011 Mapping Flood Events Using Remote Sensing in an Arid Zone River, Central Australia 19th International Congress on Modelling and Simulation, Perth, Australia, 12 - 16 December 2011 The University of Melbourne, pp 2030 - 2036 Các trang web http://ffw.mrcmekong.org http://modis.gsfc.nasa.gov http://reverb.echo.nasa.gov http://www.thoitietnguyhiem.net 40 PHỤ CHƢƠNG 41 Phụ lục 1: mực nƣớc thủy văn trạm Tân Châu Châu Đốc (Nguồn: http://ffw.mrcmekong.org/) Ngày 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 0,08 0,24 0,37 0,73 0,98 0,99 0,91 0,91 0,63 0,37 -0,09 -0,27 -0,21 -0,2 -0,2 -0,16 0,04 0,3 0,54 0,82 1,02 1,02 1,0 0,76 0,16 0,12 0,36 0,51 0,52 0,55 Sồ liệu trạm Châu Đốc Tháng 0,55 1,5 0,63 1,62 0,83 1,71 1,09 1,83 1,23 1,87 1,23 1,85 1,12 1,79 0,94 1,82 0,68 1,8 0,57 1,79 0,31 1,86 0,28 1,9 0,26 1,97 0,4 1,99 0,33 2,04 0,3 2,08 0,44 2,18 0,68 2,26 0,95 2,31 1,22 2,23 1,45 2,16 1,53 2,04 1,43 2,0 1,19 1,94 1,12 1,95 1,17 1,97 1,25 1,99 1,29 2,02 1,32 2,06 1,34 2,07 1,4 2,16 42 2,25 2,33 2,34 2.34 2,3 2,27 2,25 2,22 2,25 2,3 2,32 2,33 2,34 2,35 2,43 2,49 2,55 2,58 2,54 2,63 2,65 2,64 2,7 2,79 2,84 2,92 3,03 3,14 3,25 3,36 10 3,49 3,59 3,65 3,69 3,69 3,72 3,7 3,71 3,72 3,71 3,7 3,68 3,66 3,62 3,58 3,57 3,56 3,54 3,56 3,55 3,56 3,54 3,52 3,47 3,41 3,39 3,37 3,31 3,24 3,2 3,15 Ngày 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 0,17 0,32 0,43 0,72 0,93 0,91 0,86 0,77 0,39 0,02 -0,02 -0,11 -0,11 -0,1 -0,07 -0,04 0,13 0,38 0,58 0,82 0,97 0,94 0,91 0,76 0,43 0,34 0,65 0,77 0,8 0,81 Sồ liệu trạm Tân Châu Tháng 0,83 2,03 0,87 2,17 1,01 2,27 1,16 2,43 1,21 2,47 1,21 2,51 1,1 2,53 0,91 2,54 0,76 2,54 0,64 2,53 0,57 2,54 0,5 2,58 0,5 2,63 0,54 2,65 0,53 2,7 0,5 2,75 0,64 2,81 0,87 2,84 1,11 2,85 1,28 2,78 1,46 2,7 1,59 2,59 1,64 2,54 1,62 2,48 1,55 2,47 1,64 2,5 1,75 2,54 1,78 2,6 1,79 2,65 1,8 2,68 1,89 2,75 43 2,82 2,85 2,85 2,84 2,82 2,79 2,74 2,75 2,78 2,81 2,83 2,82 2,82 2,84 2,88 2,92 2,99 3,03 3,04 3,14 3,22 3,28 3,36 3,43 3,52 3,62 3,75 3,9 4,01 4,12 10 4,21 4,31 4,33 4,31 4,27 4,24 4,22 4,22 4,2 4,19 4,16 4,13 4,1 4,06 3,99 3,94 3,92 3,92 3,95 3,96 3,97 3,95 3,9 3,85 3,77 3,74 3,7 3,62 3,54 3,48 3,43 Phụ lục 2: Các ảnh MODIS thu thập STT 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 Tên ảnh MOD09A1.A2013185.h28v07.005.2013194140558.hdf MOD09A1.A2013185.h28v08.005.2013194140841.hdf MOD09A1.A2013193.h28v07.005.2013204212301.hdf MOD09A1.A2013193.h28v08.005.2013204213035.hdf MOD09A1.A2013201.h28v07.005.2013213180209.hdf MOD09A1.A2013201.h28v08.005.2013213180236.hdf MOD09A1.A2013209.h28v07.005.2013218193510.hdf MOD09A1.A2013209.h28v08.005.2013218193813.hdf MOD09A1.A2013217.h28v07.005.2013226124406.hdf MOD09A1.A2013217.h28v08.005.2013226123701.hdf MOD09A1.A2013225.h28v07.005.2013234131215.hdf MOD09A1.A2013225.h28v08.005.2013234130341.hdf MOD09A1.A2013233.h28v07.005.2013242172031.hdf MOD09A1.A2013233.h28v08.005.2013242172709.hdf MOD09A1.A2013241.h28v07.005.2013252120946.hdf MOD09A1.A2013241.h28v08.005.2013252120319.hdf MOD09A1.A2013249.h28v07.005.2013260221328.hdf MOD09A1.A2013249.h28v08.005.2013260221928.hdf MOD09A1.A2013257.h28v07.005.2013266131142.hdf MOD09A1.A2013257.h28v08.005.2013266125854.hdf MOD09A1.A2013265.h28v07.005.2013277212100.hdf MOD09A1.A2013265.h28v08.005.2013277213701.hdf MOD09A1.A2013273.h28v07.005.2013301134442.hdf MOD09A1.A2013273.h28v08.005.2013301134404.hdf MOD09A1.A2013281.h28v07.005.2013303134618.hdf MOD09A1.A2013281.h28v08.005.2013303132704.hdf MOD09A1.A2013289.h28v07.005.2013314140454.hdf MOD09A1.A2013289.h28v08.005.2013314140653.hdf MOD09A1.A2013297.h28v07.005.2013316135122.hdf MOD09A1.A2013297.h28v08.005.2013316134618.hdf MOD09A1.A2013305.h28v07.005.2013316145109.hdf MOD09A1.A2013305.h28v08.005.2013316145252.hdf MOD09A1.A2013313.h28v07.005.2013322112659.hdf MOD09A1.A2013313.h28v08.005.2013322113435.hdf MOD09A1.A2013321.h28v07.005.2013332221752.hdf MOD09A1.A2013321.h28v08.005.2013332223127.hdf MOD09A1.A2013329.h28v07.005.2013338220345.hdf MOD09A1.A2013329.h28v08.005.2013338220228.hdf 44 Thời gian thu ảnh 04/07/2013 04/07/2013 12/07/2013 12/07/1203 20/07/2013 20/07/2013 28/01/2012 28/01/2012 05/08/2013 05/08/2013 13/08/2013 13/08/2013 21/08/2013 21/08/2013 29/08/2013 29/08/2013 06/09/2013 06/09/2013 14/09/2013 14/09/2013 22/09/2013 22/09/2013 30/09/2013 30/09/2013 08/10/2013 08/10/2013 16/10/2013 16/10/2013 24/10/2013 24/10/2013 01/11/2013 01/11/2013 09/11/2013 09/11/2013 17/11/2013 17/11/2013 25/11/2013 25/11/2013 Trong đặc điểm ảnh đƣợc xác định nhƣ sau: Ví dụ nhƣ: MOD09A1.A2013185.h28v07.005.2013194140558.hdf MOD09A1 : tên ngắn gọn A2013185 : Ngày bắt đầu quét chụp (04/07/2013) h28v07 : Kí hiệu khu vực chụp 005 : Thế hệ 2013194140558.hdf : Thời gian xuất ảnh (14: 05: 55 ngày 12/07/2013) Phụ lục 3: Tƣơng quan diện tích ngập không che mây ĐBSCL với số liệu quan trắc trạm thủy văn Châu Đốc Tân Châu Ngày 4/7/2009 12/7/2009 20/7/2009 28/7/2009 5/8/2009 13/8/2009 21/8/2009 29/8/2009 6/9/2009 14/9/2009 22/9/2009 30/9/2009 8/10/2009 16/10/2009 24/10/2009 TÂN CHÂU STT Mực nƣớc 1,16 0,50 1,28 1,78 2,47 2,63 2,70 2,65 2,79 10 2,84 11 3,28 12 4,12 13 4,22 14 3,94 15 3,85 Diện tích 1.473,89 1.306,18 1.359,69 1.170,71 1.275,72 1.187,73 1.517,64 1.692,80 1.684,79 1.892,43 1.823,68 1.962,75 2.368,57 1.962,75 2.691,39 Ngày 4/7/2009 12/7/2009 20/7/2009 28/7/2009 5/8/2009 13/8/2009 21/8/2009 29/8/2009 6/9/2009 14/9/2009 22/9/2009 30/9/2009 8/10/2009 16/10/2009 24/10/2009 45 CHÂU ĐỐC STT Mực nƣớc 1,09 0,28 1,22 1,29 1,87 1,97 2,16 2,06 2,27 10 2,35 11 2,64 12 3,36 13 3,71 14 3,57 15 3,47 Diện tích 1.473,89 1.306,18 1.359,69 1.170,71 1.275,72 1.187,73 1.517,64 1.692,80 1.684,79 1.892,43 1.823,68 1.962,75 2.368,57 1.962,75 2.691,39 Phụ lục 4: Tƣơng quan diện tích ngập có che mây ĐBSCL với số liệu quan trắc trạm thủy văn Châu Đốc Tân Châu Ngày 4/7/2009 12/7/2009 20/7/2009 28/7/2009 5/8/2009 13/8/2009 21/8/2009 29/8/2009 6/9/2009 14/9/2009 22/9/2009 30/9/2009 8/10/2009 16/10/2009 24/10/2009 TÂN CHÂU STT Mực nƣớc 1,16 0,50 1,28 1,78 2,47 2,63 2,70 2,65 2,79 10 2,84 11 3,28 12 4,12 13 4,22 14 3,94 15 3,85 Diện tích 1.306,18 1.359,69 1.170,71 1.275,72 1.187,73 1.517,64 1.692,80 1.684,79 1.892,43 1.823,68 1.962,75 2.368,57 1.962,75 2.691,39 2.573,54 Ngày 4/7/2009 12/7/2009 20/7/2009 28/7/2009 5/8/2009 13/8/2009 21/8/2009 29/8/2009 6/9/2009 14/9/2009 22/9/2009 30/9/2009 8/10/2009 16/10/2009 24/10/2009 46 CHÂU ĐỐC STT Mực nƣớc 1,09 0,28 1,22 1,29 1,87 1,97 2,16 2,06 2,27 10 2,35 11 2,64 12 3,36 13 3,71 14 3,57 15 3,47 Diện tích 1.306,18 1.359,69 1.170,71 1.275,72 1.187,73 1.517,64 1.692,80 1.684,79 1.892,43 1.823,68 1.962,75 2.368,57 1.962,75 2.691,39 2.573,54 [...]... ngập lũ thay đổi theo không gian và thời gian ảnh hƣởng nhƣ thế nào đến quá trình sản xuất nông nghiệp và cơ cấu mùa vụ ở vùng ĐBSCL khá quan trọng Xuất phát từ vấn đề trên, đề tài nghiên cứu: “Giải đoán tình hình lũ lụt ở đồng bằng sông Cửu Long giai đoạn năm 2013 bằng ảnh MODIS đã đƣợc thực hiện nhằm mục đích sau: - Theo dõi quá trình diễn tiến lũ theo thời gian ở đồng bằng sông Cửu Long năm 2013. .. điều tra khảo sát đƣợc 1.2.2 iện trạng ngập lũ ở đồng bằng Sông Cửu Long Lũ lụt là thiên tai lớn nhất đe dọa nƣớc ta trong nhiều năm qua gây nhiều thiệt hại về ngƣời và tài sản Theo Trần Tiễn Khanh (2001), ở Miền Nam những năm gần đây đã có những trận lũ lụt lớn chƣa từng thấy xảy ra tại vùng đồng bằng sông Cửu Long (ĐBSCL) Đặc biệt là cơn lũ lụt từ tháng 9 đến tháng 11 năm 2000 đƣợc gọi cơn lũ thế 4... một tiến trình lũ Theo Lê Anh Tuấn (2004), lũ đƣợc phân biệt thành các loại sau: - Lũ nhỏ: là loại lũ có đỉnh lũ thấp hơn mức đỉnh lũ trung bình nhiều năm - Lũ vừa: là loại lũ có đỉnh lũ đạt mức đỉnh lũ nhiều năm - Lũ lớn: là loại lũ có đỉnh lũ cao hơn mức đỉnh lũ nhiều năm - Lũ đặc biệt lớn: là loại lũ cao đỉnh lũ cao hiếm thấy trong thời kỳ quan trắc - Lũ lịch sử: là loại lũ có đỉnh lũ cao nhất trong... 04/07 /2013) 24 3.2 Kết quả tạo ảnh LSWI khu vực nghiên cứu (ngày 04/07 /2013) 25 3.3 Kết quả tạo ảnh DVEL khu vực nghiên cứu (ngày 04/07 /2013) 26 3.4 Kết quả phân loại đối tƣợng dựa trên kết hợp các chỉ số EVI, LSWI, DVEL (ngày 04/07 /2013) 27 3.5 Bản đồ ngập lũ từng giai đoạn ĐBSCL năm 2013 không che mây 28 3.6 Bản đồ ngập lũ từng giai đoạn ĐBSCL năm 2013 có che mây 28 3.7 Diện tích ngập lũ Đồng Bằng Sông. .. trong khí đó lũ vừa và nhỏ nhiều hơn (Viện Khoa học thủy lợi miền Nam, 2013) 1.2.3 Nguyên nhân gây ra lũ ở đồng bằng sông Cửu Long Mƣa lớn và kéo dài (do bão lớn) là nguyên nhân chính gây ra lũ lụt, ngoài ra ở vùng đồng bằng cửa sông tiếp giáp với biển, triều cƣờng là một nhân tố làm lũ lụt trầm trọng hơn (theo Lê Anh Tuấn, 2004) Ngoài ra, còn một số yếu tố khác ảnh hƣởng đến khả năng xuất hiện lũ lớn và... ngập lũ ở Bangladesh - Ứng dụng của ảnh MODIS ở Việt Nam Sử dụng tƣ liệu ảnh vệ tinh MODIS nghiên cứu mùa vụ cây trồng, lập bảng đồ hiện trạng và biến động lớp phủ vùng đồng bằng Sông Hồng giai đoạn 2008 - 2010 (Vũ Hữu Long và ctv, 2011) Bộ Tài Nguyên và Môi trƣờng nghiên cứu ứng dụng ảnh vệ tinh MODIS để phòng chống ngập lụt và theo dõi môi trƣờng ở ĐBSCL (Lê Minh và ctv, 2005) Ứng dụng ảnh MODIS. .. của ảnh MODIS đƣợc sử dụng cho tính toán chỉ số thực vật Các kênh phổ của ảnh MODIS đƣợc sử dụng cho tính toán 3.3 Diện tích không ngập có che mây và không che mây ở ĐBSCL 3.4 Bảng Diện tích phân vùng ngập lũ ĐBSCL năm 2013 3.1 xi Trang 13 23 24 29 33 DANH SÁCH HÌNH Tiêu đề Hình Trang 1.1 Bản đồ các tỉnh Đồng Bằng Sông Cửu Long 2 1.2 Đồ thị diễn tả một tiến trình lũ 4 1.3 Một số trận lũ lớn khu vực Đồng. .. d liệu - Ảnh MODIS – MOD09A1 - Bản đồ hành chính ĐBSCL 2013 - Số liệu mực nƣớc Châu Đôc, Tân Châu 2013 - Bản đồ hiện trạng sử dụng đất ở ĐBSCL năm 2013 (Nguyễn Thị KimYến, 2013) Xử lý d liệu - Tiền xử lý ảnh viễn thám - Tạo ảnh EVI, tạo chuỗi EVI - Tạo ảnh LSWI - Tạo ảnh DVEL, tạo chuỗi DVEL Bản đồ hiện trạng sử dụng đất năm 2013 (Nguyễn Thị Kim Yến, 2013) Chuỗi ảnh DVEL Ảnh DVEL - Phân cấp ảnh DVEL... tháng 7 Trong những năm gần đây, lũ thể xảy ra muộn hơn nhƣ năm 2011, và một số năm không có lũ hay lũ rất nhỏ, chẳng hạn nhƣ 1998, 2012 Lũ đồng bằng lên và xuống chậm, theo đúng bản chất của lũ lƣu vực Nói chung, cƣờng suất lũ từ 2 - 3 cm đến 10 - 15 cm/ngày Lũ đầu vụ có cƣờng suất lớn hơn lũ chính vụ Trong thời gian qua, lũ trên đồng bằng đang có những biến đổi khác với trƣớc đây, lũ lớn dƣờng nhƣ xuất... xạ hàng năm là 140 - 150 Kcal/cm2 /năm 2 Vùng ĐBSCL có mạng lƣới sông khá phức tạp, trong đó chủ yếu là sông Cửu Long và các chi lƣu của nó Ở thƣợng nguồn trên Tân Châu (sông Tiền) và Châu Đốc (sông Hậu) khi chảy vào đồng bằng, sông có bề rộng khoảng từ 60 m đến 300 m và dần dần mở rộng khi chảy về dƣới hạ lƣu, bề rộng sông khoảng 2 km khi ra đến biển, đoạn lớn nhất là cửa sông Hậu, bề rộng sông lên ... NHIÊN BỘ MÔN TÀI NGUYÊN ĐẤT ĐAI - DANH XÀ RÂY KHẢO SÁT TÌNH HÌNH LŨ LỤT Ở ĐỒNG BẰNG SÔNG CỬU LONG GIAI ĐOẠN NĂM 2013 BẰNG ẢNH MODIS LUẬN VĂN KỸ SƢ NGÀNH QUẢN LÝ ĐẤT ĐAI Mã ngành: 52850103... 04/07 /2013) 27 3.5 Bản đồ ngập lũ giai đoạn ĐBSCL năm 2013 không che mây 28 3.6 Bản đồ ngập lũ giai đoạn ĐBSCL năm 2013 có che mây 28 3.7 Diện tích ngập lũ Đồng Bằng Sông Cửu Long năm 2013 30... lũ lụt đồng sông Cửu Long giai đoạn năm 2013 ảnh MODIS đƣợc thực nhằm mục đích sau: - Theo dõi trình diễn tiến lũ theo thời gian đồng sông Cửu Long năm 2013 - xây dựng đồ trạng ngập lũ đồ phân

Ngày đăng: 13/11/2015, 11:10

Từ khóa liên quan

Tài liệu cùng người dùng

  • Đang cập nhật ...

Tài liệu liên quan