Skip to content

Materialize ds before snow-index calcs to fix read-only buffer error #83

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.

Sign up for GitHub

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

Closed
ZachHoppinen wants to merge 2 commits into from
Closed

Materialize ds before snow-index calcs to fix read-only buffer error #83

Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
5422d0f
load dataset before snow-index calcs to avoid read-only buffer error
ZachHoppinen May 9, 2026
6667b7b
deep-copy after .load() so already-loaded read-only buffers become wr…
ZachHoppinen May 9, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
9 changes: 9 additions & 0 deletions spicy_snow/retrieval.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -164,6 +164,15 @@ def retrieve_snow_depth(aoi: shapely.geometry.Polygon,
# clip outlier values of backscatter to overall mean
ds = s1_clip_outliers(ds)

# The snow-index functions do in-place .loc[time=...] writes; without
# forcing fresh writable in-memory arrays, the second orbit's write
# fails with "assignment destination is read-only".
# .load() alone is insufficient: it materializes lazy/dask arrays but
# does NOT deep-copy already-in-memory arrays whose buffers are
# read-only (observed on fairbanks 2024_2025). Combine with
# .copy(deep=True) to guarantee writable backing.
ds = ds.load().copy(deep=True)

# -- Calulating Snow Index -- #
log.info("Calculating snow index")
# calculate delta CR and delta VV
Expand Down
Loading