在R中添加大型柵格時，@ data @值會發生什麼？

Question

簡單的算術運算工作與小柵格預期，例如：

> r1 <- raster(ncol=5, nrow=5) 
> r1[] <- rep(1, 5 * 5) 
> class(r1@data@values) 
[1] "numeric" 
> class((r1 + r1)@data@values) 
[1] "numeric" 

但是，當我用更大的柵格，奇怪的事情發生了：

> r1 <- raster(ncol=5000, nrow=5000) 
> r1[] <- rep(1, 5000 * 5000) 
> r2 <- r1 + r1 
> class(r2@data@values) 
[1] "logical" 

爲什麼r2@data@values成爲logical？ 該所得目的是在隨後的計算中，不可用的 和一般它的行爲出現異常，例如：

> summary(r1@data@values) 
    Min. 1st Qu. Median Mean 3rd Qu. Max. 
     1  1  1  1  1  1 
> summary(r2@data@values) 
    Mode NA's 
logical  0 
> object.size(r1) 
200011576 bytes 
> object.size(r2) 
12112 bytes 

基本上r2看起來破碎。

作爲一種變通方法，加入@data@values似乎很好地工作：

> r3 <- r1 
> r3@data@values <- r1@data@values * 2 
> class(r3@data@values) 
[1] "numeric" 

但你怎麼能解釋這種行爲，只有擁有較大的柵格發生什麼呢？

Answer 1

logical是因爲raster對象指向存儲實際值的磁盤上的臨時文件，而不是將25,000,000個數值存儲在內存中。這些值完全可用於後續的算術方法。使用@data@values插槽來設置值需要麻煩。使用提供的算術方法或使用setValues()函數。

這在?Arith-methods幫助頁面中詳細說明：

如果輸出柵格*的值不能在內存中舉行，他們將被保存到一個臨時文件。您可以使用選項來設置默認文件格式，數據類型和進度欄。

鑑於上面的對象r2，您可以像存儲在內存中的對象那樣使用它。這是raster所做的很多事情之一 - 它可以與存儲在磁盤上的文件無縫協作，就像它們是常規的內存中對象一樣。

r3 <- r2^2 
r3 
#class  : RasterLayer 
#dimensions : 5000, 5000, 2.5e+07 (nrow, ncol, ncell) 
#resolution : 0.072, 0.036 (x, y) 
#extent  : -180, 180, -90, 90 (xmin, xmax, ymin, ymax) 
#coord. ref. : +proj=longlat +datum=WGS84 
#data source : /private/var/folders/ck/1x5j5jvx5mq17clf5r32gg540000gn/T/R_raster_tmp/Simon/raster_tmp_2014-01-13_214043_710579.grd 
#names  : layer 
#values  : 4, 4 (min, max) 

注層的值被正確地返回爲4，即2^2。