updated write.tgrass, after change in raster

R/tgrass.R

@@ -25,18 +25,19 @@ read.tgrass = function(fname, localName = TRUE, useTempDir = TRUE, isGeoTiff = T
 
 write.tgrass = function(obj, fname, ...) {
 	# require(raster)
+	if (!requireNamespace("raster", quietly = TRUE))
+		stop("package raster required for write.tgrass")
 	if (is(obj, "STFDF")) {
 		end.time = obj@endTime
-		obj = as(obj, "RasterStack")
+		#obj = as(obj, "RasterStack")
+		obj = as(obj, "RasterBrick")
 	} else
 		end.time = NULL
 	n = paste(names(obj), ".tif", sep = "")
 	olddir = getwd()
 	dir = tempdir() # gets the same as read.tgrass got!
 	setwd(dir)
 	# write .tifs:
-	if (!requireNamespace("raster", quietly = TRUE))
-		stop("package raster required for write.tgrass")
 	for (i in 1:(raster::nlayers(obj)))
 		raster::writeRaster(raster::raster(obj, layer=i), n[i], ...)
 	# write proj.txt:

TODO

@@ -14,3 +14,7 @@ subsetting STF, STS, STI with [
 kriging without STFDF's, kriging with STIDF's (I think this is supported now)
 
 single location prediction.
+
+Sat Oct 17 22:40:29 CEST 2015
+Meng: coerce STIDF -> STFDF: do NOT set endTime to full interval?
+  alternative: (a) zero interval width, (b) mean of endTime of each time slice?