Estoy escribiendo una utilidad simple para recortar lotes de archivos de trama geotiff multibanda en la misma área (más pequeña). Con gdalwarp, puedo recortar fácilmente un archivo usando un archivo de forma de recorte de un solo polígono: gdalwarp -cutline clipper.shp -crop_to_cutline input.tif...
Tengo una pregunta con respecto a la agregación espacial en R. Lo que estoy tratando de hacer es agregar un conjunto de datos de puntos a una cuadrícula. Sin embargo, no estoy seguro de cómo hacerlo, ya que tengo poca experiencia con este tipo de cosas. Esperaba que cualquiera de ustedes pudiera...
¿Cómo puedo calcular sumas, promedios, etc. de puntos ráster (multibanda) por polígono de una capa vectorial? Me dijeron que esto se llama "estadísticas zonales". Lo intenté con QGIS primero. Hay una manera de hacerlo, pero eso es demasiado lento (convertir el ráster en vector, intersecar con la...
Tengo dos rásteres: un modelo de elevación digital (DEM) y una banda Landsat. La resolución de la cuadrícula es de 30 m, pero las cuadrículas no coinciden (ver imagen: el ráster oscuro en la parte superior es el DEM; el ráster claro debajo es la banda Landsat). Quiero realizar alguna...
Tengo una tabla de polígonos que representan las isócronas del tiempo de viaje en días particulares. Para cada punto de origen, hay cinco geometrías de isócronas (almacenadas en filas separadas). Para cada punto de origen, quiero rasterizar las cinco isocronas (un NULL binario o 1), y luego...
Soy muy nuevo en el uso de datos SIG y solo tengo una modesta experiencia con R. He estado leyendo sobre cómo analizar datos espaciales utilizando el libro PDF de espacial-analyst.net, así que no estoy completamente perdido, pero pensé que podría describirlo. mi problema y la gente pueden sugerir...
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 datos ráster en una tabla de postgres que quiero ingresar a python como una matriz numpy. Estoy usando psycopg2 para conectarme a la base de datos. Puedo descargar los datos pero vuelven como una cadena (probablemente un binario serializado). ¿Alguien sabe cómo tomar esta cadena y...
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...
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
Cuando proceso un ráster en QGIS, crea una segunda banda. Parece que la banda no tiene datos útiles. Me gusta usar la herramienta de colores de trama y no aceptará una trama con nada más que en la banda. ¿Cómo elimino la banda
Al leer la literatura SIG, a veces se siente que los autores usan los dos términos con un significado específico en mente. Siempre trato los términos indistintamente. ¿Me equivoco? Lea la primera línea en el enlace . Vea cómo dice que un servicio de imágenes proporciona acceso a ráster (e imagen),...
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...
Tengo un ráster tiff en mi base de datos PostGIS 2.0. Quiero obtener el valor del ráster en un punto lat / lon especificado. Yo tengo: postgresql-9.1 postgis-2.0-svn con soporte raster un ráster importado en PostGIS su SRID es
De vez en cuando, utilizamos datos LiDAR para modelar en ArcGIS. Anteriormente, se nos proporcionaría un .xyz o .ascii de los datos, pero más recientemente, también se nos proporcionan archivos .las además de .xyz y ascii. Estaba leyendo este artículo que dice "Tener lidar en formato LAS puede ser...
Estaba tratando de hacer una imagen de trama de una base de datos de puntos espaciados irregularmente. Los datos se ven como- > head(s100_ras) x y z 1 267573.9 2633781 213.29545 2 262224.4 2633781 69.78261 3 263742.7 2633781 51.21951 4 259328.4 2633781 301.98413 5 264109.8 2633781 141.72414 6...
¿Existe una forma directa de muestrear aleatoriamente un ráster para que la salida del proceso sea un ráster? Estoy usando un ejemplo que encontré en la r-sig-geo lista y también probé la sampleRandomfunción en el rasterpaquete. Ambos enfoques producen una salida que no estoy seguro de cómo...
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 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 . Tengo un gran...
Estoy tratando de usar QGIS v2.0.1 para lograr las mismas tareas que puedo hacer con ArcMap v10.1 y me encuentro con algunos problemas para encontrar herramientas equivalentes. En ArcMap estoy usando la herramienta Focal Statistics (mean) Spatial Analyst y no sé si GRASS o SAGA tiene un equivalente...