Vincent Sarago

IMO we should do this at the same time as https://github.com/mapbox/rasterio/issues/1774 so we can add types 🤷‍♂️

@glostis indeed, this is interesting you can also totally remove any transform and you'll get the same result ```python path = "s3://sentinel-s1-l1c/GRD/2020/3/17/IW/DV/S1B_IW_GRDH_1SDV_20200317T173959_20200317T174024_020735_0274FA_1C70/measurement/iw-vv.tiff" session = AWSSession(requester_pays=True) with rasterio.Env(session): with rasterio.open(path) as...

other interesting notes, the transform resulting from gdalwarp is not equal to the one we get with `from_gcps` `gdalwarp -ot float32 -multi -t_srs EPSG:4326 -tps iw-vv.tiff iw-vv.vrt` ```python with rasterio.open("iw-vv.vrt")...

thanks @glostis I've updated my previous comment. so with the file I'm using the transform created from gcps has `78.20` value while with gdalwarp it's `75.61` when checking the gcps,...

I think this is specific to your GDAL install (/home/a001673/mambaforge/envs/satpy/share/gdal). I have no issue using the file on either my laptop (Mac OS with home-brew gdal) or other env (docker)....

@howff the COG driver doesn't allow direct `write` but only `COPY` ```python import rasterio from rasterio.shutil import copy from rasterio.io import MemoryFile import numpy as np from pathlib import Path...

@howff oops, yep it was missing `driver="GTiff"` ```python import rasterio from rasterio.shutil import copy from rasterio.io import MemoryFile import numpy as np from pathlib import Path pth = Path('/tmp/') /...

seems similar to https://github.com/rasterio/rasterio/issues/2339

Go ahead @sgillies 👌 I'm making baby step on this subject. About the mask, I think I've got a solution explained in https://github.com/OSGeo/gdal/pull/742#issuecomment-463346403, the main problem for me is that...

Ref of past discussion : https://github.com/mapbox/rasterio/issues/1191#issuecomment-398171026 we decided in the past to materialized the mask to an alpha band because it was how gdal worked ;-)