# calculate polygon coordinates for each plot centroid. proj4string=CRS(as.character("+proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.9999079 +x_0 ...