small improvement for saveDataset() in rgdal
Seb, Thanks for a helpful example and patch - committed in SVN revision 411. Best wishes, Roger
On Fri, 14 Dec 2012, Sebastian P. Luque wrote:
Hi,
Some time ago I spent several hours trying to figure out why the
following was failing?:
library(rgdal); library(adehabitatHR)
data(puechabonsp)
locs <- puechabonsp$relocs
## Estimation of UD for four animals
(ud <- kernelUD(locs[, 1]))
locs.gtiff <- create2GDAL(ud[[1]],
drivername="GTiff")
locs.aaigrid <- copyDataset(locs.gtiff,
driver="AAIGrid")
saveDataset(locs.aaigrid, "~/tmp/test.asc")
which prints this:
Error in .local(.Object, ...) : Dataset copy failed
finding that this resulted from failure to expand the 'filename'
argument in saveDataset(). Would it make sense to patch this function
with:
---<--------------------cut here---------------start------------------->---
*** /home/sluque/Scripts/R/src/rgdal/R/gdal.R 2012-12-14 15:05:25.000000000 -0600
--- /home/sluque/Scripts/R/src/gdal.R 2012-12-14 15:05:11.000000000 -0600
***************
*** 217,223 ****
new.obj <- new(new.class,
handle = .Call('RGDAL_CopyDataset',
dataset, getDriver(dataset),
! FALSE, options, filename, PACKAGE="rgdal"))
if (returnNewObj) return(new.obj)
invisible(GDAL.close(new.obj))
--- 217,223 ----
new.obj <- new(new.class,
handle = .Call('RGDAL_CopyDataset',
dataset, getDriver(dataset),
! FALSE, options, normalizePath(filename), PACKAGE="rgdal"))
if (returnNewObj) return(new.obj)
invisible(GDAL.close(new.obj))
---<--------------------cut here---------------end--------------------->---
? It tests fine on my system:
R> sessionInfo()
R version 2.15.2 (2012-10-26)
Platform: x86_64-pc-linux-gnu (64-bit)
locale:
[1] LC_CTYPE=en_CA.UTF-8 LC_NUMERIC=C LC_TIME=en_CA.UTF-8 LC_COLLATE=en_CA.UTF-8 LC_MONETARY=en_CA.UTF-8 LC_MESSAGES=en_CA.UTF-8
[7] LC_PAPER=C LC_NAME=C LC_ADDRESS=C LC_TELEPHONE=C LC_MEASUREMENT=en_CA.UTF-8 LC_IDENTIFICATION=C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] adehabitatHR_0.4.4 adehabitatLT_0.3.7 CircStats_0.2-4 boot_1.3-7 MASS_7.3-22 adehabitatMA_0.3.4 ade4_1.5-1 deldir_0.0-21 rgdal_0.7-25
[10] sp_1.0-2 slmisc_0.9.3 latticeExtra_0.6-24 RColorBrewer_1.0-5 lattice_0.20-10
loaded via a namespace (and not attached):
[1] compiler_2.15.2 grid_2.15.2 tools_2.15.2
Cheers,
--
Seb
? https://stat.ethz.ch/pipermail/r-sig-ecology/2012-March/002895.html
_______________________________________________ R-sig-Geo mailing list R-sig-Geo at r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Roger Bivand Department of Economics, NHH Norwegian School of Economics, Helleveien 30, N-5045 Bergen, Norway. voice: +47 55 95 93 55; fax +47 55 95 95 43 e-mail: Roger.Bivand at nhh.no