我有这个来自某人的mat文件,我想把它改成R中的栅格砖。
导入正常,我可以使用R.matlab
包
library(raster)
library(rgdal)`enter code here`
library(R.matlab)
SM2RF<-readMat("/..../Nile_DISTR.mat")
str(SM2RF)
List of 3
$ lonubn: num [1:231, 1] 34.6 34.6 34.6 34.6 34.6 ...
$ latubn: num [1:231, 1] 10.12 9.88 9.62 9.38 9.12 ...
$ Psim : num [1:8766, 1:231] NaN NaN NaN NaN NaN NaN NaN N NaN ...
- attr(*, "header")=List of 3
..$ description: chr "MATLAB 5.0 MAT-file, Platform: PCWIN64, Created
on: Fri Dec 19 10:10:41 2014 "
..$ version : chr "5"
..$ endian : chr "little"
列表中包括lat,long和Psim。
我曾经使用栅格砖功能转换为栅格砖。但是,我不知道matlab中栅格的结构,它应该是24年来每日降雨网格的时间序列。但我不知道为什么,我无法识别此mat文件中的图层数量。我尝试过像数组一样的砖,比如
r <-brick(SM2RF$Psim, xmn = min(SM2RF[[1]]),xmx=max(SM2RF[[1]]), ymn = min(SM2RF[[2]]), ymx=max(SM2RF[[2]]),crs=CRS("+init=epsg:20136"))
但是,我发现了这个错误。
Error in .local(x, ...) : cannot coerce a matrix to a RasterBrick
无论如何可以帮助我吗? 最好的问候
答案 0 :(得分:1)
您可以尝试以下方式:
x <- cbind(SM2RF$lonubn, SM2RF$latubn, t(SM2RF$Psim))
b <- rasterFromXYZ(x)