Geometric statistics with subspace structure preservation for SPD matrices