Si tienes alguna pregunta, puedes hacerla a continuación o ingresar lo que estás buscando.
De forma predeterminada, rasterio.plot.show_hist(…) establece el título del histograma como “Histograma”, la etiqueta del eje x como “DN” (para Número Digital) y la etiqueta del eje y como “Frecuencia”. Existe un parámetro para el ‘título’, pero nada para la etiqueta del eje x y la etiqueta del eje y. ¿Cómo . . . Read more
Me gustaría convertir una imagen (.tiff) en puntos Shapely. Hay 45 millones de píxeles, necesito una forma de lograr esto sin un bucle (actualmente lleva más de 15 horas). Por ejemplo, tengo un archivo .tiff que, cuando se abre, es una matriz de 5000×9000. Los valores son valores de píxeles . . . Read more
Estoy intentando descargar el paquete rasterio para analizar conjuntos de datos raster en Python dentro de las disciplinas de análisis geoespacial. He descargado las dependencias de rasterio, incluyendo gdal, pyproj y fiona a través de los archivos “wheel” de esta fuente de binarios de análisis geoespacial: https://www.lfd.uci.edu/~gohlke/pythonlibs/ Instalé estos paquetes . . . Read more
Actualmente estoy usando este código para convertir un archivo de raster a un geodataframe: import rasterio from rasterio.features import shapes mask = None with rasterio.open(#INSERTAR ARCHIVO TIF AQUÍ) as src: image = src.read(1) # primer banda, aún no estoy seguro de cómo hacerlo con varias bandas results = ( {‘properties’: . . . Read more
Tengo un script que solía funcionar para calcular estadísticas zonales (mediana), pero ahora obtengo el AttributeError: el objeto ‘DatasetReader’ no tiene el atributo ‘afín’. Aquí está mi código: con rasterio.open(f'{project_data}/ras.tif’) como raster: matriz = raster.read(1) afín = raster.afín estadística = zonal_stats(f'{project_data}/Lila.shp’, matriz, affine=afín, stats=[‘mediana’], geojson_out=True) resultado = {“type”: “FeatureCollection”, “características”: . . . Read more