私が使用してラスタ画像のヒストグラムやグラフのように、値の分布図を作成したいR

Question

ラスタのヒストグラムを作成します。

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

plot()は私のデータを検証するために使用され、私にすべての緑を示し画像。おそらく3のバンド1 しかし、私は他のバンドを見ることができないのですか？明らかに、ヒストグラム内の分布は画像ファイル内の補間値を表さない。 ARCgisで作成されたヒストグラムはここにhistです。これは真の値を表していると信じています。

imageのような実際の値のヒストグラムを作成する方法に関する提案。

ベスト、 マティアス

Answer 1

あなたは試みることができる

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)) 

ggplot(melt(DEM), 
     aes(value, fill=as.factor(Var3))) + 
    geom_histogram(position="dodge") 

---

それとも、あなたの更新に関して

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)) 

Answer 2

がHIST（）関数を呼び出したりlukeAからggplotコマンドを使用し、その後、背景色にNAvalueを設定してみてくださいとして：

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

これは次のプロットになります：