{ "cells": [ { "cell_type": "markdown", "id": "a61446d2", "metadata": {}, "source": [ "# Landsat Time Series" ] }, { "cell_type": "code", "execution_count": 1, "id": "630ba918", "metadata": {}, "outputs": [], "source": [ "import os\n", "from pcxarray import pc_query, prepare_timeseries\n", "from distributed import Client # dask distributed client\n", "import matplotlib.pyplot as plt\n", "from IPython.display import HTML\n", "from matplotlib.animation import FuncAnimation\n", "from matplotlib.colors import LinearSegmentedColormap\n", "from shapely import Polygon\n", "import geopandas as gpd\n", "import numpy as np" ] }, { "cell_type": "code", "execution_count": 2, "id": "bd186f11", "metadata": {}, "outputs": [], "source": [ "geom = Polygon([\n", " (-91.3011932, 33.852768),\n", " (-90.9098053, 33.8413567),\n", " (-90.9194183, 33.5069017),\n", " (-91.2957001, 33.509193),\n", " (-91.3005066, 33.8579026),\n", "])\n", "\n", "gdf = gpd.GeoDataFrame(geometry=[geom], crs=\"EPSG:4326\")\n", "gdf = gdf.to_crs(32615) # UTM 15N\n", "selected_geom = gdf.geometry.union_all().envelope" ] }, { "cell_type": "code", "execution_count": 3, "id": "a4d300a7", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Found 5439 items\n" ] } ], "source": [ "items_gdf = pc_query(\n", " collections=\"landsat-c2-l2\",\n", " geometry=selected_geom,\n", " datetime=\"1984/2024\",\n", " crs=gdf.crs,\n", " # query={\"eo:cloud_cover\": {\"lt\": 5}} # OPTIONAL: filter by cloud cover\n", ")\n", "print(f\"Found {len(items_gdf)} items\")" ] }, { "cell_type": "code", "execution_count": null, "id": "a568e33f", "metadata": {}, "outputs": [], "source": [ "client = Client(n_workers=os.cpu_count(), memory_limit=\"4GB\")" ] }, { "cell_type": "code", "execution_count": 5, "id": "68bae4ac", "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "Constructing dask computation graph: 92%|█████████▏| 2520/2740 [05:01<00:26, 8.17timestep/s]/home/dh2306/pcxarray/src/pcxarray/io.py:92: UserWarning: Retrying loading raster from COG to error: '/vsicurl/https://landsateuwest.blob.core.windows.net/landsat-c2/level-2/standard/oli-tirs/2022/024/036/LC08_L2SP_024036_20220616_20220629_02_T1/LC08_L2SP_024036_20220616_20220629_02_T1_SR_B3.TIF?st=2025-07-08T14%3A11%3A42Z&se=2025-07-09T14%3A56%3A42Z&sp=rl&sv=2024-05-04&sr=c&skoid=9c8ff44a-6a2c-4dfb-b298-1c9212f64d9a&sktid=72f988bf-86f1-41af-91ab-2d7cd011db47&skt=2025-07-09T05%3A09%3A40Z&ske=2025-07-16T05%3A09%3A40Z&sks=b&skv=2024-05-04&sig=omp2Y4M1akRciljB2%2BUBT7NBODgMSQPcevDEFtrqDL4%3D' not recognized as being in a supported file format.. Attempt 1/5\n", " warn(f\"Retrying loading raster from COG to error: {e}. Attempt {_retries + 1}/{max_retries}\")\n", "/home/dh2306/pcxarray/src/pcxarray/io.py:92: UserWarning: Retrying loading raster from COG to error: '/vsicurl/https://landsateuwest.blob.core.windows.net/landsat-c2/level-2/standard/oli-tirs/2022/024/036/LC08_L2SP_024036_20220616_20220629_02_T1/LC08_L2SP_024036_20220616_20220629_02_T1_SR_B3.TIF?st=2025-07-08T14%3A11%3A42Z&se=2025-07-09T14%3A56%3A42Z&sp=rl&sv=2024-05-04&sr=c&skoid=9c8ff44a-6a2c-4dfb-b298-1c9212f64d9a&sktid=72f988bf-86f1-41af-91ab-2d7cd011db47&skt=2025-07-09T05%3A09%3A40Z&ske=2025-07-16T05%3A09%3A40Z&sks=b&skv=2024-05-04&sig=omp2Y4M1akRciljB2%2BUBT7NBODgMSQPcevDEFtrqDL4%3D' not recognized as being in a supported file format.. Attempt 2/5\n", " warn(f\"Retrying loading raster from COG to error: {e}. Attempt {_retries + 1}/{max_retries}\")\n", "Constructing dask computation graph: 93%|█████████▎| 2544/2740 [05:04<00:25, 7.78timestep/s]/home/dh2306/pcxarray/src/pcxarray/io.py:92: UserWarning: Retrying loading raster from COG to error: '/vsicurl/https://landsateuwest.blob.core.windows.net/landsat-c2/level-2/standard/oli-tirs/2022/024/037/LC08_L2SP_024037_20220702_20220707_02_T2/LC08_L2SP_024037_20220702_20220707_02_T2_QA_PIXEL.TIF?st=2025-07-08T14%3A11%3A42Z&se=2025-07-09T14%3A56%3A42Z&sp=rl&sv=2024-05-04&sr=c&skoid=9c8ff44a-6a2c-4dfb-b298-1c9212f64d9a&sktid=72f988bf-86f1-41af-91ab-2d7cd011db47&skt=2025-07-09T05%3A09%3A40Z&ske=2025-07-16T05%3A09%3A40Z&sks=b&skv=2024-05-04&sig=omp2Y4M1akRciljB2%2BUBT7NBODgMSQPcevDEFtrqDL4%3D' not recognized as being in a supported file format.. Attempt 1/5\n", " warn(f\"Retrying loading raster from COG to error: {e}. Attempt {_retries + 1}/{max_retries}\")\n", "/home/dh2306/pcxarray/src/pcxarray/io.py:92: UserWarning: Retrying loading raster from COG to error: '/vsicurl/https://landsateuwest.blob.core.windows.net/landsat-c2/level-2/standard/oli-tirs/2022/024/036/LC08_L2SP_024036_20220616_20220629_02_T1/LC08_L2SP_024036_20220616_20220629_02_T1_SR_B3.TIF?st=2025-07-08T14%3A11%3A42Z&se=2025-07-09T14%3A56%3A42Z&sp=rl&sv=2024-05-04&sr=c&skoid=9c8ff44a-6a2c-4dfb-b298-1c9212f64d9a&sktid=72f988bf-86f1-41af-91ab-2d7cd011db47&skt=2025-07-09T05%3A09%3A40Z&ske=2025-07-16T05%3A09%3A40Z&sks=b&skv=2024-05-04&sig=omp2Y4M1akRciljB2%2BUBT7NBODgMSQPcevDEFtrqDL4%3D' not recognized as being in a supported file format.. Attempt 3/5\n", " warn(f\"Retrying loading raster from COG to error: {e}. Attempt {_retries + 1}/{max_retries}\")\n", "/home/dh2306/pcxarray/src/pcxarray/io.py:92: UserWarning: Retrying loading raster from COG to error: '/vsicurl/https://landsateuwest.blob.core.windows.net/landsat-c2/level-2/standard/oli-tirs/2022/024/037/LC08_L2SP_024037_20220702_20220707_02_T2/LC08_L2SP_024037_20220702_20220707_02_T2_QA_PIXEL.TIF?st=2025-07-08T14%3A11%3A42Z&se=2025-07-09T14%3A56%3A42Z&sp=rl&sv=2024-05-04&sr=c&skoid=9c8ff44a-6a2c-4dfb-b298-1c9212f64d9a&sktid=72f988bf-86f1-41af-91ab-2d7cd011db47&skt=2025-07-09T05%3A09%3A40Z&ske=2025-07-16T05%3A09%3A40Z&sks=b&skv=2024-05-04&sig=omp2Y4M1akRciljB2%2BUBT7NBODgMSQPcevDEFtrqDL4%3D' not recognized as being in a supported file format.. Attempt 2/5\n", " warn(f\"Retrying loading raster from COG to error: {e}. Attempt {_retries + 1}/{max_retries}\")\n", "Constructing dask computation graph: 94%|█████████▎| 2568/2740 [05:07<00:20, 8.20timestep/s]/home/dh2306/pcxarray/src/pcxarray/io.py:92: UserWarning: Retrying loading raster from COG to error: '/vsicurl/https://landsateuwest.blob.core.windows.net/landsat-c2/level-2/standard/oli-tirs/2022/024/037/LC08_L2SP_024037_20221107_20221115_02_T1/LC08_L2SP_024037_20221107_20221115_02_T1_QA_PIXEL.TIF?st=2025-07-08T14%3A11%3A42Z&se=2025-07-09T14%3A56%3A42Z&sp=rl&sv=2024-05-04&sr=c&skoid=9c8ff44a-6a2c-4dfb-b298-1c9212f64d9a&sktid=72f988bf-86f1-41af-91ab-2d7cd011db47&skt=2025-07-09T05%3A10%3A14Z&ske=2025-07-16T05%3A10%3A14Z&sks=b&skv=2024-05-04&sig=V%2BGV6yKA6t1RJXaOGR%2BEUPCzbxrbWQ3AdgMvrRVuQcs%3D' not recognized as being in a supported file format.. Attempt 1/5\n", " warn(f\"Retrying loading raster from COG to error: {e}. Attempt {_retries + 1}/{max_retries}\")\n", "/home/dh2306/pcxarray/src/pcxarray/io.py:92: UserWarning: Retrying loading raster from COG to error: '/vsicurl/https://landsateuwest.blob.core.windows.net/landsat-c2/level-2/standard/oli-tirs/2022/024/037/LC08_L2SP_024037_20220702_20220707_02_T2/LC08_L2SP_024037_20220702_20220707_02_T2_QA_PIXEL.TIF?st=2025-07-08T14%3A11%3A42Z&se=2025-07-09T14%3A56%3A42Z&sp=rl&sv=2024-05-04&sr=c&skoid=9c8ff44a-6a2c-4dfb-b298-1c9212f64d9a&sktid=72f988bf-86f1-41af-91ab-2d7cd011db47&skt=2025-07-09T05%3A09%3A40Z&ske=2025-07-16T05%3A09%3A40Z&sks=b&skv=2024-05-04&sig=omp2Y4M1akRciljB2%2BUBT7NBODgMSQPcevDEFtrqDL4%3D' not recognized as being in a supported file format.. Attempt 3/5\n", " warn(f\"Retrying loading raster from COG to error: {e}. Attempt {_retries + 1}/{max_retries}\")\n", "/home/dh2306/pcxarray/src/pcxarray/io.py:92: UserWarning: Retrying loading raster from COG to error: '/vsicurl/https://landsateuwest.blob.core.windows.net/landsat-c2/level-2/standard/oli-tirs/2022/024/037/LC08_L2SP_024037_20221107_20221115_02_T1/LC08_L2SP_024037_20221107_20221115_02_T1_QA_PIXEL.TIF?st=2025-07-08T14%3A11%3A42Z&se=2025-07-09T14%3A56%3A42Z&sp=rl&sv=2024-05-04&sr=c&skoid=9c8ff44a-6a2c-4dfb-b298-1c9212f64d9a&sktid=72f988bf-86f1-41af-91ab-2d7cd011db47&skt=2025-07-09T05%3A10%3A14Z&ske=2025-07-16T05%3A10%3A14Z&sks=b&skv=2024-05-04&sig=V%2BGV6yKA6t1RJXaOGR%2BEUPCzbxrbWQ3AdgMvrRVuQcs%3D' not recognized as being in a supported file format.. Attempt 2/5\n", " warn(f\"Retrying loading raster from COG to error: {e}. Attempt {_retries + 1}/{max_retries}\")\n", "Constructing dask computation graph: 95%|█████████▍| 2592/2740 [05:10<00:18, 8.09timestep/s]/home/dh2306/pcxarray/src/pcxarray/io.py:92: UserWarning: Retrying loading raster from COG to error: '/vsicurl/https://landsateuwest.blob.core.windows.net/landsat-c2/level-2/standard/oli-tirs/2022/024/036/LC08_L2SP_024036_20220616_20220629_02_T1/LC08_L2SP_024036_20220616_20220629_02_T1_SR_B3.TIF?st=2025-07-08T14%3A11%3A42Z&se=2025-07-09T14%3A56%3A42Z&sp=rl&sv=2024-05-04&sr=c&skoid=9c8ff44a-6a2c-4dfb-b298-1c9212f64d9a&sktid=72f988bf-86f1-41af-91ab-2d7cd011db47&skt=2025-07-09T05%3A09%3A40Z&ske=2025-07-16T05%3A09%3A40Z&sks=b&skv=2024-05-04&sig=omp2Y4M1akRciljB2%2BUBT7NBODgMSQPcevDEFtrqDL4%3D' not recognized as being in a supported file format.. Attempt 4/5\n", " warn(f\"Retrying loading raster from COG to error: {e}. Attempt {_retries + 1}/{max_retries}\")\n", "/home/dh2306/pcxarray/src/pcxarray/io.py:92: UserWarning: Retrying loading raster from COG to error: '/vsicurl/https://landsateuwest.blob.core.windows.net/landsat-c2/level-2/standard/oli-tirs/2022/024/037/LC08_L2SP_024037_20221107_20221115_02_T1/LC08_L2SP_024037_20221107_20221115_02_T1_QA_PIXEL.TIF?st=2025-07-08T14%3A11%3A42Z&se=2025-07-09T14%3A56%3A42Z&sp=rl&sv=2024-05-04&sr=c&skoid=9c8ff44a-6a2c-4dfb-b298-1c9212f64d9a&sktid=72f988bf-86f1-41af-91ab-2d7cd011db47&skt=2025-07-09T05%3A10%3A14Z&ske=2025-07-16T05%3A10%3A14Z&sks=b&skv=2024-05-04&sig=V%2BGV6yKA6t1RJXaOGR%2BEUPCzbxrbWQ3AdgMvrRVuQcs%3D' not recognized as being in a supported file format.. Attempt 3/5\n", " warn(f\"Retrying loading raster from COG to error: {e}. Attempt {_retries + 1}/{max_retries}\")\n", "/home/dh2306/pcxarray/src/pcxarray/io.py:92: UserWarning: Retrying loading raster from COG to error: '/vsicurl/https://landsateuwest.blob.core.windows.net/landsat-c2/level-2/standard/oli-tirs/2022/024/037/LC08_L2SP_024037_20220702_20220707_02_T2/LC08_L2SP_024037_20220702_20220707_02_T2_QA_PIXEL.TIF?st=2025-07-08T14%3A11%3A42Z&se=2025-07-09T14%3A56%3A42Z&sp=rl&sv=2024-05-04&sr=c&skoid=9c8ff44a-6a2c-4dfb-b298-1c9212f64d9a&sktid=72f988bf-86f1-41af-91ab-2d7cd011db47&skt=2025-07-09T05%3A09%3A40Z&ske=2025-07-16T05%3A09%3A40Z&sks=b&skv=2024-05-04&sig=omp2Y4M1akRciljB2%2BUBT7NBODgMSQPcevDEFtrqDL4%3D' not recognized as being in a supported file format.. Attempt 4/5\n", " warn(f\"Retrying loading raster from COG to error: {e}. Attempt {_retries + 1}/{max_retries}\")\n", "Constructing dask computation graph: 95%|█████████▌| 2616/2740 [05:13<00:15, 7.97timestep/s]/home/dh2306/pcxarray/src/pcxarray/io.py:92: UserWarning: Retrying loading raster from COG to error: '/vsicurl/https://landsateuwest.blob.core.windows.net/landsat-c2/level-2/standard/oli-tirs/2022/024/037/LC08_L2SP_024037_20221107_20221115_02_T1/LC08_L2SP_024037_20221107_20221115_02_T1_QA_PIXEL.TIF?st=2025-07-08T14%3A11%3A42Z&se=2025-07-09T14%3A56%3A42Z&sp=rl&sv=2024-05-04&sr=c&skoid=9c8ff44a-6a2c-4dfb-b298-1c9212f64d9a&sktid=72f988bf-86f1-41af-91ab-2d7cd011db47&skt=2025-07-09T05%3A10%3A14Z&ske=2025-07-16T05%3A10%3A14Z&sks=b&skv=2024-05-04&sig=V%2BGV6yKA6t1RJXaOGR%2BEUPCzbxrbWQ3AdgMvrRVuQcs%3D' not recognized as being in a supported file format.. Attempt 4/5\n", " warn(f\"Retrying loading raster from COG to error: {e}. Attempt {_retries + 1}/{max_retries}\")\n", "Constructing dask computation graph: 96%|█████████▋| 2640/2740 [05:17<00:12, 7.72timestep/s]/home/dh2306/pcxarray/src/pcxarray/processing.py:319: UserWarning: Failed to read item LC08_L2SP_024036_20220616_02_T1 with error: Failed to load raster from COG after 5 attempts: '/vsicurl/https://landsateuwest.blob.core.windows.net/landsat-c2/level-2/standard/oli-tirs/2022/024/036/LC08_L2SP_024036_20220616_20220629_02_T1/LC08_L2SP_024036_20220616_20220629_02_T1_SR_B3.TIF?st=2025-07-08T14%3A11%3A42Z&se=2025-07-09T14%3A56%3A42Z&sp=rl&sv=2024-05-04&sr=c&skoid=9c8ff44a-6a2c-4dfb-b298-1c9212f64d9a&sktid=72f988bf-86f1-41af-91ab-2d7cd011db47&skt=2025-07-09T05%3A09%3A40Z&ske=2025-07-16T05%3A09%3A40Z&sks=b&skv=2024-05-04&sig=omp2Y4M1akRciljB2%2BUBT7NBODgMSQPcevDEFtrqDL4%3D' not recognized as being in a supported file format.. Skipping this item.\n", " warn(f\"Failed to read item {item_series['id']} with error: {e}. Skipping this item.\")\n", "Constructing dask computation graph: 97%|█████████▋| 2664/2740 [05:20<00:10, 7.26timestep/s]/home/dh2306/pcxarray/src/pcxarray/processing.py:319: UserWarning: Failed to read item LC08_L2SP_024037_20220702_02_T2 with error: Failed to load raster from COG after 5 attempts: '/vsicurl/https://landsateuwest.blob.core.windows.net/landsat-c2/level-2/standard/oli-tirs/2022/024/037/LC08_L2SP_024037_20220702_20220707_02_T2/LC08_L2SP_024037_20220702_20220707_02_T2_QA_PIXEL.TIF?st=2025-07-08T14%3A11%3A42Z&se=2025-07-09T14%3A56%3A42Z&sp=rl&sv=2024-05-04&sr=c&skoid=9c8ff44a-6a2c-4dfb-b298-1c9212f64d9a&sktid=72f988bf-86f1-41af-91ab-2d7cd011db47&skt=2025-07-09T05%3A09%3A40Z&ske=2025-07-16T05%3A09%3A40Z&sks=b&skv=2024-05-04&sig=omp2Y4M1akRciljB2%2BUBT7NBODgMSQPcevDEFtrqDL4%3D' not recognized as being in a supported file format.. Skipping this item.\n", " warn(f\"Failed to read item {item_series['id']} with error: {e}. Skipping this item.\")\n", "/home/dh2306/pcxarray/src/pcxarray/io.py:92: UserWarning: Retrying loading raster from COG to error: '/vsicurl/https://landsateuwest.blob.core.windows.net/landsat-c2/level-2/standard/oli-tirs/2022/024/037/LC08_L2SP_024037_20220616_20220629_02_T1/LC08_L2SP_024037_20220616_20220629_02_T1_QA_PIXEL.TIF?st=2025-07-08T14%3A11%3A42Z&se=2025-07-09T14%3A56%3A42Z&sp=rl&sv=2024-05-04&sr=c&skoid=9c8ff44a-6a2c-4dfb-b298-1c9212f64d9a&sktid=72f988bf-86f1-41af-91ab-2d7cd011db47&skt=2025-07-09T05%3A09%3A40Z&ske=2025-07-16T05%3A09%3A40Z&sks=b&skv=2024-05-04&sig=omp2Y4M1akRciljB2%2BUBT7NBODgMSQPcevDEFtrqDL4%3D' not recognized as being in a supported file format.. Attempt 1/5\n", " warn(f\"Retrying loading raster from COG to error: {e}. Attempt {_retries + 1}/{max_retries}\")\n", "/home/dh2306/pcxarray/src/pcxarray/io.py:92: UserWarning: Retrying loading raster from COG to error: '/vsicurl/https://landsateuwest.blob.core.windows.net/landsat-c2/level-2/standard/oli-tirs/2022/024/037/LC08_L2SP_024037_20220616_20220629_02_T1/LC08_L2SP_024037_20220616_20220629_02_T1_QA_PIXEL.TIF?st=2025-07-08T14%3A11%3A42Z&se=2025-07-09T14%3A56%3A42Z&sp=rl&sv=2024-05-04&sr=c&skoid=9c8ff44a-6a2c-4dfb-b298-1c9212f64d9a&sktid=72f988bf-86f1-41af-91ab-2d7cd011db47&skt=2025-07-09T05%3A09%3A40Z&ske=2025-07-16T05%3A09%3A40Z&sks=b&skv=2024-05-04&sig=omp2Y4M1akRciljB2%2BUBT7NBODgMSQPcevDEFtrqDL4%3D' not recognized as being in a supported file format.. Attempt 2/5\n", " warn(f\"Retrying loading raster from COG to error: {e}. Attempt {_retries + 1}/{max_retries}\")\n", "Constructing dask computation graph: 98%|█████████▊| 2688/2740 [05:24<00:07, 7.19timestep/s]/home/dh2306/pcxarray/src/pcxarray/processing.py:319: UserWarning: Failed to read item LC08_L2SP_024037_20221107_02_T1 with error: Failed to load raster from COG after 5 attempts: '/vsicurl/https://landsateuwest.blob.core.windows.net/landsat-c2/level-2/standard/oli-tirs/2022/024/037/LC08_L2SP_024037_20221107_20221115_02_T1/LC08_L2SP_024037_20221107_20221115_02_T1_QA_PIXEL.TIF?st=2025-07-08T14%3A11%3A42Z&se=2025-07-09T14%3A56%3A42Z&sp=rl&sv=2024-05-04&sr=c&skoid=9c8ff44a-6a2c-4dfb-b298-1c9212f64d9a&sktid=72f988bf-86f1-41af-91ab-2d7cd011db47&skt=2025-07-09T05%3A10%3A14Z&ske=2025-07-16T05%3A10%3A14Z&sks=b&skv=2024-05-04&sig=V%2BGV6yKA6t1RJXaOGR%2BEUPCzbxrbWQ3AdgMvrRVuQcs%3D' not recognized as being in a supported file format.. Skipping this item.\n", " warn(f\"Failed to read item {item_series['id']} with error: {e}. Skipping this item.\")\n", "/home/dh2306/pcxarray/src/pcxarray/io.py:92: UserWarning: Retrying loading raster from COG to error: '/vsicurl/https://landsateuwest.blob.core.windows.net/landsat-c2/level-2/standard/oli-tirs/2022/024/037/LC08_L2SP_024037_20220616_20220629_02_T1/LC08_L2SP_024037_20220616_20220629_02_T1_QA_PIXEL.TIF?st=2025-07-08T14%3A11%3A42Z&se=2025-07-09T14%3A56%3A42Z&sp=rl&sv=2024-05-04&sr=c&skoid=9c8ff44a-6a2c-4dfb-b298-1c9212f64d9a&sktid=72f988bf-86f1-41af-91ab-2d7cd011db47&skt=2025-07-09T05%3A09%3A40Z&ske=2025-07-16T05%3A09%3A40Z&sks=b&skv=2024-05-04&sig=omp2Y4M1akRciljB2%2BUBT7NBODgMSQPcevDEFtrqDL4%3D' not recognized as being in a supported file format.. Attempt 3/5\n", " warn(f\"Retrying loading raster from COG to error: {e}. Attempt {_retries + 1}/{max_retries}\")\n", "Constructing dask computation graph: 99%|█████████▉| 2712/2740 [05:26<00:03, 7.65timestep/s]/home/dh2306/pcxarray/src/pcxarray/io.py:92: UserWarning: Retrying loading raster from COG to error: '/vsicurl/https://landsateuwest.blob.core.windows.net/landsat-c2/level-2/standard/oli-tirs/2022/024/037/LC08_L2SP_024037_20220616_20220629_02_T1/LC08_L2SP_024037_20220616_20220629_02_T1_QA_PIXEL.TIF?st=2025-07-08T14%3A11%3A42Z&se=2025-07-09T14%3A56%3A42Z&sp=rl&sv=2024-05-04&sr=c&skoid=9c8ff44a-6a2c-4dfb-b298-1c9212f64d9a&sktid=72f988bf-86f1-41af-91ab-2d7cd011db47&skt=2025-07-09T05%3A09%3A40Z&ske=2025-07-16T05%3A09%3A40Z&sks=b&skv=2024-05-04&sig=omp2Y4M1akRciljB2%2BUBT7NBODgMSQPcevDEFtrqDL4%3D' not recognized as being in a supported file format.. Attempt 4/5\n", " warn(f\"Retrying loading raster from COG to error: {e}. Attempt {_retries + 1}/{max_retries}\")\n", "Constructing dask computation graph: 100%|██████████| 2740/2740 [05:29<00:00, 8.31timestep/s]\n", "/home/dh2306/pcxarray/src/pcxarray/processing.py:319: UserWarning: Failed to read item LC08_L2SP_024037_20220616_02_T1 with error: Failed to load raster from COG after 5 attempts: '/vsicurl/https://landsateuwest.blob.core.windows.net/landsat-c2/level-2/standard/oli-tirs/2022/024/037/LC08_L2SP_024037_20220616_20220629_02_T1/LC08_L2SP_024037_20220616_20220629_02_T1_QA_PIXEL.TIF?st=2025-07-08T14%3A11%3A42Z&se=2025-07-09T14%3A56%3A42Z&sp=rl&sv=2024-05-04&sr=c&skoid=9c8ff44a-6a2c-4dfb-b298-1c9212f64d9a&sktid=72f988bf-86f1-41af-91ab-2d7cd011db47&skt=2025-07-09T05%3A09%3A40Z&ske=2025-07-16T05%3A09%3A40Z&sks=b&skv=2024-05-04&sig=omp2Y4M1akRciljB2%2BUBT7NBODgMSQPcevDEFtrqDL4%3D' not recognized as being in a supported file format.. Skipping this item.\n", " warn(f\"Failed to read item {item_series['id']} with error: {e}. Skipping this item.\")\n", "/home/dh2306/pcxarray/src/pcxarray/processing.py:643: UserWarning: Failed to process a group with error: Could not read any items that supposedly had partial overlap with the geometry. Try checking the geometry and its CRS.. Skipping this group.\n", " warn(f\"Failed to process a group with error: {e}. Skipping this group.\")\n" ] }, { "data": { "text/html": [ "
<xarray.DataArray 'reproject-536c3f7c2189476a9b7a18f4ba4fea16' (time: 2739,\n",
" band: 3,\n",
" y: 1290, x: 1209)> Size: 51GB\n",
"dask.array<rechunk-merge, shape=(2739, 3, 1290, 1209), dtype=float32, chunksize=(64, 1, 1290, 1209), chunktype=numpy.ndarray>\n",
"Coordinates:\n",
" * time (time) datetime64[ns] 22kB 1984-04-20 1984-04-29 ... 2024-12-31\n",
" * band (band) <U8 96B 'green' 'nir08' 'qa_pixel'\n",
" * y (y) float64 10kB 3.748e+06 3.748e+06 ... 3.709e+06 3.709e+06\n",
" * x (x) float64 10kB 6.572e+05 6.572e+05 ... 6.934e+05 6.934e+05\n",
" spatial_ref int32 4B 32615\n",
"Attributes:\n",
" AREA_OR_POINT: Point\n",
" scale_factor: 1.0\n",
" add_offset: 0.0\n",
" _FillValue: nan