From 29db69ea2cc516eeaed8e5f182d0e95309b23fed Mon Sep 17 00:00:00 2001 From: zhaojinghao Date: Thu, 10 Nov 2022 16:17:29 +0800 Subject: [PATCH] change to par --- AOD_NetCDF_to_GeoTIFF.py | 2 +- main.py | 10 ++++++++++ 2 files changed, 11 insertions(+), 1 deletion(-) diff --git a/AOD_NetCDF_to_GeoTIFF.py b/AOD_NetCDF_to_GeoTIFF.py index 4a2a6a6..e028b16 100644 --- a/AOD_NetCDF_to_GeoTIFF.py +++ b/AOD_NetCDF_to_GeoTIFF.py @@ -14,7 +14,7 @@ def NC_to_tiffs(data, Output_folder): # 读取变量的时候，会自动根据scale factor对数值进行还原，但是Nodata的栅格会存储为-32768 # 无论是日数据还是小时数居，变量名都是"SWR" - AOD_arr = np.asarray(nc_data_obj.variables["SWR"]) # 将SWR数据读取为数组 + AOD_arr = np.asarray(nc_data_obj.variables["PAR"]) # 将SWR数据读取为数组 # 这个循环将所有Nodata的值（即-32768）全部改为0 for i in range(len(AOD_arr)): diff --git a/main.py b/main.py index 064c4cc..26c837a 100644 --- a/main.py +++ b/main.py @@ -1,6 +1,8 @@ import os import Download_AOD as Daod import datetime as dt +import glob +import AOD_NetCDF_to_GeoTIFF as trans today = dt.date(2022, 7, 22) print(today) @@ -96,6 +98,14 @@ if __name__ == "__main__": # Analysis_Path = dst_filePath + "/PAR_Minutes_Analysis" # if not os.path.exists(Analysis_Path): # os.makedirs(Analysis_Path) + data_list = glob.glob(Download_Path + "\\*.nc") + + # for循环完成解析 + for i in range(len(data_list)): + data = data_list[i] + trans.NC_to_tiffs(data, Analysis_Path) + print(data + "-----转tif成功") + today = today + dt.timedelta(days=1) else: print("选择错误！")