**The Pan-and-Tilt Hyperspectral Radiometer System (PANTHYR) for Autonomous Satellite Validation Measurements—Prototype Design and Testing**

**Dieter Vansteenwegen 1,\*, Kevin Ruddick 2, André Cattrijsse 1, Quinten Vanhellemont <sup>2</sup> and Matthew Beck <sup>2</sup>**


Received: 17 May 2019; Accepted: 3 June 2019; Published: 6 June 2019

**Abstract:** This paper describes a system, named "pan-and-tilt hyperspectral radiometer system" (PANTHYR) that is designed for autonomous measurement of hyperspectral water reflectance. The system is suitable for deployment in diverse locations (including offshore platforms) for the validation of water reflectance derived from any satellite mission with visible and/or near-infrared spectral bands (400–900 nm). Key user requirements include reliable autonomous operation at remote sites without grid power or cabled internet and only limited maintenance (1–2 times per year), flexible zenith and azimuth pointing, modularity to adapt to future evolution of components and different sites (power, data transmission, and mounting possibilities), and moderate hardware acquisition cost. PANTHYR consists of two commercial off-the-shelf (COTS) hyperspectral radiometers, mounted on a COTS pan-and-tilt pointing system, controlled by a single-board-computer and associated custom-designed electronics which provide power, pointing instructions, and data archiving and transmission. The variable zenith pointing improves protection of sensors which are parked downward when not measuring, and it allows for use of a single radiance sensor for both sky and water viewing. The latter gives cost reduction for radiometer purchase, as well as reduction of uncertainties associated with radiometer spectral and radiometric differences for comparable two-radiance-sensor systems. The system is designed so that hardware and software upgrades or changes are easy to implement. In this paper, the system design requirements and choices are described, including details of the electronics, hardware, and software. A prototype test on the Acqua Alta Oceanographic Tower (near Venice, Italy) is described, including comparison of the PANTHYR system data with two other established systems: the multispectral autonomous AERONET-OC data and a manually deployed three-sensor hyperspectral system. The test established that high-quality hyperspectral data for water reflectance can be acquired autonomously with this system. Lessons learned from the prototype testing are described, and the future perspectives for the hardware and software development are outlined.

**Keywords:** Hyperspectral reflectance; validation; autonomous measurements; ground-truth data; system design

#### **1. Introduction**

The objective of this paper is to describe the motivation, design, and prototype testing for an autonomous system of hyperspectral radiometers suitable for validation of satellite-derived water reflectance.

#### *1.1. Motivation and Objective*

Satellite imagery of marine, coastal, and inland water reflectance is now routinely used for measuring parameters such as the concentrations of chlorophyll a, a proxy for phytoplankton biomass, and suspended particulate matter, important for sediment transport applications. The satellite data are used for regulatory monitoring of the aquatic environment, e.g., via the European Union Water Framework and Marine Strategy Framework Directives [1,2], and for providing a scientific basis for coastal zone decision-making, e.g., via the assessment of impacts of human activities and constructions [3]. However, end-users of the data require reliable information on data quality, and validation of the satellite data at the level of water reflectance is particularly crucial. This is because of the large errors that may occur during the data calibration and processing, particularly during the atmospheric correction steps [4].

This validation is best achieved by a "match-up" of in situ measurements of water (surface) reflectance made at the same time as the satellite measurement [5], and experience over the last 10 years showed that only autonomous in situ systems can provide sufficient data for this purpose. In particular, AERONET-OC [6], a federated network of multispectral robotically pointed radiometers on offshore platforms all over the world, proved to be the main source of validation data [7] for spaceborne optical missions such as ENVISAT/MERIS, MODIS/AQUA, VIIRS, Sentinel-3/OLCI, Sentinel-2/MSI, Landsat-8/OLI, etc. However, the radiometer adopted within AERONET-OC is only multispectral and cannot adequately cover the spectral bands of all recent and future optical spaceborne missions without spectral interpolation/extrapolation/modelling [8] and associated uncertainties. The WATERHYPERNET network is, therefore, being developed, based closely on the concept of the successful AERONET-OC network [6] but with the essential advantage of a hyperspectral radiometer, thus enabling the validation of all visible and near-infrared bands of all present and future satellite missions providing water reflectance data.

The objective of the present paper is to describe the measurement system, called PANTHYR (pan-and-tilt hyperspectral radiometer) that was developed for use within the WATERHYPERNET network. This measurement system consists of two commercial off-the-shelf (COTS) hyperspectral radiometers, mounted on a COTS pan-and-tilt (PT) pointing system, controlled by a single-board-computer and associated custom-designed electronics which provide power, pointing instructions, and data archiving and transmission.

#### *1.2. Measurement System Requirements*

The PANTHYR system was developed to fit the following user requirements:


• Pointing accuracy of at least 5◦ azimuth and 1◦ zenith.

#### *1.3. Precursor Autonomous Systems*

A few autonomous systems with pointable radiometers already exist and are briefly described here. The most used autonomous system for measurement of water reflectance is the Seaprism version of the CIMEL CE318 multiband photometer system. This system was originally commercialized in the 1990s for the measurement of aerosol properties and includes direct sun, near sun, and principal plane and almucantar sky radiance measurements, and it became the unique instrument of the AERONET network [9]. The Seaprism version was developed in the early 2000s [10] by reprogramming of the pointing system, especially to perform downward pointing measurements of upwelling radiance from water. The system evolved over the years to include improvements of the optical and electronic components, as well as development of new versions, e.g., including the possibility of nighttime lunar measurements [11]. The system consists essentially of an optical head with two fore-optics protected by collimators, containing a filter wheel for multiband optical measurements with a very wide dynamic range (from direct sun to dark water), a robotic pointing system, a control box providing pointing and measurement instructions and managing power and data, and associated auxiliary equipment for power generation (solar panels) and data transmission (METEOSAT satellite network uplink or cellular link). The system is extremely robust, giving reliable maintenance-free performance for long-term (~1 year) deployments in a very wide range of environments including land and water sites from the tropics to the polar regions.

The OSPREY system [12] was designed to be commercialized by Biospherical Instruments as a very-high-performance modular system of radiometers on pointing systems. The radiometers include a hyperspectral spectrometer and numerous single-band microradiometers within a single thermally controlled casing and a filter wheel in front of the spectrometer allowing polarimetric, direct sun, straylight-corrected, and dark measurements. The system, described in detail in Reference [12], has a very high dynamic range and high accuracy pointing, and is capable of water, sun, sky, and moon radiance measurements.

The RFLEX system [13] consists of three hyperspectral TRIOS/RAMSES radiometers measuring downwelling irradiance, downwelling sky radiance, and upwelling water radiance mounted at fixed zenith angles (0◦, 40◦, and 140◦, respectively) on an azimuthally rotating platform with associated control software. This system was designed principally for deployment on moving ships, where the azimuthally rotating platform allows achieving an optimal 90◦ or 135◦ relative azimuth to sun for any ship heading. Instructions for system hardware construction using low-cost components (except for the COTS radiometers) were made publicly available via SourceForge (https://sourceforge.net/projects/rflex/) and source code for communication with TRIOS/RAMSES was made available via GitHub (https: //github.com/StefanSimis/PyTRIOS). The system was used operationally on two "ferryboxes" mounted on ships of opportunity in the Baltic Sea.

The DALEC system [14] is commercially available from Insitu Marine Optics and consists of three hyperspectral spectrometers measuring downwelling irradiance, and downwelling (sky) and upwelling (water) radiance embedded at fixed zenith angles (180◦, 140◦, and 40◦, respectively) within a compact azimuthally rotating body on a gimballed mount. Control and data are managed by a standard personal computer (PC) using software supplied by the manufacturer. This system was designed principally for deployment on moving ships but can also be used from fixed platforms. The system is deployed on research vessels and an example of data usage is described by Reference [15].

The SAS Solar Tracker system from Seabird Electronics consists of three hyperspectral Seabird/HyperSAS radiometers measuring downwelling irradiance, downwelling sky radiance, and upwelling water radiance mounted at fixed zenith angles (180◦, 140◦, and 40◦, respectively) on an azimuthally rotating platform, with compass-based sun-tracking to ensure optimal relative azimuth angle with respect to sun. Power and data are managed via a deck interface unit and a standard PC. The system is used routinely on a few research vessels and ships of opportunity [16].

The WISP-3 is a handheld system, commercially available from Water Insight, with three embedded radiometers measuring water reflectance from downwelling irradiance, sky radiance, and water radiance. Data for monitoring parameters such as chlorophyll a concentration are logged internally and can be downloaded to a PC as processed data. A variant called WISPstation is under development for autonomous operation. An example of data usage and inter-comparison with other systems is described in Reference [17].

The current PANTHYR design will provide hyperspectral reflectance data not available from the Seaprism system and has advantages over the RFLEX, DALEC, Suntracker, and WISP designs because the zenith pointing flexibility allows sky and water radiance measurements to be made by the same radiometer and allows downward pointing to protect instruments when not measuring. Use of a single radiance sensor provides cost saving and reduces the uncertainties associated with spectral and radiometric calibration differences of the two-radiance-sensor systems. The disadvantage of the PANTHYR sequential measurement of water and sky radiance as compared to the simultaneous measurements possible with the two-radiance-sensor systems is considered to be minimal for the good clear sun and sky illumination conditions needed for satellite validation. If illumination conditions are suboptimal, this will be detected by the PANTHYR replicate measurements and will lead to data rejection. The PANTHYR system does not pretend to achieve the very high performance expected of the OSPREY system, but represents a much lower-cost alternative which should be ideal for deployment at multiple sites worldwide by organizations with moderate budgets.

All of these systems, including PANTHYR, are developing within a context where hyperspectral spectrometers are rapidly evolving with reduction of size, power, and cost (thanks to mass production for medical and industrial applications), and where pointing systems are becoming more affordable and more easily available as COTS items (thanks to mass production for applications such as video surveillance). Other technologies facilitating PANTHYR development include the massive improvements and cost reductions in data transmission (thanks to mobile phone networks) and in microcomputers.

#### *1.4. Overview of Paper*

In this paper, the design choices are described in Section 2, together with details of the control electronics, where most of the original developments were made, and the data acquisition protocol and data processing steps, which are strongly based on precursor systems. A first seaborne test of the system on Acqua Alta Oceanographic Tower (AAOT) is described in Section 3, where data acquisition and processing are also described for two precursor systems used for comparison, a manually supervised three-sensor hyperspectral system and the autonomous multispectral AERONET-OC system. Results from the prototype tests are described in Section 4. General system performance is evaluated with description of lessons learned for system improvements. A comparison is made between data acquired with PANTHYR and data from the two precursor systems. Finally, conclusions from these prototype tests and future perspectives for system refinement are outlined in Section 5.

#### **2. System Design**

#### *2.1. Top-Level Design Choices*

To stay within the design and future purchasing budget, we decided to use mostly COTS components. Self-design of all hardware would allow for greater optimization, but design and manufacturing of all components would add cost and complexity beyond the scope of our goals. Only when commercial hardware was not available or sufficient did we create our own. As availability of the chosen components will change in the next few years, special care was taken to implement a modular design in hardware and software. Individual component changes will require updated software/hardware for only that part, while most of the system remains the same. This also facilitates adding extra hardware and capabilities in the future.

The TRIOS/RAMSES COTS hyperspectral radiometer was chosen because of its relatively high performance, highly robust and mature low-power design, and moderate price. This radiometer provides measurements for a spectral range of 320–950 nm and was used for satellite validation since 2002 [18] by many teams worldwide. The performance of this class of radiometer was extensively studied and characterized with regard to straylight [19], thermal sensitivity [20], polarization sensitivity [21], angular response of the irradiance units [22], and non-linearity. A low-cost, portable, LED-based FieldCAL device is available for rapid checking of radiometric sensitivity in the field.

The FLIR PTU-D48E COTS pointing system was chosen for the PANTHYR system because of its high performance (including good pan and tilt speeds and pointing accuracy), robust design, and moderate price.

While COTS rugged PCs are widely available, a key early decision was to use a small low-power embedded computer board with self-design of additional electronics for managing power and connectivity to components such as the radiometers, GNSS (Global Navigation Satellite System) receiver, and auxiliary sensors. This required more development time, but provided more control over power management and greater flexibility in connecting and supporting additional external devices.

The Linux operating system was chosen to facilitate software portability to next-generation hardware platforms and to improve reliability and security. Linux is open source, easy to customize, and ideal for an offshore system where logistics make on-site intervention difficult and/or expensive.

Python was chosen as the main programming language as it is available under an open-source license, is well known in the scientific community, and produces clear, readable code.

#### *2.2. System Overview and Key Components*

The overall PANTHYR system is outlined in Figure 1, showing the main hardware elements and the associated power and data connections. A small embedded computer board controls all components and forms the heart of the system. After start-up, a GNSS receiver provides UTC (Coordinated Universal Time), as well as location, allowing the system to calculate the position of the sun. An optional IP network camera can take still pictures from the measurement target areas, allowing users to check any suspect data for unusual conditions such as floating vegetation or debris, boats, birds, and other obstructions in the field of view. During a measurement cycle, the controller calculates the head position for each measurement step, points the instrument in the required zenith and azimuth angles, and makes a measurement. At the end of a measurement cycle, the head moves to a "park" position where the instruments are pointed downward to prevent fouling. The system then goes into a sleep state, conserving power while waiting for the next cycle.

A myriad of single-board computers with ever increasing performance became available on the market in the last few years. These small (around credit-card size) and cheap (<100 €) boards manage to run an operating system and have on-board storage, as well as a network connection, while consuming as little as 0.5 W of power in standby. For our application, the BeagleBone Black has the following advantages:


To connect and interface the controller to the rest of the system, we designed two electronic boards which are plugged on top of the BeagleBone. The first board translates the voltage levels of the serial ports to RS-232-compatible levels so that the controller can communicate with the instruments and GNSS receiver. The second board provides protection and filtering on the power lines, as well as allowing power saving by cutting the power to external devices when they are not in use. Five protected power outputs with a maximum power draw of 2 A each can be controlled from software. An external solid-state relay is used to control power to the head, which is the only component that needs a switched 24 V direct current supply.

**Figure 1.** Pan-and-tilt hyperspectral radiometer system (PANTHYR) system diagram.

A number of commercially available PT heads were evaluated. We looked at price, resolution, power consumption, and ruggedness. The FLIR PTU-D48 E matched our requirements most closely. This unit was developed as a rugged solution for demanding applications and has the necessary specifications and certifications such as temperature range, ingress protection 67, and salt spray protection (MIL-810G). Its maximum payload of 7.5 kg leaves margin for extra hardware.

To allow the full pan (+/−174◦) and tilt (+90◦/−30◦) movement, a free loop of cable is required between the rotating instruments and the base of the head. We found no good solution to prevent this loop of cables from getting in front of the instruments at certain positions. The next version of PANTHYR will use the same head but with built-in slip rings, providing internal connections in the head assembly to replace outside cables. This option also adds 360◦ azimuth movement capability to the head. The limited tilt range of the PTU-D48 E means that the irradiance sensor which needs to be vertical during measurements cannot be parked lower than 30◦ below the horizon; it is thought that this will provide adequate protection from the elements, but longer-term testing is required to confirm this. The limited tilt range (120◦) of the PTU-D48 E also precludes an arrangement with parallel radiance and irradiance sensors since the radiance sensor would need to be tilted through 140◦ to go from vertical (zenith 180◦) to the necessary water-viewing angle (zenith 40◦). The radiance sensor was, therefore, fixed at an angle of 40◦ to the irradiance sensor, giving a zenith angle range for the irradiance sensor of 180◦ (downwelling irradiance measurement) to 60◦ (parked) and a zenith angle range for the radiance sensor of 140◦ (sky radiance measurement) through 40◦ (water radiance measurement) to 20◦ (parked).

On average, the PT unit consumes 12.3 W during start-up, 6.2 W during hold, and up to 20 W during combined axial movements. This demonstrates the necessity of keeping the system in a low-power sleep mode as much as possible.

The InSYS MRO-L200 gateway/switch connects all Ethernet devices, as well as providing the gateway to a 4G cellular network. This industrial device can set up and manage a VPN connection to the onshore server, and has a serial port and digital output, which could be used to remotely power cycle PANTHYR in case of problems. Its high power consumption (2.5 W typical) is the main drawback. A low-power network switch and a separate 4G gateway that is switched on only during short intervals would mitigate this, but also remove our option to connect to the installation in case of problems.

Two mechanical structures were designed to mount to existing structures at the measurement sites and install the instruments. A versatile mounting system accommodates different measurement sites, while remaining simple. We started with a 50-mm steel pipe as a range of clamps are commercially available for this standard diameter. Welded on top of this pipe is a triangular plate. The head is

mounted on a second similar plate. Three threaded rods connect both plates atop each other. Leveling of the head is achieved by changing the spacing along the three rods. The result is a simple system with a variety of mounting options (Figure 2a).

On top of the head, both sides of a U-shaped bracket accept adapter plates that hold the radiometers. Bolting holes for fixed inclinations of 20◦ and 60◦, as well as a slot that allows a variable 10◦ to 70◦ inclination, are available (Figure 2b).

**Figure 2.** Engineering drawings of mounting structures. (**a**) Support providing a flexible connection between a platform structure, e.g., horizontal railings or a vertical post fixed through the swivel couplings, and pan/tilt unit to be mounted on uppermost horizontal surface; (**b**) sensor mounting bracket for attaching radiometers to top of pan/tilt unit.

#### *2.3. Software and Usage*

The main Python script checks system and user settings and controls all actions. Data and settings are stored in an SQLite database. Software will be released under the GNU GPLv3 (GNU's Not Unix General Public License version 3) at the GitHub page: www.github.com/hypermaq/panthyr.

PANTHYR has one SQLite database which contains the measurement protocol, system settings, logs, and measurement results, as well as a task queue. The user can use a (remotely) connected laptop to access a webpage where the contents of the database can be viewed, changed, and exported (Figure 3). This serves as the main user interface for system configuration.

A "worker script" regularly checks the "queue" database table for tasks that need to be performed. These tasks range from executing a measurement cycle to setting up station parameters.

To perform a measurement cycle, the system gets the necessary settings from the "settings" table and read the first line in the "protocol" table. Each line in this table contains the parameters (instrument, zenith angle, azimuth offset, number of repetitions) describing one sub-cycle in the measurement cycle. Measurement results are written to the "measurements" table. When all scans in the protocol are finished, the system resumes a standby state where it regularly checks for new tasks. In the event of a failure, a log entry is created in the "logs" table. After three failed attempts, a task is ignored to prevent it from blocking the system.

Apart from the database frontend, PANTHYR also hosts an FTP (File Transfer Protocol) server to download log files and pictures, as well as an SSH (Secure SHell) server to allow low-level system maintenance. To prevent provider restrictions enforced on cellular networks, PANTHYR initiates an outgoing OpenVPN connection to an onshore server. Users can connect to the same server and communicate with PANTHYR as if it were next to them.


**Figure 3.** Screenshot of SQLite database interface.

#### *2.4. PANTHYR Data Acquisition Protocol*

The measurement protocol presented here is based on various precursors including References [6,25,26] and, as usual for above-water radiometry, gives a water-leaving radiance or reflectance for an off-nadir geometry, which can be matched to satellite viewing geometries by appropriate models of a bidirectional reflectance distribution function. Research to determine an optimal measurement protocol for PANTHYR operations within the WATERHYPERNET network is in progress and some aspects may change in the future, e.g., potential additional zenith and/or azimuth angles and/or a different number of replicate measurements, especially for water viewing. However, the current protocol is considered as already sufficient for demonstrating performance of the hardware.

The PANTHYR system performs automated measurements every 20 min from sunrise until sunset. Each cycle consists of measurements with a 90◦, 135◦, 225◦, and/or 270◦ relative azimuth to sun. In general, and depending on the installation location, platform geometry, time of day (sun location), and associated platform shading of the water target, only one or two (or sometimes zero) of these azimuth angles are appropriate for measurement of water reflectance; other azimuth angles will be contaminated by platform shading or even direct obstruction of the water target as defined from the instrument field of view. A selection of acceptable azimuth angles is made a priori based on expert judgement (used here for prototype testing; see Section 3) or, better, a study of platform shading effects by modeling or experimentation [27]. The measurement of unacceptable azimuth angles, defined by a "no-go zone", can then be avoided to save power.

For each measurement cycle, the system performs a sub-cycle for each of the configured relative azimuth angles as described in Table 1**.** Based on the AERONET-OC protocol [6], but with repetition of the *Ed* and *Ld* replicates, each azimuthal measurement sub-cycle consists of 2 × 3 replicate scans each of downwelling irradiance, *Ed*, and sky radiance, *Ld*, and 11 replicate scans of upwelling radiance, *Lu*, where "(spectral) scan" refers to acquisition of a single instantaneous spectrum. Firstly, the irradiance sensor is pointed upward, with the radiance sensor offset by 40◦, and three replicates of *Ed* followed by three replicates of *Ld* are measured. The radiance sensor is then moved to a 40◦ downward viewing angle to make 11 replicate *Lu* scans. The irradiance and radiance sensors are then repositioned to make

three more replicate scans each of *Ld* and *Ed*. Scans are stored as digital counts (DCs), including the embedded dark pixel counts, with instrument serial number and integration time in the metadata.

**Table 1.** Pan-and-tilt hyperspectral radiometer system (PANTHYR) basic data acquisition cycle. Azimuth is measured clockwise with respect to sun with a relative azimuth of 0◦, meaning that the radiance sensor is pointing toward the sun, and a relative azimuth of 90◦, meaning that the radiance sensor is pointing to the right of sun. A zenith angle of 180◦ means that the instrument is pointing vertically upward and, hence, measuring downwelling (ir)radiance. Measurements within a predefined azimuth no-go zone will be skipped.


#### *2.5. PANTHYR Data Processing*

Each of the scans within a sub-cycle are converted from DCs to calibrated (ir)radiances, as described in Equations (1–3) of Reference [20], using the calibration files appropriate for the instrument serial number and date as follows:


Incomplete spectra are removed, as well as *Ld* and *Lu* scans with >25% difference between neighboring scans at 550 nm. *Ed* scans are removed using the same criterion after normalizing *Ed* by *Cos*(θ0), where θ<sup>0</sup> is the sun zenith angle. The sub-cycle is further processed if sufficient scans pass the quality control: *Lu* 9/11, and *Ed* and *Ld* 5/6 each. The remaining *Ed* and *Ld* measurements are then grouped and mean averaged to *Ed* and *Ld*. For each *Lu* scan, the water-leaving radiance, *Lw*, is computed by removing skyglint radiance as follows:

$$L\_{\mathfrak{w}} = L\_{\mathfrak{u}} - \rho \overline{L\_{\mathfrak{u}}} \tag{1}$$

where ρ*<sup>F</sup>* is the "wind-roughened Fresnel" coefficient that represents the fraction of incident skylight that is reflected back toward the water-viewing sensor at the air–water interface and is given by the look-up table (LUT) of Reference [25]. This LUT describes ρ*<sup>F</sup>* as function of viewing and sun geometry and wind speed. Wind speed is retrieved from ancillary data files or set to a user-defined default value if wind speed data are unavailable. The data in the LUT are linearly interpolated to the current observation geometry and wind speed.

The *Lw* scans are then converted into ("uncorrected") water-leaving radiance reflectance scans, ρ*wu*, as follows:

$$
\rho\_{\text{new}} = \pi \frac{L\_w}{\overline{E\_d}}.\tag{2}
$$

A "near infrared (NIR) similarity spectrum" correction is then applied to remove any white error from inadequate skyglint correction as shown in Equation (3). A spectrally flat measurement error, ε, is estimated using two wavelengths in the NIR [28], where λ<sup>1</sup> = 780 nm and λ<sup>2</sup> = 870 nm.

$$
\varepsilon = \frac{\alpha \rho\_{\rm un}(\lambda\_2) - \rho\_{\rm un}(\lambda\_1)}{\lambda\_2 - \lambda\_1},
\tag{3}
$$

where α is the similarity spectrum [26] ratio for the bands used; here, α(780, 870) = 1/0.523 = 1.912. Per scan, ε is subtracted from the ρ*<sup>w</sup>* at all wavelengths to give an NIR-corrected water reflectance,

$$
\rho\_{\rm w}(\lambda) = \rho\_{\rm wu}(\lambda) - \varepsilon,\tag{4}
$$

and all ρ*<sup>w</sup>* scans are mean-averaged to give the final NIR-corrected water reflectance, ρ*w*. This NIR correction is optional and may be turned off, e.g., for extremely turbid waters where it is not valid and/or for situations where the satellite measurement may include a similar constraint in the atmospheric correction which needs to be validated independently.

The final quality control to retain or reject the ρ*<sup>w</sup>* is performed according to Reference [26]; measurements are rejected when *Ld*/*Ed* > 0.05 sr−<sup>1</sup> at 750 nm (indicating clouds either in front of the sun or in the sky-viewing direction), and when the coefficient of variation (CV, standard deviation divided by the mean) of the ρ*<sup>w</sup>* scans is >10% at 780 nm.

#### **3. Prototype Testing**

#### *3.1. July 2018 Deployment at the Acqua Alta Oceanographic Tower*

After extensive testing of individual components and the full system in the laboratory and outdoors on land, a first prototype system test was performed on the Acqua Alta Oceanographic Tower (AAOT) in July 2018 during the international Field Inter-Comparison Exercise (FICE) organized by the FRM4SOC (Fiducial Reference Measurements for Satellite Ocean Colour) project. AAOT is an oceanographic platform in the Northern Adriatic Sea 15 km offshore from Venice in water of 16 m depth. It was used extensively over the last 20 years for oceanographic data collection including optical oceanography [29] and was used for a number of multi-partner optical radiometry inter-comparison exercises [30,31]. It is the location of the first AERONET-OC site and provides such data since 2002. The above-water platform structure was entirely reconstructed in 2018, but the new structure closely follows the design and layout of the original platform.

The PANTHYR system was deployed on the east corner of the top deck of the platform (Figure 4), with the irradiance sensor collector 2 m above top deck floor and, hence, about 14 m above sea level. The system ran for six days continuously from 12 to 17 July 2018, under the supervision of the developer for the first three days and with non-specialized supervision for the remaining days. Some test conditions were not the most challenging that can be expected in the future; grid power was available, and data were not transmitted over a 4G link but stored onboard. However, the basic functions of autonomous pointing and data acquisition were fully and successfully tested.

This deployment for prototype testing was also installed only at moderate height and not above the height of the top deck masts for practical reasons; for future operational deployments, the PANTHYR system should ideally be positioned above the height of all such masts and structures or with a careful analysis of any shading of downwelling irradiance. For subsequent data comparison, only PANTHYR measurements at relative azimuth angles of 225◦ or 270◦ were used.

Using the angular field-of-view data (7◦ full-width half-maximum) for the radiance instrument supplied by the manufacturer suggests that the water-viewing radiometer sees an elliptic patch of water with about 3 m diameter. Typical integration time for the water-viewing measurement is ≤512 ms during this experiment. Spatio-temporal variability caused by the surface wave field is, thus, mainly averaged and not resolved.

The FICE campaign further involved the deployment of multiple supervised radiometer systems, including seven above-water radiometer systems, one shipborne free-fall underwater radiometer system, and one platform-deployed vertical profiling system. Inter-comparison of all these systems will be the focus of a separate paper. In the current paper, results of the PANTHYR system are compared only to one supervised above-water radiometer system, the manually supervised RBINS (Royal Belgian Institute of Natural Sciences) hyperspectral above-water three-radiometer TRIOS/RAMSES system, hereafter called the "M3TRIOS" system, and the autonomous multispectral AERONET-OC system operated by the Joint Research Centre (JRC) of the European Commission at Ispra. The five TRIOS/RAMSES sensors used here (two for PANTHYR, three for M3TRIOS) were radiometrically calibrated at the Tartu Observatory laboratory just before the FICE campaign, thus minimizing calibration-related differences.

Only measurements from the five days with best sky conditions (13–17 July) are used here. For those days and during the period of day used for comparison in the present study, sun zenith angles ranged from 23.4◦ to 46.7◦ and wind speed ranged from 0.1 to 4.3 m/s.

**F G** 

**Figure 4.** Photographs of PANTHYR prototype deployment at Acqua Alta Oceanographic Tower (AAOT). (**a**) Irradiance and sky radiance measurement position; (**b**) southwest face of the Acqua Alta Oceanographic Tower; (**c**) upwelling (water + skyglint) radiance measurement position; (**d**) park position.

#### *3.2. Manually Supervised M3TRIOS System Used for Data Comparison*

The M3TRIOS system deployed during the FRM4SOC/FICE consists of three TRIOS RAMSES radiometers, one measuring downwelling irradiance, *Ed*, and two measuring radiance with fixed skyand water-viewing zenith angles of 140◦ and 40◦ for *Ld* and *Lu*, respectively. The system was deployed since 2001 at many locations from both large and small ships, and platforms including AAOT [31]. Details of the instrumentation and standard measurement protocol, based on Reference [25], and data processing and quality control can be found in Reference [26]. This is quite similar to the PANTHYR processing described in Section 2.5, except for the following elements:


Although the instrument has wider spectral range, data are here limited to the range 400–900 nm where quality was checked by previous inter-comparison exercises. For this specific implementation, the downwelling irradiance sensor was mounted on the top deck of AAOT on a telescopic mast at height 5 m above deck. The sky- and water-viewing radiance sensors were mounted at the southwest face of the platform on a frame tailor-made by Tartu Observatory and Plymouth Marine Laboratory to accommodate many radiometers and to ensure 40◦ and 140◦ zenith-angle viewing. On the other hand, 90◦, 135◦, or 270◦ relative azimuth angle to sun was achieved by manual rotation of the structure before each measurement. Measurements were made typically every 30 min during daytime between 8:00 a.m. and 1:00 p.m. UTC.

#### *3.3. AERONET-OC Data Used for Comparison*

The CIMEL/Seaprism system providing data for AAOT through the AERONET-OC network was installed on the west corner of AAOT on a purpose-built jetty to minimize any shading effects of the platform on the water being measured. Full details on the instrumentation and measurement protocol and data processing and quality control can be found in Reference [6].

Level 1.5 normalized water-leaving radiance data, *Lwn*, were used without *f*/*Q* correction [32] as it is not applied to the PANTHYR andM3TRIOS data. This data were downloaded from the AERONET-OC (AOC) website for the "Venise" site located on AAOT [6]. Level 1.5 data are automatically cloud cleared but do not have final post-deployment calibration applied and are, hence, not fully quality assured. The matchup data used in this paper do, however, pass quality control and are likely candidates for level 2.0 data if post-deployment calibration is acceptable [Zibordi, pers. comm. 18 March 2019]. *Lwn* values were converted to water reflectance, ρ*w*, by

$$
\rho\_{\text{i}\nu} = \pi \frac{L\_{\text{awn}}}{F\_0} \, \text{s} \tag{5}
$$

where *F*<sup>0</sup> is the extraterrestrial irradiance [33] resampled to gaussian band-passes with 10-nm full-width half-maximum centered on the reported exact center wavelengths of the CIMEL instrument. Each of the level 1.5 AOC measurements was matched to the closest PANTHYR measurement within 20 min, and the PANTHYR data were resampled to the same band.

#### **4. Results**

#### *4.1. System Performance*

The system tested at AAOT performed very well and demonstrated its capability of making automated measurements following a predefined protocol and interval. After installation and configuration, there were no blocking issues while the system worked autonomously, and no intervention was required during the six days of operation while the system performed as expected.

Since the system was still in a prototype phase when tested on the AAOT, some functionality is yet to be tested or demonstrated in a real-world scenario. The system was powered from the mains grid and, as the platform was visited daily, there was no need for 3G/4G data transmission. Most measurement sites will require one or both of these functions as mains power or cabled internet may not be available.

With the defined measurement protocol and a measurement interval of 20 min, it became obvious that the system speed needed improvements. The measurement protocol described in Section 2.4 took about 12 min to finish, depending on sun azimuth (since this affects required movements and how many measurements are skipped because they fall within the "no-go" zone). Post-FICE optimization of the control libraries for the pan/tilt head since resulted in faster movements. Additional speed can be gained in the control of the instruments, as well as by implementing multi-threading. Before each measurement cycle, some of the devices need time to start up and/or calibrate. The pan/tilt head, for example, needs to run a calibration routine to find the center point of both axes, and the IP camera needs time to boot and warm up. Implementation of multithreading would allow us to do both at the same time. Shorter measurement cycles result in more low-power standby time in between measurement cycles and lower power consumption, thus achieving better overall efficiency.

Some existing systems such as the CIMEL SeaPRISM hardware rely on the sun to calibrate the system azimuth during installation. While a high accuracy can be achieved with this method, it is only usable on sunny days. An electronic compass that achieves the required accuracy on metal structures, and that has a low power consumption, a marine-grade enclosure, and low cost is yet to be found. The current mode of operation for aligning the system zenith is achieved by placing a spirit level on the top bracket. This allows alignment as accurate as the used spirit level. After this is done, a digital inclinometer is placed on top of the instruments to verify their absolute inclination for different angles. Azimuth calibration is similar to the CIMEL system; the user can order the head to rotate toward the sun, after which the shadows cast by the top bracket serve as a reference while manually performing the alignment.

The images taken as part of each measurement cycle showed an unexpected interest from the local wildlife. Even though the system was (almost) continuously moving, birds seemed not to be deterred but showed a rather close interest in the new technology, both contaminating data by blocking the instrument field of view completely or subtly (Figure 5) and by potentially contaminating the instrument itself by fecal deposits. Major blocking of a radiometer field of view will be easily identified as bad data because the spectral signature will be different from the expected water, sky, and full sky irradiance targets. However, partial blocking of the instrument field of view (Figure 5b) may contaminate data in a way that cannot be automatically detected and rejected. Discussion with experts on bird life suggests that typical visible or audible "scarecrow" devices are not effective deterrents and the best approach is the use of spikes to prevent comfortable resting spots. We, therefore, strapped cable ties around the moving parts (Figure 4c), leaving the ends uncut as deterrent spikes, being careful to avoid the irradiance sensor field of view.

The integration of an IP camera alongside the radiometers proved particularly useful in identifying these unforeseen problems and may be useful for other causes of data contamination such as boats, floating debris/vegetation, or other unusual conditions in the radiometer field of view (Figure 5) since the radiometer gives only spatially integrated information. However, image analysis is presently subjective and is not easy to automate.

**Figure 5.** Sample photos taken from the camera aligned with the radiance sensor (for azimuth angles not used for measurements in this paper). (**a**) Bird blocking field of view; (**b**) more subtle measurement contamination; (**c**) floating vegetation and platform shadow on water; (**d**) unexpected water target contamination.

#### *4.2. Water Reflectance Spectra—Mean over Time*

For the three clear sky days analyzed here, 30 matchups were recorded between the PANTHYR and M3TRIOS systems and 10 between the PANTHYR and AERONET-OC systems. All 30 water reflectance spectra for PANTHYR and M3TRIOS systems are plotted in Figure 6, with the mean average of the 10 matchup AERONET-OC spectra. Differences between spectra observed on different days/times by a single system in this figure combine both possible temporal variability of the target itself (including possible BRDF effects, which are uncorrected for the PANTHYR and M3TRIOS systems), the random uncertainties associated with the measurement protocol and associated data processing, and certain instrument artefacts (e.g., imperfect irradiance sensor angular response). Since the intra-spectra differences are rather small, we can conclude that the water target was itself rather stable during the three days of analyzed data (spread over a six-day period). One outlier is clearly visible for the M3TRIOS system; detailed analysis of data for that measurement, including the manual log-book and various photos, suggests that this may result from an oily film visible at the water surface within the M3TRIOS system field of view. After exclusion of this outlier, the mean average water reflectance spectrum over the 29 matchups is very similar for PANTHYR and M3TRIOS systems. The mean average water reflectance spectrum over the 10 matchups is also very similar for the PANTHYR and the AERONET-OC systems. For a more detailed understanding of system data quality, including the impact of random uncertainties and removing the differences associated with water target temporal variability, the individual matchups were analyzed as described in the sections below.

#### *4.3. Data Comparison with M3TRIOS System—Matchup Analysis*

The 29 matchups between the PANTHYR and M3TRIOS systems, interpolated to typical multispectral wavelengths, are shown via scatterplots of *Ed*(λ), *Ld*(λ), *Lu*(λ), and ρ*w*(λ), and associated linear regression statistics in Figure 7, including root-mean-square difference (RMSD) and mean average relative difference (MARD).

**Figure 6.** Water reflectance spectra for each of the 30 stations (gray lines) from (left) the three-radiometer TRIOS/RAMSES system (M3TRIOS) system and (right) the PANTHYR system. Mean spectra over all 30 stations are superimposed for the PANTHYR (orange dashed line) and the M3TRIOS system (blue dashed lines) and for all 10 measurements from the AERONET-OC system (green dotted line joining dots where the multispectral data exist).

The MARD of *Ed*(λ) calculated over all wavelengths was 3.1% and the scatterplot (Figure 7a) suggests some spectral variability of differences with lower PANTHYR *Ed*(λ) for 410 nm and 440 nm and higher PANTHYR *Ed*(λ) at 620 nm and 675 nm. A full uncertainty analysis is not yet available for these instruments/systems as deployed at AAOT, and no corrections were applied for instrument artefacts such as imperfect cosine response, straylight, thermal sensitivity, non-linearity or polarization effects, or deployment imperfections, e.g., for optical perturbations from higher mast structures. Operational deployment of the PANTHYR system for fiducial reference measurements will require such an uncertainty analysis but is beyond the scope of this technology-proving paper. However, the observed differences for *Ed*(λ) are not a cause for concern at this stage of testing. The reader is referred also to Reference [34] for a more detailed discussion of differences between measurements of *Ed*(λ) from different instruments in typical field conditions.

The scatterplot (Figure 7b) of *Ld*(λ) sky radiance measurements shows strong correlation but for two distinct groups of data. Analysis of the metadata recorded for both systems revealed that these distinct groups of data correspond to different azimuth angles of the PANTHYR and M3TRIOS systems; they are pointing at very different portions of the sky and, in some cases, the PANTHYR system is pointing at a much brighter sky with cirrus clouds. Since both of the obvious two groups of data corresponded to cases where PANTHYR and M3TRIOS were measuring at different azimuth angles because the systems were deployed on different sides of the AAOT platform and, hence, pointed at the water in different azimuth angles, there is no cause for concern that neither group of points has a slope of one. There were only two matchups where PANTHYR and M3TRIOS were measuring at the same azimuth angle, but this was considered insufficient for detailed analysis.

**Figure 7.** Scatterplot for selected wavelengths comparing PANTHYR data with M3TRIOS system data for (**a**) *Ed*, (**b**) *Ld*, (**c**) *Lu*, and (**d**) ρ*<sup>w</sup>* for 29 matchup stations. Linear regression statistics are given, including the root-mean-square difference (RMSD) and the mean average relative difference (MARD).

The scatterplot (Figure 7c) of *Lu*(λ) water (plus reflected skyglint) radiance measurements shows strong correlation (*r*<sup>2</sup> = 0.992) with low systematic differences (slope = 0.97). The MARD of 6.5% was affected by the different azimuth angles of the two systems since this parameter includes reflected skyglint, typically 2.5–3.0% of *Ld*(λ).

The scatterplot (Figure 7d) of ρ*w*(λ) water reflectance (after skyglint correction) shows high correlation between the PANTHYR and M3TRIOS systems (*r*<sup>2</sup> = 0.996) with low systematic differences (slope = 0.97). The MARD of 7.1% was dominated by the contribution for data at 870 nm, where the water reflectance was very low and absolute differences, e.g., arising from skyglint correction, translated into very large relative differences.

The spectral variations of RMSD and MARD between PANTHYR and M3TRIOS for *Ed*(λ) and ρ*w*(λ) over these 29 matchups are shown in Figure 8. The short wavelength scale variability for the RMSD of *Ed*(λ) in Figure 8a was related to the similar spectral variability of *Ed*(λ) itself, e.g., the atmospheric oxygen absorption feature at 762 nm or the Frauenhofer lines related to solar photosphere absorption at 516–518 nm, and the way the slightly different central wavelengths of each *Ed* instrument under-resolved these features. The MARD of *Ed*(λ) between PANTHYR and M3TRIOS (Figure 8b) was <5% for the full spectral range 400–900 nm. The common radiometric calibration of the two *Ed* instruments just prior to this experiment helped limit MARD; however, a more detailed characterization of these two instruments, particularly including an analysis of their "cosine" angular response and perhaps straylight and non-linear responses, may improve results in the future.

The spectral variation of RMSD between PANTHYR and M3TRIOS for ρ*w*(λ) shown in Figure 8c suggests an overall spectral shape typical of *Lu*(λ) spectra, i.e., upwelling radiance from green water and reflected "blue" skyglint. Possible causes of this difference are multiple and certainly include imperfect correction of skyglint, but may also include different BRDF at the different azimuth angles, propagation of *Ed*(λ) differences, *L* radiometer calibration and characterization (straylight, non-linearity, polarization, thermal sensitivity, etc.), etc. A detailed analysis is beyond the scope of this technology-proving paper, which concludes merely that there is no major cause for concern at present.

The MARD between PANTHYR and M3TRIOS shown in Figure 8d for ρ*w*(λ) was <5% for the spectral range 410–580 nm, with higher MARD around 600 nm. The much higher MARD for wavelengths higher than 700 nm are clearly related to the very low ρ*w*(λ) for these waters, showing some similarity with the pure water absorption coefficient spectrum for 700–900 nm [26] and with the phytoplankton absorption coefficient at 660–680 nm. In such conditions, the MARD has no practical relevance because a satellite validation analysis would generally consider absolute differences between PANTHYR and satellite or simply not use such data.

**Figure 8.** Root-mean-square difference (RMSD) and the mean average relative difference (MARD) between PANTHYR and M3TRIOS systems over 29 stations for (**a**,**b**) *Ed*, and (**c**,**d**) ρ*w*. The 5% MARD is shown as a horizontal dashed line.

*4.4. Data Comparison with AERONET-OC System—Matchup Analysis*

The 10 matchups between the PANTHYR and AERONET-OC systems, with resampling of PANTHYR data to gaussian spectral response functions on the AERONET-OC central wavelengths, are shown via scatterplots of ρ*w*(λ), and associated linear regression statistics in Figure 9. There is a strong correlation between these datasets (*r*<sup>2</sup> = 0.996) with low systematic differences (slope = 0.95). The MARD over all wavelengths appeared high (18.3%); however, as seen in the spectral variation of MARD shown in Figure 8d, it was dominated by the data at 870 nm, where RMSD was low but ρ*w*(λ) itself was very low. For the wavelengths from 413 nm to 555 nm, MARD between PANTHYR and AERONET-OC was between 3.5% and 6.5%, except at 530 nm, where it reached 9.0%. Possible causes of the differences are multiple and may include imperfect correction of skyglint, propagation of *Ed*(λ) differences (not available online for AERONET-OC), radiometer calibration and characterization (straylight and/or spectral response functions, non-linearity, polarization, thermal sensitivity, etc.), etc. A detailed analysis is beyond the scope of this technology-proving paper, which concludes merely that there is no major cause for concern at present.

**Figure 9.** Scatterplot for selected wavelengths comparing PANTHYR data with AERONET-OC Venise Level 1.5 data for ρ*w* for 10 matchups. Linear regression statistics are given, including the root-mean-square difference (RMSD) and the mean average relative difference (MARD). All PANTHYR data presented in this figure were measured at the 270◦ relative azimuth. The black line is the reduced major axis regression. Statistics are plotted here for all wavelengths together. For per-band RMSD and MARD, see Figure 8c,d.

#### **5. Conclusions and Future Perspectives**

To gather enough high-quality ground-truth data for validation of water reflectance derived from satellite missions, an automated system is necessary. The AERONET-OC federated network of autonomous instrument systems [6] is now the main source of such validation data for all satellite missions measuring water reflectance but provides only multispectral data. The PANTHYR system described here aims to provide hyperspectral water reflectance data for satellite validation.

The design of the PANTHYR system was described in detail here. Two COTS hyperspectral radiometers (one radiance, one irradiance) and an IP camera are mounted on a COTS PT pointing system with original development of control electronics and software providing a low-cost, low-power, but robust, modular, and extendable design. Date/time/location information is received from a GNSS receiver and data are transmitted by a 4G gateway. A flexible mounting bracket was designed to easily fit the system to any suitable platform structure. The prototype was tested during an inter-comparison exercise organized at AAOT, an offshore platform in the Adriatic Sea. The system functioned autonomously over the six-day period without intervention. Data acquired from the system was compared with an established manually deployed hyperspectral system and with the automated AERONET-OC system data deployed simultaneously at AAOT. Data quality from the PANTHYR system was good; water reflectances compared to the two other systems had per-band MARDs within 5.5% for the spectral range 410–550 nm in these good measurement conditions.

Imagery from the IP camera was surprisingly useful and revealed unexpected conditions including birds within the camera field of view, which could contaminate radiometer data in significant (easily detected) or subtle (potentially undetected) ways.

The prototype tests described here were made with mains power supply and without autonomous data transmission over the internet. Use of an autonomous power supply (typically solar panels) with power supply monitoring and automated data transmission over 4G need to be tested in future work.

An upgrade of the PT unit to a version with slip rings is also planned to avoid the risk of cable snags possible with this first prototype.

Future mounting hardware will be machined from 5083 aluminum alloy instead of stainless steel to save weight. Lower overall weight results in lower shipping cost, and limiting the weight on top of the PT head reduces power consumption.

The system can also be extended to function on moving platforms, such as ships or buoys, via the addition of heading and inclination sensors and/or movement compensation mechanisms.

The data acquisition and processing described here are based strongly on precursor work and provide a robust starting point. However, future research may take advantage of the flexible pointing capability to investigate improved data acquisition protocols, e.g., with multiple zenith angle radiance measurements. The skyglint correction for above-water reflectance measurements is also considered to be a major source of measurement uncertainty, and potential improvements were investigated in many recent studies. When considered mature, such improvements can easily be incorporated in the processing software described here.

In conclusion, the PANTHYR system prototype was successfully tested in a basic configuration (without autonomous power supply and without data transmission over the internet) giving robust performance and high-quality hyperspectral data. The system prototype meets the requirements for future worldwide deployment in a network for hyperspectral validation of water reflectance data from satellites.

**Author Contributions:** Conceptualization, D.V., K.R., and A.C.; methodology, D.V., K.R., and A.C.; software, D.V. and Q.V.; validation, Q.V. and M.B.; formal analysis, K.R., Q.V., and M.B.; writing—original draft preparation, K.R. and D.V.; writing—review and editing, A.C., Q.V., M.B., and D.V.

**Funding:** This research was funded by BELSPO (Belgian Science Policy Office) in the framework of the STEREO III program project HYPERMAQ (SR/00/335), and by the European Space Agency PRODEX/HYPERNET-OC project. The AAOT fieldwork was partially funded by the European Space Agency FRM4SOC project.

**Acknowledgments:** The authors would like to thank the Plymouth Marine Laboratory (PML) and Consiglio Nazionale delle Ricerche (CNR), Dr. Gavin Tilstone, Tartu Observatory, and the FRM4SOC project for organizing the FICE at AAOT and for providing mounting frames for the M3TRIOS system radiometers. Tartu Observatory provided the radiometric calibration of all five RAMSES instruments used here. Dr. Giuseppe Zibordi is acknowledged for providing the AERONET-OC data and for kind advice. TRIOS Gmbh is acknowledged for providing details of the serial protocol communication for the RAMSES instruments, allowing us to write the necessary data acquisition software. For the latter, the RFLEX code provided by Stefan Simis provided some useful ideas on how to proceed. Michiel T'Jampens provided appreciated technical feedback. We thank the three anonymous reviewers for many constructive comments on the first version of this paper.

**Conflicts of Interest:** The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

#### **References**


© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
