在R中创建光栅直方图

我想使用以下方法创建光栅图像的值分布图，如直方图或图形：

library(raster)
library(sp)
library(rgdal)
DEM <- raster("NR.tif")
hist(DEM)
plot(DEM)

库（光栅）
图书馆（sp）
图书馆（rgdal）
你可以试试看

download.file("https://www.dropbox.com/s/t279m5ojners7fl/NR.tif?dl=1", 
              tf <- tempfile(fileext = ".tif"), mode="wb")
library(raster)
library(tiff)
library(ggplot2)
library(reshape2)
DEM <- readTIFF(tf)
plot(as.raster(DEM))



或者，关于你的更新

r <- as.raster(DEM)
tab <- as.data.frame(sort(table(r)))
ggplot(subset(tab, !r %in% c("#F0F0F0", "#000000")), 
              aes(x=r, y=Freq, fill=I(r))) + 
         geom_bar(stat="identity") + 
  theme(axis.text.x = element_text(angle=90))

r尝试将NavValue设置为背景色，然后调用hist（）函数或使用lukeA的ggplot命令：

library(raster)
ras <- stack("Downloads/NR.tif")
NAvalue(ras) <- 240
hist(ras)

库（光栅）
ras我不确定我是否理解了您想要的内容：levelplot{rasterVis}（）是否有帮助？不确定，如果我理解正确，但是：尝试
hist（DEM）
或
plot（density（DEM））
我希望编辑澄清了问题。谢谢，这确实是一个改进。但是直方图计算了灰色背景，因此不代表计算值。（这是垂直水柱中的硝酸盐梯度。）使用ARCmap创建的直方图具有明显不同的分布（）。是的，这是一个经典的RGB直方图，可以在每个图像处理工具中找到。Imho没有要求进行地理空间分析。但是你应该能够从上面得到你想要的地方（例如，你可以排除一个灰色的背景色。）谢谢。我想说得更具体些，我想我得更努力些。感谢您的精彩回答。另请参阅：
plotRGB（ras）
library(raster)
ras <- stack("Downloads/NR.tif")
NAvalue(ras) <- 240
hist(ras)