[R-sig-Geo] problem with writeRaster and GeoTIFF format

lancelot renaud.lancelot at cirad.fr
Sat May 11 01:45:35 CEST 2013


Thank you very much, Thiago. That did the trick.

All the best,

Renaud

Le 10/05/2013 23:17, Thiago V. dos Santos a écrit :
> Hi,
>
> I had the same issue in the past and I managed to solve it by:
>
> a) removing rgdal
> b) reinstalling gdal and libtiff libraries
> c) reinstalling rgdal
>
> Hope it helps.
> Greetings,
> --
> Thiago V. dos Santos
> PhD student
> Land and Atmospheric Science
> University of Minnesota
> http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/index.htm
> Phone: (612) 323 9898
>
> ----- Original Message -----
> From: lancelot <renaud.lancelot at cirad.fr>
> To: r-sig-geo <r-sig-geo at r-project.org>
> Cc:
> Sent: Friday, May 10, 2013 1:57 PM
> Subject: [R-sig-Geo] problem with writeRaster and GeoTIFF format
>
> Dear all,
>
> I have a raster
>
>  > r
> class : RasterLayer
> dimensions : 1639, 821, 1345619 (nrow, ncol, ncell)
> resolution : 0.00887, 0.00833 (x, y)
> extent : 43.21786, 50.50013, -25.60209, -11.94922 (xmin, xmax,
> ymin, ymax)
> coord. ref. : +proj=longlat +ellps=WGS84
> data source : in memory
> names : layer
> values : 15.3815, 38.32809 (min, max)
>
> I cannot write it to a GeoTIFF file:
>
>  > writeRaster(r, filename = "r.tif")
> Error in .local(.Object, ...) : Unable to create dataset
>
> The same thing happens after converting the raster to
> SpatialGridDataFrame and using writeGDAL
>
>  > r.sgdf <- as(r, Class = "SpatialGridDataFrame")
>  > writeGDAL(dataset = r.sgdf, fname = "r.tif", drivername = "GTiff")
> Error in .local(.Object, ...) : Unable to create dataset
>
> I can write the raster to other formats (e.g. ESRI .hdr labelled). I can
> also reproduce the .tif example in ?writeRaster.
>
> Any hint?
>
> Best wishes,
>
> Renaud
>
>  > sessionInfo()
> R version 3.0.0 (2013-04-03)
> Platform: x86_64-w64-mingw32/x64 (64-bit)
>
> locale:
> [1] LC_COLLATE=French_France.1252 LC_CTYPE=French_France.1252
> LC_MONETARY=French_France.1252 LC_NUMERIC=C
> [5] LC_TIME=French_France.1252
>
> attached base packages:
> [1] grid stats graphics grDevices utils datasets methods
> base
>
> other attached packages:
> [1] maptools_0.8-23 lattice_0.20-15 foreign_0.8-53 mapdata_2.2-2
> maps_2.3-2 raster_2.1-25 rgdal_0.8-9
> [8] sp_1.0-9
>
> loaded via a namespace (and not attached):
> [1] compiler_3.0.0 fortunes_1.5-1 tools_3.0.0
>  >
>
>
> --
> Renaud Lancelot
> Directeur adjoint / Deputy director
> CIRAD, UMR15, Campus International de Baillarguet TA A-DIR / B
> F34398 Montpellier
>
> EDENext Project, coordinator: http://www.edenext.eu/
>
> Tel. +33 4 67 59 37 17 - Fax +33 4 67 59 37 98
> Secr. +33 4 67 59 37 37 - Cell. +33 6 77 52 08 69
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org <mailto:R-sig-Geo at r-project.org>
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo

-- 
Renaud Lancelot
Directeur adjoint / Deputy director
CIRAD, UMR15, Campus International de Baillarguet TA A-DIR / B
F34398 Montpellier

EDENext Project, coordinator: http://www.edenext.eu/

Tel.  +33 4 67 59 37 17  -  Fax  +33 4 67 59 37 98
Secr. +33 4 67 59 37 37  - Cell. +33 6 77 52 08 69



More information about the R-sig-Geo mailing list