Hello everybody:

 

I´ve got rasters, imported using the function raster() from the package
raster grouped in stacks.

 

I also have SpatialPolygonsDataFrame, which are shapefiles imported using
the function readShapePoly() of maptools.

 

I need to obtain the frequency distribution for the values of each raster
within the stack, but only for the area falling within the polygon.

 

How can I select the area of the raster within the polygon? Or can I clip
this area (how?) and then use freq(), of the raster package? 

 

I would appreciate any help.

 

Thank you,

 

Frederico Mestre

 


	[[alternative HTML version deleted]]

