Ứng dụng mô hình oryza2000 kết hợp ảnh viễn thám MODIS thành lập bản đồ năng suất lúa tỉnh sóc trăng

92 576 7
Ứng  dụng  mô  hình  oryza2000  kết  hợp  ảnh viễn thám MODIS thành lập bản đồ năng suất lúa tỉnh sóc trăng

Đ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

Ứng dụng mô hình oryza2000 kết hợp ảnh viễn thám MODIS thành lập bản đồ năng suất lúa tỉnh sóc trăng Ứng dụng mô hình oryza2000 kết hợp ảnh viễn thám MODIS thành lập bản đồ năng suất lúa tỉnh sóc trăng Ứng dụng mô hình oryza2000 kết hợp ảnh viễn thám MODIS thành lập bản đồ năng suất lúa tỉnh sóc trăng Ứng dụng mô hình oryza2000 kết hợp ảnh viễn thám MODIS thành lập bản đồ năng suất lúa tỉnh sóc trăng Ứng dụng mô hình oryza2000 kết hợp ảnh viễn thám MODIS thành lập bản đồ năng suất lúa tỉnh sóc trăng Ứng dụng mô hình oryza2000 kết hợp ảnh viễn thám MODIS thành lập bản đồ năng suất lúa tỉnh sóc trăng Ứng dụng mô hình oryza2000 kết hợp ảnh viễn thám MODIS thành lập bản đồ năng suất lúa tỉnh sóc trăng Ứng dụng mô hình oryza2000 kết hợp ảnh viễn thám MODIS thành lập bản đồ năng suất lúa tỉnh sóc trăng Ứng dụng mô hình oryza2000 kết hợp ảnh viễn thám MODIS thành lập bản đồ năng suất lúa tỉnh sóc trăng Ứng dụng mô hình oryza2000 kết hợp ảnh viễn thám MODIS thành lập bản đồ năng suất lúa tỉnh sóc trăng Ứng dụng mô hình oryza2000 kết hợp ảnh viễn thám MODIS thành lập bản đồ năng suất lúa tỉnh sóc trăng Ứng dụng mô hình oryza2000 kết hợp ảnh viễn thám MODIS thành lập bản đồ năng suất lúa tỉnh sóc trăng Ứng dụng mô hình oryza2000 kết hợp ảnh viễn thám MODIS thành lập bản đồ năng suất lúa tỉnh sóc trăng Ứng dụng mô hình oryza2000 kết hợp ảnh viễn thám MODIS thành lập bản đồ năng suất lúa tỉnh sóc trăng

TRƯỜNG ĐẠI HỌC CẦN THƠ KHOA MÔI TRƯỜNG VÀ TÀI NGUYÊN THIÊN NHIÊN -- NGÔ HỮU LỢI ỨNG DỤNG MÔ HÌNH ORYZA2000 KẾT HỢP ẢNH VIỄN THÁM MODIS THÀNH LẬP BẢN ĐỒ NĂNG SUẤT LÚA TỈNH SÓC TRĂNG LUẬN VĂN TỐT NGHIỆP CAO HỌC NGÀNH QUẢN LÝ ĐẤT ĐAI Năm 2015 TRƯỜNG ĐẠI HỌC CẦN THƠ KHOA MÔI TRƯỜNG VÀ TÀI NGUYÊN THIÊN NHIÊN -- NGÔ HỮU LỢI ỨNG DỤNG MÔ HÌNH ORYZA2000 KẾT HỢP ẢNH VIỄN THÁM MODIS THÀNH LẬP BẢN ĐỒ NĂNG SUẤT LÚA TỈNH SÓC TRĂNG LUẬN VĂN TỐT NGHIỆP CAO HỌC NGÀNH QUẢN LÝ ĐẤT ĐAI MÃ SỐ: 60850103 HƯỚNG DẪN KHOA HỌC Gs Ts NGÔ NGỌC HƯNG Năm 2015 CHẤP THUẬN LUẬN VĂN CỦA HỘI ĐỒNG Luận văn này, với đề tựa là“Ứng dụng mô hình Oryza2000 kết hợp ảnh viễn thám MODIS thành lập đồ suất lúa tỉnh Sóc Trăng” học viên Ngô Hữu Lợi thực theo hƣớng dẫn Gs Ts Ngô Ngọc Hƣng Luận văn báo cáo đƣợc Hội đồng chấm luận văn thông qua ngày…………………… Uỷ viên Thƣ ký Phản biện Phản biện Cán hƣớng dẫn Chủ Tịch Hội Đồng Gs Ts NGÔ NGỌC HƢNG i LỜI CẢM TẠ Trãi qua thời gian học tập trƣờng Đại học Cần Thơ giúp hoàn thiện thêm kiến thức chuyên ngành Để đƣợc kết nhƣ ngày hôm nổ lực thân, giúp đỡ tận tình quý thầy cô, gia đình, bạn bè Tôi xin chân thành cảm ơn đến thầy Ngô Ngọc Hƣng tận tình hƣớng dẫn tạo điều kiện tốt để hoàn thành đề tài luận văn tốt nghiệp cách thuận lợi; thầy Võ Quang Minh 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, toàn thể quý thầy cô giảng dạy lớp Cao học Quản lý Đất đai K20 truyền đạt kiến thức sâu rộng cho làm hành trang sau phục vụ cho xã hội, cho đất nƣớc Các bạn lớp Cao học Quản lý Đất đai K20 động viên, tạo niềm tin giúp đỡ suốt thời gian qua Cuối kính gởi lòng biết ơn sâu sắc đến Cha, Mẹ ngƣời thân động viên nhắc nhở suốt trình học tập tạo điều kiện tốt để có kết nhƣ ngày hôm Chân thành cảm ơn! Ngô Hữu Lợi ii LÝ LỊCH KHOA HỌC I LÝ LỊCH SƠ LƢỢC Họ tên: Ngô Hữu Lợi Giới tính: Nam Ngày, tháng, năm sinh: 26/10/1990 Nơi sinh: Cần Thơ Quê quán: Cờ Đỏ - Cần Thơ Dân tộc: Kinh Thƣờng trú: số 48, ấp Thới Hữu, xã Thới Đông, huyện Cờ Đỏ, TP.Cần Thơ Điện thoại di động: 0982 866 255 E-mail: hloiproct90@gmail com II QUÁ TRÌNH ĐÀO TẠO Cao đẳng Hệ đào tạo: Chính quy Thời gian đào tạo: 2008 - 2011 Nơi học: Cao đẳng Kinh tế - Kỹ thuật Cần Thơ Ngành học: Quản lý Đất đai Tên môn thi tốt nghiệp: Quy hoạch sử dụng đất, Đo đạc địa Ngƣời hƣớng dẫn: Nguyễn Doanh Quyền, Phạm Văn Bé Đại học Hệ đào tạo: Chính quy Thời gian đào tạo: 2011 – 2013 Nơi học: Đại học Cần Thơ Ngành học: Quản lý Đất đai Tên luận văn: Đánh giá hiệu mô hình sử dụng đất xã Thới Đông, huyện Cờ Đỏ, thành phố Cần Thơ Ngƣời hƣớng dẫn: Th.S Nguyễn Thị Song Bình Cao học Hệ đào tạo: Chính quy Thời gian đào tạo: 2013 - 2015 Nơi học: Đại học Cần Thơ Ngành học: Quản lý Đất đai Tên luận văn: Ứng dụng ảnh Viễn thám Modis kết hợp mô hình Oryza2000 thành lập đồ suất lúa tỉnh Sóc Trăng Ngƣời hƣớng dẫn: Gs.Ts Ngô Ngọc Hƣng Trình độ ngoại ngữ : Anh văn trình độ B quốc gia, trình độ B1 khung tham chiếu Châu Âu iii LỜI CAM ĐOAN Tôi xin cam đoan công trình nghiên cứu thân Các số liệu, kết trình bày luận văn trung thực chƣa đƣợc công bố công trình nghiên cứu trƣớc Tác giả luận văn Ngô Hữu Lợi iv TÓM TẮT Nghiên cứu sử dụng phƣơng pháp giải đoán ảnh viễn thám Modis (MOD09Q1) độ phân giải 250m, chu kỳ lặp ngày để theo dõi biến động cấu mùa vụ tỉnh Sóc Trăng từ tháng 9/2012 đến tháng 12/2013 Kết giải đoán đƣợc kiểm tra thực địa nhằm đánh giá khả sử dụng ảnh Modis nghiên cứu thành lập đồ cấu mùa vụ nhƣ qua theo dõi giai đoạn sinh trƣởng lúa nghiên cứu cho thấy khả ứng dụng ảnh viễn thám Modis theo dõi sinh trƣởng lúa với độ xác cao Việc giải đoán ảnh sau đƣợc kết hợp với mô hình dự báo suất lúa Oryza2000 để dự báo sản lƣợng lúa Với điểm thu mẫu trạng suất thu hoạch kết đƣợc đánh giá, so sánh với số liệu thống kê địa phƣơng có tƣơng quan cao Điều cho thấy khả ứng dụng ảnh viễn thám Modis (MOD09Q1) kết hợp với mô hình Oryza2000 để theo dõi tình hình sản xuất lúa phạm vi rộng lớn, cấp vùng cấp quốc gia Kết nghiên cứu làm tảng thực nghiên cứu vào công tác theo dõi diện tích, đề sách sản xuất lúa cách hiệu tỉnh Sóc Trăng nói riêng ĐBSCL nói chung Từ khóa: Oryza2000, ảnh viễn thám Modis, sinh trƣởng lúa v ABSTRACT This research method is used MODIS remote sensing image interpretation (MOD09Q1) resolution of 250m, cycle repeats days to track seasonal changes structure in Soc Trang from month 9/2012 to month 12/2013 Results interpretation has been field testing to assess usability research MODIS mapping crop structure and thereby track the growth stage of the rice plant Research shows that the applicability of MODIS remote sensing in monitoring the growth of rice plants with high accuracy The image interpretation then be combined with predictive models to yield Oryza2000 forecast paddy output With 9-point sampling status and yield results are evaluated, compared with the statistics locally high correlation This suggests that the ability MODIS remote sensing applications (MOD09Q1) model combined with Oryza2000 to monitor rice production on a large scale, regional or national level Research results can implement the platform further studies on the monitoring area, set rice production policies effective in Soc Trang province in particular and the region in general Key words: Oryza2000, sensing image MODIS, growing rice vi MỤC LỤC CHẤP THUẬN LUẬN VĂN CỦA HỘI ĐỒNG i LỜI CẢM TẠ ii LÝ LỊCH KHOA HỌC iii LỜI CAM ĐOAN iv TÓM TẮT v ABSTRACT vi MỤC LỤC vii Danh sách bảng x Danh sách hình xi Danh mục từ viết tắt xiv CHƢƠNG MỞ ĐẦU 1.1 Đặt vấn đề 1.2 Mục tiêu phạm vi nghiên cứu 1.2.1 Mục tiêu tổng quát 1.2.2 Mục tiêu cụ thể 1.2.3 Phạm vi nghiên cứu 1.3 Ý nghĩa nghiên cứu 1.4 Thời gian thực CHƢƠNG LƢỢC KHẢO TÀI LIỆU 2.1 Định nghĩa đất đai 2.2 Vai trò đất đai 2.3 An ninh lƣơng thực 2.4 Mô hình Oryza2000 2.5 Tổng quan ảnh viễn thám 10 2.5.1 Khái niệm ảnh viễn thám 10 2.5.2 Nguyên lý viễn thám 11 2.5.3 Ảnh viễn thám Modis 12 2.5.4 Ảnh số thực vật (Normal Different Vegetation Index - NDVI) 12 2.6 Các giai đoạn sinh trƣởng lúa 13 2.7 Một số nghiên cứu ứng dụng nƣớc 15 vii 2.8 Đặc điểm vùng nghiên cứu 17 2.8.1 Vị trí địa lý 17 2.8.2 Đặc điểm địa hình 18 2.8.3 Khí hậu 18 2.8.4 Tài nguyên thiên nhiên 19 2.8.5 Hệ thống sông ngòi 19 2.8.6 Tình hình sản xuất lúa tỉnh Sóc Trăng 19 CHƢƠNG PHƢƠNG TIỆN, PHƢƠNG PHÁP 21 3.1 Phƣơng tiện 21 3.2 Phƣơng pháp 21 3.2.1 Phƣơng pháp kế thừa số liệu 21 3.2.2 Phƣơng pháp thu thập ảnh MODIS 21 3.2.3 Phƣơng pháp chọn mẫu nghiên cứu 23 3.2.4 Phƣơng pháp xác định suất thực tế đồng 25 3.2.5 Thành lập đồ cấu mùa vụ phân bố suất 27 3.2.6 Phƣơng pháp kiểm tra 28 CHƢƠNG KẾT QUẢ THẢO LUẬN 29 4.1 Xử lý ảnh Modis 29 4.2 Thành lập đồ trạng cấu mùa vụ canh tác lúa tỉnh Sóc Trăng 2013 36 4.3 Mô suất lúa mô hình Oryza2000 44 4.3.1 Dữ liệu đầu vào Oryza2000 44 4.3.2 Các bƣớc sử dụng Oryza2000 46 4.3.3 Kết mô suất 47 4.3.4 Năng suất thu thực tế 50 4.4 Thành lập đồ phân bố suất lúa tỉnh Sóc Trăng 51 4.4.1 Phân bố suất vụ Đông Xuân 51 4.4.2 Phân bố suất vụ Xuân Hè 53 4.4.3 Phân bố suất vụ Hè Thu 55 4.5 Đánh giá độ tin cậy 57 viii - Nguyễn Quang Tuấn, 2008 Hƣớng dẫn sử dùng phần mềm xử lý tƣ liệu viễn thám Trƣờng Đại học Huế - Phạm Thị Lệ Huyền, 2014 Ứng dụng mô hình ORYZA mô suất lúa tỉnh Sóc Trăng theo hƣớng phục vụ bảo hiểm nông nghiệp Luận văn tốt nghiệp Cao học Quản lý đất đai, trƣờng Đại học Cần Thơ - Sóc Trăng sau 20 năm tái lập - Một chặng đƣờng phát triển, 2011 Nhà xuất Cục thống kê Sóc Trăng - Trần Thị Hiền, 2010 Nghiên cứu ảnh viễn thám MODIS theo dõi tiến độ xuống giống vùng đất trồng lúa đồng song Cửu Long Luận văn Thạc sĩ ngành Khoa học đất, trƣờng Đại học Cần Thơ - Văn kiện Đại hội đại biểu toàn quốc lần thứ X, 2006 Nhà xuất trị Quốc gia - Văn kiện Đại hội đại biểu toàn quốc lần thứ XI, 2011 Nhà xuất trị Quốc gia - Võ Quang Minh, 2005 Giáo trình Viễn thám 1, Khoa Nông nghiệp SHƢD, trƣờng Đại học Cần Thơ - Võ Quang Minh, 2010 Giáo trình kỹ thuật Viễn thám Nhà xuất Đại học Cần Thơ Tiếng Anh - Bouman, 2001 ORYZA2000: modeling lowland rice - C Wikarmpapraharn and E Kositsakulchai (2010) Evaluation of ORYZA2000 and CERES-Rice Models under Potential Growth Condition in the Central Plain of Thailand Evaluation of ORYZA2000 and CERES-Rice Models under Potential Kasetsart University, Kamphaeng Saen Campus, Nakhon Pathom 73140, Thailan Thai Journal of Agricultural Science 2010, 43(1): 17-29 - Decagon Devices, 2014 AccuPAR and Ceptometer are registered trademarks of Decagon Devices - Ebrahim Amiri And Mojtaba Rezaei (2010) Evaluation of Water-Nitrogen Schemes for Rice in Iran, Using ORYZA2000 Model Department of Agriculture, Islamic Azad University, Lahijan Branch, Lahijan, Iran Communications in Soil Science and Plant Analysis, 41: 20, 2459 — 2477 - Mojtaba Esmaeilzad Limoudehi, Ebrahim Amiri, Majid Ashori (2014) Evaluating ORYZA2000 model in shifting Irrigation conditions of Bahar variety 62 rice Journal of Applied Science and Agriculture, 9(7) Special 2014, Pages: 2779-2784 - Parida, B.R.B Oinam, N.R Patel, N Sharma, R Kandwal and M.K Hazarika.2008 Land surface temperature variation in relation to vegetation type using MODIS satellite data in Gujarat state of India International Jounal of Remote Sensing, Vol 29, pp, 4219-4235 - Song, C., Woodcock, and C E., Seto (2001) Classification and change detection using landsat TM data: when and how to correct atmospheric effects?, Remote Sensing of Environment 75, pp 230-244 - Thomas M.Lillesand & Ralph W.Kiefer (1996), Remote sensing and image interpretation, University of Wisconsin, Madison, USA Địa Web - Theo dõi lúa Truy cập trang: http://theodoilua.blogspot.com/ - Cổng thông tin điện tử tỉnh Sóc Trăng Truy cập trang http://www.soctrang.gov.vn/wps/portal/ - An ninh lƣơng thực ―Không số lƣợng‖ Truy cập trang: http://baochinhphu.vn/Quoc-te/An-ninh-luong-thuc-Khong-chi-la-soluong/177100.vgp 63 PHỤ LỤC Phụ lục 1: Phiếu điều tra Tên chủ hộ: ; Mô hình sản xuất: Thâm niên canh tác:……(năm) I ĐẶC ĐIỂM THỬA RUỘNG - Vị trí đất tại: ấp: .; xã :…… .; huyện:… ; tỉnh: - Hiện trạng canh tác - Diện tích: - Chiều cao lúa:…………………………………………………………………… - Mật độ gieo sạ…………………………………………………… - Mực nƣớc ruộng: ……………………………………………………………… - Dịch hại: …………………………………………………………………… … - Phân bón: ………………………………………………………………… …… - Tọa độ (4 góc ruộng): x: ; y: x: ; y: x: ; y: x: ; y: II NHỮNG THAY ĐỔI TRONG QUÁ TRÌNH SỬ DỤNG ĐẤT III LỊCH THỜI VỤ (Dƣơng lịch) Tháng Vụ Vụ Vụ 64 10 11 12 Phụ lục 2: Địa Website tải Modis 2013 http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.06.26/MOD09Q1.A2013177.h28v08.005.2013186125708.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.06.26/MOD09Q1.A2013177.h28v07.005.2013186130828.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.07.04/MOD09Q1.A2013185.h28v08.005.2013194140843.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.07.04/MOD09Q1.A2013185.h28v07.005.2013194140559.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.07.12/MOD09Q1.A2013193.h28v07.005.2013204212301.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.07.12/MOD09Q1.A2013193.h28v08.005.2013204213037.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.07.20/MOD09Q1.A2013201.h28v07.005.2013213180210.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.07.20/MOD09Q1.A2013201.h28v08.005.2013213180237.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.07.28/MOD09Q1.A2013209.h28v07.005.2013218193510.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.07.28/MOD09Q1.A2013209.h28v08.005.2013218193814.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.08.05/MOD09Q1.A2013217.h28v07.005.2013226124407.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.08.05/MOD09Q1.A2013217.h28v08.005.2013226123701.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.08.13/MOD09Q1.A2013225.h28v08.005.2013234130341.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.08.13/MOD09Q1.A2013225.h28v07.005.2013234131217.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.08.21/MOD09Q1.A2013233.h28v07.005.2013242172033.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.08.21/MOD09Q1.A2013233.h28v08.005.2013242172710.hdf 65 http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.08.29/MOD09Q1.A2013241.h28v07.005.2013252120946.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.08.29/MOD09Q1.A2013241.h28v08.005.2013252120319.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.09.06/MOD09Q1.A2013249.h28v07.005.2013260221328.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.09.06/MOD09Q1.A2013249.h28v08.005.2013260221929.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.09.14/MOD09Q1.A2013257.h28v07.005.2013266131143.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.09.14/MOD09Q1.A2013257.h28v08.005.2013266125854.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.09.22/MOD09Q1.A2013265.h28v07.005.2013277212100.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.09.22/MOD09Q1.A2013265.h28v08.005.2013277213705.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.09.30/MOD09Q1.A2013273.h28v07.005.2013301134442.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.09.30/MOD09Q1.A2013273.h28v08.005.2013301134404.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.10.08/MOD09Q1.A2013281.h28v07.005.2013303134618.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.10.08/MOD09Q1.A2013281.h28v08.005.2013303132704.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.10.16/MOD09Q1.A2013289.h28v07.005.2013314140454.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.10.16/MOD09Q1.A2013289.h28v08.005.2013314140653.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.10.24/MOD09Q1.A2013297.h28v08.005.2013316134620.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.10.24/MOD09Q1.A2013297.h28v07.005.2013316135122.hdf 66 http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.11.01/MOD09Q1.A2013305.h28v08.005.2013316145252.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.11.01/MOD09Q1.A2013305.h28v07.005.2013316145109.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.11.09/MOD09Q1.A2013313.h28v07.005.2013322112659.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.11.09/MOD09Q1.A2013313.h28v08.005.2013322113436.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.11.17/MOD09Q1.A2013321.h28v07.005.2013332221752.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.11.17/MOD09Q1.A2013321.h28v08.005.2013332223127.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.11.25/MOD09Q1.A2013329.h28v08.005.2013338220229.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.11.25/MOD09Q1.A2013329.h28v07.005.2013338220345.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.12.03/MOD09Q1.A2013337.h28v08.005.2013346114908.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.12.03/MOD09Q1.A2013337.h28v07.005.2013346120914.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.12.11/MOD09Q1.A2013345.h28v08.005.2013354123420.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.12.11/MOD09Q1.A2013345.h28v07.005.2013354125151.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.12.19/MOD09Q1.A2013353.h28v08.005.2013363072142.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.12.19/MOD09Q1.A2013353.h28v07.005.2013363075536.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.12.27/MOD09Q1.A2013361.h28v07.005.2014006135646.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 3.12.27/MOD09Q1.A2013361.h28v08.005.2014006141402.hdf 67 http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 4.01.01/MOD09Q1.A2014001.h28v08.005.2014012032932.hdf http://e4ftl01.cr.usgs.gov//MODIS_Composites/MOLT/MOD09Q1.005/201 4.01.01/MOD09Q1.A2014001.h28v07.005.2014012032755.hdf 68 Phụ lục 3: Lịch lấy mẫu giá trị LAI qua mùa Điểm Đông Xuân Ngày lấy mẫu (ngày/tháng/năm) D1 D2 D3 D4 D5 Xuân Hè LAI (m2 lá/m2 đất) Rep1 Rep2 Rep3 9/11/2012 1,2 1,05 0,81 26/11/2012 3,32 3,16 11/12/2012 4,09 8/10/2012 Ngày lấy mẫu (ngày/tháng/năm) Hè Thu LAI (m2 lá/m2 đất) Rep3 Ngày lấy mẫu (ngày/tháng/năm) LAI (m2 lá/m2 đất) Rep1 Rep2 Rep1 Rep2 Rep3 18/03/2013 0,96 0,75 0,7 6/7/2013 1,57 0,41 0,71 3,05 2/4/2013 5,43 5,88 6,06 23/07/2013 4,58 3,97 2,07 3,44 3,29 18/04/2013 3,93 5,86 3,67 8/8/2013 5,25 4,49 0,26 0,42 0,31 29/01/2013 0,27 0,13 0,11 5/6/2013 0,14 0,18 0,2 23/10/2012 1,51 1,94 2,11 13/02/2013 0,82 0,79 0,55 21/06/2013 1,19 2,73 1,68 9/11/2012 1,84 2,34 2,2 2/3/2013 1,45 3,8 2,5 6/7/2013 4,14 4,29 4,71 18/03/2013 2,83 2,35 2,81 23/07/2013 4,83 4,44 5,31 23/10/2012 1,33 1,51 1,24 13/02/2013 0,38 0,26 0,96 21/06/2013 0,45 0,32 0,57 9/11/2012 3,66 3,93 3,13 2/3/2013 2,13 1,58 1,71 6/7/2013 2,17 1,72 2,67 26/11/2012 5,19 6,28 5,86 18/03/2013 3,08 2,94 2,74 23/07/2013 2,93 4,96 2,26 2/4/2013 3,66 2,46 2,55 8/8/2013 3,12 4,24 2,72 26/11/2012 4,89 2,18 4,88 20/05/2013 0,1 0,14 0,11 11/12/2012 4,1 4,31 5,05 5/6/2013 0,54 0,98 0,69 28/12/2012 6,22 7,75 6,82 21/06/2013 1,51 2,06 2,04 6/7/2013 5,21 3,96 4,88 23/07/2013 6,49 5,38 6,83 26/11/2012 1,78 1,56 1,1 6/7/2013 0,71 0,69 0,99 11/12/2012 2,47 2,96 2,64 23/07/2013 0,71 1,24 1,18 69 Điểm Đông Xuân Ngày lấy mẫu (ngày/tháng/năm) 28/12/2012 D6 D7 D8 D9 Xuân Hè LAI (m2 lá/m2 đất) Rep1 Rep2 Rep3 4,21 3,44 3,08 Hè Thu LAI (m2 lá/m2 đất) Ngày lấy mẫu (ngày/tháng/năm) Rep1 Rep2 Rep3 Ngày lấy mẫu (ngày/tháng/năm) LAI (m2 lá/m2 đất) Rep1 Rep2 Rep3 8/8/2013 1,59 2,54 1,52 24/08/2013 4,61 4,51 23/10/2012 2,12 2,38 2,54 1/2/2013 0,64 0,91 0,78 21/06/2013 0,68 0,61 0,61 9/11/2012 3,66 3,59 3,76 13/02/2013 1,79 1,9 1,51 6/7/2013 1,03 1,43 1,09 26/11/2012 2,37 1,75 3,43 2/3/2013 4,31 5,23 4,08 23/07/2013 4,86 4,64 3,63 18/03/2013 4,21 2,77 3,18 8/8/2013 4,79 4,46 3,64 8/10/2012 0,94 1,11 0,51 1/2/2013 2,31 2,94 2,45 21/06/2013 1,08 0,82 0,71 23/10/2012 2,45 3,03 2,85 13/02/2013 4,16 3,56 3,18 6/7/2013 1,18 1,64 0,99 9/11/2012 3,55 2,95 2,79 2/3/2013 3,65 2,98 2,73 23/07/2013 3,45 2,53 3,33 23/10/2012 0,74 1,33 1,05 23/07/2013 2,64 3,03 3,47 9/11/2012 3,1 2,04 2,07 8/8/2013 3,8 4,97 4,42 26/11/2012 3,68 3,86 4,23 24/08/2013 4,23 3,44 3,36 9/11/2012 1,51 1,11 1,07 8/8/2013 1,54 1,44 1,34 26/11/2012 3,93 2,94 4,58 24/08/2013 2,08 1,43 1,39 11/12/2012 4,3 4,8 4,28 9/9/2013 3,12 2,91 3,28 70 Phụ lục 4: Nƣớc tƣới đạm qua mùa Đông Xuân Tƣới Điểm D1 D2 Xuân Hè Tƣới Phân bón Ngày/tháng/năm Số lƣợng (mm) 15/10/2012 21/10/2012 26/10/2012 02/11/2012 04/11/2012 16/11/2012 23/11/2012 30/11/2012 07/12/2012 14/12/2012 21/12/2012 28/12/2012 23/09/2012 6/10/2012 16/10/2012 26/10/2012 05/11/2012 15/11/2012 25/11/2012 05/12/2012 15/12/2012 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 Ngày bón (ngày sau sạ) Số lƣợng phân N (kg/ha) 11 23 29 41 32 38 38 37 32 01/03/2013 18/03/2013 02/04/2013 18/04/2013 25 45 52 25 50 25 25/01/2013 13/02/2013 01/03/2013 18/03/2013 Ngày/tháng/năm 71 Phân bón Ngày Số bón Số lƣợng lƣợng (ngày phân N (mm) sau (kg/ha) sạ) 50 27 38 90 32 42 60 40 31 60 55 37 50 50 100 20 27 51 61 25 50 25 Hè Thu Phân bón Ngày bón (ngày sau sạ) Số lƣợng phân N (kg/ha) 20 25 33 48 59 32 42 31 37 26 14 25 35 55 12 25 38 25 Đông Xuân Tƣới Điểm Xuân Hè Tƣới Phân bón Ngày bón (ngày sau sạ) Số lƣợng phân N (kg/ha) Ngày/tháng/năm Số lƣợng (mm) D3 15/10/2012 21/10/2012 26/10/2012 2/11/2102 9/11/2012 16/11/2012 23/11/2012 30/11/2012 7/12/2012 14/12/2012 24/12/2012 50 50 50 50 50 50 50 50 50 50 50 10 20 45 80 23 28 28 23 23 D4 14/11/2012 29/11/2012 9/12/2012 19/12/2012 29/12/2012 50 50 50 50 50 D5 7/11/2012 14/11/2012 21/11/2012 28/11/2012 5/12/2012 50 50 50 50 50 18 33 48 63 78 19 34 58 18 26 23 29 20 36 32 32 30 30 Ngày/tháng/năm 22/01/2013 02/02/2013 24/02/2013 72 Phân bón Ngày Số bón Số lƣợng lƣợng (ngày phân N (mm) sau (kg/ha) sạ) 50 23 50 10 34 50 18 28 25 32 45 23 60 23 Hè Thu Phân bón Ngày bón (ngày sau sạ) Số lƣợng phân N (kg/ha) 25 55 32 32 25 35 60 70 25 37 37 25 22 35 45 37 45 52 20 Đông Xuân Tƣới Điểm D6 D7 D8 Xuân Hè Tƣới Phân bón Ngày/tháng/năm Số lƣợng (mm) 12/12/2012 19/12/2012 26/12/2012 3/1/2013 10/1/2013 17/01/2013 8/10/2012 28/10/2012 7/11/2012 27/11/2012 50 50 50 50 50 50 30 30 30 30 20/09/2012 19/10/2012 30/10/2012 9/11/2012 19/11/2012 6/10/2012 21/10/2012 26/10/2012 50 50 50 50 50 50 70 50 Ngày bón (ngày sau sạ) Số lƣợng phân N (kg/ha) 10 20 43 32 32 23 16 26 45 38 38 30 30 34 60 25 50 25 Ngày/tháng/năm 18/01/2013 28/01/2013 08/02/2013 18/02/2013 28/02/2013 08/03/2013 18/03/2013 28/03/2013 02/01/2013 17/01/2013 08/02/2013 73 Phân bón Ngày Số bón Số lƣợng lƣợng (ngày phân N (mm) sau (kg/ha) sạ) 30 30 30 30 30 30 30 30 50 50 50 Hè Thu Phân bón Ngày bón (ngày sau sạ) Số lƣợng phân N (kg/ha) 10 20 43 32 32 23 10 20 35 43 50 26 10 24 17 26 42 38 38 30 30 12 30 40 50 38 24 24 30 14 25 35 12 25 38 Đông Xuân Tƣới Điểm D9 Xuân Hè Tƣới Phân bón Ngày/tháng/năm Số lƣợng (mm) 10/11/2012 15/11/2012 30/11/2012 15/12/2012 23/10/2012 13/11/2012 23/11/2012 30/11/2012 10/12/2012 20/12/2012 30/12/2012 70 70 70 70 50 70 60 80 70 50 80 Ngày bón (ngày sau sạ) 11 20 27 32 72 Số lƣợng phân N (kg/ha) Ngày/tháng/năm 18 37 18 23 23 74 Phân bón Ngày Số bón Số lƣợng lƣợng (ngày phân N (mm) sau (kg/ha) sạ) Hè Thu Phân bón Ngày bón (ngày sau sạ) Số lƣợng phân N (kg/ha) 55 25 16 24 46 23 32 43 36 Phụ lục 5: Dữ liệu quản lý cấy trồng Vụ Điểm Đẻ nhánh Trổ Chín sáp Chín hoàn toàn Tình hình dịch hại Ốc bƣu vàng Ốc bƣu vàng Ốc bƣu vàng Ốc bƣu vàng, muối hạt, nhiều hạt lép, dễ rụng Ốc bƣu vàng Ốc bƣu vàng, chuột Ốc bƣu vàng Ốc bƣu vàng Ốc bƣu vàng, cỏ dại OM8576 OM108-5 OM10636 ST5 100 100 95 120 Cấy Cấy Sạ Sạ 25/10/2012 23/09/2012 28/09/2012 25/10/2012 28/10/2012 26/10/2012 28/10/2012 25/11/2012 20/12/2012 18/11/2012 30/11/2012 17/01/2013 11/1/2013 05/12/2012 20/12/2012 27/01/2013 21/01/2013 13/12/2012 28/12/2012 01/02/2013 OM8017 OM8595 OM6976 OM4900 95 95 100 100 100 Sạ Sạ Sạ Sạ Sạ 30/10/2012 28/09/2012 10/09/2012 26/09/2012 13/10/2012 01/12/2012 28/10/2012 10/10/2012 30/10/2012 25/10/2012 01/01/2013 25/11/2012 20/11/2012 6/12/2012 23/12/2012 20/01/2013 15/12/2012 10/12/2012 26/12/2012 13/01/2012 01/02/2013 20/12/2012 20/12/2012 05/01/2012 23/01/2012 D1 D2 OM6976 OM108-5 100 100 Cấy Cấy 01/03/2013 10/03/2013 25/04/2013 15/01/2013 30/01/2013 10/03/2013 D3 OM6976 100 Sạ 15/01/2013 15/02/2013 27/03/2013 D6 OM9584-1 95 Sạ 10/01/2013 10/02/2013 15/03/2013 D7 OM8017 95 Sạ 30/12/2012 28/01/2013 01/03/2013 D1 D2 D3 D4 Đông Xuân D5 D6 D7 D8 D9 Xuân Hè Giống Thời Hình gian thức Ngày sạ/cấy sinh gieo (ngày/tháng/năm) trƣởng trồng (ngày) OM6976-41 75 15/05/2013 20/05/2013 Ốc bƣu vàng 30/03/2013 15/04/2013 Ốc bƣu vàng; sâu lá; nhiễm mặn đầu vụ, cuối vụ 18/04/2013 28/04/2013 Ốc bƣu vàng, sâu lá, cỏ dại 05/04/2013 13/04/2013 Ốc bƣu vàng, chuột, rầy nâu, nhiễm phèn mặn cuối vụ 20/03/2013 25/03/2013 Ốc bƣu vàng, nhiễm phèn mặn cuối vụ Vụ Hè Thu Điểm Giống Thời Hình gian thức Ngày sạ/cấy sinh gieo (ngày/tháng/năm) trƣởng trồng (ngày) D1 D2 D3 OM6976 OM108-5 OM3673 100 100 95 Cấy Cấy Cấy D4 ST5 120 Cấy D5 ST5 120 D6 D7 D8 D9 OM9584-1 95 95 100 100 OM3673 OM6976 OM4900 14/06/2013 26/05/2013 25/06/2013 09/06/2013 10/05/2013 Đẻ nhánh Trổ Chín sáp Chín hoàn toàn 29/07/2013 09/08/2013 10/09/2103 18/09/2013 01/07/2013 03/08/2013 20/08/2013 24/08/2013 09/07/2013 03/08/2013 27/08/2013 03/09/2013 30/05/2013 13/08/2013 22/08/2013 09/09/2013 Sạ 12/06/2013 02/07/2013 10/09/2013 28/09/2013 05/10/2013 Sạ Sạ Sạ Sạ 02/06/2013 29/06/2013 19/06/2013 16/06/2013 26/07/2013 05/07/2013 25/07/2013 26/08/2013 30/08/2013 10/09/2013 12/10/2013 04/08/2013 03/08/2013 12/08/2013 15/09/2013 76 03/09/2013 09/09/2013 18/09/2013 16/10/2013 Tình hình dịch hại Ốc, sâu Ốc Ốc, Sâu Ốc, chuột; Lá lúa có nhiều đốm đỏ; Lúa có nhiều hạt lép sâu lá; lúa bị vàng, có đốm đỏ Sâu Sâu lá, đạo ôn, lem lép hạt [...]... thực tế của vùng trồng lúa Sóc Trăng trong việc áp dụng công nghệ vào việc quản lý và dự báo sản lƣợng lúa để đƣa ra các chính sách phù hợp Từ đó, đề tài Ứng dụng mô hình Oryza2000 kết hợp ảnh viễn thám MODIS thành lập bản đồ năng suất lúa tỉnh Sóc Trăng đƣợc thực hiện 1.2 Mục tiêu và phạm vi nghiên cứu 1.2.1 Mục tiêu tổng quát Sử dụng mô hình Oryza2000 kết hợp với công nghệ Viễn thám và GIS để giám... giám sát quá trình sinh trƣởng, phát triển của cây lúa qua đó dự báo năng suất và sản lƣợng lúa khu vực tỉnh Sóc Trăng nhằm đảm bảo vấn đề an ninh lƣơng thực 1.2.2 Mục tiêu cụ thể - Thành lập bản đồ hiện trạng cơ cấu mùa vụ tỉnh Sóc Trăng; - Kết hợp mô hình Oryza2000 và ảnh viễn thám MODIS thành lập bản đồ phân bố năng suất lúa; 1 - Dự báo đƣợc sản lƣợng lúa làm cơ sở hoạch định chính sách an ninh lƣơng... đầu của mô hình Oryza2000 46 Hình 4.21 Biểu đồ tƣơng quan năng suất lúa 49 Hình 4.22 Bảng tính năng suất thực tế 50 Hình 4.23 Bản đồ diện tích gieo trồng vụ Đông Xuân 51 Hình 4.24 Bản đồ phân bố năng suất vụ Xuân Hè 54 Hình 4.25 Bản đồ diện tích gieo trồng và phân bố năng suất vụ Hè Thu 55 xii Hình 4.26 Đánh giá độ tin cậy trên ENVI 57 Hình 4.27 Biểu đồ so sánh diện tích thực tế và giải đoán 58 Hình 4.28... nghiên cứu 21 Bảng 4.1 Bảng tiêu chuẩn sử dụng để phân loại cơ cấu sử dụng đất 31 Bảng 4.2 Lịch thời vụ lúa 1 vụ 38 Bảng 4.3 Lịch thời vụ lúa 2 vụ 39 Bảng 4.4 Lịch thời vụ lúa 3 vụ 40 Bảng 4.5 Kết quả mô phỏng năng suất lúa của các điểm (kg/ha) 41 Bảng 4.6 Kết quả năng suất mô phỏng và năng suất lúa thực tế (kg/ha) 42 Bảng 4.7 So sánh sản lƣợng ƣớc đoán và sản lƣợng thống kê vụ Đông Xuân 45 Bảng 4.8 So... hệ thống FSE 7 Hình 2.2 Các thành phần cơ bản của hệ thống viễn thám 11 Hình 2.3 Phổ điện từ 11 Hình 2.4 Các giai đoạn phát triển của cây lúa 14 Hình 2.5 Tƣơng quan giữa NDVI và sự phát triển của lúa ở vụ Đông xuân và Hè Thu 15 Hình 2.6 Bản đồ hành chính tỉnh Sóc Trăng 18 Hình 3.1 Trang chủ Dowload ảnh MODIS 21 Hình 3.2 Chọn vùng cần load ảnh và thời gian chụp ảnh 22 Hình 3.3 File ảnh MODIS cho phép... 22 Hình 3.4 Vị trí các điểm thu mẫu tại tỉnh Sóc Trăng Hình 3.5 Máy đo chỉ số diện tích lá (Ceptometer LP-80) Hình 3.6 Diện tích khu vực lấy mẫu thu hoạch 26 Hình 3.7 Dụng cụ thu mẫu năng suất 26 Hình 3.8 Các bƣớc xây dựng bản đồ phân bố năng suất 27 Hình 4.1 Ảnh Modis ĐBSCL trƣớc và sau khi ghép lại 29 Hình 4.2 Ảnh đƣợc hiệu chỉnh về đúng toạ độ 30 Hình 4.3 Ảnh trƣớc và sau khi che 31 Hình 4.4 Ảnh. .. chứng minh tính ƣu việt của công nghệ viễn thám trong dự báo năng suất, sản lƣợng lúa ở nƣớc ta Phạm Thị Lệ Huyền (2013), Nghiên cứu này thực hiện nhằm mô phỏng năng suất lúa ở Sóc Trăng, kết hợp với quan sát và so sánh với thực tế Đánh giá 16 khả năng ứng dụng của mô hình này tại Đồng bằng sông Cửu Long Mô hình đƣợc thực hiện với 3 vụ sản xuất: Đông Xuân, Xuân Hè và Hè Thu tại 9 điểm canh tác lúa Kết. .. 38 Hình 4.13 Bản đồ vùng canh tác lúa 38 Hình 4.14 Vùng không trồng lúa 39 Hình 4.15 Bản đồ hiện trạng cơ cấu mùa vụ Sóc Trăng 2013 39 Hình 4.16 Vùng trồng lúa vụ Xuân Hè 40 Hình 4.17 Biểu đồ thể hiện sự thay đổi giá trị NDVI vùng trồng lúa 1 vụ 41 Hình 4.18 Biểu đồ thể hiện sự thay đổi giá trị NDVI vùng trồng lúa 2 vụ 42 Hình 4.19 Biểu đồ thể hiện sự thay đổi giá trị NDVI vùng trồng lúa 3 vụ 43 Hình. .. cứu Truy tìm sự khởi đầu của nó từ một mô hình đơn giản cho sự tăng trƣởng và sản xuất, về sau đã trở thành một công cụ toàn diện mô hình lúa đƣợc áp dụng để phân tích kịch bản khác nhau Mô hình có khả năng mạnh mẽ về ƣớc đoán hạn chế thời tiết, tăng trƣởng và năng suất lúa - tiềm năng tăng trƣởng và năng suất Nó có khả năng tốt về ƣớc lƣợng sự tăng trƣởng và năng suất cây trồng trong điều kiện thực tế... còn xây dựng đƣợc bản đồ thay đổi hiện trạng cơ cấu mùa vụ mùa vụ cho cả ĐBSCL năm 2011, xây dựng quy trình xử lý ảnh Viễn thám MODIS để thành lập bản đồ cơ cấu mùa vụ cho cấp vùng và bản đồ cơ cấu mùa vụ cho cả nƣớc Nghiên cứu của C Wikarmpapraharn and E Kositsakulchai (2010) báo cáo kết quả so sánh hai mô hình tăng trƣởng cây lúa CERES và ORYZA2000 ở đồng bằng miền Trung Thái Lan ORYZA2000 là công ... trồng lúa Sóc Trăng việc áp dụng công nghệ vào việc quản lý dự báo sản lƣợng lúa để đƣa sách phù hợp Từ đó, đề tài Ứng dụng mô hình Oryza2000 kết hợp ảnh viễn thám MODIS thành lập đồ suất lúa tỉnh. .. ĐẠI HỌC CẦN THƠ KHOA MÔI TRƯỜNG VÀ TÀI NGUYÊN THIÊN NHIÊN -- NGÔ HỮU LỢI ỨNG DỤNG MÔ HÌNH ORYZA2000 KẾT HỢP ẢNH VIỄN THÁM MODIS THÀNH LẬP BẢN ĐỒ NĂNG SUẤT LÚA TỈNH SÓC TRĂNG LUẬN VĂN TỐT NGHIỆP... Năm 2015 CHẤP THUẬN LUẬN VĂN CỦA HỘI ĐỒNG Luận văn này, với đề tựa là Ứng dụng mô hình Oryza2000 kết hợp ảnh viễn thám MODIS thành lập đồ suất lúa tỉnh Sóc Trăng học viên Ngô Hữu Lợi thực theo

Ngày đăng: 21/12/2015, 01:06

Từ khóa liên quan

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

Tài liệu liên quan