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