switch to morecantile3

[vincentsarago]

ref #417

[vincentsarago]

🤔 recent changes in rasterio break tests 😭

[vincentsarago]

💥 I can't get CI to run because we need rio-cogeo, which also pin morecantile to >=2.1,<2.2 but in order to make CI run in rio-cogeo (if I update it to morecantile 3) I'll need rio-tiler 3 🤦

🔁 dependencies 🕳️

[vincentsarago]

Tests are failing because for an unknown reason the bounds returned by rasterio==1.2.8 are not the same as the one for rasterio==1.2.6 and it now considers the tile/point to be outside the bounds.

from rasterio.crs import CRS
from rasterio.warp import transform_bounds
in_crs = CRS.from_wkt('PROJCS["WGS_1984_Lambert_Azimuthal_Equal_Area",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]],PROJECTION["Lambert_Azimuthal_Equal_Area"],PARAMETER["latitude_of_center",45],PARAMETER["longitude_of_center",-100],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["metre",1],AXIS["Easting",EAST],AXIS["Northing",NORTH]]')
out_crs = CRS.from_epsg(4326)
bbox = (-4410000.000000002, 1453684.1641738773, -1506678.966595509, 4309999.999999999)
print(transform_bounds(in_crs, out_crs, *bbox, densify_pts=21))

# rasterio 1.2.6
>> (-179.91600712033699, 42.116151939674644, 179.04299285803737, 76.46962486562524)

# rasterio 1.2.8
>> (-160.73152298505175, 32.22950673496283, -108.2439562098202, 71.16432886506686)


#weirder: results from pyproj
from pyproj import CRS, Transformer
in_crs = CRS.from_wkt('PROJCS["WGS_1984_Lambert_Azimuthal_Equal_Area",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]],PROJECTION["Lambert_Azimuthal_Equal_Area"],PARAMETER["latitude_of_center",45],PARAMETER["longitude_of_center",-100],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["metre",1],AXIS["Easting",EAST],AXIS["Northing",NORTH]]')
out_crs = CRS.from_epsg(4326)
tr = Transformer.from_crs(in_crs, out_crs, always_xy=True)
bbox = (-4410000.000000002, 1453684.1641738773, -1506678.966595509, 4309999.999999999)
print(tr.transform_bounds(*bbox, densify_pts=21))

### Pyproj 3.2 and 3.1
>> (166.11052567893609, 42.116151939674644, -124.50200312011823, 76.46962486562524)