Hi, Answer in English although I understand French (after all an english speaking list here)
You can: Instead of vsistr = "/vsicurl?pc_url_signing=yes&pc_collection=sentinel-2-l2a&url=" red_href = vsistr + planetary_computer.sign(item).assets["B04"].href cogtif = gdal.Open( red_href ) use red_href = f"/vsicurl/{planetary_computer.sign(item).assets["B04"].href }&pc_url_signing=yes&pc_collection=sentinel-2-l2a" cogtif = gdal.Open( red_href ) Regards Thomas Gratier Le mer. 28 sept. 2022 à 10:12, Loïc Lozac'h <loic.loz...@inrae.fr> a écrit : > Bonjour, > > Est-il possible d'utiliser la signature de la librairie python > Planetary_Computer pour lire des COG depuis le catalogue STAC Microsoft ? > > J'ai essayé plusieurs méthode en utilisant /vsicurl/ mais les extended > filenames associé à la signature ne sont pas reconnu. > > Est-il possible de lire ces données à la manière de RasterIO ? > > J'utilise gdal version 3.5.2 > > code python en erreur : > vsistr = "/vsicurl?pc_url_signing=yes&pc_collection=sentinel-2-l2a&url=" > red_href = vsistr + planetary_computer.sign(item).assets["B04"].href > cogtif = gdal.Open( red_href ) > > response : > > Warning 6: Unsupported option: se > Warning 6: Unsupported option: sp > Warning 6: Unsupported option: sv > Warning 6: Unsupported option: sr > Warning 6: Unsupported option: skoid > Warning 6: Unsupported option: sktid > Warning 6: Unsupported option: skt > Warning 6: Unsupported option: ske > Warning 6: Unsupported option: sks > Warning 6: Unsupported option: skv > Warning 6: Unsupported option: sig > ERROR 11: HTTP response code: 403 > > code python ok : > with fiona.open(shp, "r") as shapefile: > shapes = [feature["geometry"] for feature in shapefile] > with rasterio.open(planetary_computer.sign(item).assets["B04"].href ) as > ds: > out_image, out_transform = rasterio.mask.mask(ds, shapes, crop=True) > out_meta = ds.meta > out_meta.update({"driver": "GTiff", > "height": out_image.shape[1], > "width": out_image.shape[2], > "transform": out_transform}) > with rasterio.open(out, "w", **out_meta) as dest: > dest.write(out_image) > > Merci pour vos réponses > > Bien cordialement, > > Loïc > _______________________________________________ > gdal-dev mailing list > gdal-dev@lists.osgeo.org > https://lists.osgeo.org/mailman/listinfo/gdal-dev >
_______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/gdal-dev