ncsuseal / bayesian_lsp Goto Github PK
View Code? Open in Web Editor NEWA Bayesian hierarchical model that quantifies long-term annual land surface phenology from sparse time series of vegetation indices.
License: Other
A Bayesian hierarchical model that quantifies long-term annual land surface phenology from sparse time series of vegetation indices.
License: Other
There are cases in which users want to fit BLSP to the 6-parameter double-logistic function (e.g., Fisher et al 2006) as their VI data may not from temperate deciduous forests. Actually, the BLSP model can be generalized to fit any curve functions. So, the package needs to be adapted for users to easily customize their curve function and prior information.
For ian
FitBLSP
for the fitted
outputWe are interested in the phenological information extraction method you presented. However, when we change the data of id from 1 to 3 or 5, the following message is displayed. What should we do?
Initialize model...
Estimate phenometrics...
Error in cbind(m1, c(samp[[1]][, paste0("m1", "[", i, "]")], samp[[2]][, :
object 'samp' not found
library(blsp)
library(data.table)
sample_file <- fread("sampleData.csv")
single_ts <- sample_file[id == 0, .(date, evi2)]
single_ts <- na.omit(single_ts)
head(single_ts)
avg_fit <- FitAvgModel(
date_vec = single_ts$date,
vi_vec = single_ts$evi2
)
results <- FitBLSP(
date_vec = single_ts$date,
vi_vec = single_ts$evi2,
initValues = avg_fit,
verbose = TRUE
)
head(results$phenos)
I would like to consult you on this issue, please reply to me at your convenience and thank you for your excellent work.
Recommend to add some kind of counter, either progressbar or something else to FitBLSP() because otherwise for users we're not sure how long the code needs to run
Some users commented that they want to have more phenometrics rather than the default SOS and EOS. Since the annual VI curves are estimated, more detailed phenometrics can be derived from the curves, such as Greenup, MidGreenup, Maturity, Peak, etc (following MCD12Q2). Provide a function in the package for users to specify whether they want SOS & EOS or the MCD12Q2-like phenometrics.
I will be adding more things as I think of them.
These questions are for the test data.
FitBLSP
, it doesn't seem like the QA column is inherently needed?Other
FitBLSP
I feel like it might be easier for users to find the usage of the package if we add examples to README.md. But for now, I'm not sure how to write R Markdown in the README.md, or it just can't be done that way and we have to run an R Markdown to get results and pictures and put them in the README.md manually.
I'm adding recommendations here to keep track of them. I don't want to implement these myself without asking for input.
I'll add to this list as we get closer to having a final package and documentation.
test_ts
to something like landsatEVI
or testData
or something? Maybe not necessary, I was just thinking of having the name be more informative since ts
is an abbreviation that not everyone usestest_ts
only needs columns date
, renameall_evi
to evi
, rename all_qa
to pixelQA
, rename all_from
to source
, and snow
.FitBLSP()
when having the plot, so might be easier for rendering the vignette to just have maybe 10 years worth of data and then at the end I can add a png of what the full thing looks like.
Fit_BLSP
where weights_vec
can be TRUE/FALSEweights_vec
is not numeric or is outside of range [0,1]weights_vec
is not the same length as vi_vec
A declarative, efficient, and flexible JavaScript library for building user interfaces.
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
An Open Source Machine Learning Framework for Everyone
The Web framework for perfectionists with deadlines.
A PHP framework for web artisans
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
Some thing interesting about web. New door for the world.
A server is a program made to process requests and deliver data to clients.
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
Some thing interesting about visualization, use data art
Some thing interesting about game, make everyone happy.
We are working to build community through open source technology. NB: members must have two-factor auth.
Open source projects and samples from Microsoft.
Google ❤️ Open Source for everyone.
Alibaba Open Source for everyone
Data-Driven Documents codes.
China tencent open source team.