Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

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

Open
765302995 opened this issue Jun 26, 2024 · 2 comments
Open

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

765302995 opened this issue Jun 26, 2024 · 2 comments

Comments

@765302995
Copy link

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
Copy link

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
Copy link
Author

765302995 commented Jun 29, 2024

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.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants