2017-02-20 10 views
0

ラスタデータからマップのプロットを作成しようとしています。私はこのコードを使用しています:RColorBrewerをプロットするときのエラー

library(raster) 
library(rgdal) 
library(classInt) 
library(RColorBrewer) 

NDII = raster("G:\\Sheyenne\\image_differencing\\NDII\\differenced.tif") 
value.vector = values(NDII) 
breaks.qt = classIntervals(value.vector, n = 6, style = "jenks", intervalClosure = "right") 
print (plot(NDII, breaks=breaks.qt$brks, col = brewer.pal(6, "Set1"))) 

をこれが返されます。

Error in print(plot(NDII, breaks = breaks.qt$brks, col = brewer.pal(6, : error in evaluating the argument 'x' in selecting a method for function 'print': Error in .asRaster(x, col, breaks, zrange, colNA, alpha = alpha) : could not find function "brewer.pal"

答えて

2

をあなたは何の再現性の例を提供していない、と私はあなたのエラーを再現することはできません。 あなたがラスタとしてRロゴの再現性の例を使用するのと同じ(とNDII[]を使用しての代わりに、変数にvalues(NDII)を格納する近道で)で次のコードを、私のためだけで正常に動作...

library(raster) 
library(rgdal) 
library(classInt) 
library(RColorBrewer) 

NDII = raster(system.file("external/rlogo.grd", package="raster")) 
# next line is really slow, I'd advise to run crop(NDII, extent(0,20,0,20)) 
# before to make quick tests 
breaks.qt = classIntervals(NDII[], n = 6, style = "jenks", 
          intervalClosure = "right") 
plot(NDII, breaks=breaks.qt$brks, col = brewer.pal(6, "Set1")) 

このコードでエラーを再現していますか?多分新しい、新鮮なセッションから始めることができますか?私は今、私の答えで述べた

あなたはそれより完全な答えにするために、あなたのコードとOPの間の差異を明確にすることをお勧めします

enter image description here

+0

... – Benjamin

+0

まあ、何の違いが実際にありません... 。 – ztl

+0

あなたの例は私のために再現可能です。なぜ私のコードがうまくいかないのかまだ分かりませんが、新しいセッションではそれを修正することはできませんが、それは私の最終的な目的です。 –

関連する問題