im new to python and have been doing some work on geospatial data as part of Uni research. I'm using a collection of modules but mainly Pandas and Geopandas if that helps
I have made a function that works perfectly on what I want it to do
def HydrogenStorageCavernPoints(Rock_depth, Rock_thickness, Cave_spacing, Facies, Infrastructure):
Stage1 = gpd.clip(Rock_depth, Rock_thickness) #Clip Rock depth and thickness params
Stage2 = gpd.clip(Stage1, Facies) # Put facies into the mix
Stage3 = gpd.clip(Stage2, Infrastructure) #clips for infrastructure
Stage4 = Stage3.explode() #This expands the Ggpd df
Stage5 = Stage4[Stage4.geom_type == 'Polygon'] #Makes it only polys in the GPD df
Stage6 = gpd.clip(Cave_spacing, Stage5) #more clips
Stage7 = add_cave_volumetrics(Stage6)
Return Stage7
Here's the bit i'm struggling with.
Now for my input variables, i have multiple shape files I.e 5 shape files for different iterations of Rock_depth, 4 shape files for Iterations of Rock_Thickness, 6 for cave spacing
Im trying to determine how the changing of these varibles, for which i have different shapefiles for change with one another. So i need to run this function on all possible combinations of input shape files and hopefully save them somewhere as an unique varible i could call up to look at or something similar.
I'm hoping i've made my newby ness clear and would be very grateful if someone could help me out. Having googled for a couple of hours i can't quite figure out the best way to do this
Thanks in advance!!!