diff --git a/notebooks/Sentinel-2/Sentinel-2_L1C_data_access.ipynb b/notebooks/Sentinel-2/Sentinel-2_L1C_data_access.ipynb index 17f3a15..4ab918c 100644 --- a/notebooks/Sentinel-2/Sentinel-2_L1C_data_access.ipynb +++ b/notebooks/Sentinel-2/Sentinel-2_L1C_data_access.ipynb @@ -84,12 +84,13 @@ "output_type": "stream", "text": [ "Discovering Sentinel-2 products from STAC catalog...\n", - "Found 18 Sentinel-2 L1C products\n" + "Found 2 Sentinel-2 L1C products\n" ] } ], "source": [ "import pystac_client\n", + "from datetime import date, timedelta\n", "\n", "STAC_URL = \"https://stac.core.eopf.eodc.eu\"\n", "catalog = pystac_client.Client.open(STAC_URL)\n", @@ -97,8 +98,8 @@ "# Define area of interest\n", "bbox = [11.2, 45.5, 11.3, 45.6]\n", "\n", - "# Look for all products starting from the first date until the present\n", - "temporal_extent = [\"2025-10-01\", None]\n", + "# Look for all products starting from a month ago until the present\n", + "temporal_extent = [str(date.today() - timedelta(days=30)), None]\n", "\n", "print(\"Discovering Sentinel-2 products from STAC catalog...\")\n", "s2_items = list(\n", @@ -106,7 +107,7 @@ " collections=[\"sentinel-2-l1c\"],\n", " bbox=bbox,\n", " datetime=temporal_extent,\n", - " query={\"eo:cloud_cover\": {\"lt\": 20}},\n", + " query={\"eo:cloud_cover\": {\"lt\": 60}},\n", " ).items()\n", ")\n", "print(f\"Found {len(s2_items)} Sentinel-2 L1C products\")" @@ -506,7 +507,7 @@ "Dimensions: ()\n", "Data variables:\n", " *empty*\n", - "Attributes: (3)
<xarray.Dataset> Size: 13GB\n", - "Dimensions: (x: 10980, y: 10980)\n", - "Coordinates:\n", - " * x (x) int64 88kB 600005 600015 600025 ... 709775 709785 709795\n", - " * y (y) int64 88kB 5099995 5099985 5099975 ... 4990215 4990205\n", - " spatial_ref int64 8B ...\n", - "Data variables: (12/13)\n", - " b02 (y, x) float64 964MB ...\n", - " b03 (y, x) float64 964MB ...\n", - " b04 (y, x) float64 964MB ...\n", - " b08 (y, x) float64 964MB ...\n", - " b05 (y, x) float64 964MB ...\n", - " b06 (y, x) float64 964MB ...\n", - " ... ...\n", - " b11 (y, x) float64 964MB ...\n", - " b12 (y, x) float64 964MB ...\n", - " b8a (y, x) float64 964MB ...\n", - " b01 (y, x) float64 964MB ...\n", - " b09 (y, x) float64 964MB ...\n", - " b10 (y, x) float64 964MB ...\n", - "Attributes: (26)" - ], - "text/plain": [ - "
<xarray.Dataset> Size: 310MB\n",
- "Dimensions: (time: 14, y: 1133, x: 813)\n",
- "Coordinates:\n",
- " * time (time) datetime64[ns] 112B 2025-10-02T10:18:51.025000 ... 20...\n",
+ " * x (x) float32 44kB 6e+05 6e+05 6e+05 ... 7.098e+05 7.098e+05\n",
+ " * y (y) float32 44kB 5.1e+06 5.1e+06 5.1e+06 ... 4.99e+06 4.99e+06\n",
" spatial_ref int64 8B ...\n",
- " * x (x) float64 7kB 6.716e+05 6.716e+05 ... 6.797e+05 6.797e+05\n",
- " * y (y) float64 9kB 5.052e+06 5.052e+06 ... 5.041e+06 5.041e+06\n",
- "Data variables:\n",
- " b02 (time, y, x) float64 103MB dask.array<chunksize=(1, 1133, 813), meta=np.ndarray>\n",
- " b03 (time, y, x) float64 103MB dask.array<chunksize=(1, 1133, 813), meta=np.ndarray>\n",
- " b04 (time, y, x) float64 103MB dask.array<chunksize=(1, 1133, 813), meta=np.ndarray>\n",
- "Attributes: (4)<xarray.Dataset> Size: 12GB\n", - "Dimensions: (x: 10980, y: 10980)\n", - "Coordinates:\n", - " * x (x) int64 88kB 600005 600015 600025 ... 709775 709785 709795\n", - " * y (y) int64 88kB 5099995 5099985 5099975 ... 4990215 4990205\n", - " spatial_ref int64 8B ...\n", - " band int64 8B ...\n", - "Data variables: (12/15)\n", - " b02 (y, x) float64 964MB ...\n", - " b03 (y, x) float64 964MB ...\n", - " b04 (y, x) float64 964MB ...\n", - " b08 (y, x) float64 964MB ...\n", - " b01 (y, x) float64 964MB ...\n", - " b05 (y, x) float64 964MB ...\n", - " ... ...\n", - " b12 (y, x) float64 964MB ...\n", - " b8a (y, x) float64 964MB ...\n", - " cld (y, x) uint8 121MB ...\n", - " snw (y, x) uint8 121MB ...\n", - " scl (y, x) uint8 121MB ...\n", - " b09 (y, x) float64 964MB ...\n", - "Attributes: (34)" - ], - "text/plain": [ - "
<xarray.Dataset> Size: 276MB\n",
- "Dimensions: (time: 12, y: 1133, x: 813)\n",
- "Coordinates:\n",
- " * time (time) datetime64[ns] 96B 2025-10-07T10:18:09.024000 ... 202...\n",
+ " * x (x) float32 44kB 6e+05 6e+05 6e+05 ... 7.098e+05 7.098e+05\n",
+ " * y (y) float32 44kB 5.1e+06 5.1e+06 5.1e+06 ... 4.99e+06 4.99e+06\n",
" spatial_ref int64 8B ...\n",
- " * x (x) float64 7kB 6.716e+05 6.716e+05 ... 6.797e+05 6.797e+05\n",
- " * y (y) float64 9kB 5.052e+06 5.052e+06 ... 5.041e+06 5.041e+06\n",
- "Data variables:\n",
- " b02 (time, y, x) float64 88MB dask.array<chunksize=(1, 1133, 813), meta=np.ndarray>\n",
- " b03 (time, y, x) float64 88MB dask.array<chunksize=(1, 1133, 813), meta=np.ndarray>\n",
- " b04 (time, y, x) float64 88MB dask.array<chunksize=(1, 1133, 813), meta=np.ndarray>\n",
- " scl (time, y, x) uint8 11MB dask.array<chunksize=(1, 1133, 813), meta=np.ndarray>\n",
- "Attributes: (4)<xarray.Dataset> Size: 255kB\n",
- "Dimensions: (time: 11, lat: 38, lon: 38)\n",
+ "<xarray.Dataset> Size: 232kB\n",
+ "Dimensions: (time: 10, lat: 38, lon: 38)\n",
"Coordinates:\n",
" spatial_ref int64 8B 0\n",
" * lon (lon) float64 304B 11.2 11.2 11.21 11.21 ... 11.3 11.3 11.3\n",
" * lat (lat) float64 304B 45.6 45.6 45.6 45.59 ... 45.51 45.5 45.5\n",
- " * time (time) datetime64[ns] 88B 2025-12-07T09:51:06.812310 ... 2...\n",
+ " * time (time) datetime64[ns] 80B 2026-02-22T09:54:51.548068 ... 2...\n",
"Data variables:\n",
- " oa01_radiance (time, lat, lon) float64 127kB dask.array<chunksize=(1, 38, 38), meta=np.ndarray>\n",
- " oa02_radiance (time, lat, lon) float64 127kB dask.array<chunksize=(1, 38, 38), meta=np.ndarray>\n",
- "Attributes: (13)