[R-sig-Geo] writeRaster() writing float despite stating INT2U

[Dutrieux, Loic]

Hi Augustin,

Could it be because of the typo in the code ("dataype")?

Cheers,
Loïc

-----Original Message-----
From: r-sig-geo-bounces at r-project.org [mailto:r-sig-geo-bounces at r-project.org] On Behalf Of Agustin Lobo
Sent: jeudi 11 juillet 2013 15:19
To: r-sig-geo
Subject: [R-sig-Geo] writeRaster() writing float despite stating INT2U

I'm having this strange problem: I write a raster as INT2U but get a geotif with Float64. Is perhaps a gdal problem?


> r <- raster(matrix(1:100,nrow=10,ncol=10))
> writeRaster(r, 
> file="r.tif",format="GTiff",dataype="INT2U",overwrite=TRUE,NAflag=0)
rgdal: version: 0.8-10, (SVN revision 478) Geospatial Data Abstraction Library extensions to R successfully loaded Loaded GDAL runtime: GDAL 1.10.0, released 2013/04/24 Path to GDAL shared files: /usr/share/gdal/1.10 Loaded PROJ.4 runtime: Rel. 4.8.0, 6 March 2012, [PJ_VERSION: 480] Path to PROJ.4 shared files: (autodetected)
class       : RasterLayer
dimensions  : 10, 10, 100  (nrow, ncol, ncell) resolution  : 0.1, 0.1  (x, y)
extent      : 0, 1, 0, 1  (xmin, xmax, ymin, ymax)
coord. ref. : NA
data source : /media/alobo/Iomega_HDD/PIVOT/RPivot/r.tif
names       : r
values      : 1, 100  (min, max)

> GDALinfo("r.tif")
rows        10
columns     10
bands       1
lower left origin.x        0
lower left origin.y        0
res.x       0.1
res.y       0.1
ysign       -1
oblique.x   0
oblique.y   0
driver      GTiff
projection  NA
file        r.tif
apparent band summary:
   GDType hasNoDataValue NoDataValue blockSize1 blockSize2
1 Float64           TRUE           0         10         10
apparent band statistics:
  Bmin Bmax Bmean      Bsd
1    1  100  50.5 29.01149

> sessionInfo()
R version 2.15.2 (2012-10-26)
Platform: x86_64-pc-linux-gnu (64-bit)

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C
LC_TIME=es_ES.UTF-8        LC_COLLATE=en_US.UTF-8
LC_MONETARY=es_ES.UTF-8
 [6] LC_MESSAGES=en_US.UTF-8    LC_PAPER=C                 LC_NAME=C
               LC_ADDRESS=C               LC_TELEPHONE=C
[11] LC_MEASUREMENT=es_ES.UTF-8 LC_IDENTIFICATION=C

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base

other attached packages:
[1] rgdal_0.8-10  raster_2.1-37 sp_1.0-9

loaded via a namespace (and not attached):
[1] grid_2.15.2     lattice_0.20-13 tools_2.15.2

Agus

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo at r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo