Multiscale analysis technique can describe an image from different resolutions and has been widely used for extracting features and modelling of remotely sensed images. The subsampled wavelet transforms are commonly employed for establishing the multiscale representation of an image. However, the wavelet-based features cannot describe patterns with long spatial spans and often result in noisy classification results. By contrast, the object-based image analysis can create classification maps composed of compact land objects. However, features extracted from a single scale still cannot provide discriminating information for the land cover classification. To improve the classification accuracy and alleviate noisy thematic maps, a regional multiscale classification method is proposed.The proposed method consists of two main blocks, including establishing a regional multiscale image representation and classifying on the basis of Markov Random Field (MRF). In the first block, the mean shift segmentation method is employed to create initial over-segmented regions. Thereafter, a rule combining the grey values in the regions and the shared boundary lengths among regions is designed to extract the low-frequency part of an image. The current high-frequency part is obtained by subtracting the original image and the current low-frequency portion. A regional multiscale representation can be iteratively established by replacing the original image with the low-frequency part and repeating the segmenting and decomposing process. In the second block, the classification result obtained from the original image is considered the prior of the label field in the first decomposed level, and the high-frequency part in the first level models the feature field of MRF. The classification result in the first level is obtained by solving an objective function consisting of feature and label energies. By iteratively projecting the current classification result to the next level and modelling the feature field with the high-frequency part, the final classification map can be acquired in the coarsest scale.To examine the effectiveness of the proposed method, three groups (six images in total), such as Prague textures, synthetic remote sensing textures, and QuickBird multispectral images, are adopted. The proposed regional multiscale MRF (RMRF) is compared with the Iterated Conditional Model (ICM), Graph Cut (GC), and Wavelet-based MRF classification methods (WMRF). Visual inspection and quantitative measurement are employed for the comparison. As GC fails to integrate the spatial information among neighbouring pixels, the classification results are dominated by the pixel spectral values. Textures or land cover patterns with complex spectral heterogeneity cannot be properly captured. WMRF provides a pixel-level multiscale representation, which is also insufficient for describing texture patterns across a large spatial span. ICM, GC and WMRF may create noisy thematic maps with different extensions when the spectral heterogeneity of a given land cover pattern is extreme. On the contrary, RMRF has a flexible framework to extract and model spatial information. RMF can capture the essential features of land objects with large spectral heterogeneity, resulting in maps with less noise.This study proposes an RMRF for classification. Experiments demonstrate that the established multiscale representation can efficiently describe texture patterns with wide spatial spans, such as land objects with complex structures. By combining the multiscale representation with MRF model, RMRF can achieve a high-precision semantic segmentation of ground objects. The adaptive estimation of the decomposition layer number and parameters of the algorithm are the on-going works.