[R-sig-Geo] (no subject)

[Tim Keitt]

We need R wrappers for:

http://www.vividsolutions.com/jts/jtshome.htm

or

http://geos.refractions.net/

My preference would be the latter.

THK

On Thu, 2006-04-06 at 16:38 -0400, Andy Bunn wrote:
> List,
> 
> I want to clip (subset) an sp object SpatialGridDataFrame using a polygon.
> For instance in this example I want to create a new object of class
> SpatialGridDataFrame that is clipped to the area inside the polygon on the
> map.
> 
>   data(meuse.grid)
>   coordinates(meuse.grid) <- ~x+y
>   gridded(meuse.grid) <- T
>   x = as(meuse.grid, "SpatialGridDataFrame")
>   clip <- cbind(x=c(180000,180000,180500,180500,180000),
>                 y=c(330500,331000,331000,330500,330500))
>   image(x, "dist")
>   polygon(clip)
>   # y = x[...clip?]
> 
> 
> How can I go about it?
> 
> Thanks in advance, Andy
> 
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
-- 
Timothy H. Keitt
Assistant Professor
http://www.keittlab.org/
http://www.utexas.edu/directory/index.php?q=Keitt