Displaying 3 results from an estimated 3 matches for "rsts1".
Did you mean:
rsts
2018 Mar 06
1
raster time series statistics
...39;-',unlist(dt['month']),"-01")
strptime(timelst,format="%Y-%m-%d")->t1
rts(s,time=as.yearmon(t1))->rsts
subset(rsts,'2017')->r2017
class(r2017 at time)
class(rsts at time)
apply.monthly(rsts,mean) # this creates error
rts(s,time=as.Date(t1))->rsts1
apply.monthly(rsts1,mean) # this creates output
-----Original Message-----
From: Jim Lemon [mailto:drjimlemon at gmail.com]
Sent: Tuesday, 6 March 2018 11:40 AM
To: Herr, Alexander (L&W, Black Mountain) <Alexander.Herr at csiro.au>
Cc: David Winsemius <dwinsemius at comcast.net>;...
2018 Mar 06
0
raster time series statistics
I can't test that at the moment as I don't have the libraries. Perhaps later.
Jim
On Tue, Mar 6, 2018 at 11:36 AM, <Alexander.Herr at csiro.au> wrote:
> Last line in the following (updated) code produces the error
> require(raster)
> require(rts)
> require(stringr)
> r <- raster(ncol=100, nrow=100)
> values(r) <- runif(ncell(r))
> stack(r)->s
>
2018 Mar 06
2
raster time series statistics
Last line in the following (updated) code produces the error
require(raster)
require(rts)
require(stringr)
r <- raster(ncol=100, nrow=100)
values(r) <- runif(ncell(r))
stack(r)->s
r->rs
for(i in 1:23){
rs[]<-r[]*i
addLayer(s,rs)->s
print(nlayers(s))
}
dt<-list(ID=seq(1:24),month=rep(formatC(1:12,flag=0,width=2),2),
year=sort(rep(2016:2017,12)))