mirror of
https://github.com/qgis/QGIS.git
synced 2025-03-03 00:02:25 -05:00
121 lines
4.4 KiB
Python
121 lines
4.4 KiB
Python
# -*- coding: utf-8 -*-
|
|
|
|
"""
|
|
***************************************************************************
|
|
Gridify.py
|
|
---------------------
|
|
Date : May 2010
|
|
Copyright : (C) 2010 by Michael Minn
|
|
Email : pyqgis at michaelminn dot com
|
|
***************************************************************************
|
|
* *
|
|
* This program is free software; you can redistribute it and/or modify *
|
|
* it under the terms of the GNU General Public License as published by *
|
|
* the Free Software Foundation; either version 2 of the License, or *
|
|
* (at your option) any later version. *
|
|
* *
|
|
***************************************************************************
|
|
"""
|
|
from builtins import str
|
|
|
|
__author__ = 'Michael Minn'
|
|
__date__ = 'May 2010'
|
|
__copyright__ = '(C) 2010, Michael Minn'
|
|
|
|
# This will get replaced with a git SHA1 when you do a git archive
|
|
|
|
__revision__ = '$Format:%H$'
|
|
|
|
from qgis.core import (QgsFeature,
|
|
QgsGeometry,
|
|
QgsPoint,
|
|
QgsWkbTypes,
|
|
QgsApplication)
|
|
from processing.core.GeoAlgorithm import GeoAlgorithm
|
|
from processing.core.GeoAlgorithmExecutionException import GeoAlgorithmExecutionException
|
|
from processing.core.parameters import ParameterVector
|
|
from processing.core.parameters import ParameterTableField
|
|
from processing.core.outputs import OutputVector
|
|
|
|
from processing.tools import dataobjects, vector
|
|
|
|
|
|
class HubLines(GeoAlgorithm):
|
|
HUBS = 'HUBS'
|
|
HUB_FIELD = 'HUB_FIELD'
|
|
SPOKES = 'SPOKES'
|
|
SPOKE_FIELD = 'SPOKE_FIELD'
|
|
OUTPUT = 'OUTPUT'
|
|
|
|
def icon(self):
|
|
return QgsApplication.getThemeIcon("/providerQgis.svg")
|
|
|
|
def svgIconPath(self):
|
|
return QgsApplication.iconPath("providerQgis.svg")
|
|
|
|
def group(self):
|
|
return self.tr('Vector analysis tools')
|
|
|
|
def name(self):
|
|
return 'Hub lines'
|
|
|
|
def displayName(self):
|
|
return self.tr('Hub lines')
|
|
|
|
def defineCharacteristics(self):
|
|
self.addParameter(ParameterVector(self.HUBS,
|
|
self.tr('Hub layer')))
|
|
self.addParameter(ParameterTableField(self.HUB_FIELD,
|
|
self.tr('Hub ID field'), self.HUBS))
|
|
self.addParameter(ParameterVector(self.SPOKES,
|
|
self.tr('Spoke layer')))
|
|
self.addParameter(ParameterTableField(self.SPOKE_FIELD,
|
|
self.tr('Spoke ID field'), self.SPOKES))
|
|
|
|
self.addOutput(OutputVector(self.OUTPUT, self.tr('Hub lines'), datatype=[dataobjects.TYPE_VECTOR_LINE]))
|
|
|
|
def processAlgorithm(self, feedback):
|
|
layerHub = dataobjects.getObjectFromUri(
|
|
self.getParameterValue(self.HUBS))
|
|
layerSpoke = dataobjects.getObjectFromUri(
|
|
self.getParameterValue(self.SPOKES))
|
|
|
|
fieldHub = self.getParameterValue(self.HUB_FIELD)
|
|
fieldSpoke = self.getParameterValue(self.SPOKE_FIELD)
|
|
|
|
if layerHub.source() == layerSpoke.source():
|
|
raise GeoAlgorithmExecutionException(
|
|
self.tr('Same layer given for both hubs and spokes'))
|
|
|
|
writer = self.getOutputFromName(self.OUTPUT).getVectorWriter(
|
|
layerSpoke.fields(), QgsWkbTypes.LineString, layerSpoke.crs())
|
|
|
|
spokes = vector.features(layerSpoke)
|
|
hubs = vector.features(layerHub)
|
|
total = 100.0 / len(spokes)
|
|
|
|
for current, spokepoint in enumerate(spokes):
|
|
p = spokepoint.geometry().boundingBox().center()
|
|
spokeX = p.x()
|
|
spokeY = p.y()
|
|
spokeId = str(spokepoint[fieldSpoke])
|
|
|
|
for hubpoint in hubs:
|
|
hubId = str(hubpoint[fieldHub])
|
|
if hubId == spokeId:
|
|
p = hubpoint.geometry().boundingBox().center()
|
|
hubX = p.x()
|
|
hubY = p.y()
|
|
|
|
f = QgsFeature()
|
|
f.setAttributes(spokepoint.attributes())
|
|
f.setGeometry(QgsGeometry.fromPolyline(
|
|
[QgsPoint(spokeX, spokeY), QgsPoint(hubX, hubY)]))
|
|
writer.addFeature(f)
|
|
|
|
break
|
|
|
|
feedback.setProgress(int(current * total))
|
|
|
|
del writer
|