Liver cirrhosis, as the terminal phase of chronic liver disease fibrosis, is associated with high morbidity and mortality. Traditional methods for assessing liver function, such as clinical scoring systems, offer only a global evaluation and may not accurately reflect regional liver function variations. This study aimed at evaluating the diagnostic potential of whole-liver histogram analysis of gadobenate dimeglumine (Gd-BOPTA)-enhanced magnetic resonance imaging (MRI) for predicting the progression of cirrhosis. In this retrospective study, 265 consecutive patients with cirrhosis admitted to the Department of Radiology, Shandong Provincial Hospital Affiliated to Shandong First Medical University from August 2012 to September 2019 were enrolled. After the exclusion criteria were applied, 117 patients (84 males and 33 females) were divided into Child-Pugh A cirrhosis (n=43), Child-Pugh B cirrhosis (n=49), and Child-Pugh C cirrhosis (n=25). After correction for liver signal intensity with the spleen was completed, 19 histogram features of the whole liver were extracted and modeled to evaluate liver function, with the Child-Pugh class being incorporated as a clinical parameter. Receiver operating characteristic (ROC) curves were used to assess the diagnosis capability and determine the optimal cutoffs after a mean follow-up of 42.3±19.1 (range, 8-93) months. The association between significant histogram features and the cumulative incidence of hepatic insufficiency was analyzed with the adjusted Kaplan-Meier curve model. Among 117 patients (12%), 14 developed hepatic insufficiency through a period of follow-up. Five features, including the median (P<0.01), 90th percentile (P<0.01), root mean squared (P<0.01), mean (P<0.01), and 10th percentile (P<0.05), were significantly different between the groups with and without hepatic insufficiency according to the Kruskal-Wallis test; in the ROC curve analysis, the area under the curve (AUC) of these features was 0.723 [95% confidence interval (CI): 0.653-0.793], 0.722 (95% CI: 0.652-0.792), 0.722 (95% CI: 0.652-0.792), 0.721 (95% CI: 0.651-0.791), and 0.674 (95% CI: 0.600-0.748) after correction, respectively (all P values <0.05). Median, 90th percentile, root mean squared, and mean were found to be significant factors in predicting liver insufficiency. The adjusted Kaplan-Meier curves revealed that patients with a feature level less than the cutoff, as compared to those with a level above the cutoff, showed a statistically shorter progression-free survival and higher incidences of hepatic insufficiency for significant features of median (cutoff =26.001; 21.28% versus 5.71%; P=0.02), 90th percentile (cutoff =86.263; 20.41% versus 5.88%; P<0.01), root mean squared (cutoff =1,028.477; 19.15% versus 7.14%; P=0.049), and mean (cutoff =27.484; 19.15% versus 7.14%; P=0.049). Patients with a 10th percentile less than -39.811 also showed a higher cumulative incidence of hepatic insufficiency than did those with a value higher than the cutoff (0.18% versus 7.46%; P=0.22). Whole-liver histogram analysis of Gd-BOPTA-enhanced MRI may serve as a noninvasive analytical method to predict hepatic insufficiency in patients with cirrhosis.