为时间序列中的每个值创建 netCDF 文件

问题描述 投票:0回答:0

我有投影坐标系中特定位置的流量时间序列数据 (692662.97,4156787.11)。我还有一张包含以下信息的栅格地图:

class      : RasterLayer 
dimensions : 1136, 624, 708864  (nrow, ncol, ncell)
resolution : 9.369195, 9.369195  (x, y)
extent     : 687763.9, 693610.3, 4151329, 4161972  (xmin, xmax, ymin, ymax)
crs        : NA 
source     : case01.dem 
names      : case01

时间序列是这样的:

        Time     value
1 0.00000000  226.5348
2 0.01666667  566.3369
3 0.03333333 3964.3586
4 0.05000000 5097.0325
5 0.06666667 6088.1221
6 0.08333333 7928.7172

我的数据是:

 df <- data.frame(Time = c(0, 0.016666667, 0.033333333, 0.05, 0.066666667, 0.083333333, 0.1, 0.116666667, 0.133333333, 0.15),
                     value = c(226.534776, 566.33694, 3964.35858, 5097.03246, 6088.122105, 7928.71716, 8183.568783, 8070.301395, 7928.71716, 7702.182384))

所以,作为输出,我想要 6 个不同的矩阵(或 nc 文件),尺寸为 1136x624,每个文件将显示栅格地图中特定单元格(692662.97,4156787.11)的这些值。

r raster netcat
© www.soinside.com 2019 - 2024. All rights reserved.