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

Complete the workflow

Browse Source
This commit is contained in:
Alireza Adli 2024-08-27 13:30:33 -04:00
parent e58710e892
commit bd2a523f46
1 changed files with 30 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

36
varennes_single_processes/add_height_to_nrcan.py
View File

@ -37,11 +37,35 @@ varennes_nrcan.delete_field('min_x')
varennes_nrcan.delete_field('min_y') varennes_nrcan.delete_field('min_y')
varennes_nrcan.delete_field('max_x') varennes_nrcan.delete_field('max_x')
varennes_nrcan.delete_field('max_y') varennes_nrcan.delete_field('max_y')
varennes_nrcan.delete_field('comment')
varennes_nrcan.layer.startEditing()
# varennes_nrcan_duplicate.delete_field('heightmin') # Add a new field for the height
# varennes_nrcan_duplicate.delete_field('heightmax') weighted_height = QgsField('building_height', QVariant.Double)
# varennes_nrcan_duplicate.delete_field('elevmin') varennes_nrcan.layer.dataProvider().addAttributes([weighted_height])
# varennes_nrcan_duplicate.delete_field('elevmax') varennes_nrcan.layer.updateFields()
# varennes_nrcan_duplicate.delete_field('bldgarea')
# varennes_nrcan_duplicate.delete_field('comment') height_min_idx = varennes_nrcan.layer.fields().indexOf('heightmin')
height_max_idx = varennes_nrcan.layer.fields().indexOf('heightmax')
weighted_height_idx = varennes_nrcan.layer.fields().indexOf('ave_height')
# Iterate through each feature and calculate the average weighted height
for feature in varennes_nrcan.layer.getFeatures():
height_min = feature[height_min_idx]
height_max = feature[height_max_idx]
if height_min is not None and height_max is not None:
average_value = (2/3 * height_min + 1/3 * height_max)
else:
average_value = None
# Update the new field with the calculated average
varennes_nrcan.layer.changeAttributeValue(
feature.id(), weighted_height_idx, average_value)
# Commit the changes
varennes_nrcan.layer.commitChanges()
print('Average weighted height calculation (ave_height)'
' for buildings completed successfully!')