def extract_Sentinel_temp_list(sen_directory, sen_shape_path, daytime_S3=None):
"""
returns the temperature value for each pixel of each station for each SENTINEL scene
prints an array for the scenes according to the stations
:return:
returns three arrays (1st array displayed for each station is the temperature array)
returns one temperature array for each station if a[0][0]
"""
import_list = import_polygons(shape_path=sen_shape_path)
sentinel_file_list = extract_files_to_list(path_to_folder=sen_directory, datatype=".tif")
Sen_station_time_series = []
for i, polygons in enumerate(import_list):
## Initialize empty analysis lists
Sen_final_mean = []
for j, tifs in enumerate(sentinel_file_list):
# if daytime_S3 in str(sentinel_file_list[j]):
src1 = rio.open(sentinel_file_list[j])
mask = rio.mask.mask(src1, [import_list[0][i]], all_touched=True, crop=True, nodata=np.nan)
Sen_temperature_array = mask[0][0]
## Calculate mean ##
mean_Sen = np.nanmean(Sen_temperature_array)
Sen_final_mean.append(mean_Sen)
Sen_station_time_series.append(Sen_final_mean)
return Sen_station_time_series
def eliminate_SENTINEL_cloudy_data(sen_directory, shape_path):
"""
eliminates the scenes which are not cloud_free
:return:
"""
# use this function to delete data which is not overlapping with ROI
eliminate_nanoverlap(sen_directory=sen_directory, shape_path=shape_path)
# create new cloud_free directory, overwrite if already exits
cloud_free = sen_directory + "/selected/cloud_free"
if os.path.exists(cloud_free):
shutil.rmtree(cloud_free)
os.mkdir(sen_directory + "/selected/cloud_free")
selected_tifs = extract_files_to_list(path_to_folder=sen_directory + "/selected", datatype=".tif")
# import polygons from shapefile
import_list = import_polygons(shape_path=shape_path)
# loop through all .tif files in folder
for i, tifs in enumerate(selected_tifs):
src1 = rio.open(selected_tifs[i])
bool_list = []
flag = 0 #set boolean flag to 0
# loop through all weatherstations for each .tif file
for j, polygons in enumerate(import_list):
out_image1, out_transform1 = rio.mask.mask(src1, [import_list[0][j]], all_touched=1, crop=True,
nodata=np.nan)
# extract minimal temperature at all stations for each .tif file to eliminate error values
min_temp = np.min(out_image1[0])
# create boolean values for error values and store in list (-40° is lowest remotely realistic temperature)
if min_temp < -40:
bool_list.append(False)
else:
bool_list.append(True)
# check the list for flags of error values and only save tif files without any error values at station locations
for boolean in bool_list:
if boolean == False:
flag = 1
break
if flag == 0:
shutil.copy(selected_tifs[i], cloud_free)
def eliminate_nanoverlap(sen_directory, shape_path):
"""
eliminates the scenes which would'nt match with all weather stations
:return:
"""
import_list = import_polygons(shape_path=shape_path)
file_list = extract_files_to_list(path_to_folder=sen_directory, datatype=".tif")
tifs_selected = sen_directory + "/selected"
if os.path.exists(tifs_selected):
shutil.rmtree(tifs_selected)
os.mkdir(sen_directory + "/selected")
for i, files in enumerate(file_list):
src1 = rio.open(file_list[i])
try:
for j, polygons in enumerate(import_list):
rio.mask.mask(src1, [import_list[0][j]], all_touched=0, crop=True, nodata=np.nan)
shutil.copy(file_list[i], tifs_selected)
except ValueError:
pass
def extract_MODIS_temp_list(mod_directory, mod_shape_path, daytime_MODIS=None):
"""
returns list of lists
:return:
returns
"""
import_list = import_polygons(shape_path=mod_shape_path)
modis_file_list = extract_files_to_list(path_to_folder=mod_directory, datatype=".tif")
Mod_station_time_series = []
for i, polygons in enumerate(import_list):
## Initialize empty analysis lists
Mod_final_mean = []
for j, tifs in enumerate(modis_file_list):
#if daytime_MODIS in str(modis_file_list[j]):
src1 = rio.open(modis_file_list[j])
mask = rio.mask.mask(src1, [import_list[0][i]], all_touched=True, crop=True, nodata=np.nan)
Mod_temperature_array = mask[0][0]
mean_Mod = np.nanmean(Mod_temperature_array)
Mod_final_mean.append(mean_Mod)
Mod_station_time_series.append(Mod_final_mean)
return Mod_station_time_series
def analyze_SENTINEL_temperature(sen_directory, sen_shape_path, daytime_S3):
"""
returns the temperature value for each pixel of each station for each SENTINEL scene
prints an array for the scenes according to the stations
:return:
returns three arrays (1st array displayed for each station is the temperature array)
returns one temperature array for each station if a[0][0]
"""
import_list = import_polygons(shape_path=sen_shape_path)
sentinel_file_list = extract_files_to_list(path_to_folder=sen_directory, datatype=".tif")
print(import_list)
print(sentinel_file_list)
print("#################### SENTINEL RESULTS ####################")
Sen_station_mean = []
Sen_station_median = []
Sen_station_stdev = []
for i, polygons in enumerate(import_list):
scenes = i + 1
## Initialize empty analysis lists
Sen_final_mean = []
Sen_final_median = []
Sen_final_stdev = []
Sen_final_variance = []
Sen_final_percentile = []
Sen_final_range = []
# print(scenes, ". weather station")
print("{}.{}".format(i +1, station_names[i]))
for j, tifs in enumerate(sentinel_file_list):
if daytime_S3 in str(sentinel_file_list[j]):
stations = j + 1
# print(stations, ". scene")
src1 = rio.open(sentinel_file_list[j])
mask = rio.mask.mask(src1, [import_list[0][i]], all_touched=True, crop=True, nodata=np.nan)
Sen_temperature_array = mask[0][0]
# print(Sen_temperature_array)
## Calculate mean ##
mean_Sen = np.nanmean(Sen_temperature_array)
# print("Mean =" + " " + str(mean_Sen))
Sen_final_mean.append(mean_Sen)
Station_mean = np.nanmean(Sen_final_mean)
## Calculate median ##
median_Sen = np.nanmedian(Sen_temperature_array)
# print("Median =" + " " + str(median_Sen))
Sen_final_median.append(median_Sen)
Station_median = np.nanmedian(Sen_final_median)
## Calculate stdev ##
stdev_Sen = np.nanstd(Sen_temperature_array)
# print("Stdev =" + " " + str(stdev_Sen))
Sen_final_stdev.append(stdev_Sen)
Station_stdev = np.nanmedian(Sen_final_stdev)
## Calculate variance ##
var_Sen = np.nanvar(Sen_temperature_array)
# print("Variance =" + " " + str(var_Sen))
Sen_final_variance.append(var_Sen)
## Calculate percentile ##
percentile_Sen = np.nanpercentile(Sen_temperature_array, 10)
# print("Percentile =" + " " + str(percentile_Sen))
Sen_final_percentile.append(percentile_Sen)
## Calculate range ##
range_Sen = np.nanmax(Sen_temperature_array) - np.nanmin(Sen_temperature_array)
# print("Range =" + " " + str(range_Sen))
Sen_final_range.append(range_Sen)
else:
pass
print("Station " + str(i + 1) + " mean for all scenes =" + " " + str(Station_mean))
# print("Station " + str(i + 1) + " median for all scenes =" + " " + str(Station_median))
# print("Station " + str(i + 1) + " stdev for all scenes =" + " " + str(Station_stdev))
Sen_station_mean.append(Station_mean)
# Sen_station_median.append(Station_median)
# Sen_station_stdev.append(Station_stdev)
# Plot multiple means; order of scenes is fundamental; plots.py (line 118-122)
# Sen_station_mean.append(Sen_final_mean)
# Sen_station_median.append(Sen_final_median)
return Sen_station_mean
def analyze_MODIS_temperature(mod_directory, mod_shape_path, daytime_MODIS):
"""
returns the temperature value for each pixel of each station for each MODIS scene
:return:
returns one temperature array for each station
"""
import_list = import_polygons(shape_path=mod_shape_path)
modis_file_list = extract_files_to_list(path_to_folder=mod_directory, datatype=".tif")
print(len(modis_file_list))
Mod_station_mean = []
Mod_station_median = []
Mod_station_stdev = []
print("#################### MODIS RESULTS ####################")
for i, polygons in enumerate(import_list):
scenes = i+1
## Initialize empty analysis lists
Mod_final_mean = []
Mod_final_median = []
Mod_final_stdev = []
Mod_final_variance = []
Mod_final_percentile = []
Mod_final_range = []
# print(scenes, ". weather station")
print("{}.{}".format(i + 1, station_names[i]))
for j, tifs in enumerate(modis_file_list):
if daytime_MODIS in str(modis_file_list[j]):
stations = j + 1
# print(stations, ". scene")
src1 = rio.open(modis_file_list[j])
mask = rio.mask.mask(src1, [import_list[0][i]], all_touched=True, crop=True, nodata=np.nan)
Mod_temperature_array = mask[0][0]
# print(Mod_temperature_array)
## Mittelwert berechnen ##
mean_Mod = np.nanmean(Mod_temperature_array)
# print("Mean =" + " " + str(mean_Mod))
Mod_final_mean.append(mean_Mod)
Station_mean = np.nanmean(Mod_final_mean)
## Median berechnen ##
median_Mod = np.nanmedian(Mod_temperature_array)
# print("Median =" + " " + str(median_Mod))
Mod_final_median.append(median_Mod)
Station_median = np.nanmedian(Mod_final_median)
## Calculate stdev ##
stdev_Mod = np.nanstd(Mod_temperature_array)
# print("Stdev =" + " " + str(stdev_Mod))
Mod_final_stdev.append(stdev_Mod)
Station_stdev = np.nanmedian(Mod_final_stdev)
## Calculate variance ##
var_Mod = np.nanvar(Mod_temperature_array)
# print("Variance =" + " " + str(var_Mod))
Mod_final_variance.append(var_Mod)
## Calculate percentile ##
percentile_Mod = np.nanpercentile(Mod_temperature_array, 10)
# print("Percentile =" + " " + str(percentile_Mod))
Mod_final_percentile.append(percentile_Mod)
## Calculate range ##
range_Mod = np.nanmax(Mod_temperature_array) - np.nanmin(Mod_temperature_array)
# print("Range =" + " " + str(range_Mod))
Mod_final_range.append(range_Mod)
else:
pass
### activate this for old functionality ###
# print("Station " + str(i+1) + " mean for all scenes =" + " " + str(Station_mean))
# print("Station " + str(i + 1) + " median for all scenes =" + " " + str(Station_median))
# print("Station " + str(i + 1) + " stdev for all scenes =" + " " + str(Station_stdev))
Mod_station_mean.append(Station_mean)
# Mod_station_median.append(Station_median)
# Mod_station_stdev.append(Station_stdev)
# Plot multiple means; order of scenes is fundamental; plots.py (line 118-122)
# Mod_station_mean.append(Mod_final_mean)
return Mod_station_mean