a_adli/hydroquebec_archetype_gispy
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Add pairwise clip of property assessment with several overlays

Browse Source
This commit is contained in:
Alireza Adli 2024-02-26 11:17:51 -05:00
parent b0e94b114b
commit c5aade8853
1 changed files with 40 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

50
mtl_buildings_workflow.py
View File

@ -117,11 +117,11 @@ splitting_intervals = clipped_nrcan_length // num_layers
for each in range(num_layers):
output_layer_path = f'C:/Users/a_adli/PycharmProjects/hydroquebec_archetype_gispy/data/output_data/splitted_clipped_nrcan/layer_{each}/layer_{each}.shp'
params = {'INPUT': clipped_nrcan_layer,
'EXPRESSION': f'$id >= {each * splitting_intervals} AND $id < {(each + 1) * splitting_intervals}\r\n',
'OUTPUT': output_layer_path}
splitting_params = {'INPUT': clipped_nrcan_layer,
'EXPRESSION': f'$id >= {each * splitting_intervals} AND $id < {(each + 1) * splitting_intervals}\r\n',
'OUTPUT': output_layer_path}
processing.run("native:extractbyexpression", params)
processing.run("native:extractbyexpression", splitting_params)
# Creating spatian index for each new layer
create_spatial_indext_params = {'INPUT': output_layer_path, 'OUTPUT': 'Output'}
@ -129,12 +129,42 @@ for each in range(num_layers):
# Putting the remaining features in one last layer (So, overall we're going to have [num_layers + 1] layers).
remaining_features = num_layers
os.makedirs(f'C:/Users/a_adli/PycharmProjects/hydroquebec_archetype_gispy/data/aEndeavor/divided_all/layer_{remaining_features}')
output_layer_path = f'C:/Users/a_adli/PycharmProjects/hydroquebec_archetype_gispy/data/aEndeavor/divided_all/layer_{remaining_features}/layer_{remaining_features}.shp'
params = {'INPUT': clipped_nrcan_layer,
'EXPRESSION': f'$id >= {num_layers * splitting_intervals}\r\n',
'OUTPUT': output_layer_path}
os.makedirs(f'C:/Users/a_adli/PycharmProjects/hydroquebec_archetype_gispy/data/output_data/splitted_clipped_nrcan/layer_{remaining_features}')
output_layer_path = f'C:/Users/a_adli/PycharmProjects/hydroquebec_archetype_gispy/data/output_data/splitted_clipped_nrcan/layer_{remaining_features}/layer_{remaining_features}.shp'
last_splitting_params = {'INPUT': clipped_nrcan_layer,
'EXPRESSION': f'$id >= {num_layers * splitting_intervals}\r\n',
'OUTPUT': output_layer_path}
processing.run("native:extractbyexpression", params)
processing.run("native:extractbyexpression", last_splitting_params)
# Create spatial index for the last partition (layer)
create_spatial_indext_params = {'INPUT': output_layer_path, 'OUTPUT': 'Output'}
processing.run("native:createspatialindex", create_spatial_indext_params)
# Now we clip the property assessment with each of the partitions (layers)
# from the last process. The process outputs the pairwise clipped version of
# the Property Assessment data in num_layers (the number is 21 in the first run) number of individual layers
# First we make the folders for each clipping process
clipped_layers_folder = 'C:/Users/a_adli/PycharmProjects/hydroquebec_archetype_gispy/data/output_data/clipped_property_assessment_partitions'
num_clipped_layers = 21
create_folders(clipped_layers_folder, num_clipped_layers)
# Through a for loop, we clip the Property Assessment layer by all the individual overlays
for each in range(num_clipped_layers):
clipping_layer = f'C:/Users/a_adli/PycharmProjects/hydroquebec_archetype_gispy/data/output_data/splitted_clipped_nrcan/layer_{each}/layer_{each}.shp'
clipped_layer = f'C:/Users/a_adli/PycharmProjects/hydroquebec_archetype_gispy/data/output_data/clipped_property_assessment_partitions/layer_{each}/layer_{each}.shp'
pairwise_clipping_params = {'INPUT': property_assessment_layer, 'OVERLAY': clipping_layer, 'FILTER_EXPRESSION': '', 'FILTER_EXTENT': None, 'OUTPUT': clipped_layer}
processing.run("native:clip", pairwise_clipping_params)
# Creating spatian index for each new layer
create_spatial_indext_params = {'INPUT': clipped_layer, 'OUTPUT': 'Output'}
processing.run("native:createspatialindex", create_spatial_indext_params)
print("Pairwise Clipping is completed.")
# Finally, we merge all the clipped layers to make the
# pairwise clipped Property Assessment layer
qgs.exitQgis()