As a new method for measuring the spatial distribution of Bose-Einstein condensates, the magnetic resonance imaging (MRI) method is proposed and studied in detail. The basic concepts, the resolution limit and the formalism of the MRI method are presented. It is expected that a resolution higher than that in optical imaging methods can be obtained by using the MRI method. Results of simulation of expected MRI signals for Bose-Einstein condensates containing dark solitons are also presented.