Method of magnetic resonance imaging comprising the steps of: exciting magnetic resonant (MR) active nuclei in a region of interest; creating magnetic field gradients in a phase-encode direction for spatially encoding the excited MR active nuclei; the number of phase-encode gradients producing a field of view corresponding to the region of interest; receiving radio frequency data from the region of interest using an array of at least two radio frequency receive coils, and producing from the radio frequency data an image by combining the images S.sub.i obtained from the radio frequency signals of the single coils C.sub.i of the array by weighting the said images S.sub.i of each coil C.sub.i with the sensitivity B.sub.i of the corresponding coil C.sub.i; an estimated sensitivity B.sub.i (Fourier domain) of each coil C.sub.i of the array being used, which is determined by applying a smoothing filter f.sub.i to the k-space data s.sub.i received from each coil C.sub.i. According to the invention the said f.sub.i smoothing filter being estimated from the spatial frequencies content of the k-space sensitivity map of each coil, and the filter f.sub.i being a k-space filter with a triangular or pyramidal shape, depending on the fact that a single dimensional profile or a two dimensional image is reconstructed, the triangular filter having a width and the pyramidal filter having a base which is double than the maximum spatial frequency of the relative coil field.