---------- Forwarded message ---------- From: Agustin Lobo <Agustin.Lobo at ija.csic.es> Date: 2010/12/16 Subject: Re: Problem with writing parameters in overlay() To: "Robert J. Hijmans" <r.hijmans at gmail.com>, Roger.Bivand at nhh.no, alobolistas at gmail.com Robert, The weird problem at writing raster layers in some cases is also solved once I remove gdal1.8dev:
require(rgdal)
Loading required package: rgdal Geospatial Data Abstraction Library extensions to R successfully loaded Loaded GDAL runtime: GDAL 1.7.3, released 2010/11/10 Path to GDAL shared files: /usr/share/gdal17 Loaded PROJ.4 runtime: Rel. 4.7.1, 23 September 2009 Path to PROJ.4 shared files: (autodetected) ?require(raster)
bN3v2 <- overlay(bN3, bSM3, fun=fun1,filename="bN3v2",overwrite=TRUE,format="EHdr", NAflag=0)
works fine now. Considering the bug fixed in gdal1.8 that I forwarded you, the problem should also be solved even using the newest gdal1.8, but I rather do not test until gdal1.8 is available in ubuntugis. You might want to keep an eye on this, though. Agus
On 11/28/2010 01:33 AM, Robert J. Hijmans wrote:
Is it only with overlay, or do other functions have the same problem? ?a<- calc(r2, fun=function(x)x*2, filename="test2", format='GTiff',overwrite=T) ?a<- aggregate(r2, 10, filename="test2", format='GTiff',overwrite=T) Robert On Sat, Nov 27, 2010 at 12:34 PM, Robert J. Hijmans<r.hijmans at gmail.com> ?wrote:
Hi Agus, This is very strange. Somehow the filetype and filename get mixed up. I cannot reproduce this. Does the below all work for you? (it does for me) library(raster) r<- raster(ncol=10, nrow=10) r<- writeStart(r, filename="test", format="GTiff", overwrite=TRUE) filename(r) r<- writeValues(r, 1:100, 1) filename(r) r<- writeStop(r) filename(r) r<- raster(ncol=10, nrow=10) r[]<- 1:100 r2<- r a<- overlay(r, r2, fun=sum, filename="test.tif", overwrite=TRUE, NAflag=0) b<- overlay(r, r2, fun=sum, filename="test.img",overwrite=TRUE, NAflag=0) setOptions(todisk=TRUE) a<- overlay(r, r2, fun=sum, filename="test1", format='GTiff', overwrite=TRUE, NAflag=0) r<- writeRaster(r, filename='r1', overwrite=TRUE) r2<- writeRaster(r2, filename='r1', overwrite=TRUE) a<- overlay(r, r2, fun=sum, filename="test2", format='GTiff', overwrite=TRUE, NAflag=0) setOptions(todisk=FALSE) Thanks for your help, Robert On Sat, Nov 27, 2010 at 3:37 AM, Agustin Lobo<Agustin.Lobo at ija.csic.es> ?wrote:
Robert, There is still this problem: These work:
bN3v2<- overlay(bN3, bSM3, fun=fun1,filename="bN3v2",NAflag=0) bN3v2<- overlay(bN3, bSM3, fun=fun1,filename="bN3v2",overwrite=TRUE,NAflag=0)
These do not work (is not a parsing problem?) (Not critical for me, just trying to help out testing):
bN3v2<- overlay(bN3, bSM3, fun=fun1,filename="bN3v2",overwrite=TRUE, format="EHdr", NAflag=0)
Error in .rasterObjectFromFile(x, band = band, objecttype = "RasterLayer", ?: ?file: /home/alobo/EHdr does not exist Calls: overlay ... raster -> ?raster -> ?.local -> ?.rasterObjectFromFile
bN3v2<- overlay(bN3, bSM3, fun=fun1,filename="bN3v2", format="EHdr", NAflag=0)
Error in .rasterObjectFromFile(x, band = band, objecttype = "RasterLayer", ?: ?file: /home/alobo/EHdr does not exist Calls: overlay ... raster -> ?raster -> ?.local -> ?.rasterObjectFromFile
bN3v2<- overlay(bN3, bSM3, fun=fun1,filename="bN3v2",overwrite=TRUE, format="EHdr", NAflag=0)
Error in .rasterObjectFromFile(x, band = band, objecttype = "RasterLayer", ?: ?file: /home/alobo/EHdr does not exist Calls: overlay ... raster -> ?raster -> ?.local -> ?.rasterObjectFromFile
bN3v2<- overlay(bN3, bSM3, fun=fun1,filename="bN3v2", NAflag=0, format="GTiff")
Error in .rasterObjectFromFile(x, band = band, objecttype = "RasterLayer", ?: ?file: /home/alobo/GTiff does not exist Calls: overlay ... raster -> ?raster -> ?.local -> ?.rasterObjectFromFile
sessionInfo()
R version 2.12.0 (2010-10-15) Platform: i486-pc-linux-gnu (32-bit) locale: ?[1] LC_CTYPE=en_US.UTF-8 ? ? ? ? ?LC_NUMERIC=C LC_TIME=en_US.UTF-8 ?[4] LC_COLLATE=en_US.UTF-8 ? ? ? ?LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8 ?[7] LC_PAPER=en_US.UTF-8 ? ? ? ? ?LC_NAME=en_US.UTF-8 LC_ADDRESS=en_US.UTF-8 [10] LC_TELEPHONE=en_US.UTF-8 ? ? ?LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=en_US.UTF-8 attached base packages: [1] stats ? ? graphics ?grDevices utils ? ? datasets ?methods ? base other attached packages: [1] rgdal_0.6-27 ?raster_1.6-22 sp_0.9-72 ? ? rkward_0.5.4 loaded via a namespace (and not attached): [1] grid_2.12.0 ? ? lattice_0.19-13 tools_2.12.0 Agus -- Dr. Agustin Lobo Institut de Ciencies de la Terra "Jaume Almera" (CSIC) LLuis Sole Sabaris s/n 08028 Barcelona Spain Tel. 34 934095410 Fax. 34 934110012 email: Agustin.Lobo at ija.csic.es http://www.ija.csic.es/gt/obster