Kiedy punkty się pokrywają, istnieje ta właściwość, która pozwala automatycznie wyświetlić wiele z nich osobno wokół ich miejsca, zwana „przesunięciem punktu”. Ale to nie działa dla linii, mimo to wydaje mi się całkiem wykonalne koncepcyjnie, aby osiągnąć coś takiego:

Absolutnie muszę zobaczyć różne linie, które w rzeczywistości wszystkie są w tym samym miejscu (pracuję w sieciach telekomunikacyjnych). Jedyny sposób, jaki widzę na razie, to naprawdę tworzenie różnych linii, jak na powyższym obrazku, co powoduje błędy przestrzenne.

Używam QGIS 2.14.

Proponuję podejście, które powtarza się tylko w przypadku generatora geometrii i funkcji niestandardowej.

Przed rozpoczęciem chcę podkreślić, że skupię uwagę na wyjaśnieniu minimalnych czynności, które należy wykonać, aby odtworzyć pożądany rezultat: oznacza to, że niektóre inne drobne parametry (takie jak rozmiary, szerokości itp.) Powinny być łatwo dostosowane przez ciebie dla lepszego dopasowania do twoich potrzeb.

Dlatego to rozwiązanie działa zarówno dla Geograficznego, jak i Projekcyjnego Systemu Referencyjnego: poniżej założyłem, że użyję rzutowanego CRS (tj. Jednostki miary to metry), ale możesz je zmienić zgodnie ze swoim CRS.

Kontekst

Załóżmy, że zaczniemy od tej wektorowej warstwy linii przedstawiającej przewody (etykiety reprezentują liczbę nakładających się (zbieżnych) drutów):

Rozwiązanie

Najpierw przejdź do, Layer Properties | Stylea następnie wybierz Single symbolrenderer.

W Symbol selectoroknie dialogowym wybierz Geometry generatortyp warstwy symbolu i Linestring / MultiLinestringtyp geometrii. Następnie kliknij Function Editorkartę:

Następnie kliknij New filei wpisz draw_wiresjako nazwę nowej funkcji:

Zobaczysz, że nowa funkcja została utworzona i jest wymieniona po lewej stronie okna dialogowego. Teraz kliknij nazwę funkcji i zastąp domyślną @qgsfunctionnastępującym kodem (nie zapomnij dodać wszystkich dołączonych tutaj bibliotek):

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

Gdy to zrobisz, kliknij Loadprzycisk, a będziesz mógł zobaczyć funkcję z Custommenu Expressionokna dialogowego.

Teraz wpisz to wyrażenie (patrz obrazek poniżej jako odniesienie):

draw_wires(40, 0.3, $currentfeature, @layer_name)

Właśnie uruchomiłeś funkcję, która w wyobrażony sposób mówi:

„W przypadku bieżącej warstwy ( @nazwa_layera ) i bieżącej funkcji ( $ currentfeature ) wyświetl przewody razem, używając początkowego maksymalnego otwarcia 40 stopni i ze zmianą kierunku w odległości 0,3-krotnej długości bieżącego segmentu.”

Jedyną rzeczą, którą musisz zmienić, jest wartość pierwszych dwóch parametrów, jak chcesz, ale oczywiście w rozsądny sposób (pozostaw pozostałe parametry funkcji, jak podano).

Na koniec kliknij Applyprzycisk, aby zastosować zmiany.

Zobaczysz coś takiego:

zgodnie z oczekiwaniami.

EDYTOWAĆ

Zgodnie ze szczególnym żądaniem zgłoszonym przez PO w komentarzu:

„Czy byłoby możliwe utworzenie tego wzoru tylko między początkiem i końcem każdej polilinii zamiast między każdym wierzchołkiem?”

Lekko zredagowałem kod. Następująca funkcja powinna zwrócić oczekiwany wynik:

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