From e71199db931894a32c855554c1410b4305c18926 Mon Sep 17 00:00:00 2001 From: Guilhem Moulin Date: Wed, 27 Sep 2023 18:39:26 +0200 Subject: =?UTF-8?q?Without=20--geometry,=20set=20extent=20on=20the=20first?= =?UTF-8?q?=20--point=20=C2=B1=20--margin.?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- gis-observation-map | 11 ++++++++++- 1 file changed, 10 insertions(+), 1 deletion(-) (limited to 'gis-observation-map') diff --git a/gis-observation-map b/gis-observation-map index 19dabd3..e2740c7 100755 --- a/gis-observation-map +++ b/gis-observation-map @@ -36,7 +36,9 @@ from qgis.core import ( QgsLayerDefinition, QgsLayerTreeModel, QgsMapLayer, + QgsPointXY, QgsProject, + QgsRectangle, QgsUnitTypes, QgsVectorLayer, QgsWkbTypes @@ -152,7 +154,6 @@ if projectInstance is not None: #layerTreeView.setCurrentLayer(layer) canvas = QgsMapCanvas() canvas.setExtent( layer.extent() ) - canvas.refresh() def gpkg_intersects(path, geometries): ds = gpkg_drv.Open(path, update=0) @@ -266,6 +267,14 @@ for path in args.geometry: geometries.append(geometry) feature = layer.GetNextFeature() +if len(args.geometry) == 0: + r = geometries[0].GetGeometryRef(0) + pt0X, pt0Y = r.GetPoint_2D(0) + pt1X, pt1Y = r.GetPoint_2D(2) + canvas = QgsMapCanvas() + extent = QgsRectangle(QgsPointXY(pt0X, pt0Y), QgsPointXY(pt1X, pt1Y)) + canvas.setExtent(extent) + def geometricFilterPolygon(polygon): polygon2 = [] if polygon.GetGeometryType() != ogr.wkbPolygon: -- cgit v1.2.3