A multisensor approach for the 2016 Amatrice earthquake damage assessment ANNALS OF GEOPHYSICS, 59, Fast Track 5, 2016; DOI: 10.4401/ag-7185 1 A multisensor approach for the 2016 Amatrice earthquake damage as- sessment VITO ROMANIELLO*1, ALESSANDRO PISCINI1, CHRISTIAN BIGNAMI1, ROBERTA ANNIBALLE2, AND SALVATORE STRAMONDO1 1 Istituto Nazionale di Geofisica e Vulcanologia 2 DIET, Sapienza University of Rome *vito.romaniello@ingv.it Abstract This work proposes methodologies aimed at evaluating the damage occurred in the Amatrice town by us- ing optical and Synthetic Aperture Radar (SAR) change features obtained from satellite images. The ob- jective is to achieve a damage map employing the satellite change features in a classifier algorithm, namely the Features Stepwise Thresholding (FST) method. The main novelties of the proposed analysis concern the estimation of derived features at object scale and the exploitation of the unsupervised FST algorithm. A segmentation of the study area into several buildings blocks has been done by considering a set of poly- gons, over the Amatrice town, extracted from the open source Open Street Map (OSM) geo-database. The available satellite dataset is composed of several optical and SAR images, collected before and after the seismic event. Regarding the optical data, we selected the Normalised Difference Index (NDI), and two quantities com- ing from the Information Theory, namely the Kullback-Libler Divergence (KLD) and the Mutual Infor- mation (MI). In addition, for the SAR data we picked out the Intensity Correlation Difference (ICD) and the KLD parameter. The exploitation of these features in the FST algorithm permits to obtain a plausible damage map that is able to indicate the most affected areas. I. INTRODUCTION n earthquake damage map, available right after a seismic event (from few hours up to few days, depending on the satellite data availability), can guide the rescue teams interventions towards the most affected areas. Satellite data can be very useful for this purpose, thanks to the wide coverage and the high spatial resolution, especially to map dam- age in large regions, where villages cannot eas- ily accessed. Satellite damage assessment is based on detection techniques that are capable to identify changes of an object by observing it at different times (Ingram (1981); Jenson (1983)). Both optical and radar sensors can be exploited for change detection purposes and several examples are reported in the literature A ANNALS OF GEOPHYSICS, 59, Fast Track 5, 2016; DOI: 10.4401/ag-7185 2 (Matsuoka and Yamazaki (2004); Hoffman et al. (2004); Yonezawa and Takeuchi (2001); Stramondo et al. (2006); Chini et al. (2013)). By exploiting both optical and radar sen- sors, and using change features achieved from the statistical analysis, a more accurate and re- liable damage mapping can be obtained. Big- nami et al. (2004) have analysed the possible advantages of combining radar and optical satellite data. Regarding the SAR data, the Intensity Correlation Difference (ICD) demonstrated to be a very good damage proxy (Stramondo et al. (2006)). KLD and MI (Brunner et al. (2010)) are also suitable features that can contribute to damage estimation. As far as the optical data is concerned, the most significative results are related to the Normalised Difference Index (NDI), Kullback- Libler Divergence (KLD) (Kullback and Leiber (1951)), and Mutual Information (MI) indexes (Xie et al. (2003)). All these features show a good sensitivity to the collapse ratio. In this work, we used the NDI, KLD and MI indexes for the optical data, and the ICD and the KLD parameters for the SAR data, fol- lowing the outcomes carried out during EU funded project APhoRISM (Advanced PRoce- dures for volcanic and Seismic Monitoring, www.aphorism-project.eu). The main novelties of the proposed anal- ysis concern the estimation of derived features at object scale and the use of an unsupervised algorithm for damage estimation purposes. In order to exploit the object-oriented approach, a segmentation of the study area into several re- gions was performed. In particular, testing an operational approach, we have segmented and then generated the damage map of Amatrice by considering a set of polygons extracted from the open source Open Street Map (OSM) geo-database. Finally, an unsupervised algo- rithm to estimate the damage level by satellite features, namely in the paper Features Step- wise Thresholding (FST), has been used (Roan- iello et al. (2016)). II. DATA AND METHODS The available dataset is made up of two Senti- nel-2 optical images (1 pre- and 1 post- seismic), three COSMO-Sky SAR images (2 pre- and 1 post-seismic) and a buildings foot- print layer extracted by the Open Street Map service (see Table 1 for EO data characteris- tics). Both optical and SAR data have a spatial resolution of 10 m. Datatype Satellite Acquisition Optical Sentinel-2 2016/08/14 Optical Sentinel-2 2016/09/04 SAR COSMO-Sky 2016/07/19 SAR COSMO-Sky 2016/08/20 SAR COSMO-Sky 2016/08/28 Table 1: EO data list The Amatrice footprint layer extracted from OSM, which is at single building scale, was modified to obtain polygons surrounding more than one building. In this way, there are more pixels associated to each polygon leading a better estimation of change features over the polygon itself. The resulting layer consist of 112 polygons, inside an area of about 1 km2. Previous studies, carried out during the European APhoRISM project, allowed us to identify the features having a good sensitivity to damage at object scale (Romaniello et al. (2016)). Following the project outcomes, in the present work we used the NDI, and two quan- tities coming from the Information Theory, the KLD and the MI. In addition, for the SAR data, we picked out the ICD parameter. ANNALS OF GEOPHYSICS, 59, Fast Track 5, 2016; DOI: 10.4401/ag-7185 3 The NDI parameter is defined as: 𝑁𝐷𝐼𝑖 = π‘ƒπ‘‚π‘†π‘‡π‘–βˆ’π‘ƒπ‘…πΈπ‘– 𝑃𝑂𝑆𝑇𝑖+𝑃𝑅𝐸𝑖 (1) Where PREi and POSTi indicate the mean values of intensity, respectively for pre- and post- seismic images, associated to i-th poly- gon. The NDI parameter was computed, con- sidering pre- and post-seismic Sentinel-2 im- ages, for each polygon obtained from OSM service. Fig. 1 shows the resulting NDI map. Figure 1: NDI map over Amatrice town obtained from pre- and post-seismic Sentinel-2 image Actual NDI values range from -0.080 and 0.311, with the highest values centred on the historical centre of Amatrice town, where most severe damage were registered (Azzaro et al. 2016). The KLD parameter is defined as: 𝐾𝐿𝐷𝑖 = (π‘ƒπ‘…πΈπ‘–βˆ’π‘ƒπ‘‚π‘†π‘‡π‘–) 2+π‘‰π‘Žπ‘Ÿ(𝑃𝑅𝐸𝑖)+π‘‰π‘Žπ‘Ÿ(𝑃𝑂𝑆𝑇𝑖) 2 Γ— ( 1 π‘‰π‘Žπ‘Ÿ(𝑃𝑅𝐸𝑖) + 1 π‘‰π‘Žπ‘Ÿ(𝑃𝑂𝑆𝑇𝑖) ) βˆ’ 2 (2) Where PREi and POSTi are the same pa- rameters in the Eq. 1, and Var(PREi) and Var(POSTi) are their variances within the i-th polygon. Nominally, the KLD parameter has the same behaviour of NDI: KLD increasing values correspond to increasing damage level. The KLD values, calculated from the two Sen- tinel-2 optical images, range from 0.0 and 8.5. The MI index measures the correlation be- tween pre- and post-seismic images (see Equa- tion 3). 𝑀𝐼𝑖 = βˆ’π‘™π‘›(1 βˆ’ π‘Ÿπ‘– 2) βˆ— 0.5 (3) Where ri2 is the correlation between the pre- and post-seismic pixels within each poly- gon. MI is related to the spatial arrangement of the pixels, and it is inversely proportional to the damage grade. For the present case study, its values range from 0.0 to 3.0. As far as the SAR data is concerned, the KLD feature is derived in the same way as the optical data, while the ICD can be computed on the base of the Pearson Correlation coeffi- cient (ρi) estimated on the pre-seismic SAR im- age pair ICpre and on the co-seismic SAR im- age pair ICcos. These two intermediate output are then used to obtain the ICD: 𝐼𝐢𝐷𝑖 = πΌπΆπ‘π‘Ÿπ‘’ βˆ’ πΌπΆπ‘π‘œπ‘  (4) Each feature provides a qualitative as- sessment of the damage distribution within the town. In order to achieve a damage map useful in the rescue operations, the unsupervised al- gorithm FST (Romaniello et al. (2016)) and based on these satellite features, has been em- ployed; the scheme of the algorithm is shown in the Fig. 2. The algorithm is able to use all change features, both from optical and SAR data, analyzing the distribution of values of each parameter. The entire ranges of values of each parameter, obtained excluding the outliers, are divided into three intervals corresponding to three damage levels. This classification is repeated n-times (in the present analysis, n is ANNALS OF GEOPHYSICS, 59, Fast Track 5, 2016; DOI: 10.4401/ag-7185 4 equal to 21) varying the width of the three intervals at each iteration. Figure 2: FST method block scheme So, we obtain several damage level estimations for each polygon (specifically, 21 iterations x 5 features). After this operation, in order to assign the final damage level (low, medium and high) to a polygon, the modal value of all previous estimations is considered for that polygon. The output is a map of the damage level occurred in Amatrice. It is worth noticing that the damage clas- sification (low, medium and high) is based on a relative scale. Romaniello et al. (2016) have shown the correlation between this kind of damage classification and the collapse ratio, defined as the ratio between the number of col- lapsed buildings and the total number of buildings within a polygon. Therefore, high damage can indicate many collapses within a city block. III. RESULTS A damage-retrieved map can be obtained by using both optical and SAR satellite images as input for the FST method. The data fusion ap- proach showed the best damage assessment. The resulting damage map is reported in the Fig. 3. From the classification, 86 polygons are associated to the lowest damage grade, 12 to the medium and 14 to the highest. Figure 3: Damage classification by using all change fea- tures (both optical and SAR) over Amatrice town The most affected area is the historical centre of the town and in particular the West side as visible in Fig. 4. Figure 4: Damage classification in the historical centre. ANNALS OF GEOPHYSICS, 59, Fast Track 5, 2016; DOI: 10.4401/ag-7185 5 The data fusion approach has overcome false alarms problems occurred in the single sensor approach for several buildings. In particular, the analysis performed using only optical features has led to false positive cases recognized by the comparison with a very high resolution VHR image acquired after the seismic event (Fig. 5). Figure 5: Amatrice post- earthquake damage. In the classification maps reported in Fig. 3 and Fig. 4 the most affected area is the West side of the historical center with several col- lapsed buildings. This confirms the soundness of the FST method. Nevertheless, a quantita- tive assessment will be performed in the next future to draw much more detailed conclusion on the accuracy of the derived products. IV. DISCUSSION The proposed method is a quite robust ap- proach being based on unsupervised classifica- tion and employing several change features, both from optical and SAR data. Furthermore, this method, being based on city blocks areas, is reliable with respect to changes not related to earthquake, because we can guess that modifications caused by tem- poral changes are small (within the blocks) if compared with the changes due to such cata- strophic event. In other words, these sources of errors are averaged within the blocks. Moreover, it is worth noticing that the proposed damage classification is on a relative scale. Indeed, the method identifies always three classes of damage associated with the low, medium and high grades of changes for the specific regions we are analyzing. On the one hand, this characteristic can be considered as a limitation of the algorithm, but from an operational point of view, we think that during rescue activities it is important to know where to intervene promptly, independently from the absolute value of damage. Finally, a comparison with the macro- seismic survey will be done, as soon as a com- plete survey will be available. ACKNOWLEDGEMENTS The research leading to these results has re- ceived funding from the European Union’s Seventh Framework Programme (FP7/2007- 2013) under grant agreement number 606738, APhoRISM project. REFERENCES Azzaro, R., and the QUEST Working Group (2016). "The 24 August 2016 Amatrice earth- quake: macroseismic survey in the damage ar- ea and EMS intensity assessment", Ann. Ge- ophys., vol. 59, no. 5. Bignami, C., Chini, M., Pierdicca, N. and Stra- mondo, S. (2004). Comparing and combining the capability of detecting earthquake damage in urban areas using SAR and optical data. Proceedings of IEEE IGARSS, 1, pp. 55–58. Brunner, D., Lemoine, G. and Bruzzone, L. (2010). "Earthquake damage assessment of buildings using VHR optical and SAR image- ANNALS OF GEOPHYSICS, 59, Fast Track 5, 2016; DOI: 10.4401/ag-7185 6 ry", IEEE Trans. Geosci. Remote Sens., vol. 48, no. 5, pp. 2403-2420. Chini, M., Piscini, A., Cinti, F.R., Amici, S., Nappi, R., De Martini, P.M. (2013). The 2011 Tohoku (Japan) Tsunami inundation and liq- uefaction investigated through optical, ther- mal, and SAR data. IEEE Trans. Geosci. Remote Sens., 10, 347–351. Hoffman, J., Roth, A. and Voigt, S. (2004). In- terferometric displacement and damage as- sessment for the 2003 Bam earthquake. In Pro- ceedings of the 2004 Envisat and ERS Symposi- um, Salzburg, Austria, 6–10 September 2004, ESA SP-572. Ingram K., Knapp E., and Robinson J. W. (1981). Change detection technique develop- ment for improved urbanized area delineation, technical memorandum CSCITM-81/6087, Computer Sciences Corporation, Silver Springs, Maryland, U.S.A. Jenson J. R. (editor) (1983). Urban/suburban land use analysis. Manual of Remote Sensing, Vol. 2, second edition (Falls Church, Virginia: American Society of Photogrammetry), pp: 1571-1666. Kullback, S.; Leibler, R. A. (1951). On Infor- mation and Sufficiency. Ann. Math. Statist. 22, no. 1, 79--86. doi:10.1214/aoms/1177729694. http://projecteuclid.org/euclid.aoms/1177729 694. Matsuoka, M. and Yamazaki, F. (2004). Use of satellite SAR intensity imagery for detecting building areas damaged due to earthquakes. Earthquake Spectra, 20, pp. 975–994. Romaniello, V., Piscini, A., Bignami, C., Anni- balle, R. and Stramondo, S. (2016). Earthquake damage mapping by using remotely sensed data: the Haiti case study, Journal of Applied Remote Sensing, (submitted). Stramondo, S. Bignami, C., Chini, M., Pierdic- ca, N., Tertulliani, A. (2006). Satellite radar and optical remote sensing for earthquake damage detection: results from different case studies. International Journal of Remote Sensing, 27, 4433 - 4447. H. Xie, L. E. Pierce, and F. T. Ulaby (2003). β€œMutual information based registration of SAR images,” in Proc. IEEE IGARSS, Toulouse, France, Jul. 2003, pp. 4028–4031. Yonezawa, C. and Takeuchi, S. (2001). Decorre- lation of SAR data by urban damage caused by the 1995 Hoyogoken-Nanbu earthquake. Inter- national Journal of Remote Sensing, 22, pp. 1585– 1600.