Given a co-location pattern consisting of spatial features, the prevalent region mining process identifies local areas in which these features are co-located with a high probability. Many approaches have been proposed for co-location mining due to its key role in public safety, social-economic development and environmental management. However, traditionally, most of the solutions focus on itemsets mining and results outputting in a textual format, which fail to adequately treat all the spatial nature of the underlying entities and processes. In this paper, we propose a new co-location analysis approach to find the prevalent regions of a pattern. The approach combines kernel density estimation and polygons clustering techniques to specifically consider the correlation, heterogeneity and contextual information existing within complex spatial interactions. A kernel density estimation surface is created for each feature and subsequently the generated multiple surfaces are combined into a final surface with cell attribute representing the pattern prevalence measure value. Polygons consisting of cells are then extracted according to the predefined threshold. Through adding appended environmental data to the polygons, an outcome of similar groups is achieved using polygons clustering approach. The effectiveness of our approach is evaluated using Points-of-Interest datasets in Shenzhen, China.