Next Article in Journal
Regional Models for High-Resolution Retrieval of Chlorophyll a and TSM Concentrations in the Gorky Reservoir by Sentinel-2 Imagery
Next Article in Special Issue
An Effectiveness Evaluation Model for Satellite Observation and Data-Downlink Scheduling Considering Weather Uncertainties
Previous Article in Journal
Current Activity of the Long Point Fault in Houston, Texas Constrained by Continuous GPS Measurements (2013–2018)
Previous Article in Special Issue
Perceptual Quality Assessment of Pan-Sharpened Images
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Lower Atmosphere on the Radiometric Quality of Unmanned Aerial Vehicle Imagery

Institute of Geodesy, Faculty of Civil Engineering and Geodesy, Military University of Technology, 01-476 Warsaw, Poland
*
Author to whom correspondence should be addressed.
Remote Sens. 2019, 11(10), 1214; https://doi.org/10.3390/rs11101214
Submission received: 30 April 2019 / Accepted: 20 May 2019 / Published: 22 May 2019

Abstract

:
Unmanned aerial vehicle (UAV) imagery has been widely used in remote sensing and photogrammetry for some time. Increasingly often, apart from recording images in the red-green-blue (RGB) range, multispectral images are also recorded. It is important to accurately assess the radiometric quality of UAV imagery to eliminate interference that might reduce the interpretation potential of the images and distort the results of remote sensing analyses. Such assessment should consider the influence of the atmosphere and the seasonal and weather conditions at the time of acquiring the imagery. The assessment of the radiometric quality of images acquired in different weather conditions is crucial in terms of improving the interpretation potential of the imagery and improving the accuracy of determining the indicators used in remote sensing and in environmental monitoring. Until now, the assessment of radiometric quality of UAV imagery did not consider the influence of meteorological conditions at different times of year. This paper presents an assessment of the influence of weather conditions on the quality of UAV imagery acquired in the visible range. This study presents the methodology for assessing image quality, considering the weather conditions characteristic of autumn in Central and Eastern Europe. The proposed solution facilitates the assessment of the radiometric quality of images acquired in the visible range. Using the objective indicator of quality assessment developed in this study, images were classified into appropriate categories, allowing, at a later stage, to improve the results of vegetation indices. The obtained results confirm that the proposed quality assessment methodology enables the objective assessment of the quality of imagery acquired in different meteorological conditions.

Graphical Abstract

1. Introduction

With the intensive development of the unmanned aerial vehicle (UAV) technology, the use of low-altitude imagery is becoming increasingly popular. UAV imagery has been applied in environmental remote sensing, climate change research, and atmospheric composition studies. Intensive research in recent years has involved applying UAV remote sensing techniques in environmental and topographic monitoring [1], forest studies [2], and ecohydrology [3]. Using UAVs with a compact multispectral digital camera has reduced the cost of acquiring image data [4]. Using UAV photogrammetry, it is possible to quickly obtain extremely high spatial resolution images for topography mapping [1], three-dimensional (3D) modelling [5], or point cloud classification [6,7]. However, studies have not yet considered the broader aspect of the influence of atmospheric and meteorological conditions on the quality of UAV imagery. Until now, spectral measurements in environmental remote sensing studies were acquired on the basis of data from sensors on-board satellites and manned aircraft. With the advent of UAV imagery technologies, new types of sensors, including remote sensors, have emerged. However, studies have not always considered the radiometric quality of UAV imagery appropriate for applying this type of data in advanced spectral analyses. Differentiating the imaging sensors installed on UAVs requires the development of standardized procedures for assessing the quality of images acquired with these sensors. This is vital, especially for quantitative and spectral remote sensing approaches. An accurate assessment of radiometric quality is also important when considering the complexity of geometric and radiometric corrections necessary for accurate spectroscopy-focused environmental remote sensing.
Low-altitude remote sensing enables acquiring spectral information by measuring the radiation emitted (e.g., measuring the spectral reflectance characteristics of agricultural crops). However, the accuracy of low-altitude spectral information is influenced by image quality, which in turn depends on the lighting and atmospheric conditions, as well as on the type of measurement protocol and data-processing procedure applied [8]. Thus, developing objective radiometric quality indicators is crucial for eliminating incorrect data, which can significantly impact the results of the spectral analyses conducted and their later application for understanding Earth’s biological processes [9]. Due to the ever-increasing popularity of small spectral cameras mounted on UAVs, e.g., Tetracam MCA (Multiple Camera Array), Parrot Sequoia, Micasense, etc., there is a need for developing reliable indicators for classifying the images to given quality groups. The quality of UAV imagery data depends above all on the sensor characteristics, lighting geometry, and meteorological conditions. Therefore, digital number (DN) values are not truly representative of the surface reflectance, as they depend on the lighting conditions and the characteristics of the imaging sensor [10]. Thus, multispectral UAV imagery data should be assessed in terms of quality, considering the influence of lighting and atmospheric conditions, camera noise, and topographic conditions. Only based on an accurate assessment of the radiometric quality of UAV images and their subsequent correction can raw data be processed into physical units of reflectance that are used to develop spatial-temporal spectral representations of the objects studied, e.g., chlorophyll content in plants.
By considering the influence of atmospheric conditions on image quality, it is possible to increase the reliability of absolute radiometric calibration, during which a raw DN is transformed into an absolute reflectance value for every spectral band. In the case of relative radiometric calibration for the entire image, normalized transformation values are applied with the assumption that the sensor’s focal plane is irradiated with a uniform radiance field [10]. In other cases, DN values are transformed into spectral reflectance characteristics on the basis of the real values of spectral reflectance characteristics measured directly in the field with a spectroradiometer. The spectral range of such measures should be similar to the ranges of the UAV-mounted sensor. The images should include photographs of representative surfaces (reference panels) with known spectral characteristics. Due to the small Instantaneous Field of View (IFOV), the atmosphere is assumed to be homogenous throughout the entire image. When calibrating such images, the atmospheric profile should be measured on the operating altitude; for UAV remote sensing, it is usually between 50 and 300 m. The aim of the radiometric correction of UAV imagery is to eliminate any lower atmosphere interferences, as sunrays reflected from the photographed objects are often non-Lambertian, which leads to intensifying the hot-spot effect. Radiometric correction models can be divided into two types: physical and empirical. So far, several different methods of implementing physical models into radiometric correction of UAV imagery have been tried, e.g., MODTRAN (MODerate resolution atmospheric TRANsmission) [11,12]; however, these have not always worked well. When assessing UAV imagery, the empirical models are more practical, for which model parameters can be calculated. Common empirical methods include histogram matching methods, dark object methods, and empirical line methods [13,14]. Several remote sensing applications require the elimination of the influence of object reflectance anisotropy. For anisotropic surfaces, a change in any of the parameters leads to a change in surface luminance [15,16]. The characteristics of light reflectance for the given surface is described by the Bidirectional Reflectance Distribution Function (BRDF), first defined by Nicodemus [17].
Increasingly often, UAV imagery is acquired in different weather and lighting conditions, which is why new methods for quality assessment need to be developed for the purposes of radiometric correction. This has been attempted previously [18,19,20,21,22]. Assessing the quality of the images and classifying them into the right groups is an indispensable stage before applying radiometric correction.
Researchers have increasingly often mentioned radiometric correction methods based on the empirical approach [23,24,25]. The idea of these methods is to find the connection between known reflectance targets in relation to the DN value measured on the image. The best solution with this approach is to convert all DN values for every spectral band of the image to surface reflectance values. New methods of radiometric correction have been proposed, allowing the generation of an orthomosaics reflectance image without knowing the empirical coefficients. Such solutions are possible due to introducing a priori calibration in multispectral cameras, such as Parrot Sequoia® (Parrot Drone SAS, Paris, France) and MicaSense RedEdge® (MicaSense Inc., Seattle, WA, USA) [26]. Also, with the use of commercial software, e.g., Pix4D (Pix4D, Lausanne, Switzerland), it is possible to conduct color correction and histogram matching, or implement the image-information-based radiometric block adjustment algorithm. These solutions may not return accurate and repeatable results in the remote sensing context, especially when processing images acquired in different seasons and different lighting conditions. That is why it is difficult to select the right normalizing parameter. Additionally, considering the increased variability in the measured spectrum, the reflectance coefficient observed at smaller pixel sizes creates a significant problem for effective radiometric normalization. As a result, low accuracy normalization, when performing remote sensing UAV measurements, effectively limits the possibility of comparing images from the same area acquired at different times and in different seasonal conditions [25]. Therefore, it is crucial to develop objective indicators of image quality that can be normalized depending on the time of year before implementing further stages of radiometric correction.

1.1. Related Works—Classical Quality Assessment Indexes

Image quality can be defined by several different parameters, such as radiometric resolution and its accuracy, represented by the noise level, and geometric resolution and sharpness, described by the Modulation Transfer Function (MTF) [27]. The radiometric quality of the digital image can be defined as a detailed mapping of local irradiation changes recorded by the imaging system, while maintaining a continuum of brightness adequate for the mapped scene. The internal quality digital image radiometric is formed by the local image contrast, tonal range, random noise, and radiometric resolution [28,29]. Image radiometric quality is influenced by the sensor characteristics defined by sharpness, contrast, and resolution [30,31], but also by the atmospheric composition and meteorological conditions at a given time.
In digital image processing, image quality is a multi-aspect problem depending on the application of the processed image. Images represented as matrices are subject to objective and subjective quality assessments. The majority of traditional image assessment methods are comparative methods based on the assumption that the assessed image is compared to the original image. Studies are currently underway to find reliable image quality assessment indexes. At the moment, objective quantitative indexes are usually used for image quality assessment based on comparing the differences between the original image and the image reconstructed after compression [32,33,34]. One of the most popular and simplest objective indexes is the Mean Square Error (MSE) index, i.e., the value of the mean square error for each pixel of the image. To determine MSE, the value of the difference between the reference image and the distorted image is calculated [35]. This index was further developed as the Root Mean Square Error (RMSE)—the square root of MSE that is used to calculate signal-to-noise ratio (SNR) and peak signal-to-noise ratio (PSNR) [34]. This index is described by [36]:
PSNR = 10 log 255 2 N 2 i e i 2 [ dB ]
where 255 is the dynamic range of the signal, N is the number of pixels, and ei is the difference between the reconstructed and the original image calculated as the ith pixel.
Although PSNR is simple and widely used, it is not resistant to lighting conditions, so the results may differ from the subjective assessment and may not correlate with human perception. Another quality assessment index is SNR, which calculates the size of the error with respect to the reference image. This index, in turn, served as the basis for developing the signal-to-quantization-noise ratio (SQNR) measure, which determines the impact of quantizing signal quality. This index can be described by [34]:
SQNR = 10 log 10 signal   power quantization   noise
where quantization noise represents the difference between the quantized and reference images.
Another important image quality index is the noise quality measure (NQM), a subjective quality measure for images that cause distortions and noise in a video sequence. These two degradation factors are considered as independent and are assessed with two quality measures: a distortion measure (DM) for frequency distortions, and a measure of the impact of noise on the video sequence quality (NQM). The NQM index considers variations in contrast sensitivity as a function of resolution, variations in the local luminance mean, interactions between contrast and spatial frequencies of the image, and contrast masking effects. The DM index is calculated in three stages. The first step consists of finding frequency distortions in the degraded image. Next, the deviation value of this frequency distortion is calculated. At the final stage, the frequency distortion deviation value is weighted so that it fits the model of the frequency response of the human visual system [37,38]. Image quality is also assessed with the use of the universal image quality index (UIQI), developed for computing the loss of structural information on images with low radiometric quality. This index presents image quality as a combination of three factors: correlation distortion, luminance distortion, and contrast distortion. This coefficient is referred to as universal because it is independent of the type of distortions, image scale, or type of sensor used. Studies have demonstrated its superiority to the MSE index [39]. Image signal quality is often assessed using the structural similarity index (SSIM), which considers three types of distortions: image luminance, contrast, and structure. The coefficient can be described as:
SSIM = ( 2 μ x μ y + C 1 ) ( 2 σ x y + C 2 ) ( μ x 2 + μ y 2 + C 1 ) ( σ x 2 + σ y 2 + C 2 )
where µx is the mean value on a good-quality image (acquired with no overcast and on flat terrain); µy is the mean value on the compared image; σx and σy stand for the variance of the reference and the compared image, respectively; σxy is the covariance between the reference and the compared images; and C1 and C2 are constant values to prevent the equation from being meaningless. In the experiments, C1 and C2 have been assigned values of 0.01 and 0.03, respectively, and the dynamic range of the image was from 0 to 255.
The SSIM coefficient is characterized by contrast, brightness, and structural similarity of an image, taking values from 0 to 1. The closer the value to 1, the higher the image quality [40].

1.2. Research Purpose

This article presents a new objective index for assessing the quality of UAV imagery, considering the impact of air humidity and the sun zenith angle. The radiometric assessment of UAV images should be the first step in the process of acquiring images with remote sensors, especially in cases where we want to obtain data such as, for example, reflection coefficients necessary for radiometric correction of images. Therefore, it is important to develop an objective indicator of radiometric quality for UAV imagery, which considers the influence of the low atmosphere and direction of sunlight during acquisition of image data. The effectiveness of the developed indicator was verified using two independent test sets.
The aim of this study was to determine the impact of atmospheric fronts, including air humidity and the direction of illumination of the photographed scene, on the quality of low-altitude images. Based on specific relationships between the atmospheric conditions and lighting direction, and the differences in the reflectance values determined on the basis of dedicated radiometric calibration panels, a new image quality indicator was developed. This indicator will be useful in radiometric correction methods, as well as in remote sensing analyses.
This paper discusses the impact of meteorological conditions on the radiometric quality of UAV images, presenting the results of the proposed new quality index and an accuracy analysis based on reference image data.
The paper is structured as follows. In Section 2, the test data are introduced. In Section 3, the research method is explained. Section 4 presents experimental results. In Section 5, the results are discussed. Finally, Section 6 provides a brief summary of this work.

2. Materials

2.1. Test Area

The correction method was designed based on a series of data obtained in the lakeside village of Mieruniszki (54°10′7.15″N, 22°34′12.18″E; 193 m a.s.l. (above sea level) on September 13, 2018. This area is located on the northern end of Poland and is characterized by high latitude, so the angle of sunlight is lower, the day is shorter, and the photos acquired in the morning and evening may be of worse radiometric quality than photos obtained at the same time under the same conditions for higher latitudes. The proper selection of test data to assess the correctness of the proposed method was particularly important. The first test set was acquired in the Mieruniszki area on November 22, 2018. The second set of data was recorded in the mountain village of Kościelisko (49°17′32.63″N, 19°52′26.28″E; 941 m a.s.l.) in Southern Poland (Figure 1) on September 26, 2018. Kościelisko is located in the south of Poland and its latitude is about 5° lower than that of Mieruniszki. Both regions have their own microclimate. In the Polish mountains during the climatic autumn, strong winds are often blowing. The biggest impact on the climate is polar air from the west, which, during the winter, causes warming and an increase in cloud cover and precipitation on the northern slopes of the mountains. In summer, it causes cooling and an increase in cloud cover. In turn, Mieruniszki is located in the Polish cold pole, where the air temperature is generally lower than in the south of the country. Winds from western directions prevail, with winds from the south-east. The most windless days fall between July and September. The local microclimate and current atmospheric fronts create the synoptic situation of the region. For the test areas, the synoptic situation was investigated and described in Section 2.2. This selection of test data allowed us to verify the proposed solution for different climatic conditions and different heights of the sun.
Test areas were covered with low grass vegetation with single trees and shrubs. They were both characterized by a low degree of urbanization, with single-family houses, and road and technical infrastructure. Test flights were conducted between September and November 2018. In the climate zone of Central and Eastern Europe, this period constitutes autumn.

2.2. Data Acquisition

RGB-range image data were obtained with the DJI Phantom 4 Pro (DJI, Shenzhen, China) platform equipped with four electric motors. The platform is fitted with a high-resolution digital camera providing images in the visible range. The systems used for navigation are GPS/GLONASS (Global Positioning System/ Global Navigation Satellite System) and an optical positioning system, as well as two IMU (Inertial Measurement Unit) systems. The platform is steered with an remote control (RC) controller operating on the 2.4 GHz frequency. A 1-inch 20-megapixel complementary metal oxide semiconductor (CMOS) sensor with mechanical shutter was mounted on a three-axes stabilized gimbal. The focal length of the lens was 24 mm (full-frame). The data were stored in 24-bit JPG format. Camera sensitivity was set at ISO200 for all images with aperture 4.0 and shutter times ranging between 1/25 and 1/800 s. All images were georeferenced with the on-board GPS. Image data were acquired at altitudes between 50 and 300 m. For each altitude, 70 images of a similar texture (flat areas, medium level of urbanization, fields, meadows) were acquired.
For the purposes of calibration and determining quality coefficients, six radiometric calibration panels with greyscales ranging from black to white were developed, as per the method suggested by Deng et al. [22]. The reflectance value of the calibration panels (Figure 2) was measured with the use of the ASD (Analytical Spectral Devices) FieldSpec spectrometer (Malvern Panalytical Ltd., Malvern, UK) and was, respectively, from 5% (P1 panel) to 95% (P6 panel) (Table 1) for all the spectral bands corresponding to the DJI (Dow Jones Indices) camera bands.
Meteorological data used for developing the parameters of radiometric quality assessment were acquired from the AGAT20 ground-level measurement station and the SR10 atmosphere probing system. The SR10 system, designed for areological measurements, enables the measuring of physical parameters of the atmosphere—air temperature, dew point, humidity, pressure, and wind speed and direction—from ground level to the height of approximately 30 km, with a frequency of 1 Hz. The system is comprised of an SR10 receiver, M10 radiosondes (Figure 3), a GNSS (Global Navigation Satellite System) ground control set for the calibration of radiosondes and based on initialization (12 bands). Additionally, the system includes an omnidirectional antenna (400 MHz), a GPS antenna (400 MHz), and a tripod. The SR10 system operating range is minimum 350 km. The receiver operates in the 400–406 MHz frequency range.
AGAT20 (Figure 4) is an autonomous mobile surface weather station for conducting measurements in the field. The station measures the following meteorological elements: atmospheric pressure at the station level, temperature and humidity at 2 m above ground, and wind at 4 m above ground.
The measurement system was comprised of a wind speed and direction measuring head, air temperature and relative humidity measuring head, pressure measuring head, a measuring stand with mast, a measurement station laptop with MetNet software, and power and signal distribution box. The measuring station should be situated away from any natural or artificial terrain obstacles impeding the free flow of air.

2.3. Meteorological Conditions

2.3.1. Synoptic Situation in the Mieruniszki Test Area

On September 13, 2018, the area of Mieruniszki was under a low-activity cold front associated with a barometric low from Northern Scandinavia (Figure 5). A maritime polar air mass flowed over the region from the west, at a velocity of 50–40 km/h, and with stable equilibrium that changed to conditionally unstable as the day passed.
At the time of the flights in Mieruniszki, between the hours of 04:00 UTC and 12:00 UTC, the sky coverage was at 8-5/8 Ci, Ac, As, Cu, Sc, with the lowest cloud bases at the height of 400–600 m, rising after 07:00 UTC to 600–1000 m. The prevailing wind direction was 250–270° and the wind speed was 2–4 kt. Near the end of the considered time interval, there were light rain showers. Prevailing visibility was above 10 km, with fog only in the early hours of the morning, limiting the visibility to 6 km. During the flights, the maximum temperature was 15–17°C. The 0°C isotherm was located at 2700–3000 m (Figure 6).
On November 22, 2018, the area of Mieruniszki was in a high-pressure field associated with the barometric high from Kaliningrad (Figure 7). A continental polar air mass flowed over the region from the southwest, at a velocity of 20–30 km/h and with stable equilibrium.
During the flights in Mieruniszki, between the hours of 06:00 and 12:00 UTC, the sky coverage was at 6-8/8 with the lowest cloud bases at the height of 300–450 or 600 m. The wind was moderate and changeable (prevalently north-easterly). Prevailing visibility was above 10 km. During the flights, the maximum temperature was –3 °C. The 0 °C isotherm was located at 0 m. Within the range of the clouds, there was light to moderate icing. Numerous inversions, which had developed after the sunset on the previous evening and during the night, lasted throughout the day, which prevented any changes in the size and height of cloud bases, as well as any increases in temperature (Figure 8).

2.3.2. Synoptic Situation in the Kościelisko Test Area

On September 26, 2018, the area of Kościelisko was in a stationary high-pressure field coming in from Southern Europe (Figure 9). A maritime polar air mass flowed over the area from the northwest at the velocity of 50–45 km/h and with stable equilibrium.
During the flights in Kościelisko between the hours of 07:00 and 17:00 UTC, the sky coverage was 2-4/8 Ci, Cs, Ac, with lowest cloud bases at the height of 2400–3000 m, with periodic changes to 2-4/8 Cu, Sc, cloud bases at 1800–2400 m. The wind was moderate and changeable, with visibility over 10 km (Figure 10). The maximum temperature was 11–13°C. The 0°C isotherm was located at 3000–3400 m.

3. Proposed Methodology for Radiometric Quality Assessment

At the first stage of this study, based on the previously presented concept [44,45], every image acquired from the flights was divided into three bands: red, green, and blue. For each band in each image, the average brightness value of the pixel and the standard deviation of the brightness value of the pixel were computed. Initially, the images were divided into 100 equal fragments. For each fragment, the standard deviation of the brightness value of the pixel for each of the bands was determined. Based on the visual analysis of the images and the distribution of standard deviation values, a pattern was observed that had not been visible in studies that used traditional image quality measures. We found that images acquired in good lighting and weather conditions have high standard deviations values for the red band. For the images acquired in poor lighting in poor weather conditions, in turn, the standard deviation values for the blue band were noticeably higher than the standard deviations for the red and green bands for the same image [44]. This occurred because heavy overcast and poor lighting produce a greater dispersion of light on particles (precipitation) larger than the wavelength (Mie dispersion) [44,45]. We developed an original formula for determining the value of the radiometric quality assessment coefficient for images acquired in high-pressure areas and in autumn conditions for Central and Eastern Europe. The weight values in the WKW formula component were based on the properties of luminance, which can be defined as follows [44]:
L = 0.299 R + 0.587 G + 0.114 B
where L is the image luminance, R is the red band, G is the green band, and B is the blue band.
The weights for the particular bands can be expressed by means of coefficients as follows [44]:
W K W = ( μ R σ R ) w R + ( μ G σ G ) w G + ( μ B σ B ) w B
where WKW is the image quality assessment index value; µi is the average DN pixel value in a given band (DN in an 8-bit scale, 0–255); σ i is the pixel brightness standard deviation value in a given band; w R is 0.299, which is the weight value of the red band determined empirically based on the luminance value; w G is 0.587, which is weight value of the green band determined empirically based on the luminance value; and w B is 0.114, which is the weight value of the blue band determined empirically based on the luminance value.
The index WKW (Equation (5)) is not able to consider the influence of humidity and the solar zenith angle on the image quality. Therefore, even though atmospheric parameters and solar zenith angle do effect the pixel values, WKW was not able to measure this influence, because it is only based on the statistical characteristics of the acquired data by the simple weighted sum of three ratio. Precisely for these reasons, the improved variant is proposed. In this way, these factors are explicitly considered.
We studied the influence of atmospheric conditions on UAV imagery radiometric quality, which then served as the basis for the proposed solution for the assessment of such data. During the day (from dawn till dusk, when photogrammetric flights could be conducted), a number of measurement sessions were completed at different altitudes, while monitoring the synoptic situation. The radiosonde was also registering information on the atmospheric conditions. The analysis of the synoptic situation showed that the biggest changes occurred with humidity. Therefore, we then focused mostly on assessing the impact of humidity on radiometric quality. Due to the data gathered, it was possible to determine the functions describing the change in humidity and the change in the reflectance coefficient for calibration tests during the day. The correlation coefficient for the function describing the changes in humidity and the reflectance coefficient during the day computed for each flight altitude of the UAV platform were considered as measures of the relationship between the functions. Changes in the reflectance coefficient were determined on the basis of reference panels with known spectral characteristics, as the difference between the theoretical and practical values. The studies used panels from P1 to P6, where the theoretical reflection coefficients ranged from 2 to 95%. The correlation coefficient was calculated separately for each UAV flight altitude, as presented in Figure 11.
As demonstrated in the figures, regardless of altitude, a strong correlation exists between humidity and differences in reflectance. The correlation coefficient value is higher than 0.90, with most of the values even higher than 0.95. This high correlation between air humidity and the difference in the reflection coefficient points to the significant impact of humidity on UAV imagery radiometric quality. This problem is important especially in the case of photographs with low radiometric quality, where the correct light reflection properties are disturbed (Figure 12).
It was, therefore, crucial to develop a model for radiometric correction, which consider the impact of air humidity on pixel brightness and enable a correct imaging of objects with different levels of light reflectance. Due to the large numbers of photographs in the photogrammetric blocks, we propose selecting the images that would need correction due to the influence of air humidity. Thus, we propose using a special QA indicator for assessing the quality of low-altitude photographs.
The position of the sun and the direction of the radiation on the test surface is determined by basic angle parameters (Figure 13).
The θz angle is measured between the angle of the rays falling directly on the given surface and the normal to this plane, and is expressed as follows [46,47,48]:
cos ( θ z ) = sin ( δ )   [ sin   ( )   cos ( β )   cos ( )   sin ( β )   cos ( γ ) ] + + cos ( δ )   [ cos ( )   cos ( β )   cos ( ω )   +   sin ( )   sin ( β )   cos ( γ ) cos   ( ω ) + + sin ( β )   sin ( γ )   sin ( ω ) ]
where β is the angle of inclination of the considered surface in relation to the Earth’s surface, φ is the latitude of the photographed location, and ω is the hour angle, i.e., the deviation of the sunrise or sunset from the local meridian as a result of Earth’s rotation about its axis. Values of the hour angles for every hour of photography (hours of the solar day; h) are expressed with:
ω = 15 ( h 12 ) ( ° )
The δ angle (°) is the solar declination, i.e., the angle between the sun and the plane of the celestial equator. The value of declination for the given day of the year (n) is computed with Cooper’s formula:
δ = 23.45 sin ( 360 284 + n 365 )
In the case of a horizontal surface, when β = 0, the equation for the solar zenith angle θz (during the day, when the sun is over the horizon, the angle must be in the range between 0 and 90°) will be as follows:
cos ( θ z ) = sin ( δ ) sin ( φ ) + cos ( δ ) cos ( φ ) cos ( ω )
Considering the influence of the atmospheric conditions (humidity) and solar zenith angle θz, we developed an original formula, shown in Equation (9), enabling the calculation of the coefficient of the radiometric quality assessment for images obtained in a high-pressure area and in autumn conditions for Central and Eastern Europe. The solar zenith angle was calculated based on the geolocation and the time of UAV flight of study area. The solar zenith angle was in the range between 2 and 40°:
Q A = W K W h sin ( θ z )
where WKW is the quality assessment index for RGB images from low altitudes, h is the humidity, and sinθz is the solar zenith angle. The minimum theoretical value of the Q A indicator occurs when one of the factors in the nominator equals 0. This situation is only possible in the ideal case when the photograph is taken with perfect exposure parameters (WKW = min) or when the theoretical humidity is at 0%, which is not possible in real life. Then, Q A =   w k w   ×   0 sin ( θ Z ) = 0. In reality, however, both WKW and humidity can approach 0, and then Q A =   w k w   ×   1 sin ( θ Z )   0, which corresponds to images with good radiometric quality.
The denominator of the proposed equation cannot be equal to 0, and, in practice, it is impossible to take any photograph with the sun at zero height. When lim θ Z 0 ( sin θ Z )   = 0 , then Q A =   w k w   ×   h   sin ( θ Z )   This describes a situation corresponding to obtaining photos with low radiometric quality. The higher the WKW and humidity values, the closer the   Q A to infinity.
In the physical interpretation, we consider a deviation from the ideal model. The higher the value of this deviation, the lower the image quality.
Similarly to the experiments conducted in Kedzierski and Wierzbicki [44] on the basis of the visual analysis of images and their histograms for individual RGB bands for the obtained WKW coefficient values, we determined the intervals and classified images as having good radiometric quality, medium radiometric quality, or bad radiometric quality; the latter needed to be subjected to atmospheric correction with the use of a dedicated model, for example, as implemented in the ATmospheric CORrection (ATCOR) [49].
Selection parameters for the images were supported by the analysis of the images themselves, and of the QA indicator value for the reference block of data acquired at a low altitude in good lighting conditions.
Based on empirical studies and statistical analyses, we determined value intervals for the QA coefficient (Table 2).
The limits of the QA coefficient intervals were assigned based on the value of the average and the multiple of the standard deviation value for the coefficient studied for a sample of reference data. The upper limit of good-quality images was set as the average value increased by twice the maximum standard deviation for the tested blocks, which was 6.00. For the medium quality category, the upper limit was the average value increased by three-times the maximum standard deviation, which was 7.65. Images with bad radiometric quality had an indicator higher than 7.65. The higher the coefficient value, the lower the image quality. As the indicator approaches the maximum value, the image approaches the so-called dark image, which is blurry and out of focus.
When examining the properties of the proposed QA indicator, we conducted case studies of photographs captured in different conditions. We adopted the assumption that all photographs were taken with exposure parameters chosen so that the intensity of light falling on detector arrays was always the same. This guaranteed the same values of the WKW factor for both photo sets; for example, this value for good-quality photographs can equal 2 [44]. Our proposed solution was verified by assessing the radiometric quality of photographs obtained during separate flights in both good and bad weather conditions. A series of flights were conducted in early autumn (September 13 and 26, 2018) with high and uniform insolation, and in late autumn with a high degree of cloud cover (November 22, 2018). The data were registered at three different altitudes at three times of day, so with different solar angles and air humidity. The QA indicator was calculated for all the images.
QA can be considered a way that integrates WKW index in order to explicitly take into consideration humidity and solar zenith angle on the quality of the image. The two additional terms of QA compensate the circumstance that WKW is not able to evaluate the influence (that exists) of humidity and solar zenith angle on the acquired data. Therefore, even though the correction proposed is empirical and very simple, the experiments carried out have proved that it is effective.

4. Results

Consider the case of photographs taken in September on 50° N latitude, just after dawn at 06:00 (sun altitude approx. 5°), 07:00 (sun altitude approx. 14°), and at 12:00 (sun altitude approx. 38°) with air humidity at 80% and 40%.
As Equation (10) and the case study presented in Table 3 demonstrate, higher-quality images can be sometimes be more readily obtained in the morning, with no clouds and low humidity, than at noon with an overcast sky, and thus high air humidity. According to the proposed classification, the photographs capture right after dawn have low radiometric quality. This is consistent with common knowledge and experience; when the angle of the sunrays is too low, there are many shadows and the image is noisy. In turn, photographs taken with the sun at 14° and low humidity have a similar quality to those taken at noon with high humidity. The quality of photographs captured in the morning strongly depend on air humidity and the quality changes dynamically along with the rapid change in the sun’s vertical position in the sky after sunrise. A parallel pattern applies to images obtained in the evening, close to sunset. The highest quality imagery is captured with photographs at noon with low humidity when the QA coefficient value approaches 0. The proposed indicator closely corresponds to the actual atmospheric conditions during photogrammetric flights and can be used to predict the quality of images obtained on a given day, at a given time, and with the known forecast humidity levels.
Figure 14 demonstrates the distribution of the QA indicator values for reference images acquired at three different altitudes in the early morning, late morning, and afternoon.
The average coefficient values for reference images ranged from 1.4 to 5.4, with the standard deviation between 0.2 and 1.5. Maximum values do not exceed 9. The minimum value of the coefficient was 0 and is a theoretical value that the coefficient would reach for photographs taken in perfect conditions, i.e., with zero humidity and a 90° sunlight angle. Such conditions are unfeasible, so the indicator will always have values higher than 0.
In the Table 4 we presented results of the statistical parameters of the indicator for each photo set.
In the case of data acquired in good conditions, QA values never exceeded 10, independent of the UAV flight altitude. With data obtained before noon, and thus with highest insolation and lowest humidity, all photos were categorized as good-quality. In the sets of photos from the morning and afternoon, there appear some photos had a coefficient exceeding 7.65, which suggests low image quality. This was caused by local radiation fogging. The relatively low SD values demonstrate a narrow distribution of values around the average, which means that the majority of images are characterized by good radiometric quality (especially with photos taken in the morning and afternoon). In turn, for data acquired with high overcast in the morning and afternoon, QA values significantly exceeded 7.65 and were categorized as bad radiometric quality. Only the photos taken before noon, when the lighting was more even, were classified as medium quality, and few were classified as good quality. The relatively high SD values in the morning and afternoon photo sets suggested the occurrence of both good- and bad-quality images; therefore, analyzing the QA coefficient for each photo enables selection of photos in need of radiometric correction.
Figure 15 presents sample images from both data sets. Visual analysis confirmed that the coefficient works correctly. Photos taken in good conditions are sharp and have a low degree of noise, which is reflected by low QA values. Images acquired in poor conditions, in turn, are blurry, hazy, and noisy, as supported by low coefficient values.
Reflectance panels were also photographed in the test images to verify the accuracy of the proposed QA indicator. Figure 16 and Figure 17 depict diagrams of air humidity and residuals from the theoretical value of the reflection coefficient for test with minimum (2%) and maximum (95%) reflectance for two chosen flight altitudes, first for Kościelisko, and then for the Mieruniszki test area.
Figure 16 and Figure 17 provide the results for sets of data obtained in good weather conditions. The cloudless sky guaranteed a uniform illumination of the entire test area. As shown above, the photographs were clear and free of noise. The absolute values of reflectance residuals for the dark test (2% reflectance) did not exceed 0.2, and for the light test (95% reflectance), they did not exceed 0.4 in every RGB band. Such discrepancies between the theoretical and practical values indicate high radiometric quality and correct registering of light reflecting from objects, which is consistent with the QA assessment.
Figure 18 depicts the correspondent analyses of data obtained in poor weather conditions. The most noticeable difference is the increase in the residual for data obtained in poor weather conditions. The absolute value of these residuals increases by approximately 0.2 regardless of the UAV flight altitude, reaching values of 0.4 and 0.6 for the dark and light tests, respectively. Due to clouds, the photographed area was unevenly illuminated, which caused differences in the registered light reflection from the calibration test. Due to the high degree of clouding for the entire day, air humidity underwent only minor changes. There was a noticeable correlation between humidity and the reflection residual size. As the humidity dropped, the difference between the theoretical and practical reflectance coefficient value increased. The correlation might be slightly weaker for images registered in cloudy conditions, but strengthened with photographs captured with little or no cloudy. The conducted analyses confirm that the QA indicator facilitates the selection of images in need of special radiometric correction.

5. Discussion

Fast and accurate quality assessment of UAV imagery, especially for the purposes of quantitative and qualitative analyses, is an area of active research in remote sensing technologies. Studies conducted to date have not considered the impact of lower-altitude atmospheric conditions on the radiometric quality of UAV imagery. Some studies only mentioned quality assessment based on laboratory experiments, which are usually difficult to replicate in real life when poorer atmospheric conditions have to be considered.
In this study, we proposed a new image radiometric quality assessment indicator that accounts for the influence of air humidity and solar zenith angle. The conducted experiments confirmed the impact of humidity (especially radiation fogs) on the quality of low-altitude images. Atmospheric correction for low-altitude UAVs in environmental monitoring has been generally neglected to date. However, in our experiment, similar to that conducted by Yu et al. [50], we demonstrated that image quality is closely linked to Mie dispersion and Rayleigh scattering, as well as other factors, such as radiation fog [51]. We observed that for images captured during high humidity, there is a greater dispersion of light on water vapor molecules, which is dependent on the wavelength (Mie dispersion). The study results also confirm the conclusions reported by Yang et al. [52], where the relationship between image quality and lighting direction was analyzed. We proved that for areas with no cloud and low humidity, as well as uniform lighting of the entire imaging area, the images are characterized by good radiometric quality, and atmospheric damping does not occur. This was verified by examining the differences between the actual reflectance values and the values read from the images. Absolute reflectance residual values for the dark test (P1, 5% reflectance) did not exceed 0.2, and for the light test (P6, 95% reflectance) they did not exceed 0.4 for any RGB band. Such discrepancies between the theoretical and practical values indicate high radiometric quality and correct registering of light reflecting from objects, which is consistent with the QA assessment. With the humidity level exceeding 50%, the absolute values of reflectance residuals for the dark test (P1, 2% reflectance) exceeded 0.4, and for the light test, they were even higher than 0.6 (P6). Images captured in such atmospheric conditions and humidity were out of focus, blurry, and noisy, which was reflected in the low QA values. Therefore, it is crucial to apply atmospheric correction when processing UAV imagery data.
As demonstrated by the results of the verification studies, the methodology adopted for assessing the quality of images obtained under different atmospheric conditions can also be applied to image data acquired in different areas. The proposed methodology considers air humidity for a given altitude, as well as the solar zenith angle. As shown by the correlation coefficient (CC), there is a close link (high correlation) between humidity at a given altitude and the difference in reflection. In many cases, the correlation was higher than 0.8, with the majority of values exceeding 0.9. High correlation between humidity and the reflection coefficient difference indicates a significant impact of humidity on the radiometric quality of UAV images. The results of our experiment show that when assessing the quality of UAV imagery in the process of radiometric correction, it is vital to consider the impact of humidity for the given altitude, as this atmospheric factor distorts the correct reflectance values. The methodology proposed for assessing the impact of humidity and solar zenith angle might constitute a basis for a new model of atmospheric correction using dedicated reflectance panels. We proved that our approach is universal for the most commonly used heights in UAV remote sensing, i.e., heights from 75 to 300 m, lighting conditions, and humidity levels. The presented research demonstrates that performing this type of assessment and classifying the images into different quality groups is a legitimate and necessary approach. The conducted experiments indicate a new avenue of research to develop new models for image radiometric correction, dedicated solely to UAV imagery.
The main limitation of the proposed indicator is that it can only be used for images obtained in the visible range. Its efficiency depends on external measurements of humidity at the given altitude. The studies were conducted during the autumn period in Central and Eastern Europe, and so the thresholds set for classifying the image quality are dedicated for data obtained in similar conditions. New quality classification limits might need to be designated after similar studies are conducted during winter, and particularly in the summer, when there are significant changes not only to air humidity, but especially to the position of the sun. Thus, the indicator might approach higher values and produce different statistics. Similar problems may arise for data registered in other areas of the world. The studies were conducted in Europe, at approximately 50 °N latitude. The height of the sun also changes with the location of the imaged area, which influences the value range of the radiometric quality indicator.

6. Conclusions

Based on the conducted studies, we developed a new quality assessment indicator for images obtained in different atmospheric conditions and with uniform lighting. The indicator considers the impact of air humidity on the quality of UAV images. The proposed quality assessment method was studied for altitudes from 50 to 300 m. This is a useful operating range for UAVs registering image data for remote sensing. This new, no-reference indicator enables the objective assessment of the quality of images obtained in the visible range, accounting for the impact of humidity and solar angle at the time of obtaining the images. The indicator was developed in two stages. In the first stage, statistical parameters of the images were studied, as well as the relationship between relative luminance and image quality. The second stage considered the impact of air humidity and the direction of lighting of the photographed scene. We also verified the developed indicator on two data sets of images obtained in different atmospheric and lighting conditions. Images classified on the basis of the proposed indicator can be used in remote sensing applications, which will allow the development of new haze models dedicated to UAV image data.

Author Contributions

All authors contributed to the experimental design and participated in the collection of UAV data. All authors provided editorial advice and participated in the review process. Conceptualization, M.K. Methodology, D.W. and A.S., who also interpreted the results and wrote the paper. Data acquisition, P.W, A.F., and J.S.

Funding

This research was supported by the Ministry of National Defense Republic of Poland, Grant No. GB/1/2018/205/2018/DA-990.

Acknowledgments

This study was supported by the Military University of Technology, the Faculty of Civil Engineering and Geodesy, Department of Remote Sensing and Photogrammetry.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gonçalves, A.J.; Henriques, R. UAV photogrammetry for topographic monitoring of coastal areas. ISPRS J. Photogramm. Remote Sens. 2015, 104, 101–111. [Google Scholar] [CrossRef]
  2. Qiu, Z.; Feng, Z.-K.; Wang, M.; Li, Z.; Lu, C. Application of UAV Photogrammetric System for Monitoring Ancient Tree Communities in Beijing. Forests 2018, 9, 735. [Google Scholar] [CrossRef]
  3. Vivoni, E.R.; Rango, A.; Anderson, C.A.; Pierini, N.A.; Schreiner-McGraw, A.P.; Saripalli, S.; Laliberte, A.S. Ecohydrology with unmanned aerial vehicles. Ecosphere 2014, 5, art130. [Google Scholar] [CrossRef] [Green Version]
  4. Colomina, I.; Molina, P. Unmanned aerial systems for photogrammetry and remote sensing: A review. ISPRS J. Photogramm. Remote Sens. 2014, 92, 79–97. [Google Scholar] [CrossRef] [Green Version]
  5. Nex, F.; Remondino, F. UAV for 3D mapping applications: A review. Appl. Geomat. 2014, 6, 1–15. [Google Scholar] [CrossRef]
  6. Xiao, J.; Gerke, M.; Vosselman, G. Building extraction from oblique airborne imagery based on robust façade detection. ISPRS J. Photogramm. Remote Sens. 2012, 68, 56–68. [Google Scholar] [CrossRef]
  7. Fryskowska, A. Improvement of 3D Power Line Extraction from Multiple Low-Cost UAV Imagery Using Wavelet Analysis. Sensors 2019, 19, 700. [Google Scholar] [CrossRef] [PubMed]
  8. Aasen, H.; Honkavaara, E.; Lucieer, A.; Zarco-Tejada, P. Quantitative remote sensing at ultra-high resolution with uav spectroscopy: A review of sensor technology, measurement procedures, and data correction workflows. Remote Sens. 2018, 10, 1091. [Google Scholar] [CrossRef]
  9. Schaepman, M.E. Spectrodirectional remote sensing: From pixels to processes. Int. J. Appl. Earth Obs. Geoinf. 2007, 9, 204–223. [Google Scholar] [CrossRef] [Green Version]
  10. Honkavaara, E.; Arbiol, R.; Markelin, L.; Martinez, L.; Cramer, M.; Bovet, S.; Chandelier, L.; Ilves, R.; Klonus, S.; Marshal, P.; et al. Digital Airborne Photogrammetry—A New Tool for Quantitative Remote Sensing?—A State-of-the-Art Review On Radiometric Aspects of Digital Photogrammetric Images. Remote Sens. 2009, 1, 577–605. [Google Scholar] [CrossRef]
  11. Berk, A.; Anderson, G.P.; Acharya, P.K.; Dothe, H.; AdlerGolden, S.M.; Richtsmeier, S.C.; Hoke, M.L. MODTRAN4 Radiative Transfer Modeling for Atmospheric Correction; International Society for Optics and Photonics: Denver, CO, USA, 1999; Volume 3756, pp. 348–353. [Google Scholar]
  12. Anderson, G.P.; Berk, A.; Acharya, P.K.; Dothe, H.; Adlergolden, S.M.; Ratkowski, A.J.; Felde, G.W.; Gardner, J.A.; Hoke, M.L.; Richtsmeier, S.C. MODTRAN4: Radiative transfer modeling for remote sensing. In Proceedings of the Algorithms for Multispectral, Hyperspectral, and Ultraspectral Imagery VI, Orlando, FL, USA, 23 August 2000; pp. 1029–1032. [Google Scholar]
  13. Mishra, V.D.; Sharma, J.K.; Singh, K.K.; Thakur, N.K.; Kumar, M. Assessment of different topographic corrections in AWiFS satellite imagery of Himalaya terrain. J. Earth Syst. Sci. 2009, 118, 11–26. [Google Scholar] [CrossRef]
  14. Richter, R. A spatially adaptive fast atmospheric correction algorithm. Int. J. Remote Sens. 1996, 17, 1201–1214. [Google Scholar] [CrossRef]
  15. Yun, S.I.; Kim, K.S. Sky Luminance Measurements Using CCD Camera and Comparisons with Calculation Models for Predicting Indoor Illuminance. Sustainability 2018, 10, 1556. [Google Scholar] [CrossRef]
  16. Yang, G.; Li, C.; Wang, Y.; Yuan, H.; Feng, H.; Xu, B.; Yang, X. The DOM Generation and Precise Radiometric Calibration of a UAV-Mounted Miniature Snapshot Hyperspectral Imager. Remote Sens. 2017, 9, 642. [Google Scholar] [CrossRef]
  17. Nicodemus, F.E. Reflectance nomenclature and directional reflectance and emissivity. Appl. Opt. 1970, 9, 1474–1475. [Google Scholar] [CrossRef] [PubMed]
  18. Aasen, H.; Burkart, A.; Bolten, A.; Bareth, G. Generating 3D hyperspectral information with lightweight UAV snapshot cameras for vegetation monitoring: From camera calibration to quality assurance. ISPRS J. Photogramm. Remote Sens. 2015, 108, 245–259. [Google Scholar] [CrossRef]
  19. Schowengerdt, R.A. Remote Sensing, Models and Methods for Image Processing, 3rd ed.; Academic Press Inc.: San Diego, CA, USA, 2007. [Google Scholar]
  20. Hruska, R.; Mitchell, J.; Anderson, M.; Glenn, N.F. Radiometric and geometric analysis of hyperspectral imagery acquired from an unmanned aerial vehicle. Remote Sens. 2012, 4, 2736–2752. [Google Scholar] [CrossRef]
  21. Honkavaara, E.; Eskelinen, M.A.; Polonen, I.; Saari, H.; Ojanen, H.; Mannila, R.; Holmlund, C.; Hakala, T.; Litkey, P.; Rosnell, T.; et al. Remote Sensing of 3-D Geometry and Surface Moisture of a Peat Production AreaUsing Hyperspectral Frame Cameras in Visible to Short-Wave Infrared Spectral Ranges Onboard a Small Unmanned Airborne Vehicle (UAV). IEEE Trans. Geosci. Remote Sens. 2016, 54, 5440–5454. [Google Scholar] [CrossRef]
  22. Deng, L.; Yan, Y.; Gong, H.; Duan, F.; Zhong, R. The effect of spatial resolution on radiometric and geometric performances of a UAV-mounted hyperspectral 2D imager. ISPRS J. Photogramm. Remote Sens. 2018, 144, 298–314. [Google Scholar] [CrossRef]
  23. Wang, C.; Myint, S.W. A Simplified Empirical Line Method of Radiometric Calibration for Small UnmannedAircraft Systems-Based Remote Sensing. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2015, 8, 1876–1885. [Google Scholar] [CrossRef]
  24. Berni, J.A.J.; Zarco-Tejada, P.J.; Suárez, L.; Fereres, E. Thermal and Narrowband Multispectral Remote Sensing for Vegetation Monitoring from an Unmanned Aerial Vehicle. IEEE Trans. Geosci. Remote Sens. 2009, 47, 722–738. [Google Scholar] [CrossRef]
  25. Tu, Y.-H.; Phinn, S.; Johansen, K.; Robson, A. Assessing Radiometric Correction Approaches for Multi-Spectral UAS Imagery for Horticultural Applications. Remote Sens. 2018, 10, 1684. [Google Scholar] [CrossRef]
  26. Parrot (Ed.) Application Note: Pixel Value to Irradiance Using the Sensor Calibration Model; Parrot: Paris, France, 2017; Volume SEQ-AN-01. [Google Scholar]
  27. Crespi, M.; De Vendictis, L. A Procedure for High Resolution Satellite Imagery Quality Assessment. Sensors 2009, 9, 3289–3313. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Honkavaara, E.; Arbiol, R.; Markelin, L.; Martínez, L.; Cramer, M.; Bovet, S.; Chandelier, L.; Ilves, R.; Klonus, S.; Marshal, P.; et al. The EuroSDR project Radiometric aspects of digital photogrammetric images- Results of the empirical phase. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2011, 38, 123–130, Part 4/W19. (on CD-ROM). [Google Scholar] [CrossRef]
  29. Pyka, K. The use of wavelets for evaluation of loss in radiometric quality in the orthophoto mosaicking process AGH UWND. Ser. Rozpr. Monogr. 2005, 152, 95. (In Polish) [Google Scholar]
  30. Wierzbicki, D.; Fryskowska, A.; Kedzierski, M.; Wojtkowska, M.; Delis, P. Method of radiometric quality assessment of NIR images acquired with a custom sensor mounted on an unmanned aerial vehicle. J. Appl. Remote Sens. 2018, 12, 015008a. [Google Scholar] [CrossRef]
  31. Burdziakowski, P. UAV design and construction for real time photogrammetry and visual navigation. In 2018 Baltic Geodetic Congress (BGC Geomatics); IEEE: Olsztyn, Poland, 2018; pp. 368–372. [Google Scholar] [CrossRef]
  32. Huynh-Thu, Q.; Ghanbari, M. Scope of validity of PSNR in image/video quality assessment. Electron. Lett. 2008, 44, 800. [Google Scholar] [CrossRef]
  33. He, L.; Gao, F.; Hou, W.; Hao, L. Objective image quality assessment: A survey. Int. J. Comput. Math. 2014, 9, 2374–2388. [Google Scholar] [CrossRef]
  34. Jayasankar, U.; Thirumal, V.; Ponnurangam, D. A New Objective Image Quality Assessment Metric: For Color and Grayscale Images. 3d Res. 2018, 93, 28. [Google Scholar] [CrossRef]
  35. Eskicioglu, A.M.; Fisher, P.S. Image quality measures and their performance. IEEE Trans. Commun. 1995, 43, 2959–2965. [Google Scholar] [CrossRef]
  36. Ohm, J. Multimedia Communication Technology: Representation, Transmission and Identification of Multimedia Signals; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
  37. Damera-Venkata, N.; Kite, T.D.; Geisler, W.S.; Evans, B.L.; Bovik, A.C. Image quality assessment based on a degradation model. IEEE Trans. Image Process. 2000, 9, 636–650. [Google Scholar] [CrossRef]
  38. Wang, Y. Survey of Objective Video Quality Measurements; Tech Report; Worcester Polytechnic Insitute: Worcester, MA, USA, 2006. [Google Scholar]
  39. Wang, Z.; Bovik, A.C. A universal image qualityindex. IEEE Signal Process. Lett. 2002, 9, 81–84. [Google Scholar] [CrossRef]
  40. Wang, Z.; Bovik, A.C.; Sheikh, H.R.; Simoncelli, E.P. Image quality assessment: From error measurement to structural similarity. IEEE Trans. Image Process. 2004, 13, 600–612. [Google Scholar] [CrossRef] [PubMed]
  41. Jose, S. Topography in Europe, map, based on the Generic Mapping Tools and ETOPO2. Available online: https://pl.wikipedia.org/wiki/Plik:Europe_topography_map.png (accessed on 20 May 2019).
  42. Synoptic Maps. Available online: http://www.pogodynka.pl/polska/mapa_synoptyczna (accessed on 20 May 2019).
  43. The Hydrometeorological Service of the Polish Armed Force, Poland area military significant weather chart (low level). Available online: http:// www.mete.mil.pl (accessed on 20 May 2019).
  44. Kedzierski, M.; Wierzbicki, D. Radiometric quality assessment of images acquired by UAV’s in various lighting and weather conditions. Measurement 2015, 76, 156–169. [Google Scholar] [CrossRef]
  45. Wierzbicki, D.; Kedzierski, M.; Fryskowska, A.; Jasinski, J. Quality Assessment of the Bidirectional Reflectance Distribution Function for NIR Imagery Sequences from UAV. Remote Sens. 2018, 10, 1348. [Google Scholar] [CrossRef]
  46. Benford, F.; Back, J.E. A time analysis of Sunshine. In Trans. of American; Illumination Engineering Society: New York, NY, USA, 1939; Volume 34, pp. 200–218. [Google Scholar]
  47. Jacobson, M.Z. Fundamentals of Atmospheric Modeling; Cambridge University Press: Cambridge, UK, 2005; p. 317. ISBN 9780521548656. [Google Scholar]
  48. Chwieduk, D.; Bogdanska, B. Some recommendations for inclinations and orientations of building elements under solar radiation in Polish conditions. Renew. Energy 2004, 29, 1569–1581. [Google Scholar] [CrossRef]
  49. Richter, R.; Schläpfer, D. Geo-atmospheric processing of airborne imaging spectrometry data. Part 2: Atmospheric/topographic correction. Int. J. Remote Sens. 2002, 23, 2631–2649. [Google Scholar] [CrossRef]
  50. Yu, X.; Liu, Q.; Liu, X.; Liu, X.; Wang, Y. A physical-based atmospheric correction algorithm of unmanned aerial vehicles images and its utility analysis. Int. J. Remote Sens. 2017, 38, 3101–3112. [Google Scholar] [CrossRef]
  51. Song, Y.; Luo, H.; Lu, R.; Ma, J. Dehazed image quality assessment by haze-line theory. J. Phys. Conf. Ser. 2017, 844, 012045. [Google Scholar] [CrossRef]
  52. Yang, Y.; Lin, Z.; Liu, F. Stable imaging and accuracy issues of low-altitude unmanned aerial vehicle photogrammetry systems. Remote Sens. 2016, 8, 316. [Google Scholar] [CrossRef]
Figure 1. Localization of study area (left image based on Reference [41]).
Figure 1. Localization of study area (left image based on Reference [41]).
Remotesensing 11 01214 g001
Figure 2. Placement of the calibration panels (from P1 to P6) during image acquisition.
Figure 2. Placement of the calibration panels (from P1 to P6) during image acquisition.
Remotesensing 11 01214 g002
Figure 3. M10 radiosonde.
Figure 3. M10 radiosonde.
Remotesensing 11 01214 g003
Figure 4. AGAT20L weather station.
Figure 4. AGAT20L weather station.
Remotesensing 11 01214 g004
Figure 5. Europe synoptic map, September 13, 2018, at 12:00 UTC (based on Reference [42]).
Figure 5. Europe synoptic map, September 13, 2018, at 12:00 UTC (based on Reference [42]).
Remotesensing 11 01214 g005
Figure 6. Graphical forecast for Poland for September 13, 2018, at 06:00 and 12:00 UTC (±3 h) (based on Reference [43]).
Figure 6. Graphical forecast for Poland for September 13, 2018, at 06:00 and 12:00 UTC (±3 h) (based on Reference [43]).
Remotesensing 11 01214 g006
Figure 7. Europe synoptic map November 22, 2018, at 12:00 UTC (based on Reference [42]).
Figure 7. Europe synoptic map November 22, 2018, at 12:00 UTC (based on Reference [42]).
Remotesensing 11 01214 g007
Figure 8. Graphical forecast for Poland for November 22, 2018, at 06:00 and 12:00 UTC (±3 h) (based on Reference [43]).
Figure 8. Graphical forecast for Poland for November 22, 2018, at 06:00 and 12:00 UTC (±3 h) (based on Reference [43]).
Remotesensing 11 01214 g008
Figure 9. Europe synoptic map September 26, 2018, at 12:00 UTC (based on Reference [42]).
Figure 9. Europe synoptic map September 26, 2018, at 12:00 UTC (based on Reference [42]).
Remotesensing 11 01214 g009
Figure 10. Graphical forecast for Poland for September 26, 2018, at 06:00, 12:00, and 18:00 UTC (±3 h) (based on Reference [43]).
Figure 10. Graphical forecast for Poland for September 26, 2018, at 06:00, 12:00, and 18:00 UTC (±3 h) (based on Reference [43]).
Remotesensing 11 01214 g010
Figure 11. Correlation coefficient (CC) for humidity change function and electromagnetic wave reflection coefficient during the day at different unmanned aerial vehicle (UAV) flight altitudes for the calibration test with known spectral characteristics in the Mieruniszki test area (obtained September 13, 2018).
Figure 11. Correlation coefficient (CC) for humidity change function and electromagnetic wave reflection coefficient during the day at different unmanned aerial vehicle (UAV) flight altitudes for the calibration test with known spectral characteristics in the Mieruniszki test area (obtained September 13, 2018).
Remotesensing 11 01214 g011
Figure 12. Scheme of low atmosphere scattering model for UAV imaging.
Figure 12. Scheme of low atmosphere scattering model for UAV imaging.
Remotesensing 11 01214 g012
Figure 13. Scheme of solar radiation falling on a horizontal surface; α is the elevation angle of the sun, θz is the solar zenith angle, and γS is the azimuthal angle of the sun.
Figure 13. Scheme of solar radiation falling on a horizontal surface; α is the elevation angle of the sun, θz is the solar zenith angle, and γS is the azimuthal angle of the sun.
Remotesensing 11 01214 g013
Figure 14. Distribution of the QA indicator values for reference data acquired (left to right): morning, noon, and evening, at (a) 50 m, (b) 100 m, and (c) 150 m.
Figure 14. Distribution of the QA indicator values for reference data acquired (left to right): morning, noon, and evening, at (a) 50 m, (b) 100 m, and (c) 150 m.
Remotesensing 11 01214 g014
Figure 15. Sample photos obtained (a) September 26, 2018, in Kościelisko, and (b) November 22, 2018, in Mieruniszki.
Figure 15. Sample photos obtained (a) September 26, 2018, in Kościelisko, and (b) November 22, 2018, in Mieruniszki.
Remotesensing 11 01214 g015
Figure 16. Air humidity and residuals from the theoretical value of reflectance for 2 and 95% tests for data obtained at (a) 125 m and (b) 200 m in good weather conditions (Kościelisko area, September 26, 2018).
Figure 16. Air humidity and residuals from the theoretical value of reflectance for 2 and 95% tests for data obtained at (a) 125 m and (b) 200 m in good weather conditions (Kościelisko area, September 26, 2018).
Remotesensing 11 01214 g016
Figure 17. Air humidity and residuals from the theoretical value of reflectance for 2 and 95% tests for data obtained at (a) 125 m and (b) 200 m in good weather conditions (Mieruniszki area, September 13, 2018).
Figure 17. Air humidity and residuals from the theoretical value of reflectance for 2 and 95% tests for data obtained at (a) 125 m and (b) 200 m in good weather conditions (Mieruniszki area, September 13, 2018).
Remotesensing 11 01214 g017
Figure 18. Air humidity and residuals from the theoretical value of reflectance for 2% and 95% tests for data obtained at (a) 125 m and (b) 200 m in poor weather conditions (Mieruniszki area, November 22, 2018).
Figure 18. Air humidity and residuals from the theoretical value of reflectance for 2% and 95% tests for data obtained at (a) 125 m and (b) 200 m in poor weather conditions (Mieruniszki area, November 22, 2018).
Remotesensing 11 01214 g018
Table 1. Reflectance measurements for each band of the DJI (DJI, Shenzhen, China) camera for panels P1 to P6.
Table 1. Reflectance measurements for each band of the DJI (DJI, Shenzhen, China) camera for panels P1 to P6.
DJP1P2P3P4P5P6
Blue0.0220.0860.1390.2090.2890.363
Green0.020.0770.1250.1880.2630.342
Red0.0270.1050.1680.2470.3310.41
Table 2. Classification of radiometric quality by QA coefficient values.
Table 2. Classification of radiometric quality by QA coefficient values.
Image Radiometric QualityQA Coefficient Intervals
Good radiometric quality Q A [ 0.00 ;   6.00 )
Medium radiometric quality Q A [ 6.00 ;   7.65 )
Bad radiometric quality Q A   7.65
Table 3. Case study of image quality depending on solar zenith angle and humidity.
Table 3. Case study of image quality depending on solar zenith angle and humidity.
Solar Zenith Angle
Humidity14°38°
80%18.46.62.6
40%9.23.31.3
Table 4. QA indicator statistics for the verification sets of UAV imagery data.
Table 4. QA indicator statistics for the verification sets of UAV imagery data.
Area NameAltitude (Above Ground Level) (m)Morning 06:00–07:00 UTCBefore Noon 09:00–12:00 UTCAfter Noon 13:00–14:00 UTC
MeanSDMinMaxMeanSDMinMaxMeanSDMinMax
Mieruniszki 09/13/2018754.51.52.18.81.50.50.73.92.60.61.84.4
1004.00.82.46.11.80.40.83.42.80.42.23.9
1504.40.72.75.22.10.40.92.82.90.61.42.9
2006.51.72.28.12.30.80.84.13.20.80.55.5
3006.42.31.89.22.50.61.13.63.31.01.46.2
Kościelisko 09/26/2018755.72.71.17.71.30.90.95.61.21.21.36.9
1006.42.80.47.81.40.90.44.11.20.80.27.7
1507.71.50.59.51.61.10.33.61.31.90.45.7
2007.62.11.39.31.61.30.24.41.30.90.94.9
3005.41.70.98.71.40.50.63.41.40.80.34.3
Mieruniszki 11/22/20187524.119.40.674.16.14.50.522.232.18.70.946.5
10025.523.50.389.36.77.10.954.633.617.91.358.2
15029.919.42.278.87.36.71.521.738.832.75.8112.0
20033.116.71.257.27.93.10.915.543.623.13.1152.3
30047.527.21.4167.49.33.90.823.757.019.95.2144.2

Share and Cite

MDPI and ACS Style

Kedzierski, M.; Wierzbicki, D.; Sekrecka, A.; Fryskowska, A.; Walczykowski, P.; Siewert, J. Influence of Lower Atmosphere on the Radiometric Quality of Unmanned Aerial Vehicle Imagery. Remote Sens. 2019, 11, 1214. https://doi.org/10.3390/rs11101214

AMA Style

Kedzierski M, Wierzbicki D, Sekrecka A, Fryskowska A, Walczykowski P, Siewert J. Influence of Lower Atmosphere on the Radiometric Quality of Unmanned Aerial Vehicle Imagery. Remote Sensing. 2019; 11(10):1214. https://doi.org/10.3390/rs11101214

Chicago/Turabian Style

Kedzierski, Michal, Damian Wierzbicki, Aleksandra Sekrecka, Anna Fryskowska, Piotr Walczykowski, and Jolanta Siewert. 2019. "Influence of Lower Atmosphere on the Radiometric Quality of Unmanned Aerial Vehicle Imagery" Remote Sensing 11, no. 10: 1214. https://doi.org/10.3390/rs11101214

APA Style

Kedzierski, M., Wierzbicki, D., Sekrecka, A., Fryskowska, A., Walczykowski, P., & Siewert, J. (2019). Influence of Lower Atmosphere on the Radiometric Quality of Unmanned Aerial Vehicle Imagery. Remote Sensing, 11(10), 1214. https://doi.org/10.3390/rs11101214

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop