空间数据和内存

我想加起来geotiffs，但遇到内存问题。 R为使用所有32GB的下列R-错误...空间数据和内存

In writeValues(y, x, start = 1) : 
    Reached total allocation of 32710Mb: see help(memory.size)

我还检查R的性能，这是64位和目标是...

“C：\ Program Files文件\ r \ R-3.3.0 \ BIN \ 64 \ Rgui.exe”

的版本是

R.Version() 
$platform 
[1] "x86_64-w64-mingw32" 

$arch 
[1] "x86_64" 

$os 
[1] "mingw32" 

$system 
[1] "x86_64, mingw32" 

$status 
[1] "" 

$major 
[1] "3" 

$minor 
[1] "3.0" 

$year 
[1] "2016" 

$month 
[1] "05" 

$day 
[1] "03" 

$`svn rev` 
[1] "70573" 

$language 
[1] "R" 

$version.string 
[1] "R version 3.3.0 (2016-05-03)" 

$nickname 
[1] "Supposedly Educational"

所以看起来正在使用R.我最大内存我试过用bigmemory打包在R中。所以在下面的代码中，我尝试了改变th e矩阵转换为big.matrix，但失败并且在尝试写入输出文件时发生错误。任何关于尝试修改代码的建议，以便使用更少的内存或尝试在ff或bigmemory软件包中工作？

############ LOOP THROUGH AGE MAPS TO COMPILE THE NUMBER OF TIMES A CELL BURNS DURING A GIVEN SPAN OF TIME #################### 
## Empirical Fires 
print("1 of 3: 2010-2015") 

burn.mat<- matrix(0,nrow,ncol) #create matrix of all zero's, the dimension of your landscape (row, col) 

# Read in Historical Fire maps 
for (j in 2010:2015){ #Year Loop 

    age.tmp<- as.matrix(raster(paste('fr',j,'.tif',sep=''))) #read in Age Map 

    burn.mat<- burn.mat+(age.tmp==1) #when something has burned in ALFRESCO empirical fire history files, AGE=1. (age.tmp==0) is a 'logic' cmd, returning a 0,1 map for True/False 

    #Write the data to a geotiff 
    out <- raster(burn.mat,xmn=-1692148,xmx= 1321752, ymn = 490809.9, ymx = 2245610, crs = '+proj=aea +lat_1=55 +lat_2=65 +lat_0=50 +lon_0=-154 +x_0=0 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=m +no_defs') 
    writeRaster(out,filename=paste(outdir,'/burn.mat.hist.1950-2007.tif',sep=''),format = 'GTiff',options='COMPRESS=LZW',datatype='FLT4S',overwrite=T) 
}

来源

2016-08-23 Jen

为什么代码每次通过循环写入相同的文件名？ – mdsumner

的问题可能会自动消失，如果你使用的光栅*的对象，而不是矩阵。喜欢的东西

library(raster) 
r <- raster('fr2010.tif') 

burn.mat <- setValues(r, 0) 

for (j in 2010:2015) { 
    age.tmp <- raster(paste0('fr', j, '.tif')) 
    burn.mat <- burn.mat + (age.tmp==1) 
# if age.tmp only has values of 0 and 1 use this instead: 
# burn.mat <- burn.mat + age.tmp 
} 

# write the results outside of the loop 
writeRaster(burn.mat, filename=file.path(outdir, 'burn.mat.hist.1950-2007.tif'), options='COMPRESS=LZW',datatype='FLT4S',overwrite=TRUE)

更直接的方式没有环

files <- paste0('fr', 2010:2015, '.tif')) 
s <- stack(files) 
burn <- sum(s)

或者

burn <- sum(s == 1)

或者写入到一个文件中的一个步骤

b <- calc(s, sum, filename=file.path(outdir, 'burn.mat.hist.1950-2007.tif'), options='COMPRESS=LZW', datatype='FLT4S', overwrite=TRUE)

来源

2016-08-24 06:05:23 RobertH

空间数据和内存

回答

相关问题