Hi Rainer,

Is the Earth Movers Distance used to compare two distributions (e.g. the
result of hist function) or two exactly raster grids? actually when I
compared two raster grids(I pasted my script here):

b <-raster(paste("tavg.3MIN.1981-2010_yearmean_NE30mask.nc
",sep=""),varname='tavg')
a <-raster(paste("prcp.3MIN.1981-2010_yearmean_NE.nc
",set=""),varname='prcp')
a <- as.matrix(a)
b <- as.matrix(b)
emd2d(a, b)

This function runs forever to finish (and it used up the memory).

Thanks,

Ping

-- 
Ping Yang, Ph.D.
CUNY Environmental Crossroads Initiative
The City College of New York - CUNY

	[[alternative HTML version deleted]]

