Photometric characterization of Lucideon and Avian Technologies color standards including application for calibration of the Mastcam-Z instrument on the Mars 2020 rover

Abstract. Several commercially available color standards exist, generated by a variety of manufacturers including LabSphere, Lucideon, and Avian Technologies. Previous work has characterized the photometric properties of LabSphere Spectralon targets. Here, we measure the visible and shortwave infrared (VSWIR; 0.4 to 2.5  μm) reflectance at multiple angles and determine the photometric properties of materials manufactured by Lucideon and Avian Technologies for potential use as calibration target materials for the Mars 2020 Mastcam-Z instrument. The Lucideon black, gray 33, green, and cyan samples are found to be significantly forward scattering. The yellow, red, and gray 70 samples are found to be weakly forward scattering. The Avian Technologies AluWhite98 sample was found to be weakly backward scattering. We characterize the absorptions observable and note the occurrence of wavelength-dependent photometric properties. The reflectance and photometric data collected and released here enable the use of these color standards for calibration of data from Mastcam-Z and other Mars-2020 rover instruments as well as provide key information for many other imaging and spectroscopy applications that require the calibration of data from multiple lighting or viewing geometries.


Introduction
Calibration targets are essential in many applications for accurate radiometric calibration of data. Here, we describe the spectral and photometric properties of color standards from Lucideon and Avian Technologies, which we characterized prior to use on the Mars-2020 (M2020) rover mission. M2020 includes the Mastcam-Z (MCZ) instrument, which consists of a pair of multispectral, visible/nearinfrared [VNIR, adjustable-zoom (focusable)] cameras. The science goals of the MCZ instrument, in brief, are to characterize the geomorphology and geology at the rover field site, assess current atmospheric and astronomical phenomena, and provide operational support and scientific context for the M2020 rover. 1 The MCZ instrument includes two fixed calibration targets (primary and secondary) each composed of a reference set of color and grayscale materials of known spectral and photometric properties 1 (Fig. 1). Illumination conditions vary on the surface of Mars according to the local solar time and the aerosol opacity of the atmosphere at the time of observation. A calibration target is necessary for verification of the preflight calibration of the camera, monitoring the stability of the instrument throughout the mission, and conversion of the radiance measured in data numbers to radiance factor or Lambert albedo to enable comparison with laboratory reflectance spectra of Earth materials (e.g., Ref. 2). In addition, the MCZ calibration targets may be useful for calibration of other onboard instruments such as the VNIR (400 to 900 nm) and the shortwave infrared (SWIR; 1300 to 2600 nm) passive spectra from SuperCam, as well as engineering cameras and cameras mounted on the rover arm. Furthermore, the targets serve as a platform of known spectral characteristics to monitor dust deposition and removal throughout the mission (e.g., Refs. [3][4][5]. For the calibration targets to enable these tasks, a detailed understanding of the photometric properties of the materials that compose the calibration targets is imperative. This is the goal of this study. This study replicates some of the techniques and analyses conducted for the common reflectance standard, Spectralon ® SRS-99 (Labsphere, North Sutton, New Hampshire). [6][7][8] The materials used here, manufactured by Lucideon (Staffordshire, United Kingdom) and Avian Technologies LLC (Sunapee, New Hampshire), are commercially available color standards. Therefore, we expect that the data acquired in this study will be of utility to others who employ these standards for imaging or spectroscopy applications or are evaluating them for potential use. We assess the wavelength-dependent spectroscopic properties and photometric behavior of the calibration materials (see data in Supporting Online Materials). This paper supplements detailed characterization of the actual assembled calibration targets and the derived bidirectional reflectance function (BRDF) for the calibration materials.
2 Mastcam-Z and its Calibration Targets MCZ's Bayer-pattern CMOS detector and 12 filters [ Fig. 1(a)] are used to measure the wavelength range of 445 to 1013 nm (preflight estimate) to discriminate ironbearing minerals and certain water-bearing minerals by detection of their characteristic VNIR absorptions and/or spectral shapes. 9,10 For each camera, an additional solar filter with reduced transmission enables observation of the sun to determine atmospheric opacity. The primary MCZ calibration target is an ∼8-cm-diameter disk with eight color/grayscale chips surrounding four concentric grayscale rings [ Fig. 1(b)]. A gnomon is in the center of this calibration target. The gnomon casts a shadow across the calibration target that allows comparison of direct versus diffuse components of the solar and sky irradiance (e.g., Refs. 2 and 11). The calibration target is similar to those flown on the Mars Exploration Rovers (MERs) Spirit and Opportunity and on the Mars Science Laboratory (MSL) Curiosity rover with modifications based on lessons learned 1 and based on requirements for use of materials caused by instrumentation on board. Because the payload holds very sensitive instruments for remote detection of organic materials, silicone rubber-based materials such as those used as color references in previously flown calibration targets were avoided. The primary MCZ calibration target is mounted horizontally and located on the rover at a position where the MCZ instrument has an unobstructed view, where it will be fully illuminated during surface operations, and where there is a minimum of diffuse or reflected light from other structures on the rover. The target is viewed from the rover mast at an angle of 58 deg from vertical (emission angle ¼ 58 deg), and geological materials on the surface are viewed at a wide range of angles. The secondary calibration target [ Fig. 1(c)] is mounted on a vertical surface on the rover for monitoring of calibration targets at different angles and with different accumulation of dust.
Calibration procedures for the Mars2020 rover will be similar to those implemented on the MER and MSL missions. For each multispectral observation of the surface, a corresponding multispectral observation of the primary calibration target under similar illumination conditions will also be acquired. 2,12 From these calibration target images, regions of interest (ROIs) of uniform dust cover will be selected from each of the four grayscale rings and from the eight relatively dust-free central spots of the eight circles around the periphery. The ROI mean radiance values are plotted against laboratory-measured directional-hemispherical reflectance values of the calibration targets corrected for illumination and viewing geometry using their photometric properties. 5 To correct for the presence of dust, a two-layer scattering model is used, closely following that of Ref. 5 with a key difference in that all irradiance is assumed to come from the direction of the sun. 2 The result of this model is a conversion factor, which can be used to correct observed radiances from the Martian scene to radiance factor (I/F). 5 Lastly, the calibrated radiance factor can be divided by the cosine of the solar incidence angle to determine the relative reflectance, also known as RÃ (REFF; see Sec. 3). [12][13][14] Given the range of possible illumination geometries, full knowledge of the photometric properties of calibration target materials is required to perform this analysis.

Methods
To characterize the photometric properties of the candidate MCZ calibration target materials, it was necessary to obtain reflectance measurements at a series of geometries with a stable and well-calibrated detector and light source. We used a multiwavelength spectrogoniometer developed at Caltech [ Fig. 2(a)] that is capable of measuring at all azimuth (Az) angles (0 deg to 180 deg) and at emission (e) and incidence (i) angles from −70 deg to 70 deg with a minimum phase angle of ∼12 deg. Incidence and emission angles are measured as positive angles from the vertical, and azimuth is measured positive clockwise from the horizontal projection of the incidence angle [ Fig. 2(b)]. The phase angle (g), the angle between the incidence and emission angle, is calculated as follows: E Q -T A R G E T ; t e m p : i n t r a l i n k -; s e c 3 ; 6 3 ; 5 8 4 g ¼ arccos½cosðiÞ Ã cosðeÞ þ sinðiÞ Ã sinðeÞ Ã cosðAzÞ: The measurements were taken with an Analytical Spectral Devices FieldSpec 3 ® , which allows for reflectance measurements from 350 to 2500 nm through the use of a 512-element silicon array detector in the VNIR and two graded InGaAs photodiodes in the SWIR. These measurements are comparable to those taken with the Reflectance Experiment Laboratory (RELAB) instrumentation at Brown University. 15,16 The sample is illuminated with a halogen source with all other room lights turned off. Baffling ensures that the signal comes from a <7 mrad cone, corresponding to diameters on sample of <5 mm at e ¼ 0 deg and <2 cm at e ¼ 70 deg. Formally, this is a conical-conical reflectance (case 5, in Ref. 17) with small angles nearing the ideal BRDF. For a white reference standard, we used a Spectralon ® target that is approximately Lambertian. [7][8] The Caltech goniometer was specifically designed for the Spectralon calibration target to lie coincident with the axes of rotation of the goniometer. To match this measurement plane, the MCZ calibration target reference materials were lifted using thin washers.
Prior to measurements, we optimized the gain on each detector with the Spectralon ® present using a geometry of i ¼ 12 deg and e ¼ 0 deg. We then positioned the goniometer at the specified geometry and obtained a spectral average of 400 white reference measurements. (This includes a dark current measurement also using a 400 spectral average.) The reflectance of each candidate target material was then measured at this geometry, also using a 400 spectral average. Every 10 min the white reference was remeasured, and its spectrum saved.
To account for the differences in each detector, we used a routine available with the ASD FieldSpec 3 ® software. This routine uses five channels from each detector at the filter boundary to calculate a bias value for the VNIR and longer wavelength SWIR detectors. This bias value matches the offset with the shorter wavelength SWIR detector at the splice point.
Following Refs. 18 and 19, we output the reflectance as the reflectance factor (REFF) that is defined as follows: E Q -T A R G E T ; t e m p : i n t r a l i n k -; s e c 3 ; 3 2 6 ; 4 9 9 REFFði; e; g; λÞ ¼ L samp ði; e; g; λÞ − L dark ðλÞ L id ði; e; g; λÞ − L dark ðλÞ ; where L samp ði; e; g; λÞ is the radiance from the target measured at geometry ði; e; gÞ at a given wavelength, λ, L id ði; e; g; λÞ is the radiance from an ideal, perfectly diffuse standard surface measured at geometry ði; e; gÞ, and L dark ðλÞ is the radiance measured with no incident light. where REFF meas is the measured wavelength-dependent reflectance relative to the L std , the Spectralon standard measured in the specified geometry. And R std is the National Institute of Standards and Technology-certified value for the wavelength-dependent reflectance of Spectralon ® .
Because the purpose of this study was to obtain an understanding of the photometric properties of the target materials, we obtained the REFF measurements at 76 geometries, illustrated in Fig. 3. To evaluate the scattering properties of the samples, we observed the change in REFF in-and out-of-plane at multiple geometries with e ¼ 58 deg and e ¼ 30 deg as well as for in-plane geometries at e ¼ 0 deg and i ¼ 30 deg (Fig. 3). The emission angle of 58 deg was chosen because it is the approximate emission angle at which the MCZ detector will view the calibration target panel; e ¼ 30 deg is a standard photometric study observation point and is similar to the path from the instrument to the secondary target (e ≅ 32 deg). The remaining angles were chosen to provide a full-suite of measurement geometries enabling a characterization of the photometric properties. We assume symmetry in the photometric properties for azimuth angles 180 deg to 360 deg. Each measurement in this study includes the full spectral wavelength range of 350 to 2500 nm. The data were used to determine the wavelength-dependent reflectance of the calibration targets and their photometric properties, including characterization of wavelength-dependent changes in photometric behavior, which has been observed in prior studies of geological materials (e.g., Ref. 16). Modeling of Hapke parameters following Ref. 3 was attempted but produced poorly constrained fits and is not further discussed. Follow-up studies will obtain measurements of the targets at finer angular resolution but at limited wavelengths to determine the functional form of the BRDF.
To ensure that the sample holder material did not influence the spectra collected, we measured the sample holder material directly and noted its spectral properties (which were not in common with any of those of the target material and are not further discussed). Lastly, we confirmed absorptions and peaks in the 400-to 1100-nm range by measuring the samples using an Ocean Optics HR200CG ultraviolet-NIR Spectrometer at i and e of 0 deg.

Samples
Eight samples are candidate calibration target reference materials based on their ability to withstand the rover environment and their photometric properties. These samples are shown in Fig. 4 with the Spectralon ® white reference also shown. The red, yellow, green, cyan, gray 33, gray 70, and black samples are manufactured by Lucideon while the AluWhite 98 (from now on referred to as white) target is manufactured by Avian Technologies LLC. The Lucideon targets are 5 cm × 5 cm glazed and matted aluminum silicate tiles. The white target is a 5.08 × 5.08 cm matte tile composed of sintered alumina.

Data Quality Evaluation and Known Instrument Artifacts
The Caltech goniometer data were examined in Ref. 16, compared to that acquired in RELAB, and matched within 10%. Comparison of our measurements on 5%, 20%, and 99% reflectance standards with Spectralon ® calibration files provided by the manufacturer indicated correspondence within 5%, indicating no significant contribution from the light source, linear instrument response, and adequate instrument calibration. During our measurements, white reference spectra were acquired every 10 min to verify that the apparatus was stable. Data for wavelengths short of 400 nm and longer than 2500 nm were evaluated for quality and excluded due to low signal-to-noise ratio. The standard deviation REFF between 100 successive measurements (over ∼5 min) of Spectralon ® is ≤0.012, indicating stability in the light source and instrument (Fig. 5). The 400-to 2500-nm wavelength range is reliable with the exceptions below where caution in interpreting data is needed.
In the raw data from the spectrometer and also occasionally in the corrected reflectance data, we observe spectral features at ∼1100 and ∼1320 nm (dashed gray lines in Fig. 5). According to correspondence with the makers of the ASD instrument, these absorptions are related to varying proportions of on-axis and off-axis signal entering the spectrometer, depending on brightness. These artifacts affect several of our spectra and are more pronounced at certain geometries. Consequently, the wavelength ranges 1100 AE 30 nm and 1320 AE 30 nm should be ignored in our REFF data. We also observe weak absorptions at ∼1400 nm and ∼1900 nm (dashed blue lines in Fig. 5), which can be caused by vibrational absorptions related to water in the instrument path length. In most REFF data, this wavelength region is reliable. There are strong 1400 and 1900 nm absorptions in the samples themselves due to water in the materials; however, smaller, superimposed fluctuations should be interpreted with caution.

Target Spectral and Photometric Properties
A summary of the spectral and photometric properties of the target materials is given in Table 1. The most illustrative BRDF measurement results are given in Figs. 6 and 7; all data and normalization plots are available in the supplementary materials. All targets demonstrate some absorptions in the 400-2400 nm range (Fig. 6). All absorptions shortward of 1000 nm are broad absorptions that control the target color properties. At longer wavelengths, all samples except the gray and black samples exhibit absorptions at 1900 nm due to water. The white, yellow, and red targets exhibit 1400 nm features due to water, and the yellow and red even have small overtones at 940 and 1120 nm related to molecular water. The white and four color targets also exhibit an absorption at 2200 nm due to Si-OH or Al-OH. Additionally, the white sample has a 694-nm peak [ Fig. 6(h)] that was confirmed with the UV-NIR spectrometer and corresponds to a welldocumented fluorescence in chromium, 20 a common substitution in Al 2 O 3 , of which the white target is composed. It may be related to a Cr for Al substitution in the structure.
The black, gray 33, green, and cyan samples are significantly forward scattering at 800 nm, i.e., >200% change in reflectance over the angles measured. The yellow, red, and gray 70 samples are quasi-isotropic but weakly forward scattering at 800 nm. The white sample is quasi-isotropic but weakly back scattering at 800 nm (Fig. 6). All the samples exhibited wavelength-dependent photometric properties, some of which were quite significant ( Fig. 7; Table 1). The red, yellow, green, blue, and black samples were more forward scattering at short wavelengths and become significantly less forward scattering at longer wavelengths.
The gray 33 and gray 70 samples were found to be less forward scattering at shorter wavelengths. For all samples except white, measurements made at phase angles ∼103 deg and ∼112 deg [arrows in Fig. 7(a)] deviate from the general upward trend of greater normalized REFF with larger phase angle (these measurements were verified by repeat measurements).

Convolution to MCZ Filter Set
We approximated the MCZ camera multispectral data of the calibration target reference materials (Fig. 8) by integrating the full-resolution spectra over MCZ filter transmission profiles [ Fig. 1(a)]. At the time of this work, the MCZ filters had not been manufactured and verified, and thus the MCZ profiles were estimated as the replicas of the Mastcam filter transmission profiles. 21 The profile of the Mastcam 805 nm filter was duplicated and shifted to create hypothetical bandpasses for the ∼600 and ∼975 nm filters that are new to MCZ. Overall, the MCZ filter set reproduces the general spectral shape of the targets.

Slope as a Function of Viewing Geometry
There is a significant change in spectral slope with viewing geometry for certain measurement suites (Fig. 9). The degree of this change is variable across samples and measurement suites. As an example, Fig. 9 demonstrates two measurement suites for the gray 33 and yellow samples. The yellow sample appears to have little to no change in spectral slope for either measurement suite. The gray 33 sample, however, has significant change in spectral slope at the Az ¼ 180, e ¼ 58 measurement suite, i.e., at higher phase angle. This effect is still observed, though to a lesser degree at the e ¼ 30   Optical Engineering 027108-9 February 2019 • Vol. 58 (2) suite. This is likely due to wavelength-dependent photometric properties indicated by the data in Fig. 7. We examined the whole dataset, and the green, black, and gray 33 samples exhibit significant changes in slope with viewing geometry, whereas the rest exhibit this effect to a lesser extent. There is no consistent change in slope with viewing geometry observed in the white sample. The reader may visualize these with data available in the supplementary materials.

Change in Absorption Strength as a Function of Viewing Geometry
Two absorptions were analyzed for changes in band depth with viewing geometry: the 1470-nm absorption band on the green sample and the 660-nm absorption band on the gray 33 sample (Fig. 10). As an example, on the green sample, the band depth decreases with increasing phase angle. The sample, however, is strongly forward scattering as previously mentioned, so there is also an increase in REFF with phase angle. There is no observable trend with absorption strength for the gray 33 sample. All samples with wavelength-dependent photometric properties (Fig. 10) can exhibit these band depth changes, which can be examined in detail with data in the supplementary material.

Discussion
For some applications where photometry is important, such as in observations of materials on Mars with the Mars2020 rover, this study shows that the characterization of the angular dependence of sample reflectance for these candidate calibration target materials is critical. The suite of targets do not all behave in the same manner with weakly backward scattering (white), weakly forward scattering (yellow, red, and gray 70), and strongly forward scattering (cyan, green, gray 33, and black) behavior, depending on the sample and the wavelength of measurement. Just as with the RTV materials (pigmented synthetic silicone rubber) used in calibration targets flown on earlier missions to Mars, these calibration target materials are not Lambertian. Several of these vary significantly from isotropic behavior, but their photometric properties vary in consistent ways that can be described. Additionally, none of the targets show significant specular behavior. When compared with the calibration targets found on the MER rover calibration target and MSL Curiosity's Mastcam calibration targets, 2,11 these candidate materials from commercial vendors are spectrally flatter in the 400-to 700nm range, with the exception of the green target. Sharp absorptions should be avoided in calibration target materials, especially if they occur near the boundaries of the bandpass filters, but these targets have smoothly varying properties in the relevant 400-to 1100-nm range for Mars2020/MCZ and for any other VNIR imaging and spectroscopy investigations. The change in photometric properties with wavelength for each target provides an explanation for why changes in slope and absorption strength are observed in the samples with changes in viewing geometry. For analysis of spectra, this highlights why correction for calibration target photometric properties are crucial as the sun angle on the calibration target varies.
The white emission peak at 698 nm is an interesting, unexpected exception, but as it occurs within the center of the bandpass and is relatively small, it should not affect  calibration. Sharp absorptions are present at longer wavelengths in the red, yellow, green, cyan, and white calibration materials tested, so users for SWIR multispectral imaging applications should consider these when selecting calibration materials and bandpass cutoffs. The emission peak in white and the sharp absorptions at larger wavelengths may provide additional information to characterize wavelength position. On the other hand, because several sharp absorptions are due to bound water in the sample, it is possible that they will weaken or shift a few 10's of nm due to changing absorption shape with reduced water content under low Martian pressure or vary as a function of time of day (due to changes in relative humidity, Refs. 22 and 23). For observations at longer wavelengths than the MCZ instrument, e.g., by SuperCam on Mars, the narrow absorptions confound use as uniform surfaces but may prove useful for monitoring dust buildup (by absorption shallowing), complementing short wavelength observations. The downward trending REFF measurements at phase angles ∼103 and ∼112 may be due to internal instrument effects (i.e., related to on and off axis light) or may be manifestations of actual scattering properties of the targets. With this data set, it is not possible to distinguish between these scenarios.

Conclusions and Future Work
The measurement suite captures the effect of viewing geometry on REFF for one Avian Technology target (AluWhite98) and seven Lucideon targets (red, yellow, cyan, green, gray 70, gray 33, and black), characterizing their photometric properties. Some of the samples are quasi-isotropic: the white target  is slightly backward scattering and the red, yellow, and gray 70 targets are slightly forward scattering. The cyan, green, gray, and black targets are strongly forward scattering. The forward scattering behavior is most significant at shorter wavelengths (i.e., in the MCZ measurement range) and for all samples except the gray samples, decreases at longer wavelengths. This wavelength dependence causes changes in spectral slope and absorption band depth depending on the viewing geometry. The scattering properties have been quantified here and can thus be corrected for. Our data alone are too sparse for a full BRDF characterization as was done in Ref. 11, but this will be future work at a more restricted set of wavelengths. Additional measurements will be taken at the Bloomsberg University Goniometer, and the data will be integrated into a scattering model covering all geometries.  Optical Engineering 027108-12 February 2019 • Vol. 58 (2) assistance in development and implementation of attempted photometric models. MBM gratefully acknowledges supporting grant CF16-0981, which covered salary for KK but also purchase of the calibration target references used in this study. This work was funded by a Mars-2020, Mastcam-Z grant to BLE, and a NASA Earth and Space Sciences Fellowship #NNX15AQ95H partially supported JB.