library(rPAex)
#
# for image "tiff" format, see rPAex document
# r <- fileRaster(tiff)
#
crop <- system.file("external/corn.csv", package="rPAex")
s<-read.csv(crop,header = TRUE)
r <- raster::rasterFromXYZ(s)
raster::image(r,useRaster=FALSE,main="Corn crop\nArea outside the experimental plot")
#P<-locator(3)
P<-list(x=c(287705.9, 287711.8, 287712.8),y=c(8664203.15, 8664204.8, 8664193.67))
Q<-fourPoint(P)
polygon(Q,border = "blue")
text(Q)
models<-borderPoint(Q)
out<-externalPoint(models,r)
points(out[,1],out[,2],pch=".")
Run the code above in your browser using DataCamp Workspace