Method of MR imaging of an object positioned in an examination volume of a MR device, the method comprising the steps of: - subjecting the object to an imaging sequence of RF pulses and switched magnetic field gradients, which imaging sequence is a zero echo time sequence comprising: i) setting a readout magnetic field gradient having a readout direction and a readout strength, both staying substantially constant during a readout cycle; ii) radiating a RF pulse in the presence of the readout magnetic field gradient; iii) acquiring a FID signal in the presence of the readout magnetic field gradient, wherein the FID signal represents a radial k-space sample; iv) gradually varying the readout direction; v) sampling a spherical volume in k-space by repeating readout cycles formed by steps i) through iv) a number of times, with the readout strength being varied between repetitions, and the number of repetitions being such that the readout direction is varied until the spherical volume is covered with sufficient density; wherein the readout strength is varied such that individual positions in k-space are sampled at least two times, each time with a different value of the readout strength, such that said k-space position is sampled at two or more different sampling times and - reconstructing a MR image from the acquired FID signals, wherein signal contributions of two or more chemical species to the acquired FID signals are separated,wherein - the signal contributions of the two or more chemical species to the FID signals are derived from phase differences of the acquired FID signals induced by the variation of the readout strength and - the separation of the signal contributions is performed on the basis of a signal model including at least the MR spectrum of each of the chemical species. - subjecting the object to an imaging sequence of RF pulses and switched magnetic field gradients, which imaging sequence is a zero echo time sequence comprising: i) setting a readout magnetic field gradient having a readout direction and a readout strength, both staying substantially constant during a readout cycle; ii) radiating a RF pulse in the presence of the readout magnetic field gradient; iii) acquiring a FID signal in the presence of the readout magnetic field gradient, wherein the FID signal represents a radial k-space sample; iv) gradually varying the readout direction; v) sampling a spherical volume in k-space by repeating readout cycles formed by steps i) through iv) a number of times, with the readout strength being varied between repetitions, and the number of repetitions being such that the readout direction is varied until the spherical volume is covered with sufficient density; wherein the readout strength is varied such that individual positions in k-space are sampled at least two times, each time with a different value of the readout strength, such that said k-space position is sampled at two or more different sampling times and i) setting a readout magnetic field gradient having a readout direction and a readout strength, both staying substantially constant during a readout cycle; ii) radiating a RF pulse in the presence of the readout magnetic field gradient; iii) acquiring a FID signal in the presence of the readout magnetic field gradient, wherein the FID signal represents a radial k-space sample; iv) gradually varying the readout direction; v) sampling a spherical volume in k-space by repeating readout cycles formed by steps i) through iv) a number of times, with the readout strength being varied between repetitions, and the number of repetitions being such that the readout direction is varied until the spherical volume is covered with sufficient density; wherein - reconstructing a MR image from the acquired FID signals, wherein signal contributions of two or more chemical species to the acquired FID signals are separated,wherein - the signal contributions of the two or more chemical species to the FID signals are derived from phase differences of the acquired FID signals induced by the variation of the readout strength and - the separation of the signal contributions is performed on the basis of a signal model including at least the MR spectrum of each of the chemical species.