Quality Assessment of the Bidirectional Reﬂectance Distribution Function for NIR Imagery Sequences from UAV

: Imaging from low altitudes is nowadays commonly used in remote sensing and photogrammetry. More and more often, in addition to acquiring images in the visible range, images in other spectral ranges, e.g., near infrared (NIR), are also recorded. During low-altitude photogrammetric studies, small-format images of large coverage along and across the ﬂight route are acquired that provide information about the imaged objects. The novelty presented in this research is the use of the modiﬁed method of the dark-object subtraction technique correction with a modiﬁed Walthall’s model for correction of images obtained from a low altitude. The basic versions of these models have often been used to radiometric correction of satellite imagery and classic aerial images. However, with the increasing popularity of imaging from low altitude (in particular in the NIR range), it has also become necessary to perform radiometric correction for this type of images. The radiometric correction of images acquired from low altitudes is important from the point of view of eliminating disturbances which might reduce the capabilities of image interpretation. The radiometric correction of images acquired from low altitudes should take into account the inﬂuence of the atmosphere but also the geometry of illumination, which is described by the bidirectional reﬂectance distribution function (BRDF). This paper presents a method of radiometric correction for unmanned aerial vehicle (UAV) NIR images. The study presents a method of low-altitude image acquisition and a fusion of the method of the dark-object subtraction technique correction with a modiﬁed Walthall’s model. The proposed solution performs the radiometric correction of images acquired in the NIR range with the root mean square error (RMSE) value not exceeding 10% with respect to the original images. The obtained results conﬁrm that the proposed method will provide effective compensation of radiometric disturbances in UAV images.


Introduction
The determination of the radiometric quality of images acquired from low altitudes using the bidirectional reflectance distribution function (BRDF) is now a challenge in remote sensing and digital-image processing. Contemporary in situ measurements of the BRDF are costly and ineffective. So far, radiometric corrections were performed mainly on satellite-imagery data and classical aerial images. It is one of the most important issues in modern photogrammetry and remote sensing, especially in its many practical applications related to the natural environment. Many studies are often based on imagery data from different periods of time. Due to the impact of time and possible differences in seasons during which imagery was acquired, the objects within the scene may differ greatly in terms of pixel brightness. So far, many models of radiometric correction have been developed, but are applicable only to satellite imagery. In addition, these models are very complex and allow only to achieve approximate results. Therefore, it is necessary to carry out research on the suitability of modified models of radiometric correction, which until now were used to correct satellite imagery and classic aerial images, to correct images obtained from a low altitude. The main difference between satellite and aerial images, and low altitude images is the quality of the imaging sensor used and the imaging altitude. The application of UAVs for the acquisition of low-altitude images for remote sensing is already common. However, the radiometric quality of images acquired from low altitudes or the impact of the BRDF on UAV image-brightness distribution has not been considered in a broader sense, until now. So far, it has been mentioned in References [1,2]. Moreover, the radiometric and geometric calibration of cameras installed on UAVs often requires accurate information about the BRDF of ground reference targets, which may be natural or artificial [3,4]. In a different approach, the methods of radiometric correction of images may be grouped into empirical or semiempirical corrections. The methods provide a relative correction between images based on the DN (pixel value of image) information of the images. The empirical corrections are based on a statistical approach using mainly equalizing histograms. The other group of semiempirical corrections takes into account the most important radiometric effects like the BRDF or models of the atmosphere [5].
So far, the application of general models for correcting the brightness of objects in the entire image simultaneously, often resulted in unwanted changes of brightness of small objects in the image [6,7]. Usually, compact digital cameras that acquire images in the RGB and NIR ranges are installed on board low-cost UAVs. They are usually equipped with a constant focal length lens of a conical field of view that results in variable angle of observation in the image plane [8]. Usually, a classical photogrammetric block of images acquired from low altitudes has both a longitudinal and latitudinal overlay of 75%. It may cause some objects to be imaged in even several images. The same object imaged in the specific images will have various values of brightness, especially in the NIR-range images. Moreover, in such images the distribution of the BRDF will, to a large extent, depend on atmospheric conditions, like cloud cover, haze, or mist, and the angle of incident solar radiance. The radiometric correction of images acquired from low altitudes should be a basic process before the application of proper remote-sensing processing, especially in cases when the task is to acquire such data, e.g., the reflectance coefficients necessary for the radiometric correction of images. So far, the radiometric and geometric calibration of the sensor in laboratory conditions was one of the most important phases of the procedure. Nowadays, more and more research is conducted in order to apply methods of self-calibration or calibration in situ with the simultaneous acquisition of imagery data and orienting the images, DTM/DSM processing, radiometric correction of images, and final products in the form of orthophoto mosaic [1].
Research concerning the determination of the BRDF for images used in remote sensing has already been conducted for several years. The researchers observed in their experiments [9][10][11][12] that the BRDF model of an object might be determined with reasonable accuracy using only viewing zenith and viewing azimuth angles as inputs.
It is thus important to develop a model of radiometric correction of images acquired from low altitudes that will take into account the impact of the so-called low atmosphere and the angle of incidence of solar radiation during the acquisition of the imagery data. The conducted research proposes an approach that takes into account the impact of the haze effect and the geometry of illumination described by the BRDF.

Characteristic Radiometric Correction with the BRDF Function
When considering a beam of light incident on a surface and its intensity, it may be stated that it depends on the properties of the surface, the location, and type of the light source as well as the conditions of observation. In case of anisotropic surfaces, a change of any of these parameters results in a change of the surface luminance [13,14]. The characteristics of the light reflection for a specific surface are described by BDRF defined by Nicodemus [15]. The function may be defined for a specific wavelength as the ratio of radiance dLs of the light reflected from a surface element ds in the direction to the irradiance dEi of the surface element from the direction, as in the figure below ( Figure 1) [16]: Figure 1. Geometry of the light beams as in the bidirectional reflectance distribution function (BRDF) definition [16].
In practice, the intensity Ps is determined of the reflected light from a discrete detector in the direction of the surface element S encompassing radiation from a solid angle Ω. Then, the BRDF may be determined using the following approximation [16]: where: P i -the power of the radiation incident on the studied surface area. In photogrammetry and remote sensing, the BRDF describes the function of incidence and reflection of the light from a terrain surface by means of a physical model. Moreover, the BRDF plays an important role in the radiometric correction of satellite images. In the case of classic aerial images and images acquired from low altitudes, corrections using this function are still being studied. However, new methods based on the radiometric correction of UAV images have recently been proposed [1,17]. Low-altitude imaging is more and more often done in various weather and lighting conditions. Therefore, radiometric correction requires the development of new methods. In their work they proposed that, for low-altitude images, the central radiation elements incident into the camera sensor are the surface of reflected sunlight (L su (λ)) and surface diffuse radiance (L su (λ)) [17][18][19]: where ρ(λ, θ i , ϕ i , θ r , ϕ r ) is the BRDF and ρ(λ, 2π, θ r , ϕ r ) is the reflectance distribution function for diffuse light, τ v (λ) and τ s (λ) are the atmospheric transmittance in the view and solar paths, respectively, E 0 (λ) is the spectral irradiance on top of the atmosphere, E d (λ) is the spectral irradiance on the surface due to diffuse illumination, and θ is the solar incidence angle on a surface. θ i and θ r are the illumination and reflected light zenith angles, and ϕ i and ϕ r are the azimuth angles, respectively [17]. Another aspect of research concerning radiometric correction is the phenomenon that, for UAV images, very often the main source of heterogeneousness is the effect of a "hot spot"-part of the image that seems to be brighter when the direction of image observation gets close to the direction of sunlight [1,20]. This effect may be corrected using the BRDF model with the assumption that at the same time the ground reflectance and atmospheric effects coefficients are anisotropic. Very often, in the case of low-altitude images the reflection of sunlight from the imaged objects is non-Lambertian and it enhances the "hot-spot" effect. The BRDF models may be grouped in two sets of types: physical and empirical. In the case of low-altitude images, the empirical model is more practical since the model parameters may be calculated for it.

Empirical Models of Radiometric Correction
In the case of a lack of data from external measurements, the impact of the atmosphere on the quality of low-altitude images may be determined by means of statistical methods based on the DN value for each channel of the image. The empirical models include [21]: Simple Dark Pixel Subtraction:-the model takes into account the shift of histograms in the specific channels depending on the angle of image acquisition [22]. It is caused by reflected light from outside the field of view of the sensor, which reaches the field of view of the sensor even if the ground reflection coefficient is equal to zero.

Modified Chavez Method:
In some cases when the image content was not homogeneous but it contained urban areas, forests or flat areas, excessive correction was observed in the RED and NIR channels [23]. In such cases, the λ-κ rule was proposed for the atmospherically scattered radiance. The value of the κ coefficient was in the range from 0.5 for haze in the atmosphere to 4.0 for a clear Rayleigh-type atmosphere. Due to the fact that the blue band shift represents the greatest impact of the atmosphere, it may be expected that the values determined for the channel will be the most accurate ones. The determined value of the shift (calibrated) enables the determination of the value of κ. The greater the shift in the DN values of the pixels, the thicker the haze of the atmosphere. Moreover, the κ-coefficient rule depends on the altitude of flight. In the case of images acquired from low altitudes, the useful range of altitudes will usually be within 100 to 300 m. The values of the shift for other channels including NIR are presented in Reference [24].
Modified Walthall Method: Numerous physical and empirical models of the BRDF have been developed in a few dozen years. However, for the purpose of the study of low-altitude images, it is sufficient to apply modified empirical models. A model developed using the Walthall method [25] is one of such models. Walthall proposed in his solution a linear semiempirical model that took into account the view zenith and the relative azimuth. Additionally, the model was improved by Nilson and Kuusk, who took into account the solar-zenith dependence [20,26,27]. The model may be applied to each channel separately. The correction may be calculated for all pixels: where: ρ-modeled reflectance θ i -incident illumination zenith angle θ r -reflection view zenith angle ϕ-relative azimuth angle a, b, c, d-free parameters

Proposed Methodology of Radiometric Correction of UAV NIR Images
Having in mind the above described premises, the application of dark-object subtraction technique (DOS) was proposed for the first phase of image correction. The process consists in subtracting a defined value of DN of reflection from dark objects. The method assumes a simplified model of radiation scattering in the atmosphere. For images acquired from UAVs, it is assumed that there is no scattered radiation reflected from the object. The reflection depends only on the radiation directly reflected from the object and radiation scattered in the atmosphere that was not reflected from the Earth's surface ( Figure 2). In the next phase, a modified version of Walthall's empirical model of the BRDF is proposed as a model adjusted to the needs of the correction of images from UAVs. It is assumed that in the case of low-altitude images, the flight time is usually up to 40 minutes; therefore, the illumination direction may be assumed as constant for the entire block. The adapted Walthall model for radiometric correction then has the form: In this case, the recent value of the constant d determined empirically is dependent on the channel for which it is applied. The proposed method of radiometric correction of images is applied independently for each channel using the empirically determined model components.

Test Area
The test area was located in the vicinity of Tylicz ( directly reflected from the object and radiation scattered in the atmosphere that was not reflected from the Earth's surface ( Figure 2). In the next phase, a modified version of Walthall's empirical model of the BRDF is proposed as a model adjusted to the needs of the correction of images from UAVs. It is assumed that in the case of low-altitude images, the flight time is usually up to 40 minutes; therefore, the illumination direction may be assumed as constant for the entire block. The adapted Walthall model for radiometric correction then has the form: In this case, the recent value of the constant d determined empirically is dependent on the channel for which it is applied. The proposed method of radiometric correction of images is applied independently for each channel using the empirically determined model components.

Test Area
The test area was located in the vicinity of Tylicz (49⁰23′45′′, 21⁰01′25′′) (Poland) (Figure 3) at elevations of over 600 m above mean sea level indicating mountainous terrain. The area's altitude is much more differentiated (there are greater altitude differences) as compared with the quarry. The acquired images covered a fragment of a hill including a part of a ski lift (southern part) and an area of a town with dense settlement (northern part). The northern and southern parts are separated by the Muszynka river flowing across Tylicz. Along the river bed on the southern side there are forests. The area of the hill at the time of flight was covered with low The area's altitude is much more differentiated (there are greater altitude differences) as compared with the quarry. The acquired images covered a fragment of a hill including a part of a ski lift (southern part) and an area of a town with dense settlement (northern part). The northern and southern parts are separated by the Muszynka river flowing across Tylicz. Along the river bed on the southern side there are forests. The area of the hill at the time of flight was covered with low grass, while the area with the ski slope was covered with snow. The northern part of the study area is urbanized to a low degree. There are detached houses, roads and technical infrastructure, single trees, bushes, and mainly grass. The photogrammetric flights were done in March 2016 and February 2017.

Data Acquisition
The NIR imagery data were acquired by means of the Trimble UX5 platform. It is an unmanned aerial vehicle adapted for the acquisition of imagery data of high-terrain resolution. It is made mainly of light EPP foam and it is equipped with electric propulsion. UX5 is a fixed-wing UAV. It is capable of entirely autonomous flying at an indicated altitude with indicated longitudinal and latitudinal coverage of the images due to its onboard Global Positioning System/Inertial Navigation System (GPS/INS) system. The flight route may be monitored by means of the flight controller in real time. The UX5 platform enables automatic triggering of the camera shutter to acquire images of the Earth's surface. The aerial vehicle starts from a special mechanical launcher. The system may operate in wind up to 18 m/s and weather conditions not more severe than light rain. Sony NEX-5T (Sony, Tokyo, Japan) was mounted on board the Trimble UX5 for the acquisition of imagery data. It is a compact digital camera equipped with a complementary metal-oxide-semiconductor (CMOS) matrix of 16.1 megapixels (maximum image resolution of 4912 × 3264 pixels). The images were recorded in the JPEG format with dissipative data compression. The camera was equipped with Voigtlander lenses of a constant focal length of 15 mm and a maximum value of aperture of f/4.5. The NEX-5T camera was modified to acquire images in its full range of matrix sensitivity. The filter for cutting off radiation of wavelengths higher than 690 nm, originally installed in front of the camera matrix, was removed. The maximum wavelength recorded by the sensor is now about 1050 nm. Photogrammetric flights were conducted in the same test area in a year's interval in various weather and lighting conditions at altitudes of 150 m and 300 m. The images were acquired with the NIR-enabled camera with a black IR-only longpass filter (B + W 092 Tiffen/Wratten #89B, the Tiffen Company, Hauppauge, New York, NY, USA), which cut on at 695 nm. The application of the filter caused the imagery to contain NIR only, where the blue pixels (band 3) were again used to record nothing but pure NIR (roughly 800 nm-1050 nm), while the red band (band 1) in the imagery obtained with this filter was in fact the red edge, roughly from 690 nm through 770 nm ( Figure 4). All acquired images were recorded with radiometric resolution of 8 bits/channel [28][29][30].  The process of the radiometric correction of mosaics took into account developing DSM and orthoimages and then the application of DOS technique to the mosaics. The lowest values of DN for the specific images channels were subtracted at the DOS correction application. The operation was done for each channel separately. In order to determine the value of the dark object in the area of Tylicz, Regions of Interest (ROI) of 300 pixels were defined for both sets of test data. The ROIs were located in the area of the river (northern part of the mosaic) where the average values of DN were determined for the specific channels of the images. For the first data set (Tylicz 2016) for Band 1, the DN values were shifted by 23 digital counts, for Band 2 by 3 digital counts, and for Band 3 by 15 digital counts. For the second data set (Tylicz 2017), for Band 1 the DN values were shifted by 19 digital counts, for Band 2 by 2 digital counts, and for Band 3 by 14 digital counts.
The analysis of the images indicated that the greatest impact of the atmosphere is observed in the short waves range (Band 2). As a result, the image histogram after correction shows relocation where the pixels have nonzero values.
The analysis of the histograms (Figures 5-8) shows that the DOS method has a limitation resulting from the relocation of the histogram independently for each band. Figures 6 and 7 show that the greatest relocation is observed for the RED and NIR channels because the haze effect diminishes for longer wavelength channels. Therefore, it is probable that, in some special cases for the NIR channel, the determined value of the DN difference may be too big in comparison with the real one, especially in cases when the dark pixels are rare.    In the further phase of processing, it is proposed to take into account the BRDF radiometric correction based on the modified Walthall's model. The BRDF model parameters were determined for each channel-Band 1, Band 2, and Band 3-using the initial values of coefficients proposed by    In the further phase of processing, it is proposed to take into account the BRDF radiometric correction based on the modified Walthall's model. The BRDF model parameters were determined for each channel-Band 1, Band 2, and Band 3-using the initial values of coefficients proposed by    In the further phase of processing, it is proposed to take into account the BRDF radiometric correction based on the modified Walthall's model. The BRDF model parameters were determined for each channel-Band 1, Band 2, and Band 3-using the initial values of coefficients proposed by    In the further phase of processing, it is proposed to take into account the BRDF radiometric correction based on the modified Walthall's model. The BRDF model parameters were determined for each channel-Band 1, Band 2, and Band 3-using the initial values of coefficients proposed by In the further phase of processing, it is proposed to take into account the BRDF radiometric correction based on the modified Walthall's model. The BRDF model parameters were determined for each channel-Band 1, Band 2, and Band 3-using the initial values of coefficients proposed by Walthall. When determining them, the impact of the zenith view angle, the azimuth angle difference of the sun, and viewing directions, as well as the empirical value of reflectance at nadir view were taken into account for each of the data sets. Figures 9 and 10 present the results of the radiometric correction of orthophotomosaics. After the application of the radiometric correction, the orthophotomosaics are uniform with respect to radiometry and they take into account the impact of the zenith view angle, the azimuth angle difference of the sun, and viewing directions. Walthall. When determining them, the impact of the zenith view angle, the azimuth angle difference of the sun, and viewing directions, as well as the empirical value of reflectance at nadir view were taken into account for each of the data sets. Figures 9 and 10 present the results of the radiometric correction of orthophotomosaics. After the application of the radiometric correction, the orthophotomosaics are uniform with respect to radiometry and they take into account the impact of the zenith view angle, the azimuth angle difference of the sun, and viewing directions.

Results of the Image Quality Assessment
There are many approaches to image quality assessment presented in publications. There is, however, no universal procedure. Image quality assessment after radiometric correction is often done using the relative values of mean errors (R-RMSE) and correlation coefficient (CC) [32], correlation coefficients, standard deviation values, the differences in pixels brightness between images [33][34][35]. In the case of digital-images analysis, the relation may be expressed for each channel in the following form: where: Walthall. When determining them, the impact of the zenith view angle, the azimuth angle difference of the sun, and viewing directions, as well as the empirical value of reflectance at nadir view were taken into account for each of the data sets. Figures 9 and 10 present the results of the radiometric correction of orthophotomosaics. After the application of the radiometric correction, the orthophotomosaics are uniform with respect to radiometry and they take into account the impact of the zenith view angle, the azimuth angle difference of the sun, and viewing directions.

Results of the Image Quality Assessment
There are many approaches to image quality assessment presented in publications. There is, however, no universal procedure. Image quality assessment after radiometric correction is often done using the relative values of mean errors (R-RMSE) and correlation coefficient (CC) [32], correlation coefficients, standard deviation values, the differences in pixels brightness between images [33][34][35]. In the case of digital-images analysis, the relation may be expressed for each channel in the following form: where:

Results of the Image Quality Assessment
There are many approaches to image quality assessment presented in publications. There is, however, no universal procedure. Image quality assessment after radiometric correction is often done using the relative values of mean errors (R-RMSE) and correlation coefficient (CC) [32], correlation coefficients, standard deviation values, the differences in pixels brightness between images [33][34][35]. In the case of digital-images analysis, the relation may be expressed for each channel in the following form: where: DN-the pixel value i-the pixel number in band k k-number of the spectral band N-number of observations The relative measure of RMSE is applied to avoid large absolute values of deviations for large values of the DN.
The correlation coefficient is determined separately for each channel using the following form: where: The standard deviation (SD) describes the distribution of data around the mean value of the pixel in the image. Covariance is an un-normalized measure of the relationship between the analyzed data. Normalization is achieved by dividing the covariance by the product of standard deviations of the processed data. This processing results in the correlation coefficient. It is a normalized measure between the original image and the image after applying the DOS and BRDF corrections. The correlation coefficient values are within the range [−1, 1]. The higher the value of the correlation coefficient is, the more correlated the data are, and, therefore, the higher the similarity between the original image and the corrected one is.
The determined RMSE (Table 3) for both sets of data do not exceed 10% in all spectral channels. It means that the selected coefficients of the BRDF correction model are reasonable. They do not degrade the content of the image. The smallest value of the error was obtained for Band 2 (Green), which means good maintaining of the spectral information, especially for vegetation. For the NIR channel, the error values were 7.6% and 9.7%, respectively. The obtained values may indicate degradation of radiometric correction efficiency with increasing UAV flight altitude. It is also related to the changes of the atmospheric model and changing conditions of radiation scattering with altitude. After the radiometric correction, the correlation values (Table 4) are close to 1. This means that the global statistics of the image have not changed, while the absolute radiometric quality of the image has improved.

Discussion
A method of radiometric correction of NIR images acquired from low altitudes was proposed as a result of the research. The proposed method is based on a fusion of the dark-object subtraction technique correction and the modified Walthall's BRDF model. The proposed correction method enables the reduction of radiometric differences in the images included in a block due to haze (fog) and the angle of incident solar radiation. The results of the quality assessment of the proposed method are comparable with other research. The determined RMSE do not exceed 10% in all spectral channels. It means that the selected coefficients of the BRDF correction model are reasonable. In order to prove the efficiency of the developed radiometric correction method, additional comparison analyses were conducted for NIR images acquired from low altitudes for a sample of 898 images obtained within a measurement campaign in October 2017. The images were acquired from 200 m altitude by means of the UX-5 platform in an area in southeastern Poland. The images were used to produce an orthophotomosaic. A correction according to the proposed method was done for each channel.
The obtained values of correlation for channels Band 1 and Band 3 were 0.98, and for the channel Band 2 0.99. The results confirm that the global statistics of the image have not changed. The values of RMSE (%) were determined for each channel for additional control. They were 6.3% (Band 1), 4.2% (Band 2), and 8.2% (Band 3). The results of the verification research show that the developed method of radiometric correction of images acquired in the NIR range is also applicable for other imagery data from other regions, other seasons of the year, and other altitudes. The proposed method of radiometric correction of images acquired from low altitudes should always take into account the solar zenith and azimuth angle. Images without any initial corrections, taking into account e.g., the impact of vignetting errors, were used in the conducted research. It is vital information because, often in studies of UAV data, the images are initially corrected by applied software that removes their primitive spectral information. The study of the approach using two different sets of data and a third verification set shows that radiometric correction of NIR images is necessary. As it was proved in the calculated values of correlation, the global statistics of the images have not changed and its absolute radiometric quality has improved. The improvement, even at the level of a few percent, is extremely important because it equalizes the spectral properties of the objects read from the images and it improves the interpretation of the images. Similar results were achieved in the research of Reference [36]. The authors proposed the radiometric correction of hyperspectral images for the assessment of yield quantity and quality of a grass sward. The obtained RMSE value for a corrected nadir mosaic was a bit over 13% at the correlation coefficient of 0.985. In other research [37], the authors proposed the application of a modified BRDF correction model for colorized aerial lidar point clouds. In their research, the RMSE value did not exceed 4% and it was the smallest for the green channel, which is similar to our research. The proposed method of initial radiometric correction of UAV images may be a new introduction to the precise radiometric correction of images with the application of dedicated reflection standards. The presented approach is based on modified models of correction of aerial and satellite images that were adapted for low-altitude images. The research results prove the universality of the presented approach for various altitudes and for various lighting conditions (various seasons). The presented research shows that the application of this type of correction is justified and it indicates the necessity of continuing its improvement for developing correction coefficients, taking into account the impact of the lowest parts of the atmosphere. The conducted experiments also indicate a new research direction for developing new radiometric correction models that would be dedicated exclusively to UAV images.

Conclusions
The paper presents the results of research concerning the method of radiometric correction and its assessment for NIR images acquired from low altitudes. The authors propose a correction method taking into account the impact of the atmosphere and the direction and angle of incident solar radiation during image acquisition. The image processing is realized in two phases-the first is based on DOS technique correction. The applied method of initial correction based on DOS has, in some cases, limitations due to its inaccuracy caused by erroneous determination of zero reflection of dark objects, and due to its inaccuracy in the assessment of DN values in the NIR range, where the impact of atmospheric absorption of radiation increases. An empirical BRDF model was applied in the second phase of correction. Moreover, the proposed method has a lower computational cost when compared to the traditional method of BRDF determination for the radiometric quality of images. The images processed according to the method may be used in remote-sensing applications and developing new haze models.
Future research will focus on developing haze models (especially for radiation haze) and developing radiometric correction for UAV images at altitudes higher than 500 m.
Author Contributions: All authors contributed to the experimental design and participated in the collection of the UAV data. All authors provided editorial advice and participated in the review process. Conceptualization, methodology, and results interpretation by D.W. and M.K. Data analysis by A.F. and J.J.