Two Dimensional Spectral Factorization

Hello,

Does anyone know a function or a code to perform a 2D Spectral Factorization (of a complex signal) if it is possible.

It is quite easy when it is one-dimensional but it is not as easy when it is multidimensional.

Thanks in advance, Ignacio