Slope stability in rock masses can be controlled by weak zones associated with discontinuities, fracturing or joints. Therefore, it is necessary to establish the number of joint sets,their orientation, opening, persistence, spacing and intensity. Since fracture intensity on rock slope surfaces can be complex to estimate, conventional methods of analysis assume homogeneous distributions by setting a constant value determined from compass and tape measurements in some sampled sectors of the slope. To account for the non-homogeneous nature of the rock mass on the entire slope surface, we propose to combine spatial statistics, data collection with short-range photogrammetry and three-dimensional image analysis techniques to estimate and map fracture intensity. To do this, each joint set, defined as a set oflines, is represented by a spatial point pattern. Where the intensity function varies spatially,the non-parametric kernel is used to estimate fracture intensity and the results are mappedfor validation. The methodology is applied on two rock masses near Bogota´, Colombia, withpromising results for fracture intensity slope zoning. These findings are highly relevant for application in road infrastructure and mining projects.