The aim of this doctoral dissertation was to develop a method for quantitative rockfall susceptibility assessment at slope scale using techniques of remote sensing. For the development of this method 4 pilot areas were chosen, each with different structural-geological setting. These pilot areas are rock slopes above the town of Omiš, rock slopes and rock cuts Brljan in the Krka National Park, rock cuts Špičunak and Lokve in the Gorski kotar area. High-resolution and precision digital models of the investigated rock slopes were created based on terrestrial laser scanning and UAV based photogrammetry. Identification and mapping of the rock mass properties (degree of weathering and rock type) and the geometrical properties of the discontinuities (orientation and sets, spacing and trace length) was carried out by remote sensing on the digital models using manual, semi-automated and automated techniques and methods. Also, for the Lokve rock cut pilot area two different rock types were determined by semi-automated analysis of intensity data, derived from terrestrial laser scanner. Based on analysed rock mass and discontinuity data, pilot areas were divided into 65 quasi-homogenic engineering geological zones and each zone was associated with certain rock mass and discontinuity data that can be found in that zone. Engineering geological zones, along with rock mass and discontinuity properties, were used as input dana for quantitative rockfall susceptibility assessment. In order to quantify susceptibility, 4 susceptibility indices were developed: Plane Failure, Wedge Failure and Toppling Susceptibility Indices, from which Rockfall Susceptibility Index is derived based on the probability of union. Susceptibility indices were calculated for each 3D rock slope model cell based on the probability analysis, spatial kinematic analysis and modified Kinematic Hazard Index. This calculation processes were entirely automated by developed algorithm ROCKS in MATLAB programming platform. Rockfall susceptibility models were created by interpolation of the susceptibility indices on the digital rock slope models for each pilot area. In total 16 rockfall susceptibility models were created. Verification of the resulting rockfall susceptibility models has indicated that high values of Rockfall Susceptibility Index highlight potentially unstable rock blocks and unfavourably orientated discontinuities. Also, zones with higher density of high Rockfall Susceptibility Index values point out to certain quasi-homogenic engineering geological zones which should be priority in rockfall hazard and risk reduction by efficient design of rock slope protection measures. Developed method for quantitative rockfall susceptibility assessment has shown that it can be used in detailed engineering geological and geotechnical investigation of rock slopes, and rock slope remediation and protection assessments. Also, it represents a precise and valuable input data for rockfall hazard and risk assessment at the slope scale.