From 9c0fffb820323388002428c0d4d9bc578bc68ba4 Mon Sep 17 00:00:00 2001 From: rodrigo <> Date: Thu, 25 Feb 2021 10:15:44 +1000 Subject: [PATCH 1/2] - removed that was cleaning all links - added as a fallback to - Removed from position that was unused --- stactools_landsat/stactools/landsat/utils.py | 21 ++++++++++++-------- 1 file changed, 13 insertions(+), 8 deletions(-) diff --git a/stactools_landsat/stactools/landsat/utils.py b/stactools_landsat/stactools/landsat/utils.py index 8e97a299..10a9d834 100644 --- a/stactools_landsat/stactools/landsat/utils.py +++ b/stactools_landsat/stactools/landsat/utils.py @@ -1,3 +1,4 @@ + import datetime import dateutil @@ -76,7 +77,6 @@ def transform_stac_to_stac(item: Item, item.ext.enable("eo") # Add and update links - item.links = [] if self_link: item.links.append(Link(rel="self", target=self_link)) if source_link: @@ -100,18 +100,23 @@ def transform_stac_to_stac(item: Item, if enable_proj: try: # If we can load the blue band, use it to add proj information - blue_asset = item.assets["SR_B2.TIF"] - blue = rasterio.open(blue_asset.href) - shape = [blue.height, blue.width] - transform = blue.transform - crs = blue.crs.to_epsg() + if item.assets.get("SR_B2.TIF"): + asset = item.assets["SR_B2.TIF"] + # SR_B10 is a fallback for SR_B2 + elif item.assets.get("SR_B10.TIF"): + asset = item.assets["SR_B10.TIF"] + else: + raise ValueError('Asset SR_B2.TIF or SR_B2.TIF required') + + opened_asset = rasterio.open(asset.href) + shape = [opened_asset.height, opened_asset.width] + transform = opened_asset.transform + crs = opened_asset.crs.to_epsg() # Now we have the info, we can make the fields item.ext.enable("projection") item.ext.projection.epsg = crs - new_assets = {} - for name, asset in item.assets.items(): if asset.media_type == "image/vnd.stac.geotiff; cloud-optimized=true": item.ext.projection.set_transform(transform, asset=asset) From 7cf5c9dddbc8ce3ccbe66ed8a663832d6c350b2f Mon Sep 17 00:00:00 2001 From: rodrigo <> Date: Thu, 25 Feb 2021 11:11:13 +1000 Subject: [PATCH 2/2] - Formatted --- stactools_landsat/stactools/landsat/utils.py | 1 - 1 file changed, 1 deletion(-) diff --git a/stactools_landsat/stactools/landsat/utils.py b/stactools_landsat/stactools/landsat/utils.py index 10a9d834..80001810 100644 --- a/stactools_landsat/stactools/landsat/utils.py +++ b/stactools_landsat/stactools/landsat/utils.py @@ -1,4 +1,3 @@ - import datetime import dateutil