Dear R help contributors,
I'd like to plot ground temperature with time on X-axis and depth on Y-axis
on this datasets ( http://r.789695.n4.nabble.com/file/n3301033/NEdaily.csv
NEdaily.csv ), and to do so I use the following commands:
library(RSEIS)
xNE <- seq(1, as.numeric(as.Date(max(NEdaily[[1]])) -
as.Date(min(NEdaily[[1]]))), 1)
yNE <- rev(c(-0.3, -0.5, -0.7, -0.9, -1.1, -1.4, -1.7, -2, -2.5, -3, -4,
-5, -7, -9, -10))
zNE <-
mirror.matrix(as.matrix(NEdaily[1:(nrow(NEdaily)-1),2:length(NEdaily)]))
filled.contour(xNE,yNE,zNE
, col = myPal(20)
, zlim = c(-20,20)
, ylab = "Depth [m]",
, xlab = paste("Days since ", as.Date(min(NEdaily[[1]]), format
="%d.%m.%Y"))
)
contour(xNE,yNE,zNE, lty = 3, add = T)
contour(xNE,yNE,zNE, nlevels = 1, level = 0, add = T, lwd = 1.5)
I get this graph ( http://r.789695.n4.nabble.com/file/n3301033/NEdaily.png
NEdaily.png ) and don't understand why filled.contour and contour plots are
no set on the same dimensions and why they don't exactly overlay. Does
anyone have an idea and a solution ?!
Thanks in advance,
Xavier
--
View this message in context:
http://r.789695.n4.nabble.com/Using-filled-contour-and-contour-functions-together-tp3301033p3301033.html
Sent from the R help mailing list archive at Nabble.com.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
