Commit dce89778 authored by Gbodjo Yawogan Jean Eudes's avatar Gbodjo Yawogan Jean Eudes
Browse files

adding new files to repo

parent d9f48542
import gdal, ogr, osr
import rasterio
import numpy as np
import os
import sys
import pandas as pd
import geopandas
import datetime
def ExtractPatch(plots_file,lstTS,dates_file,outPath,outPattern,area,plots_info=False):
'''
'''
print(outPattern,plots_info)
df = geopandas.read_file(plots_file)
lstDates = []
f = open(dates_file)
for d in f :
lstDates.append(d[:8])
size = int(np.sqrt(area * 10000))
ts_array = None
dic = {}
for j in range(len(lstTS)) :
with rasterio.open(lstTS[j]) as ds:
# band = os.path.basename(lstTS[j]).split('_',1)[0]
count = ds.count
cellsize = ds.profile['transform'][0]
originX = ds.profile['transform'][2]
originY = ds.profile['transform'][5]
crs = df.crs
epsg = int(crs['init'].split(":")[-1])
Proj = osr.SpatialReference()
Proj.ImportFromEPSG(epsg)
tmp_array = None
for index, row in df.iterrows():
centroid = row.geometry.centroid
ulx = centroid.x - size/2
uly = centroid.y + size/2
xOffset = int((ulx - originX) / cellsize)
yOffset = int((uly - originY) / -cellsize)
col_nb=row_nb=int(size/cellsize)
window = rasterio.windows.Window(xOffset,yOffset,col_nb,row_nb)
if xOffset < ds.width and xOffset+col_nb < ds.width and yOffset < ds.height and yOffset+row_nb < ds.height :
if j == 0 and plots_info :
dic.setdefault('feat_index',[]).append(index)
dic.setdefault('ID',[]).append(row['ID'])
dic.setdefault('Projet',[]).append(row['Projet'])
dic.setdefault('Biom_f',[]).append(row['Biomasse_f'])
dic.setdefault('Biom_s',[]).append(row['Biomasse_s'])
dic.setdefault('Rdt_f',[]).append(row['Rdt_f'])
dic.setdefault('Rdt_s',[]).append(row['Rdt_s'])
dic.setdefault('Sowing',[]).append(row['Semi'])
dic.setdefault('Harvest',[]).append(row['Recolte'])
ts_values = None
for c in range(count) :
if ts_values is None :
ts_values = ds.read(c+1,window=window).reshape(col_nb*row_nb)
else :
ts_values = np.hstack((ts_values,ds.read(c+1,window=window).reshape(col_nb*row_nb)))
if tmp_array is None:
tmp_array = ts_values
else:
tmp_array = np.vstack((tmp_array,ts_values))
if ts_array is None:
ts_array = tmp_array
else :
ts_array = np.hstack((ts_array,tmp_array))
outPath = "./patchs/{}Ha".format(area)
if not os.path.exists(outPath):
os.makedirs(outPath)
if plots_info :
ptrn = outPattern.split('_')[0]+'_'+outPattern.split('_')[1]
outCSV = os.path.join(outPath,'{}_plots.csv'.format(ptrn))
outdf = pd.DataFrame.from_dict(dic)
outdf = outdf.drop(columns='feat_index')
outdf.to_csv(outCSV,index=False)
np.save(os.path.join(outPath,'{}.npy'.format(outPattern)),ts_array)
if __name__=="__main__":
# 2017
# Niakhar
# Radar
plots_file = "/media/je/LACIE/SENEGAL/RENDEMENTS/Rdt_Mil_Niakhar_2017.shp"
outPath = "/media/je/SATA_1/lab5"
area=1
lstTS = ["/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-1/ASC/2017/OUTPUT/CONCAT/VH_ASC_CONCAT_S1.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-1/ASC/2017/OUTPUT/CONCAT/VV_ASC_CONCAT_S1.tif"]
dates_file = "/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-1/ASC/2017/OUTPUT/CONCAT/dates_ASC.txt"
outPattern = "niakhar_2017_radar"
ExtractPatch(plots_file,lstTS,dates_file,outPath,outPattern,area,plots_info=True)
# Optical
lstTS = ["/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2017/OUTPUT/CONCAT_GAPF/B2_NIAKHAR_2017_CONCAT_S2_GAPF.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2017/OUTPUT/CONCAT_GAPF/B3_NIAKHAR_2017_CONCAT_S2_GAPF.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2017/OUTPUT/CONCAT_GAPF/B4_NIAKHAR_2017_CONCAT_S2_GAPF.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2017/OUTPUT/CONCAT_GAPF/B8_NIAKHAR_2017_CONCAT_S2_GAPF.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2017/OUTPUT/CONCAT_GAPF/B5_NIAKHAR_2017_CONCAT_S2_GAPF-10m.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2017/OUTPUT/CONCAT_GAPF/B6_NIAKHAR_2017_CONCAT_S2_GAPF-10m.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2017/OUTPUT/CONCAT_GAPF/B7_NIAKHAR_2017_CONCAT_S2_GAPF-10m.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2017/OUTPUT/CONCAT_GAPF/B8A_NIAKHAR_2017_CONCAT_S2_GAPF-10m.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2017/OUTPUT/CONCAT_GAPF/B11_NIAKHAR_2017_CONCAT_S2_GAPF-10m.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2017/OUTPUT/CONCAT_GAPF/B12_NIAKHAR_2017_CONCAT_S2_GAPF-10m.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2017/OUTPUT/INDICES/NDVI_NIAKHAR_2017_CONCAT_S2_GAPF.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2017/OUTPUT/INDICES/NDWI_NIAKHAR_2017_CONCAT_S2_GAPF.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2017/OUTPUT/INDICES/EVI_NIAKHAR_2017_CONCAT_S2_GAPF.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2017/OUTPUT/INDICES/MSAVI2_NIAKHAR_2017_CONCAT_S2_GAPF.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2017/OUTPUT/INDICES/GDVI_NIAKHAR_2017_CONCAT_S2_GAPF.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2017/OUTPUT/INDICES/CIGreen_NIAKHAR_2017_CONCAT_S2_GAPF.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2017/OUTPUT/INDICES/CIRedEdge_NIAKHAR_2017_CONCAT_S2_GAPF.tif"]
dates_file = "/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2017/OUTPUT/CONCAT_GAPF/dates.txt"
outPattern = "niakhar_2017_opt_gapf"
ExtractPatch(plots_file,lstTS,dates_file,outPath,outPattern,area)
# 2018
# Niakhar
# Radar
outPath = "/media/je/SATA_1/lab5"
area=1
lstTS = ["/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-1/ASC/2018/OUTPUT/CONCAT/VH_ASC_CONCAT_S1.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-1/ASC/2018/OUTPUT/CONCAT/VV_ASC_CONCAT_S1.tif"]
dates_file = "/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-1/ASC/2018/OUTPUT/CONCAT/dates_ASC.txt"
plots_file = "/media/je/LACIE/SENEGAL/RENDEMENTS/Rdt_Mil_Niakhar_2018.shp"
outPattern = "niakhar_2018_radar"
ExtractPatch(plots_file,lstTS,dates_file,outPath,outPattern,area,plots_info=True)
plots_file = "/media/je/LACIE/SENEGAL/RENDEMENTS/Rdt_Mil_Niakhar-SERENA_2018.shp"
outPattern = "niakhar-serena_2018_radar"
ExtractPatch(plots_file,lstTS,dates_file,outPath,outPattern,area,plots_info=True)
plots_file = "/media/je/LACIE/SENEGAL/RENDEMENTS/Rdt_Mil_Niakhar-SIMCO_2018.shp"
outPattern = "niakhar-simco_2018_radar"
ExtractPatch(plots_file,lstTS,dates_file,outPath,outPattern,area,plots_info=True)
# Optical
lstTS = ["/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2018/OUTPUT/CONCAT_GAPF/B2_NIAKHAR_2018_CONCAT_S2_GAPF.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2018/OUTPUT/CONCAT_GAPF/B3_NIAKHAR_2018_CONCAT_S2_GAPF.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2018/OUTPUT/CONCAT_GAPF/B4_NIAKHAR_2018_CONCAT_S2_GAPF.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2018/OUTPUT/CONCAT_GAPF/B8_NIAKHAR_2018_CONCAT_S2_GAPF.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2018/OUTPUT/CONCAT_GAPF/B5_NIAKHAR_2018_CONCAT_S2_GAPF-10m.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2018/OUTPUT/CONCAT_GAPF/B6_NIAKHAR_2018_CONCAT_S2_GAPF-10m.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2018/OUTPUT/CONCAT_GAPF/B7_NIAKHAR_2018_CONCAT_S2_GAPF-10m.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2018/OUTPUT/CONCAT_GAPF/B8A_NIAKHAR_2018_CONCAT_S2_GAPF-10m.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2018/OUTPUT/CONCAT_GAPF/B11_NIAKHAR_2018_CONCAT_S2_GAPF-10m.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2018/OUTPUT/CONCAT_GAPF/B12_NIAKHAR_2018_CONCAT_S2_GAPF-10m.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2018/OUTPUT/INDICES/NDVI_NIAKHAR_2018_CONCAT_S2_GAPF.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2018/OUTPUT/INDICES/NDWI_NIAKHAR_2018_CONCAT_S2_GAPF.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2018/OUTPUT/INDICES/EVI_NIAKHAR_2018_CONCAT_S2_GAPF.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2018/OUTPUT/INDICES/MSAVI2_NIAKHAR_2018_CONCAT_S2_GAPF.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2018/OUTPUT/INDICES/GDVI_NIAKHAR_2018_CONCAT_S2_GAPF.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2018/OUTPUT/INDICES/CIGreen_NIAKHAR_2018_CONCAT_S2_GAPF.tif",
"/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2018/OUTPUT/INDICES/CIRedEdge_NIAKHAR_2018_CONCAT_S2_GAPF.tif"]
dates_file = "/media/je/LACIE/SENEGAL/NIAKHAR/SENTINEL-2/2018/OUTPUT/CONCAT_GAPF/dates.txt"
plots_file = "/media/je/LACIE/SENEGAL/RENDEMENTS/Rdt_Mil_Niakhar_2018.shp"
outPattern = "niakhar_2018_opt_gapf"
ExtractPatch(plots_file,lstTS,dates_file,outPath,outPattern,area)
plots_file = "/media/je/LACIE/SENEGAL/RENDEMENTS/Rdt_Mil_Niakhar-SERENA_2018.shp"
outPattern = "niakhar-serena_2018_opt_gapf"
ExtractPatch(plots_file,lstTS,dates_file,outPath,outPattern,area)
plots_file = "/media/je/LACIE/SENEGAL/RENDEMENTS/Rdt_Mil_Niakhar-SIMCO_2018.shp"
outPattern = "niakhar-simco_2018_opt_gapf"
ExtractPatch(plots_file,lstTS,dates_file,outPath,outPattern,area)
# 2018
# Nioro
# Radar
plots_file = "/media/je/LACIE/SENEGAL/RENDEMENTS/Rdt_Mil_Nioro_2018.shp"
outPath = "/media/je/SATA_1/lab5"
area=1
lstTS = ["/media/je/LACIE/SENEGAL/NIORO/SENTINEL-1/ASC/2018/OUTPUT/CONCAT/VH_ASC_CONCAT_S1.tif",
"/media/je/LACIE/SENEGAL/NIORO/SENTINEL-1/ASC/2018/OUTPUT/CONCAT/VV_ASC_CONCAT_S1.tif"]
dates_file = "/media/je/LACIE/SENEGAL/NIORO/SENTINEL-1/ASC/2018/OUTPUT/CONCAT/dates_ASC.txt"
outPattern = "nioro_2018_radar"
ExtractPatch(plots_file,lstTS,dates_file,outPath,outPattern,area,plots_info=True)
# Optical
lstTS = ["/media/je/LACIE/SENEGAL/NIORO/SENTINEL-2/2018/OUTPUT/CONCAT_GAPF/B2_NIORO_2018_CONCAT_S2_GAPF.tif",
"/media/je/LACIE/SENEGAL/NIORO/SENTINEL-2/2018/OUTPUT/CONCAT_GAPF/B3_NIORO_2018_CONCAT_S2_GAPF.tif",
"/media/je/LACIE/SENEGAL/NIORO/SENTINEL-2/2018/OUTPUT/CONCAT_GAPF/B4_NIORO_2018_CONCAT_S2_GAPF.tif",
"/media/je/LACIE/SENEGAL/NIORO/SENTINEL-2/2018/OUTPUT/CONCAT_GAPF/B8_NIORO_2018_CONCAT_S2_GAPF.tif",
"/media/je/LACIE/SENEGAL/NIORO/SENTINEL-2/2018/OUTPUT/CONCAT_GAPF/B5_NIORO_2018_CONCAT_S2_GAPF-10m.tif",
"/media/je/LACIE/SENEGAL/NIORO/SENTINEL-2/2018/OUTPUT/CONCAT_GAPF/B6_NIORO_2018_CONCAT_S2_GAPF-10m.tif",
"/media/je/LACIE/SENEGAL/NIORO/SENTINEL-2/2018/OUTPUT/CONCAT_GAPF/B7_NIORO_2018_CONCAT_S2_GAPF-10m.tif",
"/media/je/LACIE/SENEGAL/NIORO/SENTINEL-2/2018/OUTPUT/CONCAT_GAPF/B8A_NIORO_2018_CONCAT_S2_GAPF-10m.tif",
"/media/je/LACIE/SENEGAL/NIORO/SENTINEL-2/2018/OUTPUT/CONCAT_GAPF/B11_NIORO_2018_CONCAT_S2_GAPF-10m.tif",
"/media/je/LACIE/SENEGAL/NIORO/SENTINEL-2/2018/OUTPUT/CONCAT_GAPF/B12_NIORO_2018_CONCAT_S2_GAPF-10m.tif",
"/media/je/LACIE/SENEGAL/NIORO/SENTINEL-2/2018/OUTPUT/INDICES/NDVI_NIORO_2018_CONCAT_S2_GAPF.tif",
"/media/je/LACIE/SENEGAL/NIORO/SENTINEL-2/2018/OUTPUT/INDICES/NDWI_NIORO_2018_CONCAT_S2_GAPF.tif",
"/media/je/LACIE/SENEGAL/NIORO/SENTINEL-2/2018/OUTPUT/INDICES/EVI_NIORO_2018_CONCAT_S2_GAPF.tif",
"/media/je/LACIE/SENEGAL/NIORO/SENTINEL-2/2018/OUTPUT/INDICES/MSAVI2_NIORO_2018_CONCAT_S2_GAPF.tif",
"/media/je/LACIE/SENEGAL/NIORO/SENTINEL-2/2018/OUTPUT/INDICES/GDVI_NIORO_2018_CONCAT_S2_GAPF.tif",
"/media/je/LACIE/SENEGAL/NIORO/SENTINEL-2/2018/OUTPUT/INDICES/CIGreen_NIORO_2018_CONCAT_S2_GAPF.tif",
"/media/je/LACIE/SENEGAL/NIORO/SENTINEL-2/2018/OUTPUT/INDICES/CIRedEdge_NIORO_2018_CONCAT_S2_GAPF.tif"]
dates_file = "/media/je/LACIE/SENEGAL/NIORO/SENTINEL-2/2018/OUTPUT/CONCAT_GAPF/dates.txt"
outPattern = "nioro_2018_opt_gapf"
ExtractPatch(plots_file,lstTS,dates_file,outPath,outPattern,area)
\ No newline at end of file
This source diff could not be displayed because it is too large. You can view the blob instead.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
import pandas as pd
import numpy as np
import os
import matplotlib.pyplot as plt
def norm_patch (radar_file, opt_file, plots_file, area, norm="meanstd", radar_bands=["VH","VV"],
opt_bands=["B2","B3","B4","B8","B5","B6","B7","B8A","B11","B12"],
opt_indices = ["NDVI","NDWI","EVI","MSAVI2","GDVI","CIGreen","CIRedEdge"]):
size = int(int(int(np.sqrt(area * 10000))/10)**2)
print ("Size:",size)
# Radar
radar_array = np.load(radar_file)
n_bands = len(radar_bands)
n_timestamps = int(radar_array.shape[1]/(size*n_bands))
norm_radar_array = None
for b in range(n_bands) :
ts_values = radar_array[:,b*size*n_timestamps:b*size*n_timestamps+size*n_timestamps]
if norm == "meanstd":
ts_values = (ts_values - ts_values.mean()) / ts_values.std()
elif norm == "minmax":
ts_values = (ts_values - ts_values.min()) / (ts_values.max() - ts_values.min())
if norm_radar_array is None :
norm_radar_array = ts_values
else :
norm_radar_array = np.hstack((norm_radar_array,ts_values))
radar_patchs = []
for i in range(n_timestamps):
lst = []
for j in range(n_bands):
lst.append(norm_radar_array[:,i*size+j*n_timestamps*size:i*size+j*n_timestamps*size+size])
radar_patchs.append(np.stack(lst, axis=-1))
radar_patchs = np.stack(radar_patchs,axis=-1)
print (radar_patchs.shape)
# Opt Bands
opt_array = np.load(opt_file)
n_bands = len(opt_bands)
n_timestamps = int(opt_array.shape[1]/(size*(len(opt_bands) + len(opt_indices))))
norm_opt_array = None
for b in range(len(opt_bands)) :
ts_values = opt_array[:,b*size*n_timestamps:b*size*n_timestamps+size*n_timestamps]
if norm == "meanstd":
ts_values = (ts_values - ts_values.mean()) / ts_values.std()
elif norm == "minmax":
ts_values = (ts_values - ts_values.min()) / (ts_values.max() - ts_values.min())
if norm_opt_array is None :
norm_opt_array = ts_values
else :
norm_opt_array = np.hstack((norm_opt_array,ts_values))
opt_patchs = []
for i in range(n_timestamps):
lst = []
for j in range(n_bands):
lst.append(norm_opt_array[:,i*size+j*n_timestamps*size:i*size+j*n_timestamps*size+size])
opt_patchs.append(np.stack(lst, axis=-1))
opt_patchs = np.stack(opt_patchs,axis=-1)
print (opt_patchs.shape)
# Opt Indices
indices_array = np.load(opt_file)
n_bands = len(opt_indices)
n_timestamps = int(opt_array.shape[1]/(size*(len(opt_bands) + len(opt_indices))))
norm_indices_array = None
for b in range(len(opt_bands),len(opt_bands)+len(opt_indices)) :
ts_values = indices_array[:,b*size*n_timestamps:b*size*n_timestamps+size*n_timestamps]
if norm == "meanstd":
ts_values = (ts_values - ts_values.mean()) / ts_values.std()
elif norm == "minmax":
ts_values = (ts_values - ts_values.min()) / (ts_values.max() - ts_values.min())
if norm_indices_array is None :
norm_indices_array = ts_values
else :
norm_indices_array = np.hstack((norm_indices_array,ts_values))
indices_patchs = []
for i in range(n_timestamps):
lst = []
for j in range(n_bands):
lst.append(norm_indices_array[:,i*size+j*n_timestamps*size:i*size+j*n_timestamps*size+size])
indices_patchs.append(np.stack(lst, axis=-1))
indices_patchs = np.stack(indices_patchs,axis=-1)
print (indices_patchs.shape)
outPath = "./patchs_data/{}Ha".format(area)
if not os.path.exists(outPath):
os.makedirs(outPath)
ptrn = os.path.basename(radar_file).split('_')[0]+"_"+os.path.basename(radar_file).split('_')[1]
np.save(os.path.join(outPath,"{}_radar_patchs.npy".format(ptrn)),radar_patchs)
np.save(os.path.join(outPath,"{}_opt_patchs.npy".format(ptrn)),opt_patchs)
np.save(os.path.join(outPath,"{}_indices_patchs.npy".format(ptrn)),indices_patchs)
plots_df = pd.read_csv(plots_file)
rdt = plots_df[["Biom_f","Biom_s","Rdt_f","Rdt_s"]].values
print (rdt.shape)
np.save(os.path.join(outPath,"{}_yields.npy".format(ptrn)),rdt)
if __name__ == '__main__' :
area=1
# Niakhar 2017
radar_file = "./patchs/%sHa/niakhar_2017_radar.npy"%area
opt_file = "./patchs/%sHa/niakhar_2017_opt_gapf.npy"%area
plots_file = "./patchs/%sHa/niakhar_2017_plots.csv"%area
norm_patch(radar_file,opt_file,plots_file,area,norm="minmax")
# Niakhar 2018
radar_file = "./patchs/%sHa/niakhar_2018_radar.npy"%area
opt_file = "./patchs/%sHa/niakhar_2018_opt_gapf.npy"%area
plots_file = "./patchs/%sHa/niakhar_2018_plots.csv"%area
norm_patch(radar_file,opt_file,plots_file,area,norm="minmax")
radar_file = "./patchs/%sHa/niakhar-serena_2018_radar.npy"%area
opt_file = "./patchs/%sHa/niakhar-serena_2018_opt_gapf.npy"%area
plots_file = "./patchs/%sHa/niakhar-serena_2018_plots.csv"%area
norm_patch(radar_file,opt_file,plots_file,area,norm="minmax")
radar_file = "./patchs/%sHa/niakhar-simco_2018_radar.npy"%area
opt_file = "./patchs/%sHa/niakhar-simco_2018_opt_gapf.npy"%area
plots_file = "./patchs/%sHa/niakhar-simco_2018_plots.csv"%area
norm_patch(radar_file,opt_file,plots_file,area,norm="minmax")
# Nioro 2018
radar_file = "./patchs/%sHa/nioro_2018_radar.npy"%area
opt_file = "./patchs/%sHa/nioro_2018_opt_gapf.npy"%area
plots_file = "./patchs/%sHa/nioro_2018_plots.csv"%area
norm_patch(radar_file,opt_file,plots_file,area,norm="minmax")
\ No newline at end of file
This diff is collapsed.
This diff is collapsed.
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment