[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