Description: In remote sensing change detection, Markov random field (MRF) has been used successfully to model the prior probability using class-labels dependencies. MRF has played an important role in the detection of complex urban changes using optical images. However, the preservation of details in urban change analysis turns out to be a highly complex task if multitemporal SAR images with their speckle are to be used. Here, the ability of MRF to preserve geometric details and to combat speckle effect at the same time becomes questionable. Blob-region phenomenon and fine structures removal are common consequences of the application of traditional MRF-based change detection algorithm. To overcome these limitations, the iterated conditional modes (ICM) framework for the optimization of the maximum a posteriori (MAP-MRF) criterion function is extended to include a nonlocal probability maximization step. This probability model, which characterizes the relationship between pixels’ class-labels in a nonlocal scale, has the potential to preserve spatial details and to reduce speckle effects. Two multitemporal SAR datasets were used to assess the proposed algorithm. Experimental results using three density functions [i.e., the log normal (LN), generalized Gaussian (GG), and normal distributions (ND)] have demonstrated the efficiency of the proposed approach in terms of detail preservation and noise suppression. Compared with the traditional MRF algorithm, the proposed approach proved to be less-sensitive to the value of the contextual parameter and the chosen density function. The proposed approach has also shown less sensitivity to the quality of the initial change map when compared with the ICM algorithm.