Multispectral Sensor Calibration and Characterization for sUAS Remote Sensing

This paper focuses on the calibration of multispectral sensors typically used for remote sensing. These systems are often provided with "factory" radiometric calibration and vignette correction parameters. These parameters, which are assumed to be accurate when the sensor is new, may change as the camera is utilized in real-world conditions. As a result, regular calibration and characterization of any sensor should be conducted. An end-user laboratory method for computing both the vignette correction and radiometric calibration function is discussed in this paper. As an exemplar, this method for radiance computation is compared to the method provided by MicaSense for their RedEdge series of sensors. The proposed method and the method provided by MicaSense for radiance computation are applied to a variety of images captured in the laboratory using a traceable source. In addition, a complete error propagation is conducted to quantify the error produced when images are converted from digital counts to radiance. The proposed methodology was shown to produce lower errors in radiance imagery. The average percent error in radiance was -10.98%, -0.43%, 3.59%, 32.81% and -17.08% using the MicaSense provided method and their "factory" parameters, while the proposed method produced errors of 3.44%, 2.93%, 2.93%, 3.70% and 0.72% for the blue, green, red, near infrared and red edge bands, respectively. To further quantify the error in terms commonly used in remote sensing applications, the error in radiance was propagated to a reflectance error and additionally used to compute errors in two widely used parameters for assessing vegetation health, NDVI and NDRE. For the NDVI example, the ground reference was computed to be 0.899 +/- 0.006, while the provided MicaSense method produced a value of 0.876 +/- 0.005 and the proposed method produced a value of 0.897 +/- 0.007. For NDRE, the ground reference was 0.455 +/- 0.028, MicaSense method produced 0.239 +/- 0.026 and the proposed method produced 0.435 +/- 0.038.