Wie verwende ich interpn auf SciPy, um upsample upsample eine 4D fMRI .mat Datei?

Question

Ich habe eine 4D .mat-Datei, die X, Y, Z und T für eine fMRI-Datei gibt, die Scans bei 1,5, 2 und 3hz dauert. Ich möchte die Bilder mit der niedrigeren Abtastrate interpolieren, so dass sie alle bei 3 Hz sind. Ich habe versucht, interpn in scipy zu verwenden, aber ich verstehe die Argumente dafür nicht. Seine Abmessungen betragen 80x64x33xn, abhängig von der Abtastrate.

Answer 1

Da Sie nur die Zeitachse interpolieren, sind Sie wahrscheinlich besser mit interp1d bedient. Angenommen, Ihre 4d Daten sind in D und Ihre verschiedenen Zeitraster sind T15, T2 und T3 es sollte so einfach sein wie

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) 

wo Dprime die aufwärtsgesampelten Daten sind.