Hello all. I have a bunch of 3D scans of the same subject over time joined in a 4D scan using (by `JoinSeries`

) and I have a vector with the phase weights. To create a 3D approximation of this 4D object you might want to multiply with the (normalized) phase weights and then take the sum over that axis.

Is there a one-liner for this in (Simple)ITK or should I resort to turning the vector into an image with the same dimensions, multiply and sum over the axis? The ProjectionImageFilter’s do not seem to support weights.

As an followup question, I’d like to know if there is somehow a way to assign these non-uniform “temporal spacings” in ITK? Currently I use SetMetaData where I set the phase as an additional key, but having non-uniform spacings would be convenient. I see that basically the xarray implementation between ITK <-> xarray basically sets a position *per slice*. Is there such ITK support?