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'}
Estoy usando rasterio 1.0a10 - ¿alguna idea de lo que está pasando aquí?