Lokaalstatistikute kasutamine rohumaade ja metsade kaugseires

Väitekirja elektrooniline versioon ei sisalda publikatsiooneKäesolev doktoritöö analüüsib lokaalstatistikute kasutamist rohumaade ja metsade kaugseires. Töö esimene osa käsitleb rohumaade monitoorimist tehisava-radari (synthetic aperture radar (SAR)) abil ning teine osa metsade kaugseiret kasutades optilisi sensoreid. Analüüsides rohumaade niitmise ja C- laineala tehisava-radari interferomeetrilise koherentsuse seoseid leiti, et selle parameetri kasutamisel on potentsiaali niitmise tuvastamise algoritmide ja rakenduste väljaarendamiseks. Tulemused näitavad, et pärast niitmist on VH ja VV polarisatsiooni 12-päeva interferomeetrilise koherentsuse mediaan väärtused statistiliselt oluliselt kõrgemad võrreldes niitmise eelse olukorraga. Koherentsus on seda kõrgem, mida väiksem on ajaline vahe niitmise ja pärast seda üles võetud esimese interferomeetrilise mõõtmise vahel. Hommikune kaste, sademed, põllutööde teostamine, näiteks külv või kündmine, kõrgelt niitmine ja kiire rohu kasv pärast niitmist vähendavad koherentsust ja raskendavad niitmise sündmuste eristamist. Selleks, et eelpoolnimetatud mõjusid leevendada tuleks tulevikus uurida 6-päeva koherentsuse ja niitmise sündmuste vahelisi seoseid. Käesolevas doktoritöös esitatud tulemused loovad siiski tugeva aluse edasisteks uuringuteks ja arendusteks eesmärgiga võtta C-laineala tehisava-radari andmed niitmise tuvastamisel ka praktikas kasutusele. Lisaks näidati, et ortofotodel põhinevate metsa kaugseire hinnangute andmisel on abi lokaalstatistikute kasutamisest. Analüüsides kaugseire hinnangut riigimetsa takseerandmete (national forest inventory) kohta leiti, et näidistel põhinev järeldamine (case-based reasoning (CBR)) sobib hästi selliste kaugseire ülesannete empiirilisteks lahendusteks, kus sisendandmetena on kasutatavad väga paljud erinevad andmeallikad. Leiti, et klasteranalüüsi saab kasutada kaugseire tunnuste eelvaliku meetodina. Võrreldes erinevaid tekstuuri statistikuid näidati, et lokaalselt arvutatud keskväärtus on kõige väärtuslikum tunnus. Järeldati, et nii statistiliste kui ka struktuursete lokaalstatistikute kasutamisega saab lisada pikslipõhistele kaugseire hinnangutele olulist andmestikku.This thesis studies approaches for remote sensing of grasslands and forests based on local statistics. The first part of the thesis focuses on monitoring of grasslands with SAR and the second part to monitoring of forests with optical sensors. It is shown that there is potential to develop mowing detection algorithms and applications using C-band SAR temporal interferometric coherence. The results demonstrate that after a mowing event, median VH and VV polarisation 12-day interferometric coherence values are statistically significantly higher than those from before the event. The sooner after the mowing event the first interferometric acquisition is taken, the higher the coherence. Morning dew, precipitation, farming activities, such as sowing or ploughing, high residual straws after the cut and rapid growth of grass are causing the coherence to decrease and impede the distinction of a mowing event. In the future, six-day interferometric coherence should also be analysed in relation to mowing events to alleviate some of these factors. Nevertheless, the results presented in this thesis offer a strong basis for further research and development activities towards the practical use of spaceborne C-band SAR data for mowing detection. Further, it was shown that local statistics can be useful for estimation of forest parameters from ortophotos and they could also provide helpful ancillary information to conduct a photo-interpretation tasks over forested areas. It was demonstrated that cluster analysis can be used as pre-selection method for the reduction of remote sensing features. Additionally, it was shown that case-based reasoning (a machine learning method) is well suited for empirical solutions of remote sensing tasks where there are many different data sources available. It was concluded that the use of local statistics adds valuable data to pixel-based remote sensing estimations

Improved forest height estimation by fusion of simulated GEDI Lidar data and TanDEM-X InSAR data

Interferometric Synthetic Aperture Radar (InSAR) and lidar are increasingly used active remote sensing techniques for forest structure observation. The TanDEM-X (TDX) InSAR mission of German Aerospace Center (DLR) and the upcoming Global Ecosystem Dynamics Investigation (GEDI) of National Aeronautics and Space Administration (NASA) together may provide more accurate estimates of global forest structure and biomass via their synergic use. In this paper, we explored the efficacy of simulated GEDI data in improving height estimates from TDX InSAR data. Our study sites span three major forest types: a temperate forest, a mountainous conifer forest, and a tropical rainforest. The GEDI lidar coverage was simulated for the full nominal two-year mission duration, under both cloud-free and 50%-cloud conditions. We then used these GEDI data to parameterize the Random Volume over Ground (RVoG) model driven by TDX imagery. In particular, we explored the following three strategies for forest structure estimation: 1) TDX data alone; 2) TDX + GEDI-derived digital terrain model (DTM); and 3) TDX + GEDI DTM + GEDI canopy height. We then validated the retrieved forest heights against wall-to-wall airborne lidar measurements. We found relatively large biases at 90 [m] spatial resolution, from 4.2–11.9 [m], and root mean square errors (RMSEs), from 7.9–12.7 [m] when using TDX data alone under constrained RVoG assumptions of a fixed extinction coefficient (σ) and a zero ground-to-volume amplitude ratio (μ = 0). Results improved significantly with the aid of a DTM derived from GEDI data which enabled estimation of spatially-varying σ values (vs. fixed extinction) under a μ = 0 assumption, with biases reduced to 1.7–4.2 [m] and RMSEs to 4.9–8.6 [m] across cloudy and cloud-free cases. The best agreement was achieved in the third strategy by also incorporating information of GEDI-derived canopy height to further enhance the RVoG parameters. The improved model, when still assuming μ = 0, reduced biases to less than or close to 1 m and further reduced RMSEs to 4.0–6.7 [m]. Finally, we used GEDI data to estimate spatially-varying μ in the RVoG model. We found biases of between −0.7–0.9 [m] and RMSEs in the range from 2.6–7.1 [m] over the three sites. Our results suggest that use of GEDI data improves height inversion from TDX, providing heights at more accuracy than can be achieved by TDX alone, and enabling wall-to-wall height estimation at much finer spatial resolution than can be achieved by GEDI alone

Hemiboreaalsete metsade kaardistamine interferomeetrilise tehisava-radari andmetelt

Väitekirja elektrooniline versioon ei sisalda publikatsioone.Käesolev doktoritöö uurib tehisavaradari (SAR) kasutusvõimalusi metsa kõrguse hindamiseks hemiboreaalsete metsade vööndis. Uurimistöö viidi läbi Tartu Üli¬kooli, Tartu Observatooriumi, Aalto Ülikooli, Euroopa Kosmoseagentuuri (ESA) kaugseire keskuse ESRIN ja Reach-U koostöös. Uurimistöös kasutatud satelliidi¬andmed on pärit Saksa Kosmosekeskuse (DLR) kõrglahutusega bistaatilise X-laineala tehisavaradari TanDEM-X satelliidipaarilt. Sagedasti uuenevad satelliidiandmed, nende globaalne katvus ja kõrge ruumi¬line lahutus võimaldavad tehisavaradari abil kaardistada metsi ning nendes toimu¬vaid muutusi suurtel maa-aladel. Radari abil on võimalik saada kõrge lahutusvõimega pilte, mis on tundlikud taimestikule, maapinna karedusele ja dielektrilistele omadustele. Sünkroonis lendava radaripaari samaaegselt tehtud pildid elimineerivad võimalikud ajalised muutused taimestikus ning tänu sellele on radariandmetest võimalik tuletada metsade vertikaalset struktuuri ja kõrgust. Uurimistöös käsitletakse tehisavaradari interferomeetrilise koherentsuse tund¬likkust metsa kõrguse suhtes ning analüüsitakse, millised keskkonna ja klimaati¬lised tingimused ning satelliidi orbiidiga seotud parameetrid mõjutavad radari¬piltidelt erinevate puuliikide kõrguse hindamise täpsust. Lisaks keskendub väitekiri interferomeetrilisele koherentsusele tuginevate mudelite analüüsi¬misele ning nende täpsuse hindamisele operatiivse metsa kõrguse kaardistamise raken-duseks. Vaatluse alla on võetud kolm testala, mis asuvad Soomaa rahvuspargis, Võrtsjärve idakaldal Rannus ja Peipsiveere looduskaitsealal ning katavad kokku 2291 hektarit metsa. 23 TanDEM-X satelliidipildi koherentsuspilte võrreldakse samadel testaladel aerolaserskaneerimise (LiDAR) abil mõõdetud puistute kõrgu¬sega, mis on omakorda jagatud kolme rühma (kuused, männid ja laia¬lehised segametsad). RVoG (Random Volume over Ground) taimekatte mudel ning sellest tule¬tatud lihtsamad pooleempiirilised mudelid sobituvad olemasolevate TanDEM-X koherentsuse ning LiDARi metsa puistute kõrgusandmetega hästi. Töö tule¬mused kinnitavad, et tulevikus on suurte ja erinevatest metsatüüpidest koosne¬vate metsade kõrguse kosmosest kaardistamisel otstarbekas kasutusele võtta esmalt just soovitatud lihtsamad ja universaalsemad mudelid.This thesis presents research in the field of radar remote sensing and contributes to the forest monitoring application development using space-borne synthetic aperture radar (SAR). Satellite data is particularly useful for large-scale forestry applications making high revisit monitoring of the state of forests worldwide possible. The sensitivity of SAR to the dielectric and geometrical properties of the targets, penetration capacity and coherent imaging properties make it a unique tool for mapping and monitoring forest biomes. SAR satellites are also capable of retrieving additional information about the structure of the forest, tree height and biomass estimates as an essential input for monitoring the changes in the carbon stocks. Interferometric SAR (InSAR) is an advanced SAR imaging technique that allows the retrieval of forest parameters while working in nearly all weather conditions, independently of daylight and cloud cover. This research concen¬trates on assessing the impact of different variables affecting hemiboreal forest height estimation from space-borne X-band interferometric SAR coherence data. In particular, the research analyses the changes in coherence dynamics related to seasonal conditions, tree species and imaging properties using a large collection of interferometric SAR images from different seasons over a four-year period. The study is carried out over three test sites in Estonia using the extensive multi-temporal dataset of 23 TanDEM-X images, covering 2291 hectares of forests to describe the relation between the interferometric SAR coherence mag¬nitude and forest parameters. The work demonstrates how the correlation of interferometric coherence and Airborne LiDAR Scanning (ALS)-derived forest height varies for pine and deciduous tree species, for summer (leaf-on) and winter (leaf-off) conditions and for flooded forest floor. A simple semi-empirical modelling approach is proposed as being suitable for wide area forest mapping with limited a priori information under a range of seasonal and environ¬¬mental conditions. A Random Volume over Ground (RVoG) model and three semi-empirical models are compared and validated against a large dataset of coherence magnitude and ALS-measured data over hemiboreal forests in Estonia. The results show that all proposed models perform well in describing the relationship between hemiboreal forest height and interferometric coherence, allowing in future to derive forest stand height with an accuracy suitable for a wide range of applications

Fusing simulated GEDI, ICESat-2 and NISAR data for regional aboveground biomass mapping

Accurate mapping of forest aboveground biomass (AGB) is critical for better understanding the role of forests in the global carbon cycle. NASA's current GEDI and ICESat-2 missions as well as the upcoming NISAR mission will collect synergistic data with different coverage and sensitivity to AGB. In this study, we present a multi-sensor data fusion approach leveraging the strength of each mission to produce wall-to-wall AGB maps that are more accurate and spatially comprehensive than what is achievable with any one sensor alone. Specifically, we calibrate a regional L-band radar AGB model using the sparse, simulated spaceborne lidar AGB estimates. We assess our data fusion framework using simulations of GEDI, ICESat-2 and NISAR data from airborne laser scanning (ALS) and UAVSAR data acquired over the temperate high AGB forest and complex terrain in Sonoma County, California, USA. For ICESat-2 and GEDI missions, we simulate two years of data coverage and AGB at footprint level are estimated using realistic AGB models. We compare the performance of our fusion framework when different combinations of the sparse simulated GEDI and ICEsat-2 AGB estimates are used to calibrate our regional L-band AGB models. In addition, we test our framework at Sonoma using (a) 1-ha square grid cells and (b) similarly sized irregularly shaped objects. We demonstrate that the estimated mean AGB across Sonoma is more accurately estimated using our fusion framework than using GEDI or ICESat-2 mission data alone, either with a regular grid or with irregular segments as mapping units. This research highlights methodological opportunities for fusing new and upcoming active remote sensing data streams toward improved AGB mapping through data fusion.</p

FUSING GEDI LIDAR AND TANDEM-X INSAR OBSERVATIONS FOR IMPROVED FOREST STRUCTURE AND BIOMASS MAPPING

The upcoming NASA’s Global Ecosystem Dynamics Investigation (GEDI) mission presents an unprecedented opportunity to advance current global biomass estimates. However, gaps are expected between GEDI’s ground tracks, requiring the development of fusion-based methodologies to contiguously map forest biomass at satisfactory resolutions and accuracies. This dissertation is built on the complementary advantages of observations from GEDI and DLR’s TerraSAR-X/TanDEM-X (TDX)) Interferometric Synthetic Aperture Radar (InSAR) mission. To meet the goal of mapping forest structure and biomass contiguously and accurately, three types of fusion strategies have been investigated. First, a simulated GEDI-derived digital terrain model (DTM) was utilized to improve height estimation from TDX. Forest heights were initially derived from TDX coherence alone as a baseline using the widely used Random Volume over Ground (RVoG) scattering model. Here, assumptions about RVoG parameters – extinction coefficient (σ) and ground-to-volume amplitude ratio (µ) – were made. Using an external DTM derived from simulated GEDI lidar data, RVoG model was used to calculate spatially varied σ values and derived forest heights with better accuracy. TDX forest height estimation was further improved with the aid of simulated GEDI-derived DTM and canopy heights. The additional use of simulated GEDI canopy heights as RVoG input not just refined σ but also enabled the estimation of µ. Based on these parameters, forest heights were improved across three different forest types; biases were reduced from 1.7–3.8 m using only simulated GEDI DTMs to -0.9–1.1 m by using both simulated GEDI DTMs and canopy heights. Finally, wall-to-wall TDX heights were used to improve biomass estimates from simulated GEDI data over three contrasting forest types. When using simulated GEDI sampled observations alone, uncertainties were estimated statistically to be 9.0–19.9% at 1 km. These were improved to 5.2–11.7% at the same resolution by upscaling simulated GEDI footprint biomass with TDX heights. The GEDI/TDX data fusion also enabled the generation of biomass maps at a fine spatial resolution of 100 m, with uncertainties estimated to be 6.0–14.0%. Through the exploration of these fusion strategies, it has been demonstrated that a fusion-based mapping method could realize the generation of forest biomass products from GEDI with unprecedented resolutions and accuracies, while taking advantage of global seamless observations from TDX

The SAR Handbook: Comprehensive Methodologies for Forest Monitoring and Biomass Estimation

This Synthetic Aperture Radar (SAR) handbook of applied methods for forest monitoring and biomass estimation has been developed by SERVIR in collaboration with SilvaCarbon to address pressing needs in the development of operational forest monitoring services. Despite the existence of SAR technology with all-weather capability for over 30 years, the applied use of this technology for operational purposes has proven difficult. This handbook seeks to provide understandable, easy-to-assimilate technical material to remote sensing specialists that may not have expertise on SAR but are interested in leveraging SAR technology in the forestry sector

Elevation and Deformation Extraction from TomoSAR

3D SAR tomography (TomoSAR) and 4D SAR differential tomography (Diff-TomoSAR) exploit multi-baseline SAR data stacks to provide an essential innovation of SAR Interferometry for many applications, sensing complex scenes with multiple scatterers mapped into the same SAR pixel cell. However, these are still influenced by DEM uncertainty, temporal decorrelation, orbital, tropospheric and ionospheric phase distortion and height blurring. In this thesis, these techniques are explored. As part of this exploration, the systematic procedures for DEM generation, DEM quality assessment, DEM quality improvement and DEM applications are first studied. Besides, this thesis focuses on the whole cycle of systematic methods for 3D & 4D TomoSAR imaging for height and deformation retrieval, from the problem formation phase, through the development of methods to testing on real SAR data. After DEM generation introduction from spaceborne bistatic InSAR (TanDEM-X) and airborne photogrammetry (Bluesky), a new DEM co-registration method with line feature validation (river network line, ridgeline, valley line, crater boundary feature and so on) is developed and demonstrated to assist the study of a wide area DEM data quality. This DEM co-registration method aligns two DEMs irrespective of the linear distortion model, which improves the quality of DEM vertical comparison accuracy significantly and is suitable and helpful for DEM quality assessment. A systematic TomoSAR algorithm and method have been established, tested, analysed and demonstrated for various applications (urban buildings, bridges, dams) to achieve better 3D & 4D tomographic SAR imaging results. These include applying Cosmo-Skymed X band single-polarisation data over the Zipingpu dam, Dujiangyan, Sichuan, China, to map topography; and using ALOS L band data in the San Francisco Bay region to map urban building and bridge. A new ionospheric correction method based on the tile method employing IGS TEC data, a split-spectrum and an ionospheric model via least squares are developed to correct ionospheric distortion to improve the accuracy of 3D & 4D tomographic SAR imaging. Meanwhile, a pixel by pixel orbit baseline estimation method is developed to address the research gaps of baseline estimation for 3D & 4D spaceborne SAR tomography imaging. Moreover, a SAR tomography imaging algorithm and a differential tomography four-dimensional SAR imaging algorithm based on compressive sensing, SAR interferometry phase (InSAR) calibration reference to DEM with DEM error correction, a new phase error calibration and compensation algorithm, based on PS, SVD, PGA, weighted least squares and minimum entropy, are developed to obtain accurate 3D & 4D tomographic SAR imaging results. The new baseline estimation method and consequent TomoSAR processing results showed that an accurate baseline estimation is essential to build up the TomoSAR model. After baseline estimation, phase calibration experiments (via FFT and Capon method) indicate that a phase calibration step is indispensable for TomoSAR imaging, which eventually influences the inversion results. A super-resolution reconstruction CS based study demonstrates X band data with the CS method does not fit for forest reconstruction but works for reconstruction of large civil engineering structures such as dams and urban buildings. Meanwhile, the L band data with FFT, Capon and the CS method are shown to work for the reconstruction of large manmade structures (such as bridges) and urban buildings

Establishing the sensitivity of Synthetic Aperture Radar to above-ground biomass in wooded savannas

Radar for biomass estimation has been widely investigated for temperate, boreal and tropical forests, yet tropical savanna woodlands, which generally form non-continuous cover canopies or sparse woodlands, have been largely neglected in biomass studies. This thesis evaluates the capability of Synthetic Aperture Radar (SAR) for estimating the above-ground biomass of the woody vegetation in a savanna in Belize, Central America. This is achieved by evaluating (i) polarimetric Synthetic Aperture Radar (SAR) backscatter and (ii) single-pass shortwave interferometric SAR (InSAR) as indicators of above-ground biomass. Specifically, the effect on SAR backscatter of woody vegetation structure such as canopy cover, basal area, vegetation height and above-ground biomass is evaluated. Since vegetation height is often correlated to above-ground biomass, the effectiveness of vegetation height retrieval from InSAR is evaluated as an indicator of above-ground biomass. The study area, situated in Belize, is representative of Central American savannas. Radar data used are AIRSAR fully polarimetric L- and P-band SAR, and AIRSAR C-band InSAR, Intermap Technologies STAR-3i X-band InSAR, and Shuttle Radar Topography Mission (SRTM) C-band InSAR. The field data comprise accurately georeferenced three-dimensional measurements for 1,133 trees and shrubs and 75 palmetto clumps and thickets in a transect of 800 m x 60 m which spans the main savanna vegetation strata of the study area. An additional 2,464 ground points were observed. Results show that savanna woodlands present a challenge for radar remote sensing methods due to the sparse and heterogeneous nature of savanna woodlands. Long-wave SAR backscatter is dominated not only by high biomass areas, but also by areas of leafy palmetto which have low vegetative biomass. Retrieved woodland canopy heights from X- and C-band InSAR are indicative of the general patterns of tree height, although retrieved heights are underestimated. The amount of underestimation is variable across the different canopy conditions. Of these two methods, the shortwave InSAR data give a better indication of the spatial distribution of the above-ground biomass of the woody vegetation in the savannas than SAR backscatter. These results have implications for new and planned future global biomass estimation missions, such as ALOS PALSAR, ESA’s planned P-band BIOMASS and TanDEM-X. Without appropriate mediation, SAR backscatter methods might overestimate above-ground biomass of the woody vegetation of savannas while InSAR height retrieval methods might underestimate biomass estimates. Some possible mediating approaches are discussed

Radar and multispectral remote sensing data accurately estimate vegetation vertical structure diversity as a fire resilience indicator

The structural complexity of plant communities contributes to maintaining the ecosystem functioning in fire-prone landscapes and plays a crucial role in driving ecological resilience to fire. The objective of this study was to evaluate the resilience to fire off several plant communities with reference to the temporal evolution of their vertical structural diversity (VSD) estimated from the data fusion of C-band synthetic aperture radar (SAR) backscatter (Sentinel-1) and multispectral remote sensing reflectance (Sentinel-2) in a burned landscape of the western Mediterranean Basin. We estimated VSD in the field 1 and 2 years after fire using Shannon’s index as a measure of vertical heterogeneity in vegetation structure from the vegetation cover in several strata, both in burned and unburned control plots. Random forest (RF) was used to model VSD in the control (analogous to prefire scenario) and burned plots (1 year after fire) using as predictors (i) Sentinel-1 VV and VH backscatter coefficients and (ii) surface reflectance of Sentinel-2 bands. The transferability of the RF model from 1 to 2 years after wildfire was also evaluated. We generated VSD prediction maps across the study site for the prefire scenario and 1 to 4 years postfire. RF models accurately explained VSD in unburned control plots (R2 = 87.68; RMSE = 0.16) and burned plots 1 year after fire (R2 = 80.48; RMSE = 0.13). RF model transferability only involved a reduction in the VSD predictive capacity from 0.13 to 0.20 in terms of RMSE. The VSD of each plant community 4 years after the fire disturbance was significantly lower than in the prefire scenario. Plant communities dominated by resprouter species featured significantly higher VSD recovery values than communities dominated by facultative or obligate seeders. Our results support the applicability of SAR and multispectral data fusion for monitoring VSD as a generalizable resilience indicator in fire-prone landscapes.SIEuropean Regional Development Fund (ERDF)Spanish Ministry of Economy and CompetitivenessRegional Government of Castilla and LeónBritish Ecological Societ