[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Next Article in Journal
Monitoring Sea Level and Topography of Coastal Lagoons Using Satellite Radar Altimetry: The Example of the Arcachon Bay in the Bay of Biscay
Previous Article in Journal
A Side Scan Sonar Image Target Detection Algorithm Based on a Neutrosophic Set and Diffusion Maps
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Remote Sensing
Volume 10
Issue 2
10.3390/rs10020296
remotesensing-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Novel Unsupervised Classification of Collapsed Buildings Using Satellite Imagery, Hazard Scenarios and Fragility Functions

by
Luis Moya
1,*,†,
Luis R. Marval Perez
2,†,
Erick Mas
1,†,
Bruno Adriano
1,†,
Shunichi Koshimura
1,† and
Fumio Yamazaki
3,†
1
International Research Institute of Disaster Science, Tohoku University, Aoba 468-I-E301, Aramaki, Aoba-ku, Sendai 980-0845, Japan
2
Graduate School of Information Science, Tohoku University, 6-6-05 Aramaki Aza Aoba, Aoba-ku, Sendai 980-8579, Japan
3
Department of Urban Environment Systems, Chiba University, 1-33 Yayoi-cho, Inage-ku, Chiba 263-8522, Japan
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Remote Sens. 2018, 10(2), 296; https://doi.org/10.3390/rs10020296
Submission received: 4 January 2018 / Revised: 30 January 2018 / Accepted: 12 February 2018 / Published: 14 February 2018
Browse Figures
Graphical abstract
"> Figure 1
<p>The functions <math display="inline"> <semantics> <mrow> <mi>g</mi> <mrow> <mo>(</mo> <msub> <mi>h</mi> <mi mathvariant="bold-italic">θ</mi> </msub> <mrow> <mo>(</mo> <mi mathvariant="bold-italic">x</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>=</mo> <mo>−</mo> <mi>ln</mi> <mrow> <mo>(</mo> <msub> <mi>h</mi> <mi mathvariant="bold-italic">θ</mi> </msub> <mrow> <mo>(</mo> <mi mathvariant="bold-italic">x</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> </mrow> </semantics> </math> and <math display="inline"> <semantics> <mrow> <mi>g</mi> <mrow> <mo>(</mo> <msub> <mi>h</mi> <mi mathvariant="bold-italic">θ</mi> </msub> <mrow> <mo>(</mo> <mi mathvariant="bold-italic">x</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <mo>=</mo> <mo>−</mo> <mi>ln</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>−</mo> <msub> <mi>h</mi> <mi mathvariant="bold-italic">θ</mi> </msub> <mrow> <mo>(</mo> <mi mathvariant="bold-italic">x</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> </mrow> </semantics> </math> that contribute to the elements of the summation in Equation (<a href="#FD5-remotesensing-10-00296" class="html-disp-formula">5</a>).</p> "> Figure 2
<p>Study area. (<b>a</b>) Location of the coastal area of Miyagi Prefecture (red rectangle) within the Tohoku region of Japan; (<b>b</b>) RGB color composite of the TerraSAR-X images acquired on 13 March 2011 (red) and 21 October 2010 (green and blue); (<b>c</b>) Inundation depth map of the Great East Japan Earthquake and Tsunami. The inundation values are given in units of meters.</p> "> Figure 3
<p>(<b>a</b>) Scatter plot of the bi-dimensional dataset composed of <span class="html-italic">r</span> and <span class="html-italic">d</span> values. The colored marks denote the densities at the corresponding points; (<b>b</b>) Empirical fragility functions of buildings that collapse due to a tsunami event as proposed by Koshimura et al. [<a href="#B26-remotesensing-10-00296" class="html-bibr">26</a>] (solid line) and by Suppasri et al. [<a href="#B28-remotesensing-10-00296" class="html-bibr">28</a>] (dashed line).</p> "> Figure 4
<p>Discriminant functions obtained using the fragility function of Koshimura et al. [<a href="#B26-remotesensing-10-00296" class="html-bibr">26</a>]. (<b>a</b>–<b>g</b>) Scatter plots of the bi-dimensional dataset separated by damage state (DS), together with the obtained discriminant functions; (<b>h</b>) Variations of the cost function (<math display="inline"> <semantics> <mrow> <mi>J</mi> <mo>(</mo> <mi mathvariant="bold-italic">θ</mi> <mo>)</mo> </mrow> </semantics> </math>) throughout the iterative gradient descent algorithm.</p> "> Figure 5
<p>Discriminant functions obtained using the fragility function of Suppasri et al. [<a href="#B28-remotesensing-10-00296" class="html-bibr">28</a>]. (<b>a</b>–<b>g</b>) Scatter plots of the bi-dimensional dataset separated by DS, together with the obtained discriminant functions; (<b>h</b>) Variations of the cost function (<math display="inline"> <semantics> <mrow> <mi>J</mi> <mo>(</mo> <mi mathvariant="bold-italic">θ</mi> <mo>)</mo> </mrow> </semantics> </math>) throughout the iterative gradient descent algorithm.</p> "> Figure 6
<p>Non-linear discriminant functions calculated from the fragility functions of (<b>a</b>) Koshimura et al. [<a href="#B26-remotesensing-10-00296" class="html-bibr">26</a>] and (<b>b</b>) Suppasri et al. [<a href="#B28-remotesensing-10-00296" class="html-bibr">28</a>].</p> "> Figure 7
<p>Parallel coordinate plot of seven normalized features. Red and blue marks denote samples classified as collapsed and non-collapsed buildings, respectively. The black ticks delimit the range of <math display="inline"> <semantics> <mrow> <mo>[</mo> <mi>a</mi> <mi>v</mi> <mi>e</mi> <mi>r</mi> <mi>a</mi> <mi>g</mi> <mi>e</mi> <mo>−</mo> <mi>s</mi> <mi>t</mi> <mi>d</mi> <mo>,</mo> <mi>a</mi> <mi>v</mi> <mi>e</mi> <mi>r</mi> <mi>a</mi> <mi>g</mi> <mi>e</mi> <mo>+</mo> <mi>s</mi> <mi>t</mi> <mi>d</mi> <mo>]</mo> </mrow> </semantics> </math> for each normalized feature.</p> ">
Versions Notes

Abstract

:
Although supervised machine learning classification techniques have been successfully applied to detect collapsed buildings, there is still a major problem that few publications have addressed. The success of supervised machine learning strongly depends on the availability of training samples. Unfortunately, in the aftermath of a large-scale disaster, training samples become available only after several weeks or even months. However, following a disaster, information on the damage situation is one of the most important necessities for rapid search-and-rescue efforts and relief distribution. In this paper, a modification of the supervised machine learning classification technique called logistic regression is presented. Here, the training samples are replaced with probabilistic information, which is calculated from the spatial distribution of the hazard under consideration and one or more fragility functions. Such damage probabilities can be collected almost in real time for specific disasters such as earthquakes and/or tsunamis. We present the application of the proposed method to the 2011 Great East Japan Earthquake and Tsunami for collapsed building detection. The results show good agreement with a field survey performed by the Ministry of Land, Infrastructure, Transport and Tourism, with an overall accuracy of over 80%. Thus, the proposed method can significantly contribute to a rapid estimation of the number and locations of collapsed buildings.

Graphical Abstract">

Graphical Abstract

1. Introduction

With the recent advent of remote sensing technologies, the estimation of affected areas in the aftermath of a large-scale disaster has considerably improved [1,2,3,4,5]. Interest in the use of remote sensing in disaster management is well motivated because remote sensing products can cover a wide affected area; multiple sensors equipped on satellites are constantly recording the Earth’s surface; and active sensors function independently of the time and climate. Because of this robustness, several authors have addressed the problem of retrieving information on infrastructural damage with an acceptable accuracy.
The main challenge in damage extraction using remote sensing information is to establish a discriminant function of features extracted from satellite imagery with which to properly classify damaged and non-damaged buildings. Several features have been proposed for damage detection. Among them, the most common fundamental basis for analysis is related to the detection of changes between a pair of images acquired before and after a disaster [2,6]. However, efforts have also been made to extract features from post-event images only [7,8,9,10]. A common approach to the establishment of discriminant functions is through calibration using previous events recorded by specific sensors, such as X-band Synthetic Aperture Radar (SAR) imagery [2,11], L-band SAR imagery [12] and C-band SAR imagery [13]. However, it has been reported that such methods are subject to limitations regarding the transferability of the discriminant functions to imagery recorded by other sensors [14,15]. Furthermore, these discriminant functions are limited to a feature dataset of fixed dimensionality. Consequently, more robust and versatile calibration methods are needed. To this end, machine learning techniques have been successfully applied [16,17,18,19,20,21,22]. Calibration using supervised machine learning is performed based on training samples, that is a set of features for which the corresponding damage states are known in advance. The performance of machine learning techniques is strongly related to the number of training samples used. The more training samples there are, the better is the calibration of the discriminant function that can be achieved.
In the aftermath of a large-scale disaster, it is crucial to identify damaged buildings as early as possible to facilitate efficient search-and-rescue activities. The need for early response is related to the life expectancy of occupants trapped under collapsed buildings. For instance, after the 1995 Hyogoken-Nanbu earthquake in Japan, the life spans of the victims ranged from 5–6 min for suffocating occupants to several days for uninjured, but starving or dehydrated occupants [23]. It is therefore necessary to optimize all steps in the process of building damage detection. When machine learning techniques are employed for building damage classification, significant time is spent on retrieving training data. According to the experience of the authors, such training data are usually obtained from field surveys, which might require significant time to establish the logistics, conduct the survey and digitize the data. Another common approach is the visual interpretation of optical images, which may be influenced by the experience of the person performing the analysis. Official reports from the government provide some of the most reliable information that can be used as training data; however, such reports usually become available only after several weeks or months. To reduce the time required to retrieve training data, Wieland et al. [16] performed a quantitative evaluation of the influence of the number of training samples on the accuracy of damage classification using the support vector machine (SVM) technique, and these authors concluded that a minimum number of training samples from a small study area can be used to map building damage in a much wider region. Furthermore, remote sensing is defined as the acquisition of information about an object or phenomenon without making physical contact with that object or phenomenon. Thus, from this perspective, the fact that the leading state-of-the-art techniques for damage detection require field surveys is not completely satisfactory.
In the field of risk analysis, a probabilistic estimate of damage is expressed as a function of an engineering demand parameter (EDP) [24,25,26,27,28]. The EDP refers to a property of the hazard that is causing the damage of interest (e.g., the peak ground acceleration for an earthquake or the inundation depth for a tsunami). Unlike the use of remote sensing technologies, the application of probability theory for risk analysis emerged several decades ago [29]. Since then, outstanding progress has been achieved [30]. A fragility function is a well-known relationship between an EDP and the probability that a certain damage level is reached or exceeded [25]. In recent years, several frameworks for building damage mapping in terms of probabilities have been proposed [31,32,33,34,35]. Based on the probabilities of damage occurrence under a given hazard scenario and a geocoded building inventory, the numbers of damaged buildings within individual grid cells/blocks can be estimated. Such estimation is currently possible due to advances in hazard quantification. The spatial distributions of natural hazards, such as earthquakes and tsunamis, can be measured by device networks [36,37,38,39] and/or from theoretical models [26,40]. However, this approach provides only local aggregates, and thus, the allocation of damage states to buildings is not possible.
In this paper, we propose a new approach for the classification of collapsed buildings. This approach enables the allocation of damage states to buildings without any training samples. Therefore, it is able to overcome some of the difficulties encountered in the two approaches discussed above. Our proposed approach, called IHF, requires the following inputs: satellite imagery, the spatial distribution of the hazard and a set of fragility functions for the buildings to be evaluated. IHF is a hybrid method that combines traditional machine learning methods with probabilistic damage mapping. This manuscript presents a follow-up study to one previously published by Moya et al. [41]. That study introduced the fundamental basis of the approach. Now, after major modifications to the framework, a generalization of the approach and the development of the equations are presented. This paper is structured as follows: The next section introduces the mathematical formulations underlying the IHF method. Section 3 presents the application of the IHF method to the 2011 Great East Japan Earthquake and Tsunami. Finally, in Section 4, conclusions are drawn, and future research is discussed.

2. The IHF Classification Method

2.1. Fundamentals

The proposed method is a modification of the logistic regression method [42,43]. Here, we focus on the two-class problem, in which the two classes correspond to collapsed ( C 1 ) and non-collapsed ( C 2 ) buildings. For a given set of features D = { x i , p i } i = 1 N , where x i ∈ R m and p i ∈ [ 0 , 1 ] , p i represents the probability that x i belongs to C 1 . Thus, the probability that x i belongs to C 2 is 1 − p i . x i is a vector that contains the features used for classification. It is intended that x i will be calculated from satellite imagery, whereas p i will be calculated from the fragility functions and the hazard distribution. The main purpose is to find a discriminant function f such that any value of x i can be correctly classified as belonging to either C 1 or C 2 . A linear discriminant, where the decision surface is a hyperplane, can be represented as a linear combination of the elements of the input vector:
f ( x ) = ∑ i = 1 m θ i x i + θ 0
An input vector x is assigned to class C 1 if f ( x ) ≥ 0 and to class C 2 otherwise. The discriminant function can be vectorized as follows:
f ( x ) = θ T x
where θ = ( θ 0 , θ 1 , … , θ m ) and x is redefined as x = ( 1 , x 1 , … , x m ) . Thus, the decision boundary is defined by the relation f ( x ) = 0 . From the logistic function perspective, a threshold classifier ( h θ ( x ) ) can be expressed as follows:
h θ ( x ) = 1 1 + e − θ T x
where x belongs to C 1 if h θ ( x ) ≥ 0.5 and to C 2 otherwise. The unknown vector θ is obtained by solving the following equation:
min θ ( J ( θ ) )
J ( θ ) = − 1 N ∑ i = 1 N p i ln h θ ( x i ) + ( 1 − p i ) ln ( 1 − h θ ( x i )
Each element of the summation in Equation (5) is the contribution of one sample feature and is composed of two weighted terms of the forms ln h θ ( x i ) and ln ( 1 − h θ ( x i ) ) , where h θ ( x i ) is always a positive scalar with a value less than one. Figure 1 plots these terms when the factor of −1 in front of the summation is included. There are several numerical procedures for solving Equation (4). It has been proven elsewhere that Equation (5) is concave and hence has a unique minimum [43]. One of the simplest and best-known algorithms for this purpose is the gradient descent algorithm [42], in which the following expression is applied until convergence:
θ τ + 1 = θ τ − α ∇ J ( θ )
where τ denotes the iteration number and α is a learning rate parameter. If α is too small, the convergence may be slow, whereas an excessively large value may cause oscillations and even divergence. The gradient of the error function J ( θ ) with respect to θ is given by:
∇ J ( θ ) = ∑ n = 1 N ( h θ ( x i ) − p i ) x i

2.2. Discussion of the IHF Method

We begin by discussing the similarities and differences between the IHF method and the standard logistic regression method. To this end, let us quickly review the fundamentals of logistic regression. Consider a set of training samples Z = { z i , r i } i = 1 K , where r i = 1 if z i ∈ C 1 and r i = 0 if z i ∈ C 2 . Under the assumption that the class-conditional densities, P ( z | C i ) , are normal, our estimator of P ( C i | z ) is a logistic function [42]:
y = P ^ ( C 1 | z ) = 1 1 + e − θ T z
It is assumed that r i has a Bernoulli distribution. Then, the likelihood of θ given sample set Z is:
l ( θ | Z ) = ∏ i = 1 K ( y i ) r i ( 1 − y i ) 1 − r i
The vector θ is obtained by maximizing Equation (9), which is equivalent to minimizing E = − l o g ( l ( θ | Z ) ) . Then, after rearranging terms, the following expression is obtained:
E ( θ | Z ) = − ∑ i = 1 K r i ln y i + ( 1 − r i ) ln ( 1 − y i )
Equations (5) and (10) are very similar. However, Equation (10) is calculated from training data, where r i is either one or zero. This means that for a sample with r = 0 , the first term in the corresponding element of the sum cancels out, whereas when r = 1 , the second term cancels out. In other words, the contribution of every training sample to Equation (10) is either ln y i or ln ( 1 − y i ) . A training sample never contributes in both terms. By contrast, Equation (5) is calculated from a set of samples whose classes are unknown. Equation (5) can be interpreted as follows: the contribution of each sample is a weighted combination of both ln h θ ( x i ) and ln ( 1 − h θ ( x i ) ) terms, where the weights, p i and ( 1 − p i ) , are estimated from the fragility curve.
To convince ourselves that solving Equation (4) will yield a consistent discriminant function, let us perform an exercise. For a feature vector x b whose probabilities of collapse and no collapse are 0.9 and 0.1, respectively, it is expected that a good discriminant function should classify the sample x b as collapsed. Now, let us consider two vector candidates: θ 1 , which classifies x b as belonging to C 1 (collapsed), and θ 2 , which classifies x b as belonging to C 2 (non-collapsed). The contributions of the sample x b to the sum in Equation (5), after the introduction of the factor −1, for both candidates are as follows:
T θ 1 = − 0.9 ln h θ 1 ( x b ) − 0.1 ln ( 1 − h θ 1 ( x b ) )
T θ 2 = − 0.9 ln h θ 2 ( x b ) − 0.1 ln ( 1 − h θ 2 ( x b ) )
Considering the assumption that θ 1 and θ 2 classify x b as collapsed and non-collapsed, respectively, the inequalities below must follow:
h θ 1 ( x b ) ≥ 0.5
h θ 2 ( x b ) < 0.5
From Equations (13) and (14) and Figure 1, it is easy to conclude that T θ 1 < T θ 2 , and therefore, the solution of Equation (4) will produce a discriminant function that will classify the majority of samples with a high probability of collapse as collapsed. Likewise, most samples with a low probability of collapse will be classified as non-collapsed.

3. Application of the IHF Method to the 2011 Great East Japan Earthquake and Tsunami

3.1. A Problem of Bi-Dimensional Dataset Features Using a Linear Threshold

In this section, the proposed classification method is assessed by comparing its results with ground-truth data on damaged buildings provided by the Japanese Ministry of Land, Infrastructure, Transportation and Tourism (MLIT) [44]. The experiments were conducted on the same dataset used in [41]. The dataset was prepared from two TerraSAR-X images of the coastal area of Miyagi Prefecture (Figure 2b) and a geocoded building footprint inventory. The images were recorded on 12 October 2010 and 13 March 2011, and thus, this dataset can be used to extract infrastructural damage caused by the Great East Japan Earthquake and Tsunami, which occurred on 11 March 2011. The dataset contains 31,262 samples (N = 31,262), each representing a building located in the affected area, with two features per sample ( m = 2 ). The first feature is the average difference in backscattering (d) between the two TerraSAR-X images within a rectangular box that contains the building’s footprint. To ensure the inclusion of layover and shadowing effects, a minimum distance of 5 m was established between the edges of the rectangular box and the building footprint. The second feature is the correlation coefficient (r) between the pixels of the two images located inside the same rectangular box. Thus, each vector sample has the form x i = ( 1 , d , r ) . Figure 3a shows a scatter plot of the bi-dimensional dataset, in which the colored marks represent the sample density.
For the estimation of the probability of collapse for each sample ( p i ), the spatial distribution of the hazard, i.e., the inundation depth (Figure 2c), and a fragility curve for building collapse (Figure 3b) were used. Two cases were analyzed in this study. In the first case, the fragility curve proposed by Koshimura et al. [26] was employed. In the second case, the fragility curve proposed by Suppasri et al. [28] was used. Based on the geolocation of each sample x i , its inundation depth was extracted from Figure 2c. Then, p i was calculated using one of the fragility functions (Figure 3b).
For the calibration of the discriminant function using the IHF classification method, it is necessary to use samples that exhibit a uniform distribution with respect to the hazard. Thus, the samples were grouped according to their inundation depth in ranges corresponding to multiples of 50 cm. Then, from each of the first 14 groups, 373 samples were randomly extracted. Thus, a total of 5222 samples was used for discriminant function calibration. The discriminant functions obtained for the first case (i.e., with the fragility function of Koshimura et al. [26]) are displayed in Figure 4. The entire dataset (31,262 samples) was grouped in accordance with the damage state as reported in the survey conducted by MLIT [44]. Seven damage state (DS) levels, from DS0 (undamaged) to DS6 (washed away), were defined in this field survey. Figure 4a–g show the scatter plots of the samples separated by DS together with the resulting discriminant functions. As previously, the colored marks denote the density of the samples. To investigate the effects of the random sampling, the IHF method was performed several times, and the differences among the results are barely noticeable, as seen from Figure 4a–g, where the discriminant functions from four different tests are shown. Thus, it was concluded that the effect of the random sampling was insignificant. The gradient descent algorithm was employed to solve Equation (4). A learning rate ( α ) of 0.05 was used. For every test, the initial vector θ was composed of zeros, that is, θ = ( 0 , 0 , 0 ) . Figure 4h shows the variations of J ( θ ) throughout the successive iterations. Convergence was reached after approximately 3000 iterations. From a qualitative evaluation, it is evident that the discriminant function separates the majority of the DS6 samples from the rest. To provide a quantitative evaluation, Table 1 presents the confusion matrix. To calculate the overall accuracy (OA), all samples with damage states from DS0–DS5 were merged. The results indicate an OA of 82.2% and a Cohen’s kappa coefficient of 0.62. The producer accuracies (PAs) for samples with different DSs were also calculated. A high PA is observed for samples with damage states from DS0–DS3. By contrast, the lowest PAs are observed for samples of DS4–DS5.
Similarly, the results obtained using the fragility function of Suppasri et al. [28] are displayed in Figure 5 and Table 2. The same parameters and original conditions as the previous case were employed here. Again, the effect of the random sampling was irrelevant. However, slightly faster convergence was observed. An OA of 87.5% and a Cohen’s kappa of 0.69 were achieved. In addition, high PA values with a minimum of 78.3% (for DS5) were achieved.

3.2. Generalization of the Discriminant Function

To prove that the IHF method retains the benefits of a traditional machine learning technique in terms of independence of dimensionality and discriminant function shape, two additional cases are presented. A modification of the discriminant function to consider non-linear terms can be easily implemented. For instance, let us add quadratic terms to our discriminant function (Equations (1) and (2)):
f ( x ) = θ 0 + θ 1 x 1 + θ 2 x 2 + θ 3 x 1 2 + θ 4 x 2 2
The resulting non-linear threshold is depicted in Figure 6, in which all samples with damage states from DS0–DS5 are presented as “non-collapsed buildings”. The confusion matrices are shown in Table 3 and Table 4, from which the same level of accuracy is observed as for the previous results.
In a similar way, the number of features considered can be increased. Figure 7 presents the parallel coordinate plot of a dataset composed of seven features: the correlation coefficient (r), the mean difference (d), the median difference, the mode of the differences, the standard deviation of the difference, the maximum value of the differences and the minimum value of the differences. The additional five features were calculated using the same TerraSAR-X images and the same rectangular boxes for the samples. With these features, the discriminant function is modified as follows:
f ( x ) = θ 0 + θ 1 x 1 + θ 2 x 2 + θ 3 x 3 + θ 4 x 4 + θ 5 x 5 + θ 6 x 6 + θ 7 x 7
In Figure 7, the red marks correspond to the samples classified as collapsed buildings. Likewise, the blue marks denote the samples classified as non-collapsed buildings. Note that all features have been normalized with respect to the maximum absolute value. The confusion matrix is shown in Table 5.

3.3. Discussion of the Case Study

In this section, additional comments regarding the results are addressed. In the experimental analysis, the effect of the fragility curve was observed by comparing the results obtained using the fragility curve of Koshimura et al. [26] and that of Suppasri et al. [28]. The reason why the fragility curve of Suppasri et al. [28] resulted in better performance is attributed to the fact that it was created specifically for Japanese buildings. By contrast, Koshimura et al. [26] proposed fragility curves for buildings in Indonesia. Thus, it can be seen that to achieve the best performance, it is necessary to use a fragility function designed for the specific target area. This issue might represent an important pitfall. However, at present, there is a significant number of publications concerning building damage functions for specific regions. For instance, tsunami fragility functions have been proposed for the building stocks in Sri Lanka [45], Chile [27], Thailand [46] and Samoa [47,48]. Likewise, earthquake fragility functions have been reported for Japan [24] and South America [49]. Another issue is the number of fragility functions used in the case study. According to the building inventory data within the affected area, 23,767 (89.7%) wooden buildings, 1753 (6.6%) steel-frame buildings and 975 (3.7%) reinforced concrete buildings were identified. Furthermore, buildings with less than two floors represent 99% of the total. Based on these aggregates, it was assumed that the majority of the affected buildings were wooden buildings. Therefore, only one fragility curve was employed in each test. However, in most possible applications, a large variety of building types might be found, and it is essential to point out that only one fragility function may be insufficient.
With respect to the calculation of the features, it was decided to use a bigger region than the building footprint in order to consider the shadowing and layover effects. As a secondary effect, some pixels of neighboring buildings were included, as well. However, the features r and d are aggregate values of the rectangular box used for their calculations. Thus, it should not affect the results significantly. This assumption is confirmed from Figure 4 and Figure 5, in which the collapsed and non-collapsed buildings are grouped in different regions.
Regarding the low accuracy observed for the samples classified as DS5 (Figure 4 and Figure 5 and Table 1, Table 2, Table 3, Table 4 and Table 5), according to MLIT, a building was classified as DS5 when its main structure had been compromised. Thus, it is very likely that the DS5 samples might include both collapsed and non-collapsed buildings. However, in the accuracy evaluation, all DS5 samples were assumed to be non-collapsed buildings. This assumption might be the reason for the low accuracy on the DS5 samples, which, in turn, was the main reason for the low user accuracy (UA). Another source of misclassification in the results is the presence of speckle noise in the SAR images. The speckle noise is associated with constructive and/or destructive interference of the reflected electromagnetic waves from different objects within a pixel. The TerraSAR-X images have been speckle filtered; however, it cannot be completely removed. Its effect is reflected in the features r and d, as can be shown in Figure 4, where some samples of collapsed and non-collapsed buildings are overlapped.
A comparison with classification methods based on supervised machine learning is desirable here. The results of Wieland et al. [16] and Bai et al. [22] were chosen for this purpose. Both research groups used the same TerraSAR-X images as in the present study. Moreover, the same information from MLIT was used to create the training data. However, the features (i.e., the elements of the vector x i ) were calculated following a different procedure. Thus, the comparison here should be considered only for reference. It cannot be concluded that one method is superior to the other. In addition to the images used in our case study, Wieland et al. [16] used three further SAR images: a TerraSAR-X image acquired on 19 June 2011 and two ALOS PALSAR images acquired on 5 October 2010 and 7 April 2011. Wieland et al. [16] performed several tests to assess the effects of the number of images, the number of features, etc. For the present comparison, the results they obtained using the same images as in our study are used. Table 6 shows the evaluation of the classification performance. Here, the user accuracy (UA), producer accuracy (PA) and a combination of both scores ( F 1 = ( U A − 1 + P A − 1 ) − 1 ) are presented as accuracy measures [50]. The last row (Total) lists the average scores for the collapsed and non-collapsed samples. Our results show the same level of accuracy as the results of Wieland et al. [16]. Bai et al. [22] addressed a more challenging problem, namely the detection of damaged buildings based on post-event imagery only. Moreover, three classes were defined: washed-away buildings, collapsed buildings and slightly damaged buildings. Deep neural networks were used for classification. The training samples were based on tiles, to which class labels were allocated in accordance with the label of the majority of the buildings in each tile. An overall accuracy of 74.8% and a Cohen’s kappa coefficient of 0.60 were achieved. Once again, our results reach a comparable level of accuracy (Table 1, Table 2, Table 3, Table 4 and Table 5). Before concluding these comparisons, it is worthwhile to recall that the training data in both studies, those of Wieland et al. [16] and Bai et al. [22], were prepared using the information from the MLIT field survey [44], of which the first report was released on 23 August 2011, five months after the disaster event.
Regarding run-time performance, the run-time strongly depends on the specific details of the hardware and compiler used for analysis. For this paper, the IHF method was implemented in the Python programming language, Version 2.7 [51], compiled with Microsoft Visual C++ 2008 for 32-bit systems. The Numerical Python (NumPy) library [52] was used to vectorize the calculations. The classifications reported in this case study were performed on an HP Z240 SFF Workstation (3.70 GHz). The calibration of the linear threshold function from the bi-dimensional dataset as presented in Section 3.1 took 8.98 s; the calibration of the non-linear threshold function from the bi-dimensional dataset (Section 3.2) took 9.05 s; and the calibration of the linear threshold function from the seven-dimensional dataset (Section 3.2) took 11.54 s. These observed run-times are negligible compared with the time required for the acquisition of the SAR imagery. However, the run-time can still be significantly reduced, and this will be addressed in a future publication. For instance, Python, being an interpreted language, exhibits lower performance than a compiled language, such as C++ or Fortran. Moreover, note that the run-times shown represent the time required for 10,000 iterations (Figure 4h and Figure 5h), which would not be necessary if the convergence were to be controlled in the algorithm. Furthermore, the gradient descent method for minimizing J ( θ ) could be replaced with a more complex algorithm, in which, for instance, the learning rate ( α in Equation (6)) could be manipulated during run-time to achieve faster convergence.

4. Conclusions and Prospects

A novel unsupervised building damage classification method called IHF has been defined and implemented. The IHF method is able to estimate which buildings have collapsed during a disaster event without any training samples and therefore can be used to achieve a quick response. The method is able to overcome the difficulties encountered when using supervised machine learning techniques while maintaining some of their advantages, such as the capacity to adapt to n-dimensional problems and the possibility of using non-linear discriminant functions. In our approach, training data are replaced with probabilistic information for each sample. The probability of collapse is estimated from the spatial distribution of the hazard under consideration and one or more fragility functions. The approach is based on a modification of the well-known logistic regression technique, a standard machine learning algorithm.
The IHF method was tested using data from the 2011 Great East Japan Earthquake and Tsunami. The features used in the analysis were collected from a pair of TerraSAR-X satellite images recorded before and after the event. The probability of collapse for each sample was calculated from the spatial distribution of the inundation depth and fragility curves for building collapse published in previous studies. The collapsed building estimates obtained using the calibrated discrimination function displayed good agreement with the results of a field survey performed by the Ministry of Land, Infrastructure, Transportation and Tourism. The overall accuracy and the Cohen’s kappa coefficient were greater than 81% and 0.61, respectively.
In principle, the IHF method can be embedded into a quick response framework for large-scale disasters. To this end, the calculation of the spatial distribution of the hazard must be implemented, as well. Currently, such hazard distributions can be calculated via instrumentation, numerical simulations or a combination of the two. Therefore, this requirement does not represent a limitation. Furthermore, a geocoded building database with a corresponding set of fragility functions must be available in advance and ready for use. It is intended that the IHF method must be performed for every new event. That is, the discriminant function calibrated in the case study cannot be used in other tsunami cases.
In future studies, several issues must be addressed. It will be necessary to prove the reliability of the IHF method in other scenarios, such as earthquakes. Moreover, it will be necessary to perform a sensitivity analysis to analyze the effects of uncertainties in the features, the fragility functions and the hazard distribution. Furthermore, the use of a more robust non-linear discriminant function using kernels will be investigated.

Acknowledgments

This research was supported by the Japan Science and Technology Agency (JST) through the SICORPproject “Increasing Urban Resilience to Large Scale Disaster: The Development of a Dynamic Integrated Model for Disaster Management and Socio-Economic Analysis (DIM2SEA)” and the JST (CREST) Project (Grant Number JP-MJCR1411).

Author Contributions

Luis Moya, Erick Mas and Bruno Adriano conceived and designed the methodology; Luis Moya and Luis R. Marval Perez contributed in the development of the equations; Fumio Yamazaki and Shunichi Koshimura supervised and guided this research troughout every stage. Luis Moya wrote the paper.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Maruyama, Y.; Yamazaki, F.; Matsuzaki, S.; Miura, H.; Estrada, M. Evaluation of Building Damage and Tsunami Inundation Based on Satellite Images and GIS Data Following the 2010 Chile Earthquake. Earthq. Spectra 2012, 28, S165–S178. [Google Scholar] [CrossRef]
  2. Liu, W.; Yamazaki, F.; Gokon, H.; Koshimura, S. Extraction of Tsunami-Flooded Areas and Damaged Buildings in the 2011 Tohoku-Oki Earthquake from TerraSAR-X Intensity Images. Earthq. Spectra 2013, 29, S183–S200. [Google Scholar] [CrossRef]
  3. Hashemi-Parast, S.O.; Yamazaki, F.; Liu, W. Monitoring and evaluation of the urban reconstruction process in Bam, Iran, after the 2003 M w 6.6 earthquake. Nat. Hazards 2017, 85, 197–213. [Google Scholar] [CrossRef]
  4. Nakmuenwai, P.; Yamazaki, F.; Liu, W. Automated Extraction of Inundated Areas from Multi-Temporal Dual-Polarization RADARSAT-2 Images of the 2011 Central Thailand Flood. Remote Sens. 2017, 9, 78. [Google Scholar] [CrossRef]
  5. Moya, L.; Yamazaki, F.; Liu, W.; Yamada, M. Detection of collapsed buildings due to the 2016 Kumamoto, Japan, earthquake from Lidar data. Nat. Hazards Earth Syst. Sci. Discuss. 2017, 2017, 1–20. [Google Scholar] [CrossRef]
  6. Dong, L.; Shan, J. A comprehensive review of earthquake-induced building damage detection with remote sensing techniques. ISPRS J. Photogramm. Remote Sens. 2013, 84, 85–99. [Google Scholar] [CrossRef]
  7. Balz, T.; Liao, M. Building-damage detection using post-seismic high-resolution SAR satellite data. Int. J. Remote Sens. 2010, 31, 3369–3391. [Google Scholar] [CrossRef]
  8. Chen, Q.; Li, L.; Jiang, P.; Liu, X. Building collapse extraction using modified freeman decomposition from post-disaster polarimetric SAR image. In Proceedings of the 2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Beijing, China, 10–15 July 2016; pp. 5769–5772. [Google Scholar]
  9. Dong, H.; Xu, X.; Gui, R.; Song, C.; Sui, H. Metric learning based collapsed building extraction from post-earthquake PolSAR imagery. In Proceedings of the 2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Beijing, China, 10–15 July 2016; pp. 4742–4745. [Google Scholar]
  10. Bai, Y.; Adriano, B.; Mas, E.; Gokon, H.; Koshimura, S. Object-based Building Damage Assessment Methodology Using Only Post Event ALOS-2/PALSAR-2 Dual Polarimetric SAR Intensity Images. J. Disaster Res. 2017, 12, 259–271. [Google Scholar] [CrossRef]
  11. Gokon, H.; Koshimura, S.; Matsuoka, M. Object-Based Method for Estimating Tsunami-Induced Damage Using TerraSAR-X Data. J. Disaster Res. 2016, 11, 225–235. [Google Scholar] [CrossRef]
  12. Matsuoka, M.; Nojima, N. Building Damage Estimation by Integration of Seismic Intensity Information and Satellite L-band SAR Imagery. Remote Sens. 2010, 2, 2111–2126. [Google Scholar] [CrossRef]
  13. Karimzadeh, S.; Mastuoka, M. Building Damage Assessment Using Multisensor Dual-Polarized Synthetic Aperture Radar Data for the 2016 M 6.2 Amatrice Earthquake, Italy. Remote Sens. 2017, 9, 330. [Google Scholar] [CrossRef]
  14. Gokon, H.; Koshimura, S.; Meguro, K. Verification of a Method for Estimating Building Damage in Extensive Tsunami Affected Areas Using L-Band SAR Data. J. Disaster Res. 2017, 12, 251–258. [Google Scholar] [CrossRef]
  15. Liu, W.; Yamazaki, F. Extraction of collapsed buildings in the 2016 Kumamoto earthquake using multi-temporal PALSAR-2 data. J. Disaster Res. 2017, 12, 241–250. [Google Scholar] [CrossRef]
  16. Wieland, M.; Liu, W.; Yamazaki, F. Learning Change from Synthetic Aperture Radar Images: Performance Evaluation of a Support Vector Machine to Detect Earthquake and Tsunami-Induced Changes. Remote Sens. 2016, 8, 792. [Google Scholar] [CrossRef]
  17. Gokon, H.; Post, J.; Stein, E.; Martinis, S.; Twele, A.; Mück, M.; Geiß, C.; Koshimura, S.; Matsuoka, M. A Method for Detecting Buildings Destroyed by the 2011 Tohoku Earthquake and Tsunami Using Multitemporal TerraSAR-X Data. IEEE Geosci. Remote Sens. Lett. 2015, 12, 1277–1281. [Google Scholar] [CrossRef]
  18. Jia, L.; Li, M.; Wu, Y.; Zhang, P.; Chen, H.; An, L. Semisupervised SAR Image Change Detection Using a Cluster-Neighborhood Kernel. IEEE Geosci. Remote Sens. Lett. 2014, 11, 1443–1447. [Google Scholar] [CrossRef]
  19. Frank, J.; Rebbapragada, U.; Bialas, J.; Oommen, T.; Havens, T.C. Effect of Label Noise on the Machine-Learned Classification of Earthquake Damage. Remote Sens. 2017, 9, 803. [Google Scholar] [CrossRef]
  20. Janalipour, M.; Mohammadzadeh, A. A Fuzzy-GA Based Decision Making System for Detecting Damaged Buildings from High-Spatial Resolution Optical Images. Remote Sens. 2017, 9, 349. [Google Scholar] [CrossRef]
  21. Bai, Y.; Adriano, B.; Mas, E.; Koshimura, S. Machine Learning Based Building Damage Mapping from ALOS-2/PALSAR-2 SAR Imagery: Case Study of 2016 Kumamoto Earthquake. J. Disaster Res. 2017, 12, 646–655. [Google Scholar] [CrossRef]
  22. Bai, Y.; Gao, C.; Singh, S.; Koch, M.; Adriano, B.; Mas, E.; Koshimura, S. A Framework of Rapid Regional Tsunami Damage Recognition from Post-event TerraSAR-X Imagery Using Deep Neural Networks. IEEE Geosci. Remote Sens. Lett. 2018, 15, 43–47. [Google Scholar] [CrossRef]
  23. Ohta, Y.; Murakami, H.; Watoh, Y.; Koyama, M. A Model for Evaluating Life Span Characteristics of Entrapped Occupants by an Earthquake. In Proceedings of the 13th World Conference on Earthquake Engineering, Vancouver, BC, Canada, 1–6 August 2004; p. 9. [Google Scholar]
  24. Yamazaki, F.; Murao, O. Vulnerability Functions for Japanese Buildings Based on Damage Data from the 1995 Kobe Earthquake. In Implications of Recent Earthquakes on Seismic Risk; Imperial College Press: London, UK, 2000; Volume 2, pp. 91–102. [Google Scholar]
  25. Porter, K.; Kennedy, R.; Bachman, R. Creating Fragility Functions for Performance-Based Earthquake Engineering. Earthq. Spectra 2007, 23, 471–489. [Google Scholar] [CrossRef]
  26. Koshimura, S.; Oie, T.; Yanagisawa, H.; Imamura, F. Developing Fragility Functions For Tsunami Damage Estimation Using Numerical Model And Post-Tsunami Data From Banda Aceh, Indonesia. Coast. Eng. J. 2009, 51, 243–273. [Google Scholar] [CrossRef]
  27. Mas, E.; Koshimura, S.; Suppasri, A.; Matsuoka, M.; Matsuyama, M.; Yoshii, T.; Jimenez, C.; Yamazaki, F.; Imamura, F. Developing Tsunami fragility curves using remote sensing and survey data of the 2010 Chilean Tsunami in Dichato. Nat. Hazards Earth Syst. Sci. 2012, 12, 2689–2697. [Google Scholar] [CrossRef]
  28. Suppasri, A.; Mas, E.; Charvet, I.; Gunasekera, R.; Imai, K.; Fukutani, Y.; Abe, Y.; Imamura, F. Building damage characteristics based on surveyed data and fragility curves of the 2011 Great East Japan tsunami. Nat. Hazards 2013, 66, 319–341. [Google Scholar] [CrossRef]
  29. Cornell, C.A. Engineering seismic risk analysis. Bull. Seismol. Soc. Am. 1968, 58, 1583–1606. [Google Scholar]
  30. McGuire, R. Seismic Hazard and Risk Analysis; Engineering Monographs on Miscellaneous Earthquake Engineering Topics; Earthquake Engineering Research Institute: Oakland, CA, USA, 2004. [Google Scholar]
  31. Martinelli, A.; Cifani, G.; Cialone, G.; Corazza, L.; Petracca, A.; Petrucci, G. Building vulnerability assessment and damage scenarios in Celano (Italy) using a quick survey data-based methodology. Soil Dyn. Earthq. Eng. 2008, 28, 875–889. [Google Scholar] [CrossRef]
  32. Zülfikar, A.C.; Fercan, N.Ö.Z.; Tunç, S.; Erdik, M. Real-time earthquake shake, damage, and loss mapping for Istanbul metropolitan area. Earth Planets Space 2017, 69, 9. [Google Scholar] [CrossRef]
  33. Karimzadeh, S.; Feizizadeh, B.; Matsuoka, M. From a GIS-based hybrid site condition map to an earthquake damage assessment in Iran: Methods and trends. Int. J. Disaster Risk Reduct. 2017, 22, 23–36. [Google Scholar] [CrossRef]
  34. Frolova, N.I.; Larionov, V.I.; Bonnin, J.; Sushchev, S.P.; Ugarov, A.N.; Kozlov, M.A. Seismic risk assessment and mapping at different levels. Nat. Hazards 2017, 88, 43–62. [Google Scholar] [CrossRef]
  35. Moya, L.; Mas, E.; Koshimura, S. Evaluation of tsunami fragility curves for building damage level allocation. Res. Rep. Tsunami Eng. 2017, 34, 33–41. [Google Scholar]
  36. Shabestari, K.T.; Yamazaki, F. Near-fault spatial variation in strong ground motion due to rupture directivity and hanging wall effects from the Chi-Chi, Taiwan earthquake. Earthq. Eng. Struct. Dyn. 2003, 32, 2197–2219. [Google Scholar] [CrossRef]
  37. Kato, T.; Terada, Y.; Nishimura, H.; Nagai, T.; Koshimura, S. Tsunami records due to the 2010 Chile Earthquake observed by GPS buoys established along the Pacific coast of Japan. Earth Planets Space 2011, 63, e5–e8. [Google Scholar] [CrossRef]
  38. Ozawa, S.; Nishimura, T.; Suito, H.; Kobayashi, T.; Tobita, M.; Imakiire, T. Coseismic and postseismic slip of the 2011 magnitude-9 Tohoku-Oki earthquake. Nature 2011, 471, 373–376. [Google Scholar] [CrossRef] [PubMed]
  39. Matsuoka, M.; Yamamoto, N. Web-based Quick Estimation System of Strong Ground Motion Maps Using Engineering Geomorphologic Classification Map and Observed Seismic Records. In Proceedings of the 15th World Conference on Earthquake Engineering, Lisbon, Portugal, 24–28 September 2012; p. 10. [Google Scholar]
  40. Irikura, K.; Miyakoshi, K.; Kamae, K.; Yoshida, K.; Somei, K.; Kurahashi, S.; Miyake, H. Applicability of source scaling relations for crustal earthquakes to estimation of the ground motions of the 2016 Kumamoto earthquake. Earth Planets Space 2017, 69, 10. [Google Scholar] [CrossRef]
  41. Moya, L.; Mas, E.; Adriano, B.; Koshimura, S.; Yamazaki, F.; Liu, W. An Integrated Method to Extract Collapsed Buildings from Satellite Imagery, Hazard Distribution and Fragility Curves. Int. J. Disaster Risk Reduct. 2018. submitted. [Google Scholar]
  42. Alpaydin, E. Introduction to Machine Learning, 3rd ed.; The MIT Press: Cambridge, MA, USA, 2014; Chapter 10. [Google Scholar]
  43. Bishop, C.M. Pattern Recognition and Machine Learning (Information Science and Statistics); Springer-Verlag New York, Inc.: Secaucus, NJ, USA, 2006. [Google Scholar]
  44. Ministry of Land, Infrastructure, Transport and Tourism (MLIT). Results of the Survey on Disaster Caused by the Great East Japan Earthquake (First Report). Available online: http://www.mlit.go.jp/report/press/city07_hh_000053.html (accessed on 22 December 2017).
  45. Murao, O.; Nakazato, H. Vulnerability functions for buildings based on damage survey data in Sri Lanka after the 2004 Indian Ocean Tsunami. In Proceedings of the International Conference on Sustainable Built Environment (ICSBE-2010), Kandy, Sri Lanka, 13–14 December 2010; pp. 371–378. [Google Scholar]
  46. Suppasri, A.; Koshimura, S.; Imamura, F. Developing tsunami fragility curves based on the satellite remote sensing and the numerical modeling of the 2004 Indian Ocean tsunami in Thailand. Nat. Hazards Earth Syst. Sci. 2011, 11, 173–189. [Google Scholar] [CrossRef]
  47. Reese, S.; Bradley, B.A.; Bind, J.; Smart, G.; Power, W.; Sturman, J. Empirical building fragilities from observed damage in the 2009 South Pacific tsunami. Earth-Sci. Rev. 2011, 107, 156–173. [Google Scholar] [CrossRef]
  48. Gokon, H.; Koshimura, S.; Imai, K.; Matsuoka, M.; Namegaya, Y.; Nishimura, Y. Developing fragility functions for the areas affected by the 2009 Samoa earthquake and tsunami. Nat. Hazards Earth Syst. Sci. 2014, 14, 3231–3241. [Google Scholar] [CrossRef]
  49. Villar-Vega, M.; Silva, V.; Crowley, H.; Yepes, C.; Tarque, N.; Acevedo, A.B.; Hube, M.A.; Gustavo, C.D.; María, H.S. Development of a Fragility Model for the Residential Building Stock in South America. Earthq. Spectra 2017, 33, 581–604. [Google Scholar] [CrossRef]
  50. Fawcett, T. An introduction to ROC analysis. Pattern Recognit. Lett. 2006, 27, 861–874. [Google Scholar] [CrossRef]
  51. Python Software Foundation. Python Language Reference, Version 2.7. Available online: https://www.python.org/ (accessed on 22 December 2017).
  52. Oliphant, T.E. Guide to NumPy, 2nd ed.; CreateSpace Independent Publishing Platform: Austin, TX, USA, 2015. [Google Scholar]
Remotesensing 10 00296 g001 550
Figure 1. The functions g ( h θ ( x ) ) = − ln ( h θ ( x ) ) and g ( h θ ( x ) ) = − ln ( 1 − h θ ( x ) ) that contribute to the elements of the summation in Equation (5).
Figure 1. The functions g ( h θ ( x ) ) = − ln ( h θ ( x ) ) and g ( h θ ( x ) ) = − ln ( 1 − h θ ( x ) ) that contribute to the elements of the summation in Equation (5).
Remotesensing 10 00296 g001
Remotesensing 10 00296 g002 550
Figure 2. Study area. (a) Location of the coastal area of Miyagi Prefecture (red rectangle) within the Tohoku region of Japan; (b) RGB color composite of the TerraSAR-X images acquired on 13 March 2011 (red) and 21 October 2010 (green and blue); (c) Inundation depth map of the Great East Japan Earthquake and Tsunami. The inundation values are given in units of meters.
Figure 2. Study area. (a) Location of the coastal area of Miyagi Prefecture (red rectangle) within the Tohoku region of Japan; (b) RGB color composite of the TerraSAR-X images acquired on 13 March 2011 (red) and 21 October 2010 (green and blue); (c) Inundation depth map of the Great East Japan Earthquake and Tsunami. The inundation values are given in units of meters.
Remotesensing 10 00296 g002
Remotesensing 10 00296 g003 550
Figure 3. (a) Scatter plot of the bi-dimensional dataset composed of r and d values. The colored marks denote the densities at the corresponding points; (b) Empirical fragility functions of buildings that collapse due to a tsunami event as proposed by Koshimura et al. [26] (solid line) and by Suppasri et al. [28] (dashed line).
Figure 3. (a) Scatter plot of the bi-dimensional dataset composed of r and d values. The colored marks denote the densities at the corresponding points; (b) Empirical fragility functions of buildings that collapse due to a tsunami event as proposed by Koshimura et al. [26] (solid line) and by Suppasri et al. [28] (dashed line).
Remotesensing 10 00296 g003
Remotesensing 10 00296 g004 550
Figure 4. Discriminant functions obtained using the fragility function of Koshimura et al. [26]. (a–g) Scatter plots of the bi-dimensional dataset separated by damage state (DS), together with the obtained discriminant functions; (h) Variations of the cost function ( J ( θ ) ) throughout the iterative gradient descent algorithm.
Figure 4. Discriminant functions obtained using the fragility function of Koshimura et al. [26]. (a–g) Scatter plots of the bi-dimensional dataset separated by damage state (DS), together with the obtained discriminant functions; (h) Variations of the cost function ( J ( θ ) ) throughout the iterative gradient descent algorithm.
Remotesensing 10 00296 g004
Remotesensing 10 00296 g005 550
Figure 5. Discriminant functions obtained using the fragility function of Suppasri et al. [28]. (a–g) Scatter plots of the bi-dimensional dataset separated by DS, together with the obtained discriminant functions; (h) Variations of the cost function ( J ( θ ) ) throughout the iterative gradient descent algorithm.
Figure 5. Discriminant functions obtained using the fragility function of Suppasri et al. [28]. (a–g) Scatter plots of the bi-dimensional dataset separated by DS, together with the obtained discriminant functions; (h) Variations of the cost function ( J ( θ ) ) throughout the iterative gradient descent algorithm.
Remotesensing 10 00296 g005
Remotesensing 10 00296 g006 550
Figure 6. Non-linear discriminant functions calculated from the fragility functions of (a) Koshimura et al. [26] and (b) Suppasri et al. [28].
Figure 6. Non-linear discriminant functions calculated from the fragility functions of (a) Koshimura et al. [26] and (b) Suppasri et al. [28].
Remotesensing 10 00296 g006
Remotesensing 10 00296 g007 550
Figure 7. Parallel coordinate plot of seven normalized features. Red and blue marks denote samples classified as collapsed and non-collapsed buildings, respectively. The black ticks delimit the range of [ a v e r a g e − s t d , a v e r a g e + s t d ] for each normalized feature.
Figure 7. Parallel coordinate plot of seven normalized features. Red and blue marks denote samples classified as collapsed and non-collapsed buildings, respectively. The black ticks delimit the range of [ a v e r a g e − s t d , a v e r a g e + s t d ] for each normalized feature.
Remotesensing 10 00296 g007
Table
Table 1. Accuracy evaluation of the classification results for collapsed and non-collapsed buildings obtained using a linear discriminant function based on the fragility function of Koshimura et al. [26]. NC: non-collapsed; C: collapsed; DS: damage state; PA: producer accuracy; UA: user accuracy.
Table 1. Accuracy evaluation of the classification results for collapsed and non-collapsed buildings obtained using a linear discriminant function based on the fragility function of Koshimura et al. [26]. NC: non-collapsed; C: collapsed; DS: damage state; PA: producer accuracy; UA: user accuracy.
DS0DS1DS2DS3DS4DS5DS0–DS5DS6TotalUA
NC16682190565548161179201817,52666418,19096.3
C189388994120451616144905814013,04562.4
Total18572578664960201695363222,431880431,235
PA89.884.985.180.069.655.678.192.582.2
Table
Table 2. Accuracy evaluation of the classification results for collapsed and non-collapsed buildings obtained using a linear discriminant function based on the fragility function of Suppasri et al. [28].
Table 2. Accuracy evaluation of the classification results for collapsed and non-collapsed buildings obtained using a linear discriminant function based on the fragility function of Suppasri et al. [28].
DS0DS1DS2DS3DS4DS5DS0–DS5DS6TotalUA
NC17932448626455511475284420,375186122,23691.6
C6413038546922078820566943899977.2
Total18572578664960201695363222,431880431,235
PA96.695.094.292.287.078.390.878.987.5
Table
Table 3. Accuracy evaluation of the classification results for collapsed and non-collapsed buildings obtained using a non-linear discriminant function based on the fragility function of Koshimura et al. [26].
Table 3. Accuracy evaluation of the classification results for collapsed and non-collapsed buildings obtained using a non-linear discriminant function based on the fragility function of Koshimura et al. [26].
DS0DS1DS2DS3DS4DS5DS0–DS5DS6TotalUA
NC16492158553447191153192617,13961517,75496.5
C2084201115130154217065292818913,48160.7
Total18572578664960201695363222,431880431,235
PA88.883.783.278.468.053.076.493.081.1
Table
Table 4. Accuracy evaluation of the classification results for collapsed and non-collapsed buildings obtained using a non-linear discriminant function based on the fragility function of Suppasri et al. [28].
Table 4. Accuracy evaluation of the classification results for collapsed and non-collapsed buildings obtained using a non-linear discriminant function based on the fragility function of Suppasri et al. [28].
DS0DS1DS2DS3DS4DS5DS0–DS5DS6TotalUA
NC17882419620554921444273820,086161221,69892.6
C6915944452825189423457192953775.4
Total18572578664960201695363222,431880431,235
PA96.393.893.391.285.275.489.581.787.3
Table
Table 5. Accuracy evaluation of the classification results for collapsed and non-collapsed buildings obtained from a 7-dimensional dataset using a non-linear discriminant function based on the fragility function of Suppasri et al. [28].
Table 5. Accuracy evaluation of the classification results for collapsed and non-collapsed buildings obtained from a 7-dimensional dataset using a non-linear discriminant function based on the fragility function of Suppasri et al. [28].
DS0DS1DS2DS3DS4DS5DS0–DS5DS6TotalUA
NC18112484636056341484288220,655220722,86290.3
C469428938621175017766597837378.8
Total18572578664960201695363222,431880431,235
PA97.596.495.793.687.679.492.174.987.2
Table
Table 6. Performance comparison between the IHF method and that of Wieland et al. [16]. IHF, imagery, hazard and fragility.
Table 6. Performance comparison between the IHF method and that of Wieland et al. [16]. IHF, imagery, hazard and fragility.
Class  Wieland et al.IHF Method
Koshimura et al.Suppasri et al.
UAPAF1UAPAF1UAPAF1
C0.750.770.760.620.920.750.770.790.78
NC0.760.740.750.960.780.860.920.910.91
Total0.760.760.760.790.850.800.840.850.85

Share and Cite

MDPI and ACS Style

Moya, L.; Marval Perez, L.R.; Mas, E.; Adriano, B.; Koshimura, S.; Yamazaki, F. Novel Unsupervised Classification of Collapsed Buildings Using Satellite Imagery, Hazard Scenarios and Fragility Functions. Remote Sens. 2018, 10, 296. https://doi.org/10.3390/rs10020296

AMA Style

Moya L, Marval Perez LR, Mas E, Adriano B, Koshimura S, Yamazaki F. Novel Unsupervised Classification of Collapsed Buildings Using Satellite Imagery, Hazard Scenarios and Fragility Functions. Remote Sensing. 2018; 10(2):296. https://doi.org/10.3390/rs10020296

Chicago/Turabian Style

Moya, Luis, Luis R. Marval Perez, Erick Mas, Bruno Adriano, Shunichi Koshimura, and Fumio Yamazaki. 2018. "Novel Unsupervised Classification of Collapsed Buildings Using Satellite Imagery, Hazard Scenarios and Fragility Functions" Remote Sensing 10, no. 2: 296. https://doi.org/10.3390/rs10020296

APA Style

Moya, L., Marval Perez, L. R., Mas, E., Adriano, B., Koshimura, S., & Yamazaki, F. (2018). Novel Unsupervised Classification of Collapsed Buildings Using Satellite Imagery, Hazard Scenarios and Fragility Functions. Remote Sensing, 10(2), 296. https://doi.org/10.3390/rs10020296

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop