Je propose une approche qui ne revient qu'à un générateur de géométrie et une fonction personnalisée.
Avant de commencer, je tiens à souligner que je concentrerai l'attention sur l'explication des choses minimales à faire pour reproduire le résultat souhaité: cela signifie que certains autres paramètres mineurs (comme les tailles, les largeurs, etc.) devraient être facilement ajustés par vous pour mieux répondre à vos besoins.
Par conséquent, cette solution fonctionne à la fois pour les systèmes de référence géographiques et projetés: dans ce qui suit, j'ai supposé utiliser un CRS projeté (c'est-à-dire que les unités de mesure sont des mètres), mais vous pouvez les modifier en fonction de votre CRS.
Le contexte
Supposons que nous partions de cette couche de vecteur de ligne linéaire représentant les fils (les étiquettes représentent le nombre de fils qui se chevauchent (coïncidents)):
Solution
Tout d'abord, accédez à Layer Properties | Style
, puis choisissez le Single symbol
moteur de rendu.
Dans la Symbol selector
boîte de dialogue, choisissez un Geometry generator
type de calque de symbole et un type Linestring / MultiLinestring
de géométrie. Cliquez ensuite sur l' Function Editor
onglet:
Ensuite, cliquez sur New file
et tapez draw_wires
comme nom de la nouvelle fonction:
Vous verrez qu'une nouvelle fonction a été créée et elle est répertoriée sur le côté gauche de la boîte de dialogue. Maintenant, cliquez sur le nom de la fonction et remplacez la valeur @qgsfunction
par défaut par le code suivant (n'oubliez pas d'ajouter toutes les bibliothèques attachées ici):
from qgis.core import *
from qgis.gui import *
from math import sin, cos, radians
@qgsfunction(args='auto', group='Custom')
def draw_wires(angle, percentage, curr_feat, layer_name, feature, parent):
def wires(polyline, new_angle, percentage):
for x in range(0, len(polyline)-1):
vertices = []
first_point = polyline[x]
second_point = polyline[x +1]
seg = QgsGeometry.fromPolyline([first_point, second_point])
len_feat = seg.length()
frac_len = percentage * len_feat
limb = frac_len/cos(radians(new_angle))
tmp_azim = first_point.azimuth(second_point)
angle_1 = radians(90 - (tmp_azim+new_angle))
dist_x, dist_y = (limb * cos(angle_1), limb * sin(angle_1))
point_1 = QgsPoint(first_point[0] + dist_x, first_point[1] + dist_y)
angle_2 = radians(90 - (tmp_azim-new_angle))
dist_x, dist_y = (limb * cos(angle_2), limb * sin(angle_2))
point_2 = QgsPoint(second_point[0] - dist_x, second_point[1] - dist_y)
tmp_azim = second_point.azimuth(first_point)
angle_3 = radians(90 - (tmp_azim+new_angle))
dist_x, dist_y = (limb * cos(angle_3), limb * sin(angle_3))
point_3 = QgsPoint(second_point[0] + dist_x, second_point[1] + dist_y)
angle_4 = radians(90 - (tmp_azim-new_angle))
dist_x, dist_y = (limb * cos(angle_4), limb * sin(angle_4))
point_4 = QgsPoint(first_point[0] - dist_x, first_point[1] - dist_y)
vertices.extend([first_point, point_1, point_2, second_point, point_3, point_4, first_point])
tempGeom = QgsGeometry.fromPolyline(vertices)
num.append(tempGeom)
return num
layer = QgsMapLayerRegistry.instance().mapLayersByName(layer_name)[0]
all_feats = {}
index = QgsSpatialIndex()
for ft in layer.getFeatures():
index.insertFeature(ft)
all_feats[ft.id()] = ft
first = True
tmp_geom = curr_feat.geometry()
polyline = tmp_geom.asPolyline()
idsList = index.intersects(tmp_geom.boundingBox())
occurrences = 0
for id in idsList:
test_feat = all_feats[id]
test_geom = test_feat.geometry()
if tmp_geom.equals(test_geom):
occurrences += 1
if occurrences & 0x1:
num = [tmp_geom]
else:
num = []
rapp = occurrences/2
i=2
new_angle = angle
while i <= occurrences:
draw=wires(polyline, new_angle, percentage)
i += 2
new_angle -= new_angle/rapp
first = True
for h in num:
if first:
geom = QgsGeometry(h)
first = False
else:
geom = geom.combine(h)
return geom
Une fois que vous avez fait cela, cliquez sur le Load
bouton et vous pourrez voir la fonction dans le Custom
menu de la Expression
boîte de dialogue.
Maintenant, tapez cette expression (voir l'image ci-dessous comme référence):
draw_wires(40, 0.3, $currentfeature, @layer_name)
Vous venez d'exécuter une fonction qui dit, de manière imaginaire:
"Pour la couche actuelle ( @layer_name ) et la fonction actuelle ( $ currentfeature ), affichez les fils ensemble en utilisant une ouverture maximale initiale de 40 degrés et avec un changement de direction à une distance de 0,3 fois la longueur du segment actuel."
La seule chose que vous devez changer est la valeur des deux premiers paramètres comme vous le souhaitez, mais évidemment de manière raisonnable (laissez les autres paramètres de fonction tels qu'ils sont fournis).
Enfin, cliquez sur le Apply
bouton pour appliquer les modifications.
Vous verrez quelque chose comme ceci:
comme prévu.
ÉDITER
Selon une demande spécifique soulevée par le PO dans un commentaire:
"Serait-il possible de créer ce motif uniquement entre le début et la fin de chaque polyligne plutôt qu'entre chaque sommet?"
J'ai légèrement modifié le code. La fonction suivante doit retourner le résultat attendu:
from qgis.core import *
from qgis.gui import *
from math import sin, cos, radians
@qgsfunction(args='auto', group='Custom')
def draw_wires(angle, percentage, curr_feat, layer_name, feature, parent):
def wires(polyline, new_angle, percentage):
vertices = []
len_feat = polyline.length()
frac_len = percentage * len_feat
limb = frac_len/cos(radians(new_angle))
tmp_azim = first_point.azimuth(second_point)
angle_1 = radians(90 - (tmp_azim+new_angle))
dist_x, dist_y = (limb * cos(angle_1), limb * sin(angle_1))
point_1 = QgsPoint(first_point[0] + dist_x, first_point[1] + dist_y)
angle_2 = radians(90 - (tmp_azim-new_angle))
dist_x, dist_y = (limb * cos(angle_2), limb * sin(angle_2))
point_2 = QgsPoint(second_point[0] - dist_x, second_point[1] - dist_y)
tmp_azim = second_point.azimuth(first_point)
angle_3 = radians(90 - (tmp_azim+new_angle))
dist_x, dist_y = (limb * cos(angle_3), limb * sin(angle_3))
point_3 = QgsPoint(second_point[0] + dist_x, second_point[1] + dist_y)
angle_4 = radians(90 - (tmp_azim-new_angle))
dist_x, dist_y = (limb * cos(angle_4), limb * sin(angle_4))
point_4 = QgsPoint(first_point[0] - dist_x, first_point[1] - dist_y)
vertices.extend([first_point, point_1, point_2, second_point, point_3, point_4, first_point])
tempGeom = QgsGeometry.fromPolyline(vertices)
num.append(tempGeom)
layer = QgsMapLayerRegistry.instance().mapLayersByName(layer_name)[0]
all_feats = {}
index = QgsSpatialIndex()
for ft in layer.getFeatures():
index.insertFeature(ft)
all_feats[ft.id()] = ft
first = True
tmp_geom = curr_feat.geometry()
coords = tmp_geom.asMultiPolyline()
if coords:
new_coords = [QgsPoint(x, y) for x, y in z for z in coords]
else:
coords = tmp_geom.asPolyline()
new_coords = [QgsPoint(x, y) for x, y in coords]
first_point = new_coords[0]
second_point = new_coords[-1]
polyline=QgsGeometry.fromPolyline([first_point, second_point])
idsList = index.intersects(tmp_geom.boundingBox())
occurrences = 0
for id in idsList:
test_feat = all_feats[id]
test_geom = test_feat.geometry()
if tmp_geom.equals(test_geom):
occurrences += 1
if occurrences & 0x1:
num = [polyline]
else:
num = []
rapp = occurrences/2
i=2
new_angle = angle
while i <= occurrences:
draw=wires(polyline, new_angle, percentage)
i += 2
new_angle -= new_angle/rapp
first = True
for h in num:
if first:
geom = QgsGeometry(h)
first = False
else:
geom = geom.combine(h)
return geom