我一直在使用 R 中的数据框,其中包含坐标列(其值几乎是不规则间隔的,因为我不得不将它们从http://www.arcgis.com/home/item.html?id=5771199a57cc4c29ad9791022acd7f74中的分钟-秒格式转换为十进制)和与每个坐标相关联的值列,如下所示:
Latitude Longitude Elevation MAT MWMT MCMT TD MAP MSP AHM SHM DD_0 DD5 DD_18 DD18 NFFD
54.99285 -129.9792 -2 6.8 14.4 -1.7 16.1 2473 696 6.8 20.6 283 1384 4100 30 246
54.99285 -129.9708 10 6.7 14.4 -1.8 16.1 2456 691 6.8 20.8 287 1383 4109 30 245
54.99285 -129.9625 5 6.8 14.4 -1.8 16.2 2431 686 6.9 21.0 286 1392 4097 31 246
# ... and so on.
我想构建一个与每列对应的层的 RasterBrick。最初,我认为,由于我的数据框命名为'clim_df':
clim_brick <- rasterFromXYZ(clim_df, crs = NA)
我希望您能够将每个列转换为砖块内的栅格层,但在此之后出现了以下错误(在将纬度和经度列名称更改为“y”和“x”之后):
Error in rasterFromXYZ(clim_df) : x cell sizes are not regular
从https://www.rdocumentation.org/packages/raster/versions/2.6-7/topics/rasterFromXYZ看来,这个函数需要坐标间距相等才能正常工作。
有没有什么方法可以解决这个问题?任何帮助将不胜感激。
P.S. 我以前用过Stack Overflow来解决问题,但从未直接提问。如果我格式不正确或者没有提供足够的信息,请让我知道。谢谢!