Hello,

I am attempting to extract values from each layer of a 276 layer raster
brick which should produce a one column, 276 row data frame containing EVI
values for each layer.

I am creating a SpatialPoints object from gCentroid using a shapefile with
one polygon hence the output has one point.  I immediately feed this into
extract but receive an error which I don't understand.  See below.  Any
help is MUCH appreciated!!

smPoint <- gCentroid(sm_reproj, byid=T)
smModis <- extract(stack,smPoint,df=T)

Error in `colnames<-`(`*tmp*`, value = names(x)[layers]) :
  length of 'dimnames' [2] not equal to array extent

	[[alternative HTML version deleted]]

