Sistemi informativi geografici

Domande e risposte per cartografi, geografi e professionisti GIS

3
mappatura di oltre un milione di punti usando QGIS
Sto cercando di mappare un CSV con oltre un milione di punti lat / lon negli Stati Uniti senza metadati. Non ho mai provato a mappare così tanti dati e voglio sapere quali passi dovrei prendere. So che ci sono troppi dati per QGIS, quindi ho provato a suddividerli in …
9 qgis  turf 




3
Esecuzione di una query spaziale in un ciclo in PyQGIS
Quello che sto cercando di fare: ciclo attraverso un punto di shapefile e selezionare ogni punto che cade in un poligono. Il codice seguente è ispirato a un esempio di query spaziale che ho trovato in un libro: mitte_path = r"D:\PythonTesting\SelectByLocation\mitte.shp" punkte_path = r"D:\PythonTesting\SelectByLocation\punkte.shp" polygon = QgsVectorLayer(mitte_path, 'Mitte', 'ogr') points …



1
Come faccio a creare uno SpatialLinesDataFrame da un frame di dati?
Ho creato un data.frame in R. Voglio ottenere SpatialLinesDataFrame come risultato, quindi capisco che devo convertire data.frame in linee, le linee in SpatialLines e SpatialLines in SpatialLinesDataFrame. l'ho fatto filedata.frame=data.frame(matrix(file),ncol=14, byrow=T)) file=Lines(filedata.frame) Ottengo il data.frame che voglio ma non le linee. Sono sicuro che è di base, ma perché non …
9 r  line 


1
Qual è la differenza tra i token "SHAPE @ XY" e "SHAPE @ TRUECENTROID" di accesso ai dati?
Cursore di ricerca dell'accesso ai dati . Qual è la differenza tra i token "SHAPE@XY"e "SHAPE@TRUECENTROID"? Restituirono le stesse coordinate su un poligono di prova: >>> cursor = arcpy.da.SearchCursor("test", "SHAPE@XY") >>> for row in cursor: ... print row[0] ... (559389.3838043335, 4239093.201390337) >>> cursor = arcpy.da.SearchCursor("test", "SHAPE@TRUECENTROID") >>> for row in …



8
Massimizzare l'uso della CPU
La mia sceneggiatura interseca le linee con i poligoni. È un processo lungo poiché ci sono più di 3000 linee e più di 500000 poligoni. Ho eseguito da PyScripter: # Import import arcpy import time # Set envvironment arcpy.env.workspace = r"E:\DensityMaps\DensityMapsTest1.gdb" arcpy.env.overwriteOutput = True # Set timer from datetime import …



Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.