forked from zanfire/qgis-utils
-
Notifications
You must be signed in to change notification settings - Fork 0
/
qgis_layers.py
108 lines (96 loc) · 3.75 KB
/
qgis_layers.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
from qgis.utils import iface
from qgis.core import QgsVectorLayer, QgsField, QgsMapLayerRegistry
from PyQt4.QtCore import *
from PyQt4.QtGui import *
LAYER_MAX_EXTENSION = 'Estratto_MassimaEstensione'
LAYER_VOLUMES = 'Estratto_UnitaVolumetrica-Pulizia'
LAYER_NEIGHBORS = 'Neighbors'
def create_neighbors_layer():
vl = QgsVectorLayer("Polygon", LAYER_NEIGHBORS, "memory")
pr = vl.dataProvider()
vl.startEditing()
pr.addAttributes( [ QgsField(FIELD_UUID, QVariant.String),
QgsField(FIELD_NEIGHBORS_UUID, QVariant.String),
QgsField(FIELD_SV, QVariant.Double)] )
vl.updateFields()
return vl
def get_layer(name):
layers = iface.legendInterface().layers()
for layer in layers:
if layer.name() == name:
return layer
return None
def build_neighbors():
layer_max_ext = get_layer(LAYER_MAX_EXTENSION)
layer_vols = get_layer(LAYER_VOLUMES)
layer_neighbors = create_neighbors_layer()
features_max_exts = load_feature_dict_id(layer_max_ext)
features_vols = load_feature_dict_id(layer_vols)
index = build_spatialindex(features_vols)
new_features = []
for f in features_max_exts.values():
ids = index.intersects(f.geometry().boundingBox())
neighbors = []
fet = QgsFeature()
fet.setGeometry(f.geometry())
for i in ids:
if f.geometry().contains(features_vols[i].geometry()):
neighbors.append(features_vols[i][FIELD_UUID])
#fet[FIELD_NEIGHBORS_UUID] = ','.join(neighbors)
fet.setAttributes([ neighbors[0], ','.join(neighbors), 0.0])
new_features.append(fet)
layer_neighbors.dataProvider().addFeatures(new_features)
layer_neighbors.commitChanges()
QgsMapLayerRegistry.instance().addMapLayer(layer_neighbors)
def show_neighbors(skip):
layer_neighbors = get_layer(LAYER_NEIGHBORS)
layer_vols = get_layer(LAYER_VOLUMES)
features = load_feature_dict_id(layer_neighbors)
features_uuid = load_feature_dict_UUID(layer_vols)
idx = 0
for f in features.values():
# Skip first elements.
if idx < skip:
idx += 1
continue
neighbors = f[FIELD_NEIGHBORS_UUID].split(',')
selection = []
for x in neighbors:
selection.append(features_uuid[x].id())
layer_vols.setSelectedFeatures(selection)
box = layer_vols.boundingBoxOfSelected()
iface.mapCanvas().setExtent(box)
iface.mapCanvas().refresh()
reply = QMessageBox.question(None, 'Message', 'Continue?\n (seeing feature ' + str(idx)+ ')', QMessageBox.No, QMessageBox.Yes)
if reply == QMessageBox.No:
break
idx += 1
def show_features(layer, features, neighbors):
selection = []
for x in neighbors:
selection.append(features[x].id())
layer.setSelectedFeatures(selection)
box = layer.boundingBoxOfSelected()
iface.mapCanvas().setExtent(box)
iface.mapCanvas().refresh()
def compute_SVs():
layer_neighbors = get_layer(LAYER_NEIGHBORS)
layer_vols = get_layer(LAYER_VOLUMES)
features = load_feature_dict_id(layer_neighbors)
features_uuid = load_feature_dict_UUID(layer_vols)
layer_neighbors.startEditing()
try:
for f in features.values():
neighbors = f[FIELD_NEIGHBORS_UUID].split(',')
SV = 0.0
if len(neighbors) == 1:
SV = compute_SV(features_uuid[neighbors[0]])
else:
SV = compute_SVmultiple(features_uuid, neighbors)
show_features(layer_vols, features_uuid, neighbors)
f[FIELD_SV] = SV
layer_neighbors.updateFeature(f)
except:
layer_neighbors.rollBack()
raise
layer_neighbors.commitChanges()