この問題の記述方法は不明です。私はそれが自明だと感じているが、私はそれを握ることができない。rasterToPoints抽出のx、y値を論理的に受け入れることができません
私はラスターオブジェクトのスタック(オブジェクトNDVI
)を持っています。これらのことから、私は、私はいくつか使用できなくピクセル追い出さ前処理に伴いrasterToPoints
xycoord1 <- rasterToPoints(NDVI)
xycoord <- xycoord1[,c(1:2)]
を使用してx
とy
座標を抽出してしまった:
> str(xycoord.short)
num [1:20054, 1:2] 3802292 3802523 3802755 3802987 3803218 ...
- attr(*, "dimnames")=List of 2
..$ : NULL
..$ : chr [1:2] "x" "y"
私は単に特定のx
を見つけたいんやありませんy
座標。例:
> which(xycoord.short[,1]==3802292)
integer(0)
しかし、たとえば1つの列内の値を「取得」できないようです。
> xycoord.short[,1][1]
[1] 3802292
> xycoord.short[,1][1]==xycoord.short[,1][1]
[1] TRUE
> xycoord.short[,1][1]==3802292
[1] FALSE
誰でもこの問題にお手伝いできますか?私は問題を見つけることができません。それはrasterToPoints
を通して最初の抽出と関係がありますか?ありがとう!
EDIT:私のxy座標の最初の10行
xy <- structure(c(3802291.63636448, 3802523.29272274, 3802754.94908101,
3802986.60543927, 3803218.26179754, 3803449.9181558, 3803681.57451406,
3803913.23087233, 3804144.88723059, 3804376.54358886, -49690.2888476191,
-49690.2888476191, -49690.2888476191, -49690.2888476191, -49690.2888476191,
-49690.2888476191, -49690.2888476191, -49690.2888476191, -49690.2888476191,
-49690.2888476191), .Dim = c(10L, 2L), .Dimnames = list(NULL,
c("x", "y")))
EDIT2ため dput出力: 値が明らかに丸みを帯びているように、それは理にかなってdput出力を掲示した後。 正確な数字を使用すると...
> any(xycoord.short[,1]==3802291.63636448)
[1] TRUE
Ah!問題があります。私の答えを編集させてください。 –