The physicochemical heterogeneity found on amorphous surfaces leads to a complex interaction of adsorbate molecules with topological and undercoordinated defects, which enhance the adsorption capacity and can participate in catalytic reactions. The identification and analysis of the adsorption structure observed on amorphous surfaces require novel tools that allow the segmentation of the surfaces into complex-shaped regions that contrast with the periodic patterns found on crystalline surfaces. We propose a Random Forest (RF) classifier that segments the surface into regions that can then be further analyzed and classified to reveal the dynamics of the interaction with the adsorbate. The RF segmentation is applied to the surface density map of the adsorbed molecules and employs multiple features (intensity, gradient, and the eigenvalues of the Hessian matrix) which are nonlocal and allow a better identification of the adsorption structures. The segmentation depends on a set of parameters that specify the training set and can be tailored to serve the specific purpose of the segmentation. Here, we consider an example in which we aim to separate highly heterogeneous regions from weakly heterogeneous regions. We demonstrate that the RF segmentation is able to separate the surface into a fully connected weakly heterogeneous region (whose behavior is somehow similar to crystalline surfaces and has an exponential distribution of the residence time) and a very heterogeneous region characterized by a complex residence-time distribution, which is generated by the undercoordinated defects and is responsible for the peculiar characteristics of the amorphous surface.