At least two chemical species are imaged using magnetic resonance imaging with signal separation for the two chemical species. The method includes acquiring first and second echo data at different echo times resulting in a first and second acquired complex dataset, —modelling the first and second acquired dataset, said modelling comprising a spectral signal model of at least one of the chemical species, —identifying in the first and second acquired dataset the voxels for which the modelling yields a single, unambiguous mathematical solution for the signal separation, and —resolving the ambiguity for the voxels for which the modelling yields more than one mathematical solution, if any such voxels remain.