Sentinel-5 P data usage

Đang mở

#303 mở ngày 20/3/2025 5 bình luận 0 reaction 0 người được giao Xem trên GitHub

Chưa có ai nhận issue này.

Đánh giá

Độ khó: 4/5
Thời gian dự kiến: 3-5 ngày
Mức phù hợp với người mới: 35/100

Loại issue: Lỗi
Độ rõ ràng: Cần làm rõ
Mức độ hoạt động: Đình trệ
Công nghệ: jupyter-notebook, python
Lĩnh vực: data

Hướng nghiên cứu

Bắt đầu với mã notebook tìm kiếm Planetary Computer STAC API và truyền các asset ch4 vào stackstac.stack; kiểm tra các URL blob Sentinel-5P được tạo và các phản hồi 404 được báo cáo. Xác nhận liệu các liên kết asset STAC có hợp lệ hay không và xác định hành vi dự kiến khi stack được persist và cắt mà không có lỗi khi mở asset.

Do mô hình lập chỉ mục viết ra từ nội dung của issue.

Mô tả

Hi! I've been trying to use sentinel-5 data to map CH4 emissions over my area of interest.

I replicated pretty much the code that i previously had for sentinel 1 and sentinel 2 processing:

sentinel_search_url = "https://planetarycomputer.microsoft.com/api/stac/v1"
sentinel_stac_client = pystac_client.Client.open(sentinel_search_url)

items = sentinel_stac_client.search(
    bbox=bounds,
    collections=["sentinel-5p-l2-netcdf"],
    datetime=f"{start_date}/{end_date}",
    query={"s5p:processing_mode": {"eq": "OFFL"}, "s5p:product_name": {"eq": "ch4"}},
).item_collection()

print(f"Total number of items found: {len(items)}")

scene_dates = sorted([item.datetime for item in items])

selected_scenes = [scene_dates[0]]
for i in range(1, len(scene_dates)):
    if (scene_dates[i] - selected_scenes[-1]).days >= min_time_step:
        selected_scenes.append(scene_dates[i])

items = [
    item for item in items if item.datetime in selected_scenes
]

print(f"Number of items after minimum time step filtering: {len(items)}")

print("\nSelected scenes after filtering:")
for item in items:
    print(f"Scene date: {item.datetime}")


(i had to provide a proj:epsg field as this data doesn't have it)

for item in items:
    # Access the asset metadata
    asset = item.assets['ch4']
    
    # Add proj:epsg metadata if it doesn't exist
    if 'proj:epsg' not in asset.extra_fields:
        asset.extra_fields['proj:epsg'] = 4326

sentinel_stack = stackstac.stack(
    items,
    assets=["ch4"],
    gdal_env=stackstac.DEFAULT_GDAL_ENV.updated({
        'GDAL_HTTP_MAX_RETRY': 3,
        'GDAL_HTTP_RETRY_DELAY': 5,
        'AWS_NO_SIGN_REQUEST': 'YES',
    }),
    epsg=4326,
    resolution=0.0629, 
    chunksize=(1, 1, 1000, 1000),
).to_dataset(dim='band')

# This length number represents the number  of assets (bands) that are to be extracted
# len(sentinel_stack)

sentinel_stack

# Remove attributes that are not time, y or x 
sentinel_stack = sentinel_stack.drop_vars([c for c in sentinel_stack.coords if not (c in ['time', 'y', 'x',])])
sentinel_stack

sentinel_stack_float32 = sentinel_stack.astype('float32')
sentinel_stack_float32

crs = "EPSG:4326"
crs_number = crs[5:]
sentinel_stack_float32 = sentinel_stack_float32.rio.write_crs(fr"{crs}", inplace=True)
sentinel_stack_float32

sentinel_stack_clipped = sentinel_stack_float32.rio.clip(aoi.geometry, aoi.crs, all_touched=True, drop=True)
with ProgressBar():
    sentinel_stack_clipped = sentinel_stack_clipped.persist()

Here, i get:

And when i try to save the datarray by computing the mean over this clipped area, i get a tiff file with no results

Ngôn ngữ chính: Jupyter Notebook
Star: 455
Fork: 225
Chỉ số merge pull request: Không có pull request nào được merge trong 30 ngày

Hướng dẫn đóng góp

Mở hướng dẫn đóng góp

Bắt đầu từ đâu

Đọc hết issue, rồi đọc hướng dẫn đóng góp của dự án.
Bình luận trên issue rằng bạn sẽ nhận — tránh hai người làm cùng một việc.
Fork repository và làm thay đổi trên một nhánh.
Mở pull request có tham chiếu số hiệu của issue.

Issue khác của microsoft/PlanetaryComputerExamples

which version of IMERG? Đang mở

Độ khó 1/5 Dưới một giờ Mức phù hợp với người mới 68/100

microsoft/PlanetaryComputerExamples#126 · Mở ngày 1/2/2022
Question: landsat-c2-l2 scaling & NDVI Đang mở

Độ khó 3/5 1-2 ngày Mức phù hợp với người mới 35/100

microsoft/PlanetaryComputerExamples#310 · Mở ngày 13/11/2025
Data Request: Sentinel- 2 Level 2 Angular Data Đang mở

Độ khó 5/5 Hơn một tuần Mức phù hợp với người mới 35/100

microsoft/PlanetaryComputerExamples#308 · Mở ngày 6/8/2025
Question regarding MOSAIKS Implementation Đang mở

Độ khó 4/5 3-5 ngày Mức phù hợp với người mới 30/100

microsoft/PlanetaryComputerExamples#307 · 4 bình luận · Mở ngày 31/7/2025
partition vs. fips confusion in census-data tutorial Đang mở

Độ khó 2/5 1-3 giờ Mức phù hợp với người mới 45/100

microsoft/PlanetaryComputerExamples#302 · Mở ngày 9/2/2025

Tất cả issue của microsoft/PlanetaryComputerExamples

Issue tương tự

[Datastream to BigQuery] NullPointerException in DataStreamIO when eventType is null Đang mở

bug needs triage p2

Độ khó 2/5 1-3 giờ Mức phù hợp với người mới 78/100

GoogleCloudPlatform/DataflowTemplates#4273 · 1 bình luận · Mở ngày 21/9/2026
Schema regression: derivation conversion_failure reverted Đang mở

Độ khó 1/5 Dưới một giờ Mức phù hợp với người mới 90/100

elong0527/yamaa#708 · Mở ngày 21/9/2026
Metadata correction for 2025.nlp4dh-1.26 Đang mở

correction metadata

Độ khó 1/5 Dưới một giờ Mức phù hợp với người mới 72/100

acl-org/acl-anthology#10025 · 1 bình luận · Mở ngày 20/9/2026
tl.feature_signature returns an object io.validate rejects Đang mở

bug

Độ khó 2/5 1-3 giờ Mức phù hợp với người mới 88/100

scverse/mantispy#103 · Mở ngày 20/9/2026
indexing issue in correct_cte Đang mở

Độ khó 2/5 1-3 giờ Mức phù hợp với người mới 76/100

desihub/desispec#2834 · Mở ngày 20/9/2026

Thêm issue về Data Engineering