Dear all,



I'm running on R 2.7.1 (win), and I'm trying to "clip" a psp and an owin
object, without success. Here is the code that I'm using:



############################################

library(maptools)

library(spatstat)



a<-readShapeSpatial("area.shp")

b <- as(as(a, "SpatialPolygons"), "owin")

plot(b)



g<-0:360



            pt<-runifpoint(1, b)

            pt_v<-as.data.frame(pt)



            x<-pt_v[,1]

            y<-pt_v[,2]



            dest<-gcDestination(x, y, bearing=sample(g, 1), dist=4,
dist.units = "km")

            trans <- psp(x, y, dest[,1],dest[,2],  window=owin())



            plot(trans, lwd=2, add=TRUE)



###############################################



My question is what should I do to get the information of what is inside the
clipping object (trans), having as an input object an owin (b).



Sorry if the code looks silly (and/or the question!), but I'm a nearby.



Thanks in advance for any help,



Cheers,



Alexandre Camargo Martensen

**********************************

Landscape Ecology and Conservation Lab

University of Sao Paulo

Sao Paulo – Brazil

**********************************

	[[alternative HTML version deleted]]


