# NOT RUN {
## The function is currently defined as
function (a = NIR, b = Red, c = Blue,Pixel.Depth)
{
name=EVI(a,b,c,Pixel.Depth=null)
areaXY <-c(xmin, xmax, ymin, ymax)
cropXY <- crop(name, areaXY)
plot(cropXY,lwd=4, main="EVI",xlab="easting", ylab="northing")
hist(cropXY, main="EVI",
xlab="EVI",col="red", ylab="Frequency of Pixels")
}
# }
Run the code above in your browser using DataLab