At least for the first part, if you already have a model (and a function),
you could calculate the probability of transition between classes using
raster::calc function.


s <- stack(r, r*2, sqrt(r))
# return a RasterLayer
rs1 <- calc(s, sum)


Cheers,

Roman



On Thu, Apr 7, 2011 at 12:07 AM, Jonathan Greenberg
<greenberg@ucdavis.edu>wrote:

> Folks:
>
> Does anyone know of a package (or a suggestion on how to implement) to
> calculate, for two classified raster images of the same location but
> different times, the relative probability of transitioning from one class
> to
> the other?  Additionally, once this is figured out, how to apply this
> transition matrix the one of the rasters to predict future per-class
> probabilities?  Thanks!  Ideally this should support arbitrarily large
> files.
>
> --j
>
> --
> Jonathan A. Greenberg, PhD
> Assistant Project Scientist
> Center for Spatial Technologies and Remote Sensing (CSTARS)
> Department of Land, Air and Water Resources
> University of California, Davis
> One Shields Avenue
> Davis, CA 95616
> Phone: 415-763-5476
> AIM: jgrn307, MSN: jgrn307@hotmail.com, Gchat: jgrn307
>
>        [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>



-- 
In God we trust, all others bring data.

	[[alternative HTML version deleted]]

