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 =...
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...
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....
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...
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...
Quiero crear dos polígonos. Uno es de las extensiones rectangulares de mi trama. Sé que puedo obtener la extensión de la trama usando r <- raster("band5.tif") e <- extent(r) plot(e)me da esto: ¿Cómo puedo crear UN polígono del límite del ráster como se muestra a continuación? Así es...
Quería convertir Raster en Point Vector en QGIS. ¿Es posible? QGIS tiene la opción de convertirlo en polígono, pero no encontré ninguna herramienta para convertirlo en punto. ¿Alguien puede
Tengo dos archivos LiDAR (.las), uno es original, digamos con X puntos. Y el otro es una copia del primer archivo .las pero con puntos Y, donde Y es menor que X. Ahora, quiero comparar los Modelos Digitales de Elevación (DEM) de estos dos archivos .las y visualizar cuán diferentes son. Quiero...
Tengo un Arc / Info Binary Grid --- específicamente, un ráster de acumulación de flujo ArcGIS --- y me gustaría identificar todas las celdas que tienen un valor específico (o en un rango de valores). En última instancia, me gustaría un archivo de formas de puntos que representen estas...
Estoy trabajando en un proyecto que incluirá tanto "datos ráster" (imágenes aéreas y DEM) como "datos cuadriculados", que entiendo que son otros tipos de mediciones numéricas tomadas en puntos de cuadrícula, y que podrían estar completamente no relacionados con la imagen (p. ej., densidad de...
Tengo un dataset ráster en formato ASCII. Usando Python, necesito determinar los valores miny maxdentro del conjunto de datos. Me han dicho que la información del encabezado es clave, que contiene elementos como el número de filas / columnas, el tamaño de la celda, etc. ¿No puede simplemente...
Tengo una trama de tipos de hábitat para un área específica en Escocia. Necesito crear escenarios de hábitat futuros con cambios en el hábitat para evaluar la viabilidad de la población de una especie de ave. Por ejemplo, en el futuro podría haber un 10% más de silvicultura en el área. Me gustaría...
Intentando llegar al fondo de por qué, cuando leo en un ráster de NDVI, el espacio @ data @ valores no contiene los valores reales hasta que los configure manualmente. Por ejemplo: NDVI <- raster("./filename.tif", crs="+proj=longlat +datum=WGS84") NDVI@data@values ## returns:...
Solo me pregunto, ¿siempre es mejor tener más de 3 GCP? ¿Deberían los PCG distribuirse de manera uniforme? ¿Dónde puedo encontrar algunos recursos generales sobre la teoría de esto? Además de probar un error, que simplemente me enseña qué funciona en ese momento y no qué debería funcionar en...
He creado una capa de polígono a partir de rásteres, y solo tengo una única identificación de clasificación (y acres calculados) para cada polígono. Me gustaría "disolver" todas las características conmovedoras en polígonos individuales más grandes. Las 3 áreas encerradas en rojo (abajo) contienen...
¿Hay alguna manera de realizar la misma tarea que la utilidad gdalbuildvrt usando enlaces GDAL Python? Hasta ahora no he encontrado ninguna manera de hacer esto que no sea crear un vrt de un único conjunto de datos y editar manualmente el xml. Me gustaría crear un vrt a partir de múltiples rásteres...
He leído el código fuente de varios filtros ráster QGis-1.7.4 que calculan la pendiente, el aspecto y la curvatura. Hay una fórmula en el filtro que calcula la curvatura total que me desconcierta. El archivo fuente está en la versión actual de QGis, con la siguiente ruta: qgis-1.7.4 / src /...
Tengo que convertirlo en un UTM para que el búfer funcione. wets<-readOGR(dsn=".",layer="shapefile") r.raster <- raster() extent(r.raster) <- extent(wets) res(r.raster) <- 100 wets.r <- rasterize(wet,r.raster) plot(wets.r) wetsbuf<-buffer(wets.r,width=500) Durante la creación...
Estoy tratando de convertir datos ráster en datos poligonales. Anteriormente estaba usando shapefile.py para crear datos de puntos en lat / lon individuales en el ráster. Pero eso creará una gran cantidad de datos y su consulta será un gran problema. Por lo tanto, sería mejor si pudiera convertir...
Estoy tratando de grabar un shapefile en un ráster usando RasterizeLayer de GDAL. Precreo un ráster de área de interés desde un archivo de formas diferente, dado un tamaño de píxel específico. Este AOI sirve entonces como base para todas las siguientes rasterizaciones (mismo número de columnas y...