Preguntas etiquetadas con python

Estoy sirviendo mosaicos vectoriales usando TileStache , tengo todo configurado como quiero. Mis datos se almacenan en Postgres y estoy usando el proveedor VecTiles para servir mosaicos GeoJSON . Quiero almacenar en caché todos mis mosaicos para que los mosaicos se sirvan más rápido. Estoy usando tilestache-seed.py para sembrar mi …

13 python  postgis  cache  vector-tiles  tilestache 

Tengo un geopandas que GeoDataFrame contiene cientos de formas Polygony MultiPolygongeometrías. Los polígonos se superponen en muchos lugares. Me gustaría hacer una nueva geometría que contenga un recuento de cuántas de ellas se superponen. Algo como esto: ¿Alguien tiene alguna idea sobre cómo abordar esto? Ni siquiera puedo ver una …

13 python  shapely  geopandas  pandas 

Tengo problemas para comprender el uso de índices espaciales con RTree. Ejemplo: tengo 300 puntos almacenados en búfer, y necesito saber el área de intersección de cada búfer con un archivo de forma poligonal. El archivo de forma poligonal tiene> 20,000 polígonos. Se sugirió que usara índices espaciales para acelerar …

13 python  shapely  fiona  spatial-index  rtree 

Cerrada . Esta pregunta necesita estar más centrada . Actualmente no está aceptando respuestas. ¿Quieres mejorar esta pregunta? Actualice la pregunta para que se centre en un problema solo editando esta publicación . Cerrado hace 4 años . Estoy buscando un buen recurso (MOOC, libro, etc.) para enseñar Python para …

13 python  remote-sensing 

Me gustaría llamar a la función del complemento de interpolación QGIS (método TIN) (Raster-> Interpolate) desde la consola de Python. No puedo encontrar la función correspondiente en QGIS API o en la lista de algoritmos de procesamiento. Encontré el algoritmo de triangulación SAGA, que funciona bien pero es 5-10 veces …

13 qgis  python  qgis-plugins  pyqgis 

Soy nuevo en esto. Estoy "construyendo" un complemento de Python para qgis y quiero importar módulos de GRASS y SAGA, pero no entiendo cómo hacerlo. El complemento que quiero importar es la interpolación de Kriging de SAGA o v.surf.idw de GRASS o el complemento de interpolación de qgis de ráster. …

13 qgis  python  grass  qgis-plugins  saga 

Estoy tratando de cambiar un complemento de QGIS 1.8 a 2.0. En la versión 1.8, creo una tabla de los nombres de campo como este: 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 …

13 python  pyqgis  qgis-2 

Vea este enlace para más detalles. El problema: Quiero recorrer un ráster continuo (uno que no tiene tabla de atributos), celda por celda, y obtener el valor de la celda. Quiero tomar esos valores y ejecutar condicionales sobre ellos, emulando los pasos de álgebra de mapas detallados a continuación sin …

13 python  raster  arcobjects  gdal  c# 

Logré agregar QgsVectorLayera QGIStravés de un complemento de Python, por ejemplo: vl = self.iface.addVectorLayer(uri.uri(), layerName, self.dbConn.getProviderName()) Pero, ¿cómo puedo eliminarlo?

13 qgis  python  vector  qgis-plugins 

Regularmente creo mis propios rásteres GeoTIFF con GDAL en Python, por ejemplo: 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 sin embargo, cuando el …

13 python  gdal  geotiff-tiff  statistics 

Estaba bastante insatisfecho con el cálculo de la longitud de las cadenas lineales en WGS84 en millas . Me mantuvo preguntándome si hay una manera más conveniente y pitónica de calcular la longitud de una cadena lineal WKT de acuerdo con un SRID dado. Tengo en mente algo como: srid="WGS84" …

13 python  wkt  wgs84  srid  length 

El script que estoy escribiendo reposiciona dos marcos de datos y establece sus extensiones. Al hacer esto, vuelve a dibujar toda la Vista activa 4 veces, ralentizando sustancialmente el script. Al presionar F9 o hacer clic en el botón 'Pausar dibujo' antes de ejecutar el script, se desactiva el dibujo …

13 arcgis-desktop  python  arcgis-10.0  arcmap 

Estoy usando Shapely en python y me dan un MultiLineStringcon un montón de Linestringobjetos. Puedo garantizar que todos los LineStringobjetos son líneas simples con solo 2 vértices y que todos forman parte de una sola línea (sin ramas). Quiero "conectar los puntos" y crear uno solo LineString. ¿Necesito escribir un …

13 python  shapely 

En la documentación de gdal_calc se indica Calculadora ráster de línea de comando con sintaxis numpy . Más adelante hay algunos ejemplos donde en uno de ellos: gdal_calc.py -A input.tif --outfile = result.tif --calc = "A * (A> 0)" --NoDataValue = 0 - significa establecer valores de cero e inferiores …

13 python  gdal  numpy  gdal-calc 

Estoy usando el siguiente código para encontrar un país (y, a veces, un estado) para millones de puntos GPS. El código actualmente toma aproximadamente un segundo por punto, que es increíblemente lento. El archivo de forma es de 6 MB. Leí que las geopandas usan rtrees para uniones espaciales, lo …

13 python  spatial-join  geopandas  rtree