# The Potential Scattering Model for Oil Palm Phenology Based on Spaceborne X-, C-, and L-Band Polarimetric SAR Imaging

Authors: Soni Darmawan; Ita Carolita; Rika Hernawati; Dede Dirgahayu; Agustan; Didin Agustian Permadi; Dewi Kania Sari; Widya Suryadini; Dhimas Wiratmoko; Yohanes Kunto
Journal: Journal of Sensors
 (2021)

---

## Abstract

Information about oil palm phenology is required for oil palm plantation management, but using spaceborne polarimetric radar imagery remains challenging. However, spaceborne polarimetric radar on X-, C-, and L-band is promising on structure vegetation and cloud area. This study investigates the scattering model of oil palm phenology based on spaceborne X-, C-, and L-band polarimetric Synthetic Aperture Radar (SAR) imaging. The X-, C-, and L-band polarimetric SAR are derived from spaceborne of TerraSAR-X, Sentinel-1A, and ALOS PALSAR 2. Study area is located in oil palm plantations, Asahan District, North Sumatra, Indonesia. The methodology includes data collection, preprocessing, radiometric calibration, speckle filtering, terrain correction, extraction of scattering value, and development of scattering model of oil palm phenology. The results showed different scattering characteristics for the X-, C-, and L-band polarimetric SAR of oil palm for age and found the potential of the scattering model for oil palm phenology based on the X-band on HH polarization that showed a nonlinear model withR2=0.65. The C-band on VH and VV polarization showed a nonlinear model with R2=0.56 and R2=0.89. The L-band on HV and HH polarization showed a logarithmic model with R2=0.50 and R2=0.51. In this case, the most potential of the scattering model of oil palm phenology based onR2is using C-band on VV polarization. However, the scattering model based on X-, C-, and L-band is potentially to be used and applied to identify the phenology of oil palm in Indonesia, which is the main parameter in yield estimation. For the future phenology model needs to improve accuracy by integrating multisensors, including different wavelengths on optical and microwave sensors and more in situ data.


---

## Body

## 1. Introduction

Oil palm (Elaeis guineensis) is a palm tree that is generally planted in Southeast Asia, particularly in Malaysia, Indonesia, and Thailand. Oil palm can produce cooking oil, mechanical oil, and fuel, widely utilized in daily life and industry [1, 2]. Oil palm obtained from seeds or portions in a hard mesocarp shell produces around 80% saturated fat that can be utilized as a raw material for the production of soap, cleansers, and other substances in the oleochemical industry [3–5]. Oil palm has an exceptional potential to be used as a biofuel in the future [6]. High yields and low production costs from oil palm are the reasons that encourage commercial plantation companies to develop oil palm trees for a considerable scope [7]. Because of expanding worldwide interest for nourishment and fuel, oil palm cultivation has been extended exponentially [8], and nowadays, oil palm is one of the most consumed vegetable oils in the world [9].Southeast Asia has the ideal condition for oil palm cultivation because it needs the moist equatorial condition [10]. Oil palm is adapting to the tropical atmosphere with a high precipitation rate, high intensity of global radiation, and warm temperature of 24–32°C [11]. The regular dry season in the tropical areas could significantly reduce oil palm yields [3]. The oil palms are generally planted in triangular patterns with a nine-meter interplanting distance (see Figure 1), following an industry standard to maximize yield with optimal sunlight penetration [3]. A planting density of 130–140 palms for each is the general practice. However, it varies according to the planting conditions and oil palm breed type [9, 12, 13]. Terrace planting usually facilitates a countermeasure to run off water and maximizing planting density in hilly areas [14].Figure 1
The oil palms planted in triangular patterns (source: [9, 15]).Oil palm plantations are planted per block so that all plants in each block have a uniform age. The age of oil palm is the main parameter in yield estimation. Oil palm age is one of the significant variables affecting the production of fruit bunches [3, 10]. Age information of oil palm is a good indicator of yield prediction as it influences the quality and quantity of the fresh fruit bunches [16]. Besides, age information is essential to precision agriculture to recognize anomalies among oil palms inside a specific age group to plan counteractive management practices and optimize management resources [10, 17]. In other cases, oil palm age information is required for the tax estimation, replanting time, and identification of oil palm diseases [15]. Organization or plantation management favors such useful data for maximizing oil palm yield, one of the most important influencing benefits [10]. According to Tan et al. [10], collecting age information on oil palm trees is time-consuming and costly, particularly on a large or regional scale. The age of oil palm growth affects the physical and environment of the oil palm plantation itself. According to Tan et al. [10] and Chemura et al. [16], canopy height and size of oil palms have a strong correlation with age.Spaceborne technology has been commonly applied in the agriculture and forestry sectors. These technologies are proven to provide a precise, affordable, and efficient solution for agricultural and forestry planning, monitoring, and management [18]. According to Henderson and Lewis [19], although sensors in the optical range of the electromagnetic spectrum have gotten the best consideration and have widely used, considerable effort has been invested into the utilization of radar sensors. Microwave sensors have become a promising technology in the application of remote sensing due to their cloud penetration ability and their capability to obtain data in all-weather day-night conditions. The microwave sensor also relies on its internal energy source, unlike the optical sensors, which rely on sunlight’s external energy source. Scatter radar is very sensitive to the dielectric properties (soil moisture and vegetation) and attributes of the geometric conditions (surface roughness) of objects on the surface of Earth [20].Many regions in the world (e.g., areas covered by clouds and lacking light), radar is the only sensor that can reliably provide consistent and periodic data. A radar sensor can obtain information in the electromagnetic spectrum bands K (1.1–1.7 cm), X (2.4–3.8 cm), C (3.8–7.5 cm), L (15–30 cm), and P (30–100 cm) [20], with the polarization of Horizontal to Horizontal (HH), Vertical to Vertical (VV), Horizontal to Vertical (HV), or Vertical to Horizontal (VH), which have varying ranges and azimuth resolutions [21]. Each of these wavelengths has a unique characteristic related to the reflection from forest stands [22]. The X-band interacts with leaves and canopy cover surfaces; hence, it is very suited for information on the tree canopy surface layer [23]. The C-band can penetrate through the leaves and spread to small branches and other underlying objects [24]. The L-band, which has a higher penetration capability, can penetrate the surface layer and spread to the stems and main branches [25, 26]. The P-band has the highest penetration capability and can penetrate the canopy covers [27].Many studies developed the model of oil palm phenology derived from scattering characteristics and any parameter, such as derived from biomass [27], leaf area index (LAI) [10, 17], height tree and tree diameter [10], crown projection area (CPA) [16], scattering [28], NDVI [10, 17, 29–31], and spectral band [10, 17, 28, 32]. However, to our best knowledge, information about phenology based on spaceborne X-, C-, and L-band polarimetric SAR has been limited and fragmented in Indonesia. In this study, we pioneered to investigate the scattering model of oil palm phenology in Indonesia based on spaceborne X-, C-, and L-band polarimetric SAR simultaneously.


## 2. Materials and Methods

### 2.1. Study Area

Study area in Asahan Regency, North Sumatra, Indonesia. This regency has the capital of Kisaran City and covers an area of 3702.97 km2. A part of Asahan regency was selected as the study area, which is an oil palm plantation area located in Hessa Plantation, Simpang Empat, Asahan Regency, 2.98° east longitude and 99.67° north latitude and 2.92° east longitude and 99.75° north latitude (Figure 2).Figure 2
Study Area in Asahan Regency, North Sumatra, Indonesia.The Asahan Regency is in the central part of the North Sumatra Province eastern coast. The district is mainly agricultural, and the primary land uses are oil palm and rubber plantations inland, with coconut groves and aquaculture ponds by the sea [33]. The first Indonesian oil palm plantation was founded in Asahan [34]. In 1911 during the Dutch colonial era, a Belgian company opened the first commercial oil palm plantation [35], and Asahan was a pioneer in the East Coast area of Sumatra [34, 36]. The large plantation area owned by many private companies and the state-owned Perkebunan Nusantara Company contributes a great deal to most Asahan citizens’ economic needs. The plasma plantation scheme is owned by a smallholder whose management is supported by the company (as a central plantation) [34].Asahan District has a factory capacity of 50 tons of fresh fruit bunch (FFB) per hour [37]. According to Budidarsono et al. [38], on the eastern coast of Sumatra, oil palm production (CPO) increased drastically. The area of the first farm, which was built in 1910–1914, was 2620 ha. The location biophysics is suitable for oil palm growing, with high rainfall (minimum 1600 mm/year) and a tropical climate within 10° of the equator. Land and labor, the most significant inputs, were available [38].


### 2.2. Methodology

The methodology includes data collection, preprocessing, radiometric calibration, speckle filtering, terrain correction, extracting scattering value, and development scattering model based on spaceborne X-, C-, and L-band polarimetric SAR (Figure3).Figure 3
Methodology.


#### 2.2.1. Data Collection

Collecting of variety polarimetric SAR imaging, including TerraSAR-X (X-band), Sentinel-1A (C-band), and ALOS PALSAR 2 (L-band) (Figure4). TerraSAR-X and ALOS PALSAR-2 data were obtained from the National Institute of Aeronautics and Space-Indonesia (LAPAN), while Sentinel-1A data were collected from the European Space Agency (ESA) Copernicus. The details of the SAR imaging are presented in Table 1.Figure 4
SAR imaging on TerraSAR-X on HH (a), Sentinel-1 on VV and VH (b), and ALOS PALSAR on HH and HV (c).
(a)(b)(c)Table 1
The characteristics of SAR data.
SensorTerraSAR-XSentinel-1AALOS PALSAR-2WavelengthX-bandλ=3.1 cmC-bandλ=5.6 cmL-bandλ=23.5 cmFrequency8-12 GHz4-8 GHz1-2 GHzPolarizationSingle HHDual VV/VHDual HH/HVResolutionStripmap:3×3 mInterferometric Wide Swath (IW):5×20 mStripmap: 3-10 mFrame sizeStripmap:50×30 kmIW: 250 kmStripmap:55×70 - 70×70 kmTemporal11 days12 days14 daysDate of acquisition29 August 201715 July 201817 January 2015Incidence angle range28.7°39°34.9°The topography elevation data were derived from Shuttle Radar Topographic Mission (SRTM), which has a spatial resolution of around 90 meters, and were processed by the National Aeronautics and Space Administration (NASA) and the United States Geological Survey (USGS). The data were projected in a geographic (Lat/Lon) projection, with the WGS84 horizontal datum. Field data collection was conducted in 2020 that collected planting age blocks of oil palm and some photograph in the study area (Figures5 and 6). In this study area, we used 174 blocks, which represent 11 varieties of oil palm age from 0 to 25 years. This planting age block data is used as a boundary in taking the region of interest (ROI) based on the scattering value with totally around 44,000 pixels.Figure 5
The planting age blocks of oil palm in the study area.Figure 6
Oil palm age in study area.


#### 2.2.2. Preprocessing

(1) Speckle Filtering. SAR imaging interpretation and classification principal problem is the speckle effect caused by wave coherent interference reflected from many primary scatterers [39]. Speckle appears in SAR imaging as granular noise due to waves interference reflected from many basic scatterers [40]. Speckle filtering is a procedure to increase image quality by reducing speckle. The enhanced Lee filter with a window size of 7×7 has been chosen in this study to minimize speckle effect on the images [41, 42], also to remove noise in the edge and high contrast area [40] without loss of information [43], because previous research showed reduction in speckles on SAR data using the Lee Filter before extracting the scattering value [24]. According to Lee et al., this algorithm is the best algorithm for segmentation crop filed purpose [40].(2) Terrain Correction. The preprocessed image was also terrain correction to eliminate the distortion introduced by the topographical variations. The purpose of terrain correction for SAR imaging is to minimize the SAR geometry effects (foreshortening, layover, and shadow) towards radar images [24]. The Shuttle Radar Topography Mission (SRTM) was used as the Digital Elevation Model (DEM) to provide height information [44] using Range-Doppler Terrain Correction [24]. The geometric corrections were performed on all SAR (X-, C-, and L-band) data of the study area. All images were rectified, and the data were projected in geographic (Lat/Long) projection, with the WGS84 datum.(3) Radiometric Calibration. Calibration of the scattering values is necessary for the intercomparison of radar images acquired with different sensors or images obtained by the same sensor if acquired in different modes or processed with different processors [45]. Radiometric calibration resulting in sigma naught (σ°) [45].The radiometric calibration process for X-band first calculated the radar brightnessβ°dB (beta zero, as contained in dB) derived from the polarization-specific calibration factor kS (in equation (1)) and converted from digital numbers (DN) to sigma naught (σ°) using equation (2) [39, 46]:
(1)β°dB=10×log10KSDN2.The final radiometric calibration was performed by calculating sigma naught (σ°) using
(2)σ°=β°dB+log10sinθloc.where θloc, the local incident angle image, is defined for the four scenes’ angle position, based on 2D interpolation across the performed image dimensions. The interpolation method applied is based on Delaunay’s triangulation of data using the Qhull algorithm [47]. In conclusion, the local incident image has the same dimensions as the scene image and produces specific information about each pixel’s local incident angle [39].While the scattering value of the C-band is transforming the gamma-calibrated scattering coefficient [48], the C-band product uses a radiometric calibration look-up table (LUT) to do the calibration [49]. The essential conversion of amplitude to DN and from DN to sigma naught was done automatically on SNAP, and once the sigma naught values were obtained [41, 50], the computation of scattering (σ°dB) can be performed in
(3)σ°dB=10·log10γi,where γiis the gamma-calibrated scattering coefficient of the C-band [48].The L-band was built on a 16-bit data type, and all pixels have digital numbers (DN) ranging from 0 to 65,535 [41]. The DN has to be converted to scattering (i.e., the reflected radar signals) recognized as Normalized Radar Cross Section (NRCS) and expressed in σ° in decibels (dB) [41]. The conversion of HH (DNHH) and HV (DNHV) scattering intensities into NRCS (that is σ°HH and σ°HV) [28, 51] was based on Shimada’s study [51] shown in
(4)σ°HHdB=10×log10DNHH2−CF,(5)σ°HVdB=10×log10DNHV2−CF,where σ° is the scattering coefficient and CF is the calibration factor. The CF is dependent on the processing date. In this study, CF is equal to -83.0 both for HH and HV [28].


#### 2.2.3. Extracting Scattering Characteristic

We created a region of interest (ROI) for area sampling; also the ROI was calculated to determine the scattering characteristics of X-, C-, and L-band polarimetric SAR imaging with the age of the oil palm in the forms of mean and standard deviation values using the algorithm as follows (equation (6) and (7)) [52]:
(6)x¯=∑j=1nxjn,(7)s2=∑j=1nxi−x¯2n−1,where n is the number of samples taken, xi is the sample value, x¯ is the average of the samples, and s is the standard deviation.


#### 2.2.4. The Development of Scattering Model

We analyzed the relationship between the scattering value of SAR imaging polarization (X-band on HH, C-band on VV and VH, and L-band on HH and HV) and the age of oil palm plantation. The scattering response from oil palm areas depends on the structure and density of the oil palms. All oil palm ranges, and even under those circumstances, the various structures, and densities have different scattering patterns and textures at various wavelengths [28]. Analyzing the relationship between scattering and age is using the Pearson correlation coefficient (equation (8)). Pearson’s method assumed that if the coefficient is similar to one, it implies a strong correlation between the two variables [53].
(8)R2=r=∑Xi−X¯Yi−Y¯∑Xi−X¯2∑Yi−Y¯21/2,where X¯ is the mean of the X variable and Y¯ is the mean of the Y variable. In this study, the first variable is the oil palm age, and the second variable is scattering polarization value. The coefficient of determination (R2) is an indication of the regression model [54]. The R2 value will show the percentage of variation as a regression model variable [55]. Several methods have been proposed to interpret the correlation coefficient into descriptors like “weak,” “moderate,” or “strong” relationship [56].


## 2.1. Study Area

Study area in Asahan Regency, North Sumatra, Indonesia. This regency has the capital of Kisaran City and covers an area of 3702.97 km2. A part of Asahan regency was selected as the study area, which is an oil palm plantation area located in Hessa Plantation, Simpang Empat, Asahan Regency, 2.98° east longitude and 99.67° north latitude and 2.92° east longitude and 99.75° north latitude (Figure 2).Figure 2
Study Area in Asahan Regency, North Sumatra, Indonesia.The Asahan Regency is in the central part of the North Sumatra Province eastern coast. The district is mainly agricultural, and the primary land uses are oil palm and rubber plantations inland, with coconut groves and aquaculture ponds by the sea [33]. The first Indonesian oil palm plantation was founded in Asahan [34]. In 1911 during the Dutch colonial era, a Belgian company opened the first commercial oil palm plantation [35], and Asahan was a pioneer in the East Coast area of Sumatra [34, 36]. The large plantation area owned by many private companies and the state-owned Perkebunan Nusantara Company contributes a great deal to most Asahan citizens’ economic needs. The plasma plantation scheme is owned by a smallholder whose management is supported by the company (as a central plantation) [34].Asahan District has a factory capacity of 50 tons of fresh fruit bunch (FFB) per hour [37]. According to Budidarsono et al. [38], on the eastern coast of Sumatra, oil palm production (CPO) increased drastically. The area of the first farm, which was built in 1910–1914, was 2620 ha. The location biophysics is suitable for oil palm growing, with high rainfall (minimum 1600 mm/year) and a tropical climate within 10° of the equator. Land and labor, the most significant inputs, were available [38].


## 2.2. Methodology

The methodology includes data collection, preprocessing, radiometric calibration, speckle filtering, terrain correction, extracting scattering value, and development scattering model based on spaceborne X-, C-, and L-band polarimetric SAR (Figure3).Figure 3
Methodology.


### 2.2.1. Data Collection

Collecting of variety polarimetric SAR imaging, including TerraSAR-X (X-band), Sentinel-1A (C-band), and ALOS PALSAR 2 (L-band) (Figure4). TerraSAR-X and ALOS PALSAR-2 data were obtained from the National Institute of Aeronautics and Space-Indonesia (LAPAN), while Sentinel-1A data were collected from the European Space Agency (ESA) Copernicus. The details of the SAR imaging are presented in Table 1.Figure 4
SAR imaging on TerraSAR-X on HH (a), Sentinel-1 on VV and VH (b), and ALOS PALSAR on HH and HV (c).
(a)(b)(c)Table 1
The characteristics of SAR data.
SensorTerraSAR-XSentinel-1AALOS PALSAR-2WavelengthX-bandλ=3.1 cmC-bandλ=5.6 cmL-bandλ=23.5 cmFrequency8-12 GHz4-8 GHz1-2 GHzPolarizationSingle HHDual VV/VHDual HH/HVResolutionStripmap:3×3 mInterferometric Wide Swath (IW):5×20 mStripmap: 3-10 mFrame sizeStripmap:50×30 kmIW: 250 kmStripmap:55×70 - 70×70 kmTemporal11 days12 days14 daysDate of acquisition29 August 201715 July 201817 January 2015Incidence angle range28.7°39°34.9°The topography elevation data were derived from Shuttle Radar Topographic Mission (SRTM), which has a spatial resolution of around 90 meters, and were processed by the National Aeronautics and Space Administration (NASA) and the United States Geological Survey (USGS). The data were projected in a geographic (Lat/Lon) projection, with the WGS84 horizontal datum. Field data collection was conducted in 2020 that collected planting age blocks of oil palm and some photograph in the study area (Figures5 and 6). In this study area, we used 174 blocks, which represent 11 varieties of oil palm age from 0 to 25 years. This planting age block data is used as a boundary in taking the region of interest (ROI) based on the scattering value with totally around 44,000 pixels.Figure 5
The planting age blocks of oil palm in the study area.Figure 6
Oil palm age in study area.


### 2.2.2. Preprocessing

(1) Speckle Filtering. SAR imaging interpretation and classification principal problem is the speckle effect caused by wave coherent interference reflected from many primary scatterers [39]. Speckle appears in SAR imaging as granular noise due to waves interference reflected from many basic scatterers [40]. Speckle filtering is a procedure to increase image quality by reducing speckle. The enhanced Lee filter with a window size of 7×7 has been chosen in this study to minimize speckle effect on the images [41, 42], also to remove noise in the edge and high contrast area [40] without loss of information [43], because previous research showed reduction in speckles on SAR data using the Lee Filter before extracting the scattering value [24]. According to Lee et al., this algorithm is the best algorithm for segmentation crop filed purpose [40].(2) Terrain Correction. The preprocessed image was also terrain correction to eliminate the distortion introduced by the topographical variations. The purpose of terrain correction for SAR imaging is to minimize the SAR geometry effects (foreshortening, layover, and shadow) towards radar images [24]. The Shuttle Radar Topography Mission (SRTM) was used as the Digital Elevation Model (DEM) to provide height information [44] using Range-Doppler Terrain Correction [24]. The geometric corrections were performed on all SAR (X-, C-, and L-band) data of the study area. All images were rectified, and the data were projected in geographic (Lat/Long) projection, with the WGS84 datum.(3) Radiometric Calibration. Calibration of the scattering values is necessary for the intercomparison of radar images acquired with different sensors or images obtained by the same sensor if acquired in different modes or processed with different processors [45]. Radiometric calibration resulting in sigma naught (σ°) [45].The radiometric calibration process for X-band first calculated the radar brightnessβ°dB (beta zero, as contained in dB) derived from the polarization-specific calibration factor kS (in equation (1)) and converted from digital numbers (DN) to sigma naught (σ°) using equation (2) [39, 46]:
(1)β°dB=10×log10KSDN2.The final radiometric calibration was performed by calculating sigma naught (σ°) using
(2)σ°=β°dB+log10sinθloc.where θloc, the local incident angle image, is defined for the four scenes’ angle position, based on 2D interpolation across the performed image dimensions. The interpolation method applied is based on Delaunay’s triangulation of data using the Qhull algorithm [47]. In conclusion, the local incident image has the same dimensions as the scene image and produces specific information about each pixel’s local incident angle [39].While the scattering value of the C-band is transforming the gamma-calibrated scattering coefficient [48], the C-band product uses a radiometric calibration look-up table (LUT) to do the calibration [49]. The essential conversion of amplitude to DN and from DN to sigma naught was done automatically on SNAP, and once the sigma naught values were obtained [41, 50], the computation of scattering (σ°dB) can be performed in
(3)σ°dB=10·log10γi,where γiis the gamma-calibrated scattering coefficient of the C-band [48].The L-band was built on a 16-bit data type, and all pixels have digital numbers (DN) ranging from 0 to 65,535 [41]. The DN has to be converted to scattering (i.e., the reflected radar signals) recognized as Normalized Radar Cross Section (NRCS) and expressed in σ° in decibels (dB) [41]. The conversion of HH (DNHH) and HV (DNHV) scattering intensities into NRCS (that is σ°HH and σ°HV) [28, 51] was based on Shimada’s study [51] shown in
(4)σ°HHdB=10×log10DNHH2−CF,(5)σ°HVdB=10×log10DNHV2−CF,where σ° is the scattering coefficient and CF is the calibration factor. The CF is dependent on the processing date. In this study, CF is equal to -83.0 both for HH and HV [28].


### 2.2.3. Extracting Scattering Characteristic

We created a region of interest (ROI) for area sampling; also the ROI was calculated to determine the scattering characteristics of X-, C-, and L-band polarimetric SAR imaging with the age of the oil palm in the forms of mean and standard deviation values using the algorithm as follows (equation (6) and (7)) [52]:
(6)x¯=∑j=1nxjn,(7)s2=∑j=1nxi−x¯2n−1,where n is the number of samples taken, xi is the sample value, x¯ is the average of the samples, and s is the standard deviation.


### 2.2.4. The Development of Scattering Model

We analyzed the relationship between the scattering value of SAR imaging polarization (X-band on HH, C-band on VV and VH, and L-band on HH and HV) and the age of oil palm plantation. The scattering response from oil palm areas depends on the structure and density of the oil palms. All oil palm ranges, and even under those circumstances, the various structures, and densities have different scattering patterns and textures at various wavelengths [28]. Analyzing the relationship between scattering and age is using the Pearson correlation coefficient (equation (8)). Pearson’s method assumed that if the coefficient is similar to one, it implies a strong correlation between the two variables [53].
(8)R2=r=∑Xi−X¯Yi−Y¯∑Xi−X¯2∑Yi−Y¯21/2,where X¯ is the mean of the X variable and Y¯ is the mean of the Y variable. In this study, the first variable is the oil palm age, and the second variable is scattering polarization value. The coefficient of determination (R2) is an indication of the regression model [54]. The R2 value will show the percentage of variation as a regression model variable [55]. Several methods have been proposed to interpret the correlation coefficient into descriptors like “weak,” “moderate,” or “strong” relationship [56].


## 2.2.1. Data Collection

Collecting of variety polarimetric SAR imaging, including TerraSAR-X (X-band), Sentinel-1A (C-band), and ALOS PALSAR 2 (L-band) (Figure4). TerraSAR-X and ALOS PALSAR-2 data were obtained from the National Institute of Aeronautics and Space-Indonesia (LAPAN), while Sentinel-1A data were collected from the European Space Agency (ESA) Copernicus. The details of the SAR imaging are presented in Table 1.Figure 4
SAR imaging on TerraSAR-X on HH (a), Sentinel-1 on VV and VH (b), and ALOS PALSAR on HH and HV (c).
(a)(b)(c)Table 1
The characteristics of SAR data.
SensorTerraSAR-XSentinel-1AALOS PALSAR-2WavelengthX-bandλ=3.1 cmC-bandλ=5.6 cmL-bandλ=23.5 cmFrequency8-12 GHz4-8 GHz1-2 GHzPolarizationSingle HHDual VV/VHDual HH/HVResolutionStripmap:3×3 mInterferometric Wide Swath (IW):5×20 mStripmap: 3-10 mFrame sizeStripmap:50×30 kmIW: 250 kmStripmap:55×70 - 70×70 kmTemporal11 days12 days14 daysDate of acquisition29 August 201715 July 201817 January 2015Incidence angle range28.7°39°34.9°The topography elevation data were derived from Shuttle Radar Topographic Mission (SRTM), which has a spatial resolution of around 90 meters, and were processed by the National Aeronautics and Space Administration (NASA) and the United States Geological Survey (USGS). The data were projected in a geographic (Lat/Lon) projection, with the WGS84 horizontal datum. Field data collection was conducted in 2020 that collected planting age blocks of oil palm and some photograph in the study area (Figures5 and 6). In this study area, we used 174 blocks, which represent 11 varieties of oil palm age from 0 to 25 years. This planting age block data is used as a boundary in taking the region of interest (ROI) based on the scattering value with totally around 44,000 pixels.Figure 5
The planting age blocks of oil palm in the study area.Figure 6
Oil palm age in study area.


## 2.2.2. Preprocessing

(1) Speckle Filtering. SAR imaging interpretation and classification principal problem is the speckle effect caused by wave coherent interference reflected from many primary scatterers [39]. Speckle appears in SAR imaging as granular noise due to waves interference reflected from many basic scatterers [40]. Speckle filtering is a procedure to increase image quality by reducing speckle. The enhanced Lee filter with a window size of 7×7 has been chosen in this study to minimize speckle effect on the images [41, 42], also to remove noise in the edge and high contrast area [40] without loss of information [43], because previous research showed reduction in speckles on SAR data using the Lee Filter before extracting the scattering value [24]. According to Lee et al., this algorithm is the best algorithm for segmentation crop filed purpose [40].(2) Terrain Correction. The preprocessed image was also terrain correction to eliminate the distortion introduced by the topographical variations. The purpose of terrain correction for SAR imaging is to minimize the SAR geometry effects (foreshortening, layover, and shadow) towards radar images [24]. The Shuttle Radar Topography Mission (SRTM) was used as the Digital Elevation Model (DEM) to provide height information [44] using Range-Doppler Terrain Correction [24]. The geometric corrections were performed on all SAR (X-, C-, and L-band) data of the study area. All images were rectified, and the data were projected in geographic (Lat/Long) projection, with the WGS84 datum.(3) Radiometric Calibration. Calibration of the scattering values is necessary for the intercomparison of radar images acquired with different sensors or images obtained by the same sensor if acquired in different modes or processed with different processors [45]. Radiometric calibration resulting in sigma naught (σ°) [45].The radiometric calibration process for X-band first calculated the radar brightnessβ°dB (beta zero, as contained in dB) derived from the polarization-specific calibration factor kS (in equation (1)) and converted from digital numbers (DN) to sigma naught (σ°) using equation (2) [39, 46]:
(1)β°dB=10×log10KSDN2.The final radiometric calibration was performed by calculating sigma naught (σ°) using
(2)σ°=β°dB+log10sinθloc.where θloc, the local incident angle image, is defined for the four scenes’ angle position, based on 2D interpolation across the performed image dimensions. The interpolation method applied is based on Delaunay’s triangulation of data using the Qhull algorithm [47]. In conclusion, the local incident image has the same dimensions as the scene image and produces specific information about each pixel’s local incident angle [39].While the scattering value of the C-band is transforming the gamma-calibrated scattering coefficient [48], the C-band product uses a radiometric calibration look-up table (LUT) to do the calibration [49]. The essential conversion of amplitude to DN and from DN to sigma naught was done automatically on SNAP, and once the sigma naught values were obtained [41, 50], the computation of scattering (σ°dB) can be performed in
(3)σ°dB=10·log10γi,where γiis the gamma-calibrated scattering coefficient of the C-band [48].The L-band was built on a 16-bit data type, and all pixels have digital numbers (DN) ranging from 0 to 65,535 [41]. The DN has to be converted to scattering (i.e., the reflected radar signals) recognized as Normalized Radar Cross Section (NRCS) and expressed in σ° in decibels (dB) [41]. The conversion of HH (DNHH) and HV (DNHV) scattering intensities into NRCS (that is σ°HH and σ°HV) [28, 51] was based on Shimada’s study [51] shown in
(4)σ°HHdB=10×log10DNHH2−CF,(5)σ°HVdB=10×log10DNHV2−CF,where σ° is the scattering coefficient and CF is the calibration factor. The CF is dependent on the processing date. In this study, CF is equal to -83.0 both for HH and HV [28].


## 2.2.3. Extracting Scattering Characteristic

We created a region of interest (ROI) for area sampling; also the ROI was calculated to determine the scattering characteristics of X-, C-, and L-band polarimetric SAR imaging with the age of the oil palm in the forms of mean and standard deviation values using the algorithm as follows (equation (6) and (7)) [52]:
(6)x¯=∑j=1nxjn,(7)s2=∑j=1nxi−x¯2n−1,where n is the number of samples taken, xi is the sample value, x¯ is the average of the samples, and s is the standard deviation.


## 2.2.4. The Development of Scattering Model

We analyzed the relationship between the scattering value of SAR imaging polarization (X-band on HH, C-band on VV and VH, and L-band on HH and HV) and the age of oil palm plantation. The scattering response from oil palm areas depends on the structure and density of the oil palms. All oil palm ranges, and even under those circumstances, the various structures, and densities have different scattering patterns and textures at various wavelengths [28]. Analyzing the relationship between scattering and age is using the Pearson correlation coefficient (equation (8)). Pearson’s method assumed that if the coefficient is similar to one, it implies a strong correlation between the two variables [53].
(8)R2=r=∑Xi−X¯Yi−Y¯∑Xi−X¯2∑Yi−Y¯21/2,where X¯ is the mean of the X variable and Y¯ is the mean of the Y variable. In this study, the first variable is the oil palm age, and the second variable is scattering polarization value. The coefficient of determination (R2) is an indication of the regression model [54]. The R2 value will show the percentage of variation as a regression model variable [55]. Several methods have been proposed to interpret the correlation coefficient into descriptors like “weak,” “moderate,” or “strong” relationship [56].


## 3. Results

The X-, C-, and L-band polarimetric SAR imaging was filtered using the Lee filter for each polarization. The scattering model of oil palm phenology was based on an empirical model derived from the relationship between oil palm age and scattering value of X-, C-, and L-band SAR imaging. In this study, ranges of oil palm age were divided for five years: 0-5, 5-10, 10-15, 15-20, and 20-25.


### 3.1. The Scattering Model of Oil Palm Phenology Based on X-Band

The scattering model of oil palm phenology was based on X-band with HH polarization. The value starts at the age of 0-5 years around -6.8 dB, of 5-10 years around -5.7 dB, of 10-15 years around -5 dB, and of 20-25 years around -6.4 dB (Figure7). On HH polarization, the scattering model has nonlinear regression on y=−0.0114x2+0.2741x–7.2823 with R2=0.65 (Figure 7).Figure 7
Scattering model of oil palm phenology on X-band with HH polarization.


### 3.2. The Scattering Model of Oil Palm Phenology Based on C-Band

The scattering model of oil palm plantation was based on C-band with VH and VV polarization. On VH polarization, the value starts in the age of 0-5 years around -14 dB, age of 5-10 years around -15 dB, age of 10-15 years around -14 dB, age of 15-20 years around -15 dB, and age of 20-25 years around -16 dB (Figure8(a)). On VH polarization, the scattering model has nonlinear regression y=−0.0055x2+0.1248x-15.785 with R2=0.56 (Figure 8(a)).Figure 8
Scattering model of oil palm phenology on C-band with VH polarization (a) and VV polarization (b).
(a)(b)On VV polarization, the value starts at the age of 0-5 years around -7.8 dB, at the age of 5-10 years around -7.4 dB, at the age of 10-15 years around -6.5 dB, at the age of 15-20 years around -7.4 dB, and at the age of 20-25 years around -8.0 dB (Figure8(b)). On VV polarization, the scattering model has nonlinear regression y=−0.0117x2+0.3006x-8.9534 with R2=0.89 (Figure 8(b)).


### 3.3. The Scattering Model of Oil Palm Phenology Based on L-Band

The scattering model of oil palm plantation based on L-band with HH and HV polarization. On HH polarization, the value starts at the age of 0-5 years around -12.6 dB, at the age of 5-10 years around -11 dB, at the age of 10-15 years around -11.3 dB, at the age of 15-20 years around -11.5 dB, and at the age of 20-25 years around -11.9 dB (Figure9(a)). The scattering model has logarithmic regression y=0.5293lnx–12.616 with R2=0.51 (Figure 9(a)).Figure 9
Scattering model of oil palm phenology on L-band with HH polarization (a) and HV polarization (b).
(a)(b)On HV polarization, the value starts at the age of 0-5 years around -22 dB, at the age of 5-10 years around -18 dB, at the age of 10-15 years around -19 dB, at the age of 15-20 years around -18.9 dB, and at the age of 20-25 years around -18.5 dB (Figure9(b)). The scattering model has logarithmic regression y=1.1555lnx-21.815 with R2=0.50 (Figure 9(b)).


## 3.1. The Scattering Model of Oil Palm Phenology Based on X-Band

The scattering model of oil palm phenology was based on X-band with HH polarization. The value starts at the age of 0-5 years around -6.8 dB, of 5-10 years around -5.7 dB, of 10-15 years around -5 dB, and of 20-25 years around -6.4 dB (Figure7). On HH polarization, the scattering model has nonlinear regression on y=−0.0114x2+0.2741x–7.2823 with R2=0.65 (Figure 7).Figure 7
Scattering model of oil palm phenology on X-band with HH polarization.


## 3.2. The Scattering Model of Oil Palm Phenology Based on C-Band

The scattering model of oil palm plantation was based on C-band with VH and VV polarization. On VH polarization, the value starts in the age of 0-5 years around -14 dB, age of 5-10 years around -15 dB, age of 10-15 years around -14 dB, age of 15-20 years around -15 dB, and age of 20-25 years around -16 dB (Figure8(a)). On VH polarization, the scattering model has nonlinear regression y=−0.0055x2+0.1248x-15.785 with R2=0.56 (Figure 8(a)).Figure 8
Scattering model of oil palm phenology on C-band with VH polarization (a) and VV polarization (b).
(a)(b)On VV polarization, the value starts at the age of 0-5 years around -7.8 dB, at the age of 5-10 years around -7.4 dB, at the age of 10-15 years around -6.5 dB, at the age of 15-20 years around -7.4 dB, and at the age of 20-25 years around -8.0 dB (Figure8(b)). On VV polarization, the scattering model has nonlinear regression y=−0.0117x2+0.3006x-8.9534 with R2=0.89 (Figure 8(b)).


## 3.3. The Scattering Model of Oil Palm Phenology Based on L-Band

The scattering model of oil palm plantation based on L-band with HH and HV polarization. On HH polarization, the value starts at the age of 0-5 years around -12.6 dB, at the age of 5-10 years around -11 dB, at the age of 10-15 years around -11.3 dB, at the age of 15-20 years around -11.5 dB, and at the age of 20-25 years around -11.9 dB (Figure9(a)). The scattering model has logarithmic regression y=0.5293lnx–12.616 with R2=0.51 (Figure 9(a)).Figure 9
Scattering model of oil palm phenology on L-band with HH polarization (a) and HV polarization (b).
(a)(b)On HV polarization, the value starts at the age of 0-5 years around -22 dB, at the age of 5-10 years around -18 dB, at the age of 10-15 years around -19 dB, at the age of 15-20 years around -18.9 dB, and at the age of 20-25 years around -18.5 dB (Figure9(b)). The scattering model has logarithmic regression y=1.1555lnx-21.815 with R2=0.50 (Figure 9(b)).


## 4. Discussion

Two significant characteristic categories determine the SAR scattering values are the sensor and the target characteristics [57]. The sensor category includes the frequency/wavelength of the SAR, the polarization of the transmitted and received SAR signal, the incident angle of the ground-interacting radar beam, and the sensor looking position [57]. According to Henderson and Lewis [19], the target characteristic influence varied with (1) vegetation type, (2) stand structure, and (3) canopy composition.


### 4.1. Characteristics Age of Oil Palm Based on SAR Imaging

The scattering characteristics of the forest vegetation are influenced by the thickness of the volume, density of the component particles (or scatters), size distribution of the component particles, shape distribution of the component particles, orientation distribution of the component particles, and dielectric properties of the component particles [58]. It is essential to know the polarization from which the SAR image obtained because signals at different polarizations interact differently with objects on the ground, which affects the brightness of the radar recorded in a particular polarization channel [59]. For simplicity, it is assumed that a natural scene can be described as a combination of three types of scatterers: (1) rough surface scatterers, (2) double-bounce scatterers, and (3) volume scatterers [60].


#### 4.1.1. On the X-Band

The scattering characteristics of the oil palm at the X-band of the young oil palm have the lowest scattering value as seen in Figure10 of the 4-year-old histogram which is caused by rough surface scatters. The mature oil palm has a higher scattering value as seen in Figure 10 on the 6 to 15-year-old histogram which is caused by double-bounce scattering. The old oil palm has a lower scattering value as seen in Figure 10 of the 16 to 23-year-old histogram which is caused by volume scatters.Figure 10
Histogram X-band for each age of oil palm.According to Martinis and Rieke [61], adopting the X-band, canopy attenuation, volume, and surface scattering from the top layer of the forest canopy is generally higher [62] because the X-band is sensitive to canopy surface scattering [63]. Hence, the same canopy surface seems rough for X-band wavelength causing high scattering values [64]. According to Rosenqvist [65], the X-band instrument clear relationship between canopy growth area and threshold ambiguity for trees.


#### 4.1.2. On the C-Band

The scattering characteristics of oil palm at the C-band on VV and VH polarization on the young oil palm have the lowest scattering value as seen in Figures11 and 12 on the histogram for 2 to 4-year-old histogram which is caused by rough surface scattering. The mature oil palm has a higher scattering value as seen in Figures 11 and 12 of the 7 to 14-year-old histogram which is caused by double-bounce scattering. The old oil palm has a lower scattering value as seen in Figures 11 and 12 of the 17 to 21-year-old histogram which is caused by volume scatters.Figure 11
Histogram C-band on VH for each age of oil palm.Figure 12
Histogram C-band on VV for each age of oil palm.According to Carolita et al. [15], increase and decrease in the scattering value occur because at the planting new age of trees height was relatively low and the size of the canopy is still small when the tree reached into productive age, the scattering value will rise because of the difference in the height of tree and soil. C-band data has a short wavelength, causing the scattering values for canopy density to become smaller. According to Teng et al. [6], there was also a significant contribution from the trunks to the surface volume scattering component, although this was not the dominant component. Also, at the shorter C-band wavelength, volume scattering can occur within the canopies of lower or sparse vegetation types, such as bushes, shrubs, or crops. According to Rignot et al. [66], radar signals were expected to be dominantly scattered by the foliage and the top branches and twigs of the canopy at the C-band.


#### 4.1.3. On the L-Band

The scattering characteristics of oil palm at the L-band on HH and HV polarization on the young oil palm have the lowest scattering value as seen in Figures13 and 14 on the 5-year-old histogram which is caused by in surface scattering. The mature oil palm has a higher scattering value as seen in Figures 13 and 14 of the 7 to 14-year-old histogram which is caused by double-bounce scattering. The old oil palm has a lower scattering value, as seen in Figures 13 and 14 of the 16 to 24-year-old histogram which is caused by volume scatters.Figure 13
Histogram L-band on HV for each age of oil palm.Figure 14
Histogram L-band on HH for each age of oil palm.According to Toh et al. [67], the scattering increased with age for all polarization in general for the L-band frequency. The L-band scattering coefficient showed the trend increasing with the age of oil palm more than 20 years. Besides that, according to Teng et al. [6], the scattering value for HH polarization using the L-band was increasing attenuation of waves propagating through to the soil surface and back, as the leaves, fronds, and trunks are growing larger.The scattering value for each polarization (HH and HV) among the age of oil palm was increasing; according to Darmawan et al. [28], it was hypothesized that the growth of oil palms was caused by higher trunks, more leaves, more branches, and larger canopy followed by the condition of scattering value. In any case, the scatter values derived from the L-band are depending on wavelengths, polarization, incident angle, and temporal data; environments such as land moisture and landscape [25, 28, 68]; and the structure of oil palms such as size, geometry, and orientation of leaves, trunks, branches, and aerials or stilt roots [28, 68, 69].


### 4.2. Coefficient Determination of Scattering Model of Oil Palm Phenology

The scatter model of oil palm phenology was performed separately for each radar wavelength (X-, C-, and L-band). In this study, the X-band showed the nonlinear oil palm phenology model in HH polarization with coefficient determinationR2=0.65, which means the model has moderate correlation. The C-band showed the nonlinear or quadratic model of oil palm phenology in VH polarization with coefficient determination R2=0.56, which means the model has moderate correlation, and while on VV polarization with coefficient determination R2=0.89, it means the model has strong correlation. The L-band showed the logarithmic model of oil palm phenology in HV polarization with coefficient determination R2=0.50, which means the model has moderate correlation, and while on HH polarization with coefficient determination R2=0.51, it means the model has moderate correlation. In this case, the more potential of scattering model of oil palm phenology based on R2 is using C-band on VV polarization.The scattering model with the moderate correlation also found by other researchers such as Tan et al. [10] has logarithmic regression for the correlating scattering value with oil palm age using ALOS PALSAR-2 with coefficient determination on HH polarizationR2=0.49, HV polarization withR2=0.27, and ratio of HH/HV polarization withR2=0.26, for study case in the southern part of peninsular Malaysia. Darmawan et al. [28] also used ALOS PALSAR-2 on HH and HV polarization having a logarithmic regression model with coefficient determination R2=0.62 and 0.41 for study case area in Jerantut, Pahang, Malaysia. In other cases, Avtar et al. [70] for study area near Miri City, Sarawak, Malaysia, showed a logarithmic regression model for a variety of SAR data including TerraSAR-X (X-band) on HH polarization with coefficient determination R2=0.075; Radarsat-2 (C-band) on HH, HV, and VV polarization with coefficient determination R2=0.39, 0.49, and 0.39; and ALOS PALSAR-2 (L-band) on HH and HV polarization with coefficient determination R2=0.62 and 0.77. Okarda et al. [32] studied the relationship between the age of oil palm in mineral soil and age and peatland in Central Kalimantan, Indonesia, using ALOS PALSAR-2 on HH and HV polarization. The scattering model has linear regression with coefficient determination on HH and HV between age and mineral soil with R2=0.55 and 0.37 and age and peatland on HH and HV with R2=0.36 and 0.28. For C-band data, it was studied by Carolita et al. [15] using Sentinel-1A on HH and HV polarization. The scattering model has nonlinear regression (quadratic model) with coefficient determination R2=0.68 for HH and HV with R2=0.77 for study case area in Asahan, North Sumatra Indonesia.


## 4.1. Characteristics Age of Oil Palm Based on SAR Imaging

The scattering characteristics of the forest vegetation are influenced by the thickness of the volume, density of the component particles (or scatters), size distribution of the component particles, shape distribution of the component particles, orientation distribution of the component particles, and dielectric properties of the component particles [58]. It is essential to know the polarization from which the SAR image obtained because signals at different polarizations interact differently with objects on the ground, which affects the brightness of the radar recorded in a particular polarization channel [59]. For simplicity, it is assumed that a natural scene can be described as a combination of three types of scatterers: (1) rough surface scatterers, (2) double-bounce scatterers, and (3) volume scatterers [60].


### 4.1.1. On the X-Band

The scattering characteristics of the oil palm at the X-band of the young oil palm have the lowest scattering value as seen in Figure10 of the 4-year-old histogram which is caused by rough surface scatters. The mature oil palm has a higher scattering value as seen in Figure 10 on the 6 to 15-year-old histogram which is caused by double-bounce scattering. The old oil palm has a lower scattering value as seen in Figure 10 of the 16 to 23-year-old histogram which is caused by volume scatters.Figure 10
Histogram X-band for each age of oil palm.According to Martinis and Rieke [61], adopting the X-band, canopy attenuation, volume, and surface scattering from the top layer of the forest canopy is generally higher [62] because the X-band is sensitive to canopy surface scattering [63]. Hence, the same canopy surface seems rough for X-band wavelength causing high scattering values [64]. According to Rosenqvist [65], the X-band instrument clear relationship between canopy growth area and threshold ambiguity for trees.


### 4.1.2. On the C-Band

The scattering characteristics of oil palm at the C-band on VV and VH polarization on the young oil palm have the lowest scattering value as seen in Figures11 and 12 on the histogram for 2 to 4-year-old histogram which is caused by rough surface scattering. The mature oil palm has a higher scattering value as seen in Figures 11 and 12 of the 7 to 14-year-old histogram which is caused by double-bounce scattering. The old oil palm has a lower scattering value as seen in Figures 11 and 12 of the 17 to 21-year-old histogram which is caused by volume scatters.Figure 11
Histogram C-band on VH for each age of oil palm.Figure 12
Histogram C-band on VV for each age of oil palm.According to Carolita et al. [15], increase and decrease in the scattering value occur because at the planting new age of trees height was relatively low and the size of the canopy is still small when the tree reached into productive age, the scattering value will rise because of the difference in the height of tree and soil. C-band data has a short wavelength, causing the scattering values for canopy density to become smaller. According to Teng et al. [6], there was also a significant contribution from the trunks to the surface volume scattering component, although this was not the dominant component. Also, at the shorter C-band wavelength, volume scattering can occur within the canopies of lower or sparse vegetation types, such as bushes, shrubs, or crops. According to Rignot et al. [66], radar signals were expected to be dominantly scattered by the foliage and the top branches and twigs of the canopy at the C-band.


### 4.1.3. On the L-Band

The scattering characteristics of oil palm at the L-band on HH and HV polarization on the young oil palm have the lowest scattering value as seen in Figures13 and 14 on the 5-year-old histogram which is caused by in surface scattering. The mature oil palm has a higher scattering value as seen in Figures 13 and 14 of the 7 to 14-year-old histogram which is caused by double-bounce scattering. The old oil palm has a lower scattering value, as seen in Figures 13 and 14 of the 16 to 24-year-old histogram which is caused by volume scatters.Figure 13
Histogram L-band on HV for each age of oil palm.Figure 14
Histogram L-band on HH for each age of oil palm.According to Toh et al. [67], the scattering increased with age for all polarization in general for the L-band frequency. The L-band scattering coefficient showed the trend increasing with the age of oil palm more than 20 years. Besides that, according to Teng et al. [6], the scattering value for HH polarization using the L-band was increasing attenuation of waves propagating through to the soil surface and back, as the leaves, fronds, and trunks are growing larger.The scattering value for each polarization (HH and HV) among the age of oil palm was increasing; according to Darmawan et al. [28], it was hypothesized that the growth of oil palms was caused by higher trunks, more leaves, more branches, and larger canopy followed by the condition of scattering value. In any case, the scatter values derived from the L-band are depending on wavelengths, polarization, incident angle, and temporal data; environments such as land moisture and landscape [25, 28, 68]; and the structure of oil palms such as size, geometry, and orientation of leaves, trunks, branches, and aerials or stilt roots [28, 68, 69].


## 4.1.1. On the X-Band

The scattering characteristics of the oil palm at the X-band of the young oil palm have the lowest scattering value as seen in Figure10 of the 4-year-old histogram which is caused by rough surface scatters. The mature oil palm has a higher scattering value as seen in Figure 10 on the 6 to 15-year-old histogram which is caused by double-bounce scattering. The old oil palm has a lower scattering value as seen in Figure 10 of the 16 to 23-year-old histogram which is caused by volume scatters.Figure 10
Histogram X-band for each age of oil palm.According to Martinis and Rieke [61], adopting the X-band, canopy attenuation, volume, and surface scattering from the top layer of the forest canopy is generally higher [62] because the X-band is sensitive to canopy surface scattering [63]. Hence, the same canopy surface seems rough for X-band wavelength causing high scattering values [64]. According to Rosenqvist [65], the X-band instrument clear relationship between canopy growth area and threshold ambiguity for trees.


## 4.1.2. On the C-Band

The scattering characteristics of oil palm at the C-band on VV and VH polarization on the young oil palm have the lowest scattering value as seen in Figures11 and 12 on the histogram for 2 to 4-year-old histogram which is caused by rough surface scattering. The mature oil palm has a higher scattering value as seen in Figures 11 and 12 of the 7 to 14-year-old histogram which is caused by double-bounce scattering. The old oil palm has a lower scattering value as seen in Figures 11 and 12 of the 17 to 21-year-old histogram which is caused by volume scatters.Figure 11
Histogram C-band on VH for each age of oil palm.Figure 12
Histogram C-band on VV for each age of oil palm.According to Carolita et al. [15], increase and decrease in the scattering value occur because at the planting new age of trees height was relatively low and the size of the canopy is still small when the tree reached into productive age, the scattering value will rise because of the difference in the height of tree and soil. C-band data has a short wavelength, causing the scattering values for canopy density to become smaller. According to Teng et al. [6], there was also a significant contribution from the trunks to the surface volume scattering component, although this was not the dominant component. Also, at the shorter C-band wavelength, volume scattering can occur within the canopies of lower or sparse vegetation types, such as bushes, shrubs, or crops. According to Rignot et al. [66], radar signals were expected to be dominantly scattered by the foliage and the top branches and twigs of the canopy at the C-band.


## 4.1.3. On the L-Band

The scattering characteristics of oil palm at the L-band on HH and HV polarization on the young oil palm have the lowest scattering value as seen in Figures13 and 14 on the 5-year-old histogram which is caused by in surface scattering. The mature oil palm has a higher scattering value as seen in Figures 13 and 14 of the 7 to 14-year-old histogram which is caused by double-bounce scattering. The old oil palm has a lower scattering value, as seen in Figures 13 and 14 of the 16 to 24-year-old histogram which is caused by volume scatters.Figure 13
Histogram L-band on HV for each age of oil palm.Figure 14
Histogram L-band on HH for each age of oil palm.According to Toh et al. [67], the scattering increased with age for all polarization in general for the L-band frequency. The L-band scattering coefficient showed the trend increasing with the age of oil palm more than 20 years. Besides that, according to Teng et al. [6], the scattering value for HH polarization using the L-band was increasing attenuation of waves propagating through to the soil surface and back, as the leaves, fronds, and trunks are growing larger.The scattering value for each polarization (HH and HV) among the age of oil palm was increasing; according to Darmawan et al. [28], it was hypothesized that the growth of oil palms was caused by higher trunks, more leaves, more branches, and larger canopy followed by the condition of scattering value. In any case, the scatter values derived from the L-band are depending on wavelengths, polarization, incident angle, and temporal data; environments such as land moisture and landscape [25, 28, 68]; and the structure of oil palms such as size, geometry, and orientation of leaves, trunks, branches, and aerials or stilt roots [28, 68, 69].


## 4.2. Coefficient Determination of Scattering Model of Oil Palm Phenology

The scatter model of oil palm phenology was performed separately for each radar wavelength (X-, C-, and L-band). In this study, the X-band showed the nonlinear oil palm phenology model in HH polarization with coefficient determinationR2=0.65, which means the model has moderate correlation. The C-band showed the nonlinear or quadratic model of oil palm phenology in VH polarization with coefficient determination R2=0.56, which means the model has moderate correlation, and while on VV polarization with coefficient determination R2=0.89, it means the model has strong correlation. The L-band showed the logarithmic model of oil palm phenology in HV polarization with coefficient determination R2=0.50, which means the model has moderate correlation, and while on HH polarization with coefficient determination R2=0.51, it means the model has moderate correlation. In this case, the more potential of scattering model of oil palm phenology based on R2 is using C-band on VV polarization.The scattering model with the moderate correlation also found by other researchers such as Tan et al. [10] has logarithmic regression for the correlating scattering value with oil palm age using ALOS PALSAR-2 with coefficient determination on HH polarizationR2=0.49, HV polarization withR2=0.27, and ratio of HH/HV polarization withR2=0.26, for study case in the southern part of peninsular Malaysia. Darmawan et al. [28] also used ALOS PALSAR-2 on HH and HV polarization having a logarithmic regression model with coefficient determination R2=0.62 and 0.41 for study case area in Jerantut, Pahang, Malaysia. In other cases, Avtar et al. [70] for study area near Miri City, Sarawak, Malaysia, showed a logarithmic regression model for a variety of SAR data including TerraSAR-X (X-band) on HH polarization with coefficient determination R2=0.075; Radarsat-2 (C-band) on HH, HV, and VV polarization with coefficient determination R2=0.39, 0.49, and 0.39; and ALOS PALSAR-2 (L-band) on HH and HV polarization with coefficient determination R2=0.62 and 0.77. Okarda et al. [32] studied the relationship between the age of oil palm in mineral soil and age and peatland in Central Kalimantan, Indonesia, using ALOS PALSAR-2 on HH and HV polarization. The scattering model has linear regression with coefficient determination on HH and HV between age and mineral soil with R2=0.55 and 0.37 and age and peatland on HH and HV with R2=0.36 and 0.28. For C-band data, it was studied by Carolita et al. [15] using Sentinel-1A on HH and HV polarization. The scattering model has nonlinear regression (quadratic model) with coefficient determination R2=0.68 for HH and HV with R2=0.77 for study case area in Asahan, North Sumatra Indonesia.


## 5. Conclusions

This study investigated the scattering model of oil palm phenology based on X-, C-, and L-band polarimetric SAR imaging in Asahan Regency, North Sumatra, Indonesia. We generated scattering values for different ranges of oil palm age: 0-5, 5-10, 10-15, 15-20, and 20-25 years, with different scattering characteristics for each X-, C-, and L-band polarimetric SAR imaging depending on the wavelength and age ranges of oil palm. Finally, we found the scattering model of oil palm phenology based on the X-band on HH polarization which isy=−0.0114x2+0.2741x–7.2823 with R2=0.65. The C-band on VH polarization is y=−0.0055x2+0.1248x-15.785 with R2=0.56, while that on VV polarization is y=−0.0117x2+0.3006x-8.9534 with R2=0.89. The L-band on HH polarization is y=0.5293lnx–12.616 with R2=0.51, while that on HV polarization is y=1.1555lnx–21.815 with R2=0.50. In this case, the more potential of the scattering model of oil palm phenology based on R2 is using C-band on VV polarization. However, the scattering model can be used and applied to identify the phenology of oil palm in Indonesia, which is the main parameter in yield estimation, also required for the tax estimation, replanting time, and identification of oil palm diseases. The scattering model being developed needs to improve accuracy by integrating multisource-multispectral data, including different wavelengths on optical data and microwaves and more in situ data.



---

*Source: 6625774-2021-03-08.xml*