[R] insert/copy skipped lines to output

[Muhammad Rahiz]

something like this...

x1 <- as.matrix(read.table("input.txt"))[1:6,] # reads only header
x2 < -array(1:6,dim=c(2,3))

write(x,"output.txt")
write.table(x3,"output.txt",row.names=FALSE,col.names=FALSE,append=TRUE)

But the results print;

ncols
nrows
xllcorner
yllcorner
cellsize
NODATA_value
    180.00
    290.00
-200000.00
-200000.00
   5000.00
  -9999.00
1 3 5
2 4 6




Muhammad Rahiz  |  Doctoral Student in Regional Climate Modeling
Climate Research Laboratory, School of Geography & the Environment	
Oxford University Centre for the Environment, University of Oxford
South Parks Road, Oxford, OX1 3QY, United Kingdom
Tel: +44 (0)1865-285194	 Mobile: +44 (0)7854-625974
Email: muhammad.rahiz at ouce.ox.ac.uk



David Winsemius wrote:
> 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
>
>