Magnetorelaxometry (MRX) is a well-known measurement technique which allows the retrieval of magnetic nanoparticle (MNP) characteristics such as size distribution and clustering behavior. This technique also enables the non-invasive reconstruction of the spatial MNP distribution by solving an inverse problem, referred to as MRX imaging. Although MRX allows the imaging of a broad range of MNP types, little research has been done on imaging different MNP types simultaneously. Biomedical applications can benefit significantly from a measurement technique that allows the separation of the resulting measurement signal into its components originating from different MNP types. In this paper, we present a theoretical procedure and experimental validation to show the feasibility of MRX imaging in reconstructing multiple MNP types simultaneously. Because each particle type has its own characteristic MRX signal, it is possible to take this a priori information into account while solving the inverse problem. This way each particle type's signal can be separated and its spatial distribution reconstructed. By assigning a unique color code and intensity to each particle type's signal, an image can be obtained in which each spatial distribution is depicted in the resulting color and with the intensity measuring the amount of particles of that type, hence the name multi-color MNP imaging. The theoretical procedure is validated by reconstructing six phantoms, with different spatial arrangements of multiple MNP types, using MRX imaging. It is observed that MRX imaging easily allows up to four particle types to be separated simultaneously, meaning their quantitative spatial distributions can be obtained.