We aimed to determine whether multiresolution fractal analysis of voxel-based dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) parametric maps can provide early prediction of breast cancer response to neoadjuvant chemotherapy (NACT). In total, 55 patients underwent 4 DCE-MRI examinations before, during, and after NACT. The shutter-speed model was used to analyze the DCE-MRI data and generate parametric maps within the tumor region of interest. The proposed multiresolution fractal method and the more conventional methods of single-resolution fractal, gray-level co-occurrence matrix, and run-length matrix were used to extract features from the parametric maps. Only the data obtained before and after the first NACT cycle were used to evaluate early prediction of response. With a training (N = 40) and testing (N = 15) data set, support vector machine was used to assess the predictive abilities of the features in classification of pathologic complete response versus non-pathologic complete response. Generally the multiresolution fractal features from individual maps and the concatenated features from all parametric maps showed better predictive performances than conventional features, with receiver operating curve area under the curve (AUC) values of 0.91 (all parameters) and 0.80 (K
trans
), in the training and testing sets, respectively. The differences in AUC were statistically significant (
P
< .05) for several parametric maps. Thus, multiresolution analysis that decomposes the texture at various spatial-frequency scales may more accurately capture changes in tumor vascular heterogeneity as measured by DCE-MRI, and therefore provide better early prediction of NACT response.