The hardware and bandwidth for this mirror is donated by dogado GmbH, the Webhosting and Full Service-Cloud Provider. Check out our Wordpress Tutorial.
If you wish to report a bug, or if you are interested in having us mirror your free-software or open-source project, please feel free to contact us at mirror[@]dogado.de.

DyMEP

DyMEP a Dynamic Multi-Environment Phenology Model. This phenology model allows to model crop phenology according to environmental covariates. One can chose, or ask the model, for the best combination of environmental covariates for a specific crop phenology macro-phase. The model was trained using a large dataset of crop phenology derived from the German weather service (DWD): https://opendata.dwd.de. more details can be found here: DOIs

Available Crops are: -Winter wheat -Winter Rye -Winter Barley -Spring Wheat -Spring Barley -Rapeseed -Oat -Maize -Green Peas -Green Beans

used covariates (daily resolution): - minimal temperature - maximal temperature - mean temperature - SPI (30-days standardized precipitation index) - relative humidity - global radiation - VPD

Example

This is a basic example which shows you how to solve a common problem:

library(DyMEP)
## basic example code
library(DyMEP)
## basic example code
# check what exists
available_phases <- available_crops_and_phases()
# check what covairates are implemented in the model
available_covariates <- available_environmental_covariates()


#what is the best environmental covariates for one or multiple phases?

best_DyMEP_model(env_covariates = c("tas","tasmin","VPD","SPI",
                                    "global_radiation","tasmax","RH"),
                 pheno_phases = c("sowing-emergence","jointing-heading"),
                        crop_abbrev = "WB")


phase_covariate_list <- best_DyMEP_model(env_covariates = c("tas","tasmin",
                                "VPD","SPI","global_radiation","tasmax","RH"),
                 pheno_phases = c("sowing-emergence","emergence-jointing",
                                  "jointing-heading"),
                 crop_abbrev = "WB",
                 output_list_for_prediction = T)

#create a list of wanted phases and corresponding environmental covariates
# we chose the suggested models, but we add the phase emergence-booting as well.
# (Attention, the model will not run if you provide gaps in the phenology!
#it is no problem to run a subset of phases but they need to be consecutive!)
phase_covariate_list <- list("sowing-emergence" = c("tasmin","VPD","SPI"),
                        "emergence-jointing"= c("tas","tasmin","VPD","SPI"),
                     "jointing-heading" = c("global_radiation","tas","SPI"))

#create dummy environmental data
environmental_data <- data.frame("DATE"=seq.Date(from = as.Date("2021-01-01"),
                                              to = as.Date("2022-12-31"),by=1),
                             "tas"=runif(730,min=-10,max=40),
                             "RH"=runif(730,min=0,max=100),
                             "tasmin"=runif(730,min=-10,max=40),
                             "tasmax"=runif(730,min=-10,max=40),
                             "VPD" = runif(730,min=-10,max=40),
                             "SPI"= runif(730,min=-1,max=4),
                             "global_radiation"= runif(730,min=0,max=10000))
# predict phenology, we that the crop was sown at the "2022-01-01" 

pheno_phase_prediction(phase_covariate_list = phase_covariate_list,
                       environmental_data = environmental_data,
                       phase_starting_date =as.Date("2022-01-01"),
                       crop_abbrev = "WR")


detailed_output <- pheno_phase_prediction(phase_covariate_list =
                                            phase_covariate_list,
                       environmental_data = environmental_data,
                       phase_starting_date =as.Date("2022-01-01"),
                       crop_abbrev = "WR",
                       output_type = "detailed_information")


# get overview plot of the prediction
DyMEP_prediction_visualizer(detailed_output)
# check the DRC curves of the used model
DyMEP_DRC_visualizer(detailed_output)

# get summary
summary(detailed_output)



# output are the end dates of the stages (attention in the example
# we use random generated dummy data, the resulting dates are therefore 
# not extremely meaningful)

These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.
Health stats visible at Monitor.