Создать дугу линий из линии и значения

9

Я пытаюсь воссоздать сюжет Origin-Destination так:

введите описание изображения здесь

Мне удалось объединить данные в таблицу MSOA to LAD и нарисовать карту, подобную этой, для одного из исходных MSOA.

введите описание изображения здесь

Который, как только вы учитываете (теперь смешные) расстояния, которые люди в Пиковом районе добираются до работы, близок.

Но мне очень нравится эффект, которого автор достиг, "расклеивая" линии. Очевидно, что с потоками 522 и 371 я не могу выбрать одну линию на каждого пассажира, но было бы неплохо создать пропорциональную дугу линий, чтобы показать количество людей, совершающих поездку.

Я думал, что смогу использовать Geometry Generator, но без конструкции цикла я не могу добиться прогресса.

Ян Тертон
источник
Этот инструмент ESRI может быть интересен вам или, по крайней мере, является трамплином для идей кода по созданию «клина» линий.
Хорнбидд
И когда линия символизирует, скажем, 50 (100, 200) пассажиров на линию? С некоторым кодом на Python (или с генератором геометрии, в котором я не уверен) вы можете вращать линии (х / 50) с разным количеством.
Стефан

Ответы:

5

Отличный вызов!

Этот ответ в основном использует генератор геометрии и был написан в QGIS 3.2. Сбой QGIS (без сохранения!) Сразу после того, как я сначала построил линии, и я почти сдался, но недавно использованный список выражений спас день - еще один бонус к использованию генератора геометрии

Я начал с двух наборов точек, одного источника и трех пунктов назначения. Направления помечены со счетом:

Начальные точки

Затем я сгенерировал линии, соединяющие исходную точку со всеми пунктами назначения, используя виртуальный слой, используя следующий код:

SELECT d.Count_MF, Makeline( s.geometry, d.geometry) 'geometry' 
  FROM Source AS s JOIN Destinations AS d

Связанные точки

Затем я использовал следующее выражение генератора Geometry для стилизации линий:

 intersection(
   geom_from_wkt( 
     'MULTILINESTRING ((' ||  $x_at( 0)  || ' ' || $y_at( 0)  || ', ' || 
     array_to_string(
       array_remove_at( string_to_array( regexp_replace(
             geom_to_wkt(nodes_to_points( tapered_buffer(  $geometry ,0, "Count_MF" * 200, floor("Count_MF" / 10)),true)),
             '[\\(\\)]','')),0)
     , ') , ('  ||  $x_at( 0)  || ' ' || $y_at( 0)  || ', ' )
    || '))')
    ,buffer( point_n(  $geometry ,1), $length))

Это занимает каждую строку и применяет следующие шаги:

  1. Создает конусный буфер, идущий от нулевой ширины в источнике до ширины, масштабируемой по счетчику назначения в конце назначения. Плотность точек буфера также масштабируется с помощью атрибута количества мест назначения.
  2. Вершины буферного полигона преобразуются в точки (это, возможно, излишне), а затем экспортируются в WKT, а скобки удаляются с помощью регулярного выражения перед преобразованием в массив.
  3. Затем массив расширяется обратно до строки WKT для многоканальной строки, вставляя координаты исходной точки и соответствующее форматирование - это создает отдельную линию для каждой из извлеченных вершин, соединенных с исходной точкой.
  4. WKT преобразуется обратно в геометрический объект и, наконец, пересекается с буфером исходной точки, чтобы обрезать их обратно в круг, на котором находится точка назначения (см. Выходные данные a, tapered_bufferчтобы понять, почему это необходимо)

Фанаты

При написании шагов я понимаю, что преобразование в массив и из него не нужно, и все манипуляции с WKT можно выполнять с помощью регулярных выражений. Это выражение ниже, и если tapered_arrayфункция может быть заменена другой, то это также может быть использовано в QGIS 2.18.

intersection(
   geom_from_wkt(
    'MULTILINESTRING ((' ||  $x_at( 0)  || ' ' || $y_at( 0)  || ', ' ||
  replace(
    regexp_replace(
      regexp_replace(
        geom_to_wkt(tapered_buffer(  $geometry ,0, "Count_MF" * 200, floor("Count_MF" / 10))),
      '^[^,]*,',''),
    ',[^,]*$',''),
  ',',') , ('  ||  $x_at( 0)  || ' ' || $y_at( 0)  || ', ')
  || '))')
,buffer( point_n(  $geometry ,1), $length))
Энди Харфут
источник
6

Ваш вопрос сделал меня любопытным.

Это решение работает только для QGIS 2.x в консоли Python

Как упомянуто в моем комментарии здесь моя идея создать дугу линий с Python.

У меня есть два точечных слоя:

я. Тот, кто держит капитал (id, капитал)

II. Один держит города (идентификатор, город, пассажиры)

Количество пассажиров «разделено на банкноты», и это будут линии, которые образуют дугу. Таким образом, пассажиры 371 представляют собой комбинацию 3х100, 1х50, 2х10 и 1х1 и всего 7 банкнот. После этого линии оформляются стилем, основанным на правилах.

Вот код:

from qgis.gui import *
from qgis.utils import *
from qgis.core import *
from PyQt4 import QtGui, uic
from PyQt4.QtGui import *
from PyQt4.QtCore import *

for lyr in QgsMapLayerRegistry.instance().mapLayers().values():
    if lyr.name() == "capital":
        capital_layer = lyr

for lyr in QgsMapLayerRegistry.instance().mapLayers().values():
    if lyr.name() == "town":
        town_layer = lyr

    # creating the memory layer
d_lyr = QgsVectorLayer('LineString', 'distance', 'memory')
QgsMapLayerRegistry.instance().addMapLayer(d_lyr)
prov = d_lyr.dataProvider()
prov.addAttributes( [ QgsField("id", QVariant.Int), QgsField("banknote",QVariant.Int)])

    # function to create the banknotes
def banknoteOutput(number):
    number_list = []
    number_list.append(number)
    banknote_count = []
    temp_list = []
    banknote_list = []
    for n in number_list:
        total_sum = 0
        total = int(n/100)
        total_sum = total_sum + total
        if total > 0:
            banknote_count.append([total, 100])
        n = n-(total*100)
        total = int(n/50)
        total_sum = total_sum + total
        if total > 0:
            banknote_count.append([total, 50])
        n = n-(total*50)
        total = int(n/10)
        total_sum = total_sum + total
        if total > 0:
            banknote_count.append([total, 10])
        n = n-(total*10)
        total = int(n/5)
        total_sum = total_sum + total
        if total > 0:
            banknote_count.append([total, 5])
        n = n-(total*5)
        total = int(n/1)
        total_sum = total_sum + total
        if total > 0:
            banknote_count.append([total, 1])
        for i in banknote_count:
            temp_list.append(i*i[0])
        banknote_list = [item for sublist in temp_list for item in sublist][1::2]
        return banknote_list

        # creating lines with the amount of banknotes
for capital in capital_layer.getFeatures():
    for town in town_layer.getFeatures():
        commuter_splitting = banknoteOutput(town['commuters'])
        for i,banknote in enumerate(commuter_splitting):
            angle = 2
            distance = QgsDistanceArea()
            distance.measureLine(capital.geometry().asPoint(), town.geometry().asPoint())
            vect = QgsFeature()
            vect.setGeometry(QgsGeometry.fromPolyline([capital.geometry().asPoint(), town.geometry().asPoint()]))
            vect.geometry().rotate(0+(i*angle), capital.geometry().asPoint())
            vect.setAttributes([int(town["id"]), int(banknote)])
            prov.addFeatures([vect])

d_lyr.updateExtents()
d_lyr.triggerRepaint()
d_lyr.updateFields()

Результат может выглядеть так:

введите описание изображения здесь

ОБНОВЛЕНИЕ: различие мужчина / женщина

Результаты в 4 слоя памяти.

from qgis.gui import *
from qgis.utils import *
from qgis.core import *
from PyQt4 import QtGui, uic
from PyQt4.QtGui import *
from PyQt4.QtCore import *

for lyr in QgsMapLayerRegistry.instance().mapLayers().values():
    if lyr.name() == "capital":
        capital_layer = lyr

for lyr in QgsMapLayerRegistry.instance().mapLayers().values():
    if lyr.name() == "town":
        town_layer = lyr

    # function to create the banknotes
def banknoteOutput(number):
    number_list = []
    number_list.append(number)
    banknote_count = []
    temp_list = []
    banknote_list = []
    for n in number_list:
        total_sum = 0
        total = int(n/100)
        total_sum = total_sum + total
        if total > 0:
            banknote_count.append([total, 100])
        n = n-(total*100)
        total = int(n/50)
        total_sum = total_sum + total
        if total > 0:
            banknote_count.append([total, 50])
        n = n-(total*50)
        total = int(n/10)
        total_sum = total_sum + total
        if total > 0:
            banknote_count.append([total, 10])
        n = n-(total*10)
        total = int(n/5)
        total_sum = total_sum + total
        if total > 0:
            banknote_count.append([total, 5])
        n = n-(total*5)
        total = int(n/1)
        total_sum = total_sum + total
        if total > 0:
            banknote_count.append([total, 1])
        for i in banknote_count:
            temp_list.append(i*i[0])
        banknote_list = [item for sublist in temp_list for item in sublist][1::2]
        return banknote_list

    # creating the male memory layer
cmt_male = QgsVectorLayer('LineString', 'Commuters_Male', 'memory')
QgsMapLayerRegistry.instance().addMapLayer(cmt_male)
prov_male = cmt_male.dataProvider()
prov_male.addAttributes( [ QgsField("id", QVariant.Int), QgsField("banknote",QVariant.Int)])

    # creating the male polygon memory layer
cmt_male_polygon = QgsVectorLayer('Polygon', 'Commuters_Male_Poly', 'memory')
QgsMapLayerRegistry.instance().addMapLayer(cmt_male_polygon)
prov_cmt_male_polygon = cmt_male_polygon.dataProvider()
prov_cmt_male_polygon.addAttributes( [ QgsField("id", QVariant.Int), QgsField("banknote",QVariant.Int)])

    # creating lines with the amount of banknotes
for capital in capital_layer.getFeatures():
    for town in town_layer.getFeatures():
        commuter_splitting = banknoteOutput(town['cmt_male'])
        points = []
        for i,banknote in enumerate(reversed(commuter_splitting)):
            angle = 2
            distance = QgsDistanceArea()
            distance.measureLine(capital.geometry().asPoint(), town.geometry().asPoint())
            vect = QgsFeature()
            vect.setGeometry(QgsGeometry.fromPolyline([capital.geometry().asPoint(), town.geometry().asPoint()]))
            vect.geometry().rotate(0+(i*angle), capital.geometry().asPoint())
            vect.setAttributes([int(town["id"]), int(banknote)])
            points.append(vect.geometry().asPolyline()[1])
            prov_male.addFeatures([vect])
        polygon = QgsFeature()
        points.insert(0,capital.geometry().asPoint())
        points.insert(len(points),capital.geometry().asPoint())
        polygon.setGeometry(QgsGeometry.fromPolygon([points]))
        polygon.setAttributes([1, 2])
        prov_cmt_male_polygon.addFeatures([polygon])

cmt_male.updateExtents()
cmt_male.triggerRepaint()
cmt_male.updateFields()
cmt_male_polygon.updateExtents()
cmt_male_polygon.triggerRepaint()
cmt_male_polygon.updateFields()

    # creating the female memory layer
cmt_female = QgsVectorLayer('LineString', 'Commuters_Female', 'memory')
QgsMapLayerRegistry.instance().addMapLayer(cmt_female)
prov_female = cmt_female.dataProvider()
prov_female.addAttributes( [ QgsField("id", QVariant.Int), QgsField("banknote",QVariant.Int)])

    # creating the female polygon memory layer
cmt_female_polygon = QgsVectorLayer('Polygon', 'Commuters_Female_Poly', 'memory')
QgsMapLayerRegistry.instance().addMapLayer(cmt_female_polygon)
prov_cmt_female_polygon = cmt_female_polygon.dataProvider()
prov_cmt_female_polygon.addAttributes( [ QgsField("id", QVariant.Int), QgsField("banknote",QVariant.Int)])

    # creating lines with the amount of banknotes
for capital in capital_layer.getFeatures():
    for town in town_layer.getFeatures():
        commuter_splitting = banknoteOutput(town['cmt_female'])
        points = []
        for i,banknote in enumerate(commuter_splitting):
            angle = 2
            distance = QgsDistanceArea()
            distance.measureLine(capital.geometry().asPoint(), town.geometry().asPoint())
            vect = QgsFeature()
            vect.setGeometry(QgsGeometry.fromPolyline([capital.geometry().asPoint(), town.geometry().asPoint()]))
            vect.geometry().rotate(-angle-(i*angle), capital.geometry().asPoint())
            vect.setAttributes([int(town["id"]), int(banknote)])
            points.append(vect.geometry().asPolyline()[1])
            prov_female.addFeatures([vect])
        polygon = QgsFeature()
        points.insert(0,capital.geometry().asPoint())
        points.insert(len(points),capital.geometry().asPoint())
        polygon.setGeometry(QgsGeometry.fromPolygon([points]))
        polygon.setAttributes([1, 2])
        prov_cmt_female_polygon.addFeatures([polygon])

cmt_female.updateExtents()
cmt_female.triggerRepaint()
cmt_female.updateFields()
cmt_female_polygon.updateExtents()
cmt_female_polygon.triggerRepaint()
cmt_female_polygon.updateFields()

Результат может выглядеть так:введите описание изображения здесь

Одна вещь, которая не идеальна с картографической точки зрения:

На первый взгляд размер дуги линии может вызывать раздражение, так как большая дуга может представлять больше пассажиров. Дуга может быть больше с меньшим количеством пассажиров (289 пассажиров / 11 банкнот), чем с большим количеством пассажиров (311 пассажиров / 5 банкнот).

Стефан
источник