[R-sig-Geo] parallelize projectRaster()
Yerguner
yasebaytok at gmail.com
Wed Jun 18 17:44:09 CEST 2014
Hi Jonathan,
Thanks a lot for your reply. I've already tried your gdalUtils() package but
it creates a empty .tif file (all zero values) and continues to run several
(>14) hrs. May be I'm not using the function correctly:
library(raster)
rasterOptions(tmpdir="/pegasus/yb3/modvar2000/")
library(ncdf4)
library(rgdal)
library(gdalUtils)
gdal_setInstallation()
ifile2 <- "cold00_s2.nc"
newproj <- "+proj=laea +lat_0=45 +lon_0=-100 +x_0=0 +y_0=0 +a=6370997
+b=6370997 +units=m +no_defs"
gdalwarp(srcfile= ifile2, dstfile= "ohey_file.tif", t_srs= newproj,
output_Raster=TRUE, multi=TRUE, tr=c(231.6564, 231.6563),r="near")
class : RasterBrick
dimensions : 2049, 4417, 9050433, 365 (nrow, ncol, ncell, nlayers)
resolution : 231.6564, 231.6563 (x, y)
extent : -1075057, -51831.17, 241865.9, 716529.7 (xmin, xmax, ymin,
ymax)
coord. ref. : +proj=laea +lat_0=45 +lon_0=-100 +x_0=0 +y_0=0 +a=6370997
+b=6370997 +units=m +no_defs
data source : /pegasus/yb3/daymet2000/ohey_file.tif
names : ohey_file.1, ohey_file.2, ohey_file.3, ohey_file.4,
ohey_file.5, ohey_file.6, ohey_file.7, ohey_file.8, ohey_file.9,
ohey_file.10, ohey_file.11, ohey_file.12, ohey_file.13, ohey_file.14,
ohey_file.15, ...
Thanks,
Yasemin
--
View this message in context: http://r-sig-geo.2731867.n2.nabble.com/parallelize-projectRaster-tp7586598p7586604.html
Sent from the R-sig-geo mailing list archive at Nabble.com.
More information about the R-sig-Geo
mailing list