Sto servendo tessere vettoriali usando TileStache , ho tutto impostato come voglio. I miei dati sono archiviati in Postgres e sto usando il provider VecTiles per servire i riquadri GeoJSON . Voglio mettere in cache tutte le mie tessere per rendere le tessere più veloci. Sto usando tilestache-seed.py per eseguire …
Ho un geopandas GeoDataFrame contenente centinaia di formosePolygon e MultiPolygongeometrie. I poligoni si sovrappongono in molti punti. Vorrei creare una nuova geometria che contasse il numero di quanti si sovrappongono. Qualcosa come questo: Qualcuno ha qualche idea su come affrontare questo? Non riesco nemmeno a vedere un modo per entrare. …
Ho difficoltà a comprendere l'uso degli indici spaziali con RTree. Esempio: ho 300 punti bufferizzati e ho bisogno di conoscere l'area di intersezione di ciascun buffer con un file di forma poligonale. Il file di forma poligonale ha> 20.000 poligoni. Mi è stato suggerito di utilizzare indici spaziali per accelerare …
Chiuso . Questa domanda deve essere più focalizzata . Al momento non accetta risposte. Vuoi migliorare questa domanda? Aggiorna la domanda in modo che si concentri su un problema solo modificando questo post . Chiuso 4 anni fa . Sto cercando una buona risorsa (MOOC, libro, ecc.) Per insegnare a …
Vorrei chiamare la funzione del plugin di interpolazione QGIS (metodo TIN) (Raster-> Interpolate) dalla console di Python. Non riesco a trovare la funzione corrispondente nell'API QGIS o nell'elenco degli algoritmi di elaborazione. Ho trovato l'algoritmo SAGA Triangulation, che funziona bene ma è 5-10 volte più lento e la velocità è …
Sono nuovo in questo. Sto "costruendo" un plugin Python per qgis e voglio importare moduli da GRASS e SAGA ma non capisco come farlo. Il plugin che voglio importare è l'interpolazione di Kriging da SAGA o v.surf.idw da GRASS o il plugin di interpolazione da Raster qgis. Voglio un metodo …
Sto cercando di cambiare un plugin da QGIS 1.8 a 2.0. Nella versione 1.8 creo una tabella con i nomi dei campi in questo modo: layer = vlayer provider = layer.dataProvider() feat = QGSFeatrue() allAttrs = provider.attributeIndexes() provider.select(allAttrs) fieldNames = [] while provider.nextFeature(feat) fields = layer.pendingFields() for (k, field) in …
Vedi questo link per maggiori dettagli. Il problema: Voglio passare in rassegna un raster continuo (uno che non ha una tabella di attributi), cella per cella e ottenere il valore della cella. Voglio prendere quei valori ed eseguire i condizionali su di essi, emulando i passaggi di algebra della mappa …
Sono riuscito a aggiungere QgsVectorLayerin QGISattraverso un plugin python, per esempio: vl = self.iface.addVectorLayer(uri.uri(), layerName, self.dbConn.getProviderName()) Ma come posso rimuoverlo?
Creo regolarmente i miei raster GeoTIFF con GDAL in Python, ad esempio: from osgeo import gdal from numpy import random data = random.uniform(0, 10, (300, 200)) driver = gdal.GetDriverByName('GTiff') ds = driver.Create('MyRaster.tif', 200, 300) band = ds.GetRasterBand(1) band.WriteArray(data) ds = band = None # save, close tuttavia, quando il risultato …
Non ero abbastanza soddisfatto del calcolo della lunghezza dei linestring in WGS84 in miglia . Mi ha fatto domandare se esiste un modo più conveniente, Pythonic per calcolare la lunghezza di una stringa lineare WKT secondo un dato SRID. Ho in mente qualcosa come: srid="WGS84" line="LINESTRING(3.0 4.0, 3.1 4.1)" print …
Lo script che sto scrivendo riposiziona due frame di dati e ne imposta l'estensione. In questo modo, ridisegna l'intera Visualizzazione attiva 4 volte, rallentando sostanzialmente lo script. Premendo F9 o facendo clic sul pulsante "Pausa disegno" prima di eseguire lo script, il disegno viene disabilitato e lo script viene eseguito …
Sto usando Shapely in Python e mi viene dato un MultiLineStringcon un mucchio di Linestringoggetti. Posso garantire che tutti gli LineStringoggetti sono linee semplici con solo 2 vertici e che fanno tutti parte di un'unica linea (senza rami). Voglio "collegare i punti" e creare un singolo LineString. Devo scrivere un …
Nella documentazione per gdal_calc è indicato il calcolatore raster della riga di comando con sintassi numpy . Più avanti ci sono alcuni esempi in cui in uno di essi: gdal_calc.py -A input.tif --outfile = result.tif --calc = "A * (A> 0)" --NoDataValue = 0 - significa impostare i valori di …
Sto usando il codice qui sotto per trovare un paese (e talvolta lo stato) per milioni di punti GPS. Il codice attualmente richiede circa un secondo per punto, che è incredibilmente lento. Lo shapefile è di 6 MB. Ho letto che Geopandas usa le risorse per i join spaziali, rendendole …
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.