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 > [email protected] > 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 _______________________________________________ R-sig-Geo mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-sig-geo
