1. Introduction
PV cells work under the photoelectric effect: when a light photon with appropriate frequency affects the cell surface, the silicon atoms are excited and an electron is ripped off them. In these conditions, each cell is converted into an electric generator that, grouped with the other cells of the PV module, produces the voltage and power supply to the external circuit. A single PV cell will overheat when it is inversely polarized or its connector is damaged. The reason is that all the energy circulates through a small part of the cell causing a significant temperature increase. Thus, regardless of the cause of the damage in the cell, the damaged surface will suffer a local overheating effect. Thereby, heat transfer in solar energy conversion devices is one of the most important parameters to be studied for the determination of their state and their efficiency. In addition to other techniques mainly based on electrical properties [
1,
2], infrared thermography (IRT) has been proven as an adequate technique for the evaluation of solar energy installations [
3]. The reason for the application of IRT is its capability to measure superficial temperature of objects from the radiation emitted in the thermal infrared band of the spectrum. Particularly, the detection of relative temperature differences allows de determination of thermal gradients and the detection of anomalous behaviors. Several thermographic techniques have been applied for the detection of pathologies and the evaluation of the productivity of the energy modules. These tests are mainly performed in laboratories, with limited sample size and mainly focused on the detection and precise characterization of pathologies at the level of PV cell with the aim at improving its design and performance [
4,
5,
6]. In the cases where the thermographic studies were performed in-situ, two main aspects were analyzed. The first is the characterization and parameterization of the external factors that affect the absolute precision of the thermographic inspection for the performance of quantitative studies [
7]. The second is the evaluation of energy facilities on small-scale communities [
8].
Other pathologies that can reduce the productivity of the PV facilities are related to the geometric defects in the mounting. PV panel productivity depends directly on the incident solar radiation which is maximized when the solar rays strike orthogonally the photovoltaic surface. In this way, azimuth and tilt parameters of PV panels should be carefully defined according to the geographical location of the PV power station and validated after the construction process to confirm the fulfillment of the project guidelines.
The generation of 3D models from dense 3D point clouds is well-known in the scientific community [
9,
10,
11,
12]. These products present proven quality and applicability for the development of automatic segmentation and geolocation of elements of interest. The 3D models are generated through the segmentation and parameterization of the point clouds, using techniques that are usually based on initial pre-filtering strategies complemented with more sophisticated algorithms which take advantage of the geometric and radiometric information associated with each 3D point. Segmentation techniques can be focused on different strategies such as iterative search of basic primitives (planes, lines, etc.) [
13,
14,
15,
16,
17]; evaluation of different sets of characteristics calculated from a point and its neighborhood [
18,
19,
20]; or multiclass classification techniques based on supervised machine learning algorithms [
21,
22,
23]. In this way, it is possible to locate different elements regardless the complexity of the scenario going from simple geometries such as roofs [
24] or columns [
25] to complex geometries such as trees [
26,
27], buildings [
24,
28,
29] or vehicles [
30].
The current procedure for monitoring PV stations is based on the thermographic inspection, panel-by-panel, performed in-situ by a human operator who tries to estimate the percentage of damaged panel that can affect its productivity. This procedure is not only slow and expensive but also dependent on the subjectivity and expertise of the operator. Some companies use unmanned aerial vehicles equipped with thermographic video cameras in a similar way as proposed in this studies [
2,
31,
32], but the processing is performed manually through the visualization of the data recorded by a human operator. Thus, capture time is reduced but processing time and human subjectivity are still elements that affect the quality of the final product, although some researchers have tried to solve this subjectivity by using image processing algorithms [
32]. An alternative to the use of thermographic cameras is the integration of RGB cameras in the UAV [
33,
34]; the main drawback of this option is that the type of defect detected depends on the flight conditions.
Both current established inspection protocols and more novel methodologies proposed by researchers are still affected by the subjectivity implicit in the inspection processes supervised by human operators. Although these novel methodologies present advances in data acquisition, they do not propose automatisms that allow speeding up the inspection processes. In addition, these novel methodologies do not provide a metric product that allows a geometric evaluation of the facilities in the inspection or construction processes. In order to improve the limitations of the methodologies remarked above, we developed a methodology able to segment automatically individual PV panels from a precise dataset obtained with a low-cost aerial platform equipped with RGB cameras. The segmented individual panels are evaluated through the complementary thermographic information. Particularly, this methodology provides advances in the robustness and efficiency of the classical inspection and validation procedures, which performed by a human operator. Results are obtained through unsupervised detection techniques adding semantics to the geometry contained in the point cloud. This semantic information is accurately geo-referenced reaching centimeter accuracy. Thus, the inspection procedure becomes accurate, faster and not affected by human subjectivity.
The proposed methodology consists of three complementary segmentation strategies: (i) a geometrically-based segmentation process is used to extract points belonging to PV panel clusters; (ii) an intensity-based segmentation [
35,
36,
37] is used to segment each individual PV panel through the detection of points belonging to their frames; and (iii) a thermographic-based segmentation is used for the detection and location of pathologies through a qualitative approach based on relative temperature differences. Thus, the computation of the thermophysical properties of PV modules or their precise absolute temperature is not necessary and the accurate location of the pathologies is performed through a custom 5D metric product. In this case, images from the RGB cameras are processed with photogrammetric and computer vision algorithms [
38,
39] for their orientation and generation of a 3D dense point cloud of the PV panels.
The data processing methodology is scalable to any other high-density aerial data sources such as LiDAR (Light Detection And Ranging) [
40]. The use of geo-referenced 3D dense point clouds will allow the performance of accurate analysis of areas, azimuths and tilts of the PV panels which are parameters strictly linked to the productivity of PV installations [
41,
42,
43].
The paper has been structured as follows: after this introduction,
Section 2 includes a detailed explanation of the materials and methods used for data acquisition and processing towards the automatic segmentation of PV panels and their classification according to the existence of geometric or thermal irregularities.
Section 3 is devoted to analyzing the methodology through the results of its application to a PV power station selected as case study. Finally,
Section 4 establishes the most relevant conclusions of the approach. The selected case study is a PV power station located in Gotarrendura (center of Spain) (coordinates 40°42′N, 4°44′W), which is nationally well-known as pioneer in the implementation of new energy models based on the exploitation of renewable energy sources and on self-sufficiency.
3. Experimental Results
3.1. Study Case
The proposed methodology has been selected due to the knowledge of the existence of damaged PV panels, adequate for testing the methodology. The PV installations consist of 21 fixed PV panel clusters distributed in 8 equidistant lines. Each line is composed by 30 PV panels of the same type with dimensions of 1.5 m × 0.8 m. Temporal security regulations limited the inspection to the first 6 lines from south to north. The result was an aerial survey of 16 PV clusters with an extension of 4000 m2 in a rectangular shape of 50 m × 80 m.
Regarding the RGB Flight planning (
Figure 7a), time between shots was established as 2 s for an approximated flight speed of 10 km/h, ensuring image acquisition with minimum forward and side overlaps of 70% and 30%, respectively. The focal length of the RGB sensor was fixed to 14 mm and the flight altitude over the ground was 40 m, resulting on a GSD of 1 cm. As a result, 43 RGB images were acquired covering the whole study area.
The full resolution RGB images were processed according to the photogrammetric and computer vision methodologies. RGB radiometry of each point was converted to its intensity value on the visible spectrum. The point cloud generated was preprocessed with the voxel grid algorithm to homogenize point density, obtaining as result a point cloud of 39,194,155 points, implying an approximated average resolution of 10,000 points/m2 (≈1 point/GSD2).
In order to georeference the scene and perform the quality assessment of the survey, a topographic GNSS (Global Navigation Satellite System) survey was performed, consisting of nine GCP homogeneously distributed over the study area (
Figure 7b). Because the absolute positioning accuracy of the scene should not be strict for this type of work, GCP were measured using RTK (Real Time Kinematics) positioning, ensuring that on-flight positioning precision provided by the equipment was better than 2 cm. Five of these GCP (points 1, 3, 5, 7 and 9) were used for the absolute georeferencing of the scene (
Table 2) obtaining a final deviation of 0.012 m. The remaining points (points 2, 4, 6 and 8) were used as check points to perform the geometric quality assessment of the final product (
Table 3) obtaining a final deviation of 0.033 m.
In the case of the IRT flight planning, time between shots and flight speed should be carefully correlated to ensure a high redundancy in the information to allow the supervised selection to discard blurred or radiometrically inconsistent imagery. Experience derived from previous surveys allows establishing an approximated flight speed of 10 Km/h as an appropriate parameter for the sensor used. Flight planning was designed through the definition of way points at the beginning and end of each line of PV panel clusters, guaranteeing that the vehicle overhangs the panels in the vertical of the terrain. Thermographic redundancy is guaranteed due to the high frame rate of the sensor, established in 50 fps. Flight altitude selected was 10 m over the ground, resulting on a GSD of 2.5 cm. In this case, 116 thermographic images were acquired and supervised by a human operator to select a total amount of 80 for the study of the study area. The thermographic texture was projected over the point cloud obtaining a hybrid RGB-I 3D point cloud with thermographic texture mapped over the points of the PV panels.
The 5D point cloud with intensity and thermographic texture is the input of the algorithm developed. In the first step, ground points are removed using the progressive morphological filter, obtaining as result a point cloud of the PV panels with 2,524,058 points. In this first filtering process, only geometry (X, Y, Z) is required.
The point cloud of the PV panels without ground points is introduced into the geometric segmentation process resulting in the extraction of 16 PV panel clusters (
Figure 8a) that have been automatically evaluated and classified regarding their geometric properties (azimuth and tilt) (
Table 4) without finding any construction defect limiting the productivity of the installation. These clusters are processed with the intensity-based segmentation process, resulting in the segmentation of 480 PV panels (
Figure 8b).
Each individual PV panel is evaluated through temperature radiometry and classified regarding the existence of pathologies or damaged areas, quantifying the percentage of area affected (
Figure 9). In total, nine damaged PV panels were detected (
Table 5).
The results are stored using both simple ASCII files and GIS Shapefiles, allowing the simple integration of the information derived from the proposed methodology with external data sources (
Figure 10).
3.2. Accuracy Assessment
In order to perform an accuracy assessment of the proposed methodology, fieldwork has been carried out to acquire validation information that can be established as ground truth. Regarding the PV panel segmentation, the validation has consisted in a supervised monitoring of the segmentation results. This process determines the high robustness of the segmentation process, with no lack or excess of segmentation: only the existing PV panels were segmented, with no appearance of new panels resulting from over-segmentation. The ground truth for the thermal pathologies was defined through a visual inspection. The whole case study was evaluated from the ground by an expert human operator aided with a thermographic camera. The human operator took note of all detected thermal pathologies, referencing the damaged PV panels through a predefined reference system based on their position within the facilities (line, cluster and number of panel). The proposed methodology was also able to detect all pathologies detected by the human operator in the field survey, without creating new (non-existent) ones. Regarding damaged areas of each PV panel, these were also measured by an expert operator following two different methodologies. The first, through a visual estimation of the damaged area (in percentage). This subjective approach is the usual methodology according to the established inspection protocols. The second, designed to obtain a ground truth less affected by the subjectivity of the operator and which will be used for comparison with the proposed methodology, was generated by measuring the damaged surfaces using a flexometer. Given the arbitrariness in the shape of pathologies, the areas of the thermal pathologies were calculated from their approximation to basic primitives such as rectangles (height and width measurements) and circumferences (diameter measurements), which are assumed to result in over-sized area values. The comparison between the results and the ground truth obtained (
Table 6) manifests the effects of the subjectivity on the visual inspection, showing the tendency of the inspector to quantify the damaged surface by counting only the number of cells fully affected by overheating. This results in an under-sizing of the actual damaged surface because the productivity reduction by the overheating of adjacent cells is not taken into account. The comparison of the results of the proposed methodology with the ground truth established from the direct measurement of damage dimensions (with the flexometer) shows the accuracy of the proposed workflow with discrepancies lower than 2% of the panel area (1.22 m
2), corresponding to discrepancies smaller than 0.02 m
2 for the 480 solar panels evaluated.
Finally, regarding the geometric evaluation of the facilities, the accuracy of the azimuth and tilt of the PV clusters was tested through the comparison of the parameters extracted with the proposed methodology and the values defined as ground truth by the topographic survey (
Table 7). In particular, the topographic survey was supported by the GCP and using a reflectorless total station, Trimble VX, whose precision for 3D point surveying without prism is 10 mm. Solar panels azimuth and tilt ground truth values were acquired making multiple angular and distance measurements in order to guarantee more precision and reliability. Discrepancies were lower than one degree for angular parameters, being these results completely acceptable due to the precisions established for the measurement of the GCPs used both for the absolute orientation of the point cloud and for the establishment of the reference system in the topographic survey.
3.3. Computing Efficiency Analysis
The machine used to perform the computation was a Microsoft Windows 8.1 workstation with 32 GB RAM, a 3.40 GHz Intel Core i5-3570K processor and an Nvidia Quadro 2000 GPU. The dense point cloud generation was the most computationally demanding process investing a total of 17 min and 45 s for the whole data set.
Regarding the algorithm developed for the automatic classification of segmented PV panels, the execution time performed through a single thread algorithm was 2 min and 23 s for the whole study area. However, multithreading computation has been implemented, reducing computation time proportionally to the number of concurrent running threads.
4. Discussion
The advantages of using the proposed methodology for the evaluation of the state of PV panels from PV power stations include both efficiency and robustness against traditional thermographic inspection methods.
Whilst, in this study, the required data were acquired with an UAV, the proposed methodology is valid for processing datasets captured with different airborne sensors such as LiDAR or any other RGB and thermographic sensors transported by any manned or unmanned aerial vehicle that guarantees the required resolution and precision.
The quality of the product obtained with the proposed methodology relies on the factors mentioned below. The main problem is the complexity of the automation in the detection of pathologies minimizing the appearance of false positives, mainly due to the high number of factors to consider in the thermal imaging interpretation, such as environmental conditions, solar reflections, surface state of the element under study and emissivity. This problem can be solved through the election of an appropriate survey time window and the follow-up of the guidelines specified in
Section 2.2.1 at the time of the thermographic acquisition. Geolocation accuracy is dependent on the quality of the ground control points used, in such way that the use of equipment and survey processes is necessary to guarantee the precision of the project.
The proposed methodology was able to detect 100% of the pathologies established in the ground truth, justifying the accuracy of the methodology for maintenance tasks. However, it is necessary to highlight that the use of the median as an estimator of the reference temperature value for non-damaged surfaces would not be a valid definition if the damages in a panel cluster cover more than one half of the surface. If this condition is fulfilled, the proposed methodology for the detection of thermal anomalies on the photovoltaic surface would set the temperature of damaged surfaces as the reference temperature value and therefore the entire surface would be classified as undamaged. However, if this happened, the safety systems of the PV power station would generate warnings produced both by the existence of very high temperatures and by the low productivity of the system. Therefore, this is not a typical maintenance task and the proposed methodology is designed to avoid this situation, allowing the detection and resolution of incidences at an earlier stage.
Maximum discrepancies in the measurement of damaged surfaces descend from 5% for visual inspection to 2% using the proposed methodology. The geometric evaluation of the facilities presents discrepancies regarding the ground truth lower to one degree for angular parameters (azimuth and tilt). Although these accuracy improvements in the measurement of damaged areas cannot be considered as a decisive criterion for the replacement of current evaluation protocols for the proposed methodology, improvements in processing time and the metric qualities of the product are decisive to make more efficient inspection processes, especially in large PV power stations. Regarding land-based inspection methods, where the facilities are traversed by a human operator equipped with a thermographic camera and a GPS tracking system, our inspection method drastically reduces survey time reaching the day productivity of a human operator in just a few minutes of flight. In those cases where aerial vehicles are already used for the thermographic survey, the proposed unsupervised algorithm removes the tedious and manual process of thermographic interpretation and evaluation by a human operator, removing human errors and reporting metric and objective information from the evaluated facilities. In all cases, geolocation is highly improved reaching centimeter accuracy in the geolocation of pathologies. Thus, the inspection procedure becomes accurate, faster and not affected by human subjectivity.
Trying to compare the proposed methodology against other similar techniques, no similar strategies for the evaluation of large-scale PV power stations were found.
5. Conclusions
This paper presents a methodology for the automatic processing of 5D point clouds with geometric, intensity and thermographic information for the evaluation of the state of PV power stations through a completely automatic process, which allows the detection and evaluation of pathologies on the facilities. A combined geometric, radiometric (RGB intensity-based) and statistical process allows the identification of individual PV panels and the evaluation of construction inaccuracies, which could considerably limit the productivity of the facility by not taking full advantage of available solar resources. Individual clustered PV panels are evaluated again through a combined geometric, radiometric (temperature-based) and statistical approach, to detect and quantify surfaces affected by pathologies that can drastically reduce the production capacity of the element. The resulting hybrid product provides a complete and georeferenced, thermographic and metric information of the elements in the PV power station. This information enables faster and better detection of pathologies, their spatial location and interpretation of the real state of the facilities. The high-level of automation of the procedure allows the evaluation of wide PV power stations in a fast and accurate way, establishing as only limitations the capacity of payload and autonomy of the aerial platform as well as the operational capacities of the computer equipment where the process is executed.
This project opens new trends for future work both from a sensorial and methodological point of view. Concerning the first, the accurate registration of the thermographic sensor regarding the aerial platform, through the accurate relative orientation and timestamp synchronization between imaging sensors and navigation system, would allow the automatic registration between thermographic images and 3D point clouds. Robust integration of high density LiDAR sensors on UAV platforms is still a developing field. The use of these sensors will not only reduce processing time in the photogrammetric generation of 3D point clouds avoiding the difficulties related to the reconstruction of texture-less surfaces but will also provide a new radiometric value to integrate in the process, as the intensity of the returned signal is usually at near infrared wavelengths.
From a methodological point of view, regarding the possibility to take advantage of the large amount of information derived from this methodology, it would be interesting to advance in the definition of an universal data structure allowing to work with lighter and more agile information that we can include as an additional layer on a GIS (Geographic Information System) enabling the integration of this information with other data sources (cadastral information, demographic information, urban parameters, etc.) to perform energy demand studies allowing the proper dimensioning of the installations.