Commit 5c62f9e8 authored by Laura LINDEPERG's avatar Laura LINDEPERG
Browse files

complete geomorpho code to extract all indicators

parent 18388e95
......@@ -25,10 +25,10 @@ ZHp_foldername = 'C:/Users/laura.lindeperg/Documents/DonneesLaura/Watersheds/GEO
## BDTopage
# topage_foldername = 'C:/Users/laura.lindeperg/Documents/DonneesLaura/Watersheds/GEOMORPHO/BDTopage/'
topage_foldername = 'C:/Users/laura.lindeperg/Documents/DonneesLaura/Watersheds/GEOMORPHO/BDTopage/'
## BDCarthage
# carthage_foldername = 'C:/Users/laura.lindeperg/Documents/DonneesLaura/Watersheds/GEOMORPHO/BDCarthage/'
carthage_foldername = 'C:/Users/laura.lindeperg/Documents/DonneesLaura/Watersheds/GEOMORPHO/BDCarthage/'
......@@ -49,12 +49,12 @@ for i in code_for_test:
# for i in watershed_code:
watershed_contour_i = gpd.read_file(shp_foldername+i+'.shp')
ZHp_i = rxr.open_rasterio(ZHp_foldername+i+'_ZHp.tif')
# hydro_sections_Topage_i = gpd.read_file(topage_foldername+i+'_BDTopage.shp')
# hydro_sections_Carthage_i = gpd.read_file(carthage_foldername+i+'_BDCarthage.shp')
hydro_sections_Topage_i = gpd.read_file(topage_foldername+i+'_BDTopage.shp')
hydro_sections_Carthage_i = gpd.read_file(carthage_foldername+i+'_BDCarthage.shp')
geomorpho_i = GeomorphologicProperties(i)
geomorpho_i.compute_ZHp_ratio(ZHp_i)
# geomorpho_i.extract_geomorphologic_properties(watershed_contour_i, ZHp_i, hydro_sections_Topage_i, hydro_sections_Carthage_i)
# geomorpho_i.compute_ZHp_ratio(ZHp_i)
geomorpho_i.extract_geomorphologic_properties(watershed_contour_i, ZHp_i, hydro_sections_Topage_i, hydro_sections_Carthage_i)
geomorpho_indicators = geomorpho_indicators.append(geomorpho_i.__dict__, ignore_index=True)
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