My solution:

Calculating and applying the VDMs via script externally and using the resulting unwarped time series as the base functional dataset in CONN. I changed the "indirect" pipeline by removing VDM calculation und Realign&Unwarp and inserted a simple Realign step. Not optimal, because of the additional interpolation, but perhaps the only reasonable way.

Best,
Martin