Luận án Tiến sĩ Nghiên cứu xây dựng mô hình giám sát sự bốc thoát hơi nước của lớp phủ khu vực Tây Bắc Việt Nam từ dữ liệu ảnh vệ tinh

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC MỎ - ĐỊA CHẤT LÊ HÙNG CHIẾN NGHIÊN CỨU XÂY DỰNG MƠ HÌNH GIÁM SÁT SỰ BỐC - THỐT HƠI NƯỚC CỦA LỚP PHỦ KHU VỰC TÂY BẮC VIỆT NAM TỪ DỮ LIỆU ẢNH VỆ TINH LUẬN ÁN TIẾN SĨ NGÀNH: KỸ THUẬT TRẮC ĐỊA - BẢN ĐỒ HÀ NỘI, 2022 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC MỎ - ĐỊA CHẤT LÊ HÙNG CHIẾN NGHIÊN CỨU XÂY DỰNG MƠ HÌNH GIÁM SÁT SỰ BỐC - THỐT HƠI NƯỚC CỦA LỚP PHỦ KHU VỰC TÂY BẮC VIỆT NAM TỪ DỮ LIỆU ẢNH VỆ TINH Ngành: Kỹ thuật Trắc địa - Bản đồ Mã số: 9.520503 Người hướng dẫn khoa học: PGS.TS TRẦN XUÂN TRƯỜNG PGS.TS DOÃN HÀ PHONG HÀ NỘI, 2022 i LỜI CAM ĐOAN Tôi xin cam đoan cơng trình nghiên cứu riêng tơi Q trình nghiên cứu thực nghiêm túc, khoa học Số liệu kết trình bày luận án xác, trung thực chưa cơng bố cơng trình khác Tác giả Lê Hùng Chiến ii MỤC LỤC LỜI CAM ĐOAN i MỤC LỤC ii DANH MỤC BẢNG BIỂU ix DANH MỤC HÌNH xi MỞ ĐẦU 1 Tính cấp thiết đề tài Mục tiêu nghiên cứu 3 Đối tượng nghiên cứu 4 Phạm vi nghiên cứu Nội dung nghiên cứu luận án Phương pháp nghiên cứu Ý nghĩa khoa học thực tiễn luận án Luận điểm bảo vệ Những điểm luận án 10 Cấu trúc luận án 11 Lời cảm ơn CHƯƠNG TỔNG QUAN VẤN ĐỀ NGHIÊN CỨU .9 1.1 Các khái niệm bốc thoát nước 1.1.1 Bốc nước (E) 1.1.2 Thoát nước (T) 10 1.1.3 Bốc thoát nước ET (Evaporation Transpiration) 11 1.1.4 Các yếu tố ảnh hưởng đến bốc thoát nước 11 1.1.5 Bốc thoát nước tham chiếu ET0 (Potential evaptransporation) 13 1.1.6 Thoát nước điều kiện tiêu chuẩn (ETc) 14 1.1.7 Thoát nước điều kiện không tiêu chuẩn (ETc adj) 14 1.1.8 Lượng bốc thoát thực tế ET (Actual evapotransporation) 15 1.1.9 Mơ hình ước tính giám sát lượng bốc thoát nước 15 1.2 Phương pháp xác lượng bốc thoát nước sử dụng liệu khí tượng 15 1.2.1 Các phương pháp đo trực tiếp 15 1.2.2 Các mô hình sử dụng lượng xạ mặt trời (radiaton-based models) 18 1.2.3 Các mơ hình kết hợp (combined models) 20 1.3 Các mơ hình xác ước tính lượng bốc nước từ liệu ảnh vệ tinh 22 iii 1.3.1 Mô hình cân lượng bề mặt đất SEBAL (Surface Energy Balance Algorithms for Land) 23 1.3.2 Mô hình số cân lượng bề mặt SEBI (Surface Energy Balance Index) 25 1.3.3 Mơ hình Hệ thống cân lượng bề mặt SEBS (Surface Energy Balance System) 27 1.3.4 Mô hình số cân lượng xạ bề mặt đơn giản S-SEBI (Simplified Surface Energy Balance Index) 28 1.3.5 Mô hình đồ bốc nước độ phân giải cao với hiệu chỉnh bên METRIC (Mapping ET with Internalized Calibration) 30 1.4 Các kết nghiên cứu giới liên quan đến đề tài 33 1.5 Các kết nghiên cứu nước liên quan đến lĩnh vực đề tài 36 1.6 Đánh giá chung phương pháp mơ hình ước tính lượng bốc nước từ bề mặt lớp phủ 39 1.7 Một số vấn đề thảo luận phát triển luận án 42 Tiểu kết chương 44 CHƯƠNG CƠ SỞ KHOA HỌC SỬ DỤNG DỮ LIỆU ẢNH VỆ TINH PHỤC VỤ ƯỚC TÍNH, GIÁM SÁT LƯỢNG BỐC THOÁT HƠI NƯỚC BỀ MẶT LỚP PHỦ 46 2.1 Khái quát viễn thám 46 2.1.1 Nguyên lý viễn thám 46 2.1.2 Đặc tính phản xạ thực vật 49 2.1.3 Đặc tính phản xạ phổ nước 50 2.1.4 Đặc tính phản xạ đối tượng đô thị 51 2.2 Đặc điểm ảnh vệ tinh Landsat 51 2.3 Vai trò liệu ảnh vệ tinh Landsat việc chiết xuất, tính tốn tham số phục vụ ước tính lượng bốc nước 53 2.4 Khả ứng dụng liệu ảnh vệ tinh Landsat việc ước tính lượng bốc thoát nước bề mặt lớp phủ 55 2.4.1 Tính giá trị lượng xạ rịng mặt trời (Rn) từ liệu ảnh vệ tinh Landsat 56 2.4.2 Xác định giá trị nhiệt ẩn q trình bốc nước λ từ liệu ảnh vệ tinh Landsat 58 2.4.3 Xác định số Psychrometric (γ) từ liệu ảnh vệ tinh Landsat mơ hình số độ cao DEM 58 iv 2.4.4 Sử dụng ảnh vệ tinh Landsat xác định giá trị độ dốc đường cong áp suất bão hòa (Δ) 59 2.5 Tính giá trị tham số từ ảnh vệ tinh Landsat thông tin độ cao phục vụ ước tính, giám sát lượng bốc nước bề mặt lớp phủ 59 2.5.1 Tính giá trị lượng xạ ròng Rni từ ảnh vệ tinh Landsat theo mơ hình SEBAL 59 2.5.2 Tính giá trị xạ rịng trung bình ngày Rnd từ Rni tính từ ảnh vệ tinh Landsat 66 2.5.3 Tính xạ rịng trung bình ngày Rnd từ số liệu khí tượng đo trực tiếp trạm quan trắc theo mơ hình FAO 56 – Penman - Monteith 67 2.5.4 Tính giá trị nhiệt ẩn q trình bốc nước (λ) từ liệu ảnh vệ tinh Landsat 68 2.5.5 Tính giá trị số Psychrometric (γ) từ liệu ảnh vệ tinh Landsat thông tin độ cao từ DEM 69 2.5.6 Tính giá trị độ dốc đường cong áp suất bão hòa (Δ) từ liệu ảnh vệ tinh Landsat 69 2.6 Xác định hệ số a, b mơ hình Priestley - Taylor phù hợp với địa hình khí hậu khu vực Tây Bắc Việt Nam 70 2.7 Đề xuất mơ hình, quy trình ước tính, giám sát lượng bốc nước sử dụng kết hợp mơ hình SEBAL với liệu ảnh vệ tinh Landsat mô hình Priestley Taylor 71 2.7.1 Đề xuất mơ hình ước tính giám sát lượng bốc thoát nước bề mặt lớp phủ khu vực Tây Bắc Việt Nam 71 2.7.2 Quy trình ước tính, giám sát lượng bốc thoát nước sử dụng kết hợp mơ hình SEBAL với liệu ảnh vệ tinh Landsat mơ hình Priestley - Taylor 72 2.8 Xây dựng chương trình ước tính, giám sát lượng bốc nước sử dụng kết hợp mơ hình SEBAL mơ hình Priestley-Taylor Google Earth Engine 75 2.8.1 Khái quát Google Earth Engine 75 2.8.2 Những lợi ích Google Earth Engine việc xây dựng chương trình ước tính, giám sát lượng bốc nước 75 Tiểu kết chương 77 CHƯƠNG KẾT QUẢ THỰC NGHIỆM VÀ THẢO LUẬN .78 3.1 Điều kiện tự nhiên tỉnh Hịa Bình 78 v 3.1.1 Vị trí địa lý 78 3.1.2 Địa hình, địa mạo 79 3.1.3 Điều kiện khí hậu 79 3.1.4 Thủy văn 82 3.1.5 Thực trạng cấu tài nguyên đất tỉnh Hịa Bình 83 3.1.6 Tài nguyên nước 84 3.1.7 Tài nguyên rừng 84 3.2 Dữ liệu phục vụ nghiên cứu 85 3.2.1 Dữ liệu ảnh vệ tinh Landsat mơ hình số độ cao DEM (SRTM) 85 3.2.2 Dữ liệu khí tượng 87 3.3 Thực nghiệm xác định hệ số a, b mơ hình Priestley – Taylor với điều kiện địa hình, khí hậu khu vực Tây Bắc Việt Nam từ liệu quan trắc khí tượng, thủy văn tỉnh Hịa Bình 90 3.3.1 Kết tính giá trị lượng xạ rịng Rnd từ liệu khí tượng, thủy văn đo trực tiếp trạm khí tượng thủy văn Hịa Bình theo mơ hình FAO 56 91 3.3.2 Kết tính giá trị nhiệt ẩn q trình bốc thoát nước (λ), số Psychrometric (γ), độ dốc đường cong áp suất nước bão hòa (Δ) từ liệu khí tượng thủy văn đo trực tiếp trạm khí tượng thủy văn Hịa Bình theo mơ hình FAO 56 93 3.3.3 Kết tính hệ số a, b mơ hình Priestley – Taylor với điều kiện địa hình, khí hậu tỉnh Hịa Bình, Sơn La thuộc vùng Tây Bắc Việt Nam từ liệu quan trắc khí tượng thủy văn 95 3.4 Thực nghiệm tính lượng bốc nước thực tế từ bề mặt lớp phủ tỉnh Hịa Bình khu vực Tây Bắc Việt Nam sử dụng kết hợp mơ hình SEBAL với liệu ảnh vệ tinh Landsat mơ hình Priestley – Taylor với hệ số a, b xác định thực nghiệm 97 3.4.1 Kết tính giá trị lượng xạ rịng trung bình ngày Rnd từ liệu ảnh vệ tinh Landsat (Rnd _VT) 97 3.4.2 So sánh kết tính giá trị xạ rịng trung bình ngày theo mơ hình FAO 56 (Rnd_FAO) lượng xạ ròng trung bình ngày tính từ ảnh vệ tinh Landsat (Rnd_VT) 109 3.4.3 Kết tính giá trị nhiệt ẩn q trình bốc nước (λ) với tham số nhiệt độ bề mặt tính từ ảnh Landsat 114 3.4.4 Kết tính giá trị số Psychrometric (γ) với giá trị độ cao chiết xuất từ DEM liệu ảnh vệ tinh Landsat 116 vi 3.4.5 Kết tính giá trị độ dốc đường cong áp suất bão hòa (Δ) với tham số nhiệt độ bề mặt tính từ ảnh vệ tinh Landsat 117 3.4.6 Tính lượng bốc nước thực tế ETa theo mơ hình Priestley – Taylor với tham số chiết xuất, tính tốn từ liệu ảnh vệ tinh Landsat giá trị độ cao từ DEM tỉnh Hịa Bình khu vực Tây Bắc Việt Nam 119 3.4.7 So sánh lượng bốc nước thực tế đo trạm khí tượng thủy văn ETa_Đo lượng bốc thoát nước tính sử dụng kết hợp mơ hình viễn thám SEBAL với liệu ảnh vệ tinh Landsat mô hình Priestley – Taylor ETa_VT 123 3.5 Xây dựng chương trình ước tính, giám sát lượng bốc nước bề mặt lớp phủ từ liệu ảnh vệ tinh Google Earth Engine 126 3.5.1 Sơ đồ khối chương trình ước tính giám sát lượng bốc thoát nước bề mặt lớp phủ 126 3.5.2 Các giao diện chương trình 127 Tiểu kết Chương 129 KẾT LUẬN KIẾN NGHỊ .131 Kết luận 131 Kiến nghị 132 DANH MỤC CÁC CƠNG TRÌNH KHOA HỌC, BÀI BÁO ĐÃ CÔNG BỐ 133 TÀI LIỆU THAM KHẢO .134 vii DANH MỤC KÍ HIỆU, CHỮ VIẾT TẮT Kí hiệu Tiếng Anh Tiếng Việt E Evaporation Bốc T Transpiration Thoát ETa Actual Evapotransporation Bốc thoát nước thực tế ET0 Reference Evapotranspiration Bốc thoát nước tham chiếu Rs Solar radiation Năng lượng xạ mặt trời Rn Net radiation Năng lượng xạ ròng Rns Short radiation Rnl Long radiation N Daylight hours Năng lượng xạ mặt trời sóng ngắn Năng lượng xạ mặt trời sóng dài Tổng nắng ngày Ts Surface temperature Nhiệt độ bề mặt Ta Near surface temperature NDVI SAVI Normalised Difference Vegetation Index Soil Adjusted Vegetation Index Nhiệt độ khơng khí gần bề mặt đất Chỉ số thực vật LAI Leaf Area Index Chỉ số thực vật có hiệu chỉnh phản xạ mặt đất Chỉ số diện tích α Surface albedo Suất phân sai bề mặt ε0 Broad band surface emissivity Phát xạ bề mặt dải rộng εNB Narrow band surface emissivity Phát xạ bề mặt dải hẹp λ Latent heat of vaporization γ Psychrometric constant Năng lượng nhiệt hóa tiềm ẩn Hằng số Psychrometric Δ Slope curve saturation vapour pressure Atmospheric pressure Độ dốc đường cong áp suất bão hòa Áp suất khơng khí Adjustment of wind speed to standard height Altitude above sea level Vận tốc gió độ cao tiêu chuẩn Food and Agriculture Organization of the United Nations Tổ chức Nông nghiệp Lương thực Liên hợp quốc P u2 z FAO Độ cao tuyệt đối điểm viii Landsat Landsatellite Ảnh vệ tinh NASA, Mỹ SEBAL SEBI Surface Energy Balance Algorithms for Land Surface Energy Balance Index SEBS Surface Energy Balance System METRIC Mapping ET with Internalized Calibration Simplified - Surface Energy Balance Index Google earth engine Mô hình cân lượng bề mặt đất Mơ hình số cân lượng bề mặt Mơ hình hệ thống cân lượng bề mặt Mơ hình đồ bốc thoát nước với hiệu chỉnh nội Mơ hình số cân lượng bề mặt đơn giản S- SEBI GGE RH độ e0 (T) ea (Kpa) Tmax (0C) Tmin (0C) Rnl (MJ/m2/d) Rns (MJ/m2/d) Rn FAO (MJ/m2/d) 66 5.304470 3.500950 37.7 30.2 1.291 12.341 11.049 76 4.454369 3.385321 35.5 26.2 1.160 11.643 10.483 82 4.781710 3.921002 36.4 27.8 0.832 11.161 10.328 67 4.932233 3.304596 35.9 29.4 1.143 11.176 10.033 76 4.795227 3.644373 36.4 27.9 0.774 10.137 9.363 67 5.393751 3.613813 38.5 30.0 1.205 12.225 11.020 68 4.701301 3.196885 36.2 27.4 1.341 11.869 10.528 72 5.231050 3.766356 37.4 30.0 1.146 12.458 11.312 72 3.846461 2.769452 30.2 26.4 1.154 10.630 9.476 49 5.158514 2.527672 36.9 30.0 1.774 11.900 10.127 70 3.351517 2.346062 28.4 23.5 0.830 8.997 8.167 Thời điểm ngày 07/7/2020 TT Độ vĩ φ (độ) tg (φ) ω_s (δ) sin(δ) cos(δ) sin(ω_s) cos(ω_s) sin(φ) 20.8167 0.379987 99.133475 0.393761 0.383664 0.923473 0.987460 -0.157869 0.355207 20.6500 0.376663 99.053309 0.393761 0.383664 0.923473 0.987680 -0.156488 0.352487 20.0667 0.365088 98.774305 0.393761 0.383664 0.923473 0.988430 -0.151679 0.342947 20.4833 0.373346 98.973336 0.393761 0.383664 0.923473 0.987897 -0.155110 0.349765 20.4500 0.372684 98.957383 0.393761 0.383664 0.923473 0.987940 -0.154835 0.349220 20.8167 0.379987 99.133475 0.393761 0.383664 0.923473 0.987460 -0.157869 0.355207 20.5167 0.374010 98.989344 0.393761 0.383664 0.923473 0.987854 -0.155386 0.350310 20.8833 0.381317 99.165557 0.393761 0.383664 0.923473 0.987372 -0.158421 0.356293 21.3833 0.391341 99.407427 0.393761 0.383664 0.923473 0.986694 -0.162586 0.364429 10 21.4333 0.392347 99.431713 0.393761 0.383664 0.923473 0.986625 -0.163004 0.365241 11 20.8301 0.380255 99.139927 0.393761 0.383664 0.923473 0.987442 -0.157980 0.355426 Rs (MJ/m2 /ngày) Rs/Rso cos(φ) dr Ra (MJ/m2/ ngày) Z (m) Rso (MJ/m2 /ngày) n (Giờ nắng thực) N (Giờ nắng vĩ độ) 0.934788 0.967143 39.574 22.7 29.698 9.8 13.22 24.564 0.827116 0.935817 0.967143 39.542 165.5 29.787 7.6 13.21 21.263 0.713815 0.939355 0.967143 39.430 61.1 29.621 9.3 13.17 23.779 0.802795 0.936838 0.967143 39.510 11.3 29.642 8.6 13.20 22.752 0.767562 0.937041 0.967143 39.504 41.2 29.661 9.1 13.19 23.499 0.792257 0.934788 0.967143 39.574 22.6 29.698 9.6 13.22 24.264 0.817036 0.936634 0.967143 39.517 20.1 29.653 9.4 13.20 23.951 0.807700 0.934374 0.967143 39.586 25.4 29.710 9.2 13.22 23.669 0.796665 0.931231 0.967143 39.679 642.9 30.269 10.1 13.25 25.037 0.827163 0.930913 0.967143 39.688 169.0 29.900 9.9 13.26 24.740 0.827433 0.934705 0.967143 39.576 971.9 30.451 9.0 13.22 23.367 0.767350 RH độ e0 (T) ea (Kpa) Tmax (0C) Tmin (0C) Rnl (MJ/m2/d) Rns (MJ/m2/d) Rn FAO (MJ/m2/d) 79 4.863311 3.842016 37.7 27.1 2.153 18.914 16.761 82 4.835978 3.965502 38.0 26.6 1.607 16.372 14.766 82 4.622069 3.790097 36.7 26.3 2.094 18.310 16.216 80 4.768226 3.814581 36.4 27.7 1.946 17.519 15.573 79 4.768226 3.766899 36.9 27.2 2.094 18.094 16.001 81 4.835978 3.917142 37.6 27.0 2.026 18.684 16.658 79 4.714620 3.724550 36.9 26.8 2.197 18.442 16.245 80 4.768226 3.814581 36.8 27.3 2.058 18.225 16.167 83 4.040532 3.353641 33.3 25.0 2.628 19.279 16.651 77 4.206695 3.239155 37.7 22.0 2.802 19.050 16.248 83 3.994118 3.315118 31.5 26.4 2.385 17.992 15.607 Thời điểm ngày 18/8/2021 TT Độ vĩ φ (độ) tg (φ) ω_s (δ) sin(δ) cos(δ) sin(ω_s) cos(ω_s) sin(φ) 20.8167 0.379987 94.973076 0.222199 0.220375 0.975415 0.996308 -0.085850 0.355207 20.6500 0.376663 94.929864 0.222199 0.220375 0.975415 0.996372 -0.085099 0.352487 20.0667 0.365088 94.779423 0.222199 0.220375 0.975415 0.996592 -0.082484 0.342947 20.4833 0.373346 94.886750 0.222199 0.220375 0.975415 0.996436 -0.084350 0.349765 20.4500 0.372684 94.878149 0.222199 0.220375 0.975415 0.996449 -0.084200 0.349220 20.8167 0.379987 94.973076 0.222199 0.220375 0.975415 0.996308 -0.085850 0.355207 20.5167 0.374010 94.895381 0.222199 0.220375 0.975415 0.996423 -0.084500 0.350310 20.8833 0.381317 94.990367 0.222199 0.220375 0.975415 0.996282 -0.086151 0.356293 21.3833 0.391341 95.120690 0.222199 0.220375 0.975415 0.996084 -0.088416 0.364429 10 21.4333 0.392347 95.133772 0.222199 0.220375 0.975415 0.996063 -0.088643 0.365241 11 20.8301 0.380255 94.976553 0.222199 0.220375 0.975415 0.996303 -0.085911 0.355426 Rs (MJ/m2 /ngày) Rs/Rso cos(φ) dr Ra (MJ/m2/ ngày) Z (m) Rso (MJ/m2 /ngày) n (Giờ nắng thực) N (Giờ nắng vĩ độ) 0.934788 0.977382 38.156 22.7 28.634 10.2 12.66 24.906 0.869801 0.935817 0.977382 38.156 165.5 28.743 7.5 12.66 20.844 0.725162 0.939355 0.977382 38.155 61.1 28.663 9.9 12.64 24.484 0.854207 0.936838 0.977382 38.156 11.3 28.626 9.8 12.65 24.317 0.849482 0.937041 0.977382 38.156 41.2 28.649 9.3 12.65 23.564 0.822533 0.934788 0.977382 38.156 22.6 28.634 10.0 12.66 24.605 0.859279 0.936634 0.977382 38.156 20.1 28.632 9.6 12.65 24.014 0.838705 0.934374 0.977382 38.156 25.4 28.636 9.9 12.67 24.451 0.853861 0.931231 0.977382 38.153 642.9 29.105 10.3 12.68 25.031 0.860007 0.930913 0.977382 38.152 169.0 28.743 11.4 12.68 26.682 0.928306 0.934705 0.977382 38.156 971.9 29.359 9.6 12.66 24.002 0.817533 RH độ e0 (T) ea (Kpa) Tmax (0C) Tmin (0C) Rnl (MJ/m2/d) Rns (MJ/m2/d) Rn FAO (MJ/m2/d) 81 4.467079 3.618334 37.2 24.6 2.552 19.178 16.626 83 4.467079 3.707675 36.9 24.9 1.861 16.050 14.189 84 4.328993 3.636354 35.8 24.9 2.445 18.853 16.407 81 4.122885 3.339537 36.0 23.0 2.766 18.724 15.958 84 4.087414 3.433428 36.2 22.5 2.524 18.145 15.621 82 4.429044 3.631816 37.0 24.5 2.486 18.946 16.459 84 4.194633 3.523491 35.9 23.7 2.501 18.491 15.990 83 4.316625 3.582799 36.4 24.2 2.509 18.828 16.318 74 4.052208 2.998634 32.1 26.3 3.244 19.274 16.029 74 4.674760 3.459323 35.5 27.9 3.048 20.546 17.497 80 3.381362 2.705089 28.7 23.5 3.254 18.481 15.227 Phụ lục 03: Tính hệ số a, b mơ hình Priestley - Taylor từ liệu khí tượng đo trực tiếp trạm khí tượng thủy văn Sử dụng giá trị lượng xạ rịng, Nhiệt ẩn q trình bốc nước, Hằng số Psychrometric, Độ dốc đường cong áp suất nước bão hịa tính từ liệu khí tượng thủy văn lượng bốc thoát nước thực tế đo trực tiếp trạm khí tượng thủy văn để xác định hệ số a, b mô hình Priestley -Taylor Thời điểm ngày 01/7/2015 Ma trận phương trình số hiệu chỉnh Ai Bi Li 9.194 1.0 -7.1 7.073 1.0 -5.7 7.028 1.0 -5.5 8.910 1.0 -7.2 6.559 1.0 -5.0 9.141 1.0 -7.2 8.839 1.0 -6.8 10.028 1.0 -7.9 Hệ số a, b tính theo số liệu ngày 01/7/2015 a1 0.795 b1 -0.087 Thời điểm ngày 15/6/2016 Ma trận phương trình số hiệu chỉnh Ai Bi Li 9.025 1.0 -6.7 8.505 1.0 -6.4 8.270 1.0 -6.3 8.377 1.0 -6.2 8.588 1.0 -6.6 8.925 1.0 -7.0 8.939 1.0 -6.8 9.201 1.0 -6.9 Hệ số a, b tính theo số liệu ngày 15/6/2016 a2 0.771 b2 -0.118 Thời điểm ngày 04/6/2017 Ma trận phương trình số hiệu chỉnh Ai Bi Li 8.984 1.0 -9.6 5.995 1.0 -5.7 8.878 1.0 -7 9.676 1.0 -7.8 5.995 1.0 -4.6 9.038 1.0 -9.5 9.155 1.0 -7.4 9.449 1.0 -7.3 Hệ số a, b tính theo số liệu ngày 04/6/2017 a3 0.862 b3 0.123 Thời điểm ngày 25/7/2018 Ma trận phương trình số hiệu chỉnh Ai Bi Li 5.594 1.000 -4.400 4.498 1.000 -3.500 6.340 1.000 -4.900 6.247 1.000 -4.800 5.803 1.000 -4.637 5.731 1.000 -3.200 6.848 1.000 -5.200 5.778 1.000 -4.617 Hệ số a, b tính theo số liệu ngày 25/7/2018 a4 0.767 b4 -0.081 Thời điểm ngày 09/7/2019 Ma trận phương trình số hiệu chỉnh Ai Bi Li 5.803 1.0 -4.4 5.737 1.0 -4.5 5.564 1.0 -4.9 5.370 1.0 -4.2 5.043 1.0 -4.0 5.765 1.0 -4.7 5.328 1.0 -4.4 4.726 1.0 -3.7 Hệ số a, b tính theo số liệu ngày 09/7/2019 Thời điểm ngày 07/7/2020 a5 b5 0.835 -0.171 Ma trận phương trình số hiệu chỉnh Ai Bi Li 9.001 1.0 -7.2 7.899 1.0 -6.2 8.818 1.0 -7.0 8.411 1.0 -6.7 8.632 1.0 -7.2 8.959 1.0 -7.0 9.046 1.0 -7.1 8.530 1.0 -6.8 Hệ số a, b tính theo số liệu ngày 07/7/2020 a6 0.782 b6 0.123 Thời điểm ngày 18/8/2021 Ma trận phương trình số hiệu chỉnh Ai Bi Li 9.145 1.0 -7.0 7.759 1.0 -5.9 8.615 1.0 -6.2 9.005 1.0 -6.7 8.828 1.0 -5.8 9.077 1.0 -6.8 9.445 1.0 -7.0 8.784 1.0 -6.4 Hệ số a, b tính theo số liệu ngày 18/8/2021 a7 0.730 b7 0.029 Phụ lục 04: Code chương trình ước tính, giám sát lượng bốc nước bề mặt lớp phủ từ liệu ảnh vệ tinh Landsat var poi_hoabinh = ee.Geometry.Point([105.33362975494495, 20.6666483915459]), aoi_hoabinh = ee.Geometry.Polygon( [[[104.83650913612476, 21.316023838232816], [104.83650913612476, 20.438374359233357], [105.80330601112476, 20.438374359233357], [105.80330601112476, 21.316023838232816]]], null, false), reg_vietnam = ee.FeatureCollection("users/hoatx/rg_vn_province"); var dem = ee.Image("USGS/SRTMGL1_003"); var raw1 = ee.Image('VNUF_GEE/LC08_C01_T1/LC08_127046_20150701'); var raw2 = ee.Image('VNUF_GEE/LC08_C01_T1/LC08_127046_20170604'); var raw3 = ee.Image('VNUF_GEE/LC08_C01_T1/LC08_127046_20210818'); var reg_hoabinh = reg_vietnam.filterMetadata('MATINH','equals', 17); // CÁC HÀM XỬ LÝ TRUNG GIAN -var cloudMask = function(image) { var qa = image.select('BQA'); /// Check that the cloud bit is off // See https://www.usgs.gov/land-resources/nli/landsat/landsat-collection-1-level-1quality-assessment-band var mask = qa.bitwiseAnd(1