3

ラスター オブジェクト (tempMap) に変換した後、マトリックス (temp_matrix) をプロットしたいと考えています。さらに、緯度と経度の位置が同じプロット ウィンドウで利用できるポイントをいくつか追加したいと考えています。私はいくつかのアプローチを試しましたが、利用可能なポイントは緯度/経度の特定の場所であり、取得しているラスター オブジェクトの範囲が異なるため、どれも機能しないようです。この問題で私を助けてください。以下に、問題のサンプル データを示します。

library(raster)
temp_matrix<-array(NA,c(11,11))
temp_matrix[1,]<-c(NA,NA,NA,NA,NA,NA,NA,0,0,-6,-6)
temp_matrix[2,]<-c(0,0,0,0,NA,NA,1,0,0,0,0)
temp_matrix[3,]<-c(1,0,0,-1,-1,0,0,0,1,0,0)
temp_matrix[4,]<-c(1,1,0,0,0,0,-1,-1,0,0,0)
temp_matrix[5,]<-c(1,NA,NA,NA,NA,-1,-1,-1,0,-1,-1)
temp_matrix[6,]<-c(NA,NA,NA,NA,NA,NA,-1,-1,-1,0,0)
temp_matrix[7,]<-c(NA,NA,NA,NA,NA,NA,NA,0,-1,0,0)
temp_matrix[8,]<-c( NA,NA,NA,NA,NA,NA,NA,0,0,0,-1)
temp_matrix[9,]<-c(NA,NA,NA,NA,NA,NA,NA,-1,0,0,0)
temp_matrix[10,]<-c(NA,NA,NA,NA,NA,NA,NA,NA,-1,-1,-2)
temp_matrix[11,]<-c(NA,NA,NA,NA,NA,NA,NA,NA,-2,-3,-2)
plot(raster(temp_matrix),axes = FALSE,legend=FALSE)
tempMap <- raster(temp_matrix)

# plot the points over this raster layer
point_1<-c(10,10) # should appear on 2nd row from top i.e. over temp_matrix[2,10]
point_2<-c(9,10)  # should appear on 3rd row from top i.e. over temp_matrix[3,10]
point_3<-c(1,10)  # should appear on lowermost row i.e.over temp_matrix[11,10]
4

1 に答える 1

6

デフォルト値に依存するのではなく、目的の範囲を設定する必要があります。

raster国家の助けとして

## S4 method for signature 'matrix'
raster(x, xmn=0, xmx=1, ymn=0, ymx=1, crs=NA, template=NULL)

と

  • xmn:: 最小 x 座標 (左ボーダー)

  • xmx:: 最大 x 座標 (右ボーダー)

  • ymn:: 最小 y 座標 (下の境界線)

  • ymx:: 最大 y 座標 (上の境界線)

xmn、xmx、ymnおよびymxを希望する値 (この場合は 1、11、1、11) に設定できます)

tempMap <- raster(temp_matrix, xmn = 1, xmx = 11, ymn = 1, ymx=11)
plot(tempMap,axes = FALSE,legend=FALSE)
points(c(10,9,1), c(10,10,10))

于 2013-02-05T23:38:17.087 に答える