Lookup Table Approach for Radiometric Calibration of Miniaturized Multispectral Camera Mounted on an Unmanned Aerial Vehicle

: Over recent years, miniaturized multispectral cameras mounted on an unmanned aerial vehicle (UAV) have been widely used in remote sensing. Most of these cameras are integrated with low-cost, image-frame complementary metal-oxide semiconductor (CMOS) sensors. Compared to the typical charged coupled device (CCD) sensors or linear array sensors, consumer-grade CMOS sensors have the disadvantages of low responsivity, higher noise, and non-uniformity of pixels, which make it di ﬃ cult to accurately detect optical radiation. Therefore, comprehensive radiometric calibration is crucial for quantitative remote sensing and comparison of temporal data using such sensors. In this study, we examine three procedures of radiometric calibration: relative radiometric calibration, normalization, and absolute radiometric calibration. The complex features of dark current noise, vignetting effect, and non-uniformity of detector response are analyzed. Further, appropriate procedures are used to derive the lookup table (LUT) of correction factors for these features. Subsequently, an absolute calibration coe ﬃ cient based on an empirical model is used to convert the digital number (DN) of images to radiance unit. Due to the radiometric calibration, the DNs of targets observed in the image are more consistent than before calibration. Compared to the method provided by the manufacturer of the sensor, LUTs facilitate much better radiometric calibration. The root mean square error (RMSE) of measured reﬂectance in each band (475, 560, 668, 717, and 840 nm) are 2.30%, 2.87%, 3.66%, 3.98%, and 4.70% respectively.


Introduction
Over recent years, miniaturized multispectral cameras based on cost-effective sensors have been extensively applied in the remote sensing field [1][2][3][4][5][6], such as MicaSense RedEdge-MX (MicaSense Inc., Seattle, WA, USA), MS600 pro (Changguang YuSense Information Technology and Equipment (Qingdao) Co., Ltd., Qingdao, China), and Parrot Sequoia (Parrot Drone SAS, Paris, France). Such multispectral cameras are usually mounted on unmanned aerial vehicles (UAVs), which facilitates individual researchers, small teams, or commercial industries to obtain spectral information with ultra-high spatial resolution [7]. From the perspective of a low-cost and lightweight device, consumer-grade complementary metal-oxide semiconductor (CMOS) sensors are the main choice for integrating such  Typically, MicaSense RedEdge-MX is used to generate precise and quantitative information on the vigor and health of crops. For aviation applications at low altitudes, the integration time should not be too high or too low to avoid insufficient exposure and the affect of image motion on the geometric accuracy. An integration time in the range of 1/2000-1/500 s is an appropriate and common choice [19,22,33]. Typically, MicaSense RedEdge-MX is used to generate precise and quantitative information on the vigor and health of crops. For aviation applications at low altitudes, the integration time should not be too high or too low to avoid insufficient exposure and the affect of image motion on the geometric accuracy. An integration time in the range of 1/2000-1/500 s is an appropriate and common choice [19,22,33].

Integrating Sphere System
The radiation characteristics of the sensor were observed using an integrating sphere (XTH2000, Labsphere Inc., North Sutton, NH, USA). The structure of the integrating sphere system and the measurement scheme are shown in Figure 2. The integrating sphere provided stable and uniform optical radiation in the wavelength range of 300-2400 nm. Its aperture was 20 cm, and the optical uniformity was more than 98%. Typically, MicaSense RedEdge-MX is used to generate precise and quantitative information on the vigor and health of crops. For aviation applications at low altitudes, the integration time should not be too high or too low to avoid insufficient exposure and the affect of image motion on the geometric accuracy. An integration time in the range of 1/2000-1/500 s is an appropriate and common choice [19,22,33].

Integrating Sphere System
The radiation characteristics of the sensor were observed using an integrating sphere (XTH2000, Labsphere Inc., North Sutton, NH, USA). The structure of the integrating sphere system and the measurement scheme are shown in Figure 2. The integrating sphere provided stable and uniform optical radiation in the wavelength range of 300-2400 nm. Its aperture was 20 cm, and the optical uniformity was more than 98%.

Radiometric Calibration Model
The DN recorded by the detectors of CMOS sensors looking at a target with emitting radiance L(x, y) were modeled as follows [5,[12][13][14]34]: L(x, y) = a * 1 g * t e * 2 N * where x, y was the position of the target on the image; R(x, y) was the correction factor at x, y in the LUT of response; V(x, y) was the correction factor at x, y in the LUT of vignetting effect; Dark(x, y) was the correction factor at x, y in the LUT of dark current; DN(x, y) was the raw pixel value at x, y; g and t e were the gain and integration time of sensors, respectively; 2 N was the digital bit rate (8 bit, 12 bit, 16 bit); N was the number of bits; a and b were the quantization parameters of analog-to-digital conversion, namely quantized coefficient and offset of absolute radiometric calibration, respectively. The radiometric calibration process was divided into three steps: relative radiometric calibration, data normalization, and absolute radiometric calibration. Equation (1) could be simplified by decomposing it into three parts.
Firstly, relative radiometric calibration made DN values comparable within the image. The main task of the relative radiometric calibration was to eliminate the influence of dark current, vignetting effect, and non-uniform detection efficiency [2,13,17,18,35].

of 25
Secondly, data normalization made the images, which were captured under different camera parameter settings, comparable [3,14,33].
Based on Equations (2)-(4), DN could be converted to the radiance unit only if the LUTs of Dark(x, y), V(x, y), R(x, y), a, b were derived. In the next step, we adopted reasonable solutions to obtain these parameter LUTs in experiments so that the incident radiation field for an arbitrary operational image could be determined.

Correction of Dark Current
The correction parameters of dark current were defined in terms of the LUT of non-uniformity in 2D image space rather than the mean value of dark current in the image [7,9,17]. To reduce the complexity of the dark current correction, a representative template was selected as the LUT of dark current. This LUT was an effective choice to eliminate dark current for images taken under different camera settings of gain and integration time. The dark current was characterized by gain and integration time to determine the representative template after the temperature of the device was stabilized. The LUT of dark current can be obtained as following procedure in Figure 3.

Vignetting Effect Correction
Vignetting is a phenomenon in which the center of the image is brighter and the edge area is darker due to the limitation of apertures and defects of optical devices [26,31,37]. Polynomial fitting based on radial distance is commonly utilized for vignetting correction of miniaturized multispectral cameras [25,38,39]. Due to the imperfect manufacturing process, the brightness center is not consistent with the sensor center and is difficult to determine. The vignetting effect on 2D images is more complex and may be asymmetric. In practice, there can be deviations in the polynomial fitting. The vignetting correction method based on image LUT has high correction accuracy [17,25,29]. Each pixel has a unique correction factor in the LUT of vignetting correction.
Theoretically, a vignette is a gentle gradient from the middle to the edge [40]. Due to the nonuniform response and the influence of random noise, the images of uniform objects are not smooth. Thus, in an image obtained with uniform incident light, the vignetting effect corresponds to low frequency, while random noise and response noise are high-frequency noises. Here, the Gaussian low-pass filtering method was selected to eliminate random noise and response noise and extract background brightness from the vignetting images. The filter could prevent distortion of the edge signal in the smoothing process [15,25].

Vignetting Effect Correction
Vignetting is a phenomenon in which the center of the image is brighter and the edge area is darker due to the limitation of apertures and defects of optical devices [26,31,37]. Polynomial fitting based on radial distance is commonly utilized for vignetting correction of miniaturized multispectral cameras [25,38,39]. Due to the imperfect manufacturing process, the brightness center is not consistent with the sensor center and is difficult to determine. The vignetting effect on 2D images is more complex and may be asymmetric. In practice, there can be deviations in the polynomial fitting. The vignetting correction method based on image LUT has high correction accuracy [17,25,29]. Each pixel has a unique correction factor in the LUT of vignetting correction.
Theoretically, a vignette is a gentle gradient from the middle to the edge [40]. Due to the non-uniform response and the influence of random noise, the images of uniform objects are not smooth. Thus, in an image obtained with uniform incident light, the vignetting effect corresponds to low frequency, while random noise and response noise are high-frequency noises. Here, the Gaussian low-pass filtering method was selected to eliminate random noise and response noise and extract background brightness from the vignetting images. The filter could prevent distortion of the edge signal in the smoothing process [15,25].
Then, the vignetting image was normalized through the correction factors. v(x, y) was the pixel value in the vignetting image; V max . was the maximum pixel value within the vignetting image.
Several normalized vignetting images were obtained at different incident radiation levels. The effective LUT of the vignetting correction represented the mean image of several normalized vignetting images.

Response Correction
After correction of dark current and vignetting effect, images captured under a uniform illumination retained the error of non-uniform pixel response and random noise. Under standard and uniform illumination, the photoelectric response of each detector in CMOS was inconsistent due to the imperfect manufacturing process [6,34]. Here, multiple images under the same radiance of the integrating sphere were corrected for dark current and vignetting effect. Then, these corrected images were averaged to eliminate random noise.
The correction factor of the response function was defined as the ratio of pixel response value to the average response value. The pixel response value and average response value were obtained after correcting the image for dark current and vignetting effect.
Similarly, the effective LUT of response correction represented the mean image of several corrected images obtained at different incident radiation levels.

Parameters of Absolute Radiometric Calibration
According to Equations (3) and (4), the linear fitting method was used to derive the coefficient and offset. The dataset used for linear fitting was collected at the different radiation intensities and camera settings of gain and integration time. This dataset must be corrected using relative radiometric calibration beforehand. Linear fitting indicated a linear causal relationship between two variables, i.e., where i was the number of observed value; a and b were the fitting coefficient and offset, respectively; aX i + b expressed the linear relationship between Y and X. ε i was a random variable that reflected the distribution around the straight line of statistical relationship and obeyed normal distribution. a and b were statistically analyzed based on the sample dataset. The evaluated values of a and b were α and β, which obeyed the following linear equation:Ŷ Generally, α and β should be calculated to minimize the deviation between each sample observation point (Y i , X i ) and the fitting line.

Experimental Schemes
The radiation calibration experiment was divided into three parts: dark current experiment, uniformity experiment, and absolute calibration experiment.
(1) Dark current experiment: These experiments were conducted with a lens cover to simulate the environment without incident light. The variations in the dark current as a function of gain (1×, 2×, 4×, and 8×) and integration time (0.44, 0.59, 0.78, 1.0, 1.4, 1.9, and 2.5) were examined. Based on the method described in Section 2.2.2, the correction factor LUT for the non-uniformity of dark current was obtained. (2) Uniformity experiment: The quantized values of the image captured in the uniform light field of the integrating sphere were non-uniform due to the vignetting effect and differences in pixel response. Firstly, the dark current in the images was eliminated using the corresponding LUT. Then, the vignetting effect was analyzed based on Section 2.2.3. Subsequently, the difference in pixel response was examined using the images that had already been corrected for dark current and vignetting effect based on Section 2. All the images were corrected for dark current, vignetting effect, and response according to the previous two steps. Then, the average value of the image was calculated as the effective quantized value to establish a linear relationship with incident radiance by linear regression.
The detailed experimental scheme is shown in Figure 4. In the experiment, the images were collected in 16-bit TIFF format, and the spectral radiance of the integrating sphere was recorded synchronously. The data were processed using MATLAB R2019a software (The MathWorks, Inc., Natick, MA, USA).
Remote Sens. 2020, 12, x 8 of 26 The detailed experimental scheme is shown in Figure 4. In the experiment, the images were collected in 16-bit TIFF format, and the spectral radiance of the integrating sphere was recorded synchronously. The data were processed using MATLAB R2019a software (The MathWorks, Inc., Natick, MA, USA).

Correction of Dark Current
To extract an optimal template image of dark current, dark current was analyzed in terms of gain and integration time. It is clear from Figure 5 that the DN of typical pixels did not increase with the increase in gain. Actually, there were no obvious changes in the dark current value of typical pixels for gain below 2×. The dark current values of the pixels were different when the gain was set to 4× and 8×. The value of some pixels increased, while that of other pixels decreased.
Remote Sens. 2020, 12, x 9 of 26 inconsistent the dark current value. Thus, to limit the effect of noise, the gain of the sensor should be set to a low level of 1× or 2× in practice. The influence of gain and integration time (0.44-2.5 ms) on the dark current was experimentally analyzed, and the corresponding results are shown in Figures 7 and 8, respectively. It is clear from  Figure 6 shows the statistical characteristics of dark current images. As the gain changed, the histogram of the dark current image appeared obviously different. The larger the gain, the more inconsistent the dark current value. Thus, to limit the effect of noise, the gain of the sensor should be set to a low level of 1× or 2× in practice. The influence of gain and integration time (0.44-2.5 ms) on the dark current was experimentally analyzed, and the corresponding results are shown in Figures 7 and 8, respectively. It is clear from Figure 7 that the dark current value of typical pixels in each band did not change with the increase in integration time. This suggested that the image template obtained at any integration time could be used to correct the dark current of the image captured at other integration time settings. The influence of gain and integration time (0.44-2.5 ms) on the dark current was experimentally analyzed, and the corresponding results are shown in Figures 7 and 8, respectively. It is clear from Figure 7 that the dark current value of typical pixels in each band did not change with the increase in integration time. This suggested that the image template obtained at any integration time could be used to correct the dark current of the image captured at other integration time settings.
Remote Sens. 2020, 12, x 10 of 26  The histograms of dark current images were slightly changed under different values of integration time. To clearly display the statistical features of images, the standard deviation (STD) of the dark current image was selected to quantify this change. The STD of the dark current image in the integration time range of 0.44-2.5 ms, which is typically used for a multispectral camera in aviation applications, was statistically analyzed, and the results are shown in Figure 9. The histograms of dark current images were slightly changed under different values of integration time. To clearly display the statistical features of images, the standard deviation (STD) of the dark current image was selected to quantify this change. The STD of the dark current image in the integration time range of 0.44-2.5 ms, which is typically used for a multispectral camera in aviation applications, was statistically analyzed, and the results are shown in Figure 9. From an integration time in the range of 0.44-2.5 ms, the STD was negligibly changed. This indicated there was no difference in the dark current of images captured between 0.44 and 2.5 ms.
The above experimental analysis provided guidance for selecting a representative dark current image as the correction factor LUT of Dark(x, y). Besides, the integration time of 1 ms is commonly used in practical application as it lies in the typical range used for aerial photography [19,22,33]. Consequently, the LUT of dark current Dark(x, y) in the radiometric calibration model was programmed as a typical image, which was the average of multiple images captured using gain of 1× and integration time of 1 ms.

Vignetting Effect
The vignetting effect of the image in which the effect of dark current had been eliminated is shown in Figure 10. The DN value decreased from the middle to the edge in the images captured by looking at the pupil of the integrating sphere, which output a uniform light field. Due to the vignetting effect, the edge brightness of the image was roughly 0.60-0.77 times the middle brightness in all the bands. The maximum luminance difference was 0.60 times in Band 5. From an integration time in the range of 0.44-2.5 ms, the STD was negligibly changed. This indicated there was no difference in the dark current of images captured between 0.44 and 2.5 ms.
The above experimental analysis provided guidance for selecting a representative dark current image as the correction factor LUT of Dark(x, y). Besides, the integration time of 1 ms is commonly used in practical application as it lies in the typical range used for aerial photography [19,22,33]. Consequently, the LUT of dark current Dark(x, y) in the radiometric calibration model was programmed as a typical image, which was the average of multiple images captured using gain of 1× and integration time of 1 ms.

Vignetting Effect
The vignetting effect of the image in which the effect of dark current had been eliminated is shown in Figure 10. The DN value decreased from the middle to the edge in the images captured by looking at the pupil of the integrating sphere, which output a uniform light field. Due to the vignetting effect, the edge brightness of the image was roughly 0.60-0.77 times the middle brightness in all the bands. The maximum luminance difference was 0.60 times in Band 5.
Based on Equations (5) and (6), the Gaussian filter template was created by inputting the appropriate STD σ. The convolution algorithm was used between the image and Gaussian filter ( f ilter gsussian ) template to generate the vignetting images. The fitting performance of all the bands are shown in Figure 11.
Then, the correction factor LUT of the vignetting effect was obtained by normalizing these images according to Equation (7). The normalized LUT of the vignetting effect of all the bands are displayed in Figure 12.
The correction factor LUT of the vignetting effect was used for image correction, and the results are shown in Figure 13. Compared to Figure 10, the brightness of the corrected image was even, which proved that the vignetting effect was eliminated. From an integration time in the range of 0.44-2.5 ms, the STD was negligibly changed. This indicated there was no difference in the dark current of images captured between 0.44 and 2.5 ms.
The above experimental analysis provided guidance for selecting a representative dark current image as the correction factor LUT of Dark(x, y). Besides, the integration time of 1 ms is commonly used in practical application as it lies in the typical range used for aerial photography [19,22,33]. Consequently, the LUT of dark current Dark(x, y) in the radiometric calibration model was programmed as a typical image, which was the average of multiple images captured using gain of 1× and integration time of 1 ms.

Vignetting Effect
The vignetting effect of the image in which the effect of dark current had been eliminated is shown in Figure 10. The DN value decreased from the middle to the edge in the images captured by looking at the pupil of the integrating sphere, which output a uniform light field. Due to the vignetting effect, the edge brightness of the image was roughly 0.60-0.77 times the middle brightness in all the bands. The maximum luminance difference was 0.60 times in Band 5. Based on Equations (5) and (6), the Gaussian filter template was created by inputting the appropriate STD σ. The convolution algorithm was used between the and Gaussian filter ( ) template to generate the vignetting images. The fitting performance of all the bands are shown in Figure 11.  Based on Equations (5) and (6), the Gaussian filter template was created by inputting the appropriate STD σ. The convolution algorithm was used between the and Gaussian filter ( ) template to generate the vignetting images. The fitting performance of all the bands are shown in Figure 11.  Then, the correction factor LUT of the vignetting effect was obtained by normalizing these images according to Equation (7). The normalized LUT of the vignetting effect of all the bands are displayed in Figure 12. The correction factor LUT of the vignetting effect was used for image correction, and the results are shown in Figure 13. Compared to Figure 10, the brightness of the corrected image was even, which proved that the vignetting effect was eliminated. Then, the correction factor LUT of the vignetting effect was obtained by normalizing these images according to Equation (7). The normalized LUT of the vignetting effect of all the bands are displayed in Figure 12. The correction factor LUT of the vignetting effect was used for image correction, and the results are shown in Figure 13. Compared to Figure 10, the brightness of the corrected image was even, which proved that the vignetting effect was eliminated.

Response Correction
The images obtained under the same output radiation intensity of the integration sphere and camera settings were corrected for dark current and vignetting effect. Their mean image was calculated to eliminate the random noise. The response differences that were retained in the images are shown in Figure 14. An obvious stripe noise was observed in the images of each band.

Response Correction
The images obtained under the same output radiation intensity of the integration sphere and camera settings were corrected for dark current and vignetting effect. Their mean image was calculated to eliminate the random noise. The response differences that were retained in the images are shown in Figure 14. An obvious stripe noise was observed in the images of each band. According to Equation (8), the correction factor LUT of the response function was calculated, which was then used to correct the images that were already corrected for dark current and vignetting effect. The results are shown in Figure 15. According to Equation (8), the correction factor LUT of the response function was calculated, which was then used to correct the images that were already corrected for dark current and vignetting effect. The results are shown in Figure 15.
The results suggested that DN values of images with response factor correction became smoother and more compact than before. This indicated that the difference in response between pixels had been effectively corrected. In the residual error, random noise was the main component. According to Equation (8), the correction factor LUT of the response function was calculated, which was then used to correct the images that were already corrected for dark current and vignetting effect. The results are shown in Figure 15. The results suggested that DN values of images with response factor correction became smoother and more compact than before. This indicated that the difference in response between pixels had been effectively corrected. In the residual error, random noise was the main component.

Absolute Calibration Parameters
We collected 21 groups of data for each band. The datasets were measured under three incident brightness levels (the output power of the integrating sphere was 100%, 70%, and 40%) and seven integration times (0.44, 0.59, 0.78, 1.0, 1.4, 1.9, and 2.5 ms). The gain was set to 1×. Each group of data

Absolute Calibration Parameters
We collected 21 groups of data for each band. The datasets were measured under three incident brightness levels (the output power of the integrating sphere was 100%, 70%, and 40%) and seven integration times (0.44, 0.59, 0.78, 1.0, 1.4, 1.9, and 2.5 ms). The gain was set to 1×. Each group of data included five images captured continuously under the same condition. Five images in one group were averaged as one image to eliminate the random noise. Thus, 21 effective images were used for absolute radiometric calibration processing. The process is described as follows: • Firstly, according to Equation (2), all the collected images were corrected for dark current, vignetting effect, and non-uniform pixel response with the LUTs derived above.

•
Secondly, according to Equation (3), the 21 images were normalized with respect to the integration time, gain, and digital bit rate, and the mean DN was calculated for each corrected image.

•
Finally, the parameters of absolute radiometric calibration, i.e., a and b, were obtained by least-squares fitting based on the mean DN of 21 corrected images and the spectral radiance output of the integration sphere [6,18,33].
The linear fitting results are shown in Figure 16, which proved the linear relationship between the data and radiance. The coefficient a, offset b, R-square, root mean square error (RMSE) and mean square (MS) of residual of linear fitting are displayed in Table 2. The R-squares of all the bands were remarkably close to 1, and the MS of residual of all the bands were close to 0. This indicated that the forecast precision of linear precision was good. In other words, the obtained values of coefficient a and offset b in absolute calibration were valid.

Verification Method
To verify the effectiveness of the proposed radiometric calibration model, the relative accuracy and absolute accuracy of the model were evaluated. The relative accuracy was used to validate the uniformity of the corrected image, while the absolute accuracy was used to validate the radiance in images.

Verification Method
To verify the effectiveness of the proposed radiometric calibration model, the relative accuracy and absolute accuracy of the model were evaluated. The relative accuracy was used to validate the uniformity of the corrected image, while the absolute accuracy was used to validate the radiance in images.

Verification Method
To verify the effectiveness of the proposed radiometric calibration model, the relative accuracy and absolute accuracy of the model were evaluated. The relative accuracy was used to validate the uniformity of the corrected image, while the absolute accuracy was used to validate the radiance in images.
Five observation targets were combined: calibrated reflectance panels (CRP) provided by the manufacturer, white target01 (WT01), white target02 (WT02), gray target01 (GT01), and gray target02 (GT02). As shown in Figure 17a, the observation targets were arranged on flat ground.  The reflectance of five targets was measured using a RS-8800 field spectrometer (Spectral Evolution, USA) with a spectral range of 350-2500 nm. The target was observed from multiple angles. The average reflectance obtained through multi-angle observations was considered as an effective reflectance. The reflectance of targets at the spectral bands consistent with the center wavelength of RedEdge-MX are shown in Table 3. The experiment was conducted outdoors in cloudless conditions. According to the experimental scheme shown in Figure 17b, five spectral images were captured using a Mica Rededge-MX multispectral camera fixed on a tripod with a height of 1.5 m from the target such that the viewing angle was perpendicular to the ground. The Band1 image of targets is shown in Figure 18a as an example. The diagonal pixels are extracted in Figure 18b   Further, the radiometric calibration method provided by the manufacturer of Rededge-MX [41] was also used to correct the same data and compared it with the proposed LUT-based method. The comparison results are discussed in the following subsections.

Relative Accuracy
After relative correction, the result was analyzed from three aspects: diagonal pixels of images, mean DN of targets, and STD of DN of targets. It was expected that the values of pixels should be more balanced than before in a single target and between similar targets. Firstly, the diagonal pixels of original images and the images corrected by two methods were extracted, and the corresponding results are shown in Figure 19.
It is clear from Figure 19 that the diagonal pixels of all the spectral images corrected by the two methods were smaller than that of the original spectral images because the dark current was eliminated. For one target, the pixels were flatter in the images corrected by the two methods. For similar targets, the consistency of DN values was also improved. In addition, Figure 19b shows that the effect of the two methods was roughly similar. However, the flatness of the pixel curve corrected by the LUT method was better than that of the curve corrected by the manufacturer's method for the Further, the radiometric calibration method provided by the manufacturer of Rededge-MX [41] was also used to correct the same data and compared it with the proposed LUT-based method. The comparison results are discussed in the following subsections.

Relative Accuracy
After relative correction, the result was analyzed from three aspects: diagonal pixels of images, mean DN of targets, and STD of DN of targets. It was expected that the values of pixels should be more balanced than before in a single target and between similar targets. Firstly, the diagonal pixels of original images and the images corrected by two methods were extracted, and the corresponding results are shown in Figure 19.
It is clear from Figure 19 that the diagonal pixels of all the spectral images corrected by the two methods were smaller than that of the original spectral images because the dark current was eliminated. For one target, the pixels were flatter in the images corrected by the two methods. For similar targets, the consistency of DN values was also improved. In addition, Figure 19b shows that the effect of the two methods was roughly similar. However, the flatness of the pixel curve corrected by the LUT method was better than that of the curve corrected by the manufacturer's method for the WT01 target and WT02 target and between them. This implied that the proposed LUT method facilitated better relative calibration than the manufacturer's method.
For quantitative comparison, the mean and STD of DN of each observed targets in the original images and corrected images based on the two calibration methods are shown in Figures 20 and 21.
It is evident from Figure 20 that the interval between the mean of the DN of WT01 and that of WT02 in the corrected results increased as compared to that in the original images. This may be attributed to the tilt observed by the sensor. However, the mean of DN of other similar targets became closer than before. This indicated the effectiveness of the two calibration methods. Moreover, the mean DN of all the same targets corrected by the LUT method was more uniform than that of the same targets corrected by the manufacturer's method, as shown in Figure 20. This implied that the effect of the LUT method was better than that of the manufacturer's method.
The STD of the DN of the WT02 target in the corrected images based on the two methods was larger than that of the original images. This also indicated the problem of sensor tilt similar to Figures 19a and 20. However, the STD of other targets became smaller than that in the original images. This confirmed that the uniformity of the image was improved. Besides, the STD of all the targets corrected by the LUT method was smaller than that of the targets corrected by the manufacturer's method. This indicated that the images corrected by the LUT method had better uniformity than those corrected by the manufacturer's method.
Remote Sens. 2020, 12, x 18 of 26 WT01 target and WT02 target and between them. This implied that the proposed LUT method facilitated better relative calibration than the manufacturer's method.   It is evident from Figure 20 that the interval between the mean of the DN of WT01 and that of WT02 in the corrected results increased as compared to that in the original images. This may be attributed to the tilt observed by the sensor. However, the mean of DN of other similar targets became closer than before. This indicated the effectiveness of the two calibration methods. Moreover, the mean DN of all the same targets corrected by the LUT method was more uniform than that of the same targets corrected by the manufacturer's method, as shown in Figure 20. This implied that the effect of the LUT method was better than that of the manufacturer's method. The STD of the DN of the WT02 target in the corrected images based on the two methods was larger than that of the original images. This also indicated the problem of sensor tilt similar to Figure  19a and Figure 20. However, the STD of other targets became smaller than that in the original images. This confirmed that the uniformity of the image was improved. Besides, the STD of all the targets corrected by the LUT method was smaller than that of the targets corrected by the manufacturer's  It is evident from Figure 20 that the interval between the mean of the DN of WT01 and that of WT02 in the corrected results increased as compared to that in the original images. This may be attributed to the tilt observed by the sensor. However, the mean of DN of other similar targets became closer than before. This indicated the effectiveness of the two calibration methods. Moreover, the mean DN of all the same targets corrected by the LUT method was more uniform than that of the same targets corrected by the manufacturer's method, as shown in Figure 20. This implied that the effect of the LUT method was better than that of the manufacturer's method. The STD of the DN of the WT02 target in the corrected images based on the two methods was larger than that of the original images. This also indicated the problem of sensor tilt similar to Figure  19a and Figure 20. However, the STD of other targets became smaller than that in the original images. This confirmed that the uniformity of the image was improved. Besides, the STD of all the targets corrected by the LUT method was smaller than that of the targets corrected by the manufacturer's

Absolute Accuracy
In practical applications, reflectance products are universally used for remote sensing analysis. In this study, a simple empirical model was used to convert the radiance of absolute radiometric calibration to reflectance [18,34]. Then, the reflectance of targets was compared with the actual reflectance measured by ASD to assess the absolute accuracy of the absolute radiometric calibration. CRP with better isotropy was selected as the known object. The equation was as follows: where image j cal−norm was the k-band image, which had been corrected by absolute radiometric calibration; ρ k CRP was the reflectance of CRP in k-band; CRP k cal−norm was the average radiance of CRP in image j cal−norm ; image k re f was the reflectance of k-band image. The absolute error of reflectance ∆ was defined as the difference between the actual reflectance x i of the radiation target i and that measured by the sensor x i . RMSE was used to evaluate the overall measurement accuracy of the sensors.
The average reflectance of WT01, WT02, GT01, and GT02 was extracted from the reflectance images image k re f calibrated by the LUT method and the manufacturer's method. Then, the absolute error and RSME were calculated according to Equations (10) and (11). The reflectance error of targets in the five spectral images corrected by the LUT method and manufacturer's method is shown in Figures 22 and 23, respectively.
The absolute error of reflectance ∆ was defined as the difference between the actual reflectance , of the radiation target and that measured by the sensor . RMSE was used to evaluate the overall measurement accuracy of the sensors.
The average reflectance of WT01, WT02, GT01, and GT02 was extracted from the reflectance images calibrated by the LUT method and the manufacturer's method. Then, the absolute error and RSME were calculated according to Equations (10) and (11). The reflectance error of targets in the five spectral images corrected by the LUT method and manufacturer's method is shown in Figures 22 and 23, respectively.  Comparing the data in Figures 22 and 23, it was clear that the RMSE of reflectance of targets in the five spectral images corrected by the LUT method was less than that of spectral images corrected by the manufacturer's method. The RMSE of reflectance calibrated using the proposed LUT method were 2.30%, 2.87%, 3.66%, 3.98%, and 4.70% for Band1, Band2, Band3, Band4, and Band5, respectively.
It should be noted that the normalization procedure of the calibration method provided by the manufacturer of the multispectral camera deteriorated the result, and the calibration parameter of Comparing the data in Figures 22 and 23, it was clear that the RMSE of reflectance of targets in the five spectral images corrected by the LUT method was less than that of spectral images corrected by the manufacturer's method. The RMSE of reflectance calibrated using the proposed LUT method were 2.30%, 2.87%, 3.66%, 3.98%, and 4.70% for Band1, Band2, Band3, Band4, and Band5, respectively.
It should be noted that the normalization procedure of the calibration method provided by the manufacturer of the multispectral camera deteriorated the result, and the calibration parameter of integration time was not appropriate for the global shutter. Additionally, there was a systematic deviation between the radiance images of all the bands calculated by the two methods. This was probably because the performance of the sensor declined with use, and the absolute calibration coefficients measured before the delivery were unable to meet the current requirements. These two problems are shown in Figure 24.
Compared to Figure 19, there was a large deviation between the radiance value derived by two methods after normalization and absolute radiometric calibration. This is more obvious in the images of Band1 and Band2 in Figure 24. The curve of diagonal pixels is obviously raised on the left side and decreased on the right side in Figure 24b

Discussion
The radiometric calibration process of remote sensing instruments involves many aspects and is complex, especially for 2D image sensors. According to the comparison with the method provided

Discussion
The radiometric calibration process of remote sensing instruments involves many aspects and is complex, especially for 2D image sensors. According to the comparison with the method provided by the manufacturer of the sensor, it was verified that the proposed radiometric calibration procedures based on the LUTs achieve encouraging results for Rededge-MX multispectral camera.
The dark current was characterized in terms of integration time and gain. This has rarely been attempted in previous studies. It was proved that compared to the mean value of the dark current image, the LUT of dark current is more appropriate for correcting the non-uniformity of dark current in the 2D image space. The effect of temperature was not considered in the dark current analysis as the inclusion of temperature can complicate the problem, and this will be discussed in future work. Dark current experiments were conducted after the device temperature was stabilized to room temperature. During the process of obtaining the normalized LUT of the vignetting effect, extracting the vignetting effect from images is a challenge, which is also the basis for eliminating the vignetting effect. The proposed method based on the Gaussian filter is simple and efficient, and it successfully eliminated the vignetting effect. Notably, this method is more effective for irregular or complex vignetting effect as shown in Band4 and Band5 of RedEdge-MX. The reliability of the results was verified from two aspects: relative accuracy and absolute accuracy, and comparison with the calibration method provided by the manufacturer of the spectral camera. In earlier studies, only one of the verification approaches has been commonly used [14,24,42]. This limits the wide utilization of sensors. On the other hand, the verification results of the proposed radiometric calibration model provided excellent relative accuracy and absolute accuracy.
Overall, the proposed radiometric calibration method based on LUTs provides outstanding results for complex and low-cost sensors. Because the basic principles of imaging sensors are similar, this method can also be used for sensors of the same type or for simple, advanced sensors such as the linear sensor or 2D-frame CCD sensor. This validates the universality of the proposed method. However, it may be noted that this experiment was implemented using an integration sphere in the laboratory environment. For field applications, the integration sphere may be replaced in the reflectance calibration procedures, but the experimental conditions need to be strictly controlled, including the uniformity of reflectance calibration, direction of sensor observation, interference from surrounding objects, etc. Accurate radiometric calibration of sensors is the basis of remote sensing. Accordingly, our future research will focus on radiation correction and widening the applications of the sensor.

Conclusions
In this study, a radiometric calibration model was proposed to determine the LUTs of correction factors for the dark current, vignetting effect, and non-uniform detection response and the absolute calibration parameters were also analyzed. The LUT of dark current correction was determined as the image template captured at a camera setting of 1× gain and 1 ms integration time. The LUT of the vignetting effect correction was obtained as the image with normalized brightness filtered by the Gaussian function. The LUT of the response correction was the ratio of the image value of each pixel to the mean value of the image obtained after correcting the dark current and vignetting effect. Then, the parameters of absolute calibration were derived by the linear fitting.
The LUTs of correction factors effectively improved the images. The absolute calibration parameters precisely converted the DN value to the radiance unit. The calibration accuracy of each band was evaluated by observing the reference targets and comparison with the method provided by the manufacturer of the sensor. The results revealed that the performance of the proposed radiometric calibration method based on LUTs was much better than that of the method provided by the manufacturer in terms of both relative accuracy and absolute accuracy. The RMSE of reflectance of reference targets derived from radiance images were 2.30%, 2.87%, 3.66%, 3.98%, and 4.70% for Band1, Band2, Band3, Band4, and Band5, respectively. Overall, the proposed LUT method based on the inherent characteristics of the image effectively realizes radiometric calibration of small multispectral cameras. The present results can serve as a useful reference for the improvement of multispectral cameras and boosting their practical applications.