search

microsoft / PlanetaryComputer

Issues, discussions, and information about the Microsoft Planetary Computer
https://planetarycomputer.microsoft.com/
MIT License
176 stars 6 forks source link

ODC-stac can't access Landsat-5 image piexl data #365

Open 765302995 opened 1 week ago

765302995 commented 1 week ago

code: for i in landsat_5(this is a collecion): data = odc.stac.stac_load( [i], bands=bands_of_interest, bbox=bbox_of_interest ).isel(time=0)

output information:

L5 processing: 93%|█████████▎| 644/690 [04:09<00:05, 8.69it/s]Aborting load due to failure while reading: https://landsateuwest.blob.core.windows.net/landsat-c2/level-2/standard/tm/1985/053/014/LT05_L2SP_053014_19850909_20200918_02_T1/LT05_L2SP_053014_19850909_20200918_02_T1_SR_B4.TIF?st=2024-06-25T05%3A28%3A37Z&se=2024-06-26T06%3A13%3A37Z&sp=rl&sv=2024-05-04&sr=c&skoid=9c8ff44a-6a2c-4dfb-b298-1c9212f64d9a&sktid=72f988bf-86f1-41af-91ab-2d7cd011db47&skt=2024-06-26T05%3A24%3A50Z&ske=2024-07-03T05%3A24%3A50Z&sks=b&skv=2024-05-04&sig=7HkxwNKSHv8vB5qnbvf9f5YqoQ/2aEziYGztk%2BSikVU%3D:1

L5 processing: 94%|█████████▎| 646/690 [04:09<00:05, 8.79it/s]Aborting load due to failure while reading: https://landsateuwest.blob.core.windows.net/landsat-c2/level-2/standard/tm/1985/053/014/LT05_L2SP_053014_19850824_20200918_02_T1/LT05_L2SP_053014_19850824_20200918_02_T1_SR_B4.TIF?st=2024-06-25T05%3A28%3A37Z&se=2024-06-26T06%3A13%3A37Z&sp=rl&sv=2024-05-04&sr=c&skoid=9c8ff44a-6a2c-4dfb-b298-1c9212f64d9a&sktid=72f988bf-86f1-41af-91ab-2d7cd011db47&skt=2024-06-26T05%3A24%3A50Z&ske=2024-07-03T05%3A24%3A50Z&sks=b&skv=2024-05-04&sig=7HkxwNKSHv8vB5qnbvf9f5YqoQ/2aEziYGztk%2BSikVU%3D:1 Aborting load due to failure while reading: https://landsateuwest.blob.core.windows.net/landsat-c2/level-2/standard/tm/1985/053/014/LT05_L2SP_053014_19850808_20200918_02_T1/LT05_L2SP_053014_19850808_20200918_02_T1_SR_B4.TIF?st=2024-06-25T05%3A28%3A37Z&se=2024-06-26T06%3A13%3A37Z&sp=rl&sv=2024-05-04&sr=c&skoid=9c8ff44a-6a2c-4dfb-b298-1c9212f64d9a&sktid=72f988bf-86f1-41af-91ab-2d7cd011db47&skt=2024-06-26T05%3A24%3A50Z&ske=2024-07-03T05%3A24%3A50Z&sks=b&skv=2024-05-04&sig=7HkxwNKSHv8vB5qnbvf9f5YqoQ/2aEziYGztk%2BSikVU%3D:1

TomAugspurger commented 1 week ago

Thanks for opening this issue. A few questions:

  1. Can you show a full reproducible example?
  2. Are you able to load just those assets by themselves?
  3. Is this reproducible, or was that a temporary networking issue?
765302995 commented 1 week ago

1.Yes, this is an example: import pandas as pd import numpy as np from tqdm import tqdm import pystac_client import planetary_computer import odc.stac

import matplotlib.pyplot as plt

from pystac.extensions.eo import EOExtension as eo catalog = pystac_client.Client.open( "https://planetarycomputer.microsoft.com/api/stac/v1", modifier=planetary_computer.sign_inplace, ) bbox_of_interest = [ -84.751904,39.571318 , -84.746904, 39.576318] time_of_interest = "1982-01-01/2022-12-21" data_value_L5 = pd.DataFrame(columns=['time', 'blue', 'green','red', 'nir08', 'swir16', 'swir22', "lwir11" ]) for i in items_L5: bands_of_interest = ["nir08", "red", "green", "blue", "swir16","swir22", "lwir11" ] data = odc.stac.stac_load( [i], bands=bands_of_interest, bbox=bbox_of_interest ).isel(time=0)

  1. yes, I can single load those assets by themselves.
  2. is a reproducible, you can try above example code.