Vorrei ottenere l'intersezione di più poligoni. Usando il shapelypacchetto di Python , posso trovare l'intersezione di due poligoni usando la intersectionfunzione. Esiste una funzione efficiente simile per ottenere l'intersezione di più poligoni? Ecco uno snippet di codice per capire cosa intendo: from shapely.geometry import Point coord1 = ( 0,0 ) …
Ci sono state un sacco di domande su questo, ma nessuno di loro ha risolto il mio problema. Ho una macchina con: Windows 7 x64 Python 3.4.3 x64 GDAL installata: gdalinfo --version --> GDAL 1.11.4, released 2016/01/25; gdal è stato installato da una ruota GDAL-1.11.4-cp34-none-win_amd64.whlda qui Il GDAL è installato …
Esiste un pacchetto Python che fornisce un'implementazione di un algoritmo Straight Skeleton? Sono consapevole che il progetto CGAL open source (C ++) contiene un'implementazione ma sembra che cgal-bindings non includa questo pacchetto CGAL . In ogni caso, preferirei un'implementazione Python pura che potrei modificare / estendere in base alle mie …
Vorrei aggiornare un attributo di una funzione. Tuttavia, non so aggiornarlo utilizzando l'oggetto funzione. Devo utilizzare il fornitore di dati per aggiornarlo. layers = QgsMapLayerRegistry.instance().mapLayersByName('my_line') layer = layers[0] dp = layer.dataProvider() it = dp.getFeatures() for i in range(0, dp.featureCount()): feat = it.next() attrs = { 2 : 30 } layer.dataProvider().changeAttributeValues({ …
Sto riscontrando problemi nell'uso della funzione Polygonize in Python. L'esempio del libro di cucina per questo può essere trovato qui . La parte rilevante del mio codice è: sourceRaster = gdal.Open('myraster.tif') band = sourceRaster.GetRasterBand(1) bandArray = band.ReadAsArray() outShapefile = "polygonized" driver = ogr.GetDriverByName("ESRI Shapefile") if os.path.exists(outShapefile+".shp"): driver.DeleteDataSource(outShapefile+".shp") outDatasource = driver.CreateDataSource(outShapefile+ …
Sto leggendo uno shapefile usando Fiona: with fiona.open('data/boroughs/boroughs_n.shp') as source: mpolys = [p for p in source] candidate = polys[0]['geometry'] questo mi dà un dict contenente un elenco di elenchi di coordinate e digitare "MultiPolygon" Tuttavia, non sono sicuro di come utilizzare i dati delle coordinate per creare un MultiPolygon …
Esiste un modo per eseguire la stessa attività dell'utilità gdalbuildvrt utilizzando i collegamenti GDAL Python? Finora non ho trovato alcun modo per farlo se non creare un vrt di un singolo set di dati e modificare manualmente l'xml. Vorrei creare un vrt da più raster (essenzialmente eseguendo un mosaico). È …
sto cercando un metodo per calcolare la coordinata X, Y di un punto interno (noto anche come labelPoint in arcmap) di un poligono (irregolare). Sto provando a trovare una soluzione senza l'uso del modulo ArcGIS la differenza tra un "centroide" e un labelPoint è che il centroide può trovarsi al …
[EDIT: la soluzione a questo era semplicemente usare OGR per leggere gli shapefile. Vedi l'esempio di geographika.] In uno shapefile ESRI, non vi è alcuna distinzione tra poligoni e multipoligoni. Inoltre, non vi è alcuna distinzione esplicita tra fori interni e anelli esterni (oltre alla "mano" di un determinato poligono). …
Quali driver ha il pacchetto python fiona? Quando controllo il manuale , [...] and the possible formats are enumerated in the fiona.drivers list. tuttavia, viene visualizzato quando digito Python from fiona import drivers print drivers > <function drivers at 0x108763050> Come posso guardare "in" quello?
Sto lavorando a uno strumento che converte i layer raster in array per l'elaborazione con NumPy, e idealmente mi piacerebbe poter lavorare con raster che vengono impacchettati in un .gdb senza esportarli tutti (specialmente se ciò richiede ArcGIS o ArcPy coinvolgenti ). È possibile con il driver OpenFileGDB? Da quello …
usando Gdal / OGR c'è un modo per convertire un livello di linea shp, kml o PostGIS in una rete di nodi e segmenti in modo che possa essere utilizzato in pacchetti come networkx?
Sto cercando di aprire un file di forma con il notebook Jupyter, ma non riesco ad aprirlo. Lo shapefile si trova nella stessa cartella della mia cartella di lavoro e Python non può identificarlo anche se utilizzo il percorso completo. import shapefile as shp import matplotlib.pyplot as plt sf = …
Ho un progetto QGIS che contiene un compositore. Ovviamente se vado nella GUI posso esportarlo come pdf ecc. Per il mio script make, anche se voglio farlo tramite uno script Python. Dalla lettura online ecco il mio tentativo: #!/usr/bin/env python import sys from qgis.gui import QgsMapCanvas, QgsLayerTreeMapCanvasBridge from PyQt4.QtCore import …
We use cookies and other tracking technologies to improve your browsing experience on our website,
to show you personalized content and targeted ads, to analyze our website traffic,
and to understand where our visitors are coming from.
By continuing, you consent to our use of cookies and other tracking technologies and
affirm you're at least 16 years old or have consent from a parent or guardian.