On Mar 28, 2010, at 1:11 PM, Muhammad Rahiz wrote:
Yes, it works. Just wondering if the technique can be optimized...
IN WHAT WAY????
--
David.
David Winsemius wrote:
On Mar 28, 2010, at 12:51 PM, Muhammad Rahiz wrote:
Dear Jannis,
Thanks for the tip. It works but I'd like to improve on the way I
did it.
x <- array(1:50,dim=c(10,10)) # data
h1 <- c("ncols 180") # header, line 1
h2 <- c("nrows 290") # header, line 2
h3 <- c("xllcorner -200000") # header, line 3
h4 <- c("yllcorner -200000") # header, line 4
h5 <- c("cellsize 5000") # header, line 5
h6 <- c("NODATA_value NA") # header, line 6
h0 <- rbind(h1,h2,h3,h4,h5,h6)
write(ho,"file.txt")
write
.table(x,"file.txt",row.names=FALSE,col.names=FALSE,append=TRUE)
The header should display as;
ncols 180
nrows 290
xllcorner -200000
yllcorner -200000
cellsize 5000
NODATA_value NA
Executing that code appears to deliver what you requested. So
rather than asking us to read you mind, can you explain what is
missing?
ncols 180
nrows 290
xllcorner -200000
yllcorner -200000
cellsize 5000
NODATA_value NA
1 11 21 31 41 1 11 21 31 41
2 12 22 32 42 2 12 22 32 42
3 13 23 33 43 3 13 23 33 43
4 14 24 34 44 4 14 24 34 44
5 15 25 35 45 5 15 25 35 45
6 16 26 36 46 6 16 26 36 46
7 17 27 37 47 7 17 27 37 47
8 18 28 38 48 8 18 28 38 48
9 19 29 39 49 9 19 29 39 49
10 20 30 40 50 10 20 30 40 50
David Winsemius, MD
West Hartford, CT
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.