library(rPAex)
crop <- system.file("external/cassava.csv", package="rPAex")
s<-read.csv(crop,header = TRUE)
r <- raster::rasterFromXYZ(s)
#-----
# P<-locator(3)
P<-list(x=c(287698.21, 287700.99, 287702.39), y=c(8664200.68, 8664201.57,8664190.63))
Q<-fourPoint(P)
#-----
nPl<-11; nPw<-3; long=1; wide=0.9
raster::image(r,useRaster=FALSE,main="Cassava crop\nnear infrared image")
img<-imageField(r, P=Q, nPl, nPw, long, wide, ploting = TRUE, border="blue",lwd=1)
R<-img$Qbase
head(R)
q<-agricolae::tapply.stat(R[,2:3],R[,1],mean)
text(q[,2],q[,3],q[,1],cex=1)
Run the code above in your browser using DataLab