利用GF-1和MODIS準同步光學遙感圖像反演內波參數的研究?

梅 源， 張旭東， 孫麗娜， 孟俊敏， 王 晶? ?

(1.中國海洋大學信息科學與工程學院，山東 青島 266100； 2.國家海洋局第一海洋研究所，山東 青島 266061)

利用GF-1和MODIS準同步光學遙感圖像反演內波參數的研究?

梅 源1， 張旭東1， 孫麗娜2， 孟俊敏2， 王 晶1? ?

(1.中國海洋大學信息科學與工程學院，山東 青島 266100； 2.國家海洋局第一海洋研究所，山東 青島 266061)

高分一號(GF-1)光學遙感衛星的發射為海洋內波的研究提供了精細的觀測資料，利用其圖像高空間分辨率的優勢，結合高時間分辨率的中分辨率成像光譜儀(MODIS)的匹配數據，能夠對南中國海東沙島附近的內波開展特征參量的反演。本文處理了GF-1和MODIS的準同步圖像，選取東沙島附近的相鄰時刻的同一條內波進行分析，提取其空間位移，計算出內孤立波傳播的群速度。分析GF-1和MODIS圖像上同一條內波的灰度剖面，計算了明暗條紋間距，結合當地水深和溫鹽數據，應用非線性薛定諤方程反演內孤立波的振幅，探究大陸架附近內孤立波的振幅演變。研究結果表明：內孤立波的群速度在東沙島附近向西傳播過程中逐漸變??；在相似的溫鹽條件下，遙感圖像中內波的明暗條紋間距和振幅呈現負相關，在深海區水深和振幅呈現正相關關系，在淺海區由于非線性作用增強會出現振幅隨水深變小反而增大的現象。

海洋內波；光學遙感；群速度；振幅；準同步

目前，遙感是大范圍觀測內波的重要手段。我國高分專項近期成功發射的高分一號(以下簡稱GF-1)光學遙感衛星的空間分辨率達到米級。而以中分辨率成像光譜儀(以下簡稱MODIS)為代表的中分辨率光學遙感器，有著圖像覆蓋范圍大、時間分辨率高、數據公開免費的優勢。將GF-1和MODIS兩者相結合，為深入研究內波精細結構的反演提供了可能。

衛星遙感觀測表明海洋中的內孤立波活動十分活躍，普遍存在于世界各大洋[1]。Alford等[2]根據SAR遙感圖像和實測數據，對南中國海內波的生成和傳播路徑進行長期深入研究，得出南中國海內波產生于呂宋海峽周期性的內潮，其傳播伴隨著能量流的傳遞，并且黑潮影響內波的傳播方向。針對南中國海北部海域的內波動力學特征[3]和大振幅內波能量屬性[4]的研究也從未停止。同時，大量遙感數據統計了該海域的海洋內波時空分布，并且建立模型進行信息反演[5-9]。楊頂田[10]利用中巴資源衛星研究內波在海南島附近海域的分布和傳播方向，反演了內波的振幅、相速度等參數信息。Chen等[11]利用兩層模型和多層模型下的KDV/BO方程的系數對應相等，從而得到了最合適的上下層水深及密度，再結合MODIS圖像反演出內波的振幅。Liu等[12-13]利用MODIS和SAR探測到的中國南海的同一條內波，計算了內波的相速度，并發現內波的相速度與水深成正比。次年，通過Envisat和ERS-2同步衛星相鄰30 min的三對SAR圖像追蹤了13列內孤立波的軌跡，由空間位移和時間間隔計算出內孤立波的相速度，與當地水深和月度氣候學海洋分層數據得出的理論值相吻合。黃曉冬等[14]基于MODIS圖像，對位于南海北部水深3 000 m左右的內波個例進行了研究，他們采用Zheng等[15]提出的Peak-to-Peak方法反演了內波振幅和相速度?？傊?，學者們利用SAR和MODIS等觀測不同海域內波并進行了參數反演。SAR雖然具有全天候、全天時探測能力，但重訪時間過長，而以MODIS為代表的光學遙感傳感器雖然時間分辨率較高，但最好的空間分辨率只有250 m，都無法實現內波精細結構和準連續時間演變的探測。GF-1光學遙感圖像的空間分辨率達到米級，彌補了MODIS空間分辨率的不足，我們以此為切入，結合MODIS的高時間分辨率，選取東沙島附近內波為研究對象進行觀測和參數反演。

本文采用2013年7月—2014年6月的GF-1和MODIS數據，分析南中國海東沙島(20°N～22°N，116°E～118°E)附近的同一條內波，計算內波傳播的群速度并反演內波振幅，探究內波在近海大陸架的特征參數變化規律。

1 東沙島內波光學遙感圖像收集與處理

GF-1是我國“高分系列項目”第一顆高分辨率對地觀測光學遙感衛星，于2013年4月26日發射升空，其上搭載兩臺2 m分辨率全色/8 m分辨率多光譜相機、四臺16 m分辨率多光譜相機，時間分辨率為4 d。GF-1的軌道參數表明，衛星每日于北京時間11時左右穿過南中國海區域，為我們觀測該海域內波提供了有利條件。我們選取16 m分辨率的GF-1圖像對東沙島內波進行研究。

MODIS是搭載在美國Terra和Aqua衛星上的重要光學傳感器，是觀測內波的主要數據來源。MODIS采用掃描式成像，圖像刈幅可達2 330 km，三種空間分辨率為250、500和1 000 m，時間分辨率為1～2 d。Terra衛星于10時30分降軌經過赤道，在其時間段附近可以收集覆蓋南海的大量MODIS圖像。我們選取250m分辨率的MODIS圖像分析東沙島的內波。

本文收集了2013年7月至2014年6月累計一年的東沙島地區(20°N～22°N，116°E～118°E)GF-1和MODIS圖像作為研究對象，由于光學遙感圖像受到天氣影響較大，能夠同時在高分辨率的GF-1圖像和中分辨率的MODIS圖像上追蹤到清晰的內波軌跡，獲得3對(6景)圖像。由于空間位置固定，成像時間相近，可以確定遙感圖像中為同一條內波的傳播。

對GF-1和MODIS圖像進行幾何校正和添加經緯度后，獲取了南海東沙島分別在2014年4月18日，2014年3月20日和2013年9月7日3對GF-1與MODIS圖像，追蹤了6列內波的軌跡，如圖1所示。

2 計算內波傳播的群速度

基于GF-1和MODIS獲得的東沙島內波空間位置信息(見圖1)，從波峰線的推移可以看出內波在東沙島海域從東向西傳播。由于GF-1和MODIS拍攝的時間差Δt是確定的，測量出同一條內波的空間位移Λ，就可以得出內波傳播的群速度Cg。即

Cg=Λ/Δt。

(1)

本文采用的GF-1圖像分辨率為16 m，MODIS圖像分辨率為250 m，GF-1圖像分辨率較MODIS圖像分辨率更高，讀取內波明暗條紋間距中心點經度時更精確。

對于圖1(a)(b)2014年4月18日的內波A和B圖像，GF-1上內波A1和B1成像于03:14:10 UTC，MODIS上內波A2和B2成像于02:45 UTC，(c)(d)2014年3月20日的內波C、D和E圖像，GF-1上內波C1、D1和E1成像于03:05:43 UTC，MODIS上內波C2、D2和E2成像于03:15 UTC，(e)(f)2013年9月7日的內波F圖像，GF-1上內波F1成像于03:18:15 UTC，MODIS上內波F2成像于03:25 UTC，東沙島所在緯20°N區域每個經度值間隔約104.31km，同一條內波的空間位移Λ由準同步的GF-1和MODIS圖像上明暗條紋的中心點經度差乘以一個經度對應的距離得出，計算出東沙島附近內波傳播的群速度(見表1)。

(GF-1圖像(a)(c)(e)與MODIS圖像(b)(d)(f)相匹配對應，A～F表示追蹤的6列內波軌跡。GF-1 images (a) (c) (e) match MODIS images (b) (d) (f), A～F represent paths of 6 internal waves in track.)

圖1 東沙島內波的GF-1和MODIS圖像Fig.1 GF-1 and MODIS images of internal waves near Dongsha Atoll表1 追蹤的6列東沙島附近內波傳播的群速度Table 1 The group velocity of 6 internal waves in track near Dongsha Atoll

Note: ①Internal wave；②Latitude；③The central longitude of bright-dark stripes on GF-1 images；④The central longitude of bright-dark stripes on MODIS images；⑤Spatial displacement；⑥Time difference；⑦The group velocity

由時間可以看出，GF-1和MODIS都于北京時間上午11時左右經過東沙島上空，圖像顯示南中國海內波自東向西傳播，傳播至東沙島附近分裂。按上述方法，得出在東沙島周圍，內波傳播群速度約1～1.5 m/s；東沙島的東部，內波傳播的群速度約2～2.5 m/s。在此海域范圍，Liu等[16]研究2006年5月三處CTD站點的實測數據，計算得出自東向西水深分別為1 100、650和320 m時，內波傳播的群速度依次為2.139 4、1.511 5和1.062 2 m/s的結論。得出的東沙島海域內波傳播速度與該學者的實測結果一致。

基于兩星數據的時間間隔較短，同一條內波在兩景圖像的位置接近，此時測量空間位移的誤差是傳播群速度誤差的主要來源。結合追蹤的6列內波A～F所在位置的水深數據，分析三組不同圖像，東沙島東部水深較深處，內波A、B、C、F的傳播速度接近或大于2 m/s，東沙島周圍水深較淺處，內波D、E的傳播速度小于2 m/s；分析同一組圖像，以2014年3月20日為例，東沙島東部的內波C，東沙島南部的內波D和東沙島北部的內波E水深依次變淺，其傳播速度也依次變小。因此，上述結果表明東沙島附近內波的群速度隨水深變淺而減小。

3 內波振幅的反演

采用高階完全非線性薛定諤方程描述內波[17]：

(2)

式中：α是頻散系數；β是非線性系數；并且引進了兩個修正項；α1是高階頻散系數；β1和β2是高階非線性系數，公式中第2、3項為頻散項，第4、5、6項為非線性項。

采用方程(2)描述內波傳播，結合遙感圖像的成像機理，可以得到如下反演模型[18]。當αβ<0 時，內波半寬度與遙感圖像中內波的亮暗間距的關系為

D=1.32l。

(3)

當αβ>0時，可以得出

D=1.76l。

(4)

式中：間距D代表的是能夠從遙感圖像上直接測量的內波亮暗條紋間距；l代表內波的半寬度，內波的振幅為

(5)

反演內波的振幅需要明暗條紋間距D，頻散系數α和非線性系數β。首先，α和β可由總水深，上層水深，上層平均密度和下層平均密度算出，將從地形高程數據官網Etopo1和溫鹽數據集下載所需經緯度范圍內的水深數據和溫度鹽度數據進行處理可得以上參量；其次，明暗條紋間距D可從GF-1和MODIS圖像上直接提取內波信息；最后，按如上公式(5)可得內波的振幅值A0。

在GF-1和MODIS圖像中選取一個無云的子圖像，將子圖像沿垂直于內波傳播方向取灰度剖面曲線。

((a)GF-1圖像的灰度剖面曲線；(b)MODIS圖像的灰度剖面曲線。(a) Gray level profile curve of GF-1 images; (b) Gray level profile curve of MODIS images.)

圖2 2014-04-18東沙島內波A于20°50′N的灰度剖面曲線
Fig.2 Gray level profile curve of internal wave A near Dongsha Atoll at 20°50′N on April 18, 2014

圖2是2014年4月18日東沙島內波A在緯度20°50′N處提取的灰度剖面曲線。對比GF-1和MODIS圖像的灰度剖面曲線可以看出，GF-1的灰度剖面曲線像素點密集，峰谷之間包含的像素點較多，曲線比較平滑；MODIS的灰度剖面曲線像素點稀疏，峰谷之間包含的像素點較少，曲線比較波折。因此，GF-1提取的明暗條紋間距D較MODIS更精確。圖2(a)GF-1圖像成像于太陽耀斑區，灰度剖面圖先暗后亮，2(b)MODIS圖像成像于非耀斑區，同一條內波灰度剖面圖先亮后暗[19]。

為了反演內波振幅，分別提取圖1中三組內波圖像的當地溫鹽水深數據，得出如下參數。對于圖1(a)(b)的內波A和B圖像，上層水深為40m，上層密度為1 022.2 kg/m3，下層密度為1 024 kg/m3，(c)(d)的內波C、D和E圖像，上層水深為25 m，上層密度為1 022.3 kg/m3，下層密度為1 024 kg/m3，(e)(f)的內波F圖像，上層水深為50 m，上層密度為1 022.4 kg/m3，下層密度為1 023.8 kg/m3，根據上述方法反演的東沙島附近內波振幅結果見表2。

表2 東沙島附近GF-1和MODIS圖像6列內波振幅的反演Table 2 Inversion of 6 internal waves amplitude on GF-1 and MODIS images near Dongsha Atoll

由于GF-1空間分辨率較MODIS更高，圖像灰度剖面圖下像素點更密集，峰、谷的位置更精確，而MODIS讀取的明暗條紋間距比GF-1的偏大。對于準同步的內波B1和B2、D1和D2、F1和F2，MODIS圖像上讀取的明暗條紋間距比GF-1的大100～300 m，而內波A1和A2、C1和C2、E1和E2，MODIS圖像上讀取的明暗條紋間距比GF-1的大400～600 m。分析表2中三組圖像的參數，追蹤的6列內波A～F在GF-1和MODIS圖像上位置接近，其中內波A、D、E和F在相匹配的圖像上水深近似相同，相差不超過20 m，MODIS圖像的明暗條紋間距大于GF-1圖像的明暗條紋間距，MODIS圖像的反演振幅小于GF-1圖像的反演振幅，可以看出內波的明暗條紋間距和振幅呈負相關。分析同一組圖像，以2014年4月18日為例，東沙島東部的內波A和B，明暗條紋間距相差不大，內波A位于內波B的西側，A處水深小于B處，內波A反演振幅小于內波B反演振幅；分析不同組圖像，以2014年3月20日東沙島南部的內波D和2013年9月7日東沙島東部的內波F為例，明暗條紋間距相差也不大，D處水深大于F處，內波D反演振幅卻小于內波F反演振幅?？梢妰炔ǖ恼穹退畹年P系比較復雜，反演振幅決定于頻散系數α和非線性系數β，明暗條紋間距相近情況下，在深海區水深和內波的反演振幅呈正相關，在淺海區非線性作用增強會導致水深變小振幅增大的情況。

4 結語

本文利用2013-07—2014-06的GF-1和MODIS數據，對比分析南中國海東沙島的同一條內波的參數特征。經過大數量GF-1圖像的篩選，結合同時間段的MODIS圖像，得出了東沙島內波傳播的群速度。同時，基于高階完全非線性薛定諤方程，結合內波在遙感圖像上的明暗條紋間距反演振幅，得出東沙島所在的大陸架附近內波振幅。

南中國海的東沙島內波自東向西傳播，其群速度隨水深的變淺而減小。東沙島所在大陸架附近反演的內波振幅量級都為幾十米，這與甘錫林等[20]的遙感統計數據相吻合:在東沙島附近,內波相速度在1.5～2 m/s左右,振幅在20～80 m左右?？偹钕嘟闆r下，內波的明暗條紋間距和反演振幅呈現負相關關系；明暗條紋間距相近情況下，反演振幅決定于頻散系數α和非線性系數β，在南海深水區域總水深和內波的反演振幅呈現正相關關系，在淺水區域非線性作用增強會導致水深變小振幅增大。

由GF-1和MODIS圖像數據得到的內波在東沙島傳播的群速度和反演振幅，可以清晰地看出內波在近海大陸架傳播速度和振幅的變化。由于MODIS遙感圖像的分辨率低導致反演內波振幅偏小，誤差較大。無疑，高分辨率的GF-1光學遙感為海洋內波的研究提供了精細化的觀測資料，大大提高了探測內波的精度。

[1] Jackson C. Internal wave detection using the moderate resolution imaging spectroradiometer (MODIS)[J]. Journal of Geophysical Research, 2007, 112: C11012-1- C11012-13.

[2] Alford M H. The formation and fate of internal waves in the South China Sea[J]. Nature, 2015, 521 (7550): 65-69.

[3] Guo C and Chen X. A review of internal solitary wave dynamics in the northern South China Sea[J]. Progress in Oceanography, 2014, 121: 7-23.

[4] Lien R C. Large-amplitude internal solitary waves observed in the northern South China Sea: properties and energetics[J]. Journal of Physical Oceanography, 2014, 44(4): 1095-1115.

[5] 蔡樹群, 何建玲, 謝皆爍. 近10年來南海孤立內波的研究進展[J]. 地球科學進展, 2011, 26 (7): 703-710.

Cai S Q, He J L and Xie J S. Recent decadal progress of the study on internal solitons in the South China Sea[J]. Advances in Earth Science, 2011, 26(7): 703-710.

[6] 張旭東, 王晶, 魏子章, 孫麗娜, 孟俊敏. 利用光學遙感探測南中國海內波[J]. 激光與光電子學進展, 2015, 52(4): 042802-1-042802-5.

Zhang X D, Wang J, Wei Z Z, Sun L N and Meng J M. Detection of internal solitary wave in the South China Sea using optical remote sensing[J]. Laser & Optoelectronics Progress, 2015, 52(4), 042802-1-042802-5.

[7] Song S Y, Wang J, Meng J M, et al. Nonlinear Schrodinger equation for internal waves in deep sea[J]. Acta Physics Sinica, 2010, 1123-1129.

[8] Wang J, Guo K and Meng J M. Study of the propagation model for large-amplitude internal waves in deep sea[J]. Chinese Journal of Lasers, 2012, 39 (s2): 214004-1-214004-6.

[9] 王晶, 郭凱, 孫美玲, 孟俊敏. 3層模型的內波傳播方程與參數反演[J]. 遙感學報, 2015, 19(2): 188-194.

Wang J, Guo K, Sun M L and Meng J M. Nonlinear schrodinger equation and parametric inversion of internal wave propagation using three-layer model[J]. Journal of Remote Sensing, 2015, 19(2): 188-194.

[10] Yang D T. Investigating internal waves east of the Hainan Island using optical satellite remote sensing data[J]. Journal of Tropical Oceanography, 2009, 28 (5): 29-34.

[11] Chen G Y. Interaction and generation of long‐crested internal solitary waves in the South China Sea[J]. Journal of Geophysical Research, 2011, 116: C06013-1-C06013-7.

[12] Liu B Q. Fusion of SAR and MODIS images for oceanic internal waves tracking in the South China Sea[J]. Proceedings of SPIE, 2013, 8921: 89210L.

[13] Liu B Q. Internal solitary wave propagation observed by tandem satellites[J]. Geophysical Research Letters, 2014, 41(6): 2077-2085.

[14] 黃曉冬, 趙瑋. 基于 MODIS 圖像的內孤立波信息反演——以南海北部深水區為例[J]. 中國海洋大學學報(自然科學版), 2014, 44(7): 19-23.

Huang X D and Zhao W. Information of internal solitary wave extracted from MODIS image: a case in the deep water of northern South China Sea[J]. Periodical of Ocean University of China, 2014, 44(7): 19-23.

[15] Zheng Q N. Theoretical expression for an ocean internal soliton synthetic aperture radar image and determination of the soliton characteristic half width[J]. Journal of Geophysical Research, 2001, 106 (C12): 31415-31423.

[16] Liu A K, Su F C, Hsu M K, et al. Generation and evolution of mode-two internal waves in the South China Sea[J]. Continental Shelf Research, 2013, 59: 18-27.

[17] 馬瑞玲. 深海內波傳播模型及其數值模擬的研究[D]. 青島: 中國海洋大學, 2012.

Ma R L. Study on the Propagation Model for the Deep-Sea Internal Waves and Its Numerical Simulation[D]. Qingdao: Ocean University of China, 2012.

[18] Li X Y, Wang J, Sun M L, et al. Internal wave parameter inversion at Malin Shelf edge based on the nonlinear Schrodinger equation[J]. Applied Mechanics and Materials, 2014, 441: 388-392.

[19] 康健. 光學遙感影像中耀斑區內孤立波信息提取模型[D]. 青島: 中國海洋大學, 2009.

Kang J. Internal Solitary Wave Retrieval Model from Sun Glint of Optical Remote Sensing Images[D]. Qingdao: Ocean University of China, 2009.

[20] 甘錫林, 黃韋艮, 楊勁松, 等. 利用多源遙感衛星數據研究南海內波的時空分布特征[J]. 遙感技術與應用, 2007, 22(2): 242-245.

Gan X L, Huang W G, Yang J S, et al. The study of spatial and temporal distribution characteristics of internal waves in the South China Sea from multi-satellite data[J]. Remote Sensing Technology and Application, 2007, 22(2): 242-245.

StudyonInversionofInternalWaveParameterUsingGF-1andMODISQuasi-SynchronousOpticalRemoteSensingImages

MEI Yuan1, ZHANG Xu-Dong1, SUN Li-Na2, MENG Jun-Min2, WANG Jing1

(1.College of Information Science and Engineering, Ocean University of China, Qingdao 266100, China； 2.First Institute of Oceanography, State Oceanic Administration, Qingdao 266061, China)

“GF-1” optical remote sensing satellite provides detailed data for the study of ocean internal wave. Combine the high spatial resolution images from GF-1 with the matching data of high temporal resolution from moderate-resolution imaging spectroradiometer (MODIS), inversion of the internal wave characteristic parameter near Dongsha Atoll in the South China Sea was studied. In this paper, GF-1 and MODIS quasi-synchronous images with the same internal wave in adjacent time near Dongsha Atoll were processed. Then the spatial displacement was extracted and thus the group velocity of internal solitary wave was calculated. In addition, gray level profiles of the same internal wave in GF-1 and MODIS images were analyzed, and the distance between the positive and negative peaks were calculated. Combining with the local bathymetry and data of temperature and salinity, the evolution of internal solitary wave amplitude near the continental shelf was studied with the application of nonlinear Schrodinger equation for amplitude inversion. Results showed that the group velocity becomes smaller with internal solitary wave near Dongsha Atoll propagating westward. Under the similar condition of temperature and salinity, the distance between the positive and negative peaks in remote sensing images and the amplitude present negative correlation, water depth and amplitude present positive correlation in the deep sea area. In shallow sea, however, amplitude of internal waves will increase with the decreasing depth of water because of the strong nonlinear effect.

ocean internal wave; optical remote sensing; group velocity; amplitude; quasi-synchronous

P731.24，P236

A

1672-5174(2018)02-113-07

10.16441/j.cnki.hdxb.20150373

梅源， 張旭東，孫麗娜， 等. 利用GF-1和MODIS準同步光學遙感圖像反演內波參數的研究[J]. 中國海洋大學學報(自然科學版), 2018, 48(2): 113-119.

MEI Yuan, ZHANG Xu-Dong, SUN Li-Na, et al. Study on inversion of internal wave parameter using GF-1 and MODIS quasi-synchronous optical remote sensing images[J]. Periodical of Ocean University of China, 2018, 48(2): 113-119.

國家自然科學基金項目(61171161)資助

Supported by the National Natural Science Foundation of China(61171161)

2015-06-15；

2016-02-16

梅源(1990-)，男，碩士生，主要從事光學遙感方面的研究。E-mail：yuanmei.register@gmail.com

? ? 通訊作者：E-mail:wjing@ouc.edu.cn

責任編輯 陳呈超