Sử dụng tư liệu vệ tinh MODIS đánh giá chất lượng không khí khu vực đồng bằng sông hồng

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN KHOA ĐỊA LÝ PHẠM THỊ DUYÊN SỬ DỤNG TƢ LIỆU VỆ TINH MODIS ĐÁNH GIÁ CHẤT LƢỢNG KHÔNG KHÍ KHU VỰC ĐỒNG BẰNG SÔNG HỒNG LUẬN VĂN THẠC SĨ KHOA HỌC Hà Nội – 2016 ĐẠI HỌC QUỐC GIA HÀ NỘI TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN KHOA ĐỊA LÝ PHẠM THỊ DUYÊN SỬ DỤNG TƢ LIỆU VỆ TINH MODIS ĐÁNH GIÁ CHẤT LƢỢNG KHÔNG KHÍ KHU VỰC ĐỒNG BẰNG SÔNG HỒNG Chuyên ngành: Bản đồ, Viễn thám Hệ thông tin địa lý Mã số: 66 04 02 14 LUẬN VĂN THẠC SĨ KHOA HỌC HƢỚNG DẪN KHOA HỌC CHỦ TỊCH HỘI ĐỒNG TS Bùi Quang Thành PGS.TS Đinh Thị Bảo Hoa Hà Nội – 2016 LỜI CẢM ƠN Lời cho phép đƣợc bày tỏ lòng biết ơn sâu sắc tới TS Bùi Quang Thành, ngƣời trực tiếp hƣớng dẫn, định hƣớng nhiệt tình giúp đỡ, hỗ trợ thực hiện, hoàn thành luận văn thạc sĩ Tôi xin trân trọng cảm ơn thầy cô giáo giảng dạy, truyền đạt giúp nâng cao kiến thức chuyên ngành thời gian học tập khoa Địa lý, trƣờng đại học Khoa học Tự nhiên, ĐHQG Hà Nội, đặc biệt đào tạo, giúp đỡ thầy cô môn Bản đồ - Viễn thám GIS Qua đây, xin gửi lời cảm ơn tới bạn lớp cao học K13 – Bản đồ, Viễn thám GIS ủng hộ tạo không khí vui vẻ, thoải mái trình học tập nhƣ trình làm luận văn Tôi xin chân thành cảm ơn Dự án “Nghiên cứu thủy tai biến đổi khí hậu xây dựng hệ thống thông tin nhiều bên tham gia nhằm giảm thiểu tính dễ bị tổn thương Bắc Trung Bộ Việt Nam (CPIS)” (2012-2015) đề tài VT-UD.06/16-20(thuộc Chƣơng trình Khoa học Công nghệ độc lập cấp Nhà nƣớc Công nghệ Vũ trụ) cung cấp nguồn liệu công cụ cần thiết để thực luận văn Lời cuối cùng, xin cảm ơn động viên, ủng hộ tinh thần gia đình, bạn bè đồng nghiệp tạo động lực giúp hoàn thành luận văn Hà Nội, ngày… tháng……năm 2016 Học viên cao học Phạm Thị Duyên MỤC LỤC DANH MỤC BẢNG BIỂU DANH MỤC HÌNH VẼ DANH MỤC CÁC TỪ VIẾT TẮT MỞ ĐẦU Tính cấp thiết Mục tiêu nghiên cứu: Nhiệm vụ nghiên cứu: Phƣơng pháp nghiên cứu Phạm vi nghiên cứu Ý nghĩa khoa học thực tiễn Cấu trúc luận văn CHƢƠNG 11 CƠ SỞ LÝ LUẬN VỀ PHÁT HIỆN Ô NHIỄM KHÔNG KHÍ BẰNG ĐỘ DÀY QUANG HỌC SOL KHÍ 11 1.1 Một số khái niệm 11 1.1.1 Khái niệm sol khí 11 1.1.2 Khái niệm hạt PM 12 1.2 Tác động bụi sol khí sức khỏe cộng đồng 13 1.3 Tổng quan phát ô nhiễm không khí độ dày quang học sol khí 15 CHƢƠNG PHƢƠNG PHÁP ĐÁNH GIÁ CHẤT LƢỢNG KHÔNG KHÍ 18 KHU VỰC ĐỒNG BẰNG SÔNG HỒNG BẰNG TƢ LIỆU VỆ TINH MODIS 18 2.1 Cơ sở khoa học phƣơng pháp đánh giá chất lƣợng không khí khu vực ĐBSH tƣ liệu vệ tinh 18 2.1.1 Phƣơngpháp ƣớctínhsolkhí từảnhvệtinh 18 Trang1 2.1.2 Phƣơngpháptínhtoánhồiquy 20 2.2 Khu vực nghiên cứu 23 2.3 Dữ liệu nghiên cứu 26 2.3.1 Dữ liệu ảnh sol khí MODIS 26 2.3.2 Dữ liệu PM 2.5 từ trạm quan trắc mặt đất 30 2.4 Tính toán giá trị PM 2.5 từ liệu MODIS AOD/AOT 33 2.4.1 Quy trình tổng quát 33 2.4.2 Quy trình tích hợp liệu vệ tinh quan trắc 35 2.4.3 Quy trình tính toán thành lập đồ PM 2.5 từ ảnh sol khí MODIS 36 2.5 Ứng dụng web-based GIS phục vụ giám sát bụi PM 2.5 36 CHƢƠNG 40 KẾT QUẢ NGHIÊN CỨU VÀ THẢO LUẬN 40 3.1 Thực nghiệm 40 3.1.1 Dữ liệu sol khí 40 3.1.2 Dữ liệu thu đƣợc từ trạm quan trắc 43 3.1.3 Xử lý hình học ảnh 46 3.1.4 Tính toán hồi quy AOD/AOT PM 2.5 46 3.2 Thành lập đồ số PM2.5 cho khu vực đồng sông Hồng 54 3.3 WebGis thể đồ PM 2.5 56 KẾT LUẬN VÀ KIẾN NGHỊ 58 TÀI LIỆU THAM KHẢO 60 PHỤ LỤC 62 Trang2 DANH MỤC BẢNG BIỂU Bảng 1.1 Tác động sol khí bụi đến sức khỏe ngƣời kết hợp với kim loại nặng 14 Bảng Các bệnh có tỉ lệ ngƣời mắc bệnh cao phạm vi toàn Việt Nam (Nguồn: Niên giám thống kê y tế năm 2008) 15 Bảng 2.1 Mức độ tƣơng quan thể qua R² 23 Bảng 2.2 Bƣớc sóng phổ kênh ảnh ứng dụng thực tiễn 26 Bảng 2.3 Ví dụ cấu trúc liệu PM 2.5 đƣợc lƣu trữ file exel 33 Bảng 3.1 Danh mục cảnh ảnh Terra 40 Bảng 3.2 Danh mục cảnh ảnh Aqua 41 Bảng 3.3 Danh mục trạm thực địa 43 Bảng 3.4 Tƣơng quan AOD PM2.5 trạm 53 Trang3 DANH MỤC HÌNH VẼ Hình 1.1 Kích thƣớc loại bụi PM so với đƣờng kính sợi tóc ngƣời 12 Hình 2.1 Các tia phản xạ khí đến đầu thu 19 Hình 2.2 Ví dụ mối liên hệ tuyến tính 22 Hình 2.3 Bản đồ vùng đồng sông Hồng 24 Hình 2.4 Cấu trúc nhóm band ảnh raster file HDF 29 Hình 2.5 Trạm quan trắc đặt Long Biên, Hà Nội 31 Hình 2.6 Đồ thị số liệu quan trắc từ 07 trạm số AIQ website Trung tâm quan trắc môi trƣờng 32 Hình 2.7 Quy trình tính toán chất lƣợng không khí khu vực ĐBSH 34 Hình 2.8 Quy trình tích hợp liệu vệ tinh liệu quan trắc 35 Hình 2.9 Quy trình tính toán PM 2.5 từ ảnh AOD/AOT 36 Hình 2.10 Thiết kế kỹ thuật Web-based GIS 37 Hình 2.11 Thiết kế chi tiết Web-based GIS 38 Hình 2.12 Giao diện WebGis thuộc CPIS 39 Hình 3.1 Ví dụ tƣ liệu thu thập từ ảnh số hiệu MOBAOD.A2014002.0345.006.2015072135002 từ vệ tinh Terra 41 Hình 3.2 Ví dụ tƣ liệu thu thập từ ảnh số hiệu MYBAOD.A2014001.0610.006.2014273233747từ vệ tinh Aqua 43 Hình 3.3 Vị trí 06 trạm mặt đất đồ GIS 44 Hình 3.4 Vị trí 06 trạm ảnh Modis Terra 45 Hình 3.5 Vị trí 06 trạm ảnh Modis Aqua 45 Hình 3.6 Giao diện công cụ MRT 46 Hình 3.7 Tƣơng quan AOD PM2.5 trạm Nguyễn Văn Cừ 47 Hình 3.8 Sai số mô hình hồi quy AOD PM2.5 trạm Nguyễn Văn Cừ 47 Hình 3.9 Tƣơng quan AOD PM2.5 trạm Phú Thọ 48 Hình 3.10 Sai số mô hình hồi quy AOD PM2.5 trạm Phú Thọ 48 Hình 3.11 Tƣơng quan AOD PM2.5 trạm Hạ Long 49 Hình 3.12 Sai số mô hình hồi quy AOD PM2.5 trạm Hạ Long 49 Hình 3.13 Tƣơng quan AOD PM2.5 trạm Huế 50 Hình 3.14 Sai số mô hình hồi quy AOD PM2.5 trạm Huế 50 Trang4 Hình 3.15 Tƣơng quan AOD PM2.5 trạm Đà Nẵng 51 Hình 3.16 Sai số mô hình hồi quy AOD PM2.5 trạm Đà Nẵng 51 Hình 3.17 Tƣơng quan AOD PM2.5 trạm Khánh Hòa 52 Hình 3.18 Sai số mô hình hồi quy AOD PM2.5 trạm Khánh Hòa 52 Hình 3.19 Tƣơng quan AOD PM2.5 06 trạm 53 Hình 3.20 Sai số mô hình hồi quy AOD PM2.5 tạị 06 trạm 53 Hình 3.21 Biểu đồ hệ số xác định bội R² theo trạm 54 Hình 3.22 Bản đồ giá trị PM 2.5 khu vực đồng sông Hồng ngày 20/12/2016 55 Hình 3.23 Giao diện kết nối liệu raster Geosever 56 Hình 3.24 Tạo style cho lớp đồ PM 2.5 với SLD 57 Hình 3.25 Bản đồ PM 2.5 WebGis 57 Trang5 DANH MỤC CÁC TỪ VIẾT TẮT MODIS Một loại cảm biến có độ phân giải trung bình Moderate resolution Imaging Spectrometer AQUA Một loại vệ tinh mang theo cảm biến MODIS TERRA Một loại vệ tinh mang theo cảm biến MODIS MOD Tiếp đầu ngữ file ảnh MODIS thu nhận từ vệ tinh TERRA MYD Tiếp đầu ngữ file ảnh MODIS thu nhận từ vệ tinh AQUA AOD/AOT Độ sâu quang học sol khí AQI Chỉ số chất lƣợng không khí GIS Hệ thông tin địa lý NASA Cơ quan hàng không Vũ trụ Mỹ Bộ TN&MT Bộ Tài nguyên & Môi trƣờng ĐBSH Đồng sông Hồng CENMA Trung tâm Quan trắc Phân tích Tài nguyên Môi trƣờng Hà Nội DN Digital number- Giá trị số R, r Correlation Coefficient – Hệ số tƣơng quan R², r² Coeffocient of Determination- Hệ số xác định bội PM Paticulate Matter – Bụi Trang6 MỞ ĐẦU Tính cấp thiết Môi trƣờng tất bao quanh ngƣời có ảnh hƣởng trực tiếp gián tiếp tới tồn phát triển ngƣời, không khí năm thành phần môi trƣờng Do vậy, chất lƣợng môi trƣờng không khí vấn đề đáng đƣợc quan tâm Ngày nay, với phát triển không ngừng xã hội, môi trƣờng nói chung môi trƣờng không khí nói riêng bị ô nhiễm nặng nề tàn phá ngƣời.Tổ chức Y tế Thế Giới (WHO) gần đƣa báo cáo ô nhiễm không khí năm vừa qua (năm 2014), dựa số liệu mức độ ô nhiễm 1.600 thành phố khắp 19 quốc gia Tổ chức sử dụng hệ thống đánh giá có tên PM2.5 PM10 PM2.5 đƣợc coi hệ thống tốt đƣợc dùng để đánh giá tác động ô nhiễm không khí lên sức khỏe xác định nồng độ bụi ô nhiễm có đƣờng kính từ 2,5 micromet trở xuống Những hạt bụi ô nhiễm khói, bụi bẩn, nấm mốc phấn hoa, đƣợc tổng hợp từ kim loại nặng hợp chất hữu độc hại Đây đƣợc coi mối nguy hiểm lớn cho thể ngƣời bị “tích lũy” hệ thống hô hấp Việc giám sát chất lƣợng không khí thu hút đƣợc quan tâm tất quốc gia giới, có Việt Nam Tại Việt Nam, trạm quan trắc khí tƣợng đo lƣờng số liên quan chất lƣợng không khí (cũng nhƣ số PM2.5), nhiên, số lƣợng trạm hạn chế Do vậy, việc tính toán, đánh giá chất lƣợng không khí diện rộng, toàn khu vực thực đƣợc Cùng với phát triển công nghệ vệ tinh quan sát Trái đất, khả ứng dụng tƣ liệu ảnh vệ tinh kết hợp với hệ thông tin địa lý giám sát chất lƣợng không khí đƣợc áp dụng nhiều quốc gia giới (khoảng 40 năm trở lại đây) Ảnh vệtinh cho biết thông tin toàn khu vực nghiên cứu theo cấu trúc mạng lƣới liên tục pixel kềnhau Mỗi pixel đƣợc ví nhƣlà trạm quan trắc, sốđiểm từảnh vệtinh sẽdày đặc nhiều so với trạm đo mặt đất Kết quảphân tích ảnh vệtinh sẽcho giá trịô nhiễm thểhiện pixel tùy thuộc vào độphân giải ảnh toàn vùng đồng thời vào thời điểm quan sát, mà với điều kiện khảnăng thiết bịvà trạm đo mặt đất không thểnào đạt đƣợc Trang7 Hình 19- Tương quan AOD PM2.5 06 trạm Hình 20- Sai số mô hình hồi quy AOD PM2.5 tạị 06 trạm Tổng hợp kết độ tƣơng quan liệu trạm mặt đất tƣ liệu MODIS đƣợc thể bảng dƣới đây: Bảng 4- Tương quan AOD PM2.5 trạm Số TT Trạm R² Hà Nội 0,575 Phú Thọ 0,156 Hạ Long 0,548 Huế 0,610 Đà Nẵng 0,578 Khánh Hòa 0,353 06 trạm 0,377 Trang53 0.7 0.6 0.5 0.4 0.3 0.2 0.1 Hà Nội Phú Thọ Hạ Long Huế Đà Nẵng Khánh Hòa 06 trạm Hình 21-Biểu đồ hệ số xác định bội R² theo trạm Kết luận: Dựa giá trị nồng độ bụi PM 2.5 đo đƣợc từ trạm quan trắc mặt đất, ta thấy đƣợc nồng độ PM 2.5 vùng nông thôn thành thị, tháng năm tƣơng đối khác biệt Ở khu vực đƣợc đƣa vào thực nghiệm, khí hậu có phân mùa rõ rệt (đặc biệt khu vực ĐBSH), nồng độ PM 2.5 mùa nóng (vào tầm tháng đến tháng 8) thấp hẳn so với thời điểm mùa lạnh (tầm tháng 1, tháng 2), mùa hè lƣợng mƣa nhiều kéo theo hạt bụi đƣợc rửa trôi Với kết chạy mô hình nhƣ trên,hệ số xác định bội R2 (Determination Coefficiant) dao động khoảng từ 0,156(tại Phú Thọ) tới 0,610 (tại Huế) Khu vực Phú Thọ có biến động giá trị PM2.5 lớn hẳn so với khu vực khác kéo theo thông số ảnh dao động lớn, liệu phân tán mạnh, R2có giá trị thấp.Tại khu vực khác, độ tƣơng quan liệu tƣơng đối tốt, 0,5 Nhìn chung, đại lƣợng R2 chƣa cao so với 1, bởi: - Số lƣợng điểm sử dụng để lập hàm hồi quy mỏng; - Sự phân bố điểm quan trắc mặt đất chƣa trải cách tối ƣu - Tại số điểm, giá trị PM2.5 đo trạm quan trắc dao động lớn 3.2 Thành lập đồ số PM2.5 cho khu vực đồng sông Hồng Dựa hàm hồi quy tính dựa 06 trạm, phƣơng trình hồi quy có dạng y = 22,8038 * x - 6,6905 (3.1) đó: y – giá trị PM2.5; x – giá trị AOD/AOT trích xuất từ ảnh Trang54 Với R2 = 0,377, giá trị số xác định bội không cao, nhiên so với nghiên cứu trƣớc đây, phƣơng trình đƣợc sử dụng để tính toán số PM2.5 phục vụ cho việc định hƣớng khu vực bị ô nhiễm hàm lƣợng vật chất lơ lửng không khí Từ phƣơng trình hồi quyvà từ ảnh MODIS tải vào ngày 20 /12/2016, thành lập đƣợc đồ số PM 2.5 khu vực đồng sông Hồng nhƣ sau: Hình 22- Bản đồ giá trị PM2.5 khu vực đồng sông Hồng ngày 20/12/2016 Trang55 Giá trị nồng độ PM 2.5 đƣợc phân ngƣỡng theo chuẩn quy định Bộ Tài nguyên & Môi trƣờng – ngƣỡng có giá trị 50 µg/m³ Kết cho thấy khu vực ĐBSH đồng ven biển có giá trị nồng độ PM 2.5 cao hẳn so với khu vực vùng núi cao 3.3 WebGis thể đồ PM 2.5 Sau tính toán đƣợc giá trị PM 2.5 thành lập đƣợc đồ PM 2.5 toàn khu vực ĐBSH, kết đƣợc đƣa lên WebGis (CPIS) để ngƣời dùng quan tâm nắm bắt đƣợc thông tin giá trị PM 2.5 cách trực quan chi tiết Để thực đƣợc điều này, cần tiến hành kết nối liệu PM 2.5 với CSDL hệ thống Nhƣ mô hình WebGis trình bày mục 2.4, CSDL PostGis Tuy nhiên, thực tế, liệu để chạy hệ thống WebGis CPIS chƣa nhiều chƣa đáp ứng đủ nguồn nhân lực vận hành hệ thống (một cách bản, đảm bảo quy trình) nên việc kết nối lƣu trữ liệu vào PostgreSQL/ PostGis tạm thời đƣợc bỏ qua Nhƣ vậy, ảnh PM 2.5 đƣợc kết nối trực tiếp với GeoSever nhƣ hình dƣới: Hình 23-Giao diện kết nối liệu raster Geosever Sau liệu đƣợc thêm vào GeoSever, tiến hành trình bày đồ PM 2.5 với SLD (Styled Layer Desciptor) Trang56 Hình 24- Tạo style cho lớp đồ PM 2.5 với SLD Kết đƣa đồ PM 2.5 lên WebGis đƣợc minh họa hình 37 Ngƣời dùng sử dụng chức “Truy vấn không gian” để xem thông tin giá trị PM 2.5 vị trí đồ Hình 25- Bản đồ PM 2.5 WebGis Trang57 KẾT LUẬN VÀ KIẾN NGHỊ Kết luận Qua phần lý thuyết thực nghiệm, đề tài chứng minh khả ứng dụng công nghệ viễn thám vào lĩnh vực giám sát ô nhiễm không khí Việt Nam nói chung, việc sử dụng tƣ liệu MODIS vào đánh giá chất lƣợng không khí nói riêng.Phƣơng pháp tiến hành cho kết thể tranh toàn cảnh tình trạng ô nhiễm không khí khu vực ĐBSH nhƣ toàn lãnh thổ Việt Nam Mặc dù số lƣợng điểm quan trắc thực địa mỏng không trải phạm vi nghiên cứu nhƣng kết cho thấy số R2 xác định tƣơng đối khả quan, bƣớc đầu thể tính đắn phƣơng pháp ƣớc tính.So sánh với nghiên cứu tƣơng tự khu vực giới, kết thu đƣợc tƣơng đƣơng Trong phạm vi đề tài, kết chấp nhận đƣợc nghiên cứu cải tiến độ xác mô hình tính toán tƣơng lai Kết thực nghiệm cho thấy giá trị nồng độ PM 2.5 khu vực ĐBSH đồng ven biển cao hẳn so với khu vực vùng núi, đặc biệt cao khu vực Hà Nội Thực tế, khu có nồng độ PM 2.5 cao khu vực có tốc độ đô thị hóa mạnh, tập trung nhiều khu công nghiệp, làng nghề, hoạt động xây dựng giao thông dày đặc, Việc đƣa kết đồ bụi PM 2.5 lên webgis hoàn toàn phù hợp với xu này, phƣơng pháp truyền tải thông tin đến cộng đồng cách hiệu quả, trực quan hữu hiệu Với thực trạng quan trắc môi trƣờng không khí tự động Việt Nam (số lƣợng trạm quan trắc mặt đất ít, nhiều trạm ngừng hoạt động, việc đầu tƣ xây dựng mạng lƣới trạm quan trắc dày đặc gặp nhiều khó khăn chi phí lớn), phƣơng pháp nghiên cứu đề tài góp phần bổ sung phƣơng pháp quan trắc phục vụ việc đánh giá tình hình ô nhiễm bụi khu vực đồng sông Hồng Kết thực nghiệm chứng minh, phƣơng pháp áp dụng cho toàn khu vực Việt Nam Kiến nghị Ngoài phƣơng pháp sử dụng ảnh sol khí MODIS để tính toàn nồng độ PM 2.5, tiến hành nghiên cứu thêm phƣơng pháp khác để tính toán hàm lƣợng bụi Trang58 Nghiên cứu sử dụng sản phẩm ảnh MODIS độ phân giải 10 km nên giới hạn độ chi tiết tƣơng lai sử dụng tƣ liệu ảnh có độ phân giải cao để có đƣợc đồ bụi PM 2.5 chi tiết Để nâng cao độ xác mô hình tính toán, cần tăng dày liệu quan trắc thực địa Trang59 TÀI LIỆU THAM KHẢO http://doc.edu.vn/tai-lieu/luan-van-sol-khi-va-mo-hinh-regcm-51575/ http://www.mfg.mtu.edu/cyberman/environment/air/pmintro.html http://editors.eol.org/eoearth/wiki/Impact_of_local_air_pollution Jimoda, L.A., EFFECTS OF PARTICULATE MATTER ON HUMAN HEALTH, THE ECOSYSTEM, CLIMATE AND MATERIALS: A REVIEW FACTA UNIVERSITATIS 2012 9: p 17 Kaufman, Y.J., Fraser, R.S., Ferrare, R.A, Satellite measurements of large- scale air pollution Journal of Geophysics Research, 1990 95: p 14 Hadjimitsis, D.G., Clayton, C.R.I, Determination of aerosol optical thickness through the derivation of an atmospheric correction for shortwave length Landsat TM and ASTER image data: an application to areas located in the vicinity of airports at UK and Cyprus Applied Geomatics 10 11 12 13 14 15 16 Journal, 2009 1: p Jun Wang, S.A.C., Intercomparison between satellite-derived aerosol optical thickness and PM2.5 mass GEOPHYSICAL RESEARCH LETTERS, 2003 30 Jill A Engel-Co, C.H.H., Basil W Coutant, Raymond M Hoff, Qualitative and quantitative evaluation of MODIS satellite sensor data for regional and urban scale air quality Atmospheric Environment 2004 38: p 14 Phonekeo, V., Application of Aerosol Optical Thickness Product (MOD04) to Air Quality Monitoring in Thailand Bangkok Metropolitan, 2009 Luong Chinh Ke, H.T.V.T., Tran Ngoc Tuong, Lam Dao Nguyen, Detecting air pollution in VietNam by optical satellite images The 31st Asian Conference on Remote Sensing, 2010 Trần Thị Vân, N.P.K., Hà Dƣơng Xuân Bảo, Viễn thám độ dày quang học mô phân bố bụi PM10 khu vực nội thành TPHCM Tạp chí Khoa học ĐHQGHN: Các Khoa học Trái đất Môi trƣờng, 2014 30 N I SIFAKIS, N.A.S., D K PARONIS, Quantitative mapping of air pollution density using Earth observation: a new processing method and application to an urban area Inter Jour of Remote Sensing, 1998 19(17): p 11 https://websrv1.ctu.edu.vn/coursewares/kinhte/phantichdulieu/chuong6.htm https://websrv1.ctu.edu.vn/coursewares/congnghe2/pv_toan/chuong5.htm https://mathbits.com/MathBits/TISection/Statistics2/correlation.htm http://ptit.edu.vn/wps/portal/nongthonvn/!ut/p/c4/04_SB8K8xLLM9MS SzPy8xBz9CP0os3hLizBHd1cfIwN_MyM3A08vc2cXVx83Y49AY_2CbEdFA Trang60 O8ydjg!/?WCM_GLOBAL_CONTEXT=/wps/wcm/connect/nongthonvn/nongt 17 honvn/vungnongthon/dongbangsonghong/af2fae80404cbe13a81bfa9171cb7767 Nguyễn Đức Ngữ, N.T.H., Khí hậu tài nguyên khí hậu Việt Nam NXB Nông 18 nghiệp, 2004 http://tcmt.monre.gov.vn/pages/article.aspx?item=X%C3%A2yd%E1%BB%B1ng, -.h.C.A.n.-t.E.B.n.-h.E.B.-t 19 http://www.quantracmoitruong.gov.vn/Portals/0/Images/news/2014/Than g%208/Bai%20trinh%20bay%20van%20hanh%20Tram%20Khi.pdf 20 Phạm Xuân Trƣờng, N.Q.M., Xây dựng ứng dụng Web-base Gis cho hệ thống 21 thông tin đất đai Tạp chí KHKT Mỏ-Địa chất, 2014 48: p http://danida.vnu.edu.vn:8080/congdong/composer/ 22 Bui Quang Thanh, N.Q.H., Nguyen Trung Kien, Pham Van Cu, Pham Minh and Hai, WEB-BASED PARTICIPATORY INFORMATION SYSTEM FOR VULNERABILITY MAPPING IN CENTRAL PROVINCES OF VIET NAM Trang61 PHỤ LỤC Phụ lục 1: Dữ liệu PM 2.5 đo trạm Đà Nẵng tƣơng ứng với liệu MODIS Aqua Terra đƣa vào thực nghiệm Time PM.10 PM.2.5 PM.1 1/1/2014 6:10 54.497 53.668 52.239 1/17/2014 6:10 44.5167 42.8075 41.3383 2/2/2014 6:10 21.3325 20.9892 20.4658 2/18/2014 6:10 26.335 26.27 26.0817 3/6/2014 6:10 24.7125 24.5025 23.6617 3/22/2014 6:10 18.2783 13.6817 10.53 4/7/2014 6:10 24.66 21.9067 19.9308 4/23/2014 6:10 25.0183 23.2875 21.1525 5/9/2014 6:10 40.955 21.0575 11.8375 5/25/2014 6:10 21.3792 17.1892 14.2275 6/10/2014 6:10 9.8617 7.6783 6.2975 6/26/2014 6:10 9.4325 7.8067 6.6117 7/12/2014 6:10 NA NA NA 7/28/2014 6:10 19.1242 15.555 12.9617 8/13/2014 6:10 12.2808 5.41 3.625 8/29/2014 6:10 38.1667 21.5033 16.1783 9/14/2014 6:10 29.6958 19.0942 13.7567 1/2/2014 3:45 64.531 63.908 62.943 1/18/2014 3:45 33.3217 32.2183 30.7258 2/3/2014 3:45 12.9433 12.85 12.5075 2/19/2014 3:45 1.175 0.83417 0.62417 3/7/2014 3:45 30.355 29.7092 27.8092 3/23/2014 3:45 12.38 9.4117 7.7183 4/8/2014 3:45 41.1817 34.7342 30.6183 4/24/2014 3:45 28.3542 25 21.5475 5/10/2014 3:45 30.0267 24.8467 21.3917 5/26/2014 3:45 19.6392 14.705 12.1142 6/11/2014 3:45 13.17 9.7983 8.0817 6/27/2014 3:45 11.7175 9.8242 8.3908 7/13/2014 3:45 NA NA NA 7/29/2014 3:45 18.5833 17.2492 15.3875 8/14/2014 3:45 14.2883 8.5117 6.61 8/30/2014 3:45 122.108 121.522 119.352 9/15/2014 3:45 20.9733 12.5875 8.4433 Trang62 Phụ lục 2: Dữ liệu PM 2.5 đo trạm Hạ Long tƣơng ứng với liệu MODIS Aqua Terra đƣa vào thực nghiệm Trang63 Phụ lục 3: Dữ liệu PM 2.5 đo trạm Huế tƣơng ứng với liệu MODIS Aqua Terra đƣa vào thực nghiệm Time 1/1/2014 6:10 1/17/2014 6:10 2/2/2014 6:10 2/18/2014 6:10 3/6/2014 6:10 3/22/2014 6:10 4/7/2014 6:10 4/23/2014 6:10 5/9/2014 6:10 5/25/2014 6:10 6/10/2014 6:10 6/26/2014 6:10 7/12/2014 6:10 7/28/2014 6:10 8/13/2014 6:10 8/29/2014 6:10 9/14/2014 6:10 1/2/2014 3:45 1/18/2014 3:45 2/3/2014 3:45 2/19/2014 3:45 3/7/2014 3:45 3/23/2014 3:45 4/8/2014 3:45 4/24/2014 3:45 5/10/2014 3:45 5/26/2014 3:45 6/11/2014 3:45 6/27/2014 3:45 7/13/2014 3:45 7/29/2014 3:45 8/14/2014 3:45 8/30/2014 3:45 9/15/2014 3:45 10/1/2014 3:45 PM.10 PM.2.5 PM.1 197.238 134.058 114.877 81.697 67.815 61.503 30.2583 22.5925 19.5442 66.303 47.2675 43.4017 37.6767 26.5792 22.8475 10.8083 9.6433 8.7433 51.129 36.6358 30.0667 67.046 44.805 34.815 45.3217 30.22 18.7317 63.177 28.6042 19.9633 46.0583 20.3608 15.6017 35.8617 12.0825 6.8208 46.375 15.6675 7.8083 25.63 17.09 11.2692 37.3583 10.1767 5.9042 24.4967 16.3367 8.2967 19.7433 13.1642 8.1133 86.199 77.948 75.576 84.777 67.242 58.638 51.138 35.1975 21.91 2.91333 2.61333 1.94667 42.6708 30.7942 25.1208 18.815 14.3533 12.7858 51.822 35.8267 28.4392 40.8275 23.0225 17.2217 27.5625 13.415 7.1242 31.8842 21.345 15.115 7.7158 4.16083 2.6925 21.7242 12.6283 6.1675 14.63 9.77 4.4575 15.6975 10.4733 2.11583 18.9517 7.1083 3.15833 34.0083 22.6733 13.715 26.2908 19.285 15.38 14.91 10.2725 6.8467 Trang64 Phụ lục 4: Dữ liệu PM 2.5 đo trạm Khánh Hòa tƣơng ứng với liệu MODIS Aqua Terra đƣa vào thực nghiệm Time 1/1/2014 6:10 1/17/2014 6:45 2/2/2014 0:00 2/18/2014 6:10 3/6/2014 6:10 3/22/2014 6:10 4/7/2014 6:10 4/23/2014 6:10 5/9/2014 6:10 5/25/2014 6:10 6/10/2014 6:10 6/26/2014 6:10 7/12/2014 6:10 7/28/2014 6:10 8/13/2014 6:10 8/29/2014 6:10 9/14/2014 6:10 9/30/2014 6:10 1/2/2014 3:45 1/18/2014 3:45 2/3/2014 3:45 2/19/2014 3:45 3/7/2014 3:45 3/23/2014 3:45 4/8/2014 3:45 4/24/2014 3:45 5/10/2014 3:45 5/26/2014 3:45 6/11/2014 3:45 6/27/2014 3:45 7/13/2014 3:45 7/29/2014 3:45 8/14/2014 3:45 8/30/2014 3:45 9/15/2014 3:45 10/1/2014 3:45 PM.10 PM.2.5 PM.1 41.2367 35.7608 31.0967 67.389 35.8044 23.8656 35.63 27.13 20.8133 43.3983 37.3158 33.2633 47.2083 37.4292 32.7308 38.8742 29.4408 19.685 59.842 43.4925 34.1242 40.1708 31.0083 25.1625 36.595 26.5683 20.2125 33.3 24.3575 18.365 15.9717 4.58083 2.75417 49.6333 23.605 13.7283 17.7233 12.0783 8.6325 24.6108 15 9.0442 41.5008 36.89 32.9133 28.4408 22.2942 17.2767 13.0267 12.3642 10.7783 28.0733 26.92 24.3308 69.625 64.028 58.028 26.36 19.1133 11.1883 32.2375 26.38 20.8442 39.8058 34.2475 30.0308 43.4442 38.4558 33.1808 20.1458 16.9408 11.8383 27.9067 22.1767 17.8533 41.0358 32.4833 24.9025 23.3225 17.6267 13.6075 18.1408 16.4017 14.1792 35.045 22.7058 16.7192 17.6917 14.7983 11.5342 10.9075 8.3242 6.6908 18.9617 15.1708 11.5392 28.0133 24.4925 20.4108 18.1267 11.0883 6.9533 8.025 7.555 6.7417 21.0783 16.8217 13.62 Trang65 Phụ lục 5: Dữ liệu PM 2.5 đo trạm Nguyễn Văn Cừ tƣơng ứng với liệu MODIS Aqua Terra đƣa vào thực nghiệm Time 1/1/2014 6:10 1/17/2014 6:10 2/2/2014 6:10 2/18/2014 6:10 3/6/2014 6:10 3/22/2014 6:10 4/7/2014 6:10 4/23/2014 6:10 5/9/2014 6:10 5/25/2014 6:10 6/10/2014 6:10 6/26/2014 6:10 7/12/2014 6:10 7/28/2014 6:10 8/13/2014 6:10 8/29/2014 6:10 9/14/2014 6:10 9/30/2014 6:10 1/2/2014 3:45 1/18/2014 3:45 2/3/2014 3:45 2/19/2014 3:45 3/7/2014 3:45 3/23/2014 3:45 4/8/2014 3:45 4/24/2014 3:45 5/10/2014 3:45 5/26/2014 3:45 6/11/2014 3:45 6/27/2014 3:45 7/13/2014 3:45 7/29/2014 3:45 8/14/2014 3:45 8/30/2014 3:45 9/15/2014 3:45 10/1/2014 3:45 PM.10 PM.2.5 PM.1 298.265 230.237 206.307 198.725 171.033 153.835 51.028 46.2633 42.8167 117.293 116.382 113 46.4058 40.2667 37.0083 34.1058 32.0492 30.2275 52.811 37.6517 29.7158 51.131 40.3125 28.6492 60.513 44.4508 36.9592 106.935 80.137 61.382 128.85 90.322 63.292 90.32 61.108 38.8683 77.451 53.187 37.98 19.385 17.75 16.5475 NA NA NA 12.0192 9.31 6.4117 80.368 53.694 30.8075 40.47 31.3642 25.64 232.623 186.612 167.66 186.64 168.857 156.915 36.1883 32.5633 29.8467 46.6892 38.2817 36.18 32.1808 31.36 30.2175 52.88 50.992 49.0442 23.6625 21.5617 19.4183 54.893 39.3067 22.6775 24.4592 16.6033 11.63 78.593 54.771 37.6317 92.13 76.097 49.6192 99.357 71.487 49.3683 49.8117 36.1025 26.045 9.1642 9.1483 8.9908 NA NA NA 21.68 18.7567 15.86 22.6333 18.7658 15.0067 57.589 40.3492 29.7 Trang66 Phụ lục 6: Dữ liệu PM 2.5 đo trạm Phú Thọ tƣơng ứng với liệu MODIS Aqua Terra đƣa vào thực nghiệm Time 1/1/2014 6:10 1/17/2014 6:10 2/2/2014 6:10 2/18/2014 6:10 3/6/2014 6:10 3/22/2014 6:10 4/7/2014 6:10 4/23/2014 6:10 5/9/2014 6:10 5/25/2014 6:10 6/10/2014 6:10 6/26/2014 6:10 7/12/2014 6:10 7/28/2014 6:10 8/13/2014 6:10 8/29/2014 6:10 9/14/2014 6:10 9/30/2014 6:10 1/2/2014 3:45 1/18/2014 3:45 2/3/2014 3:45 2/19/2014 3:45 3/7/2014 3:45 3/23/2014 3:45 4/8/2014 3:45 4/24/2014 3:45 5/10/2014 3:45 5/26/2014 3:45 6/11/2014 3:45 6/27/2014 3:45 7/13/2014 3:45 7/29/2014 3:45 8/14/2014 3:45 8/30/2014 3:45 9/15/2014 3:45 10/1/2014 3:45 PM.10 115.874 139.778 27.9533 39.3642 43.9575 27.6667 11.0733 23.1267 3.96833 NA NA NA NA NA 3.6525 7.9458 73.025 35.1583 194.617 127.133 22.025 20.705 23.3075 33.3217 12.7892 18.6675 0.79167 NA NA NA NA NA 33.9517 15.4533 43.0425 67.562 PM.2.5 92.876 105.557 26.8992 38.3075 42.9025 25.5208 10.4433 20.52 3.92167 NA NA NA NA NA 3.17917 6.3317 49.1433 23.4683 148.168 108.59 21.2533 20.26 22.8108 31.085 12.6292 17.5825 0.79167 NA NA NA NA NA 22.8983 11.7783 37.5775 50.075 PM.1 86.322 95.424 25.3383 37.6825 41.195 24.1867 9.6675 17.4075 3.9025 NA NA NA NA NA 2.58583 4.72417 36.7925 15.9317 133.727 101.318 19.9983 19.88 22.3767 29.9775 12.1808 15.9917 0.78333 NA NA NA NA NA 16.765 9.1483 31.3192 41.7983 Trang67 ... PHƢƠNG PHÁP ĐÁNH GIÁ CHẤT LƢỢNG KHÔNG KHÍ KHU VỰCĐỒNG BẰNG SÔNG HỒNG BẰNG TƢ LIỆU VỆ TINH MODIS 2.1 Cơ sở khoa học phƣơng pháp đánh giá chất lƣợng không khí khu vực ĐBSH tƣ liệu vệ tinh 2.1.1... sĩ: Sử dụng tư liệu vệ tinh MODIS đánh giá chất lượng không khí khu vực Đồng sông Hồng" với mong muốn góp phần giải phần vấn đề nhƣ trình bày, góp sức mục tiêu chung bảo vệ sức khỏe cộng đồng. .. khỏe cộng đồng 13 1.3 Tổng quan phát ô nhiễm không khí độ dày quang học sol khí 15 CHƢƠNG PHƢƠNG PHÁP ĐÁNH GIÁ CHẤT LƢỢNG KHÔNG KHÍ 18 KHU VỰC ĐỒNG BẰNG SÔNG HỒNG BẰNG TƢ LIỆU VỆ TINH MODIS