Retrieval of vegetation height in rice fields using polarimetric SAR interferometry with TanDEM-X data

This work presents for the first time a demonstration with satellite data of polarimetric SAR interferometry (PolInSAR) applied to the retrieval of vegetation height in rice fields. Three series of dual-pol interferometric SAR data acquired with large baselines (2–3 km) by the TanDEM-X system during its science phase (April–September 2015) are exploited. A novel inversion algorithm especially suited for rice fields cultivated in flooded soil is proposed and evaluated. The validation is carried out over three test sites located in geographically different areas: Sevilla (SW Spain), Valencia (E Spain), and Ipsala (W Turkey), in which different rice types are present. Results are obtained during the whole growth cycle and demonstrate that PolInSAR is useful to produce accurate height estimates (RMSE 10–20 cm) when plants are tall enough (taller than 25–40 cm), without relying on external reference information.This work has been supported by the Spanish Ministry of Economy and Competitiveness (MINECO) and EU FEDER under project TIN2014-55413-C2-2-P. The research leading to these results has received funding from the European Union Seventh Framework Programme (FP7/2007–2013) under grant agreement 606983, and the Land-SAF (the EUMETSAT Network of Satellite Application Facilities) project. The in-situ measurements in the Ipsala site were conducted with the funding of The Scientific and Technological Research Council of Turkey (TUBITAK, Project No.: 113Y446)

Rice Plant Height Monitoring from Space with Bistatic Interferometry

This chapter provides an overview of the possibility to derive paddy rice plant heights with spaceborne bistatic SAR interferometry (InSAR). By using the only available interferometer in space, TanDEM-X, an investigation of rice crops located in Turkey is performed. Before analyzing the main outcomes, an introduction to the generation of elevation models with InSAR is provided, with a special focus on the agricultural land cover. The processing chain and the modifications foreseen to properly produce plant elevations and a roadmap for the quality assessment are described. The results obtained, with a very high interferometric coherence supporting an accurate estimation due to a limited electromagnetic wave penetration into the canopy, support a temporal change analysis on a field-by-field basis. For the purpose, an automatic approach to segment the fields without external auxiliary data is also provided. The study is concluded with an analysis of the impact of the wave polarization in the results

Vertical Accuracy of Freely Available Global Digital Elevation Models (ASTER, AW3D30, MERIT, TanDEM-X, SRTM, and NASADEM)

https://www.mdpi.com/2072-4292/12/21/348

Integration of LIDAR and IFSAR for mapping

LiDAR and IfSAR data is now widely used for a number of applications, particularly those needing a digital elevation model. The data is often complementary to other data such as aerial imagery and high resolution satellite data. This paper will review the current data sources and the products and then look at the ways in which the data can be integrated for particular applications. The main platforms for LiDAR are either helicopter or fixed wing aircraft, often operating at low altitudes, a digital camera is frequently included on the platform, there is an interest in using other sensors such as 3 line cameras of hyperspectral scanners. IfSAR is used from satellite platforms, or from aircraft, the latter are more compatible with LiDAR for integration. The paper will examine the advantages and disadvantages of LiDAR and IfSAR for DEM generation and discuss the issues which still need to be dealt with. Examples of applications will be given and particularly those involving the integration of different types of data. Examples will be given from various sources and future trends examined

Extracting Urban Morphology for Atmospheric Modeling from Multispectral and SAR Satellite Imagery

This paper presents an approach designed to derive an urban morphology map from satellite data while aiming to minimize the cost of data and user interference. The approach will help to provide updates to the current morphological databases around the world. The proposed urban morphology maps consist of two layers: 1) Digital Elevation Model (DEM) and 2) land cover map. Sentinel-2 data was used to create a land cover map, which was realized through image classification using optical range indices calculated from image data. For the purpose of atmospheric modeling, the most important classes are water and vegetation areas. The rest of the area includes bare soil and built-up areas among others, and they were merged into one class in the end. The classification result was validated with ground truth data collected both from field measurements and aerial imagery. The overall classification accuracy for the three classes is 91 %. TanDEM-X data was processed into two DEMs with different grid sizes using interferometric SAR processing. The resulting DEM has a RMSE of 3.2 meters compared to a high resolution DEM, which was estimated through 20 control points in flat areas. Comparing the derived DEM with the ground truth DEM from airborne LIDAR data, it can be seen that the street canyons, that are of high importance for urban atmospheric modeling are not detectable in the TanDEM-X DEM. However, the derived DEM is suitable for a class of urban atmospheric models. Based on the numerical modeling needs for regional atmospheric pollutant dispersion studies, the generated files enable the extraction of relevant parametrizations, such as Urban Canopy Parameters (UCP).Peer reviewe

Estimation of RVoG Scene Parameters by Means of PolInSAR With TanDEM-X Data: Effect of the Double-Bounce Contribution

This article evaluates the effect of the double-bounce (DB) decorrelation term that appears in single-pass bistatic acquisitions, as in the TanDEM-X system, on the inversion of scene parameters by means of polarimetric SAR interferometry (PolInSAR). The retrieval of all scene parameters involved in the Random Volume over Ground (RVoG) model (i.e., ground topography, vegetation height, extinction, and ground-to-volume ratios) is affected by this term when the radar response from the ground is dominated by the DB. The estimation error in all these parameters is analyzed by means of simulations over a wide range of system configurations and scene variables for both agricultural crops and forest scenarios. Simulations demonstrate that the inclusion of the DB term, which complicates the inversion algorithm, is necessary for the angles of incidence shallower than 30° to achieve an estimation error below 10% in vegetation height and to avoid a significant underestimation in the ground-to-volume ratios. At steep incidences, this decorrelation term does not affect the estimation of vegetation height and ground-to-volume ratios. Regarding the extinction, this parameter is intrinsically not well estimated, since most retrieved values are close to the initial guesses employed for the optimization algorithm, regardless of the use or not of the DB decorrelation term. Finally, these findings are compared with the experimental results from the TanDEM-X data acquired over the rice fields in Spain for the available system parameters (baseline and incidence angle) of the acquired data set.This work was supported in part by the Spanish Ministry of Science, Innovation and Universities, the State Agency of Research (AEI), and in part by the European Funds for Regional Development (EFRD) under Project TEC2017-85244-C2-1-P. The work of Noelia Romero-Puig was supported in part by the Generalitat Valenciana and in part by the European Social Fund (ESF) under Grant ACIF/2018/204

Biomass Representation in Synthetic Aperture Radar Interferometry Data Sets

This work makes an attempt to explain the origin, features and potential applications of the elevation bias of the synthetic aperture radar interferometry (InSAR) datasets over areas covered by vegetation. The rapid development of radar-based remote sensing methods, such as synthetic aperture radar (SAR) and InSAR, has provided an alternative to the photogrammetry and LiDAR for determining the third dimension of topographic surfaces. The InSAR method has proved to be so effective and productive that it allowed, within eleven days of the space shuttle mission, for acquisition of data to develop a three-dimensional model of almost the entire land surface of our planet. This mission is known as the Shuttle Radar Topography Mission (SRTM). Scientists across the geosciences were able to access the great benefits of uniformity, high resolution and the most precise digital elevation model (DEM) of the Earth like never before for their a wide variety of scientific and practical inquiries. Unfortunately, InSAR elevations misrepresent the surface of the Earth in places where there is substantial vegetation cover. This is a systematic error of unknown, yet limited (by the vertical extension of vegetation) magnitude. Up to now, only a limited number of attempts to model this error source have been made. However, none offer a robust remedy, but rather partial or case-based solutions. More work in this area of research is needed as the number of airborne and space-based InSAR elevation models has been steadily increasing over the last few years, despite strong competition from LiDAR and optical methods. From another perspective, however, this elevation bias, termed here as the “biomass impenetrability”, creates a great opportunity to learn about the biomass. This may be achieved due to the fact that the impenetrability can be considered a collective response to a few factors originating in 3D space that encompass the outermost boundaries of vegetation. The biomass, presence in InSAR datasets or simply the biomass impenetrability, is the focus of this research. The report, presented in a sequence of sections, gradually introduces terminology, physical and mathematical fundamentals commonly used in describing the propagation of electromagnetic waves, including the Maxwell equations. The synthetic aperture radar (SAR) and InSAR as active remote sensing methods are summarised. In subsequent steps, the major InSAR data sources and data acquisition systems, past and present, are outlined. Various examples of the InSAR datasets, including the SRTM C- and X-band elevation products and INTERMAP Inc. IFSAR digital terrain/surface models (DTM/DSM), representing diverse test sites in the world are used to demonstrate the presence and/or magnitude of the biomass impenetrability in the context of different types of vegetation – usually forest. Also, results of investigations carried out by selected researchers on the elevation bias in InSAR datasets and their attempts at mathematical modelling are reviewed. In recent years, a few researchers have suggested that the magnitude of the biomass impenetrability is linked to gaps in the vegetation cover. Based on these hints, a mathematical model of the tree and the forest has been developed. Three types of gaps were identified; gaps in the landscape-scale forest areas (Type 1), e.g. forest fire scares and logging areas; a gap between three trees forming a triangle (Type 2), e.g. depending on the shape of tree crowns; and gaps within a tree itself (Type 3). Experiments have demonstrated that Type 1 gaps follow the power-law density distribution function. One of the most useful features of the power-law distributed phenomena is their scale-independent property. This property was also used to model Type 3 gaps (within the tree crown) by assuming that these gaps follow the same distribution as the Type 1 gaps. A hypothesis was formulated regarding the penetration depth of the radar waves within the canopy. It claims that the depth of penetration is simply related to the quantisation level of the radar backscattered signal. A higher level of bits per pixels allows for capturing weaker signals arriving from the lower levels of the tree crown. Assuming certain generic and simplified shapes of tree crowns including cone, paraboloid, sphere and spherical cap, it was possible to model analytically Type 2 gaps. The Monte Carlo simulation method was used to investigate relationships between the impenetrability and various configurations of a modelled forest. One of the most important findings is that impenetrability is largely explainable by the gaps between trees. A much less important role is played by the penetrability into the crown cover. Another important finding is that the impenetrability strongly correlates with the vegetation density. Using this feature, a method for vegetation density mapping called the mean maximum impenetrability (MMI) method is proposed. Unlike the traditional methods of forest inventories, the MMI method allows for a much more realistic inventory of vegetation cover, because it is able to capture an in situ or current situation on the ground, but not for areas that are nominally classified as a “forest-to-be”. The MMI method also allows for the mapping of landscape variation in the forest or vegetation density, which is a novel and exciting feature of the new 3D remote sensing (3DRS) technique. Besides the inventory-type applications, the MMI method can be used as a forest change detection method. For maximum effectiveness of the MMI method, an object-based change detection approach is preferred. A minimum requirement for the MMI method is a time-lapsed reference dataset in the form, for example, of an existing forest map of the area of interest, or a vegetation density map prepared using InSAR datasets. Preliminary tests aimed at finding a degree of correlation between the impenetrability and other types of passive and active remote sensing data sources, including TerraSAR-X, NDVI and PALSAR, proved that the method most sensitive to vegetation density was the Japanese PALSAR - L-band SAR system. Unfortunately, PALSAR backscattered signals become very noisy for impenetrability below 15 m. This means that PALSAR has severe limitations for low loadings of the biomass per unit area. The proposed applications of the InSAR data will remain indispensable wherever cloud cover obscures the sky in a persistent manner, which makes suitable optical data acquisition extremely time-consuming or nearly impossible. A limitation of the MMI method is due to the fact that the impenetrability is calculated using a reference DTM, which must be available beforehand. In many countries around the world, appropriate quality DTMs are still unavailable. A possible solution to this obstacle is to use a DEM that was derived using P-band InSAR elevations or LiDAR. It must be noted, however, that in many cases, two InSAR datasets separated by time of the same area are sufficient for forest change detection or similar applications

Integration of high-resolution, Active and Passive Remote Sensing in support to Tsunami Preparedness and Contingency Planning

In the aftermath of the Sri Lanka tsunami disaster, a stack of synoptic procedures and remote sensing techniques was chosen for satisfying the urgent mapping needs of the Government. This choice presented the undebated advantage of (a) allowing to start the work immediately (b) without relying upon ground logistics until the onset of the air campaign, (c) minimizing the duration of the work on spot, while (d) covering fast - and at an otherwise unreacheable resolution - large portions of a difficult-to-penetrate territory, (e) keeping the work sustainable and, overall, (f) allowing to carry out the work. This combination of airborne and spaceborne techniques is ready-to-use worldwide, and the techniques for flooding simulation and scenario building can be chosen at whatever level of complexity - choosing preferably robustness. It is also worth noting further that the new generation of metric resolution, X-band Radar satellite constellations (as TerraSAR-X and Cosmo-SkyMED), may allow creating LiDAR-like products avoiding airborne missions. The products of the space-and-air campaign were handed over by the Ambassador of Italy to the Minister for Disaster Management and Humanitarian Affairs on 7th December 2006, Colombo, Sri Lanka

Radarkaugseire rakendused metsaüleujutuste ja põllumajanduslike rohumaade jälgimiseks

Väitekirja elektrooniline versioon ei sisalda publikatsioone.Käesolev doktoritöö keskendub radarkaugseire rakenduste arendamisele kahes keerukas looduskeskkonnas: üleujutatud metsas ja põllumajanduslikel rohumaadel. Uurimistöö viidi läbi Tartu Observatooriumis, Tartu Ülikoolis, Ventspilsi Kõrgkoolis ja Aalto Ülikoolis. Töö esimene osa käsitleb X-laineala polarimeetrilise radarisignaali käitumist regulaarselt üleujutatavas metsas Soomaa näitel ning teine osa põllumajanduslike rohumaade seisundi ja polarimeetriliste ning interferomeetriliste tehisava-radari parameetrite vahelisi seoseid. 2012 kevadel Soomaa testalal TerraSAR-X andmetega läbi viidud eksperiment näitas, et topelt-peegeldusele tundlik HH-VV polarimeetriline kanal pakub tõesti kontrastsemat tagasihajumisepõhist üleujutatud metsa eristust üleujutamata metsast kui traditsiooniline HH polarimeetriline kanal. HH-VV kanali eelis HH kanali ees on seda suurem, mida madalam on mets ning raagus tingimustes lehtmetsas oli HH-VV kanali eelis HH kanali ees suurem kui okasmetsas. Lisaks on üleujutusele tundlik HH ja VV kanali polarimeetriline faasivahe, mida on soovitatud ka varasemates töödes kasutada täiendava andmeallikana üleujutuste kaardistamisel. Käesolevas doktoritöös mõõdeti polarimeetrilise X-laineala tehisava-radari HH/VV faasivahe suurenemine üleujutuste tõttu erineva kõrgusega okas- ja lehtmetsas. 2013 a vegetatsiooniperioodil korraldati Rannu test-alal välimõõtmistega toetatud eksperiment uurimaks X- ja C-laineala polarimeetrilise ning X-laineala interferomeetrilise tehisava-radari parameetrite undlikkust rohumaade tingimuste muutustele. Ilmnes, et ühepäevase vahega kogutud X-laineala tehisava-radari interferomeetriliste paaride koherentsus korreleerus rohu kõrgusega. Koherentsus oli seda madalam, mida kõrgem oli rohi - leitud seost on võimalik potentsiaalselt rakendada niitmise tuvastamiseks. TerraSAR-X ja RADARSAT-2 polarimeetriliste aegridade analüüsi tulemusel leiti kaks niitmisele tundlikku parameetrit: HH/VV polarimeetriline koherentsus ja polarimeetriline entroopia. Niitmise järel langes HH/VV polarimeetriline koherentsus järsult ning polarimeetriline entroopia tõusis järsult. Rohu tagasikasvamise faasis hakkas HH/VV polarimeetriline koherentsus aeglaselt kasvama ning entroopia aeglaselt kahanema. Täheldatud efekt oli tugevam TerraSARX X-laineala aegridadel kui RADARSAT-2 C-riba tehisava-radari mõõtmistel ning seda selgemini nähtav mida rohkem biomassi niitmise järgselt maha jäi. Leitud HH/VV polarimeetrilise koherentsuse ja polarimeetrilise entroopia käitumine vastas taimkatte osakestepilve radarikiirguse tagasihajumismudelile. Mudeli järgi põhjus- 60 tas eelnimetatud parameetrite iseloomulikku muutust rohukõrte kui dipoolide orientatsiooni ja korrastatuse muut niitmise tõttu, mis on kooskõlas meie välimõõtmiste andmetega.This thesis presents research about the application of radar remote sensing for monitoring of complex natural environments, such as flooded forests and agricultural grasslands. The study was carried out in Tartu Observatory, University of Tartu, Ventspils University College, and Aalto University. The research consists of two distinctive parts devoted to polarimetric analysis of images from a seasonal flooding of wetlands, and to polarimetric and interferometric analysis of a summer-long campaign covering eleven agricultural grasslands. TerraSAR-X data from 2012 were used to assess the use of the double-bounce scattering mechanism for improving the mapping of flooded forest areas. The study confirmed that the HH–VV polarimetric channel that is sensitive to double-bounce scattering provides increased separation between flooded and unflooded forest areas when compared to the conventional HH channel. The increase in separation increases with decreasing forest height, and it is more pronounced for deciduous forests due to the leaf-off conditions during the study. The phase difference information provided by the HH–VV channel may provide additional information for delineating flooded and unflooded forest areas. Time series of X-band (TanDEM-X and COSMO-SkyMed) and C-band (RADARSAT-2) data from 2013 were analyzed in respect to vegetation parameters collected during a field survey. The one-day repeat-pass X-band interferometric coherence was shown to be correlated to the grassland vegetation height. The coherence was also found to be potentially useful for detecting mowing events. The polarimetric analysis of TanDEM-X and RADARSAT-2 data identified two parameters sensitive to mowing events - the HH/VV polarimetric coherence magnitude and the H2α entropy. Mowing of vegetation consistently caused the coherence magnitude to decrease and the entropy to increase. The effect was more pronounced in case of X-band data. Additionally, the effect was stronger with more vegetation left on the ground after mowing. The effect was explained using a vegetation particle scattering model. The changes in polarimetric variables was shown to be caused by the change of orientation and the randomness of the vegetation