An estimate of solar irradiation potential over large regions requires the knowledge of the long-term spatiotemporal distribution of the solar radiation as well as the identification of the suitable surfaces where the photovoltaic (PV) installations can be built. These main components can be modelled in different ways and are thus affected by different sources of uncertainty. Thus, when estimating the exploitable potential over large regions, it is important to measure the accuracy of the entire process. In this work, we provide a generic method to estimate the solar irradiation potential of rooftops over large regions and an estimate of the corresponding uncertainties when calculating the long-term electricity generation of PV plants. This method uses satellite based solar radiation data covering a period of 22 years, with a temporal resolution of 30 min and a spatial resolution of 3.8-5.6 km. Suitable surfaces on rooftops are identified using Digital Surface Models combined with building footprints. This allows to determine the geometry of rooftops, such as slope, and orientation with a spatial resolution of 0.5 m. Finally, we calculated the electricity generation based on models which take into account all characteristics of PV system components. In order to estimate the accuracy of the model for electricity production, we compared the monthly generation of 500 PV plants in Switzerland consisting of different PV technologies (mono-crystalline, poly-crystalline and thin film) with the estimates. The validation results show a correlation coefficient (R-2) of 0.9 and a median monthly relative error between 0.28% (August) and 28.08% (December). The monthly estimates are more accurate during summer time, while spatially and technology-wise no significant differences are found.