Commit: b5df8df897105669b6e64a1d2eac721facf8bdfb - myprojects-hg-reborn (hg)

Commit MetaInfo

• Revisión: b5df8df897105669b6e64a1d2eac721facf8bdfb (tree)
• Tiempo: 2007-07-29 19:22:30
• Autor: iselllo
• Commiter: iselllo

Log Message

I added the code 3D_R_plotting.R which correctly uses the lattice
package for 3D plots in R and properly save the figure in a file.

Cambiar Resumen

• add: R-codes/3D_R_plotting.R

Diferencia incremental

diff -r 8f44f41ea645 -r b5df8df89710 R-codes/3D_R_plotting.R
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/R-codes/3D_R_plotting.R	Sun Jul 29 10:22:30 2007 +0000
@@ -0,0 +1,73 @@
+rm(list=ls())
+library(lattice)
+data_150<-matrix(scan("y_complete00100"),ncol=100,byrow=T)
+print("Ok reading the data")
+data_150<-as.matrix(data_150)
+
+
+D_mean<-read.table("D_mean_seq",header=FALSE)
+D_mean<-as.matrix(D_mean)
+
+ time<-read.table("time_grid", header=FALSE)
+time<-as.matrix(time)
+
+t_len<-length(time)
+t_seq<-seq(1:length(time))
+
+
+g <- expand.grid(x = D_mean, y = sort(time,decreasing=FALSE))
+
+lentot<-length(D_mean)*length(time)
+mydens<-seq(1,lentot)
+# dim(mydens)<-c(lentot,1)
+for (i in t_seq)
+{j<-(i-1)*length(D_mean)+1
+mydens[j:(j+length(D_mean)-1)]<-data_150[i, ]
+}
+
+g$z<-mydens
+
+# myfig_old<-wireframe(z ~ x * y, g, drape = TRUE,shade=TRUE,
+#  scales = list(arrows = FALSE),pretty=FALSE, aspect = c(1,1), colorkey = TRUE
+#  ,zoom=1, main=expression("Evolution Size Distribution"), zlab = list(expression("density"),rot =90),distance=0.0
+#  ,xlab = list(expression("log(D_p [nm])"),rot =0),
+#  ylab = list(expression("Time [s]"),rot =0),
+#  perspective=TRUE,screen = list(z = 200, x = -35,y= 00),cex.axis=2,
+#  )
+
+# trellis.device(pdf,
+#  file="evolution-3D-size-distribution_old.pdf")
+#  trellis.par.set("axis.line",
+#  list(col="transparent"))
+#  print(myfig_old)
+#  dev.off()
+
+# , at = c(5,10,30, 50, 100, 250, 500)
+
+
+
+myfig<-wireframe(z ~ x * y, g,font = 3, tck = 1,screen = list(z = 180, x = -80, y=0),colorkey = FALSE,  xlab = expression(paste(D[agg]," [nm]")),
+ ylab = expression(paste(tau ," [s]")),
+ zlab = list(expression(paste(dN/dlogD[agg]," ["*cm^-3*"]")),rot=90),
+scales = list(arrows = FALSE,cex= 0.7, col = "black",font = 1, tck = 0.8,x=list(log = TRUE, at = c(5,10,30, 50, 100, 250, 500))), shade=TRUE,
+ light.source= c(10,10,10), shade.colors = function(irr, ref,
+ height, w = 0.4)
+ grey(w * irr + (1 - w) *
+ (1 - (1 - ref)^0.4)),aspect = c(1, 0.65)) 
+
+
+
+# trellis.device(postscript,
+#  file="evolution-3D-size-distribution.ps",onefile = FALSE, paper = "special",
+#  horizontal = FALSE, width = 8, height = 6)
+ 
+trellis.device(pdf,
+  file="evolution-3D-size-distribution.pdf",onefile = FALSE, paper = "special"
+, width = 8, height = 6)
+
+trellis.par.set("axis.line",
+ list(col="transparent"))
+ print(myfig)
+ dev.off()
+
+print("So far so good")
\ No newline at end of file