r - Plotting crop calendars -

-

i have csv file (crop_calendar.csv) containing information on development stages of crop in particular region. each row has following structure:

crop_name   sowing_dat    emergence_date  flowering_date  maturity_date  harvest_date

which gives example:

winter_wheat    18.08   28.08   24.06   30.07   3.08 winter_rye      18.08   28.08   15.06   23.07   29.07 spring_wheat    27.04   10.05   1.07    4.08    7.08 spring_barley   27.04   12.05   27.06   1.08    5.08

now, i'd put information in graphic looks that: crop calendar example

any idea how lots of crop (rows) , @ different locations?

here example assuming have day.of.year() of sowing , duration (in days) of 3 periods each crop , each country.

the crop calendar

#making random numbers reproducible set.seed(12345) rawdata <- expand.grid(   crop = paste("crop", letters[1:8]),    country = paste("country", letters[10:13]) ) #day.of.year of sowing rawdata$sowing <- runif(nrow(rawdata), min = 0, max = 365) #number of days until mid season rawdata$midseason <- runif(nrow(rawdata), min = 10, max = 30) #number of days until harvest rawdata$harvest <- runif(nrow(rawdata), min = 20, max = 150) #number of days until end of harvest rawdata$harvest.end <- runif(nrow(rawdata), min = 10, max = 40)  dataset <- data.frame(crop = character(0), country = character(0), period = character(0), duration = numeric(0))  #sowing around new year last.day <- rowsums(rawdata[, c("sowing", "midseason")]) if(any(last.day >= 365)){   dataset <- rbind(     dataset,     cbind(       rawdata[last.day >= 365, c("crop", "country")],       period = "sowing",       duration = last.day[last.day >= 365] - 365     )   )   dataset <- rbind(     dataset,     cbind(       rawdata[last.day >= 365, c("crop", "country")],       period = "mid-season",       duration = rawdata$harvest[last.day >= 365]     )   )   dataset <- rbind(     dataset,     cbind(       rawdata[last.day >= 365, c("crop", "country")],       period = "harvest",       duration = rawdata$harvest.end[last.day >= 365]     )   )   dataset <- rbind(     dataset,     cbind(       rawdata[last.day >= 365, c("crop", "country")],       period = na,       duration = 365 - rowsums(rawdata[last.day >= 365, c("midseason", "harvest", "harvest.end")])     )   )   dataset <- rbind(     dataset,     cbind(       rawdata[last.day >= 365, c("crop", "country")],       period = "sowing",       duration = 365 - rawdata$sowing[last.day >= 365]     )   )   rawdata <- rawdata[last.day < 365, ] }  #mid-season around new year last.day <- rowsums(rawdata[, c("sowing", "midseason", "harvest")]) if(any(last.day >= 365)){   dataset <- rbind(     dataset,     cbind(       rawdata[last.day >= 365, c("crop", "country")],       period = "mid-season",       duration = last.day[last.day >= 365] - 365     )   )   dataset <- rbind(     dataset,     cbind(       rawdata[last.day >= 365, c("crop", "country")],       period = "harvest",       duration = rawdata$harvest.end[last.day >= 365]     )   )   dataset <- rbind(     dataset,     cbind(       rawdata[last.day >= 365, c("crop", "country")],       period = na,       duration = 365 - rowsums(rawdata[last.day >= 365, c("midseason", "harvest", "harvest.end")])     )   )   dataset <- rbind(     dataset,     cbind(       rawdata[last.day >= 365, c("crop", "country")],       period = "sowing",       duration = rawdata$midseason[last.day >= 365]     )   )   dataset <- rbind(     dataset,     cbind(       rawdata[last.day >= 365, c("crop", "country")],       period = "mid-season",       duration = 365 - rowsums(rawdata[last.day >= 365, c("sowing", "midseason")])     )   )   rawdata <- rawdata[last.day < 365, ] }   #harvest around new year last.day <- rowsums(rawdata[, c("sowing", "midseason", "harvest", "harvest.end")]) if(any(last.day >= 365)){   dataset <- rbind(     dataset,     cbind(       rawdata[last.day >= 365, c("crop", "country")],       period = "harvest",       duration = last.day[last.day >= 365] - 365     )   )   dataset <- rbind(     dataset,     cbind(       rawdata[last.day >= 365, c("crop", "country")],       period = na,       duration = 365 - rowsums(rawdata[last.day >= 365, c("midseason", "harvest", "harvest.end")])     )   )   dataset <- rbind(     dataset,     cbind(       rawdata[last.day >= 365, c("crop", "country")],       period = "sowing",       duration = rawdata$midseason[last.day >= 365]     )   )   dataset <- rbind(     dataset,     cbind(       rawdata[last.day >= 365, c("crop", "country")],       period = "mid-season",       duration = rawdata$harvest[last.day >= 365]     )   )   dataset <- rbind(     dataset,     cbind(       rawdata[last.day >= 365, c("crop", "country")],       period = "harvest",       duration = 365 - rowsums(rawdata[last.day >= 365, c("sowing", "midseason", "harvest")])     )   )   rawdata <- rawdata[last.day < 365, ] }   #no crop around new year dataset <- rbind(   dataset,   cbind(     rawdata[, c("crop", "country")],     period = na,     duration = rawdata$sowing   ) ) dataset <- rbind(   dataset,   cbind(     rawdata[, c("crop", "country")],     period = "sowing",     duration = rawdata$midseason   ) ) dataset <- rbind(   dataset,   cbind(     rawdata[, c("crop", "country")],     period = "mid-season",     duration = rawdata$harvest   ) ) dataset <- rbind(   dataset,   cbind(     rawdata[, c("crop", "country")],     period = "harvest",     duration = rawdata$harvest.end   ) ) dataset <- rbind(   dataset,   cbind(     rawdata[, c("crop", "country")],     period = na,     duration = 365 - rowsums(rawdata[, c("sowing", "midseason", "harvest")])   ) )  labels <- c("", "jan.", "feb.", "mar.", "apr.", "may", "jun.", "jul.", "aug.", "sep.", "okt.", "nov.", "dec.") breaks <- cumsum(c(0, 31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31)) ggplot(dataset, aes(x = crop, y = duration, colour = period, fill = period)) + geom_bar(stat = "identity") + facet_wrap(~country) + coord_flip() + scale_fill_manual(values = c("sowing" = "darkgreen", "mid-season" = "grey", "harvest" = "yellow")) + scale_colour_manual(values = c("sowing" = "black", "mid-season" = "black", "harvest" = "black"), guide = "none") + scale_y_continuous("", breaks = breaks, labels = labels, limits = c(0, 365)) + theme_bw() + theme(axis.text.x = element_text(hjust = 1))

Comments

Post a Comment

Popular posts from this blog

javascript - Laravel datatable invalid JSON response -

-
everything work fine on code below fetch user logs: js: <script> $(document).ready(function () { var table =$('#systemlogs').datatable({ responsive: true, processing: true, serverside: true, "language": { "url": "/datatables/media/plug-in/persian.json" }, ajax: '{!! url('/admin/systemlogs/data/systemlogsdatatable') !!}', columns: [ { data: 'name', name: 'name' }, { data: 'message_text', name: 'message_text' }, { data: 'remote_addr', name: 'remote_addr' }, { data: 'log_created_at', name: 'management_logs.created_at' }, { data: 'log_updated_at', name: 'management_logs.updated_at' }, ] }); }); </script> routes: route::get('/systemlogs/data/systemlogsdatatable','systemlogscontrolle...
Read more

java - Exception in thread "main" org.springframework.context.ApplicationContextException: Unable to start embedded container; -

-
i getting error when run spring application: exception in thread "main" org.springframework.context.applicationcontextexception: unable start embedded container; nested exception org.springframework.context.applicationcontextexception: unable start embeddedwebapplicationcontext due missing embeddedservletcontainerfactory bean. @ org.springframework.boot.context.embedded.embeddedwebapplicationcontext.onrefresh(embeddedwebapplicationcontext.java:124) @ org.springframework.context.support.abstractapplicationcontext.refresh(abstractapplicationcontext.java:476) @ org.springframework.boot.context.embedded.embeddedwebapplicationcontext.refresh(embeddedwebapplicationcontext.java:109) @ org.springframework.boot.springapplication.refresh(springapplication.java:691) @ org.springframework.boot.springapplication.run(springapplication.java:320) @ org.springframework.boot.springapplication.run(springapplication.java:952) @ org.springframework.boot.springa...
Read more

sql server 2008 - My Sql Code Get An Error Of Msg 245, Level 16, State 1, Line 1 Conversion failed when converting the varchar value '8:45 AM' to data type int -

-
msg 245, level 16, state 1, line 1 conversion failed when converting varchar value '8:45 am' data type int. script declare @specimenid bigint = 0, @script varchar(max); declare @starttime int = (select convert(int,substring(start_time,0,charindex(':', start_time)) - 12) appt appt._id=112601) --select @starttime set @script = 'select distinct ptlast_name + '', '' + pt.first_name patient_name, dbo.formatdate(pt.dob,''mm/dd/yyyy'') dob, (convert(varchar, (datediff(yy, dob, appt_date) - case when (month(dob) > month(appt_date)) or (month(dob) = month(appt_date) , day(dob) > day(appt_date)) 1 else 0 end)) + ''/'' + convert(varchar, (datediff(m, dateadd(yy, (datediff(yy, dob, appt_date) - case when (month(dob) > month(appt_date)) or (month(dob) = month(appt_date) , day(dob) > day(appt_date)) 1 else 0 end), d...
Read more