27

次の最小限の作業例のように、デカルト座標 (x、y) に対して値 (l) が指定されているデータ フレームがあります。

set.seed(2013)
df <- data.frame( x = rep( 0:1, each=2 ),
                  y = rep( 0:1,  2),
                  l = rnorm( 4 ))

df
#   x y           l
# 1 0 0 -0.09202453
# 2 0 1  0.78901912
# 3 1 0 -0.66744232
# 4 1 1  1.36061149

ラスター パッケージを使用してラスターを作成したいのですが、ドキュメンテーションを読んでも、ラスター セルにある形式でデータをロードする簡単な方法がわかりませんでした。for ループを使用してそれを行ういくつかの方法を思いつきましたが、私が見逃しているもっと直接的なアプローチがあるのではないかと思います。

4

3 に答える 3

54

より簡単な解決策は次のとおりです。

 library(raster)
 dfr <- rasterFromXYZ(df)  #Convert first two columns as lon-lat and third as value                
 plot(dfr)
 dfr                  
 class       : RasterLayer 
 dimensions  : 2, 2, 4  (nrow, ncol, ncell)
 resolution  : 1, 1  (x, y)
 extent      : -0.5, 1.5, -0.5, 1.5  (xmin, xmax, ymin, ymax)
 coord. ref. : NA 
 data source : in memory
 names       : l 
 values      : -2.311813, 0.921186  (min, max)

プロット

さらに、CRS 文字列を指定することもできます。詳細な議論はここで入手できます。

于 2016-06-09T07:11:47.490 に答える
33

これが1つのアプローチです。SpatialPixelsDataFrame

library(raster)
# create spatial points data frame
spg <- df
coordinates(spg) <- ~ x + y
# coerce to SpatialPixelsDataFrame
gridded(spg) <- TRUE
# coerce to raster
rasterDF <- raster(spg)
rasterDF
# class       : RasterLayer 
# dimensions  : 2, 2, 4  (nrow, ncol, ncell)
# resolution  : 1, 1  (x, y)
# extent      : -0.5, 1.5, -0.5, 1.5  (xmin, xmax, ymin, ymax)
# coord. ref. : NA 
# data source : in memory
# names       : l 
# values      : -0.6674423, 1.360611  (min, max)

help('raster')では、さまざまなクラスのオブジェクトからラスターを作成するいくつかの方法について説明しています。

于 2013-10-28T04:56:03.417 に答える