1 votos

Rasterio.mask devuelve el error 'las formas están fuera de los límites del raster' - pero definitivamente están dentro de los límites del raster

Estoy usando rasterio.mask.mask para enmascarar un raster con un shapefile.

import rasterio
import fiona

image_object = rasterio.open('/home/cate/image_name.img')

aoi_shapefile = fiona.open('/home/cate/AOI.shp', "r")
aoi = [polygon["geometry"] for polygon in aoi_shapefile]

image_aoi, image_aoi_transform = rasterio.mask.mask(image_object, aoi)

y me sale UserWarning: shapes are outside bounds of raster. Are they in different coordinate reference systems?

Por lo que puedo decir, la forma está definitivamente dentro de los límites de la trama..

image_object.bounds
BoundingBox(left=117.4095061135422, bottom=33.73209332058161, right=120.49755473423146, top=35.64128279392083)

aoi_shapefile.bounds
(119.82677788035676, 34.076577156414714, 120.37042708436398, 34.575927115244426)

image_object.crs
CRS({'init': 'epsg:4326'})

aoi_shapefile.crs
{'init': 'epsg:4326'}

Esta es la vista en QGIS: enter image description here

Estoy usando rasterio 1.0a10 - ¿alguna idea de lo que está pasando aquí?

2voto

Antonio Haley Puntos 2588

Cate, este es un error de rasterio que se ha corregido en la versión 1.0a12 (en PyPI ahora).

pip install --pre -U rasterio

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X