我有一个看起来像这样的data.frame：
> head(dat)
YEAR DAY LAT LON BIOMASS
2000 8 39.78333 -73.33333 32.7
2000 8 39.86667 -73.46667 6.8
2001 7 40.10000 -73.15000 9.5
2001 7 40.23333 -72.73333 11.3
2001 8 40.25000 -73.26667 3.2
2002 7 40.41667 -72.73333 159.7
我希望将这些数据转换为BIOMASS的一系列2D矩阵（LAT x LON），用于每天的组合。我希望网格分辨率为0.10度。我知道在某些情况下，网格中会有不止一个观察结果。在这些情况下，我会采用重叠值的均值。此外，会有许多没有关联值的网格点，因为我想插入NA。
我的实际date.frame相当大（50年，每年约25天）所以我需要一些有效的东西。我试图在R中找到一个可以做到这一点的包，但没有运气。有没有人有这方面的经验？我可以使用Matlab，虽然我在R中更舒服。