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 geo_index = ScrubLayer(qgis_path, input_paths['GeoIndex'], 'GeoIndex') varennes_boundary = \ ScrubLayer(qgis_path, input_paths['Varennes Boundary'], 'Varennes Boundary') property_assessment_2022 = ScrubLayer( qgis_path, input_paths['Property Assessment 2022'], 'Property Assessment 2022') 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) print('Processing the Property Assessment 2022 layer') print(property_assessment_2022) property_assessment_2022.create_spatial_index() property_assessment_2022.clip_layer( varennes_boundary.layer_path, output_paths['Varennes Clipped Property Assessment']) varennes_clipped_property_assessment = ScrubLayer( qgis_path, output_paths['Varennes Clipped Property Assessment'], 'Varennes Clipped Property Assessment') varennes_clipped_property_assessment.create_spatial_index() print(varennes_clipped_property_assessment)