A fully automatic, multiscale and multiblock fuzzy C-means (MsbFCM) classification method with intensity correction for MR images is presented in this paper. We use a bilateral filter to process MR images and to build a multiscale image series by increasing the standard deviation of spatial function and reducing the standard deviation of range function. We separate every scale image into multiple blocks and for every block a multiscale fuzzy C-means classification method is applied along the scales from the coarse to fine levels to overcome the effect of intensity inhomogeneity. The method is robust for noise MR images with intensity inhomogeneity because of its multiscale and multiblock bilateral filtering scheme. Our method was compared with the conventional FCM, a modified FCM (MFCM) and multiscale FCM (MsFCM) method on synthesized images, simulated brain MR images, and real MR images. The MsbFCM method achieved an overlap ratio of greater than 91% as validated by the ground truth even if original images have 9% noise and 40% intensity inhomogeneity. Experimental results using real MR images demonstrate the effectiveness of the proposed method. Our MsbFCM classification method is accurate and robust for various MR images.