Hi, I've been working on a set of functions for a bit more than a month. This morning I updated all my packages (but not R) and now, when I call:
outfile <- "b4.sinks.sgrd" writeRaster(myraster, outfile, overwrite = TRUE, datatype = "FLT4S", NAflag = -9999) or infile <- "b4.sinks.sdat" r <- raster(infile) from within one of my functions, I get this error: Error in if (is.null(projs) | is.na(projs) | nchar(projs) < 3) { : argument is of length zero In addition: Warning message: In is.na(projs) : is.na() applied to non-(list or vector) of type 'NULL' My function that calls writeRaster() and raster() hasn't changed in a while. It was working fine just yesterday! From the error message, it looks like .makeProj (in makeProjString.R) calls .newCRS with an unknown projection, although I never changed my raster's projections since I started working on this project. Here's the output of proj4string(myraster): "+proj=utm +zone=18 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0" Here's the output of sessionInfo(): R version 3.0.2 (2013-09-25) Platform: x86_64-w64-mingw32/x64 (64-bit) locale: [1] LC_COLLATE=French_Canada.1252 LC_CTYPE=French_Canada.1252 LC_MONETARY=French_Canada.1252 [4] LC_NUMERIC=C LC_TIME=French_Canada.1252 attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] snow_0.3-13 raster_2.1-66 devtools_1.3 stringr_0.6.2 RSAGA_0.93-6 plyr_1.8 [7] shapefiles_0.7 foreign_0.8-55 gstat_1.0-17 rgdal_0.8-11 data.table_1.8.10 sp_1.0-14 loaded via a namespace (and not attached): [1] digest_0.6.3 evaluate_0.5.1 grid_3.0.2 httr_0.2 intervals_0.14.0 lattice_0.20-23 [7] memoise_0.1 parallel_3.0.2 RCurl_1.95-4.1 spacetime_1.0-9 tools_3.0.2 whisker_0.3-2 [13] xts_0.9-7 zoo_1.7-10 I'm getting frustrated as I'm about to finish this program I've been working on for a while and now I'm stuck!! Thanks for your help, Guillaume _______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo