To produce image information from an object it is subjected to a continuous static magnetic field along a Z axis and to sets of sequences of orthogonal gradients G.sub.x, G.sub.y and G.sub.z to the magnetic field. Spins in a selected plane (the X-Z plane) are excited by selective rf pulses and an associated G.sub.y gradient and the selected spins are subjected to all three gradients of which the G.sub.z gradient provides twist or warp to each column of spins extending along the Z axis to phase-encode the columns. The spin-echo signals are read out in the presence of a G.sub.x gradient. In each set of sequences a different value of Z gradient is employed. The Fourier transformed spin-echo signals obtained from each sequence, when arranged in order of increasing G.sub.z gradient and subjected to a second Fourier transform represent the distribution of spin density in the Z direction, thus giving a two-dimensional image of the selected X-Z plane.