# NOT RUN {
{
require(raster)
# read raster data
r <- brick(system.file("extdata", "ndvi.tif", package="fieldRS"))
# read field data
data(fieldData)
# derive centroids
c <- spCentroid(fieldData)
# plot polygons and compare with centroids
plot(fieldData)
points(c, col="red")
}
# }
Run the code above in your browser using DataLab