Skip to content
Snippets Groups Projects

Enh documentation supers2

Merged Cresson Remi requested to merge enh_documentation_supers2 into main
9 files
+ 121
62
Compare changes
  • Side-by-side
  • Inline
Files
9
@@ -8,12 +8,14 @@ api = Client.open(
modifier=sign_inplace
)
bbox = [4, 42.99, 5, 44.05]
def mosa(year):
"""Return a pyotb application that perform a mosaic."""
res = api.search(bbox=bbox, datetime=[f'{year}-01-01', f'{year}-12-25'])
res = api.search(
bbox=[4, 42.99, 5, 44.05],
datetime=[f"{year}-01-01", f"{year}-12-25"],
collections=["spot-6-7-drs"]
)
urls = [f"/vsicurl/{r.assets['src_xs'].href}" for r in res.items()]
return pyotb.Mosaic({"il": urls})
@@ -23,7 +25,7 @@ def ndvi(xs):
return pyotb.BandMath({"il": [xs], "exp": "(im1b4-im1b1)/(im1b4+im1b1)"})
ndvi_22 = ndvi(mosa('2022'))
ndvi_21 = ndvi(mosa('2021'))
ndvi_22 = ndvi(mosa("2022"))
ndvi_21 = ndvi(mosa("2021"))
delta_ndvi = ndvi_22 - pyotb.Superimpose({"inr": ndvi_22, "inm": ndvi_21})
delta_ndvi.write("ndvi_loss.tif?&box=5000:5000:4096:4096")
Loading