Lorsque les points se chevauchent, il existe cette propriété qui permet d'afficher automatiquement le lot d'entre eux séparément autour de l'endroit où ils se trouvent, appelé «déplacement de point». Mais cela ne fonctionne pas pour les lignes, même si cela me semble tout à fait conceptuellement réalisable afin de réaliser quelque chose comme ça:

J'ai absolument besoin de voir les différentes lignes qui sont en réalité toutes au même endroit (je travaille dans les réseaux de télécommunications). La seule façon que je vois pour l'instant est de créer vraiment des lignes différentes comme dans l'image ci-dessus, créant ainsi des erreurs spatiales.

J'utilise QGIS 2.14.

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 symbolmoteur de rendu.

Dans la Symbol selectorboîte de dialogue, choisissez un Geometry generatortype de calque de symbole et un type Linestring / MultiLinestringde géométrie. Cliquez ensuite sur l' Function Editoronglet:

Ensuite, cliquez sur New fileet tapez draw_wirescomme 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 @qgsfunctionpar 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 Loadbouton et vous pourrez voir la fonction dans le Custommenu de la Expressionboî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 Applybouton 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