如何在SciPy上使用interpn上采样上采样4d fMRI .mat文件？

Question

我有一个4d .mat文件，它为fMRI文件提供X，Y，Z和T​​，该文件以1.5,2和3hz扫描。我想插入较低的采样率图像，因此它们都是3hz。我曾尝试在scipy中使用interpn，但我不明白它的论点。其尺寸为80x64x33xn，具体取决于采样率。

Answer 1

由于您只是插入时间轴，因此最好使用interp1d。假设您的4d数据位于D且您的各个时间网格为T15，T2和T3，那么它应该像

一样简单

from scipy.interpolate import interp1d
# next line assumes interpolation along the last dimension of D
# if your data are arranged differently, use axis kwarg
f = interp1d(T15, D)
Dprime = f(T3)

其中Dprime是上采样数据。