Yes, you are right. Thanks for the correction.
Besides, I have found a bug when add.contour=FALSE. Tomorrow I will
upload to R-Forge a new version of rasterVis to fix this bug and to
improve the labelling of the x-axis.
Best,
Oscar.
El Wed, 17 Aug 2011 15:10:25 -0700
Robert J. Hijmans
Hi Michael,
Thanks for the reply.
The data file can be rechieved from the following link:
ftp://ftp.wiley.com/public/sci_tech_med/spatio_temporal_data
SST011970_032003.datftp://ftp.wiley.com/public/sci_tech_med/spatio_temporal_data/SST011970_032003.dat
On Wed, 17 Aug 2011, Jianyun Wu wrote:
Hi Michael,
Thanks for the reply.
The data file can be rechieved from the following link:
ftp://ftp.wiley.com/public/sci_tech_med/spatio_temporal_data
From Edzer Pebesma's notes in demo/CressieWikle.R in spacetime for the
equator:
library(spacetime)
Here's an update using the two files you need at a minimum:
## read the entire matrix file
x - as.matrix(read.table(SST011970_032003.dat))
## read the coordinates
lonlat - read.table(SSTlonlat.dat)
names(lonlat) - c(Longitude, Latitude)
## choose just one lat
x1 - x[lonlat$Latitude == -29, ]
Hi Roger,
It seems really helpful, I will go through Vignettes of spacetime package
to see how it goes.
Thanks and Regards
Fred
On Wed, Aug 17, 2011 at 5:13 PM, Roger Bivand roger.biv...@nhh.no wrote:
On Wed, 17 Aug 2011, Jianyun Wu wrote:
Hi Michael,
Thanks for the reply.
The data file
Fred,
Perhaps you can use the Hovmoller function in the rasterVis package.
Robert
On Wed, Aug 17, 2011 at 12:20 AM, Jianyun Wu jianyun.fred...@gmail.comwrote:
Hi Roger,
It seems really helpful, I will go through Vignettes of spacetime package
to see how it goes.
Thanks and Regards
Fred
Hi,
You will find a hovmoller function in the rasterVis package:
http://rastervis.r-forge.r-project.org/
With your data:
library(raster)
library(rasterVis)
dat - read.table('SST011970_032003.dat', header=FALSE)
loc - read.table('SSTlonlat.dat', header=FALSE)
rasterList - lapply(dat,
Dear Experts,
I am trying to plot the Hovmueller Diagram (example below), which is time
vs. longitude in R.
