[R-sig-Geo] error exporting ESRI raster to raster native format

Herry Alexander.Herr at csiro.au
Fri Jul 25 02:43:12 CEST 2014


hi list,

after loading a large ESRI raster file with a raster attribute table I
cannot process it. A simple writeRaster results in an unusable file
importing without the RAT.
>rster<-raster('GDA94albers/clum50m0314', RAT=FALSE)

image(rster) works fine
>writeRaster(rster,'luse4BA.grd', progress='text')->rs
image(rs) results in 

"Error in setValues(outras, m) : 
  length(values) is not equal to ncell(x), or to 1
In addition: Warning messages:
1: In .readCells(x, cells, 1) : NAs introduced by coercion
2: In cbind(cells[, 1], vals[, 2:ncol(vals)]) :
  number of rows of result is not a multiple of vector length (arg 1)"

I get the same result importing the raster and doing a deratify...
info on the raster:
>rster
class       : RasterLayer 
dimensions  : 76740, 80200, 6154548000  (nrow, ncol, ncell)
resolution  : 50, 50  (x, y)
extent      : -1888000, 2122000, -4847000, -1010000  (xmin, xmax, ymin,
ymax)
coord. ref. : +proj=aea +lat_1=-18 +lat_2=-36 +lat_0=0 +lon_0=132 +x_0=0
+y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs 
data source : /data/her134/BA/GDA94albers/clum50m0314 
names       : clum50m0314 
values      : 100, 663  (min, max)

I am suspecting this is a problem with importing the ESRI raster...

Any suggestions?

Thanks and cheers
herry






--
View this message in context: http://r-sig-geo.2731867.n2.nabble.com/error-exporting-ESRI-raster-to-raster-native-format-tp7586795.html
Sent from the R-sig-geo mailing list archive at Nabble.com.



More information about the R-sig-Geo mailing list