Me gustaría realizar una superposición espacial para identificar el polígono dentro del cual cae un conjunto de puntos. Sin embargo, primero quiero amortiguar y disolver los polígonos, de modo que los puntos que caen en cualquier lugar dentro de los polígonos fusionados (pero no dentro de ningún...
Tengo un proyecto QGIS con capas, basado en consultas SQL PostGIS que usan la fecha actual. Como los datos en la base de datos cambian, estas consultas devuelven datos diferentes cada día nuevo. ¿Es posible cambiar dinámicamente el nombre de la capa en el panel del árbol de capas, por lo que...
En primer lugar, soy un novato en SIG. Tengo una nube de puntos lidar con GPS y datos de telemetría asociados. Luego, transformo las coordenadas de cada punto en un sistema de coordenadas ECEF y finalmente transformo las coordenadas angulares en un sistema de coordenadas cartesianas, utilizando la...
Quiero cargar archivos .adf en R. Los datos provienen de esta página: http://www.fao.org/geonetwork/srv/en/metadata.show?id=14057 Intenté el siguiente código que encontré después de algunas investigaciones en Internet. El problema es que en la clase RasterLayer obtengo valores negativos que no...
Estoy tratando de proyectar una trama. En R existe la projectRaster()función para esto (debajo de un ejemplo completamente reproducible): # example Raster require(raster) r <- raster(xmn=-110, xmx=-90, ymn=40, ymx=60, ncols=40, nrows=40) r <- setValues(r, 1:ncell(r)) projection(r) # project...
Así que estoy trabajando en un proyecto que ha estado en desarrollo durante más de un año, pero solo recientemente me involucré. Mi tarea consiste en convertir aproximadamente 2.500 archivos KML a otro formato, importar estos archivos en un Geoserver y luego copiar la información de estilo asociada...
Cerrada . Esta pregunta necesita detalles o claridad . Actualmente no está aceptando respuestas. ¿Quieres mejorar esta pregunta? Agregue detalles y aclare el problema editando esta publicación . Cerrado hace 5 meses . Estoy tratando de crear un búfer...
He creado un data.frame en R. Como resultado, quiero obtener SpatialLinesDataFrame, por lo que entiendo que tengo que convertir data.frame en líneas, las líneas en SpatialLines y SpatialLines en SpatialLinesDataFrame. hice esto filedata.frame=data.frame(matrix(file),ncol=14,
Trabajando con QGIS 2.4.0 Chugiak Tengo una pregunta sobre el recorte de imágenes ráster. Tengo una imagen aérea que quiero recortar usando la extensión de un archivo de forma. Hasta ahora todo bien, utilicé la herramienta raster / extracción / clipper para hacer esto. El problema es que el tamaño...
Quiero crear un conjunto de datos vectoriales de puntos de 10000 puntos (o más grandes) dentro de un ráster binario, donde los puntos deben estar restringidos a áreas donde el valor del ráster es 1. Intenté los siguientes pasos. Poligonalizar ráster QGIS: Vector -> Herramientas de...
Estoy usando el método GDAL ReadAsArray para trabajar con datos ráster usando numpy (específicamente reclasificación). Como mis rásteres son grandes, proceso las matrices en bloques, iterando en cada bloque y procesando con un método similar al ejemplo de GeoExamples . Ahora estoy buscando la...
Los barrios reina y torre son dos formas comunes de calcular estadísticas para una celda focal. También se les conoce como vecindarios Moore y (von) Neumann. Creo que también existe una "versión obispo", es decir, desde la esquina inferior izquierda a la superior derecha. Pero, ¿por qué estos dos...
Necesito encontrar un algoritmo o método que pueda detectar latitude longitude puntos atípicos en una trayectoria durante el procesamiento posterior , que luego se puede arreglar (volver a la trayectoria de la trayectoria en función de sus vecinos). Como ejemplo del tipo de puntos atípicos que me...
Intento desarrollar un complemento para agregar capas postgis a qgis desde un complemento como parte del aprendizaje de Python. El código para agregar una capa postgis se toma del libro de cocina y funciona bien si lo ejecuto desde la consola de Python dentro de qgis (si se presiona el botón OK en...
Tengo un gran proyecto que consta de alrededor de 10x archivos de proyecto QGIS. Cada proyecto contiene un compositor de impresión con un atlas configurado para generar hasta 10 mapas. ¿Cuál es la forma más fácil de automatizar la apertura de cada proyecto y crear un PDF utilizando el atlas...
Estoy ejecutando un servidor Windows Server 2012. PostgreSQL 9.3.5, 64 bits POSTGIS 2.1.3 QGIS Server 2.6.1-2 QGIS Desktop 2.8.3 y 2.12 Estoy usando una tableta Microsoft Surface Pro 4 con QGIS Desktop 2.12 para editar algunas capas almacenadas en la base de datos Postgresql. Las capas en la...
Esta es una pregunta de geometría de SQL Server aparentemente simple que pensé que tendría una solución lista para usar, pero no tengo suerte para encontrarla. Mi intención es seleccionar todos los registros dentro de una tabla que tengan polígonos anidados (contenidos) dentro de un polígono más...
Estoy usando la biblioteca R maptools para analizar el archivo de forma en una lista de polígonos. La función readShapeSpatial shp <- readShapeSpatial("<path to my shapefile>") me da un objeto SpatailPolygonDataFrame . En mi ejemplo, mi SpatialPolygonsDataFrame tiene las siguientes...
Tengo una serie de cuadrículas hexagonales de 1 km que cubren varios condados en los Estados Unidos en una base de datos postgreSQL / postGIS. Cada cuadrícula tiene el CRS EPSG: 3857, y la capa de los condados tiene EPSG: 3857. Al ver las cuadrículas con los condados en QGIS, todo parece...
Estoy tratando de crear el choropleth()mapa correcto en R usando el GISToolspaquete. Yo uso choro.legend()para mostrar mi leyenda. Sin embargo, no puedo crear la ubicación correcta de la leyenda, Y / O NO entiendo el papel del fmtparámetro en mi choro.legend()función. Me parece que fmtpuede...