Hi all, Jonathan very kindly helped me with dramatically improve my raster algebra performance (an earlier post).
I have another bottle-neck now, for which I was hoping that someone may have a solution. The stack() function takes a long time, so I moved to the quickStack function suggested on http://stackoverflow.com/questions/9929324/increase-performance-speed. Changing the raster stack to a matrix using the values() function is now the bottleneck. Does anyone know if it is possible to read a series of rasters directly into a matrix and thereby drop the values() step? I have been trying to take a look at the code for the stack() function, but can't seem to see it. How does one view the function code - or is it written in C? And how does one figure that out? An example of the code would be: in <- stack(list.files("E:\Temp")) v <- values(in) Thanks very much for any help in advance. Cheers, Alastair -- Dr. Alastair J. Potts Claude Leon Postdoctoral Fellow Botany Department Nelson Mandela Metropolitan University Cell #: 082 491-7275 Office #: 041 504-2398 (Tues, Thurs, Fri) Fax #: 086 686-4251 [[alternative HTML version deleted]] _______________________________________________ R-sig-Geo mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-sig-geo
