Commit c2078f61 authored by Laura LINDEPERG's avatar Laura LINDEPERG
Browse files

Saturated water content maps (for drainage porosity)

parent f39a2fb7
......@@ -67,6 +67,10 @@ BDLisa_shp = gpd.read_file(BDLisa_shp_path)
BRGM_geol_path = 'C:/Users/laura.lindeperg/Documents/DonneesLaura/CarteGeolBRGM/FR_vecteur/FR_vecteur/GEO001M_CART_FR_S_FGEOL_2154.shp'
BRGM_shp = gpd.read_file(BRGM_geol_path)
# SATURATED - Water retention of topsoil: saturated water content - cm3/cm3 no units
SAT_map_path = 'C:/Users/laura.lindeperg/Documents/DonneesLaura/FS_water_2016_07/ths_fao_octop.tif'
SAT_map = rxr.open_rasterio(SAT_map_path).squeeze()
# KS - Saturated hydraulic conductivity of topsoil - log10(cm/day)
KS_map_path = 'C:/Users/laura.lindeperg/Documents/DonneesLaura/FS_water_2016_07/ks_fao_octop.tif'
KS_map = rxr.open_rasterio(KS_map_path).squeeze()
......@@ -170,7 +174,7 @@ watershed_code = shp_watersheds.loc[:,'Code']
# watershed_code_IV = shp_BV_IV.loc[shp_BV_IV.loc[:,'Code'].isin(watershed_code)].loc[:, 'Code']
# Get a sample of them for test
code_for_test = watershed_code.loc[580:]
code_for_test = watershed_code.loc[0:1]
# code_for_test = ['K9341810']
# code_for_test = ['J4742020']
......@@ -198,6 +202,10 @@ for i in watershed_code:
shp_contour_i = gpd.read_file(shp_foldername+i+'.shp')
## Get hydraulic properties data within the watershed and save them in file
# SATURATED - Water retention of topsoil: saturated water content - cm3/cm3 no units
raster_SAT = SAT_map.rio.clip(shp_contour_i.geometry.apply(mapping), shp_contour_i.crs)
raster_SAT.rio.to_raster('C:/Users/laura.lindeperg/Documents/DonneesLaura/Watersheds/GEOL/JRC-ESDAC/SAT/'+i+'_SAT.tif')
# KS - Saturated hydraulic conductivity of topsoil - log10(cm/day)
raster_KS = KS_map.rio.clip(shp_contour_i.geometry.apply(mapping), shp_contour_i.crs)
raster_KS.rio.to_raster('C:/Users/laura.lindeperg/Documents/DonneesLaura/Watersheds/GEOL/JRC-ESDAC/KS/'+i+'_KS.tif')
......
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