Can I parameterize the state of a quantum system given reduced density matrices describing its subparts?

As the simplest example, consider a set of two qubits where the reduced density matrix of each qubit is known. If the two qubits are not entangled, the overall state would be given by the tensor product of the one qubit states. More generally, I could write a set of contraints on the elements of a two-qubit density matrix to guarantee the appropriate reduced description.

Is there is a way to do this more elegantly and systematically for arbitrary bi-partite quantum systems? I'm particularly interested in systems where one of the Hilbert spaces is infinite dimensional, such as a spin 1/2 particle in a harmonic oscillator.