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...
Tengo dos capas en ArcGIS: una es un punto vectorial y la otra es una capa ráster. Quiero agregar el valor de la cuadrícula de datos ráster al punto vectorial. El problema es que el punto vectorial está lejos de cualquier cuadrícula ráster existente con valor, por lo que el comando "Muestra" no...
Estoy interesado en una forma de identificar rásteres (p. Ej., Formato ESRI GRID) o archivos de forma corruptos que impidan su uso en una herramienta de análisis. Tengo un gran número en el árbol de carpetas estructuradas y me gustaría poder escribir algo que recorra el árbol de archivos y pueda...
Me gustaría importar todos los rásteres de mosaico múltiples (cada mosaico está compuesto por 7 capas) de la carpeta indicada en mi R. Luego, acceda a ellos como rásteres de múltiples capas individuales. Así que lo hice: # read all mosaics named "mos....img" in R
Tengo un DEM de una región montañosa. Me gustaría extraer las áreas del DEM por encima de las cuales se encuentra el 60% del área total del DEM. Entonces, esencialmente me gustaría extraer las áreas que representan las partes superiores de las montañas. He oído hablar de esto conocido como la...
Refiriéndome a mi pregunta anterior , ¿cómo puedo "geocodificar" / convertir, usando R, una imagen ráster generada a un formato que pueda mostrarse como una capa en un mapa? En resumen, mi objetivo es recrear los resultados del software comercial SpatialKey ( captura de pantalla ). Estoy usando...
Estoy tratando de crear un ráster usando las bibliotecas Python gdal y llegué al punto en el que se emiten los datos, pero los datos de salida se voltean en el eje x del punto de origen . Sé que debo pasar por alto algo, pero no puedo entender a dónde me estoy equivocando. ¿Algunas ideas? Al crear...
Tengo dos imágenes ráster, booleanas clasificadas (1 o 2). Se superponen en áreas determinadas. Necesito un ráster que sea solo las secciones superpuestas. Uso de Arcinfo 10. Answer solo puede usar métodos de trama, sin conversión a polígonos, puntos,
Me gustaría pasar de un ráster binario de bosque / no bosque a una resolución de 30 m a un ráster a una resolución de 240 m con el valor de cada celda como la proporción de esa celda arbolada, es decir, los valores de datos para las celdas de 240 m irían de 0 (donde todas las celdas en el ráster de...
Mi experiencia : programador experimentado, historiador aficionado, principiante completo en SIG (así que estoy usando QGIS). Sin embargo, aprendí a usar la herramienta de georreferenciación QGIS. La entrada: supongamos que tengo una imagen como esta Wikipedia: Archivo: VoyagesOfRabbanBarSauma.jpg...
Tengo un ráster GeoTIFF en ArcGIS Desktop, que se puede descargar desde los siguientes dos enlaces: http://www.rtwilson.com/downloads/MODIS_MOD04_16-06-2006_Subset.tfw http://www.rtwilson.com/downloads/MODIS_MOD04_16-06-2006_Subset.tif Si carga estos en ArcGIS Desktop, verá grandes áreas...
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...
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...
Me gustaría agregar atributos a un ráster almacenado en una base de datos postgis. Utilizo el siguiente comando para insertar el ráster en la base de datos: raster2pgsql -s 4236 -I -C -M /home/user/mysrc/outlabeled.tif public.test | psql -d template_postgis El ráster outlabeled.tifes, de hecho,...
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...
Estoy tratando de recortar un ECW grande (detalles a continuación) pero el archivo ráster es demasiado grande para ser procesado por completo. Algunos detalles de la ECW a continuación Controlador: Ondas comprimidas ECW / ERDAS (SDK 5.0) tamaño de archivo: 50 gb El tamaño es 450000, 565081...
Tengo una imagen rasterizada de elevaciones. He agregado un archivo de formas vectoriales que contiene divisiones políticas locales. Quiero mantener los colores de elevación, pero también quiero mantener las líneas negras que representan las divisiones políticas. He reducido la transparencia del...
Estoy intentando mejorar un proceso de vector / pitón extremadamente engorroso para un modelo de peligro natural. Por el momento tenemos una secuencia de comandos larga que genera líneas de distancia / rumbo desde un punto determinado para determinar: El tipo de polígono que se cruza (p. ej.,...
Compilé e instalé el maestro QGIS actual como de costumbre e importé algunos datos ráster en mi base de datos postgres 9.3 con la extensión postgis 2.0.2. Cada vez que intento importar los datos ráster en qgis usando el administrador de la base de datos, recibo este mensaje: GDAL provider: Cannot...
Tengo dos capas .ecw superpuestas de la misma área, aunque una está recortada en el medio, con un borde blanco. Cuando configuro el valor de fondo en 255,255,255, y lo defino como hueco, solo vacía ese valor RGB específico, cuando en realidad, hay otros valores "blancos" en mi capa,...