##ps_layer=vector ##plot_mean=string TRUE ##plot_mean_only=string FALSE ##date1=string 01/01/2018 ##date2=string 01/2/2021 ##title=string Landslide ##dateformat=string %d/%m/%Y ##begin_timeseries=number 6 ##xlabel=string Date ##ylabel=string Deformation (mm) ##outputfile=string /home/calderone/backup_script_R_2022/output/defo.pdf ##output_width=number 20 ##output_height=number 5 ##ymax=number 20 ##ymin=number -20 #output_plots_to_html ##showplots source("/home/calderone/backup_script_R_2022/r_functions/plotting_functions_backup.R") library(ggplot2) library(reshape) ps_layer_df=as.data.frame(ps_layer) ps_layer_df$geometry <- NULL write.csv(ps_layer_df,"/home/calderone/backup_script_R_2022/ex//df.csv") begin_timeseries <- begin_timeseries*-1 ps_layer_df <- ps_layer_df[,c(-1:begin_timeseries)] #print(ps_layer_df) if (length(ps_layer_df$name)==0){ ps_layer_df$name <- "landslide" } if (length(ps_layer_df$plot)==0){ ps_layer_df$plot <- 1 } ps_layer_df$ps <- rownames(ps_layer_df) #print ("ps_layer_df") #write.csv(ps_layer_df,"/home/calderone/backup_script_R_2022/ex//df_rstudio.csv") ps_layer_plot_frame <- melt(ps_layer_df,id.vars=c("name","plot","ps")) #print (ps_layer_plot_frame) ps_layer_plot_frame$date <- as.Date(ps_layer_plot_frame$variable,format="d%Y%m%d") plot_frame_sum<-summarySE(ps_layer_plot_frame,measurevar="value",groupvars=c("date","name"),na.rm=T) if(plot_mean=="TRUE"){ defo<-ggplot()+ geom_line(data=ps_layer_plot_frame, aes(x = as.Date(date), y=value,colour=ps))+ geom_point(data=ps_layer_plot_frame, aes(x = as.Date(date), y=value,colour=ps))+ scale_color_grey(start=0.6,end=0.6,guide = 'none')+ geom_errorbar(data=plot_frame_sum,aes(x = as.Date(date),ymin=value-ci, ymax=value+ci,linetype=name), size=1,width=10)+ geom_line(data=plot_frame_sum,aes(x = as.Date(date), value,linetype=name), size=2)+ ylab(ylabel) + theme_bw(base_size = 20) + theme(legend.justification = c(0, 1), legend.position = c(0, 1))+ ggtitle(title)+ ylim(ymin,ymax)+ xlab(xlabel)+ scale_x_date(breaks = "6 month", minor_breaks = "2 month",limits=as.Date(c(date1, date2),format=dateformat))+ theme(axis.text.x = element_text(angle = 0, hjust = 1))+ theme(panel.grid.major = element_line(colour = "darkgrey"),panel.grid.minor = element_line(colour = "darkgrey")) } if(plot_mean_only=="TRUE"){ defo<-ggplot()+ geom_errorbar(data=plot_frame_sum,aes(x = as.Date(date),ymin=value-ci, ymax=value+ci,colour=name), size=1,width=10)+ geom_line(data=plot_frame_sum,aes(x = as.Date(date), value,colour=name), size=2)+ scale_color_grey(start=0.3,end=0.6,guide = 'none')+ ylab(ylabel) + theme_bw(base_size = 20) + theme(legend.justification = c(0, 1), legend.position = c(0, 1))+ ggtitle(title)+ xlab(xlabel)+ylim(ymin,ymax)+ scale_x_date(breaks = "6 month", minor_breaks = "2 month",limits=as.Date(c(date1, date2),format=dateformat))+ theme(axis.text.x = element_text(angle = 0, hjust = 1))+ theme(panel.grid.major = element_line(colour = "darkgrey"),panel.grid.minor = element_line(colour = "darkgrey")) } defo pdf(outputfile,width=output_width, height=output_height) defo dev.off()