STT Tên trạm Vĩ độ Kinh độ Nhiệt độ không khí (oC) 11/12/2019 11h 11/12/2019 14 h 1 Bái Thượng 19.900 105.383 16.1 24.8 2 Hồi Xuân 20.367 105.117 20.0 24.0 3 Như Xuân 19.633 105.567 20.8 23.5 4 Sầm Sơn 19.750 105.900 19.5 22.3 5 Thanh Hoá 19.750 105.783 20.4 22.5 6 Tĩnh Gia 19.450 105.783 20.0 22.8 7 Yên Định 19.983 105.667 21.1 23.3
STT Tên trạm Vĩ độ Kinh độ Nhiệt độ không khí (
oC) 29/09/2019 11h 29/09/2019 14 h 1 Bái Thượng 19.900 105.383 31.5 34.4 2 Hồi Xuân 20.367 105.117 27.6 33.2 3 Như Xuân 19.633 105.567 30.5 32.8 4 Sầm Sơn 19.750 105.900 31.4 32.6 5 Thanh Hoá 19.750 105.783 30.3 32.9 6 Tĩnh Gia 19.450 105.783 31.2 33.1 7 Yên Định 19.983 105.667 30.9 33.2
Với hai dữ liệu này, nhiệt độ được thử nghiệm với thời điểm nhiệt độ thấp và nhiệt độ cao để đánh giá phương pháp.
4.2.2. Thực nghiệm tính toán nhiệt độ không khí.
1. Xử lý dữ liệu và đăng ký hệ toạ độ.
Dữ liệu DEM được chuyển đổi cắt ghép theo khu vực nghiên cứu như hình dưới đây:
Hình 4.4. Mô hình số địa hình tại khu vực đã được cắt ghép và xử lý
Dữ liệu NDVI được chuyển đổi cắt ghép theo khu vực nghiên cứu.
Hình 4.6. Chỉ số NDVI đã xử lý ngày 11/12/2019
Dữ liệu MODIS ở dạng HDF được tải về sẽ được chuyển đổi về hệ toạ độ địa lý. Với mỗi thời điểm, có 2 ảnh nên hai ảnh được ghép lại với nhau để có nhiệt độ bề mặt từng thời điểm trong khu vực nghiên cứu. Đồng thời chuyển đổi nhiệt độ K về nhiệt độ C. Công việc này được thực hiện trên phần mềm ENVI.
2. Tính toán các ảnh của các biến đầu vào.
Ảnh giá trị nghịch đảo DEM (DEM_ND) và tỉ lệ chỉ số thực vật (NDVIper)
Thực hiện nghịch đảo DEM và tỉ lệ chỉ số thực vật bằng phần mềm ENVI:
DEM_ND = 1/DEM
NDVIper = (2-(NDVI+1))/2
3. Trích xuất các giá trị lập hàm hồi quy Sử dụng phần mềm GIS trích xuất các giá trị: - Nhiệt độ không khí (từ dữ liệu trạm);
- Nhiệt độ bề mặt đất (từ ảnh nhiệt độ bề mặt tính từ MODIS); - Tham số nghịch đảo độ cao (từ ảnh nghịch đảo độ cao) - Tham số tỉ lệ chỉ số thực vật (từ ảnh tỉ lệ chỉ số thực vật)
4. Tiến hành hồi quy đa biến
Hồi quy đa biến (3 biến) tính toán tham số và tương quan. Trong quá trình tính toán phát hiện dữ liệu ở các trạm không tương quan với các biến thì tiến hành điều tra và loại bỏ nếu cần thiết. Hồi quy theo công thức:
T_air = a+ b*T_suf + c* NDVIper + d* DEM_ND T_air :nhiệt độ không khí bằng nhiệt độ đo tại trạm T_suf: nhiệt độ bề mặt từ dữ liệu MODIS.
NDVIper và DEM_ND được chiết xuất ở trên
Với thời điểm 11/12/2019 lúc 11h số liệu hồiquy như sau: