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

Puistute takseertunnuste hindamine aerolidari mõõtmisandmete põhjal hemiboreaalsetes metsades

A Thesis for applying for the degree of Doctor of Philosophy in Forestry.Forest management and planning requires up-to-date data, which commonly is acquired using field experts and ground measurements. Nowadays, more and more of data about forest stands is measured using remotely sensing methods. Most common methods include aerial photography and laser scanning from airplanes, also spectral measurements from satellites or even drone images and applications. This doctoral thesis focuses on developing applications and methods for utilising the airborne laser scanning (ALS) data that is freely available for the whole Estonia. The ALS measurements are carried out by the Estonian Land Board on a routine basis twice a year – in spring and summer. The first variable that was studied in this thesis was forest height. Based on the thesis, the most reliable method for forest height assessment was using the ALS point-cloud 80th height percentile (HP80). The small circular plot (radius of 15…30 m) and stand based studies showed high correlations with the field-measured forest heights and with great confidence it can be said, that ALS-based forest height estimations are close or even with higher accuracy, than field inspected. The second studied variable was standing wood volume. The ALS-based methods and models that were developed throughout this thesis used the idea, that standing wood volume is based on forest height and density. For this the HP80 and a threshold-based point count ratio was used (canopy cover - CC). ALS-based CC (CCALS) estimates were studied and compared with digital hemispherical photo based measurements. The results showed similar errors as were shown in other similar studies, with around 10-15% root mean square error (RMSE). The strongest correlation was shown using all echoes above a 1.3 metre threshold. Combining the CCALS and HP80 showed standing wood volume estimates with a similar error as we would receive from field measurements (<20%). The freely available multitemporal ALS data showed promising results for forest height growth monitoring and detecting small-scale disturbances. CCALS was shown to have strong predictive value, when compared with a four year difference in thinned and unthinned stands. The nation-wide ALS data can also be combined with forest height predictions from satellites, providing a faster update compared to the ALS data. Promising results were shown using the interferometric synthetic aperture radar (InSAR). Stand species maps generated using self-learning algorithms and satellite based spectral data can be used for developing species specific models of standing wood volume prediction. By combining these different datasets we can construct a nation-wide forest resource to help make better decisions for forest management and targeting fieldwork.Metsades majandamisotsuste langetamiseks ja metsamajanduslike tööde planeerimiseks on metsaomanikel vaja andmeid. Harjumuspäraselt on andmete kogumiseks tehtud metsas maapealseid mõõtmisi. Viimastel aastakümnetel on metsade inventeerimiseks üha enam aga kasutatud mittekontaktseid mõõtmisi - lennukitelt tehtavad aerofotosid, laserskaneerimist, satelliitidelt tehtavaid kiirgusmõõtmisi või viimastel aastatel ka droonidelt tehtud pilte. Antud doktoritöö on võtnud fookusesse aerolaserskaneerimise (ALS) andmete põhjal Eesti metsadesse sobilike rakenduste väljatöötamise. ALS mõõtmisi teeb Eesti Maa-amet rutiinsete lendude käigus kaks korda aastas, nii kevadel kui ka suvel. Aastast 2008 alustatud mõõtmiste tulemusel on Eesti üks väheseid riike maailmas, kus on vabalt kasutada mitmekordselt kogu riiki kattev ALS andmestik. Doktoritöö tulemusel töötati välja metsa kõrguse hindamiseks sobilikud meetodid, kasutades selleks punktipilvede kõrgusprotsentiile. Tugevamaid seoseid metsas proovitükkidel mõõdetud kõrgustega näitas punktipilve 80-protsentiil (HP80) ja uuringute põhjal võib väita, et metsa kõrguse määramine suvistelt aerolidari andmetelt on ligilähedane täpsustele, mida saadakse metsas kohapeal mõõtes. Teine oluline tunnus, mida metsade majandamise planeerimisel silmas peetakse, on kasvava metsa tagavara. Teadustöö põhjal töötati välja mudelite kujud ja metoodika, mille abil prognoositud tagavara oli sarnase veapiiriga, mis on lubatud metsas hinnanguid tegevatele taksaatoritele (<20%). Väljatöötatud ALS-põhine mudeli kuju järgib loogikat, et metsa tagavara on otseselt seotud mõõdetud kõrguse ja metsa tihedusega. Tihenduse hindamiseks aerolidari andmetelt kasutatakse nivoopõhist punktide suhtearvu ehk nn katvushinnangut (CCALS). Katvushinnangu täpsuse valideerimiseks ja tihedas metsas sobiva prognoosimeetodi väljatöötamiseks tehti välimõõtmisi kasutades poolsfäärikaameraid. Poolsfääripiltide põhjal tehtud valideerimise tulemused andsid sarnaseid veahinnanguid, mida on ka varasemates teadusuuringutes esitletud (RMSE = 10…15%). Kahe sarnasest fenoloogilisest perioodist ALS andmestiku lahutamisel uuriti ka muutuste tuvastamise võimalikkust. Uuringud andsid paljulubavaid tulemusi metsade kõrguskasvu hindamiseks ja CCALS osutus ka oluliseks tunnuseks väiksemate häiringute, nagu näiteks harvendusraie, tuvastamiseks. Kogu riiki katva ALS andmestiku kombineerimisel erinevate satelliitandmetega või näiteks spektraalsete mõõtmiste põhjal tehtud puistu liigiliste koosseisu kaartidega on võimalik antud töös välja pakutud meetodite abil anda igal aastal kogu Eesti metsaressursside ülevaade. Samuti on võimalik koostada vaid kaugseirevahendeid ja proovitükkidel lähendatud mudeleid kasutades eraldiste põhised takseerkirjeldused, mida siis taksaatorid saavad näiteks kasutada oma välitööde kavandamisel.  Publication of this thesis is supported by the Estonian University of Life Sciences

Utilization of bistatic TanDEM-X data to derive land cover information

Forests have significance as carbon sink in climate change. Therefore, it is of high importance to track land use changes as well as to estimate the state as carbon sink. This is useful for sustainable forest management, land use planning, carbon modelling, and support to implement international initiatives like REDD+ (Reducing Emissions from Deforestation and Degradation). A combination of field measurements and remote sensing seems most suitable to monitor forests. Radar sensors are considered as high potential due to the weather and daytime independence. TanDEM-X is a interferometric SAR (synthetic aperture radar) mission in space and can be used for land use monitoring as well as estimation of biophysical parameters. TanDEM-X is a X-band system resulting in low penetration depth into the forest canopy. Interferometric information can be useful, whereas the low penetration can be considered as an advantage. The interferometric height is assumable as canopy height, which is correlated with forest biomass. Furthermore, the interferometric coherence is mainly governed by volume decorrelation, whereas temporal decorrelation is minimized. This information can be valuable for quantitative estimations and land use monitoring. The interferometric coherence improved results in comparison to land use classifications without coherence of about 10% (75% vs. 85%). Especially the differentiation between forest classes profited from coherence. The coherence correlated with aboveground biomass in a R² of about 0.5 and resulted in a root mean square error (RSME) of 14%. The interferometric height achieved an even higher correlation with the biomass (R²=0.68) resulting in cross-validated RMSE of 7.5%. These results indicated that TanDEM-X can be considered as valuable and consistent data source for forest monitoring. Especially interferometric information seemed suitable for biomass estimation

Deep Learning Model Transfer in Forest Mapping Using Multi-Source Satellite SAR and Optical Images

Deep learning (DL) models are gaining popularity in forest variable prediction using Earth observation (EO) images. However, in practical forest inventories, reference datasets are often represented by plot- or stand-level measurements, while high-quality representative wall-to-wall reference data for end-to-end training of DL models are rarely available. Transfer learning facilitates expansion of the use of deep learning models into areas with sub-optimal training data by allowing pretraining of the model in areas where high-quality teaching data are available. In this study, we perform a “model transfer” (or domain adaptation) of a pretrained DL model into a target area using plot-level measurements and compare performance versus other machine learning models. We use an earlier developed UNet based model (SeUNet) to demonstrate the approach on two distinct taiga sites with varying forest structure and composition. The examined SeUNet model uses multi-source EO data to predict forest height. Here, EO data are represented by a combination of Copernicus Sentinel-1 C-band SAR and Sentinel-2 multispectral images, ALOS-2 PALSAR-2 SAR mosaics and TanDEM-X bistatic interferometric radar data. The training study site is located in Finnish Lapland, while the target site is located in Southern Finland. By leveraging transfer learning, the SeUNet prediction achieved root mean squared error (RMSE) of (Formula presented.) m and R2 of 0.882, considerably more accurate than traditional benchmark methods. We expect such forest-specific DL model transfer can be suitable also for other forest variables and other EO data sources that are sensitive to forest structure.</p

Application of the Trace Coherence to HH-VV PolInSAR TanDEM-X Data for Vegetation Height Estimation

This article investigates, for the first time, the inclusion of the operator Trace Coherence (TrCoh) in polarimetric and interferometric synthetic aperture radar (SAR) methodologies for the estimation of biophysical parameters of vegetation. A modified inversion algorithm based on the well-known Random Volume over Ground (RVoG) model, which employs the TrCoh, is described and evaluated. In this regard, a different set of coherence extrema is used as input for the retrieval stage. In addition, the proposed methodology improves the inversion algorithm by employing analytical solutions rather than approximations. Validation is carried out exploiting single-pass HH-VV bistatic TanDEM-X data, together with reference data acquired over a paddy rice area in Spain. The added value of the TrCoh and the convenience of the use of analytical solutions are assessed by comparing with the conventional polarimetric SAR interferometry (PolInSAR) algorithm. Results demonstrate that the modified proposed methodology is computationally more effective than current methods on this dataset. For the same scene, the steps required for inversion are computed in 6 min with the conventional method, while it only takes 6 s with the proposed approach. Moreover, vegetation height estimates exhibit a higher accuracy with the proposed method in all fields under evaluation. The root-mean-squared error reached with the modified method improves by 7 cm with respect to the conventional algorithm.This work was supported in part by the Spanish Ministry of Science and Innovation, in part by 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 Generalitat Valenciana and in part by the European Social Fund (ESF) under Grant ACIF/2018/204

Comparison between TanDEM-X- and ALS-based estimation of aboveground biomass and tree height in boreal forests

Interferometric Synthetic Aperture Radar (InSAR) data from TerraSAR-X add-on for Digital Elevation Measurement (TanDEM-X) were used to estimate aboveground biomass (AGB) and tree height with linear regression models. These were compared to models based on airborne laser scanning (ALS) data at two Swedish boreal forest test sites, Krycklan (64 degrees N19 degrees E) and Remningstorp (58 degrees N13 degrees E). The predictions were validated using field data at the stand-level (0.5-26.1 ha) and at the plot-level (10 m radius). Additionally, the ALS metrics percentile 99 (p99) and vegetation ratio, commonly used to estimate AGB and tree height, were estimated in order to investigate the feasibility of replacing ALS data with TanDEM-X InSAR data. Both AGB and tree height could be estimated with about the same accuracy at the stand-level from both TanDEM-X- and ALS-based data. The AGB was estimated with 17.2% and 14.6% root mean square error (RMSE) and the tree height with 7.6% and 4.1% RMSE from TanDEM-X data at the stand-level at the two test sites Krycklan and Remningstorp. The Pearson correlation coefficients between the TanDEM-X height and the ALS height p99 were r=.98 and r=.95 at the two test sites. The TanDEM-X height contains information related to both tree height and forest density, which was validated from several estimation models

Advances in Radar Remote Sensing of Agricultural Crops: A Review

There are enormous advantages of a review article in the field of emerging technology like radar remote sensing applications in agriculture. This paper aims to report select recent advancements in the field of Synthetic Aperture Radar (SAR) remote sensing of crops. In order to make the paper comprehensive and more meaningful for the readers, an attempt has also been made to include discussion on various technologies of SAR sensors used for remote sensing of agricultural crops viz. basic SAR sensor, SAR interferometry (InSAR), SAR polarimetry (PolSAR) and polarimetric interferometry SAR (PolInSAR). The paper covers all the methodologies used for various agricultural applications like empirically based models, machine learning based models and radiative transfer theorem based models. A thorough literature review of more than 100 research papers indicates that SAR polarimetry can be used effectively for crop inventory and biophysical parameters estimation such are leaf area index, plant water content, and biomass but shown less sensitivity towards plant height as compared to SAR interferometry. Polarimetric SAR Interferometry is preferable for taking advantage of both SAR polarimetry and SAR interferometry. Numerous studies based upon multi-parametric SAR indicate that optimum selection of SAR sensor parameters enhances SAR sensitivity as a whole for various agricultural applications. It has been observed that researchers are widely using three models such are empirical, machine learning and radiative transfer theorem based models. Machine learning based models are identified as a better approach for crop monitoring using radar remote sensing data. It is expected that the review article will not only generate interest amongst the readers to explore and exploit radar remote sensing for various agricultural applications but also provide a ready reference to the researchers working in this field

Forest attributes mapping with SAR data in the romanian South-Eastern Carpathians requirements and outcomes

Esta tesis doctoral se centra en la estimación de variables forestales en la zona Sureste de los Cárpatos Rumanos a partir de imágenes de radar de apertura sintética. La investigación abarca parte del preprocesado de las imágenes, métodos de generación de mosaicos y la extracción de la cobertura de bosque, sus subtipos o su biomasa. La tesis se desarrolló en el Instituto Nacional de Investigación y Desarrollo Forestal Marín Dracea (INCDS) y la Universidad de Alcalá (UAH) gracias a varios proyectos: el proyecto EO-ROFORMON del INCDS (Prototyping an Earth-Observation based monitoring and forecasting system for the Romanian forests), y el proyecto EMAFOR de la UAH (Synthetic Aperture Radar (SAR) enabled Analysis Ready Data (ARD) cubes for efficient monitoring of agricultural and forested landscapes). El proyecto EO-ROFORMON fue financiado por la Autoridad Nacional para la Investigación Científica de Rumania y el Fondo Europeo de Desarrollo Regional. El proyecto EMAFOR fue financiado por la Comunidad Autónoma de Madrid (España). El objetivo de esta tesis es el desarrollo de algoritmos para la extracción de variables forestales de uso general como la cobertura, el tipo o la biomasa del bosque a partir de imagen de radar de apertura sintética. Para alcanzar dicho propósito se analizaron posibles fuentes de sesgo sistemático que podrían aparecer en zonas de montaña (ej., normalización topográfica, generación de mosaicos), y se aplicaron técnicas de aprendizaje de máquina para tareas de clasificación y regresión. La tesis contiene ocho secciones: una introducción, cinco publicaciones en revistas o actas de congresos indexados, una pendiente de publicación (quinto capítulo) y las conclusiones. La introducción contextualiza la importancia del bosque, cómo se recoge la información sobre su estado (ej., inventario forestal) y las iniciativas o marcos legislativos que requieren dicha información. A continuación, se describe cómo la teledetección puede complementar la información de inventario forestal, detallando el contexto histórico de las distintas tecnologías, su funcionamiento, y cómo pueden ser aplicadas para la extracción de información forestal. Por último, se describe la problemática y el monitoreo del bosque en Rumanía, detallando el objetivo de la tesis y su estructura. El primer capítulo analiza la influencia del modelo digital de elevaciones (MDE) en la calidad de la normalización topográfica, analizando tres MDE globales (SRTM, AW3D y TanDEM-X DEM) y uno nacional (PNOA-LiDAR). Los experimentos se basan en la comparación entre órbitas, con un MDE de referencia, y la variación del acierto en la clasificación dependiendo del MDE empleado para la normalización. Los resultados muestran una menor diferencia ente órbitas al utilizar un MDE con una mejor resolución (ej. TanDEM-X, PNOA-LIDAR), especialmente en el caso de zonas con fuertes pendientes o formas del terreno complejas, como pueden ser los valles. En zonas de alta montaña las imágenes de radar de apertura sintética (SAR) sufren frecuentes distorsiones. Estas distorsiones dependen de la geometría de adquisición, por lo que es posible combinar imágenes adquiridas desde varias órbitas para que la cobertura sea lo más completa posible. El segundo capítulo evalúa dos metodologías para la clasificación de usos del suelo utilizando datos de Sentinel-1 adquiridos desde varias órbitas. El primer método crea clasificaciones por órbita y las combina, mientras que el segundo genera un mosaico con datos de múltiples órbitas y lo clasifica. El acierto obtenido mediante combinación de clasificaciones es ligeramente mayor, mientras que la clasificación de mosaicos tiene importantes omisiones de las zonas boscosas debido a problemas en la normalización topográfica y a los efectos direccionales. El tercer capítulo se enfoca en separar la cobertura forestal de otras coberturas del suelo (urbano, vegetación baja, agua) analizando la utilidad de las variables basadas en la coherencia interferométrica. En él se realizan tres clasificaciones de máquina vector-soporte basadas en un conjunto concreto de variables. El primer conjunto contiene las estadísticas anuales de la retrodispersión (media y desviación típica anual), el segundo añade la coherencia a largo plazo (separación temporal mayor a un año), el tercero incluye las estadísticas de la coherencia a corto plazo (mínima separación temporal). Utilizar variables basadas en la coherencia aumenta el acierto de la clasificación hasta un 5% y reduce los errores de omisión de la cobertura forestal. El cuarto capítulo evalúa la posibilidad de detectar talas selectivas utilizando datos de Sentinel-1 y Sentinel-2. Sus resultados muestran que la detección resulta muy difícil debido a la saturación de los sensores y la confusión introducida por el efecto de la fenología. El quinto capítulo se centra en la clasificación de tipos de bosque basado en una serie temporal de datos Sentinel-1. Se basa en la creación de un conjunto de modelos que describen la relación entre la retrodispersión y el ángulo local de incidencia para un determinado tipo de bosque y fecha concreta. Para cada píxel se calcula el residuo respecto al modelo de cada uno de los tipos de bosque, acumulando dichos residuos a lo largo de la serie temporal. Hecho esto, cada píxel es asignado al tipo de bosque que acumula un menor residuo. Los resultados son prometedores, mostrando que frondosas y coníferas tienen un comportamiento distintivo, y que es posible separar ambos tipos de bosque con un alto grado de acierto. El sexto capítulo está dedicado a la estimación de biomasa utilizando datos Sentinel-1, ALOS PALSAR y regresión Random Forest. Se obtiene un error similar para ambos sensores a pesar de utilizar una banda diferente (band-C vs. -L), con poca reducción en el error cuando ambas bandas se utilizan conjuntamente. Sin embargo, el ajuste de un estimador adaptado a las condiciones locales de Rumanía sí ofreció una reducción de del error al ser comparado con las estimaciones globales de biomasa