王百順，范代讀，顧君暉，許建雄

（1.國家海洋局東海信息中心，上海 200137；2.同濟(jì)大學(xué) 海洋地質(zhì)國家重點(diǎn)實(shí)驗(yàn)室，上海 200092；3.同濟(jì)大學(xué) 長江水環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室，上海 200092）

海底地形是構(gòu)造運(yùn)動(dòng)、海洋動(dòng)力、泥沙輸移、生物活動(dòng)以及近岸顯著人類活動(dòng)等多種因素綜合作用的反映，其變化過程相當(dāng)復(fù)雜.海底地形地貌信息是海洋科學(xué)研究的基礎(chǔ)資料，對(duì)海洋漁業(yè)、海底資源開發(fā)、海洋權(quán)益保障和軍事活動(dòng)等具有重要意義.目前對(duì)海底地形的研究遠(yuǎn)遠(yuǎn)落后于陸地地形.

海洋測繪是以海圖等形式提供一個(gè)可以“看得見”的海底世界.水深數(shù)據(jù)可以通過單波束、多波束測深系統(tǒng)或海圖數(shù)字化等手段獲取，這些水深數(shù)據(jù)需要經(jīng)過一定的后期數(shù)據(jù)處理、網(wǎng)格化，生成海底數(shù)字高程模型，由此進(jìn)行海底地形的可視化與各種統(tǒng)計(jì)分 析［1－7］.內(nèi)插生成格網(wǎng)數(shù)字高層模型（digital elevation model，DEM）數(shù)據(jù)的方法較多，研究者往往根據(jù)研究區(qū)域大小、海底地形復(fù)雜性和水深數(shù)據(jù)獲取方式及其分辨率等的差異，采用不同的內(nèi)插值方法.任建武等［1］提出對(duì)比較離散的水深數(shù)據(jù)（單波束或海圖數(shù)字化獲?。┓中蝺?nèi)插.針對(duì)多波束測深數(shù)據(jù)的特點(diǎn)，高金耀等［2］、馬建林等［5］和王海棟等［7］分別使用了分形布朗運(yùn)動(dòng)（fractional Brown motion，F(xiàn)BM）和張力樣條結(jié)合的插值算法、改進(jìn)的距離反比權(quán)重法和多細(xì)節(jié)層次模型（level of detail，LOD）技術(shù)以及抗差最小二乘配置算法進(jìn)行海底數(shù)字地形模型的構(gòu)建，以達(dá)到相應(yīng)的精度和效率.

海底地形的可視模型是通過一系列的圖表進(jìn)行直觀表征.任建武等［1］通過水深分層設(shè)色和多方位分層設(shè)色法建立了南黃海輻射沙洲的立體可視、可比模型及各層水面積統(tǒng)計(jì)表等.楊剛等［4］提出利用最大梯度追蹤算法檢測正、負(fù)地形與平緩地形的分界線，并利用“三態(tài)值”算法識(shí)別正、負(fù)和平緩地形.隨著海洋測繪技術(shù)的發(fā)展，水深數(shù)據(jù)時(shí)／空分辨率和海域使用與開發(fā)程度逐步提高，可視模型也將更加多樣化.

2005—2010年期間，國家海洋局東海分局承擔(dān)多項(xiàng)“九O八”專項(xiàng)調(diào)查項(xiàng)目，共完成我國近海海域12個(gè)區(qū)塊的地形測量.本文選擇DX47區(qū)塊南部調(diào)查區(qū)，即對(duì)南黃海輻射沙洲北翼的單波束測深數(shù)據(jù)進(jìn)行海底地形的可視模型分析，包括：利用殘差分析以確定適合的插值方法，利用地形各參數(shù)進(jìn)行地形分割，對(duì)脊、槽相間的海底地形單元進(jìn)行精細(xì)劃分和定量分析.研究結(jié)果為2009年6月10日國務(wù)院討論通過的《江蘇沿海地區(qū)發(fā)展規(guī)劃》［8］的實(shí)施，尤其在沙脊實(shí)施高灘圈圍和沙脊間潮流主干通道建設(shè)萬噸級(jí)以上的深水航道提供依據(jù).

1 研究區(qū)概況

南黃海輻射沙洲分布在江蘇岸外，北起射陽河口，南至長江口北部的蒿枝港，南北長200km，東西寬140km，共有70多條沙脊及沙脊之間的潮流通道組成，面積達(dá)22470km2.規(guī)模巨大、形態(tài)獨(dú)特的輻射沙脊在全球?qū)崒俸币姡缫褳槲覈Ｑ蟮刭|(zhì)學(xué)家和海岸地貌學(xué)家所關(guān)注，自20世紀(jì)60年代以來，相繼開展了各種調(diào)查［9－12］.20世紀(jì)80年代任美鍔院士領(lǐng)導(dǎo)的“江蘇省海岸帶和海涂資源綜合調(diào)查”較清晰地揭示了輻射狀沙脊群的面貌［10］.1990—1995年由南京大學(xué)王穎院士主持的國家自然科學(xué)基金“八五”重點(diǎn)課題“黃海輻射沙洲形成演變研究”，聯(lián)合河海大學(xué)、同濟(jì)大學(xué)與中國科學(xué)院海洋研究所等，開展了輻射沙洲水動(dòng)力、沉積地貌與演變、晚第四紀(jì)地層與沉積環(huán)境演變等的綜合研究［11］.

據(jù)實(shí)測地形資料分析，沙洲在不斷增長，年均增長面積80km2.其中，沙洲外圍的水下沙脊遭受沖刷，侵蝕的沉積物較粗部分被輸送到沙脊頂部，使沙脊不斷增高和變寬，出露水面，近陸一側(cè)不斷發(fā)生沙洲淤高合并；較細(xì)的物質(zhì)則隨潮流懸浮輸送至潮間帶沉積，使沿岸灘涂不斷淤漲［11］.沙洲和潮灘不斷淤漲為江蘇省提供了極為寶貴的潛在土地資源，而潮流通道沖刷且穩(wěn)定為江蘇省平原海岸建設(shè)深水港提供了特殊、有利的條件，一些水深超過15m以上的大型潮流通道可建設(shè)10萬噸級(jí)以上深水航道［13］.這些有利條件一定程度上促使一省的發(fā)展規(guī)劃升格為國家戰(zhàn)略，隨著國務(wù)院于2009年6月10日討論并通過了《江蘇沿海地區(qū)發(fā)展規(guī)劃》，輻射沙洲的海底地形地貌及其演變趨勢和動(dòng)力機(jī)制再次成為研究熱點(diǎn)［8，13－14］.根據(jù)規(guī)劃，到2020年將對(duì)在鹽城射陽河口至南通東灶港之間的大豐、東臺(tái)、如東、啟東海岸和弶港外的輻射沙洲等地進(jìn)行圍墾，形成1800平方公里左右的墾區(qū)，遠(yuǎn)期將在東沙、條子泥、腰沙等淺灘沙洲進(jìn)行圍墾，可望形成4666.7km2土地后備資源.同時(shí)，擴(kuò)大港口能力，推進(jìn)大豐港區(qū)、洋口港區(qū)深水航道升級(jí)前期研究論證工作［8］.

2 水深數(shù)據(jù)來源與測量方法

測量資料來自國家海洋局東海信息中心.按照“九O八”專項(xiàng)任務(wù)設(shè)計(jì)，于2008年9月—2009年10月在南黃海輻射沙洲北翼實(shí)施了單波束水深測量.測量區(qū)域位于江蘇省大豐市、如東縣和射陽縣岸外近海海域（見圖1），距離海岸線6.14～20.83km，南北長134.02km，東西寬44.01～71.45km，面積6550km2，測點(diǎn)共5898個(gè).主測線間距2.5km，檢測線平均間距12.5km.

測量儀器設(shè)備主要有：差分GPS導(dǎo)航定位儀，美國 Trimble／Ashtech 公司生產(chǎn)，型號(hào)DSM／BR2G，實(shí)時(shí)定位精度≤3m；精密測深儀，無錫海鷹加科生產(chǎn)，型號(hào) HY1600，測深范圍0.3～300m，測深精度0.1%；自容式潮位儀，日本ALEC生產(chǎn)，型號(hào)ACTD－CMP，分辨率0.02m，精度±0.1%滿量程.

圖1 研究區(qū)位置及測量站點(diǎn)分布示意圖（站點(diǎn)抽?。▎挝唬簁m）Fig.1 Sketch map of the study area and distribution of bathymetric surveying stations（unit：km）

潮位校對(duì)采用實(shí)測和預(yù)報(bào)相結(jié)合方法.距0m等深線20km以內(nèi)的水深訂正采用實(shí)測驗(yàn)潮資料；在小洋口港設(shè)置有臨時(shí)驗(yàn)潮站，同時(shí)，使用了周邊嵐山、連云港、燕尾、濱海、太陽沙等長期站潮位數(shù)據(jù).岸線20km以外采用國家海洋信息中心提供的預(yù)報(bào)潮位（理論深度基準(zhǔn)面），根據(jù)實(shí)測和預(yù)報(bào)潮位將實(shí)測水深訂正到理論深度基準(zhǔn)面.以往測量作業(yè)中的潮位改正、近海水深測量一般采用模擬法，沿岸一般采用直線分帶和三角分帶法.目前，無論近海還是沿岸水深測量一般都采用時(shí)差法進(jìn)行潮位訂正，在某些難以布設(shè)驗(yàn)潮站的海區(qū)，利用潮汐調(diào)和常數(shù)和預(yù)報(bào)法進(jìn)行潮位改正，這也是DX47區(qū)塊采用的方法.實(shí)測潮位、定位和測深同步進(jìn)行，測量數(shù)據(jù)包括測點(diǎn)位置和相應(yīng)的水深.

3 最佳插值法篩選

地形數(shù)字高程分析及可視化需要對(duì)測點(diǎn)進(jìn)行插值，根據(jù)測量數(shù)據(jù)特征和分析目的選擇插值方法.表現(xiàn)原始數(shù)據(jù)空間分布的插值方法主要包括反距離加權(quán)法、克里金法、最小曲率法、改進(jìn)謝別德法、徑向基函數(shù)法和最近鄰點(diǎn)法等，各種方法還包含一系列的調(diào)整參數(shù).本次實(shí)驗(yàn)運(yùn)用殘差分析來選擇該批次水深數(shù)據(jù)的最佳插值方法.首先，用各種插值方法分別計(jì)算實(shí)際測量點(diǎn)上的插值；其次，計(jì)算各個(gè)站點(diǎn)測量值與插值的殘差，并對(duì)殘差進(jìn)行統(tǒng)計(jì)分析；最后根據(jù)統(tǒng)計(jì)參數(shù)對(duì)比，選擇最合適的插值方法.

用六種插值方法對(duì)5898個(gè)測量點(diǎn)（n）進(jìn)行殘差計(jì)算，統(tǒng)計(jì)參數(shù)見表1.除改進(jìn)謝別德法外，其余五種插值方法殘差均值普遍較小，說明實(shí)際測量值與插值總體一致.殘差最大值與最小值相對(duì)于研究區(qū)水深而言偏大，原因在于局部存在地形陡變，由于測量點(diǎn)密度不夠，在地形陡變區(qū)域進(jìn)行插值易生成較大殘差.可以采用多種參數(shù)對(duì)殘差進(jìn)行評(píng)價(jià)，特別是標(biāo)準(zhǔn)誤差、平均偏差和標(biāo)準(zhǔn)偏差.克里金法、最小曲率法和最近鄰點(diǎn)法的標(biāo)準(zhǔn)誤差最小，為0.017～0.018，最近鄰點(diǎn)法殘差均值較大；最小曲率法的平均偏差較克里金法大一些.因此，就研究海域現(xiàn)有測量數(shù)據(jù)而言，最合適的插值方法是克里金法，其次是最小曲率法，改進(jìn)謝別德法殘差最大，是最先應(yīng)該舍棄的插值方法.

表1 六種插值方法殘差分析統(tǒng)計(jì)（n＝5898）Tab.1 Residual analysis on the interpolation results of six different methods

4 海底地形可視化與統(tǒng)計(jì)分析

4.1 水深分層設(shè)色表達(dá)法

由水深值進(jìn)行分層設(shè)色是表達(dá)海底地形的最直觀簡潔方法，被廣泛應(yīng)用于航海、海洋工程和公益領(lǐng)域.研究區(qū)最大水深為36.7m，平均水深為6.0m.地形總體向東北呈加深趨勢.

沙脊或潮流槽在水深圖上呈條帶狀相間分布，向外呈扇形展開（見圖2）.區(qū)內(nèi)分布六條大型海底沙脊，長度35～115km，寬度3～11km.

根據(jù)水深變化，研究區(qū)可以定性地劃分出四種地形單元：淺灘、沙脊、潮流槽和寬谷.沙脊和潮流槽分布最廣，顯示為明顯的正、負(fù)地形；淺灘位于輻射沙洲的根部，低潮時(shí)可出露，地勢平緩（圖2中大片白色相連區(qū)域）；寬谷分布在潮流槽的外側(cè)，地勢低緩，水深常大于15m，逐漸過渡到開闊的陸架.用水深劃分的地形單元之間呈過渡狀態(tài).

圖2 水深設(shè)色海底地形圖Fig.2 Bathymetric map using color bar

4.2 地形傾角表達(dá)法

傾角ST是重要的地形參數(shù)之一.在地形分析中，一般在剖面上對(duì)傾角進(jìn)行分析，但剖面上的傾角不一定就是相應(yīng)點(diǎn)的最大傾角，因?yàn)槠拭娣较虿灰欢偸谴怪钡匦巫畲髢A斜方向，且剖面難以表示整個(gè)研究區(qū)的地形傾角.根據(jù)地形傾角公式（見公式（1））［15－16］，可以計(jì)算整體每個(gè)測量點(diǎn)或插值點(diǎn)與水平面之間的最大傾角，以度表示.

式中：x，y，z為大地坐標(biāo)系的x，y，z軸，分別對(duì)應(yīng)于經(jīng)度、緯度和高程，m.

計(jì)算結(jié)果表明，研究區(qū)最小傾角為0°，最大傾角31.8°，平均傾角5.3°.分布最廣的是傾角＜5°的地形區(qū)，傾角≥20°的地形區(qū)僅在局部分布（見圖3）.根據(jù)傾角大小可將研究區(qū)劃分成三個(gè)地形單元：0°～5°為水平－微傾斜地形區(qū)，包括地勢低緩的寬谷和淺灘，以及沙脊頂部與潮流槽的谷底，占總面積的60%；5°～10°為緩傾斜地形區(qū)，主要分布在沙脊的峰轉(zhuǎn)坡處或潮流槽的谷轉(zhuǎn)坡處，占26%；≥10°為傾斜地形區(qū)，分布在沙脊與潮流槽相接的斜坡處，占14%.

4.3 地形傾向

地形傾向AT也是基本的地形參數(shù)，有時(shí)伴隨著傾角一起描述.傾向是測點(diǎn)或插值點(diǎn)最大傾斜方向，采用公式（2）［15－16］計(jì)算每個(gè)插值點(diǎn)或測量點(diǎn)傾向，以度表示，向北傾向計(jì)為0°，順時(shí)針方向測算.地形傾向從0°～360°，主要傾向?yàn)?5°～90°，其次為90°～150°，260°～360°，傾向0°～45°和150°～260°地形區(qū)欠發(fā)育（見圖4）.

4.4 地形剖面曲率

在地形分析時(shí)，一般只研究三個(gè)基本參數(shù)：深（高）度、角度、傾向.可用多種方式表達(dá)深（高）度參數(shù)在空間上的分布，如平面等值線圖、影像圖、渲染圖、線框圖、三維表面圖等，角度參數(shù)多以剖面表達(dá)，傾向參數(shù)一般不做專門分析.很少有學(xué)者分析其他地形參數(shù)，如地形曲率等，至于更高級(jí)的地形分析則僅在理論教學(xué)探討中進(jìn)行.事實(shí)上，各種地形曲率和高級(jí)地形分析對(duì)于地形分區(qū)、地形重要界線確定、地形成因等研究具有重要意義.

地形剖面曲KP率表示沿最大梯度方向傾斜的變化率.當(dāng)?shù)匦纹拭媲蕿檎禃r(shí)，表示地形上凸，為負(fù)值時(shí)，表示地形下凹.脊是一種上凸地形單元，而槽是一種下凹地形單元，對(duì)于任一個(gè)測點(diǎn)或插值點(diǎn)，可以公式（3）［15－16］求其曲率.當(dāng)曲率為正則表明該點(diǎn)位于上凸曲面上，即位于沙脊之上.當(dāng)曲率為負(fù)值則表明該點(diǎn)位于下凹曲面之上，即位于潮流槽中.當(dāng)曲率為零則表明該點(diǎn)處于曲面凸凹銜接處，即脊、槽的分界線上.由此，可以明確界定出脊、槽地形單元的界線（見圖5）.

圖5 地形剖面曲率劃分的槽與脊Fig.5 Topographic divisions of ridges and channels based on the profile curvature

剖面曲率不僅可以將原水深地形圖上的根部細(xì)微地形單元?jiǎng)澐珠_來，而且將次級(jí)的脊、槽也劃分開來.由此發(fā)現(xiàn)，沙脊上普遍發(fā)育“次級(jí)槽”，而潮流槽中多有發(fā)育“次級(jí)脊”，形態(tài)各異，關(guān)系復(fù)雜，可見脊、槽體系水動(dòng)力環(huán)境的復(fù)雜多變.

根據(jù)剖面曲率，研究區(qū)地形單元可分為三大類：①沙脊，形態(tài)各異，總體呈條帶狀，面積2721km2；②潮流槽，分布于沙脊之間，向外有次級(jí)分汊現(xiàn)象，面積2930km2；③寬谷，面積899km2.以上三種地形單元分別占總面積的41.6%，44.7%和13.7%.與水深數(shù)據(jù)劃分的地形單元相比，缺少根部的淺灘.

4.5 地形平面曲率與切線曲率分析

等值線的彎曲程度可以用地形平面曲率KH來測量，用水平方向上地形方位角的變化率來表示，見公式（4）.地形切線曲率KT是測量垂直于梯度方向的垂直面走向的變化率，或等值線切線方向的變化率，見公式（5）［15－16］.

研究區(qū)地形平面曲率為－20×10－6～23475×10－6，平均891×10－6.圖6為地形平面曲率與地形水深疊加圖，可見地形比較復(fù)雜；正值區(qū)占49.5%，負(fù)值區(qū)占50.5%.地形平面曲率與水深、傾角為依據(jù)劃分的地形單元相關(guān)性低，這表明沙脊和潮流槽的地形復(fù)雜程度相似.研究區(qū)切線曲率與平面曲率分布特征相似，僅在數(shù)值和局部存在細(xì)微差異.這兩個(gè)參數(shù)反映了相似的海洋學(xué)涵義，即負(fù)值區(qū)位于輻40射沙洲的內(nèi)側(cè)，水深較淺，輻聚流場較強(qiáng)；正值區(qū)位于輻射洲的外側(cè)，水深較大，流場較弱.

圖6 地形平面曲率分布Fig.6 Distribution pattern of plan curvature

5 結(jié)論

南黃海輻射沙洲海底地形較復(fù)雜，采用多種地形參數(shù)可以從不同角度更好地表達(dá)地形特征與變化，并進(jìn)行定量分析.殘差分析表明，克里金法是研究區(qū)現(xiàn)有水深數(shù)據(jù)的最佳插值法.通過水深分層設(shè)色法可以定性地區(qū)分出輻射沙洲的四種地形單元：淺灘、沙脊、潮流槽和寬谷.不同地形單元或同一地形單元的不同部位可以有明顯不同的坡度和傾向，其中以水平－微傾斜地形（傾角＜5°）為主，占整個(gè)沙洲面積的60%，分布在寬谷、淺灘、沙脊頂部和潮流槽的谷底等.地形剖面曲率分層設(shè)色可以更細(xì)致地刻畫脊、槽地形單元的微地貌，包括沙脊上發(fā)育的次級(jí)槽和潮流槽中的次級(jí)脊等，并且可以定量分析不同地形單元面積及其所占比例，但該方法未能識(shí)別出淺灘地形單元.研究區(qū)切線曲率與平面曲率分布特征非常相似，兩者可能反映了相似的海洋動(dòng)力特征，即正（負(fù)）值表示輻射沙洲外（內(nèi)）側(cè)較弱（強(qiáng)）的流場.

［1］任建武，孫亞梅，李海宇.海底地形及其變化可視模型的研究［J］.測繪學(xué)報(bào)，1996，25（4）：257.REN Jianwu，SUN Yamei，LI Haiyu.A study on visual models of submarine topography and its change［J］.Acta Geodaetica et Cartographica Sinica，1996，25（4）：257.

［2］高金耀，金翔龍，吳自銀.多波束數(shù)據(jù)的海底數(shù)字地形模型構(gòu)建［J］.海洋通報(bào)，2003，22（1）：30.GAO Jinyao，JIN Xianglong， WU Ziyin.Construction of submarine DTM from Raw Multibeam Data［J］.Marine Science Bulletin，2003，22（1）：30.

［3］賈俊濤，翟京生，譚冀川.海底地形數(shù)據(jù)的格網(wǎng)DEM內(nèi)插方法［J］.海洋測繪，2003，23（6）：24.JIA Juntao，ZHAI Jingsheng，TAN Jichuan.The interpolation method of building grid to the data of terrain under sea［J］.Survering and Charting，2003，23（6）：24.

［4］楊剛，杜德文，呂海龍.數(shù)字海底地形分割算法［J］.海洋科學(xué)進(jìn)展，2004，22（2）：204.YANG Gang，DU Dewen，LU Hailong.Segmentation algorithm for digital sea－floor terrain［J］.Advances in Marine Science，2004，22（2）：204.

［5］馬建林，金菁，來向華.多波束測深海底數(shù)字地形模型的建立［J］.海洋測繪，2005，25（5）：15.MA Jianlin，JI Jing，LAI Xianghua.The construction of seafloor digital terrain model based on multibeam soundings［J］.Hydrographic Survering and Charting，2005，25（5）：15.

［6］周冠華，溫珍河，姜效典，等.南海海底地形可視化分析及其地質(zhì)意義［J］.海洋地質(zhì)與第四紀(jì)地質(zhì)，2006，26（2）：139.ZHOU Guanhau， WEN Zhenhe，JIANG Xiaodian，et al.Analysis on visualization of the submarine landorm of South China Sea and its geological meanings［J］.Marine Geology and Quaternary Geology，2006，26（2）：139.

［7］王海棟，柴洪洲，趙東明.基于抗差最小二乘配置的海底地形生成研究［J］.系統(tǒng)仿真學(xué)報(bào)，2010，22（9）：2091.WANG Haidong，CHAI Hongzhou，ZHAO Dongming.Research on seabed terrain generation based on robust least－squares collocation［J］.Journal of System Simulation，2010，22（9）：2091.

［8］國家發(fā)展改革委.江蘇沿海地區(qū)發(fā)展規(guī)劃［EB／OL］.（2010－01－26）［2010－12－26］http：／／www.sdpc.gov.cn／gzdt／P020100126607833080317.pdf，2009.National Development and Reform Commission of China.The Development Plan for the Jiangsu Coastal Zone ［EB／OL］.（2010－01－26）［2010－12－26］.http：／／www.sdpc.gov.cn／gzdt／P020100126607833080317.pdf，2009.

［9］李成治，李本川.蘇北沿海暗沙成因的研究［J］.海洋與湖沼，1981，12（4）：321.LI Chengzhi，LI Benchuan.Studies on the formation of Subei sand cays［J］.Oceanologia ET Limnologia Sinica，1981，12（4）：321.

［10］任美鍔.江蘇省海岸帶和海涂資源綜合調(diào)查報(bào)告［M］.北京：海洋出版社，1986.REN Mei’e.Comprehensive investingation of coastal zone and tidal flat resources，Jiangsu Province［M］.Beijing：Ocean Press，1986.

［11］王穎.黃海陸架輻射沙脊群［M］.北京：中國環(huán)境科學(xué)出版社，2002.WANG Ying.Radiative sandy ridge field on continnetal shelf of the Yellow Sea［M］.Beijing：China Environmental Science Press，2002.

［12］Li C，Zhang J，F(xiàn)an D，et al.Holocene regression and the tidal radial sand ridge system formation in the Jiangsu coastal zone，east China［J］.Marine Geology，2001，173（1－4）：97.

［13］朱大奎.江蘇海洋環(huán)境與沿海經(jīng)濟(jì)發(fā)展［J］.水資源保護(hù)，2003（6）：19.ZHU Dakui.Jiangsu marine environment and coastal economic development［J］.Water Resource Conservation，2003（6）：19.

［14］陳可鋒，陸培東，王艷紅，等.南黃海輻射沙洲趨勢性演變的動(dòng)力機(jī)制分析［J］.水科學(xué)進(jìn)展，2010，21（2）：267.CHEN Kefeng， LU Peidong， WANG Yanhong， et al.Hydrodynamic mechanism of evolvement trends in radial sandbank of south Yellow Sea，China［J］.Advances in Water Science，2010，21（2）：267.

［15］Moore I D， Lewis A， Gallant J C. Terrain attributes：estimation methods and scale effects［C］／／Jakeman A J，Beck M B，McAleer M J，eds.Modeling Change in Environmental Systems.Chichester：John Wiley and Sons，1993：189－214.

［16］Mitasova H，Hofierka J.Interpolation by regularized spline with tension：II.application to terrain modeling and surface geometry analysis［J］.Mathematical Geology，1993，25（6）：657.