# Improved Vegetation Profiles with GOCI Imagery Using Optimized BRDF Composite

Authors: Sang-il Kim; Do-Seob Ahn; Kyung-Soo Han; Jong-Min Yeom
Journal: Journal of Sensors
 (2016)

---

## Abstract

The purpose of this study was to optimize a composite method for the Geostationary Ocean Color Imager (GOCI), which is the first geostationary ocean color sensor in the world. Before interpreting the sensitivity of each composite with ground measurements, we evaluated the accuracy of bidirectional reflectance distribution function (BRDF) performance by comparing modeled surface reflectance from BRDF simulation with GOCI-measured surface reflectance according to composite period. The root mean square error values for modeled and measured surface reflectance showed reasonable accuracy for all of composite days since each BRDF composite period includes at least seven cloud-free angular sampling for all BRDF performances. Also, GOCI-BRDF-adjusted NDVIs with four different composite periods were compared with field-observation NDVI and we interpreted the sensitivity of temporal crop dynamics of GOCI-BRDF-adjusted NDVIs. The results showed that vegetation index seasonal profiles appeared similar to vegetation growth curves in both field observations from crop scans and GOCI normalized difference vegetation index (NDVI) data. Finally, we showed that a 12-day composite period was optimal in terms of BRDF simulation accuracy, surface coverage, and real-time sensitivity.


---

## Body

## 1. Introduction

Terrestrial products such as surface albedo and vegetation index play an important role in soil moisture and energy balance and are critical for understanding the boundary layer energy mechanism between the atmosphere and land surface [1, 2]. Satellite imagery is considered an efficient tool for estimating terrestrial products, having the advantage of a worldwide spatial and temporal resolution. Surface vegetation conditions and dynamics in terrestrial products [3] have mostly been constructed from optical satellites such as the Advanced Very High Resolution Radiometer (AVHRR) [4], Moderate-Resolution Imaging Spectroradiometer (MODIS) [4–6], SPOT/VEGETATION [7], METEOSAT Second Generation (MSG) [8], Geostationary Operational Environmental Satellite (GOES) [9], and Medium Resolution Imaging Spectrometer (MERIS) [10]. Atmospheric effects and surface anisotropy scattering are the main obstacles to accurate estimates when using optical satellites to measure terrestrial products. Although atmospheric correction models for the most common optical satellites are easily accessible using commercial remote sensing software, the bidirectional reflectance distribution function (BRDF) has a complicated methodology, which may restrict its utility among user groups. However, BRDF is a critical tool for reliably reducing the anisotropy associated with surface scattering effects, one of the main fluctuating error sources of relative solar-target-sensor geometry [11].The most widely used semiempirical BRDF model for optical sensors is sensitive to the number of observations made under geometric conditions, and observations of a diversity of conditions have resulted in a higher quality BRDF product [12]. Sun- and Geosynchronous satellites have different angular sampling observations [13–15]. A typical MODIS among polar-orbiting satellites performs BRDF modeling using a 16-day composite period to provide a suitable tradeoff between the availability of sufficient angular sampling observations and the temporal stability of the surface [16, 17]. Geostationary satellites ensure repeatable sun illumination conditions during a given day, acquiring sufficient cloud-free surface reflectance to make a reliable BRDF model [18].The Geostationary Ocean Color Imager (GOCI), the first geostationary optical ocean color sensor, has the potential to accurately estimate land surface products [19–22], because it provides images with high temporal resolution and has a spectral response function (SRF) over the solar spectral wavelength range, similar to MODIS sensors. SeaWiFs (Sea-Viewing Wide Field-of-View Sensor), another ocean color sensor, has previously been used for estimating terrestrial vegetation [23, 24].It is impossible to obtain land surface spectral information under a cloudy area with optical satellite sensors, so composite days must be generated. The most commonly used 16-day composite for BRDF modeling in geostationary satellites [20, 25] shows reduced sensitivity to crop dynamics over selected rice paddy areas relative to real-time conditions. The sensitivity of BRDF models is reduced with longer composite days due to the additional time involved, despite adequate cloud-free sampling. Yeom and Kim [21], when comparing in situ Normalized Distribution Vegetation Index (NDVI) with GOCI-BRDF-adjusted NDVI, showed that the BRDF-adjusted NDVI profile was shifted relative to the in situ profile. These adjusted vegetation profiles are thus less sensitive to real-time change due to the composite period used for the BRDF model [21]. Therefore, it is critical to optimize the BRDF composite method for sensitivity to vegetation profiles, with a choice-adaptive composite period.


## 2. Materials and Methods

### 2.1. Study Area and Satellite Data

The study was conducted in the southern part of the Korean peninsula (Figure1). The Korean peninsula is characterized by having four distinct seasons. The monsoon season, called Jang-Ma, brings considerable amounts of rain to the Korean peninsula in the summer between June and August, during which the Korean peninsula receives more than half of its annual precipitation. The annual mean temperature in South Korea ranges from 8.4 to 18.6°C, except in the high mountain areas (Korea Meteorological Administration) [26].Figure 1
Geostationary Ocean Color Imager (GOCI) full disk image acquired on 11 September 2014, over Korea. The blue box marks the study area. Field-observation sites are denoted by red points (a: Kyehwa, b: Kimjae).The GOCI was launched successfully on 27 June 2010, with a 2500 km × 2500 km field of view and eight multispectral bands, including visible to near-infrared bands. GOCI has a 500 m spatial resolution eight times a day, with a reobservation time of around 1 hour from 9 a.m. to 4 p.m. The characteristics of GOCI are presented in Table1. We used bands 5 and 8 of NDVI, due to the similar spectral characteristics of these bands with respect to those of MODIS.Table 1
Detailed characteristics of the Geostationary Ocean Color Imager (GOCI) and Moderate-Resolution Imaging Spectroradiometer (MODIS) sensors used for estimating land-surface products.
Satellite
GOCI
Orbit type
Geosynchronous
Altitude
≈36,000 km
Spatial resolution
500 m at nadir
Wavelength (nm)
B1: 402–422
B5: 650–670
B2: 433–453
B6: 675–685
B3: 480–500
B7: 735–755
B4: 545–565
B8: 845–885


### 2.2. Field Observation Using a Multispectral Radiometer (MSR)

Field measurements were obtained with a handheld portable CROPSCAN multispectral radiometer (MSR) to compare with satellite-based vegetation profiles. The MSR-16 consists of 16 spectral bands over the wavelength range of 450–1750 nm and a measurement distance of 2 m from the sensor to the vegetation canopy. The MSR-16 measurements are unaffected by atmospheric conditions. To measure spectral information on rice paddies, we observed three different points within two sites under clear sky conditions (red points in Kyehwa and Kimjae, Figure1) and then calculated the average of three points at each site (Figure 2). The instrument was operated during the growing season (13 June 2014, 26 June 2014, 15 July 2014, 22 July 2014, 11 August 2014, 26 August 2014, 15 September 2014, and 3 October 2014). Field observations were compared with satellite observation pixels to geometrically match observation points (500 × 500 m).Figure 2
Detailed field-observation site images acquired on 5 October 2014, using the Landsat-8 operational land imager. The red dots are field-observation sites (a and b from Figure1), and the rectangles represent sites corresponding to one GOCI pixel.
(a)
(b)


### 2.3. Preprocessing of Satellite Data and BRDF Modeling

GOCI imagery was preprocessed before BRDF modeling. The preprocessing was included in converting the digital number to top-of-atmosphere (TOA) reflectance, cloud masking, and atmospheric correction. Cloud masking used a threshold method to reduce the probability of cloudy contamination [27]. An atmospheric correction model was applied to the TOA reflectance to estimate top-of-canopy reflectance using a look-up table based on the Second Simulation of the Satellite Signal in the Solar Spectrum [28–31].To correct for surface anisotropy effects, we used the BRDF model developed by Roujean et al. [32]; in the semiempirical BRDF model proposed, there are three scattering components: isotropic scattering, geometric scattering, and volumetric scattering from a horizontally homogeneous vegetation canopy [33]. Consider the following:(1)
ρ
θ
s
,
θ
v
,
∅
=
k
0
+
k
1
f
1
θ
s
,
θ
v
,
∅
+
k
2
f
2
θ
s
,
θ
v
,
∅
,where ρ is the observed reflectance in a given spectral channel, θ
s is the solar zenith angle (SZA), θ
v is the viewing zenith angle (VZA), ∅ is the relative azimuth angle (RAA), k
0 is isotropic coefficient of the Lambertian reflectance when θ
s
=
0 andθ
v
=
0, k
1 and k
2 are the empirical coefficients of the geometric and volumetric kernel, and f
1 and f
2 are physical geometric and volumetric scattering kernels depending on θ
s, θ
v, and ∅ between the sun and view direction of view. We applied our own fixed viewing angle and mean SZA during the composition period, without adjusting the VZA to the nadir direction based on Yeom and Kim [20].Inadequate observation sampling leads to poor surface coverage, while long composite periods reduce BRDF model sensitivity to real-time vegetation change, even with adequate cloud-free observations. To improve the vegetation profiles estimated with BRDF composite methods, several composite periods were added to enhance the accuracy and sensitivity of BRDF modeling relative to ground measurements. Here, we used 4-, 8-, 12-, and 16-day composite periods for each pixel to determine the optimized GOCI-based, BRDF-adjusted NDVI for each study area. Figure3 shows the concept of alterative composite periods for BRDF modeling. To perform kernel-driven, semiempirical BRDF modeling, at least seven observations of clear free surface reflectance are required for reliable simulation results [15, 34]. Finally, we estimated the daily NDVI product-based BRDF-adjusted surface reflectance (BAR).Figure 3
Concept of composite period for BRDF modeling.


## 2.1. Study Area and Satellite Data

The study was conducted in the southern part of the Korean peninsula (Figure1). The Korean peninsula is characterized by having four distinct seasons. The monsoon season, called Jang-Ma, brings considerable amounts of rain to the Korean peninsula in the summer between June and August, during which the Korean peninsula receives more than half of its annual precipitation. The annual mean temperature in South Korea ranges from 8.4 to 18.6°C, except in the high mountain areas (Korea Meteorological Administration) [26].Figure 1
Geostationary Ocean Color Imager (GOCI) full disk image acquired on 11 September 2014, over Korea. The blue box marks the study area. Field-observation sites are denoted by red points (a: Kyehwa, b: Kimjae).The GOCI was launched successfully on 27 June 2010, with a 2500 km × 2500 km field of view and eight multispectral bands, including visible to near-infrared bands. GOCI has a 500 m spatial resolution eight times a day, with a reobservation time of around 1 hour from 9 a.m. to 4 p.m. The characteristics of GOCI are presented in Table1. We used bands 5 and 8 of NDVI, due to the similar spectral characteristics of these bands with respect to those of MODIS.Table 1
Detailed characteristics of the Geostationary Ocean Color Imager (GOCI) and Moderate-Resolution Imaging Spectroradiometer (MODIS) sensors used for estimating land-surface products.
Satellite
GOCI
Orbit type
Geosynchronous
Altitude
≈36,000 km
Spatial resolution
500 m at nadir
Wavelength (nm)
B1: 402–422
B5: 650–670
B2: 433–453
B6: 675–685
B3: 480–500
B7: 735–755
B4: 545–565
B8: 845–885


## 2.2. Field Observation Using a Multispectral Radiometer (MSR)

Field measurements were obtained with a handheld portable CROPSCAN multispectral radiometer (MSR) to compare with satellite-based vegetation profiles. The MSR-16 consists of 16 spectral bands over the wavelength range of 450–1750 nm and a measurement distance of 2 m from the sensor to the vegetation canopy. The MSR-16 measurements are unaffected by atmospheric conditions. To measure spectral information on rice paddies, we observed three different points within two sites under clear sky conditions (red points in Kyehwa and Kimjae, Figure1) and then calculated the average of three points at each site (Figure 2). The instrument was operated during the growing season (13 June 2014, 26 June 2014, 15 July 2014, 22 July 2014, 11 August 2014, 26 August 2014, 15 September 2014, and 3 October 2014). Field observations were compared with satellite observation pixels to geometrically match observation points (500 × 500 m).Figure 2
Detailed field-observation site images acquired on 5 October 2014, using the Landsat-8 operational land imager. The red dots are field-observation sites (a and b from Figure1), and the rectangles represent sites corresponding to one GOCI pixel.
(a)
(b)


## 2.3. Preprocessing of Satellite Data and BRDF Modeling

GOCI imagery was preprocessed before BRDF modeling. The preprocessing was included in converting the digital number to top-of-atmosphere (TOA) reflectance, cloud masking, and atmospheric correction. Cloud masking used a threshold method to reduce the probability of cloudy contamination [27]. An atmospheric correction model was applied to the TOA reflectance to estimate top-of-canopy reflectance using a look-up table based on the Second Simulation of the Satellite Signal in the Solar Spectrum [28–31].To correct for surface anisotropy effects, we used the BRDF model developed by Roujean et al. [32]; in the semiempirical BRDF model proposed, there are three scattering components: isotropic scattering, geometric scattering, and volumetric scattering from a horizontally homogeneous vegetation canopy [33]. Consider the following:(1)
ρ
θ
s
,
θ
v
,
∅
=
k
0
+
k
1
f
1
θ
s
,
θ
v
,
∅
+
k
2
f
2
θ
s
,
θ
v
,
∅
,where ρ is the observed reflectance in a given spectral channel, θ
s is the solar zenith angle (SZA), θ
v is the viewing zenith angle (VZA), ∅ is the relative azimuth angle (RAA), k
0 is isotropic coefficient of the Lambertian reflectance when θ
s
=
0 andθ
v
=
0, k
1 and k
2 are the empirical coefficients of the geometric and volumetric kernel, and f
1 and f
2 are physical geometric and volumetric scattering kernels depending on θ
s, θ
v, and ∅ between the sun and view direction of view. We applied our own fixed viewing angle and mean SZA during the composition period, without adjusting the VZA to the nadir direction based on Yeom and Kim [20].Inadequate observation sampling leads to poor surface coverage, while long composite periods reduce BRDF model sensitivity to real-time vegetation change, even with adequate cloud-free observations. To improve the vegetation profiles estimated with BRDF composite methods, several composite periods were added to enhance the accuracy and sensitivity of BRDF modeling relative to ground measurements. Here, we used 4-, 8-, 12-, and 16-day composite periods for each pixel to determine the optimized GOCI-based, BRDF-adjusted NDVI for each study area. Figure3 shows the concept of alterative composite periods for BRDF modeling. To perform kernel-driven, semiempirical BRDF modeling, at least seven observations of clear free surface reflectance are required for reliable simulation results [15, 34]. Finally, we estimated the daily NDVI product-based BRDF-adjusted surface reflectance (BAR).Figure 3
Concept of composite period for BRDF modeling.


## 3. Result

Before determining the optimized composite period of BRDF modeling, the accuracy of each composite period was assessed by simulating surface reflectance using BRDF coefficient parameters at the same sensor-target-sun geometry observations. Figure4 shows the correlation between the measured channel surface reflectance of GOCI and the modeled reflectance from each of the BRDF composite periods. The results of modeled surface reflectance are shown in Figures 4(a)–4(d). Bands 5 and 8 have the lowest root mean square error (RMSE) values in all composite periods, which indicates that there was no difference in channel reflectance among the composite periods. The measured reflectance and modeled channel reflectance were not affected by composite period under clear sky conditions. We inferred that similar accuracy would be obtained for each BRDF composite period by including at least seven cloud-free angular sampling for all BRDF performances.Figure 4
Scatterplots show the modeled channel reflectance versus measured channel reflectance for a (a) 4-day, (b) 8-day, (c) 12-day, and (d) 16-day composite period.
(a)
(b)
(c)
(d)The GOCI BRDF-adjusted NDVI from each composite period was compared to assess the spatial coverage of all composite periods. Figure5 shows sample data of GOCI BRDF-adjusted NDVI from different composite periods. The 4-day composite period showed low coverage rate of usable surface reflectance pixels due to having a relatively low number of cloud-free observations (Figure 5(a)). The coverage increased with increasing composite period length due to the inclusion of more cloud-free measurements. The cloudy areas were nearly absent in Figures 5(c) and 5(d), due to adequate cloud-free sampling within the longer composite periods. In other words, longer composite days account for surface reflectance that would otherwise be obscured by cloud cover. However, longer composite days reduce sensitivity to changes in surface reflectance. Consequently, it is critical to identify an optimized BRDF composite period.Figure 5
GOCI normalized difference vegetation index- (NDVI-) based bidirectional reflectance distribution function- (BRDF-) adjusted surface reflectance, by composite period (7 September 2014): (a) 4-day, (b) 8-day, (c) 12-day, and (d) 16-day composite period. White areas correspond to the ocean and cloud cover.Figure6 shows a comparison of GOCI NDVI-based BAR and field-observation NDVI from crop scans for the two crop sites; the measured NDVI is higher than the GOCI NDVI. This discrepancy could be explained by atmospheric correction and land cover types. The issue of atmospheric correction based on satellite data can be difficult to resolve. The field-measured NDVI value is observed directly at the rice paddy area, while GOCI NDVI includes many types of land cover, such as crops, roads, and structures, despite careful selection of the study sites. GOCI NDVI describes the general vegetation phenology pattern based on satellite data; the measured NDVI and GOCI NDVI have similar vegetation growth curves. Figure 6(a) shows the NDVI time series-based BRDF-adjusted surface reflectance at the Kyehwa site with late-maturing crops; Figure 6(b) shows the data for the Kimjae site with early-maturing crops. The BRDF-adjusted NDVI profile for all composite periods is sequentially shifted to the right compared with the vegetation growth curve of field-measured NDVI in both areas. BAR NDVI with 12-day and 16-day composite periods better followed the crop dynamics than the 4-day and 8-day composite periods. Therefore, the calculated BAR NDVI on 12-day and 16-day composite periods shows a stable crop profile. It is difficult to compare BAR NDVI and field-observation NDVI one-to-one, because heavy rain between Julian days 200 and 210 in June 2014 caused a sudden drop in BAR NDVI values. However, crop blooming during the field-observation periods appeared as maximum NDVI more quickly in the 12-day composite period than in the 16-day composite period. BAR NDVI in the 12-day composite period was sensitive to real-time change in comparison with field-observation NDVI and provided reliable temporal crop dynamics.Figure 6
NDVI time series-based BRDF-adjusted surface reflectance at (a) the Kyehwa site and (b) the Kimjae site, by composite periods. The black, grey, red, and blue dots represent 4-day, 8-day, 12-day, and 16-day composite periods, respectively, and the solid line is the field-observation NDVI.
(a)
(b)


## 4. Conclusion

This study estimated vegetation profiles using an optimized BRDF composite method. We ran our simulations using various composite periods and compared the NDVI profile derived from GOCI using four composite periods with field-observation NDVI. The BRDF composite method was tested with regard to modeling accuracy, spatial coverage, and sensitivity to real-time vegetation profiles. Finally, we determined that the 12-day composite period was both more sensitive to real-time vegetation changes and similar in terms of spatial coverage to the 16-day composite period. Although 4-day and 8-day composites were more sensitive to changes in ground measurements than the 12-day composite, they produced discontinuous vegetation profiles.The optimized BRDF composite period identified in this study will improve the GOCI satellite’s ability to measure real-time terrestrial products using the BRDF composite method.



---

*Source: 7165326-2016-02-23.xml*