Skip to content
Prev 44931 / 398528 Next

adding header info to write.table

On Thu, 26 Feb 2004, femke wrote:

            
zz <- file("myraster", "w")
cat("ncols         532\nnrows         999\n", file=zz)
cat("xllcorner     510465\nyllcorner     4766375\n", file=zz)
cat("cellsize      30\nNODATA_value  -9999\n", file=zz)
write.table(data, file=zz, append=TRUE, col.names=FALSE, row.names=FALSE)
close(zz)

using connections (see help(connections)) and the append= argument in 
write.table. Is everybody writing ArcGIS ASCII rasters recently?