A theoretical framework for the generation of natural orbitals, natural spin orbitals, as well as orbital-and spin-magnetizations from multi-configurational ab-initio wavefunction calculations including spin-orbit coupling is presented. It is shown how these computational orbital and magnetization tools can be used to interpret and rationalize the magnetic properties of selected complexes containing transition metals, lanthanides, and actinides.