remote sensing
Ocean Remote Sensing with Synthetic Aperture Radar
Edited by Xiaofeng Yang, Xiaofeng Li, Ferdinando Nunziata and Alexis Mouche Printed Edition of the Special Issue Published in Remote Sensing
www.mdpi.com/journal/remotesensing Ocean Remote Sensing with Synthetic Aperture Radar
Special Issue Editors Xiaofeng Yang Xiaofeng Li Ferdinando Nunziata Alexis Mouche
MDPI • Basel • Beijing • Wuhan • Barcelona • Belgrade Special Issue Editors Xiaofeng Yang Chinese Academy of Sciences China
Xiaofeng Li National Oceanic and Atmospheric Administration USA
Ferdinando Nunziata Universita` degli Studi di Napoli Parthenope Italy
Alexis Mouche Laboratoire dOceanographie´ Physique et Spatiale, Ifremer France
Editorial Office MDPI AG St. Alban-Anlage 66 Basel, Switzerland
This edition is a reprint of the Special Issue published online in the open access journal Remote Sensing (ISSN 2072-4292) in 2017 (available at: http://www.mdpi.com/journal/ remotesensing/special issues/ocean rs SAR).
For citation purposes, cite each article independently as indicated on the article page online and as indicated below:
Lastname, F.M.; Lastname, F.M. Article title. Journal Name. Year. Article number, page range.
First Edition 2018
Image courtesy of Xiaofeng Yang, Xiaofeng Li, Ferdinando Nunziata and Alexis Mouche
ISBN 978-3-03842-720-9 (Pbk) ISBN 978-3-03842-719-3 (PDF)
Articles in this volume are Open Access and distributed under the Creative Commons Attribution (CC BY) license, which allows users to download, copy and build upon published articles even for commercial purposes, as long as the author and publisher are properly credited, which ensures max- imum dissemination and a wider impact of our publications. The book taken as a whole is c 2018 MDPI, Basel, Switzerland, distributed under the terms and conditions of the Creative Commons li- cense CC BY-NC-ND (http://creativecommons.org/licenses/by-nc-nd/4.0/). Table of Contents
About the Special Issue Editors ...... v
Preface to ”Ocean Remote Sensing with Synthetic Aperture Radar” ...... vii
Xuan Zhou, Jinsong Chong, Haibo Bi, Xiangzhen Yu, Yingni Shi and Xiaomin Ye Directional Spreading Function of the Gravity-Capillary Wave Spectrum Derived from Radar Observations doi: 10.3390/rs9040361 ...... 1
Hui Meng, Xiaoqing Wang, Jinsong Chong, Xiangfei Wei and Weiya Kong Doppler Spectrum-Based NRCS Estimation Method for Low-Scattering Areas in Ocean SAR Images doi: 10.3390/rs9030219 ...... 15
Weiya Kong, Jinsong Chong and Hong Tan Performance Analysis of Ocean Surface Topography Altimetry by Ku-Band Near-Nadir Interferometric SAR doi: 10.3390/rs9090933 ...... 37
Xiangguang Leng, Kefeng Ji, Shilin Zhou and Huanxin Zou Azimuth Ambiguities Removal in Littoral Zones Based on Multi-Temporal SAR Images doi: 10.3390/rs9080866 ...... 52
He Wang, Jingsong Yang, Alexis Mouche, Weizeng Shao, Jianhua Zhu, Lin Ren and Chunhua Xie GF-3 SAR Ocean Wind Retrieval: The First View and Preliminary Assessment doi: 10.3390/rs9070694 ...... 66
Xiuzhong Li, Biao Zhang, Alexis Mouche, Yijun He and William Perrie Ku-Band Sea Surface Radar Backscatter at Low Incidence Angles under Extreme Wind Conditions doi: 10.3390/rs9050474 ...... 78
Lanqing Huang, Bin Liu, Xiaofeng Li, Zenghui Zhang, Wenxian Yu Technical Evaluation of Sentinel-1 IW Mode Cross-Pol Radar Backscattering from the Ocean Surface in Moderate Wind Condition doi: 10.3390/rs9080854 ...... 92
Lizhang Zhou, Gang Zheng, Xiaofeng Li, Jingsong Yang, Lin Ren, Peng Chen, Huaguo Zhang and Xiulin Lou An Improved Local Gradient Method for Sea Surface Wind Direction Retrieval from SAR Imagery doi: 10.3390/rs9070671 ...... 113
Yi Yu, Xiaofeng Yang, Weimin Zhang, Boheng Duan, Xiaoqun Cao and Hongze Leng Assimilation of Sentinel-1 Derived Sea Surface Winds for Typhoon Forecasting doi: 10.3390/rs9080845 ...... 130
Boheng Duan, Weimin Zhang, Xiaofeng Yang, Haijin Dai, and Yi Yu Assimilation of Typhoon Wind Field Retrieved from Scatterometer and SAR Based on the Huber Norm Quality Control doi: 10.3390/rs9100987 ...... 147
iii Weizeng Shao, Jing Wang, Xiaofeng Li and Jian Sun An Empirical Algorithm for Wave Retrieval from Co-Polarization X-Band SAR Imagery doi: 10.3390/rs9070711 ...... 165
Pengzhen Chen, Lei Liu, Xiaoqing Wang, Jinsong Chong, Xin Zhang and Xiangzhen Yu Modulation Model of High Frequency Band Radar Backscatter by the Internal Wave Based on the Third-Order Statistics doi: 10.3390/rs9050501 ...... 179
Olga Lavrova and Marina Mityagina Satellite Survey of Internal Waves in the Black and Caspian Seas doi: 10.3390/rs9090892 ...... 194
Miao Kang, Kefeng Ji, Xiangguang Leng and Zhao Lin Contextual Region-Based Convolutional Neural Network with Multilayer Fusion for SAR Ship Detection doi: 10.3390/rs9080860 ...... 221
Yichang Chen, Gang Li, Qun Zhang and Jinping Sun Refocusing of Moving Targets in SAR Images via Parametric Sparse Representation doi: 10.3390/rs9080795 ...... 235
Tingting Jin, Xiaolan Qiu, Donghui Hu and Chibiao Ding An ML-Based Radial Velocity Estimation Algorithm for Moving Targets in Spaceborne High- Resolution and Wide-Swath SAR Systems doi: 10.3390/rs9050404 ...... 250
Ji-Wei Zhu, Xiao-Lan Qiu, Zong-Xu Pan, Yue-Ting Zhang and Bin Lei An Improved Shape Contexts Based Ship Classification in SAR Images doi: 10.3390/rs9020145 ...... 267
Yongfeng Cao, Linlin Xu and David Clausi Exploring the Potential of Active Learning for Automatic Identification of Marine Oil Spills Using 10-Year (2004–2013) RADARSAT Data doi: 10.3390/rs9101041 ...... 278
Shuangshang Zhang, Qing Xu, Quanan Zheng and Xiaofeng Li Mechanisms of SAR Imaging of Shallow Water Topography of the Subei Bank doi: 10.3390/rs9111203 ...... 298
Xiaolin Bian, Yun Shao, Wei Tian, Shiang Wang, Chunyan Zhang, Xiaochen Wang and Zhixin Zhang Underwater Topography Detection in Coastal Areas Using Fully Polarimetric SAR Data doi: 10.3390/rs9060560 ...... 319
Wensheng Wang, Martin Gade and Xiaofeng Yang Detection of Bivalve Beds on Exposed Intertidal Flats Using Polarimetric SAR Indicators doi: 10.3390/rs9101047 ...... 335
iv About the Special Issue Editors
Xiaofeng Yang received his B.S. degree in environmental science from Sichuan University, Chengdu, China, in 2005, and his Ph.D. degree in cartography and geographic information systems from the Institute of Remote Sensing Applications (IRSA), Chinese Academy of Sciences (CAS), Beijing, China, in 2010. From 2010 to 2015, he was an Associate Professor with the Institute of Remote Sensing and Digital Earth, CAS. He is currently a Full Professor with the State Key Laboratory of Remote Sensing Science, RADI, CAS. His research interests include satellite oceanography, synthetic aperture radar image processing, and marine atmospheric boundary layer process studies. Dr. Yang serves as an Editor for MDPIs Remote Sensing. He is an IEEE Senior Member, and the Secretary of Technical Committee on Space Earth Science, Chinese Society of Space Research. He has also served as the Assistant Chief Scientist of the Chinese Water Circle Observation Mission.
Xiaofeng Li, Scientist at NOAA, received his B.S. degree in optical engineering from Zhejiang Uni- versity, China in 1985 and his Ph.D. degree in physical oceanography from the North Carolina State University, USA in 1997. He is the author of more than 100 peer-reviewed publications covering the topics in remote sensing observation and theoretical/numerical model studies of various types of oceanic and atmospheric phenomena, satellite image processing, ocean surface oil spill and target detection/classification with multi-polarization SAR, and development of sea surface temperature algorithms. Dr. Li currently serves as an Associate Editor of IEEE Transactions on Geoscience and Remote Sensing, is an Associate Editor of the International Journal of Remote Sensing, and the Ocean Section Editor-in-Chief of Remote Sensing. He is an Editorial Board Member of the International Journal of Digital Earth, Big Earth Data, and CAAI Transactions on Intelligence Technology.
Ferdinando Nunziata, PhD, was born in Italy in 1982. He received his B.Sc. and M.Sc. degrees (summa cum laude) in telecommunications engineering and his Ph.D. degree (curriculum electro- magnetic fields) from the Universita` degli Studi di Napoli Parthenope, Napoli, Italy, in 2003, 2005, and 2008, respectively. Since 2010, he has been an Assistant Professor of electromagnetic fields with the Universita` degli Studi di Napoli Parthenope. He authored/coauthored more than 60 peer-reviewed journal papers that deal with applied electromagnetics.
Alexis Mouche, received his Ph.D. degree in ocean remote sensing from the Universite de Versailles Saint-Quentin, Versailles, France, in 2005. He is currently a Senior Research Scientist with the Labo- ratoire dOceanographie´ Physique et Spatiale, Institut Francais de Recherche pour lExploitation de la Mer, Plouzane, France. His research interests include the interaction of electromagnetic and oceanic waves for ocean remote sensing applications.
v
Preface to “Ocean Remote Sensing with Synthetic Aperture Radar”
The oceans covers approximately 71% of the Earth’s surface, 90% of the biosphere and contains 97% of Earth’s water. In 1978, NASA launched the first SeaSat satellite, primilary aiming at ocean observations and the microwave synthetic aperture radar (SAR) was one of four instruments. Since then, the global oceans have been observed on SAR images, which has a high resolution (<100 m) and a large swath (450 km for ScanSAR mode images). The microwave SAR can image the ocean surface in all weather conditions and day or night. An increasing number of SAR satellites have become available since the early 1990s, such as the ERS-1/-2 and Envisat satellites, the Radarsat-1/-2 satellites, the COSMO- SkyMed constellation, TerraSAR-X and TanDEM-X, the Gaofen-3, among others. Recently, the European Space Agency lauched a new generation of SAR satellites (Sentinel-1A in 2014 and Sentinel-1B in 2016). This operational SAR mission, for the first time, provides researchers with free and open SAR images necessary to carry out broader and deeper investigation of the global oceans. SAR remote sensing of ocean and coastal monitoring has become a research hotspot in geoscience and remote sensing. This book—Progress in SAR Oceanography—provides an update of the state-of-the- science research on ocean remote sensing with SAR. Overall, the book presents a variety of marine applications in ocean research topics such as, oceanic surface and internal waves characteristics studies, high-resolution sea surface wind retrieval, shallow-water bathymetry mapping, oil spill detection, coastline and inter-tidal zone classification, ship and other man-made objects detection, as well as remotely sensed data assimilation. The book is aimed at a wide audience, ranging from graduate students, university faculty members, scientists to policy makers and managers.
Xiaofeng Yang, Xiaofeng Li, Ferdinando Nunziata, Alexis Mouche Special Issue Editors
vii
remote sensing
Article Directional Spreading Function of the Gravity-Capillary Wave Spectrum Derived from Radar Observations
Xuan Zhou 1,*, Jinsong Chong 1, Haibo Bi 2,3, Xiangzhen Yu 4, Yingni Shi 5 and Xiaomin Ye 6 1 National Key Laboratory of Microwave Imaging Technology, Institute of Electronics, Chinese Academy of Sciences, Beijing 100190, China; [email protected] 2 Key Laboratory of Marine Geology and Environment, Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266061, China; [email protected] 3 Laboratory for Marine Geology, Qingdao National Laboratory for Marine Science and Technology, Qingdao 266200, China 4 Shanghai Radio Equipment Research Institute, Shanghai 200090, China; [email protected] 5 College of Information Science & Engineering, Ocean University of China, Qingdao 266100, China; [email protected] 6 National Satellite Ocean Application Service, Beijing 100081, China; [email protected] * Correspondence: [email protected]; Tel.: +86-134-666-60511
Academic Editors: Xiaofeng Li, Ferdinando Nunziata, Alexis Mouche, Raphael M. Kudela and Prasad S. Thenkabail Received: 3 December 2016; Accepted: 1 April 2017; Published: 12 April 2017
Abstract: Directional spreading function of the gravity-capillary wave spectrum can provide the high-wavenumber wave energy distribution among different directions on the sea surface. The existing directional spreading functions have been mainly developed for the low-wavenumber gravity wave with buoy data. In this paper, we use radar observations to derive the directional spreading function of the gravity-capillary wave spectrum, which is expressed as the second-order Fourier series expansion. So far the standard form of the second-order harmonic coefficient has not been proposed to correctly unify the gravity and gravity-capillary wave. Our strategy is to introduce a correcting term to replace the inaccurate gravity-capillary spectral component in Elfouhaily’s directional spreading function. The second-order harmonic coefficient at L, C and Ku band calculated by the radar observation is used to fit the correcting term to obtain one at the full gravity-capillary wave region. According to our proposed the directional spreading function, there is a spectral region between the gravity and gravity-capillary range where it signifies the negative upwind–crosswind asymmetry at low and moderate speed range. And this is not reflected by the previous models, but has been confirmed by radar observations. The Root Mean Square Difference (RMSD) of the proposed second-order harmonic coefficient versus the radar-observed one at L, C band Ku band is 0.0438, 0.0263 and 0.0382, respectively. The overall bias and RMSD are −0.0029 and 0.0433 for the whole second-order harmonic coefficient range, respectively. The result verifies the accuracy of the proposed directional spreading function at L, C band Ku band.
Keywords: directional spreading function; gravity-capillary wave; radar observations
1. Introduction The gravity-capillary wave plays an important role in air-sea interaction because it affects the mass, momentum and energy flux through the air-sea interface. The wind-induced turbulence transfers wind energy from the atmosphere to the gravity-capillary wave by the friction at the interface, and then because the phase speed of the gravity-capillary wave is less than one of the gravity waves, its energy
Remote Sens. 2017, 9, 361 1 www.mdpi.com/journal/remotesensing Remote Sens. 2017, 9, 361 is pillaged by the gravity wave and wind wave will grow [1,2]. The energy propagation is colinear with gravity crests propagation in absence of currents and internal gravity waves. The directional spectrum is used to describe the gravity-capillary waves, which can give the wave energy distribution among different directions on the sea surface. There are two main ways to observe the azimuthal behavior of the directional spectrum: in situ measurement and remote sensing measurement. Buoy and its array provided two main types of directional spreading functions of ocean wave spectrum. The cosine-shape spreading function was first proposed by Longuet-Higgines et al. [3] according to the motion of a flotation buoy. Mitsuyasu et al. [4] estimated the spreading parameter of the cosine-shape function using a cloverleaf buoy data. Hasselmann et al. [5] improved the cosine-shape functions using wave data collected by pitch-and-roll buoys. The sech-shape spreading function was advanced by Donelan et al. [6] using data collected from a 14-element wave gauge array because they found that the distribution of wave energy in the direction transverse to the main wave direction behaved like a hyperbolic secant. The above-mentioned directional spreading functions measured by buoy are suitable for the gravity wave spectrum. However, the wavelength of the gravity-capillary wave is too short for buoy and its array to measure the azimuthal behavior of the gravity-capillary wave. Radar is the important method to measure the gravity-capillary wave spectrum due to the Bragg resonance scattering between electromagnetic waves and gravity-capillary waves [7–9]. Apel [10], Caudal et al. [11] and Liu et al. [12] extended the sech-shape spreading function to the gravity-capillary domain using radar data. However, the sech-shape spreading function cannot explain the angular scattering behavior of radar data due to its noncentrosymmetric property. Guissard [13] pointed out that the directional spreading function of gravity-capillary wave spectrum should contain only even harmonics if expressed as a Fourier series. And then it was also applied to a directional spectrum by Elfouhaily et al. [14] and Hwang et al. [15]. Unfortunately, when these directional spreading functions are used to demonstrate scattering properties of the sea surface, there are obvious differences between theoretical calculations and radar observations. For example, the L-band negative upwind-crosswind asymmetry [16–18] of backscatter at low wind speeds is not explained by the existing directional spreading functions. The radar observation is related to the directional spectrum by sea surface backscatter model, such as the Two-Scale Model (TSM). However, the double integrals in TSM are very inconvenient to calculate the directional spreading function. For the VV polarization, the solutions of Small-Perturbation Method (SPM) are approximately equal to the TSM solutions in intermediate incidence angles. Therefore, the SPM is used to relate the directional spectrum to the GMFs at VV polarization in intermediate incidence angles. The new directional spreading function at the full gravity-capillary wave region is derived from the L-, C- and Ku-band one calculated by the SPM. This paper is organized as follows: Section 2 describes the L-, C- and Ku-band GMFs. The detail of methodology is given in Section 3. Section 4 validate the directional spreading function by radar observation from SMAP SAR, METOP-A ASCAT and QuikSCAT SeaWinds-1. The whole paper is discussed and concluded in Sections 5 and 6, respectively.
2. Data Description In this paper, the L-, C- and Ku-band geophysical model functions (GMFs), which empirically describe the backscattering properties of sea surface, are chosen to serve as a proxy of radar observations to derive a new directional spreading function of the gravity-capillary wave spectrum. The combination of the L, C and Ku bands provide a good coverage of the gravity-capillary wave spectrum for the wavenumber ranging from 25 to 500 rad/m. The following derivation is based on L-band GMF [18], C-band CMOD5 GMF [19] and Ku-band NSCAT2 GMF [20]. The L-band GMF and CMOD5 GMF express the NRCSs as second-order cosine harmonic functions of the radar-observed azimuthal angle with the analytical functions [21]. The NSCAT2 GMF is given as the lookup table with
2 Remote Sens. 2017, 9, 361 respect to the NRCS, the polarization, the 10-m-height wind speed, the relative wind direction and the incidence angle. ◦ Figure 1 shows the contour plots of SMAP SAR GMF, CMOD5 GMF and NSCAT2 GMF in 40 incidence angles. The contour line of each wind speed is symmetric around the wind direction. At C ◦ ◦ and Ku band, the maxima of the contour line occur in the upwind (0 ) and downwind (180 ) directions ◦ ◦ and minima in the crosswind (90 or 270 ) directions. However, the contour lines of low wind speeds overlap ones of moderate and high wind speeds at L band, which is obviously different from C- and Ku-band pattern. This is not explained by the existing directional spreading functions.
(a)
(b)
(c)
◦ Figure 1. The SMAP SAR GMF (a); CMOD5 GMF (b) and NSCAT2 GMF (c)in40 incidence angles.
3. Methodology The directional spreading function of the gravity-capillary wave spectrum is related to the radar observation through sea surface backscatter model. We first introduce the directional wave spectrum and propose a basic form of directional spreading function. And then TSM and SPM, which are two basic kinds of sea surface backscatter model, are described and compared. It is generally known that TSM is more suitable for a realistic sea surface due to introduce the double integrals to describe the sea surface tilting effect. However, the double integrals in TSM are very inconvenient to be used to calculate the directional spreading function with radar observations. Fortunately, the solutions of SPM are approximately equal to the TSM solutions in intermediate incidence angles at VV polarization.
3 Remote Sens. 2017, 9, 361
Finally, we derive the directional spreading function from SPM and calculate its parameters with L-, C- and Ku-band GMFs.
3.1. Basic Form of Directional Spreading Function The directional wave spectrum can provide the directional distribution of ocean wave energy on the sea surface. With the increase of the quality and quantity of available data, more and more directional wave spectrums have been proposed [22,23]. In most cases the directional wave spectrum Ψ(k, φ) can be described as a function of both the wave wavenumber and the wave direction relative to the wind as follows: Ψ(k, φ) = ϕ(k)·D(k, φ) (1) where k is the wave wavenumber, φ is the wave direction relative to the wind, ϕ(k) is the omnidirectional wave spectrum and D(k, φ) is the directional spreading function defined as:
Ψ(k, φ) D(k, φ) = (2) 2π Ψ( φ) φ 0 k, d
If the directional wave spectrum is expressed as a Fourier series, the directional spreading function should contain only even harmonics:
∞ ( φ) = 1 [ + ( φ)] D k, π 1 ∑ a2n cos 2n (3) 2 n=1 where a2n is the coefficient of even harmonics. In fact the Fourier series expansion is usually truncated to second order [10]: 1 D(k, φ) = [1 + Δ(k) cos(2φ)] (4) 2π where Δ(k) is the second-order harmonic coefficient and the function of both the wavenumber and the wind speed. Unfortunately, up till now the shape of the directional spreading function has been a controversial issue, and the standard form of Δ(k) has not been given to correctly unify the gravity and gravity-capillary wave [11]. Here, the form of the directional spreading function from Elfouhaily’s spectrum is used and Δ(k) is expressed as: Δ( ) = { + 2.5 + ( )2.5} k tanh a0 ap c/cp am cm/c (5)
∗ where a0 and ap are both constants, am is the function of u /cm, c is the phase speed, cp is the phase ∗ 2.5 speed of the dominant long wave and u is the friction velocity at the sea surface. ap c/cp and 2.5 am(cm/c) in Equation (5) are related to the directionality of the gravity wave and the gravity-capillary wave, respectively. However, Elfouhaily’s spectrum cannot correctly reflect the directionality of the gravity-capillary wave because radar data is excluded from its development. A correcting term is introduced to replace the gravity-capillary spectral component in the directional spreading function of Elfouhaily’s spectrum: Δ( ) = { 2.5 + δ( )} k tanh ap c/cp k, U10 (6) δ where ap is equal to 4, U10 is the 10-m-height wind speed and is a correction factor which is a function of wavenumber and wind speed.
3.2. Comparison and Selection of Sea Surface Backscatter Model Sea surface backscatter model can describe the relation between radar observation and directional wave spectrum. The SPM and TSM, which are two basic approaches to calculate ocean-surface scattering, are suitable for the small-scale surface and the tilted small-scale surface, respectively.
4 Remote Sens. 2017, 9, 361
3.2.1. Small-Perturbation Method According to electromagnetic scattering perturbation theory, the Normalized Radar Cross Section (NRCS) of the gravity-capillary wave surface without regard to the tilting effect can be calculated by the first-order SPM [24]: σ (θ) = π 4 4 θ (θ) 2Ψ( θ φ) 0pq 16 kR cos gpq 2kR sin , (7) σ where 0 is the NRCS, the indices p and q represent transmitting and receiving polarizations, = π λ λ θ respectively; kR is the radar wavenumber, kR 2 / , is the radar wavelength; is the incidence angle; gpq(θ) is the first-order scattering coefficient.
3.2.2. Two-Scale Model In fact, the gravity-capillary waves are tilted by the gravity waves of sea surface. The tilting effect θ θ modifies the incidence angle referenced to a horizontal surface as the local angle i. Accounting for the sea surface tilting effect, the NRCS is calculated by TSM [24]: