Hello,
I am trying to interpolate a variable within my dataset over a couple regions in Brazil. However, I am having a difficult time trying to get the box around my interpolated area to expand so that I can see the entire region. I am not sure if any one has had this issue before, but I would greatly appreciate some advice. I have attached my code and a screenshot of the plot I have been generating.
library(dismo)
BM <- CRS("+proj=aea +lat_1=-5 +lat_2=-42 +lat_0=-32 +lon_0=-60 +x_0=0 +y_0=0 +ellps=aust_SA +units=m +no_defs")
topsoil_BM <- spTransform(topsoil_sp, BM)
cerrado_BM <- spTransform(cerrado, BM)
v <- voronoi(topsoil_BM)
plot(v)
brazil <- aggregate(cerrado_BM)
vbrazil <- intersect(v, brazil)
spplot(vbrazil, 'd13C', col.regions=rev(get_col_regions()))