mirror of
				https://github.com/qgis/QGIS.git
				synced 2025-10-31 00:06:02 -04:00 
			
		
		
		
	Prevent changes to files that weren't changed between releases. This eases review of the changes between releases significantly.
		
			
				
	
	
		
			161 lines
		
	
	
		
			5.7 KiB
		
	
	
	
		
			Python
		
	
	
	
	
	
			
		
		
	
	
			161 lines
		
	
	
		
			5.7 KiB
		
	
	
	
		
			Python
		
	
	
	
	
	
| # -*- coding: utf-8 -*-
 | |
| 
 | |
| """
 | |
| ***************************************************************************
 | |
|     ExtentFromLayer.py
 | |
|     ---------------------
 | |
|     Date                 : August 2012
 | |
|     Copyright            : (C) 2012 by Victor Olaya
 | |
|     Email                : volayaf at gmail 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.                                   *
 | |
| *                                                                         *
 | |
| ***************************************************************************
 | |
| """
 | |
| 
 | |
| __author__ = 'Victor Olaya'
 | |
| __date__ = 'August 2012'
 | |
| __copyright__ = '(C) 2012, Victor Olaya'
 | |
| 
 | |
| import os
 | |
| from math import floor, ceil
 | |
| 
 | |
| from qgis.PyQt.QtGui import QIcon
 | |
| from qgis.PyQt.QtCore import QVariant
 | |
| 
 | |
| from qgis.core import (QgsApplication,
 | |
|                        QgsField,
 | |
|                        QgsFeatureSink,
 | |
|                        QgsGeometry,
 | |
|                        QgsFeature,
 | |
|                        QgsWkbTypes,
 | |
|                        QgsProcessing,
 | |
|                        QgsProcessingException,
 | |
|                        QgsProcessingParameterMapLayer,
 | |
|                        QgsProcessingParameterDistance,
 | |
|                        QgsProcessingParameterDefinition,
 | |
|                        QgsProcessingParameterFeatureSink,
 | |
|                        QgsFields)
 | |
| 
 | |
| from processing.algs.qgis.QgisAlgorithm import QgisAlgorithm
 | |
| 
 | |
| pluginPath = os.path.split(os.path.split(os.path.dirname(__file__))[0])[0]
 | |
| 
 | |
| 
 | |
| class ExtentFromLayer(QgisAlgorithm):
 | |
| 
 | |
|     INPUT = 'INPUT'
 | |
|     BY_FEATURE = 'BY_FEATURE'
 | |
|     ROUND_TO = 'ROUND_TO'
 | |
| 
 | |
|     OUTPUT = 'OUTPUT'
 | |
| 
 | |
|     def icon(self):
 | |
|         return QgsApplication.getThemeIcon("/algorithms/mAlgorithmExtractLayerExtent.svg")
 | |
| 
 | |
|     def svgIconPath(self):
 | |
|         return QgsApplication.iconPath("/algorithms/mAlgorithmExtractLayerExtent.svg")
 | |
| 
 | |
|     def tags(self):
 | |
|         return self.tr('polygon,vector,raster,extent,envelope,bounds,bounding,boundary,layer,round,rounded').split(',')
 | |
| 
 | |
|     def group(self):
 | |
|         return self.tr('Layer tools')
 | |
| 
 | |
|     def groupId(self):
 | |
|         return 'layertools'
 | |
| 
 | |
|     def __init__(self):
 | |
|         super().__init__()
 | |
| 
 | |
|     def initAlgorithm(self, config=None):
 | |
|         self.addParameter(QgsProcessingParameterMapLayer(self.INPUT, self.tr('Input layer')))
 | |
| 
 | |
|         round_to_parameter = QgsProcessingParameterDistance(self.ROUND_TO,
 | |
|                                                             self.tr('Round values to'),
 | |
|                                                             parentParameterName=self.INPUT,
 | |
|                                                             minValue=0,
 | |
|                                                             defaultValue=0)
 | |
|         round_to_parameter.setFlags(round_to_parameter.flags() | QgsProcessingParameterDefinition.FlagAdvanced)
 | |
|         self.addParameter(round_to_parameter)
 | |
| 
 | |
|         self.addParameter(QgsProcessingParameterFeatureSink(self.OUTPUT, self.tr('Extent'), type=QgsProcessing.TypeVectorPolygon))
 | |
| 
 | |
|     def name(self):
 | |
|         return 'polygonfromlayerextent'
 | |
| 
 | |
|     def displayName(self):
 | |
|         return self.tr('Extract layer extent')
 | |
| 
 | |
|     def processAlgorithm(self, parameters, context, feedback):
 | |
|         layer = self.parameterAsLayer(parameters, self.INPUT, context)
 | |
| 
 | |
|         round_to = self.parameterAsDouble(parameters, self.ROUND_TO, context)
 | |
| 
 | |
|         fields = QgsFields()
 | |
|         fields.append(QgsField('MINX', QVariant.Double))
 | |
|         fields.append(QgsField('MINY', QVariant.Double))
 | |
|         fields.append(QgsField('MAXX', QVariant.Double))
 | |
|         fields.append(QgsField('MAXY', QVariant.Double))
 | |
|         fields.append(QgsField('CNTX', QVariant.Double))
 | |
|         fields.append(QgsField('CNTY', QVariant.Double))
 | |
|         fields.append(QgsField('AREA', QVariant.Double))
 | |
|         fields.append(QgsField('PERIM', QVariant.Double))
 | |
|         fields.append(QgsField('HEIGHT', QVariant.Double))
 | |
|         fields.append(QgsField('WIDTH', QVariant.Double))
 | |
| 
 | |
|         (sink, dest_id) = self.parameterAsSink(parameters, self.OUTPUT, context,
 | |
|                                                fields, QgsWkbTypes.Polygon, layer.crs())
 | |
|         if sink is None:
 | |
|             raise QgsProcessingException(self.invalidSinkError(parameters, self.OUTPUT))
 | |
| 
 | |
|         try:
 | |
|             # may not be possible
 | |
|             layer.updateExtents()
 | |
|         except:
 | |
|             pass
 | |
| 
 | |
|         rect = layer.extent()
 | |
| 
 | |
|         if round_to > 0:
 | |
|             rect.setXMinimum(floor(rect.xMinimum() / round_to) * round_to)
 | |
|             rect.setYMinimum(floor(rect.yMinimum() / round_to) * round_to)
 | |
|             rect.setXMaximum(ceil(rect.xMaximum() / round_to) * round_to)
 | |
|             rect.setYMaximum(ceil(rect.yMaximum() / round_to) * round_to)
 | |
| 
 | |
|         geometry = QgsGeometry.fromRect(rect)
 | |
| 
 | |
|         minx = rect.xMinimum()
 | |
|         miny = rect.yMinimum()
 | |
|         maxx = rect.xMaximum()
 | |
|         maxy = rect.yMaximum()
 | |
|         height = rect.height()
 | |
|         width = rect.width()
 | |
|         cntx = minx + width / 2.0
 | |
|         cnty = miny + height / 2.0
 | |
|         area = width * height
 | |
|         perim = 2 * width + 2 * height
 | |
| 
 | |
|         feat = QgsFeature()
 | |
|         feat.setGeometry(geometry)
 | |
|         attrs = [
 | |
|             minx,
 | |
|             miny,
 | |
|             maxx,
 | |
|             maxy,
 | |
|             cntx,
 | |
|             cnty,
 | |
|             area,
 | |
|             perim,
 | |
|             height,
 | |
|             width,
 | |
|         ]
 | |
|         feat.setAttributes(attrs)
 | |
|         sink.addFeature(feat, QgsFeatureSink.FastInsert)
 | |
| 
 | |
|         return {self.OUTPUT: dest_id}
 |