use otbGeodesicMorphologyIterativeDecompositionFilter and create Vector image of 3*N bands for each level (leveling, convex, concave) * N. For mono scale decomposition set the number of iteration to 1.