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

PRODEM dates are not constant #46

Open
mankoff opened this issue Jun 12, 2024 · 0 comments
Open

PRODEM dates are not constant #46

mankoff opened this issue Jun 12, 2024 · 0 comments

Comments

@mankoff
Copy link
Member

mankoff commented Jun 12, 2024

**** PRODEM
#+BEGIN_SRC bash :exports both :results verbatim
MSG_OK "dh/dt"
g.mapset -c PRODEM
r.mask -r
f=$(ls ${DATADIR}/PRODEM/PRODEM??.tif | head -n1) # debug
for f in $(ls ${DATADIR}/PRODEM/PRODEM??.tif); do
y=20$(echo ${f: -6:2})
r.in.gdal -r input=${f} output=DEM_${y} band=1
# r.in.gdal -r input=${f} output=var_${y} band=2
# r.in.gdal -r input=${f} output=dh_${y} band=3
# r.in.gdal -r input=${f} output=time_${y} band=4
# r.univar -g time_2019 # mean = DOI 213 = 01 Aug
done
g.region raster=DEM_2019 -pa
#+END_SRC

Currently PRODEMs are assumed to be tied to the median date (August 1). This should be improved by using the date associated with each pixel. Lining up the dates between the earliest PRODEM and the earlier Khan DEMs should also be done more carefully.

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

1 participant