Sentinel-1 InSAR coherence to detect floodwater in urban areas: Houston and hurricane harvey as a test case

This paper presents an automatic algorithm for mapping floods. Its main characteristic is that it can detect not only inundated bare soils, but also floodwater in urban areas. The synthetic aperture radar (SAR) observations of the flood that hit the city of Houston (Texas) following the landfall of Hurricane Harvey in 2017 are used to apply and validate the algorithm. The latter consists of a two-step approach that first uses the SAR data to identify buildings and then takes advantage of the Interferometric SAR coherence feature to detect the presence of floodwater in urbanized areas. The preliminary detection of buildings is a pre-requisite for focusing the analysis on the most risk-prone areas. Data provided by the Sentinel-1 mission acquired in both Strip Map and Interferometric Wide Swath modes were used, with a geometric resolution of 5 m and 20 m, respectively. Furthermore, the coherence-based algorithm takes full advantage of the Sentinel-1 mission's six-day repeat cycle, thereby providing an unprecedented possibility to develop an automatic, high-frequency algorithm for detecting floodwater in urban areas. The results for the Houston case study have been qualitatively evaluated through very-high-resolution optical images acquired almost simultaneously with SAR, crowdsourcing points derived by photointerpretation from Digital Globe and Federal Emergency Management Agency's (FEMA) inundation model over the area. For the first time the comparison with independent data shows that the proposed approach can map flooded urban areas with high accuracy using SAR data from the Sentinel-1 satellite mission

Towards a 20m global building map from Sentinel-1 SAR Data

This study introduces a technique for automatically mapping built-up areas using synthetic aperture radar (SAR) backscattering intensity and interferometric multi-temporal coherence generated from Sentinel-1 data in the framework of the Copernicus program. The underlying hypothesis is that, in SAR images, built-up areas exhibit very high backscattering values that are coherent in time. Several particular characteristics of the Sentinel-1 satellite mission are put to good use, such as its high revisit time, the availability of dual-polarized data, and its small orbital tube. The newly developed algorithm is based on an adaptive parametric thresholding that first identifies pixels with high backscattering values in both VV and VH polarimetric channels. The interferometric SAR coherence is then used to reduce false alarms. These are caused by land cover classes (other than buildings) that are characterized by high backscattering values that are not coherent in time (e.g., certain types of vegetated areas). The algorithm was tested on Sentinel-1 Interferometric Wide Swath data from five different test sites located in semiarid and arid regions in the Mediterranean region and Northern Africa. The resulting building maps were compared with the Global Urban Footprint (GUF) derived from the TerraSAR-X mission data and, on average, a 92% agreement was obtained.Peer ReviewedPostprint (published version

Flood mapping in vegetated areas using an unsupervised clustering approach on Sentinel-1 and-2 imagery

The European Space Agency's Sentinel-1 constellation provides timely and freely available dual-polarized C-band Synthetic Aperture Radar (SAR) imagery. The launch of these and other SAR sensors has boosted the field of SAR-based flood mapping. However, flood mapping in vegetated areas remains a topic under investigation, as backscatter is the result of a complex mixture of backscattering mechanisms and strongly depends on the wave and vegetation characteristics. In this paper, we present an unsupervised object-based clustering framework capable of mapping flooding in the presence and absence of flooded vegetation based on freely and globally available data only. Based on a SAR image pair, the region of interest is segmented into objects, which are converted to a SAR-optical feature space and clustered using K-means. These clusters are then classified based on automatically determined thresholds, and the resulting classification is refined by means of several region growing post-processing steps. The final outcome discriminates between dry land, permanent water, open flooding, and flooded vegetation. Forested areas, which might hide flooding, are indicated as well. The framework is presented based on four case studies, of which two contain flooded vegetation. For the optimal parameter combination, three-class F1 scores between 0.76 and 0.91 are obtained depending on the case, and the pixel- and object-based thresholding benchmarks are outperformed. Furthermore, this framework allows an easy integration of additional data sources when these become available

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

A Human-Centered Framework for the Understanding of Synthetic Aperture Radar Images

The limited usage of SAR data in the end-user community and in applicative contexts testified the failure of the recent literature, in which the research privileged the automatic extraction of information at the expense of users' experience with data. The development of new products and processing frameworks providing user-friendly representations and extraction of the physical information is a necessary condition for the full exploitation of SAR sensors. In this Book, the necessity to restore users’ centrality in remote sensing data analysis is analyzed and achieved through the introduction of two new classes of RGB SAR products obtained via multitemporal processing, whose principal characteristics are the ease of interpretation and the possibility to be processed with simple, end-user-oriented technique. These proposed approach aims to definitely fill the gap between the academy and the applications. The rationale is to provide ready-to-use images, in which the technical expertise with electromagnetic models, SAR imaging and image processing has been absorbed in the products formation phase. In such way, the idea that SAR images are too complicated to be interpreted and processed without a high technical expertise in order to extract physical information is overcame

3D space intersection features extraction from Synthetic Aperture Radar images

The main purpose of this Thesis is to develop new theoretical models in order to extend the capabilities of SAR images space intersection techniques to generate three dimensional information. Furthermore, the study aims at acquiring new knowledge on SAR image interpretation through the three dimensional comprehension of the scene. The proposed methodologies allow to extend the known radargrammetric applications to vector data generation, exploiting SAR images acquired with every possible geometries. The considered geometries are points, circles, cylinders and lines. The study assesses the estimation accuracy of the features in terms of absolute and relative position and dimensions, analyzing the nowadays operational SAR sensors with a special focus on the national COSMO-SkyMed system. The proposed approach is original as it does not require the direct matching between homologous points of different images, which is a necessary step for the classical radargrammetric techniques; points belonging to the same feature, circular or linear, recognized in different images, are matched through specific models in order to estimate the dimensions and the location of the feature itself. This approach is robust with respect to the variation of the viewing angle of the input images and allows to better exploit archive data, acquired with diverse viewing geometries. The obtained results confirm the validity of the proposed theoretical approach and enable important applicative developments, especially in the Defence domain: (i) introducing original three dimensional measurement tools to support visual image interpretation; (ii) performing an advanced modelling of building counting only on SAR images; (iii) exploiting SAR images as a source for geospatial information and data; (iv) producing geospatial reference information, such as Ground Control Point, without any need for survey on the ground

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

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