A COMPARISON OF RANDOM FOREST AND LIGHT GRADIENT BOOSTING MACHINE FOR FOREST ABOVE-GROUND BIOMASS ESTIMATION USING A COMBINATION OF LANDSAT, ALOS PALSAR, AND AIRBORNE LIDAR DATA

Extracting damages caused by the 2008 Ms 8.0 Wenchuan earthquake from SAR remote sensing data

Satellite forest inventories are the only feasible way to map Canada’s vast, remote forest regions, such as those in the Northwest Territories (NWT). A method used to create such inventories is the k-nearest neighbour (k-NN) algorithm, which spatially extends information from forest inventory (FI) plots to the entire forest land base using wall-to-wall features typically derived from Landsat data. However, the benefits of integrating L-band synthetic aperture radar (SAR) data, strongly correlated to forest biomass, have not been assessed for Canadian northern boreal forests. Here we describe an optimized multivariate k-NN implementation of a 151,700 km2 area in southern NWT that included ca. 2007 Landsat and dual-polarized Phased Array type L-band SAR (PALSAR) data on board the Advanced Land Observing Satellite (ALOS). Five forest attributes were mapped at 30 m cells: stand height, crown closure, stand/total volume and aboveground biomass (AGB). We assessed accuracy gains compared to Landsat-based maps. To circumvent the scarcity of FI plots, we used 3600 footprints from the Geoscience Laser Altimeter System (GLAS) as surrogate FI plots, where forest attributes were estimated using Light Detection and Ranging (LiDAR) metrics as predictors. After optimization, k-NN predicted forest attribute values for each pixel as the average of the 4 nearest (k = 4) surrogate FI plots within the Euclidian space of 9 best features (selected among 6 PALSAR, 10 Landsat, and 6 environmental features). Accuracy comparisons were based on 31 National Forest Inventory ground plots and over 1 million airborne LiDAR plots. Maps that included PALSAR HV backscatter resulted in forest attribute predictions with higher goodness of fit (adj. R2), lower percent mean error (ME%), and percent root mean square error (RMSE%), and lower underestimation for larger attribute values. Predictions were most accurate for conifer stand height (RMSE% = 32.1%, adj. R2 = 0.58) and AGB (RMSE% = 47.8%, adj. R2 = 0.74), which is much more abundant in the area than mixedwood or broadleaf. Our study demonstrates that optimizing k-NN parameters and feature space, including PALSAR, Landsat, and environmental variables, is a viable approach for inventory mapping of the northern boreal forest regions of Canada.

Studies are needed to evaluate the ability of Advanced Land Observing Satellite (ALOS) Phased Array type L-band Synthetic Aperture Radar (PALSAR) for forest aboveground biomass (AGB) extraction in mountainous areas. In this article, forest biomass was estimated at plot and stand levels, and different biomass grades, respectively. Light detection and ranging (LiDAR) data with about one hit per m2 were first used for forest biomass estimation at the plot level, with R 2 of 0.77. Then the LiDAR-derived biomass, as prior knowledge, was used to investigate the relationship between ALOS PALSAR data and biomass. The results showed that at each biomass level, the range of the back-scatter coefficient in HH and HV polarization (where H and V represent horizontal and vertical polarizations, respectively, and the first of the two letters refers to the transmission polarization and the second to the received polarization) was very large and there was no obvious relationship between the synthetic aperture radar (SAR) back-scatter coefficient and biomass at plot level. At stand level and in different biomass grades, the back-scatter coefficient increased with the increase of forest biomass, and a logarithm equation can be used to describe the relationship. The main reason may be that forest structure is complex at the plot level, while the average value could partly decrease the influence of forest structure at stand level. Meanwhile, terrain radiometric correction (TRC) was investigated and found effective for forest biomass estimation.

The demarcation of linear geomorphological features due to fracture systems, underlying or exposed faults, drainage lines etc. from satellite images has been one of the most widely used applications in geological studies. Synthetic Aperture Radar (SAR) data has been recognized as an effective tool for lineament extraction due to its unique sensitivity to subtle topography, soil moisture and texture. However, the typical SAR sensor parameters such as frequency, polarization and look direction plays a vital role in enhancing the performance for a specific application. Therefore, the study aims to analyse the multi-frequency (C- and L-band), multi-polarization (like and cross) and multi-look direction SAR data for geological lineaments extraction. In this paper, Sentinel-1 and Advanced Land Observation Satellite (ALOS)-2 Phased Array type L-band Synthetic Aperture Radar (PALSAR)-2 datasets acquired over a part of the Indo-Burma ranges of Manipur region (northeastern part of India) are considered for analysis. The lineaments extracted from SAR datasets are compared with reference to lineaments provided by the Geological Survey of India (GSI) through their portal. The study observes that the like-polarization C-band SAR data provide optimal results for this application. It is also observed that the integration of lineaments extracted from multi-look SAR data can significantly improve the performance than that derived from single-look data. Overall, the study also focuses on the spatial distribution of lineaments with respect to the local structural regime for a better understanding of the regional tectonic setup.

Abstract. Providing an accurate above-ground biomass (AGB) map is of paramount importance for carbon stock and climate change monitoring. The main objective of this study is to compare the performance of pixel-based and object-based approaches for AGB estimation of temperate forests in north-eastern of New York State. Second, the capabilities of optical, SAR, and optical + SAR data were investigated. To achieve the goals, the random forest (RF) regression algorithm was used to model and predict the AGB values. Optical (i.e. Landsat 5TM, Landsat 8 OLI, and Sentinel-2), synthetic aperture radar (SAR) (Sentinel-1 and global phased array type L-band SAR (PALSAR/PALSAR-2)), and their integration have been used to estimate the AGB. It is worth mentioning that the airborne light detection and ranging (LiDAR) AGB raster has been used as a reference data for training/testing purposes. The results demonstrate that the OBIA approach enhanced the RMSE of AGB estimation about 5.32 Mg/ha, 8.9 Mg/ha, and 5.29 Mg/ha for optical, SAR, and optical + SAR data, respectively. Moreover, optical + SAR data with the RMSE of 42.63 Mg/ha and R2 of 0.72 for pixel-based and RMSE of 37.31 Mg/ha and R2 of 0.77 for object-based approach provided the best results.

Snow cover is an important parameter for hydrological modelling and climate change modelling. Various methods are available only for wet snow-cover mapping using conventional synthetic aperture radar (SAR) data. Total snow (wet + dry) cover mapping with SAR data is still a topical research area. Therefore, incoherent target decomposition theorems have been implemented on fully polarimetric SAR data to characterize the scattering of various targets. Further classification techniques – both unsupervised and supervised – have been applied for accurate mapping of total snow cover. For this purpose, Advanced Land Observing Satellite – phased array-type L-band SAR (ALOS–PALSAR) data (12 May 2007) have been analysed for snow classification of glaciated terrain in and around Badrinath region in Himalaya. An ALOS-Advanced Visible and Near Infrared Radiometer (AVNIR)-2 image (6 May 2007) was also used to provide assistance in the selection of different training classes. It has been found that the application of incoherent target decomposition theorems such as H/A/α and four-component scattering mechanism models are good for extracting the desired information of snow cover from fully polarimetric PALSAR data. Finally, based on these target decomposition theorems and the Wishart classifier, PALSAR data have been classified into snow or non-snow cover, and the user accuracy of snow classes was found to be better than the user accuracy of other classes. Hence, the application of incoherent target decomposition theorems with full polarimetric ALOS-PALSAR data is useful for snow-cover mapping.

Many supervised classification algorithms have been proposed, however, they are rarely evaluated for specific application. This research examines the performance of machine learning classifiers support vector machine (SVM), neural network (NN), Random Forest (RF) against maximum classifier (MLC) (traditional supervised classifier) in forest resources and land cover categorization, based on combination of Advanced Land Observing Satellite (ALOS) Phased Array type L-band Synthetic Aperture Radar (PALSAR) and Landsat Thematic Mapper (TM) data, in Northern Tanzania. Various data categories based on Landsat TM surface reflectance, ALOS PALSAR backscattering and their derivatives were generated for various classification scenarios. Then a separate and joint processing of Landsat and ALOS PALSAR data were executed using SVM, NN, RF and ML classifiers. The overall classification accuracy (OA), kappa coefficient (KC) and F1 score index values were computed. The result proves the robustness of SVM and RF in classification of forest resource and land cover using mere Landsat data and integration of Landsat and PALSAR (average OA = 92% and F1 = 0.7 to 1). A two sample t-statistics was utilized to evaluate the performance of the classifiers using different data categories. SVM and RF indicate there is no significance difference at 5% significance level. SVM and RF show a significant difference when compared to NN and ML. Generally, the study suggests that parametric classifiers indicate better performance compared to parametric classifier.

Recently land-use change has been the main concern for worldwide environment change and is being used by city and regional planners to design sustainable cities. Nakuru in the central Rift Valley of Kenya has undergone rapid urban growth in last decade. This paper focused on urban growth using multi-sensor satellite imageries and explored the potential benefits of combining data from optical sensors (Landsat, Worldview-2) with Radar sensor data from Advanced Land Observing Satellite (ALOS) Phased Array type L-band Synthetic Aperture Radar (PALSAR) data for urban land-use mapping. Landsat has sufficient spectral bands allowing for better delineation of urban green and impervious surface, Worldview-2 has a higher spatial resolution and facilitates urban growth mapping while PALSAR has higher temporal resolution compared to other operational sensors and has the capability of penetrating clouds irrespective of weather conditions and time of day, a condition prevalent in Nakuru, because it lies in a tropical area. Several classical and modern classifiers namely maximum likelihood (ML) and support vector machine (SVM) were applied for image classification and their performance assessed. The land-use data of the years 1986, 2000 and 2010 were compiled and analyzed using post classification comparison (PCC). The value of combining multi-temporal Landsat imagery and PALSAR was explored and achieved in this research. Our research illustrated that SVM algorithm yielded better results compared to ML. The integration of Landsat and ALOS PALSAR gave good results compared to when ALOS PAL- SAR was classified alone. 19.70 km2 of land changed to urban land-use from non-urban land-use between the years 2000 to 2010 indicating rapid urban growth has taken place. Land-use information is useful for the comprehensive land-use planning and an integrated management of resources to ensure sustainability of land and to achieve social Eq- uity, economic efficiency and environmental sustainability.

Fusion of optical, radar and waveform LiDAR observations for land cover classification

Fuel load is the key factor driving fire ignition, spread and intensity. The current literature reports the light detection and ranging (LiDAR), optical and airborne synthetic aperture radar (SAR) data for fuel load estimation, but the optical and SAR data are generally individually explored. Optical and SAR data are expected to be sensitive to different types of fuel loads because of their different imaging mechanisms. Optical data mainly captures the characteristics of leaf and forest canopy, while the latter is more sensitive to forest vertical structures due to its strong penetrability. This study aims to explore the performance of Landsat Enhanced Thematic Mapper Plus (ETM+) and Advanced Land Observing Satellite (ALOS) Phased Arrayed L-band Synthetic Aperture Radar (PALSAR) data as well as their combination on estimating three different types of fuel load—stem fuel load (SFL), branch fuel load (BFL) and foliage fuel load (FFL). We first analyzed the correlation between the three types of fuel load and optical and SAR data. Then, the partial least squares regression (PLSR) was used to build the fuel load estimation models based on the fuel load measurements from Vindeln, Sweden, and variables derived from optical and SAR data. Based on the leave-one-out cross-validation (LOOCV) method, results show that L-band SAR data performed well on all three types of fuel load (R2 = 0.72, 0.70, 0.72). The optical data performed best for FFL estimation (R2 = 0.66), followed by BFL (R2 = 0.56) and SFL (R2 = 0.37). Further improvements were found for the SFL, BFL and FFL estimation when integrating optical and SAR data (R2 = 0.76, 0.81, 0.82), highlighting the importance of data selection and combination for fuel load estimation.

It is well known that many of the nighttime acquisitions of the L-band Advanced Land Observing Satellite (ALOS) Phased Array-type L-band Synthetic Aperture Radar (PALSAR) instrument over equatorial regions show significant distortions of the image amplitude information. These distortions have the form of amplitude stripes that are roughly aligned with the local geomagnetic field. While ionospheric scintillation has been identified as the source of these distortions, the exact nature of the induced artifacts on synthetic aperture radar (SAR) image quality and SAR signal phase has not yet been studied in sufficient detail. Hence, this paper provides a quantitative analysis of equatorial scintillation effects on SAR image quality and SAR phase. We have performed a statistical analysis of ALOS PALSAR images over equatorial regions to describe the observed distortions and relate them to ionospheric parameters. An ionospheric simulator was developed and validated that is capable of simulating ionospheric distortions based on ionospheric scintillation parameters. Using this simulator, we found that ionospheric scintillation in the equatorial zone can cause significant distortions of SAR image amplitudes, image focus, and SAR signal phase. We determined threshold ionospheric environmental conditions that lead to the formation of these image distortions. Based on these thresholds, we quantified the likelihood of occurrence of ionospheric distortions for the global equatorial belt and for L-band sensors ALOS PALSAR, ALOS-2 PALSAR-2, and NASA-ISRO SAR (NISAR).

DOI:10.17014/ijog.6.2.123-140Mount Sinabung was re-activated at August 28th, 2010 after a long repose interval. The early stage of a phreatic eruption was then followed by magmatic eruptions at September 15th, 2013 for years until now. To understand the ground surface changes accompanying the eruption periods, comprehensive analyses of surface and subsurface data are necessary, especially the condition in pre- and syn-eruption periods. This study is raised to identify ground surface and topographical changes before, intra, and after the eruption periods by analyzing the temporal signature of surface roughness, moisture, and deformation derived from Synthetic Aperture Radar (SAR) data. The time series of SAR backscattering intensity were analyzed prior to and after the early eruption periods to know the lateral ground surface changes including estimated lava dome roughness and surface moisture. Meanwhile, the atmospherically corrected Differential Interferometric SAR (D-InSAR) method was also applied to know the vertical topographical changes prior to the eruptions. The atmospheric correction based on modified Referenced Linear Correlation (mRLC) was applied to each D-InSAR pair to exclude the atmospheric phase delay from the deformation signal. The changes of surface moistures on syn-eruptions were estimated by calculating dielectric constant from SAR polarimetric mode following Dubois model. Twenty-one Phased Array type L-band SAR (PALSAR) data on board Advanced Land Observing Satellite (ALOS) and nine Sentinel-1A SAR data were used in this study with the acquisition date between February 2006 and February 2017. For D-InSAR purposes, the ALOS PALSAR data were paired to generate twenty interferograms. Based on the D-InSAR deformation, three times inflation-deflation periods were observed prior to the early eruption at August 28th 2010. The first and second inflation-deflation periods at the end of 2008 and middle 2009 presented migration of magma batches and dike generations in the deep reservoir. The third inflation-deflation periods in the middle of 2010 served as a precursor signal presenting magma feeding to the shallow reservoir. The summit was inflated about 1.4 cm and followed by the eruptions. The deflation of about 2.3 cm indicated the release pressure and temperature in the shallow reservoir after the early eruption at August 28th, 2010. The last inflation-deflation period was also confirmed by the increase of the lava dome roughness size from 5,121 m2 on July to 6,584 m2 on August. The summit then inflated again about 1.1 cm after the first eruption and followed by unrest periods presented by lava dome growth and destruction at September 15th, 2013. The volcanic products including lava and pyroclastics strongly affected the moisture of surface layer. The volcanic products were observed to reduce the surface moisture within syn-eruption periods. The hot materials are presumed responsible for the evaporation of the surface moisture as well.

In this paper, we propose a new polarimetric synthetic aperture radar (SAR) data interpretation method based on a locally averaged Stokes vector. We first propose a method to extract discriminators from all three components of the averaged Stokes vector. Based on the extracted discriminators, we build four physical interpretation layers with ascending priorities, i.e., the basic structure layer, the low-coherence targets layer, the man-made targets layer, and the low-backscattering targets layer. An intuitive final image can be generated by simply stacking the four layers in the priority order. We test the performance of the proposed method over Advanced Land Observing Satellite Phased Array type L-band SAR (ALOS-PALSAR) data. Experimental results show that the proposed method has high interpretation performance, particularly for skew-aligned or randomly distributed buildings and isolated man-made targets such as bridges.

Digital Elevation Models (DEMs) contribute to geomorphological and hydrological applications. DEMs can be derived using different remote sensing-based datasets, such as Interferometric Synthetic Aperture Radar (InSAR) (e.g., Advanced Land Observing Satellite (ALOS) Phased Array type L-band SAR (PALSAR) and Shuttle Radar Topography Mission (SRTM) DEMs). In addition, there is also the Digital Surface Model (DSM) derived from optical tri-stereo ALOS Panchromatic Remote-sensing Instrument for Stereo Mapping (PRISM) imagery. In this study, we evaluated satellite-based DEMs, SRTM (Global) GL1 DEM V003 28.5 m, ALOS DSM 28.5 m, and PALSAR DEMs 12.5 m and 28.5 m, and their derived channel networks/orders. We carried out these assessments using Light Detection and Ranging (LiDAR) Digital Surface Models (DSMs) and Digital Terrain Models (DTMs) and their derived channel networks and Strahler orders as reference datasets at comparable spatial resolutions. We introduced a pixel-based method for the quantitative horizontal evaluation of the channel networks and Strahler orders derived from global DEMs utilizing confusion matrices at different flow accumulation area thresholds (ATs) and pixel buffer tolerance values (PBTVs) in both ±X and ±Y directions. A new Python toolbox for ArcGIS was developed to automate the introduced method. A set of evaluation metrics—(i) producer accuracy (PA), (ii) user accuracy (UA), (iii) F-score (F), and (iv) Cohen’s kappa index (KI)—were computed to evaluate the accuracy of the horizontal matching between channel networks/orders extracted from global DEMs and those derived from LiDAR DTMs and DSMs. PALSAR DEM 12.5 m ranked first among the other global DEMs with the lowest root mean square error (RMSE) and mean difference (MD) values of 4.57 m and 0.78 m, respectively, when compared to the LiDAR DTM 12.5 m. The ALOS DSM 28.5 m had the highest vertical accuracy with the lowest recorded RMSE and MD values of 4.01 m and −0.29 m, respectively, when compared to the LiDAR DSM 28.5 m. PALSAR DEM 12.5 m and ALOS DSM 28.5 m-derived channel networks/orders yielded the highest horizontal accuracy when compared to those delineated from LiDAR DTM 12.5 m and LiDAR DSM 28.5 m, respectively. The number of unmatched channels decreased when the PBTV increased from 0 to 3 pixels using different ATs.

We propose the use of SAR (Synthetic Aperture Radar) backscatter intensity data as a tool for monitoring coast. ALOS (Advanced Land Observing Satellite) PALSAR (Phased Array type L-band SAR) data are analyzed to investigate the temporal changes in shape and area of mouthbars in the Mekong River distributaries. River mouthbars with strong backscatter, which is surrounded by the water, are successfully extracted using a histogram thresholding algorithm. Estimated areas of river mouthbars gradually increase on an annual time scale. Seasonal variations of areas were also recognized, and they appear to be related to the mean monthly sea level, discharge, and rainfall. Moreover, it is highly likely that tidal height at the time of SAR data acquisition clearly reflect the area of river mouthbar by examining over half-day intervals for ascending and descending images. This study shows that if appropriately calibrated, the SAR data is useful to quantifying long-term changes in river mouthbar.