Estoy tratando de procesar algunos datos ráster usando ogr / gdal y parece que no puedo obtener la plena utilización de todos los núcleos de mi máquina. Cuando solo ejecuto el proceso en un solo núcleo, obtengo el 100% de utilización de ese núcleo. Cuando trato de dividirme en multinúcleo (en el...
En nuestras encuestas aéreas de vida silvestre, generalmente generamos docenas de archivos GPX, y generalmente convertimos la mayoría de ellos en archivos shape para mostrar en una aplicación SIG. Los archivos GPX son una excelente manera de mantener todos (o casi todos) sus datos GPS en un solo...
Necesito ayuda ... Tengo problemas con QGIS. Tengo algunos rásteres que parecen tener valores nulos, aunque los reclasifiqué para obtener un ráster booleano :( ¿Cómo puedo cambiar los valores NULL a
Soy un tipo de vectores en su mayoría, pero para un proyecto actual necesito trabajar con un montón de viejos mapas de papel escaneados (¡mapas de daños de bombas ww2 para Londres, si estás interesado!) Escaneamos y georreferenciamos los mapas, y ahora queremos crear una capa compuesta en mosaico...
Estoy usando ArcGIS 10.1 y quiero crear un nuevo ráster basado en dos rásteres preexistentes. El RasterToNumPyArray tiene un buen ejemplo que quiero adaptarse. import arcpy import numpy myArray = arcpy.RasterToNumPyArray('C:/data/inRaster') myArraySum = myArray.sum(1) myArraySum.shape =...
Tengo ortofotos de un área que necesito para servir como mapa de fondo para mis vectores. Lo obtuve como 2500 archivos (71,5MB cada uno) en formato TIFF sin formato con el archivo de palabra TFW correspondiente: 180GB de datos. El sistema de coordenadas es local y coincide con mis vectores (no hay...
¿Hay alguna manera de establecer todos los píxeles de valor 0 en un ráster en nodata usando GDAL? Me gustaría usar fillnodata para rellenar algunos agujeros, pero los agujeros están establecidos en 0, por lo que primero debo
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
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 2 años . Antecedentes...
Quiero proyectar este mapa en proyección robinson: library(ggmap) world <- map_data("world") ggplot() + geom_path(data = world, aes(long, lat, group = group)) Y me gustaría cambiar la proyección a "Robinson" (siguiendo los consejos de la respuesta a mi pregunta anterior: ¿Qué proyección...
Estoy analizando los caminos probables tomados por los viajeros en su camino al trabajo, en base a los datos del flujo del Censo, las distribuciones de densidad de trabajo y la red de transporte. En este momento estoy usando pgRouting y el complemento pgRouting para QGIS (descrito aquí ) para mi...
Dado un objeto espacial en R, ¿cómo recorto todos sus elementos para que queden dentro de un cuadro delimitador? Hay dos cosas que me gustaría hacer (idealmente, sabría cómo hacer ambas cosas, pero cualquiera de las dos es una solución aceptable para mi problema actual: restringir un archivo de...
Publiqué en el foro de arcgis pero nunca obtuve respuesta. ¿Alguien puede decirme cómo resolver esto? Suponga que tiene un registro de seguimiento de vehículos del 1 de octubre al 31 de diciembre con información sobre la velocidad del viento. Todos los datos se guardan en la geodatabase...
Tengo una pregunta con respecto a trazar un mapa coroplético. ¿Qué se recomienda con respecto a trazar un mapa con saltos personalizados? Y haz que se vea un poco atractivo. En mi script , he usado el comando plot y spplot para visualizar algunos datos, pero no estoy completamente satisfecho con...
¿Existe una biblioteca JavaScript de código abierto como Leaflet u OpenLayers que se puede usar en una aplicación web o móvil y carga mapas basados en vectores en lugar de imágenes? Solo necesito la información del camino y algunos PDI. La única biblioteca que encontré fue polymaps , pero aún...
He usado la gIntersectfunción del rgeospaquete para determinar qué polígonos se cruzan en dos SpatialPolygonsDataFrames (spdfs). El resultado es una matriz lógica con nombres de fila que comprende las ID de polígono de sfdf # 1 y colname las ID de polígono de spdf # 2. Me gustaría usar esta matriz...
Tenía un complemento de trabajo en QGIS 1.8 que leía datos de una base de datos MSAccess y los agregaba a una serie de capas de memoria. Hay algo de procesamiento involucrado en el ínterin, así que no creo que usar QGIS para leer directamente desde la base de datos sea una opción. Me gustaría...
Mi pregunta es sobre la extracción de valores de rásteres en la ubicación de los puntos. Con la función de extracción, esto es muy fácil, y la función me da un marco de datos con los valores de todas las variables en los puntos. Quiero tener en ese marco de datos las coordenadas de cada punto....
Cerrada . Esta pregunta está basada en la opinión . Actualmente no está aceptando respuestas. ¿Quieres mejorar esta pregunta? Actualice la pregunta para que pueda ser respondida con hechos y citas editando esta publicación . Cerrado hace 4 años . He sido...
Tengo estos dos polígonos: library(sp); library(rgeos); library(maptools) coords1 <- matrix(c(-1.841960, -1.823464, -1.838623, -1.841960, 55.663696, 55.659178, 55.650841, 55.663696), ncol=2) coords2 <- matrix(c(-1.822606, -1.816790, -1.832712, -1.822606, 55.657887, 55.646806, 55.650679,...