ALOS-2/PALSAR-2 Calibration, Validation, Science and Applications

Twelve edited original papers on the latest and state-of-art results of topics ranging from calibration, validation, and science to a wide range of applications using ALOS-2/PALSAR-2. We hope you will find them useful for your future research

Disaster debris estimation using high-resolution polarimetric stereo-SAR

AbstractThis paper addresses the problem of debris estimation which is one of the most important initial challenges in the wake of a disaster like the Great East Japan Earthquake and Tsunami. Reasonable estimates of the debris have to be made available to decision makers as quickly as possible. Current approaches to obtain this information are far from being optimal as they usually rely on manual interpretation of optical imagery. We have developed a novel approach for the estimation of tsunami debris pile heights and volumes for improved emergency response. The method is based on a stereo-synthetic aperture radar (stereo-SAR) approach for very high-resolution polarimetric SAR. An advanced gradient-based optical-flow estimation technique is applied for optimal image coregistration of the low-coherence non-interferometric data resulting from the illumination from opposite directions and in different polarizations. By applying model based decomposition of the coherency matrix, only the odd bounce scattering contributions are used to optimize echo time computation. The method exclusively considers the relative height differences from the top of the piles to their base to achieve a very fine resolution in height estimation. To define the base, a reference point on non-debris-covered ground surface is located adjacent to the debris pile targets by exploiting the polarimetric scattering information. The proposed technique is validated using in situ data of real tsunami debris taken on a temporary debris management site in the tsunami affected area near Sendai city, Japan. The estimated height error is smaller than 0.6m RMSE. The good quality of derived pile heights allows for a voxel-based estimation of debris volumes with a RMSE of 1099m3. Advantages of the proposed method are fast computation time, and robust height and volume estimation of debris piles without the need for pre-event data or auxiliary information like DEM, topographic maps or GCPs

Polarimetric Synthetic Aperture Radar

This open access book focuses on the practical application of electromagnetic polarimetry principles in Earth remote sensing with an educational purpose. In the last decade, the operations from fully polarimetric synthetic aperture radar such as the Japanese ALOS/PalSAR, the Canadian Radarsat-2 and the German TerraSAR-X and their easy data access for scientific use have developed further the research and data applications at L,C and X band. As a consequence, the wider distribution of polarimetric data sets across the remote sensing community boosted activity and development in polarimetric SAR applications, also in view of future missions. Numerous experiments with real data from spaceborne platforms are shown, with the aim of giving an up-to-date and complete treatment of the unique benefits of fully polarimetric synthetic aperture radar data in five different domains: forest, agriculture, cryosphere, urban and oceans

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

Temporal Characteristics of Boreal Forest Radar Measurements

Radar observations of forests are sensitive to seasonal changes, meteorological variables and variations in soil and tree water content. These phenomena cause temporal variations in radar measurements, limiting the accuracy of tree height and biomass estimates using radar data. The temporal characteristics of radar measurements of forests, especially boreal forests, are not well understood. To fill this knowledge gap, a tower-based radar experiment was established for studying temporal variations in radar measurements of a boreal forest site in southern Sweden. The work in this thesis involves the design and implementation of the experiment and the analysis of data acquired. The instrument allowed radar signatures from the forest to be monitored over timescales ranging from less than a second to years. A purpose-built, 50 m high tower was equipped with 30 antennas for tomographic imaging at microwave frequencies of P-band (420-450 MHz), L-band (1240-1375 MHz) and C-band (5250-5570 MHz) for multiple polarisation combinations. Parallel measurements using a 20-port vector network analyser resulted in significantly shorter measurement times and better tomographic image quality than previous tower-based radars. A new method was developed for suppressing mutual antenna coupling without affecting the range resolution. Algorithms were developed for compensating for phase errors using an array radar and for correcting for pixel-variant impulse responses in tomographic images. Time series results showed large freeze/thaw backscatter variations due to freezing moisture in trees. P-band canopy backscatter variations of up to 10 dB occurred near instantaneously as the air temperature crossed 0⁰C, with ground backscatter responding over longer timescales. During nonfrozen conditions, the canopy backscatter was very stable with time. Evidence of backscatter variations due to tree water content were observed during hot summer periods only. A high vapour pressure deficit and strong winds increased the rate of transpiration fast enough to reduce the tree water content, which was visible as 0.5-2 dB backscatter drops during the day. Ground backscatter for cross-polarised observations increased during strong winds due to bending tree stems. Significant temporal decorrelation was only seen at P-band during freezing, thawing and strong winds. Suitable conditions for repeat-pass L-band interferometry were only seen during the summer. C-band temporal coherence was high over timescales of seconds and occasionally for several hours for night-time observations during the summer. Decorrelation coinciding with high transpiration rates was observed at L- and C-band, suggesting sensitivity to tree water dynamics.The observations from this experiment are important for understanding, modelling and mitigating temporal variations in radar observables in forest parameter estimation algorithms. The results also are also useful in the design of spaceborne synthetic aperture radar missions with interferometric and tomographic capabilities. The results motivate the implementation of single-pass interferometric synthetic aperture radars for forest applications at P-, L- and C-band

Study of the speckle noise effects over the eigen decomposition of polarimetric SAR data: a review

This paper is focused on considering the effects of speckle noise on the eigen decomposition of the co- herency matrix. Based on a perturbation analysis of the matrix, it is possible to obtain an analytical expression for the mean value of the eigenvalues and the eigenvectors, as well as for the Entropy, the Anisotroopy and the dif- ferent a angles. The analytical expressions are compared against simulated polarimetric SAR data, demonstrating the correctness of the different expressions.Peer ReviewedPostprint (published version

Advanced pixel selection and optimization algorithms for Persistent Scatterer Interferometry (PSI)

Tesi amb diferents seccions retallades per dret de l'editorPremi Extraordinari de Doctorat, promoció 2018-2019. Àmbit de les TICGround deformation measurements can provide valuable information for minimization of associated loss and damage caused by natural and environmental hazards. As a kind of remote sensing technique, Persistent Scatterer Interferometry (PSI) SAR is able to measure ground deformation with high spatial resolution, efficiently. Moreover, the ground deformation monitoring accuracy of PSI techniques can reach up to millimeter level. However, low coherence could hinderthe exploitation of SAR data, and high-accuracy deformation monitoring can only be achieved by PSI for high quality pixels. Therefore, pixel optimization and identification of coherent pixels are crucial for PSI techniques. In this thesis, advanced pixel selection and optimization algorithms have been investigated. Firstly, a full-resolution pixel selection method based on the Temporal Phase Coherence (TPC) has been proposed. This method first estimates noise phase term of each pixel at interferogram level. Then, for each pixel, its noise phase terms of all interferograms are used to assess this pixel’s temporal phase quality (i.e., TPC). In the next, based on the relationship between TPC and phase Standard Deviation (STD), a threshold can be posed on TPC to identify high phase quality pixels. This pixel selection method can work with both Deterministic Scatterers (PSs) and Distributed Scatterers (DSs). To valid the effectiveness of the developed method, it has been used to monitor the Canillo (Andorra) landslide. The results show that the TPC method can obtained highest density of valid pixels among the employed three approaches in this challenging area with X-band SAR data. Second, to balance the polarimetric DInSAR phase optimization effect and the computation cost, a new PolPSI algorithm is developed. This proposed PolPSI algorithm is based on the Coherency Matrix Decomposition result to determine the optimal scattering mechanism of each pixel, thus it is named as CMD-PolPSI. CMDPolPSI need not to search for solution within the full space of solution, it is therefore much computationally faster than the classical Equal Scattering Mechanism (ESM) method, but with lower optimization performance. On the other hand, its optimization performance outperforms the less computational costly BEST method. Third, an adaptive algorithm SMF-POLOPT has been proposed to adaptive filtering and optimizing PolSAR pixels for PolPSI applications. This proposed algorithm is based on PolSAR classification results to firstly identify Polarimetric Homogeneous Pixels (PHPs) for each pixel, and at the same time classify PS and DS pixels. After that, DS pixels are filtered by their associated PHPs, and then optimized based on the coherence stability phase quality metric; PS pixels are unfiltered and directly optimized based on the DA phase quality metric. SMF-POLOPT can simultaneously reduce speckle noise and retain structures’ details. Meanwhile, SMF-POLOPT is able to obtain much higher density of valid pixels for deformation monitoring than the ESM method. To conclude, one pixel selection method has been developed and tested, two PolPSI algorithms have been proposed in this thesis. This work make contributions to the research of “Advanced Pixel Selection and Optimization Algorithms for Persistent Scatterer InterferometryLes mesures de deformació del sòl poden proporcionar informació valuosa per minimitzar les pèrdues i els danys associats causats pels riscos naturals i ambientals. Com a tècnica de teledetecció, la interferometria de dispersors persistents (Persistent Scatter Interferometry, PSI) SAR és capaç de mesurar de forma eficient la deformació del terreny amb una alta resolució espacial. A més, la precisió de monitorització de la deformació del sòl de les tècniques PSI pot arribar a arribar a nivells del mil·límetre. No obstant això, una baixa coherència pot dificultar l’explotació de dades SAR i el control de deformació d’alta precisió només es pot aconseguir mitjançant PSI per a píxels d’alta qualitat. Per tant, l’optimització de píxels i la identificació de píxels coherents són crucials en les tècniques PSI. En aquesta tesi s¿han investigat algorismes avançats de selecció i optimització de píxels. En primer lloc, s'ha proposat un mètode de selecció de píxels de resolució completa basat en la coherència temporal de fase (Temporal Phase Coherence, TPC). Aquest mètode estima per primera vegada el terme de fase de soroll de cada píxel a nivell d’interferograma. A continuació, per a cada píxel, s'utilitzen els termes de la fase de soroll de tots els interferogrames per avaluar la qualitat de fase temporal d'aquest píxel (és a dir, TPC). A la següent, basant-se en la relació entre el TPC i la desviació estàndard de fase (STD), es pot plantejar un llindar de TPC per identificar píxels de qualitat de fase alta. Aquest mètode de selecció de píxels es capaç de detectar tant els dispersors deterministes (PS) com els distribuïts (DS). Per validar l’eficàcia del mètode desenvolupat, s’ha utilitzat per controlar l’esllavissada de Canillo (Andorra). Els resultats mostren que el mètode TPC pot obtenir la major densitat de píxels vàlids, comparat amb els mètodes clàssics de selecció, en aquesta àrea difícil amb dades de SAR de banda X. En segon lloc, per equilibrar l’efecte d’optimització de fase DInSAR polarimètrica i el cost de càlcul, es desenvolupa un nou algorisme de PolPSI. Aquest algorisme proposat de PolPSI es basa en el resultat de la descomposició de la matriu de coherència per determinar el mecanisme de dispersió òptim de cada píxel, de manera que es denomina CMD-PolPSI. CMDPolPSI no necessita buscar solucions dins de l’espai complet de la solució, per tant, és molt més eficient computacionalment que el mètode clàssic de mecanismes d’igualtat de dispersió (Equal Scattering Mechanism, ESM), però amb un efecte d’optimització no tant òptim. D'altra banda, el seu efecte d'optimització supera el mètode BEST, el que te un menor cost computacional. En tercer lloc, s'ha proposat un algoritme adaptatiu SMF-POLOPT per al filtratge adaptatiu i l'optimització de píxels PolSAR per a aplicacions PolPSI. Aquest algorisme proposat es basa en els resultats de classificació PolSAR per identificar primer els píxels homogenis polarimètrics (PHP) per a cada píxel i, alhora, classificar els píxels PS i DS. Després d'això, els píxels DS es filtren pels seus PHP associats i, a continuació, s'optimitzen en funció de la mètrica de qualitat de la fase d'estabilitat de coherència; els píxels classificats com PS no es filtren i s'optimitzen directament en funció de la mètrica de qualitat de la fase DA. SMF-POLOPT pot reduir simultàniament el soroll de la fase interferomètrica i conservar els detalls de les estructures. Mentrestant, SMF-POLOPT aconsegueix obtenir una densitat molt més alta de píxels vàlids per al seguiment de la deformació que el mètode ESM. Per concloure, en aquesta tesi s’ha desenvolupat i provat un mètode de selecció de píxels, i s’han proposat dos algoritmes PolPSI. Aquest treball contribueix a la recerca en "Advanced Pixel Selection and Optimization Algorithms for Persistent Scatterer Interferometry"Postprint (published version

Cryosphere Applications

Synthetic aperture radar (SAR) provides large coverage and high resolution, and it has been proven to be sensitive to both surface and near-surface features related to accumulation, ablation, and metamorphism of snow and firn. Exploiting this sensitivity, SAR polarimetry and polarimetric interferometry found application to land ice for instance for the estimation of wave extinction (which relates to sub surface ice volume structure) and for the estimation of snow water equivalent (which relates to snow density and depth). After presenting these applications, the Chapter proceeds by reviewing applications of SAR polarimetry to sea ice for the classification of different ice types, the estimation of thickness, and the characterisation of its surface. Finally, an application to the characterisation of permafrost regions is considered. For each application, the used (model-based) decomposition and polarimetric parameters are critically described, and real data results from relevant airborne campaigns and space borne acquisitions are reported