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:...
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...
Este es un proceso que lleva solo unos segundos en el software SIG. Mi intento de hacerlo en R utiliza una gran cantidad de memoria y luego falla. ¿Hay algo mal en mi código o es algo que R no puede hacer? He leído que R puede funcionar dentro de Grass, ¿puedo usar una función Grass desde dentro de...
¿Qué tan importante es corregir Landsat para la atmósfera, la bruma y convertir el DN en reflectancia? ¿Qué impacto tiene en los resultados, por ejemplo, para muchos de los índices de
Tengo un archivo ráster donde todas las celdas tienen valores entre -3 y 7, ningún valor de datos es -9999. ¿Cómo puedo calcular el número de celdas con un valor específico, por ejemplo, 6? ¿Puedo usar la
Estoy creando algunos rásteres de punto flotante en GDAL con el controlador GeoTIFF. Cuando cargo las imágenes resultantes en QGIS o Arc, la simbología predeterminada pone el mínimo en -3.40282e38 y el máximo en 3.40282e38, por lo que el ráster se ve exactamente gris. ¿Hay alguna manera de poder...
Tengo un conjunto de rásteres (.vrt) con datos diarios de humedad del suelo. Quiero sumar los píxeles de todos los rásteres para tener una medida por mes. Sin embargo, el mundo entero no está cubierto todos los días, lo que da como resultado valores de nodatos en lugares donde existen medidas para...
He importado una imagen tif ráster en PostgreSQL con este comando: raster2pgsql -s 32643 -I -M filepath.tif -F -t 100x100 public.databassename > filepath.sql E importó el archivo SQL de salida dentro de la base de datos PostgreSQL que ejecuta esto: psql -U postgres -d databasename -f...
Estoy escribiendo un método de Python que agrega dos rásteres y genera una única salida de ráster. Por razones que escapan a mi control, la extensión de los rásteres de entrada es diferente, pero se superponen. ¿Es posible operar exclusivamente en los píxeles de ráster de entrada que se...
Para un estudio sobre la expansión de las redes de carreteras en la selva tropical, estoy tratando de extraer caminos de las imágenes de Landsat. Ya tenemos compuestos nítidos y libres de nubes en los que las carreteras son claramente visibles a simple vista, pero extraerlas en las características...
Estoy tratando de hacer una coincidencia de histograma con Python para mejorar el proceso de mosaico de múltiples rásteres superpuestos. Estoy basando mi código en el que se encuentra en: http://www.idlcoyote.com/ip_tips/histomatch.html Hasta la fecha, he logrado recortar el área superpuesta de...
He descargado varios mosaicos de trama SRTM DEM de https://dwtkns.com/srtm30m/ , y he creado versiones de sombreado a partir de los mosaicos. Sin embargo, después de una inspección más cercana, los mosaicos de sombreado ahora tienen líneas de borde que no sé cómo ocultar o eliminar (como en la...
He estado tratando de convertir archivos de velocidad del viento en rásteres sin éxito. Fue discutido en un foro de Arcgis que seguí. Esto implicó reemplazar los punto y coma con espacios y agregar un encabezado. Esperaba incluir uno de los archivos .asc originales junto con una versión .csv que...
Tengo que hacer un análisis de criterios múltiples para responder la pregunta: "cuál es el mejor lote para desarrollar". Algunos de los criterios son: distancia de la parada de autobús más cercana (capa de puntos con paradas de autobús) distancia de la tienda más cercana (capa de puntos con...
Tengo un mosaico de un DEM derivado de LiDar de 1 m. Necesito generar un subconjunto de datos como un DEM de 10 m. Actualmente estoy usando la herramienta de agregado en ARCGIS 10 para producir un valor medio para cada nuevo píxel de 10 m. ¿Algún consejo sobre si esta es la mejor técnica para tal...
Esto es 4 imágenes (NE, SE, NW, SW) secciones de imágenes. Tienen cierta superposición entre ellos. Trabajando con 1m de ortoimágenes de 3 bandas. Algunas veces se vuelan diferentes partes en diferentes momentos / fuentes, de ahí la diferencia en épocas del año. Las imágenes están en formato...
En mi trabajo diario, constantemente me piden que calcule áreas de datasets ráster globales en proyección geográfica a una resolución de 30 segundos de arco. Estos conjuntos de datos son normalmente el resultado de una operación Combinar (un ejemplo típico es una clase de vegetación combinada con...
... en lugar de mantener los rásteres como GEOTIFF almacenados localmente? Tengo un gran proyecto que requerirá poner aproximadamente 50,000, 10mb GEOTIFF como WMS a un mapa web basado en OpenLayers. Solo se mostrará 1 WMS en un momento dado, por lo que el rendimiento no es un problema. Teniendo...