Question - DWT 3D Wavelet Coefficients Organization

[ktaswell]

I'm trying to implement a soft thresholding method for 3D data which preserves scaling (or approximation) coefficients at each level, then soft thresholds all others. Unfortunately, I have been having trouble identifying indices which correspond to the scaling coefficients. I was trying to make sense of the detailn() and detailindex() functions and it makes sense to me for 1D condition, but I don't understand how to use these to create an index list for higher dimensions, like the 3D condition that I need.

Also to complicate further, I need to handle both dyadic cube and dyadic non-cube dimensions.

Could someone help me with this?

Thanks