Open access peer-reviewed chapter

Monitoring of Tsunami/Earthquake Damages by Polarimetric Microwave Remote Sensing Technique

Written By

Dong Li, Yunhua Zhang, Liting Liang, Jiefang Yang and Xun Wang

Submitted: 25 October 2019 Reviewed: 17 January 2020 Published: 20 February 2020

DOI: 10.5772/intechopen.91242

From the Edited Volume

Tsunami - Damage Assessment and Medical Triage

Edited by Mohammad Mokhtari

Chapter metrics overview

898 Chapter Downloads

View Full Metrics

Abstract

Polarization characterizes the vector state of EM wave. When interacting with polarized wave, rough natural surface often induces dominant surface scattering; building also presents dominant double-bounce scattering. Tsunami/earthquake causes serious destruction just by inundating the land surface and destroying the building. By analyzing the change of surface and double-bounce scattering before and after disaster, we can achieve a monitoring of damages. This constitutes one basic principle of polarimetric microwave remote sensing of tsunami/earthquake. The extraction of surface and double-bounce scattering from coherency matrix is achieved by model-based decomposition. The general four-component scattering power decomposition with unitary transformation (G4U) has been widely used in the remote sensing of tsunami/earthquake to identify surface and double-bounce scattering because it can adaptively enhance surface or double-bounce scattering. Nonetheless, the strict derivation in this chapter conveys that G4U cannot always strengthen the double-bounce scattering in urban area nor strengthen the surface scattering in water or land area unless we adaptively combine G4U and its duality for an extended G4U (EG4U). Experiment on the ALOS-PALSAR datasets of 2011 great Tohoku tsunami/earthquake demonstrates not only the outperformance of EG4U but also the effectiveness of polarimetric remote sensing in the qualitative monitoring and quantitative evaluation of tsunami/earthquake damages.

Keywords

  • disaster monitoring
  • damage evaluation
  • tsunami
  • earthquake
  • microwave remote sensing
  • synthetic aperture radar (SAR)
  • polarimetric SAR (PolSAR)
  • polarimetric decomposition
  • scattering model
  • unitary transformation

1. Introduction

Tsunami and earthquake seriously endanger people’s lives and properties. Efficient and accurate monitoring and assessment are of crucial importance for the fast response, management, and mitigation of the disasters [1, 2, 3]. Compared with the optical remote sensing, microwave remote sensing technology such as synthetic aperture radar (SAR) has been widely applied to monitoring natural and human-induced disasters for its all-day and all-weather working capacity [4].

Polarization is an essential property of the electromagnetic wave [5, 6, 7, 8]. The polarization state of wave will change when interacting with ground object. For example, rough natural surface such as land and water often induces the strong Bragg surface scattering, while building often presents the dominant double-bounce scattering because of the dihedral corner reflectors formed by ground and the vertical wall of building. Therefore, by analyzing the polarization of the scattering wave, we can acquire the physical and geometrical information regarding the object. This is the main task of SAR polarimetry (PolSAR) [9, 10, 11].

Tsunami is often accompanied by earthquake and flooding [1, 2, 3]. It damages and inundates the buildings and causes the collapse of the ground-wall dihedral structures as well as the enhancement of the direct surface scatterers. Therefore, by analyzing the power of double-bounce scattering and surface scattering before and after the event, we can achieve an efficient monitoring of the disasters. This simple strategy has been successfully adopted in the polarimetric microwave remote sensing of tsunami/earthquake [12, 13, 14, 15, 16, 17, 18, 19, 20, 21].

Nonetheless, the extraction of double-bounce scattering and surface scattering from PolSAR image is not so straightforward because each pixel in PolSAR is a 3×3 complex coherency matrix T with nine degrees of freedom (DoF). A widely used approach to achieve this is to decompose T on the canonical scattering models [22]. The first such decomposition was devised by Freeman and Durden [23] which expands T on the surface scattering, double-bounce scattering, and volume scattering (describes the complex scattering in vegetation area). This three-component decomposition, however, is responsible for only five DoF of T because of the symmetric reflection assumption. This assumption was tackled by Yamaguchi et al. [24] by introducing a fourth helix component and two additional models of volume scattering. The resulted four-component decomposition (Y4O) then only leaves three DoF unaccounted: the 13 element of T, i.e., T13, and the real part of the 23 element of T, i.e., ReT23. A same target will present differently by a simple rotation about the line of sight of radar. Deorientation should be first conducted on T to eliminate the influence [25]. As a result, ReT23 changes to zero and Y4O with rotation (Y4R) accounts for seven DoF [26]. Based on Y4R, Sato et al. [27] further proposed to add a new model to characterize volume scattering generated by even-bounce structure. However, Sato’s extended Y4R (S4R) still leaves T13 unaccounted. To solve this, Singh et al. [28] in 2013 proposed a general four-component decomposition (G4U) based on a special unitary matrix. G4U enables T13 included in the accounted models by conducting unitary transformation to the rotated version of T . Singh et al. [28] claimed that G4U could make full use of polarimetric parameters. As a result, in comparison with the four-component decompositions such as S4R and Y4R, G4U could enhance double-bounce scattering power over urban area while enhancing surface scattering contribution over an area where surface scattering is preferable [28]. This makes G4U very suitable to the remote sensing of tsunami/earthquake [16, 20] and establishes G4U the state-of-the-art four-component scattering power decomposition [29, 30].

This chapter is dedicated to enable an extension to G4U for better monitoring of tsunami/earthquake disaster. It is indicated that the unitary transformation in G4U adds a T13-related but redundant balance equation to the original self-contained equation system in Y4R and S4R. Then T13 is accounted for by G4U, but we obtain no exact solution to the system but the approximate ones. The general expression of the approximate solutions enables a generalized G4U (GG4U), while G4U and S4R represent two special forms. A dual G4U (DG4U) is also obtained. The general solution indicates that G4U cannot always enhance the double-bounce scattering power over urban area nor strengthen the surface scattering power over the area where surface scattering is dominant unless we adaptively integrate G4U and DG4U for an extended G4U (EG4U). Experiments on the PolSAR images of Miyagi Prefecture, Japan, acquired by the L-band spaceborne ALOS-PALSAR system before and after the March 11, 2011, Off-Tohoku 9.0 tsunami/earthquake demonstrate not only the outperformance of EG4U but also the effectiveness of polarimetric remote sensing in the monitoring of tsunami/earthquake disaster.

The remainder of this chapter is arranged as follows. Section 2 presents the basic principle of PolSAR and the polarization descriptors first. The advanced four-component scattering power decompositions are then described in Section 3 to develop the EG4U. By decomposing the ALOS-PALSAR datasets of the 2011 great Tohoku tsunami/earthquake using EG4U, Section 4 evaluates and analyzes the polarimetric monitoring of disaster damages further. The chapter is eventually concluded in Section 5.

Advertisement

2. SAR polarimetry and polarization descriptors

SAR is an active microwave remote sensing technique dedicated to acquire the large-scaled 2D coherent image of the earth’s surface reflectivity [9]. It transmits microwave pulses and receives the backscattering from the illuminated terrain to synthesize a high spatial resolution image. Such an active operation enables SAR an all-day working capacity independent of solar illumination. In addition, operating in the microwave region of electromagnetic spectrum avoids the effects of rain and clouds, which allows SAR an almost all-weather continuous monitoring of the earth surface [9].

Polarization characterizes the vector state of the electromagnetic wave. The polarization state of wave will change when interacting with a ground object. By processing and analyzing such change of polarization, we can obtain the material, roughness, shape, and orientation information regarding the object. The core of this change is the (Sinclair) scattering matrix S of the object, which transforms the incident electric filed EI into the scattered electric filed ES [31]:

ES=ejkrrSEIEHSEVS=ejkrrSHHSHVSVHSVVEHIEVIE1

where r denotes the distance from radar to ground object, k is the wave number, and subscript H or V represents the horizontal or vertical polarization. Matrix S is obtained by first transmitting H-polarized wave (EHI) and receiving scatterings in H - and V-polarization simultaneously to measure the first column SHH and SVH and then transmitting V-polarized wave (EVI) and also receiving in H- and V-polarization simultaneously for the second column SHV and SVV. In reciprocal backscattering, we have SHV=SVH and matrix S covers five DoF then.

Generally, almost all the ground scatterers are situated in the dynamically changing environment and subjected to spatial and/or temporal variations [32]. Such scatterer is called the distributed target, and we can no longer model its scattering with a determined scattering matrix S . The 3×3 coherency matrix T is then constructed as the statistical average of the acquired scatterings to describe the second-order moment of the fluctuations [9]:

T=kk=T11T12T13T21T22T23T31T32T33,k=12SHH+SVVSHHSVV2SHVE2

where · and superscript represent the operations of ensemble average and conjugate transpose and k denotes the Pauli vector. The spatial/temporal depolarization pushes the DoF of T to nine. Therefore, different from the conventional SAR image, each pixel in PolSAR image is not a complex number but a 3×3 coherency matrix T.

The coherency matrix T in Eq. (2) is expressed in the H-V polarization basis; we can also formulate it in some other orthonormal basis by simply taking the unitary transformation of T:

UnitaryT=defU3TU3E3

where U3 is the special unitary matrix that describes the relationship between H-V and the new orthonormal basis. Target deorientation is just based on the real rotation matrix [25]:

U3θ=1000cos2θsin2θ0sin2θcos2θ,2θ=12tan12ReT23T22T33E4

Combining Eq. (4) into Eq. (3), the deoriented coherency matrix T is

T=U3θTU3θ=T11T12T13T21T22jImT23T31jImT32T33.E5

Deorientation makes T23 become purely imaginary and reduces DoF from nine to eight. In order to eliminate the imaginary part further, Singh et al. developed an imaginary rotation matrix [28]:

U3φ=1000cos2φjsin2φ0jsin2φcos2φ,2φ=12tan12ImT23T22T33.E6

A coherency matrix T with zero T23 is then achieved:

T=U3φTU3φ=T11T12T13T21T220T310T33.E7
Advertisement

3. Advanced four-component scattering power decompositions

Polarimetric incoherent decomposition plays an important role in the discrimination and recognition of the distributed target [22]. It pursues the scattering mechanism of the unknown target by extracting the dominant or average target (such as the Huynen-type phenomenological dichotomies [7, 32] and the eigenvalue/eigenvector-based target decompositions [9, 33]) from T or expanding T on the canonical models (such as the model-based target decompositions [23, 24, 25, 26, 27, 28]). Among these decompositions, the four-component scattering power decompositions such as Y4R, S4R, and G4U have been a hot topic recently [29].

3.1 Y4R and S4R

Y4R and S4R decompose the target by linearly expanding matrix T on the four canonical scattering models, as illustrated in Figure 1:

Figure 1.

The canonical models involved in the four-component model-based scattering power decompositions.

T=fSTS+fDTD+fVTV+fCTCE8

where TS, TD, TV, and TC denote the surface scattering model, the double-bounce scattering model, the volume scattering model, and the helix scattering model, respectively:

TS=1β0ββ20000,TD=α2α0α10000,TV=ad0db000c,TC=1200001±0j1j.E9

Parameters fS, fD, fV, and fC in Eq. (8) represent the contributions of the four components; β and α in TS and TD are complex parameters; a, b, c, and d in TV are real constants satisfying a+b+c=1, which involve in four volume scattering models and are adaptively selected according to the branch conditions [27, 28]. Combining Eqs. (8) and (9), the S4R/Y4R scattering balance equation system on unknowns fS, fD, fV, fC, α, and β is formulated [26, 27]:

fS+fDα2+fVa=T111)fSβ+fDα+fVd=T122)fSβ2+fD+fVb+fC2=T223)±jfC2=jImT234)fVc+fC2=T335).E10

Nevertheless, we obtain no scattering balance equation on T13 in Eq. (10). Hence, there always exists a T13 -related unaccounted residue in Y4R and S4R.

3.2 G4U

To model T13, G4U uses U3φ to conduct unitary transformation to both sides of Eq. (10) first and then eliminates the influence of φ [28]. As a result, an additional balance equation is brought into G4U, and we obtain the following scattering balance equation system [30]:

fS+fDα2+fVa=T111)fSβ+fDα+fVd=T12+T13fSβ+fDα+fVd=T12T132)fSβ2+fD+fVb+fC2=T223)±jfC2=jImT234)fVc+fC2=T335).E11

Comparing Eq. (11) with Eq. (10), we can find that Eq. (11–2) gives a dichotomy to Eq. (10–2). The redundancy makes Eq. (11) have no such exact solution like Eq. (10) but some approximate ones. In G4U, Singh et al. preferred the first equation of (11–2) only.

3.3 GG4U: generalization of G4U

Obviously, Eq. (11) provides us a generalized G4U (GG4U). Here we focus on the general solution to (11) for the unknowns fS, fD, fV, fC, α, and β. Let

S=T11fVaC1=T12+T13fVdC2=T12T13fVdD=T22fVbfC2C=1+μ2C1+1μ2C2E12

where μ is a real constant. Then Eq. (11) can be rearranged as

fS+fDα2=S1)fSβ+fDα=C2)fSβ2+fD=D3)fC=2ImT234)fV=12c2T33fC5).E13

Eq. (13) comprises of five equations and six unknowns. Following Freeman-Durden [23] and Yamaguchi et al. [24], we can fix α or β in terms of the sign of SD for the superior between surface scattering and double-bounce scattering:

BC>0dominant surface scatteringα=0BC0dominant doublebounce scatteringβ=0E14

where BC=SD. Combining Eqs. (13) and (14), we can then simply obtain the scattering power of each of the four components, i.e., the surface scattering power PS, the double-bounce scattering power PD, the volume scattering power PV, and the helix scattering power PC:

PS=fS1+β2=S+C2S,BC>0SC2D,BC0PD=fD1+α2=DC2S,BC>0D+C2D,BC0PC=fCHT33ImT23PV=12c2T33PCE15

where H· denotes the Heaviside step function, which is used here to adjust the value of PC for nonnegative PV ruling [27]. It can be easily validated that PS+PD+PV+PC=T11+T22+T33. Thus GG4U gives a decomposition of scattering power.

3.4 Special decompositions

By taking appropriate value to μ, we can have some different decompositions, which are denoted as Gμ. Here we are particularly interested to the following special cases of Gμ.

Case (1): G+1G4U

C=C1=T12+T13fVd=CG4U.E16

This is just the parameter C used in G4U. GG4U changes to G4U in this case.

Case (2): G1DG4U

C=C2=T12T13fVd.E17

This acts as the complement of case (1); thus we name it the dual G4U (DG4U).

Case (3): G0S4R

C=C1+C22=T12fVd=CS4R.E18

This is the parameter C used in S4R, i.e., S4R also shows a special form of GG4U. Hence, the essential difference between S4R and G4U just lies in the different definition of parameter C in Eqs. (16) and (18). The unitary transformation is just to enable the T13 entry contained in CG4U and finally in PS and PD. Parameter C defined in Eq. (12) is a generalization of CG4U and CS4R.

3.5 Theoretical evaluation of S4R and G4U

S4R can improve Y4R by strengthening the double-bounce scattering in urban area [27]. Singh et al. [28] indicated that G4U could further improve S4R in this aspect by strengthening surface scattering in the area where surface scattering is preferable to double-bounce scattering, while increasing the double-bounce scattering in the urban area where the double-bounce scattering is preferable to surface scattering. By combining the ruling in Eq. (14), we can formulate these observations as

PSG4UPSS4R,BC>0PDG4UPDS4R,BC0.E19

In terms of the general expression of PS and PD in (15), here we give a simple validation to Eq. (19) by combining μ=0 and μ=1 into Eqs. (12) and (15):

PSG4U=S+C12SPSS4R=S+C1+C224S,BC>0;PDG4U=D+C12DPDS4R=D+C1+C224D,BC0.E20

From Eq. (20) we have

PSG4UPSS4R=2C12C1+C224S,BC>0PDG4UPDS4R=2C12C1+C224D,BC0.E21

Then Eq. (19) will hold if 2C12C1+C220. Obviously, this condition is not always tenable. Hence, despite better performance in some areas, G4U cannot improve S4R for every target area. To tackle this, the extended G4U (EG4U) will be developed in the following as an adaptive combination of G4U and DG4U.

3.6 EG4U: adaptive combination of G4U and DG4U

Combining μ=1 into Eqs. (12) and (15), DG4U surface and double-bounce scattering powers can be formulated as

PSDG4U=S+C22S,BC>0PDDG4U=D+C22D,BC0.E22

Combining Eqs. (20) and (22), after some simple deduction, we obtain

PSG4U+PSDG4U2PSS4R=C1C224S0,BC>0PDG4U+PDDG4U2PDS4R=C1C224D0,BC0E23
PSG4UPSDG4U=C12C22S,BC>0PDG4UPDDG4U=C12C22D,BC0.E24

We can immediately obtain from Eq. (23) that

maxPSG4UPSDG4UPSS4R,BC>0maxPDG4UPDDG4UPDS4R,BC0.E25

From Eq. (24) we obtain

maxPSG4UPSDG4U=PSG4U,BC>0maxPDG4UPDDG4U=PDG4U,BC0BC1>0maxPSG4UPSDG4U=PSDG4U,BC>0maxPDG4UPDDG4U=PDDG4U,BC0,BC10E26

where BC1=C1C2. Eq. (26) just lays the foundation for EG4U:

EG4UG±1=G+1=G4U,BC1>0G1=DG4U,BC10.E27

As the adaptive combination of G4U and DG4U, EG4U is also a special case of GG4U. So we denote it as G±1. By bringing μ=+1 or μ=1 into Eqs. (12) and (15) based on the branch condition BC1, we can achieve the scattering powers of four components in EG4U. Furthermore, from Eqs. (25) to (27), we have

PSEG4U=maxPSG4UPSDG4UPSS4RPSG4UPSDG4U,BC>0PDEG4U=maxPDG4UPDDG4UPDS4RPDG4UPDDG4U,BC0.E28

Compared with S4R and G4U, EG4U increases surface scattering in area where surface scattering is superior to double-bounce scattering and strengthens double-bounce scattering in area where double-bounce scattering is preferable to surface scattering. Therefore, EG4U achieves not only a nice improvement to S4R, but also an effective extension to G4U. This may make EG4U more suitable to the remote sensing of tsunami/earthquake. We will investigate this in Section 4. The procedure of EG4U is outlined in Algorithm 1.

Algorithm 1: EG4U

01: Input: T

02: Conduct deorientation to T for T

03: Compute helix power PC=2ImT23HT33ImT23

04: Calculate branch condition BC

05: Determine volume scattering model based on branch condition

06: Obtain volume scattering power PV=2T33PC/2c

07: Compute parameters S, D, C1, and C2, as well as branch condition BC1

08: Implement SPAN reservation ruling based on S+D

09: if S+D>0

10: Adaptively select between G4U and DG4U based on BC1

11: if BC1>0

12: C=C1

13: else

14: C=C2

15: end if

16: Calculate surface scattering power PS and double-bounce scattering power

PD according to BC

17: if BC>0

18: PS=S+C2/S,PD=DC2/S

19: else

20: PS=SC2/D,PD=D+C2/D

21: end if

22: Implement nonnegative PS and PD ruling

23: else

24: PS=PD=0,PV=T11+T22+T33PC

25: end if

26: Output: PS,PD,PV,PC

Advertisement

4. Monitoring of disaster by EG4U decomposition of ALOS-PALSAR images of 2011 Tohoku tsunami/earthquake

As indicated in Subsection 3.4, G4U and S4R represent two special forms of GG4U of equal status. Hence, G4U cannot fully improve S4R only if we ascend the status of G4U by combining the duality of G4U, i.e., DG4U and G4U together for EG4U. EG4U can adaptively strengthen the surface scattering and double-bounce scattering. Therefore, it may improve the competence and performance of G4U in the remote sensing of damages caused by earthquake/tsunami disaster. We demonstrate these in the following by decomposing the ALOS-PALSAR images of the 2011 great Tohoku tsunami/earthquake using EG4U.

4.1 Great Tohoku earthquake and tsunami

The great Tohoku earthquake is also known as the great Sendai earthquake or the great East Japan earthquake, which was a magnitude 9.0–9.1 (Mw) undersea megathrust earthquake off the coast of northeast Japan (the epicenter is shown in Figure 2 as “

”) that occurred on March 11, 2011, the most powerful earthquake ever recorded in Japan [34]. The earthquake triggered powerful tsunami, which swept the mainland of Japan, killed over 10,000 people (mainly through drowning), and damaged over 1,000,000 buildings (half of them are collapsed and even totally collapsed) [35].

Figure 2.

Location of the great Tohoku tsunami/earthquake epicenter () and the ALOS-PALSAR footprint of the two selected fully polarimetric datasets (red rectangle, pre-event; blue rectangle, post-event).

4.2 Datasets

The Advanced Land Observing Satellite (ALOS) was launched in 2006 by the Japanese Space Agency (JAXA). It has three remote sensing payloads, i.e., the Panchromatic Remote-sensing Instrument for Stereo Mapping (PRISM) for digital elevation mapping, the Advanced Visible and Near Infrared Radiometer type 2 (AVNIR-2) for precise land coverage observation, and the Phased Array type L-band SAR (PALSAR) for all-day/all-weather land observation [36].

To demonstrate the capability of polarimetric remote sensing for damage monitoring, we choose two quad-polarization single-look complex-level 1.1 (ascending orbit) datasets acquired around Miyagi Prefecture, Japan, before and after the earthquake/tsunami with 138 days’ temporal baseline, as summarized in Table 1. The ALOS-PALSAR footprint of the two datasets is shown in Figure 2.

Scene IDAcquire dataIncidence angle1Azimuth resolutionGround-range resolution2
ALPSRP2570907602010-11-2123.802°4.5 m23.5 m
ALPSRP2772207602011-04-0823.836°4.5 m23.5 m

Table 1.

ALOS-PALSAR datasets used in the experiment and their characteristics.

The incidence angle here indicates the incidence angle at the scene center.


The ground-range resolution is defined as the slant-range resolution/sin(incidence angle) [9], while the slant-range resolution of the two datasets is both 9.5 m.


4.3 Method

The flowchart of EG4U-based monitoring and evaluation of damages caused by tsunami/earthquake disaster is illustrated in Figure 3. We first co-register the two datasets based on the image features [37, 38, 39, 40]. The boxcar filtering [9] is then carried out to both datasets to suppress the speckles. To ensure the pixel size in both image directions comparable, the window size for ensemble average is chosen as 2 pixels in ground-range direction and 12 pixels in azimuth direction, i.e., we integrate the scattering matrix S of a total of 24 pixels for the estimation of a coherency matrix T in Eq. (2). From T we calculate the orientation angle θ according to Eq. (4) and implement the deorientation operation for the deoriented coherency matrix ⟨[T′]⟩ according to Eq. (5). Finally, EG4U is used to decompose ⟨[T′]⟩ to extract scattering powers PS, PD, PV, and PC and construct the RGB pseudo-color scattering power visualization result by encoding RGB with PDPVPS. This process is executed on each cell of the two datasets until we obtain the complete pre- and post-event scattering power images shown in Figure 4, based on which we evaluate EG4U on monitoring of the tsunami/earthquake disaster in the following.

Figure 3.

Flowchart of EG4U-based monitoring of tsunami/earthquake disaster.

Figure 4.

Color-coded scattering power image of the study area (a) before and (b) after the great Tohoku tsunami/earthquake disaster. The framed patch regions A, B, and C are extracted for particular analysis.

4.4 Evaluation and analysis

For better comparison and analysis, we also display the optical image of the study area obtained from ©Google Earth in Figure 5. Our intuitive impression of Figure 4(a) and (b) is their consistency and nice correspondence to the optical image. The blue color mainly appears in the water and land areas because of the dominant surface scattering there. The red color mainly arises in the urban area, such as the Ishinomaki City and Higashi-Matsushima City, with a large number of buildings. The ground and the vertical walls of buildings constitute the dihedral corner structures, which generally reflect the dominant double-bounce scattering. Mountain presents the green color, i.e., the dominant volume scattering. The well-developed branch and crown structures of trees on the mountain complicate the scattering process, depolarize the scattering wave, and show themselves as the complex mixed volume scattering in PolSAR image. Therefore, by color-coding the scattering powers obtained by EG4U, we can achieve a nice discrimination of the ground objects.

Figure 5.

Optical image of the study area obtained from ©Google earth. Particular attention is paid to the framed patch regions A, B, and C.

Despite the consistency, we can also observe the obvious difference between the pre- and post-event scattering power images. A lot of red pixels in Figure 4(a) change to blue pixels in Figure 4(b), particularly in the urban areas of Ishinomaki and Higashi-Matsushima, which illustrate the change from the dominant double-bounce scattering to the dominant surface scattering, denote the decrease of the dihedral structures, and indicate the collapse of buildings. Take Ishinomaki City framed in Patch A for instance; it is interesting to observe that the strong change mainly arises in the area by the seaside, while tiny change occurs in the area away from the coast. This finding is also validated by the corresponding optical images acquired before and after the event shown in Figure 6(a) and (b). Therefore, the severe damages brought by the Tohoku tsunami/earthquake are probably mainly due to the flooding rather than the earthquake. Flooding from the Onagawa Bay and the Mangokuura Sea also swept the town of Onagawa framed in Patch B, as shown in Figure 6(c) and (d) in terms of the pre- and post-event optical images. A large majority of red pixels of Patch B in Figure 4(a) change to blue pixels or even green pixels in Figure 4(b), which indicates that nearly all the buildings in Onagawa were badly damaged by the flooding except for a few buildings constructed in high elevation. The collapsed buildings not only present the dominant surface scattering here, but also the dominant volume scattering because of the complex scattering in such mountain area. The biggest change caused by flooding appears in the area along the Kitakami River. Take the town of Kamaya framed in Patch C, for example, as shown in Figure 6(e), besides several buildings, the most part of Kamaya is farmland. This area can be clearly distinguished from the Kitakami River in Figure 4(a) before the disaster. However, after the disaster, nearly all the land and buildings in Kamaya are flooded by the water from Kitakami River as shown in Figure 6(f), which present in Figure 4(b) as the wide distribution of blue pixels and show the dominant surface scattering here. Therefore, by decomposing the pre- and post-event PolSAR datasets with EG4U to construct the color-coded scattering power images, we can achieve a simple but accurate monitoring of the damages caused by tsunami/earthquake disaster.

Figure 6.

Optical images of (first row, i.e. (a) and (b)) patch A, (second row, i.e. (c) and (d)) patch B, and (third row, i.e. (e) and (f)) patch C obtained from ©Google earth (first column, i.e. (a), (c), and (e)) before and (second column, i.e. (b), (d), and (f)) after the 3.11 great Tohoku tsunami/earthquake.

From the above analysis, we can obtain that flooding which resulted from tsunami is the main contributor to the severe damages in the 3.11 great Tohoku earthquake. The flooding destroyed the buildings and inundated the lands. All these damages present themselves in the polarization domain as the change of the dominant scattering mechanism from double-bounce scattering to surface scattering and in the image domain as the change of pixel color from red to blue. The boundary condition BC has been widely used in model-based decomposition as a crucial feature to discriminate surface scattering and double-bounce scattering [23, 24, 26, 27, 28]. As expressed in Eq. (14), BC>0 indicates stronger surface scattering than double-bounce scattering, while BC0 denotes stronger double-bounce scattering than surface scattering. Therefore, besides the qualitative evaluation in terms of color, we can further achieve an quantitative evaluation of the damages by analyzing the dominant scattering according to BC. Figure 7(a) and (b) show the binary images of BC before and after the disaster, respectively. The white pixel denotes BC>0, i.e., the dominant surface scattering, which mainly occupies the water and land areas, while the black one denotes BC0, i.e., the dominant double-bounce scattering, which mainly occupies the urban and mountain areas. Before disaster, the black pixels account for 15.1641% of the whole image, while this ratio decreases to 13.0785% after the disaster, i.e., the dominant scattering mechanism of about 2.0856% area of the scene is changed from double-bounce scattering to surface scattering. As shown in Figure 7, the change mainly arises in the urban area like the Ishinomaki City and Higashi-Matsushima City, in the land area like the town of Kamaya, as well as in the water area like the Mangokuura Sea, Onagawa Bay, and Kitakami River. This further provides us a consistently quantitative evaluation of the damages. All these demonstrate the importance and value of polarimetric microwave remote sensing technique in the monitoring of tsunami/earthquake damages.

Figure 7.

Binary display of the branch condition BC extracted from (a) pre- and (b) post-event ALOS-PALSAR datasets. The white pixels correspond to BC>0, while the black pixels denote BC0.

Singh et al. [28] indicated that G4U could enhance double-bounce scattering over urban area while strengthen surface scattering contribution over water and land area. This establishes G4U the state-of-the-art four-component scattering power decomposition and enables its wide application to the remote sensing of forestry, agriculture, wetland, snow, glaciated terrain, earth surface, manmade target, environment, and damages caused by earthquake, tsunami, and landslide [29, 30]. Nevertheless, the rigorous derivation in Eq. (21) validates that G4U cannot always enhance the double-bounce scattering nor strengthen the surface scattering power unless we adaptively integrate G4U and its duality, i.e., DG4U, for EG4U based on another boundary condition BC1. As expressed in Eq. (27), G4U is selected only when BC1>0; otherwise, we should turn to DG4U. The binary images Figure 8(a) and (b) further show the pre- and post-event BC1, respectively, where the white pixels (i.e., BC1>0) indicate the area where G4U operates and the black pixels (i.e., BC10) give the area where DG4U operates. The white pixels account for 46.4260% of the pre-event image, which conveys that G4U achieves better result than S4R only for 46.4260% area. As for the rest 53.5740% area, we should resort to DG4U for improvement. The ratio of white pixels increases to 49.5247% after the disaster. Nevertheless, there are still half a little more areas where G4U will underestimate the surface or double-bounce scattering. If we adopt G4U in this area to evaluate damages caused by tsunami/earthquake, the reduced double-bounce scattering from G4U may lead to the underestimation of building scale and overestimation of damage level. EG4U can adaptively increase the surface scattering or double-bounce scattering. Hence, it definitely improves the competence and performance of G4U in the remote sensing of damages caused by earthquake/tsunami.

Figure 8.

Binary display of the branch condition BC1 extracted from (a) pre- and (b) post-event ALOS-PALSAR datasets. The white pixels correspond to BC1>0, while the black pixels denote BC10.

Advertisement

5. Conclusion

Flooding is the main contributor to the severe damages in the great Tohoku tsunami/earthquake. It destroyed the buildings and inundated the lands by the seaside. All these damages present themselves in the polarization domain as the change of the dominant scattering mechanism from double-bounce scattering to surface scattering and in the image domain as the change of pixel color from red to blue. The color-coded scattering power image is very useful and powerful in the qualitative evaluation of damages. The boundary condition BC further enables a nice quantitative evaluation of disaster. The unitary transformation in G4U adds a T13-related but redundant balance equation to the original self-contained equation system. The general solution enables a generalized G4U, while G4U just represents a special form. The strict derivation conveys that G4U cannot always strengthen the double-bounce scattering in urban area nor strengthen the surface scattering in water or land area unless we adaptively combine G4U and its duality for EG4U. Experiment on the ALOS-PALSAR datasets of 2011 great Tohoku tsunami/earthquake demonstrates not only the outperformance of EG4U but also the effectiveness of polarimetric remote sensing in the qualitative monitoring and quantitative evaluation of tsunami/earthquake damages. Efficient and accurate monitoring and assessment are of crucial importance for the fast response, management, and mitigation of the disasters. The all-day and all-weather working capacity is a significant advantage of microwave remote sensing. Polarimetric remote sensing is an effective technique in the discrimination and recognition of ground objects.

Advertisement

Acknowledgments

This work was supported in part by the National Natural Science Foundation of China under Grant No. 41871274 and No. 61971402 and by the Strategic High-Tech Innovation Fund of Chinese Academy of Sciences under Grant CXJJ19B10.

Advertisement

Conflict of interest

The authors declare no conflict of interest.

Advertisement

Notes

Sections 2 and 3 of this chapter are extracted from a journal paper of the authors submitted to IEEE Transactions on Geoscience and Remote Sensing on June 07, 2017. The paper is still under review at the time of publication of this chapter. For more details about the paper, please refer to Reference [30].

References

  1. 1. Mokhtari M, editor. Tsunami-A Growing Disaster. London: IntechOpen; 2011. p. 232. DOI: 10.5772/922
  2. 2. Morner N-A, editor. The Tsunami Threat–Research and Technology. London: IntechOpen; 2011. p. 714. DOI: 10.5772/573
  3. 3. Mokhtari M, editor. Tsunami. London: IntechOpen; 2016. p. 164. DOI: 10.5772/61999
  4. 4. Marghany M, editor. Advanced Remote Sensing Technology for Synthetic Aperture Radar Applications, Tsunami Disasters, and Infrastructure. London: IntechOpen; 2019. p. 167. DOI: 10.5772/intechopen.78525
  5. 5. Li D, Zhang Y. Adaptive model-based classification of PolSAR data. IEEE Transactions on Geoscience and Remote Sensing. 2018;56(12):6940-6955. DOI: 10.1109/TGRS.2018.2845944
  6. 6. Li D, Zhang Y. Random similarity-based entropy/alpha classification of PolSAR data. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing. 2017;10(12):5712-5723. DOI: 10.1109/JSTARS.2017.2748234
  7. 7. Li D, Zhang Y. Unified Huynen phenomenological decomposition of radar targets and its classification applications. IEEE Transactions on Geoscience and Remote Sensing. 2016;54(2):723-743. DOI: 10.1109/TGRS.2015.2464113
  8. 8. Li D, Zhang Y. Random similarity between two mixed scatterers. IEEE Geoscience and Remote Sensing Letters. 2015;12(12):2468-2472. DOI: 10.1109/LGRS.2015.2484383
  9. 9. Lee J-S, Pottier E. Polarimetric Radar Imaging: From Basics to Applications. Boca Raton: CRC Press; 2009. p. 422. DOI: 10.1201/9781420054989
  10. 10. Cloude SR. Polarisation Applications in Remote Sensing. Oxford: Oxford University Press; 2010. p. 453. DOI: 10.1063/1.3502550
  11. 11. van Zyl JJ, Kim Y. Synthetic Aperture Radar Polarimetry. Hoboken: John Wiley & Sons, Inc.; 2011. p. 312. DOI: 10.1002/9781118116104
  12. 12. Yamaguchi Y. Disaster monitoring by fully polarimetric SAR data acquired with ALOS-PALSAR. Proceedings of the IEEE. 2012;100(10):2851-2860. DOI: 10.1109/JPROC.2012.2195469
  13. 13. Sato M, Chen S-W, Satake M. Polarimetric SAR analysis of tsunami damage following the march 11, 2011 East Japan earthquake. Proceedings of the IEEE. 2012;100(10):2861-2875. DOI: 10.1109/JPROC.2012.2200649
  14. 14. Watanabe M, Motohka T, Miyagi Y, Yonezawa C, Shimada M. Analysis of urban areas affected by the 2011 off the Pacific coast of Tohoku earthquake and tsunami with L-band SAR full-polarimetric mode. IEEE Geoscience and Remote Sensing Letters. 2012;9(3):472-476. DOI: 10.1109/LGRS.2011.2182030
  15. 15. Li X, Zhang L, Guo H, Sun Z, Liang L. New approaches to urban area change detection using multitemporal RADARSAT-2 polarimetric synthetic aperture radar (SAR) data. Canadian Journal of Remote Sensing. 2012;38(3):253-255. DOI: 10.5589/m12-018
  16. 16. Singh G, Yamaguchi Y, Boerner W-M, Park S-E. Monitoring of the March 11, 2011, off-Tohoku 9.0 earthquake with super-tsunami disaster by implementing fully polarimetric high-resolution POLSAR techniques. Proceedings of the IEEE. 2013;101(3):831-846. DOI: 10.1109/JPROC.2012.2230311
  17. 17. Chen S-W, Sato M. Tsunami damage investigation of built-up areas using multitemporal spaceborne full polarimetric SAR images. IEEE Transactions on Geoscience and Remote Sensing. 2013;51(4):1985-1997. DOI: 10.1109/TGRS.2012.2210050
  18. 18. Li N, Wang R, Deng Y, Liu Y, Wang C, Balz T, et al. Polarimetric response of landslides at X-band following the Wenchuan earthquake. IEEE Geoscience and Remote Sensing Letters. 2014;11(10):1722-1726. DOI: 10.1109/LGRS.2014.2306820
  19. 19. Chen S-W, Wang X-S, Sato M. Urban damage level mapping based on scattering mechanism investigation using fully polarimetric SAR data for the 3.11 East Japan earthquake. IEEE Transactions on Geoscience and Remote Sensing. 2016;54(12):6919-6929. DOI: 10.1109/TGRS.2016.2588325
  20. 20. Ji Y, Sumantyo JTS, Chua MY, Waqar MM. Earthquake/tsunami damage level mapping of urban areas using full polarimetric SAR data. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing. 2018;11(7):2296-2309. DOI: 10.1109/JSTARS.2018.2822825
  21. 21. Zhai W, Huang C, Pei W. Two new polarimetric feature parameters for the recognition of the different kinds of buildings in earthquake-stricken areas based on entropy and eigenvalues of PolSAR decomposition. Remote Sensing. 2018;10(10):1613. DOI: 10.3390/rs10101613
  22. 22. Cloude SR, Pottier E. A review of target decomposition theorems in radar polarimetry. IEEE Transactions on Geoscience and Remote Sensing. 1996;34(2):498-518. DOI: 10.1109/36.485127
  23. 23. Freeman A, Durden SL. A three-component scattering model for polarimetric SAR data. IEEE Transactions on Geoscience and Remote Sensing. 1998;36(3):963-973. DOI: 10.1109/36.673687
  24. 24. Yamaguchi Y, Moriyama T, Ishido M, Yamada H. Four-component scattering model for polarimetric SAR image decomposition. IEEE Transactions on Geoscience and Remote Sensing. 2005;43(8):1699-1706. DOI: 10.1109/TGRS.2005.852084
  25. 25. Zhu F, Zhang Y, Li D. A novel deorientation method in PolSAR data processing. Remote Sensing Letters. 2016;7(11):1083-1092. DOI: 10.1080/2150704X.2016.1217438
  26. 26. Yamaguchi Y, Sato A, Boerner W-M, Sato R, Yamada H. Four-component scattering power decomposition with rotation of coherency matrix. IEEE Transactions on Geoscience and Remote Sensing. 2011;49(6):2251-2258. DOI: 10.1109/TGRS.2010.2099124
  27. 27. Sato A, Yamaguchi Y, Singh G, Park S-E. Four-component scattering power decomposition with extended volume scattering model. IEEE Geoscience and Remote Sensing Letters. 2012;9(2):166-170. DOI: 10.1109/LGRS.2011.2162935
  28. 28. Singh G, Yamaguchi Y, Park S-E. General four-component scattering power decomposition with unitary transformation of coherency matrix. IEEE Transactions on Geoscience and Remote Sensing. 2013;51(5):3014-3022. DOI: 10.1109/TGRS.2012.2212446
  29. 29. Li D, Zhang Y, Liang L. A Concise Survey of G4U. 2019. Available from: https://arxiv.org/ftp/arxiv/papers/1910/1910.14323.pdf [Accessed: 28 December 2019]
  30. 30. Li D, Zhang Y, Liang L. A mathematical extension to the general four-component scattering power decomposition with unitary transformation of coherency matrix. IEEE Transactions on Geoscience and Remote Sensing. 2020:1-18. Under review
  31. 31. Sinclair G. The transmission and reception of elliptically polarized waves. Proceedings of the IRE. 1950;32(2):148-151. DOI: 10.1109/JRPROC.1950.230106
  32. 32. Huynen JR. Phenomenological theory of radar targets [thesis]. Delft: Delft University of Technology; 1970
  33. 33. Cloude SR, Pottier E. An entropy based classification scheme for land applications of polarimetric SAR. IEEE Transactions on Geoscience and Remote Sensing. 1997;35(1):68-78. DOI: 10.1109/36.551935
  34. 34. Wikipedia. 2011 Tohoku Earthquake and Tsunami. 2019. Available from: https://en.wikipedia.org/wiki/2011_T%C5%8Dhoku_ earthquake_and_tsunami. [Accessed: 28 December 2019]
  35. 35. National Police Agency of Japan. Police Countermeasures and Damage Situation associated with 2011 Tohoku district-off the Pacific Ocean Earthquake. 2019. Available from: https://www.npa.go.jp/news/other/earthquake2011/pdf/higaijokyo_e.pdf [Accessed: 28 December 2019]
  36. 36. JAXA Earth Observation Research Center. About ALOS–Overview and Objectives. 2019. Available from: https://www.eorc.jaxa.jp/ALOS/en/about/about_index.htm [Accessed: 28 December 2019]
  37. 37. Li D, Zhang Y. A novel approach for the registration of weak affine images. Pattern Recognition Letters. 2012;33(12):1647-1655. DOI: 10.1016/j.patrec.2012.04.009
  38. 38. Li D, Zhang Y. A fast offset estimation approach for InSAR image subpixel registration. IEEE Geoscience and Remote Sensing Letters. 2012;9(2):267-271. DOI: 10.1109/LGRS.2011.2166752
  39. 39. Li D, Zhang Y. On the appropriate feature for general SAR image registration. Proceedings of SPIE. 2012;8536:8536X. DOI: 10.1117/12.970520
  40. 40. Li D, Zhang Y. The appropriate parameter retrieval algorithm for feature-based SAR image registration. Proceedings of SPIE. 2012;8536:8536Y. DOI: 10.1117/12.970522

Written By

Dong Li, Yunhua Zhang, Liting Liang, Jiefang Yang and Xun Wang

Submitted: 25 October 2019 Reviewed: 17 January 2020 Published: 20 February 2020