## 4 Results
### IVCI Model Results  
Initially, we resampled the resolution of the SBAS-InSAR results to 30 m \(\times\) 30 m, matching the grid size of the initial IV-CatBoost map. The results were then classified into five categories based on the natural breaks of annual deformation rates: \(<\)\(-\)14.88 mm/y, \(-\)14.88 mm/y to \(-\)9.14 mm/y, \(-\)9.14 mm/y to \(-\)5.39 mm/y, \(-\)5.39 mm/y to \(-\)2.46 mm/y, and \(>\)\(-\)2.46 mm/y. Each category was assigned a value of 5, 4, 3, 2, or 1. In order to ensure the accuracy of the model, no interpolation was performed in the blank areas. Only areas with annual subsidence rate results were refined. The data shown in Table 2 and landslide susceptibility were combined using ArcGIS to obtain the combined IVCI model. The updated results are shown in Figure 7.  
The distribution of susceptibility zones in the IVCI model is similar to the IV-CatBoost zoning. The IVCI model results (Table 4) show a reduction in the very low, low, and moderate susceptibility zones, while an increase is observed in the high and very high susceptibility zones. The differences in proportions between the zones are within 0.9%, indicating that the IV-CatBoost model's classification trend is correct, although some regions require further refinement. To visually observe the changes in the grid before and after refinement, the grids were subtracted before and after integration to facilitate the analysis of changes in different susceptibility levels. According to Table 5, after integration, there were 94,124 grids (84.71 km\({}^{2}\)) that had their susceptibility levels increased by one level, 46,992 grids (42.29 km\({}^{2}\)) increased by two levels, 19,535 grids (17.58 km\({}^{2}\)) increased by three levels, and 5133 grids (4.62 km\({}^{2}\)) increased by four levels. In total, 165,784 grids obtained a more accurate susceptibility classification under the correction of InSAR.  
Figure 9: Map of SBAS-InSAR annual subsidence rates.  
In Figure 10, there are three distinct areas where the susceptibility levels have undergone significant corrections. To provide a clearer explanation of the model's refinements and validate the accuracy of the IVCI model, these areas are labeled as Y-1, Y-2, and Y-3. Detailed discussions on these significant correction areas are presented in Section 5.1.  
\begin{table}
\begin{tabular}{c c c c c c} \hline \hline \multicolumn{6}{c}{**Landslide**} \\ \multicolumn{1}{c}{**Susceptibility Class**} & **IV-CatBoost** & \multicolumn{2}{c}{**IVCI**} & \multicolumn{2}{c}{**Landslide Susceptibility Increase**} \\ \hline
**Class** & **No. Cells** & **\%** & **NO. Cells** & **\%** & **Class** & **No. Cells** \\ \hline VL & 4,146,823 & 21.63\% & 4,071,042 & 21.24\% & 0 & 19,005,155 \\ L & 3,763,380 & 19.63\% & 3,741,168 & 19.51\% & +1 & 94,124 \\ M & 3,343,630 & 17.44\% & 3,372,395 & 17.59\% & +2 & 46,992 \\ H & 3,358,220 & 17.52\% & 3,403,205 & 17.75\% & +3 & 19,535 \\ VH & 4,558,886 & 23.78\% & 4,583,129 & 23.91\% & +4 & 5133 \\ \hline \hline \end{tabular}
\end{table}
Table 5: Corrections of the partition differences before and after the model.  
Figure 10: Map of the post-correction level differences between the IV-CatBoost and IVCI model.  
\begin{table}
\begin{tabular}{c c c c c c} \hline \hline \multicolumn{6}{c}{**Landslide Susceptibility Classes**} \\ \cline{2-6}  & **VH** & **H** & **M** & **L** & **VL** \\ \hline Area of zone (km\({}^{2}\)) & 4102.33 & 3022.40 & 3009.27 & 3387.04 & 3732.14 \\ Proportion (\%) & 23.78\% & 17.52\% & 17.44\% & 19.63\% & 21.63\% \\ Landslide ratio (\%) & 74.16\% & 16.55\% & 6.02\% & 2.21\% & 1.06\% \\ \hline \hline \end{tabular}
\end{table}
Table 4: IV-CatBoost model areas of the landslide susceptibility classes. Very-high susceptibility zones (VH); high susceptibility zones (H); moderate susceptibility zones (M); low susceptibility zones (L); very-low susceptibility zones (VL).