2016-04-08 5 views
2

スタック:ラスタを組み合わせることで、次のように私は、ファイルに読んだ私はCRU TSから3.23</p> <p>を温度データを含む複数のラスタのスタックを持って

nc = stack("cru_ts3.23.1951.1960.pre.dat.nc", varname = "pre") 
nc2 = stack("cru_ts3.23.1961.1970.pre.dat.nc", varname = "pre") 
... 
nc7 = stack("cru_ts3.23.2011.2014.pre.dat.nc", varname = "pre") 

は、私はこれらすべてのラスタを作成するために組み合わせることができます任意の方法はあります一つの大きなもの?私は降水量の傾向を見ています。彼らがすべて1つのオブジェクトに入っていると、ずっと簡単です!ありがとう

答えて

2

raster::stack(nc, nc2, nc3, nc4, nc5, nc6, nc7)トリックを行う必要があります。非常に類似した手法が、例えば、 here

はしかし、あなたが希望のオブジェクトは、同一の座標参照システム(CRS)、解像度(res)、およびextentを持っていることが要求されるstackすることに注意してください。そうしないと、stackが失敗する可能性が高くなります。上記の基準の1つ(または複数)が異なる場合は、raster::projectRasterraster::resampleを展開して、均一な空間プロパティを持つオブジェクトを作成できます。

+1

また、次のように1行でスタッキングを行うこともできます: 'stack(" cru_ts3.23.1951.1960.pre.dat.nc "、" cru_ts3.23.1951.1970.pre.dat.nc "、" cru_ts3.23.1951 .2014.pre.dat.nc "、varname =" pre ")' – maRtin

+1

True story @maRtin、私はちょうど上記の作成されたオブジェクトに固執したいと思います。 – fdetsch

関連する問題