Hello everyone,

I'm new to this list and to using R for geographical data analysis. I am
playing around with the gdistance package but have a small problem:

I am trying to use a RasterLayer to define origin and goal points for the
passage function, which I believe is possible based on the documentation:

origin
SpatialPoints, matrix or numeric object with coordinates or RasterLayer
object with origin cells set to TRUE

goal
SpatialPoints, matrix or numeric object with coordinates or RasterLayer
object with origin cells set to TRUE

Here is a basic idea of what I'm doing:

####################################
map.test <- raster(ncol=10,nrow=10)
map.test[] <- 100

origin <- raster(ncol=ncol(map.test),nrow=nrow(map.test))
origin[] <- FALSE
origin[93] <- TRUE
goal <- raster(ncol=ncol(map.test),nrow=nrow(map.test))
goal[] <- FALSE
goal[9] <- TRUE

tran.test <- transition(map.test,transitionFunction=mean,directions=8)

passage(tran.test, origin, goal,theta=0,totalNet=total)

and this gives the error message:
Error in which(getValues(origin)) : argument to 'which' is not logical
##################################

Any ideas what I'm doing wrong?

Thanks very much
Josephine

	[[alternative HTML version deleted]]

