Fix and clip NRCan, beginning Geoindex layer process

This commit is contained in:
Alireza Adli 2024-02-13 15:29:31 -05:00
parent ec0448685e
commit d1c8f24cdf

View File

@ -4,7 +4,7 @@ from services_scripts.load_layer import load_layer
import processing
# Set the path to QGIS installation
QgsApplication.setPrefixPath("C:/Program Files/QGIS 3.34.1/apps/qgis", True)
QgsApplication.setPrefixPath('C:/Program Files/QGIS 3.34.1/apps/qgis', True)
# Initialize QGIS application
qgs = QgsApplication([], False)
@ -25,15 +25,32 @@ print(f'Fixing {nrcan_name} geometries')
fixed_nrcan = 'C:/Users/a_adli/PycharmProjects/hydroquebec_archetype_gispy/data/tests/python_fixed_04/py_fixes_04.shp'
params_fixing_nrcan = {'INPUT': nrcan, 'METHOD': 0, 'OUTPUT': fixed_nrcan}
fixed_layer = processing.run("native:fixgeometries", params_fixing_nrcan)['OUTPUT']
fixed_layer_nrcan = processing.run('native:fixgeometries', params_fixing_nrcan)['OUTPUT']
fixed_nrcan, fixed_nrcan_name = load_layer(fixed_layer, 'Fixed NRCan')
fixed_nrcan, fixed_nrcan_name = load_layer(fixed_layer_nrcan, 'Fixed NRCan')
print(f'{fixed_nrcan_name} data count: {fixed_nrcan.featureCount()}')
params_create_index_nrcan = {'INPUT': fixed_layer, 'OUTPUT': 'Output'}
montreal_boundary_layer = 'C:/Users/a_adli/PycharmProjects/hydroquebec_archetype_gispy/data/montreal_boundary|layername=Montreal_boundary'
clipped_nrcan_layer = 'C:/Users/a_adli/PycharmProjects/hydroquebec_archetype_gispy/data/tests/clipped_nrcan/clipped_nrcan.shp'
params_clipping_nrcan = {'INPUT': fixed_nrcan, 'OVERLAY': montreal_boundary_layer, 'FILTER_EXPRESSION': '', 'FILTER_EXTENT': None, 'OUTPUT': clipped_layer}
processing.run("native:clip", params_clipping_nrcan)
print(f'Clipping of {fixed_nrcan_name} is completed.')
clipped_nrcan, clipped_nrcan_name = load_layer(clipped_nrcan_layer, 'Clipped NRCan')
print(f'{clipped_nrcan_name} data count: {clipped_nrcan.featureCount()}')
params_create_index_nrcan = {'INPUT': clipped_nrcan, 'OUTPUT': 'Output'}
indexed_layer = processing.run("native:createspatialindex", params_create_index_nrcan)
print(f'Creating spatial index for {clipped_nrcan_name} is completed.')
print("Creating spatial index is completed.")
geoindex_unclipped_layer = 'C:/Users/a_adli/PycharmProjects/hydroquebec_archetype_gispy/data/Geoindex_81670/mamh_usage_predo_2022_s_poly.shp'
geoindex_unclipped, geoindex_unclipped_name = load_layer(geoindex_unclipped_layer, 'GeoIndex Unclipped')
print(f'{geoindex_unclipped_name} data count: {geoindex_unclipped.featureCount()}')
qgs.exitQgis()