2024-07-23 10:26:08 -04:00
|
|
|
from scrub_layer_class import *
|
|
|
|
|
from varennes_io_paths import *
|
|
|
|
|
|
|
|
|
|
# Making folders for the output data layers
|
|
|
|
|
create_output_folders(output_paths, output_paths_dir)
|
|
|
|
|
|
|
|
|
|
# Initialize the input data layers
|
2024-07-23 15:14:51 -04:00
|
|
|
geo_index = ScrubLayer(qgis_path, input_paths['GeoIndex'], 'GeoIndex')
|
|
|
|
|
varennes_boundary = \
|
|
|
|
|
ScrubLayer(qgis_path, input_paths['Varennes Boundary'], 'Varennes Boundary')
|
2024-07-23 17:22:07 -04:00
|
|
|
|
|
|
|
|
|
|
|
|
|
print('Processing the GeoIndex layer')
|
|
|
|
|
print(geo_index)
|
|
|
|
|
geo_index.create_spatial_index()
|
|
|
|
|
|
|
|
|
|
geo_index.fix_geometries(output_paths['Fixed GeoIndex'])
|
|
|
|
|
|
|
|
|
|
# Defining a new layer for the fixed GeoIndex
|
|
|
|
|
geo_index_fixed = ScrubLayer(qgis_path, output_paths['Fixed GeoIndex'],
|
|
|
|
|
'Fixed GeoIndex')
|
|
|
|
|
geo_index_fixed.create_spatial_index()
|
|
|
|
|
print(geo_index_fixed)
|
|
|
|
|
geo_index_fixed.clip_layer(varennes_boundary.layer_path,
|
|
|
|
|
output_paths['Clipped Fixed GeoIndex'])
|
|
|
|
|
|
|
|
|
|
geo_index_clipped = \
|
|
|
|
|
ScrubLayer(qgis_path,
|
|
|
|
|
output_paths['Clipped Fixed GeoIndex'], 'Clipped Fixed GeoIndex')
|
|
|
|
|
geo_index_clipped.create_spatial_index()
|
|
|
|
|
print(geo_index_clipped)
|
