• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于HJ衛(wèi)星的近岸Ⅱ類水體葉綠素a濃度定量遙感反演研究
    ——以灤河口北部海域?yàn)槔?/h1>
    2017-04-14 09:16:58羅建美霍永偉韓曉慶
    海洋學(xué)報(bào) 2017年4期
    關(guān)鍵詞:波段葉綠素反演

    羅建美,霍永偉,韓曉慶*

    (1. 河北地質(zhì)大學(xué) 土地資源與城鄉(xiāng)規(guī)劃學(xué)院,河北 石家莊 050031;2. 中國(guó)科學(xué)院遺傳與發(fā)育生物學(xué)研究所 農(nóng)業(yè)資源研究中心 中國(guó)科學(xué)院農(nóng)業(yè)水資源重點(diǎn)實(shí)驗(yàn)室,河北 石家莊 050022; 3. 中國(guó)科學(xué)院大學(xué),北京 100049;4. 河北省國(guó)土資源利用規(guī)劃院,河北 石家莊 050051)

    基于HJ衛(wèi)星的近岸Ⅱ類水體葉綠素a濃度定量遙感反演研究
    ——以灤河口北部海域?yàn)槔?/p>

    羅建美1,2,3,霍永偉4,韓曉慶4*

    (1. 河北地質(zhì)大學(xué) 土地資源與城鄉(xiāng)規(guī)劃學(xué)院,河北 石家莊 050031;2. 中國(guó)科學(xué)院遺傳與發(fā)育生物學(xué)研究所 農(nóng)業(yè)資源研究中心 中國(guó)科學(xué)院農(nóng)業(yè)水資源重點(diǎn)實(shí)驗(yàn)室,河北 石家莊 050022; 3. 中國(guó)科學(xué)院大學(xué),北京 100049;4. 河北省國(guó)土資源利用規(guī)劃院,河北 石家莊 050051)

    為驗(yàn)證環(huán)境衛(wèi)星影像在近岸Ⅱ類水體中葉綠素a濃度反演的適宜性,獲得適合灤河口北部近岸海域的葉綠素a濃度的高精度反演算法,基于HJ-1A CCD2影像數(shù)據(jù)和現(xiàn)場(chǎng)同步實(shí)測(cè)數(shù)據(jù),構(gòu)建、驗(yàn)證并確立適合研究區(qū)葉綠素a濃度反演模型。結(jié)果表明:HJ-1A CCD2衛(wèi)星可用于水文情況較復(fù)雜的近岸Ⅱ類水體葉綠素a濃度的反演研究;最佳波段組合(B3/B1)構(gòu)建的倒數(shù)模型,可反演葉綠素a濃度。結(jié)果顯示葉綠素a濃度整體由近岸向內(nèi)海方向逐漸減小,自北向南逐漸降低,并在河口處形成向外海突出的高值區(qū);葉綠素a濃度最大值出現(xiàn)在金山嘴南側(cè),為25.0 μg/L。離岸葉綠素a濃度最小值出現(xiàn)在灤河口以東約13.2 km處,僅4.6 μg/L;近岸葉綠素a濃度最小值出現(xiàn)在灤河口附近,約為11.0 μg/L。該研究可為相關(guān)研究提供技術(shù)參考,研究成果可為當(dāng)?shù)卣茖W(xué)管理海洋環(huán)境、制定海洋政策提供決策依據(jù),為該海域可持續(xù)發(fā)展提供數(shù)據(jù)支持。

    HJ 衛(wèi)星;灤河口北部;葉綠素a濃度;模型;反演

    1 引言

    葉綠素a是浮游植物、藻類植物中最豐富的色素,是浮游生物體的重要組成成分之一,其濃度的高低與該水體中藻類的種類和數(shù)量密切相關(guān),是表征水體富營(yíng)養(yǎng)化狀況的主要指標(biāo),也是水體富營(yíng)養(yǎng)化評(píng)價(jià)的重要參數(shù),更是研究全球碳循環(huán)的主要手段之一。因此國(guó)內(nèi)外學(xué)者對(duì)其進(jìn)行了深入研究,目前葉綠素a濃度的計(jì)算方法主要有經(jīng)驗(yàn)統(tǒng)計(jì)方法[1—8]、半分析半經(jīng)驗(yàn)方法[9—13]和解析模型方法[14—20]3類。上述3種方法各有利弊,與后兩種方法相比,經(jīng)驗(yàn)統(tǒng)計(jì)方法針對(duì)某一特定區(qū)域具有較好的適用性,可較好的反演特定區(qū)域的葉綠素a濃度,便于建立業(yè)務(wù)化應(yīng)用的葉綠素a濃度提取模型,在實(shí)際應(yīng)用中具有較大優(yōu)勢(shì)。

    與常規(guī)監(jiān)測(cè)方法相比,基于遙感的葉綠素a濃度監(jiān)測(cè)具有大尺度、實(shí)時(shí)、動(dòng)態(tài)監(jiān)測(cè)的優(yōu)勢(shì),國(guó)內(nèi)外學(xué)者利用CZCS、OCTS、OCM、SeaWiFS、MODIS、MERIS、Landsat TM、Hyperion等國(guó)外多光譜、高光譜衛(wèi)星影像進(jìn)行了大量研究,成果較為突出[21—41]。應(yīng)用國(guó)產(chǎn)衛(wèi)星進(jìn)行葉綠素a濃度的研究起步較晚,但近幾年隨著我國(guó)衛(wèi)星迅速發(fā)展,傳感器觀測(cè)波段及觀測(cè)精度較以往有了很大幅度提高,國(guó)內(nèi)學(xué)者利用國(guó)產(chǎn)衛(wèi)星,如HY-1A/B、CBERS、GF-1、HJ-1A/B CCD、HJ-1A HIS等,針對(duì)葉綠素a的反演進(jìn)行了大量研究,取得了諸多富有成效的研究成果[42—52]。其中,HJ系列衛(wèi)星多光譜數(shù)據(jù)因其分辨率較高(30 m)、涵蓋常用波段(可見(jiàn)光及近紅外等4個(gè)波段)、對(duì)地刈幅寬(360 km)、重訪周期短(4 d)、可免費(fèi)獲取等特點(diǎn),在葉綠素a濃度反演中具有較強(qiáng)的優(yōu)勢(shì)。但在應(yīng)用范圍方面,HJ衛(wèi)星應(yīng)用于內(nèi)陸封閉水體較多[41,52—60],應(yīng)用于開(kāi)放海域的Ⅰ類水體較少,應(yīng)用于水質(zhì)條件復(fù)雜的近岸Ⅱ類水體則較少[61—64]。這是由于近岸Ⅱ類水體既受海水?dāng)y帶浮游植物、懸浮物、溶解有機(jī)物等影響,又受入海河流攜帶陸源懸浮泥沙、黃色物質(zhì)等影響,各組分之間彼此混合、交互作用。利用遙感影像反演Ⅱ類水體葉綠素濃度,受海水中上述物質(zhì)影響,使得適用于大洋Ⅰ類水體的經(jīng)驗(yàn)公式基本不適用于近岸Ⅱ類水體。即使同為Ⅱ類水體,因水質(zhì)組分及光譜特征差別較大,其反演算法也隨時(shí)空不同存在較大區(qū)別,建立更大尺度的普適性反演算法具有一定困難,但針對(duì)較小尺度建立的反演模型,在特定區(qū)域其反演精度卻比較精準(zhǔn),尤其是結(jié)合衛(wèi)星過(guò)境時(shí)的水質(zhì)同步實(shí)測(cè)數(shù)據(jù),其反演結(jié)果的準(zhǔn)確性和可信度則更高[16, 51, 58, 61]。

    本研究以灤河口北部近岸Ⅱ類海域?yàn)檠芯繀^(qū),選取2014年HJ衛(wèi)星遙感影像的適當(dāng)波段組合,基于SPSS軟件構(gòu)建葉綠素a濃度估算模型,結(jié)合現(xiàn)場(chǎng)同步實(shí)測(cè)數(shù)據(jù),評(píng)價(jià)模型的準(zhǔn)確性和適宜性,以期獲得適合研究區(qū)葉綠素a濃度的高精度反演算法,為國(guó)產(chǎn)衛(wèi)星在近岸海域水質(zhì)遙感動(dòng)態(tài)監(jiān)測(cè)和預(yù)測(cè)等方面的應(yīng)用提供案例支持。研究方法和過(guò)程可為同類研究提供技術(shù)參考,研究成果可為當(dāng)?shù)卣茖W(xué)管理海洋環(huán)境、制定海洋政策提供決策依據(jù),為該海域可持續(xù)發(fā)展提供數(shù)據(jù)支持。

    2 研究區(qū)域

    研究區(qū)位于河北省秦皇島市昌黎縣,灤河口以北、金山嘴以南近海海域,介于39.4°~39.8°N、119.3°~119.7°E之間,面積約800 km2,平均水深7 m。屬半濕潤(rùn)大陸性氣候,年均溫11℃,年均降水量712 mm,入海水系主要有灤河、七里海、大蒲河及北側(cè)的洋河、戴河等。大量懸浮泥沙、溶解物、生活污水及沿海養(yǎng)殖池塘排水隨徑流入海,成分復(fù)雜,其中,水體中溶解氧(DO)6.3~8.6 mg/L,pH(酸堿值)8.0~8.1,鹽度28.9~31.1,油類52.7~127.0 μg/L,水溫15.0~20.3℃,透明度0.8~1.8 m,化學(xué)需氧量(COD)1.0~1.6 mg/L,水質(zhì)指標(biāo)符合國(guó)家Ⅱ類水質(zhì)標(biāo)準(zhǔn)。本研究實(shí)測(cè)采樣點(diǎn)共24個(gè),呈4列南北向分布(圖1)。

    3 研究數(shù)據(jù)及處理

    本研究主要采用HJ衛(wèi)星遙感影像及同步實(shí)測(cè)數(shù)據(jù),基于ENVI、ArcGIS、SPSS等軟件平臺(tái),將遙感影像進(jìn)行幾何校正、大氣校正等預(yù)處理,通過(guò)分析不同波段與實(shí)測(cè)數(shù)據(jù)之間的相關(guān)性,選取適合研究區(qū)的波段組合構(gòu)建數(shù)據(jù)模型,并檢驗(yàn)其準(zhǔn)確性,最后利用數(shù)據(jù)模型計(jì)算研究區(qū)葉綠素a濃度情況。研究技術(shù)路線見(jiàn)圖2。

    3.1 遙感影像數(shù)據(jù)及處理

    3.1.1 影像數(shù)據(jù)

    本研究采用HJ-1A CCD2的影像數(shù)據(jù),景號(hào)為453-68,成像時(shí)間為2014年5月22日01時(shí)41分03秒(格林威治時(shí)間),分辨率是30 m,研究區(qū)無(wú)云,數(shù)據(jù)來(lái)源于中國(guó)資源衛(wèi)星應(yīng)用中心網(wǎng)站。

    參考影像為兩景經(jīng)過(guò)幾何精校正的Landsat5-TM影像,景號(hào)分別為121-032、121-033,成像時(shí)間分別是2010年9月27日、2010年9月11日,分辨率均為30 m,研究區(qū)無(wú)云。為便于后期處理,已將這兩景影像進(jìn)行鑲嵌處理。參考影像為西安80 坐標(biāo)系,6 度分帶20 N帶,Gauss Kruger投影,IAG75 橢球體。

    3.1.2 影像處理

    (1)輻射定標(biāo)

    為實(shí)現(xiàn)定量研究,需建立遙感傳感器的數(shù)字量化輸出值DN與其所對(duì)應(yīng)視場(chǎng)中輻射亮度值之間的定量關(guān)系,即對(duì)影像進(jìn)行輻射定標(biāo)。

    L=DN/g+L0,

    (1)

    式中,L為輻射亮度值,單位是W/(m2·sr·μm),DN為影像的數(shù)字量化值,g為絕對(duì)定標(biāo)系數(shù)增益,單位是W-1·m2·sr·μm,L0為偏移值,單位是W/(m·sr·μm)。

    在ENVI 5.1中,利用針對(duì)環(huán)境衛(wèi)星影像使用IDL語(yǔ)言開(kāi)發(fā)的Envi_hj1a1b_tools工具,通過(guò)讀取頭文件信息,進(jìn)行輻射定標(biāo)、圖層疊加等操作。

    (2)影像裁剪

    由于環(huán)境星的成像幅面為360 km,遠(yuǎn)超出研究所需范圍,為減小處理數(shù)據(jù)量,利用Subset Data from ROI模塊裁剪出研究區(qū)范圍,以便后期數(shù)據(jù)處理。在裁剪時(shí)保留一部分陸地,利于幾何校正尋找同名地物點(diǎn)操作。

    圖1 研究區(qū)范圍示意圖和采樣點(diǎn)位置Fig.1 The study area and the sampling positions

    圖2 研究技術(shù)路線圖Fig.2 The technology roadmap of the study

    (3)幾何校正

    利用ENVI的流程化模塊進(jìn)行幾何校正,該模塊通過(guò)設(shè)置自動(dòng)控制點(diǎn)的匹配方法、過(guò)濾條件、校正波段、控制點(diǎn)個(gè)數(shù)等條件,通過(guò)選取具有相同波段信息的影像圖層,軟件可以自動(dòng)尋找同名地物點(diǎn),最后只需刪除或修改自動(dòng)生成的控制點(diǎn)即可。本研究設(shè)置匹配方法為Cross Correlation,最小匹配值0.6,匹配模型為Fitting Global Transform,最大允許誤差數(shù)為1,搜索窗口為200,匹配窗口為61,匹配波段均為紅色波段,重采樣方法為最近鄰法,以最大程度減少光譜信息損失,最終校正結(jié)果的總誤差控制在0.3個(gè)像元內(nèi)。

    (4)大氣校正

    常用的大氣校正模型有暗黑像元模型[41,65]、6S模型[19,28,33,59]、FLAASH模型[13,48,52,54,56,58,62-63,66]等,根據(jù)韓曉慶等[67]的研究,在海岸帶地區(qū)FLAASH模型的校正精度要優(yōu)于其他兩種模型,杜挺等[53]的研究也證實(shí)了這一結(jié)果,因此本研究選用FLAASH模型進(jìn)行大氣校正,各輸入?yún)?shù)見(jiàn)表1。

    表1 FLAASH大氣校正模型輸入?yún)?shù)表

    (5)海陸分離

    采用ENVI的Feature Extraction(面向?qū)ο蠓诸?,?jiǎn)稱FE)模塊進(jìn)行海陸分離操作,自動(dòng)提取研究區(qū)海域[68]。該模塊通過(guò)集合臨近像元識(shí)別感興趣的光譜要素,利用遙感影像數(shù)據(jù)的空間、紋理和光譜信息進(jìn)行分割和分類的特點(diǎn),設(shè)置分割及合并閾值生成一個(gè)分割圖像,并輸出矢量等多種格式。與手工提取、監(jiān)督分類或NDWI等模型提取海域相比,該方法不但可以有效提高信息提取效率,避免人機(jī)交互提取時(shí)產(chǎn)生的人為誤差、監(jiān)督分類因樣本不典型造成分類結(jié)果的偏差及NDWI等模型因地物復(fù)雜導(dǎo)致提取信息不準(zhǔn)確等問(wèn)題,而且可以隨時(shí)預(yù)覽處理效果,如果不符合要求,只需簡(jiǎn)單調(diào)整閾值重新預(yù)覽,直至得到滿意的結(jié)果,節(jié)省了數(shù)據(jù)處理時(shí)間和數(shù)據(jù)存儲(chǔ)空間。該方法尤其適用于地物影像特征變化不大,地物邊界相對(duì)明顯的區(qū)域。本研究設(shè)置的分割閾值為45,合并閾值為95。

    3.2 同步實(shí)測(cè)數(shù)據(jù)及處理

    3.2.1 同步實(shí)測(cè)數(shù)據(jù)

    研究區(qū)為海域,與內(nèi)陸水體相比,海水的周期性運(yùn)動(dòng)使得海水中的物質(zhì)交換更加頻繁,隨海水運(yùn)移的速度更快,海域反演要求影像成像時(shí)間與數(shù)據(jù)采集時(shí)間間隔更短。為使研究結(jié)果真實(shí)可信,減少延時(shí)誤差,本研究現(xiàn)場(chǎng)數(shù)據(jù)采集時(shí)間為2014 年5 月22 日9 時(shí)至12 時(shí),與環(huán)境衛(wèi)星過(guò)境時(shí)間基本保持一致。數(shù)據(jù)采集當(dāng)天晴朗無(wú)云,風(fēng)速2.1 m/s,風(fēng)向ENE,波浪較小。

    3.2.2 數(shù)據(jù)處理

    同步采樣點(diǎn)共24 個(gè),使用5 L的卡蓋式采水器采集表層海水水樣,深度為0.5 m,每個(gè)采樣點(diǎn)取3 個(gè)樣品,在實(shí)驗(yàn)室中分別測(cè)定3個(gè)樣品后,取其均值為該采樣點(diǎn)的葉綠素濃度,并用GPS同步測(cè)量樣點(diǎn)經(jīng)緯度。

    在實(shí)驗(yàn)室采用分光光度法測(cè)定海水中的葉綠素a濃度。使用玻璃纖維濾膜(Whatman GF/C,47 mm)過(guò)濾海水樣品,然后放入貝賽離心管,加入濃度為90% 的丙酮溶液,搖勻后放置冰箱儲(chǔ)存24 h,溫度為4℃。此后,將該樣品放入離心機(jī),速度設(shè)定為4 000 r/min,離心時(shí)間為10 min。以濃度為90% 的丙酮溶液作參比,使用分光光度計(jì)分別在750 nm、664 nm、647 nm、630 nm波長(zhǎng)處測(cè)定吸光值。按公式(2)計(jì)算得出樣品中葉綠素a濃度:

    ρChla=[11.85×(E664-E750)-1.54×

    (2)

    式中,ρChla為樣品中葉綠素a濃度,單位為μg/L;E664、E750、E647、E630分別為664nm、750nm、647nm、630nm波長(zhǎng)處測(cè)定吸光值;ν為樣品中提取液體積,單位為mL;V為海水樣品實(shí)際用量,單位為L(zhǎng);L為測(cè)定池光程,單位為cm。

    研究樣品測(cè)定葉綠素a濃度范圍6.9~16.3 μg/L。

    3.3 葉綠素反演

    利用準(zhǔn)同步HJ-1A CCD2和實(shí)測(cè)葉綠素a濃度數(shù)據(jù),分析得出與實(shí)測(cè)數(shù)據(jù)相關(guān)性最大的波段組合,建立葉綠素a濃度反演模型,并對(duì)模型精度進(jìn)行對(duì)比分析。根據(jù)前人經(jīng)驗(yàn)[51],將野外實(shí)測(cè)24 個(gè)樣點(diǎn)分為兩組,其中2/3(16 個(gè))樣點(diǎn)用于模型反演,1/3(8 個(gè))樣點(diǎn)用于驗(yàn)證模型。

    3.3.1 相關(guān)分析

    利用SPSS 20計(jì)算各波段與實(shí)測(cè)葉綠素a濃度的相關(guān)系數(shù),結(jié)果表明B3(紅色)波段與葉綠素a濃度的相關(guān)性最高,達(dá)到0.568,置信度為0.043,表明為在置信度(雙側(cè))為0.05時(shí),二者顯著相關(guān)。B2(綠色)波段,相關(guān)系數(shù)為0.437,B1(藍(lán)色)波段為0.347,B4(近紅外)波段為0.311(表2)。

    表2 各波段反射率與葉綠素a濃度實(shí)測(cè)值的相關(guān)關(guān)系

    注:*為在置信度(雙側(cè))為0.05時(shí),相關(guān)性是顯著的。

    部分學(xué)者研究成果表明,波段比值在一定程度上能進(jìn)一步消除大氣對(duì)電磁波傳輸信號(hào)的影響[69],因此,本研究進(jìn)一步嘗試采用不同波段比值與實(shí)測(cè)葉綠素a濃度進(jìn)行相關(guān)分析,以便得出與研究區(qū)實(shí)測(cè)數(shù)據(jù)最吻合的波段組合。

    在ENVI中,利用Band Math功能分別對(duì)影像的4 個(gè)波段進(jìn)行比值運(yùn)算,計(jì)算結(jié)果與實(shí)測(cè)葉綠素a濃度值進(jìn)行相關(guān)分析,結(jié)果顯示B3/B1波段比值與實(shí)測(cè)葉綠素a濃度相關(guān)性最高,為0.720,置信度為0.006,表明在置信度(雙側(cè))為0.01時(shí),二者顯著相關(guān)(表3),這雖與陳莉瓊等[16]利用HJ衛(wèi)星反演葉綠素時(shí)的結(jié)論相似,但與其他學(xué)者常用的藍(lán)、綠波段比值及紅、近紅外波段組合反演水體葉綠素的算法有所不同[56,61—62]。為此,本研究進(jìn)一步分析反射率影像的光譜曲線特征,發(fā)現(xiàn)除個(gè)別點(diǎn)外,大部分反演點(diǎn)B1(藍(lán)波段)的反射率值最高,B2(綠波段)、B3(紅波段)、B4(近紅外波段)的反射率值依次降低(圖3),并非如一般海水中,B2(綠波段)反射率最高,B1(藍(lán)波段)、B2(綠波段)、B4(近紅外波段)反射率值較低。這可能由于研究區(qū)水體受石油污染(52.7~127.0 μg/L,平均85.1 μg/L),水體的影像特征發(fā)生變化,導(dǎo)致研究區(qū)葉綠素反演波段與常規(guī)反演波段不同,但該波段組合的反演結(jié)果較之其他組合更能準(zhǔn)確反應(yīng)研究區(qū)葉綠素特征。

    表3 各波段組合反射率比值與葉綠素a濃度實(shí)測(cè)值的相關(guān)關(guān)系

    續(xù)表3

    注:**為在置信度(雙側(cè))為0.01時(shí),相關(guān)性是顯著的。

    圖3 反演點(diǎn)反射率光譜特征Fig.3 The spectral characteristics of reflectivity in inverse points

    3.3.2 數(shù)據(jù)建模及精度評(píng)價(jià)

    在SPSS中,以B3/B1為自變量,以實(shí)測(cè)葉綠素a濃度為因變量,選擇線性模型等10種常用的曲線模型,構(gòu)建波段與實(shí)測(cè)葉綠素a濃度之間的函數(shù)關(guān)系式(表4)。

    表4 葉綠素a濃度反演模型及估值對(duì)比

    續(xù)表4

    為選取最佳模型,計(jì)算了各模型的相關(guān)系數(shù)(r2)、均方根誤差(RMSE)及相對(duì)誤差(RE),對(duì)比結(jié)果下:二次曲線模型和三次曲線模型的Sig.值均大于0.01,表明上述兩個(gè)模型的模擬結(jié)果在置信度(雙側(cè))為0.01時(shí),與實(shí)測(cè)葉綠素a濃度不顯著相關(guān),因此不采用上述兩個(gè)模型;其余8個(gè)模型的Sig.值均小于0.01,表明在置信度(雙側(cè))為0.01時(shí),模型計(jì)算結(jié)果與實(shí)測(cè)葉綠素a濃度顯著相關(guān)。其中,線性模型、復(fù)合模型、成長(zhǎng)模型、指數(shù)模型的Sig.值與其余4個(gè)模型相比較大,且r2較小,表明相關(guān)度比其他4個(gè)模型低,故不采用上述4個(gè)模型;在剩余的4個(gè)模型中,冪指數(shù)模型與S曲線模型的RMSE值分別為2.877 μg/L、2.853 μg/L,RE比例分別為18.536%、18.517%,與對(duì)數(shù)模型、倒數(shù)模型相比,誤差較大,不采用上述兩種模型;與對(duì)數(shù)模型相比,倒數(shù)模型的模擬值與葉綠素a濃度實(shí)測(cè)值的相關(guān)系數(shù)更大,RMSE較小,且通過(guò)進(jìn)一步對(duì)比驗(yàn)證點(diǎn)的模擬值與實(shí)測(cè)值,結(jié)果顯示倒數(shù)模型的相關(guān)系數(shù)為0.939 5,對(duì)數(shù)模型為0.884 5(圖4),表明在10種曲線模型中,倒數(shù)模型能夠較真實(shí)的模擬研究區(qū)葉綠素a濃度實(shí)際情況。因此,選擇倒數(shù)模型作為研究區(qū)海域葉綠素a濃度反演模型。

    圖4 倒數(shù)(a)、對(duì)數(shù)(b)模型驗(yàn)證點(diǎn)實(shí)測(cè)值與模擬值比較Fig.4 Comparison in values between measured and estimated by reciprocal(a) and logarithmic(b) models

    圖5 葉綠素a濃度反演圖Fig.5 The inversion map of Chl a concentration

    4 研究結(jié)果

    在ENVI中利用Band Math功能,根據(jù)構(gòu)建的反演模型,計(jì)算得出葉綠素a濃度分布情況(圖5)。

    研究結(jié)果呈現(xiàn)以下特點(diǎn):葉綠素a濃度由岸向海逐漸減小,自北向南逐漸降低;戴河口的葉綠素a濃度等值線較密,且與洋河共同形成突入海中的葉綠素a濃度高值區(qū);大蒲河口葉綠素a濃度呈長(zhǎng)條狀突入海中;七里海入??谔卣鞑幻黠@;灤河口形成錐形葉綠素高值區(qū)。戴河、洋河、大蒲河由于入海水量較大,突入海中的葉綠素a濃度呈盾型,該特征以大蒲河口最為明顯,呈長(zhǎng)條狀突入海中,這是由于河流沖淡水帶來(lái)上游大量的營(yíng)養(yǎng)鹽,在河口外形成突出的高營(yíng)養(yǎng)鹽區(qū)域,便于藻類生長(zhǎng);灤河入海水量較少,突入海中的葉綠素a濃度呈錐型;七里海內(nèi)水體隨潮水漲落而進(jìn)出河口,水體與中葉綠素a濃度與外側(cè)海水差別不大,河口處葉綠素a濃度特征不明顯。

    葉綠素a濃度最大值出現(xiàn)在金山嘴南側(cè),高達(dá)25.0 μg/L,因?yàn)樵搮^(qū)域自北向南的沿岸流,在金山嘴以南形成波影區(qū),波速較緩,海水中的葉綠素不易擴(kuò)散,形成富集,加之戴河的沖淡水帶來(lái)部分營(yíng)養(yǎng)成分,利于藻類穩(wěn)定生長(zhǎng);離岸葉綠素a濃度最小值出現(xiàn)在灤河口以東約13.2 km處,僅4.6 μg/L。近岸葉綠素a濃度最小值出現(xiàn)在南側(cè)灤河口附近,約11.0 μg/L左右。由于自21世紀(jì)初以來(lái),灤河入海水量逐漸減少,入海營(yíng)養(yǎng)鹽含量降低,葉綠素a濃度隨之減小。自北向南的沿岸流在此處達(dá)灣頂,流速較大,大量葉綠素a被沖刷到下游地區(qū)——灤河口西南側(cè)富集,也部分造成了灤河口北側(cè)葉綠素a較低;此外,研究區(qū)東側(cè)海域有一條寬度約3 km、葉綠素a濃度明顯高于左右兩側(cè)海域的縱向條帶,這是由于HJ-1A CCD2衛(wèi)星影像存在部分條帶和噪聲,在背景均一的地物如水體等位置,條帶現(xiàn)象更加明顯,導(dǎo)致反演結(jié)果出現(xiàn)部分異常。

    5 結(jié)論

    本文利用HJ-1A CCD2衛(wèi)星的30 m空間分辨率多光譜影像,結(jié)合現(xiàn)場(chǎng)同步實(shí)測(cè)數(shù)據(jù),在灤河口北部海域開(kāi)展葉綠素a濃度反演研究,主要結(jié)論如下:

    (1)通過(guò)與實(shí)測(cè)葉綠素a濃度進(jìn)行相關(guān)分析,確定(B3/B1)為反演研究區(qū)最佳波段組合,并根據(jù)該波段組合構(gòu)建倒數(shù)模型,可用于反演葉綠素a濃度。

    (2)葉綠素a濃度整體呈由岸向海逐漸減小、自北向南逐漸降低的特征,且在河口處形成向外海突出的高值區(qū)。

    (3)葉綠素a濃度最大值出現(xiàn)在金山嘴南側(cè),為25.0 μg/L。離岸葉綠素a濃度最小值出現(xiàn)在灤河口以東約13.2 km處,僅4.6 μg/L;近岸葉綠素a濃度最小值出現(xiàn)在灤河口附近,約11.0 μg/L左右。

    (4)研究結(jié)果證明HJ-1A CCD2衛(wèi)星可用于反演水文情況較復(fù)雜的沿海區(qū)域葉綠素a濃度反演研究,但影像本身的條帶噪聲導(dǎo)致反演結(jié)果出現(xiàn)部分異常。

    (5)本研究雖取得一定成果,但由于Ⅱ類水體光學(xué)特性的復(fù)雜性,不僅受浮游植物影響,還受其他物質(zhì)如懸浮物、黃色物質(zhì)等有色溶解有機(jī)物的影響。在下一步工作中,還要充分考慮上述物質(zhì)對(duì)葉綠素濃度反演的影響,進(jìn)一步分析波段組合,并修正模型參數(shù),以提高模型的普適性和反演精度。

    致謝:本研究衛(wèi)星影像由中國(guó)資源衛(wèi)星應(yīng)用中心提供,同步數(shù)據(jù)由國(guó)家海洋環(huán)境監(jiān)測(cè)中心同志指導(dǎo)測(cè)定,各位審稿專家為本文提出諸多建設(shè)性的意見(jiàn)和建議,在此一并致謝!

    [1] Moses W J, Gitelson A A, Berdnikov S, et al. Estimation of chlorophyll-aconcentration in case Ⅱ waters using MODIS and MERIS data-successes and challenges[J]. Environmental Research Letters, 2009, 4(4): 549-567.

    [2] Ekstrand S. Landsat TM based quantification of chlorophyll-aduring algae blooms in coastal waters[J]. International Journal of Remote Sensing, 1992, 13(10): 1913-1926.

    [3] Gitelson A A, Schalles J F, Hladik C M. Remote chlorophyll-aretrieval in turbid, productive estuaries: Chesapeake Bay case study[J]. Remote Sensing of Environment, 2007, 109(4): 464-472.

    [4] 馬萬(wàn)棟, 吳傳慶, 殷守敬, 等. 基于實(shí)測(cè)光譜的煙臺(tái)四十里灣水體葉綠素a遙感反演[J]. 中國(guó)環(huán)境監(jiān)測(cè), 2015, 31(4): 135-140.

    Ma Wandong, Wu Chuanqing, Yin Shoujing, et al. Chlorophyll-aretrieval of Yantai Sishili Bay based oninsituspectral characteristics[J]. Environmental Monitoring in China, 2015, 31(4): 135-140.

    [5] 蔣剛, 肖建, 鄭永康, 等. 基于統(tǒng)計(jì)學(xué)習(xí)理論的海色遙感葉綠素濃度反演方法[J]. 中國(guó)海洋大學(xué)學(xué)報(bào)(自然科學(xué)版), 2006, 36(4): 559-563.

    Jiang Gang, Xiao Jian, Zheng Yongkang, et al. Use ocean color remote sensing data to retrieve chlorophyll-aconcentration based on statistic learning theory[J]. Periodical of Ocean University of China, 2006, 36(4): 559-563.

    [6] 劉朝相. 內(nèi)陸水體葉綠素a濃度的遙感反演模型研究——以官?gòu)d水庫(kù)媯水河段為例[D]. 北京: 首都師范大學(xué), 2014.

    Liu Chaoxiang. Remote sensing retrieval model of chlorophyllaconcentration in inland waters[D]. Beijing: Capital Normal University, 2014.

    [7] 張勇, 張永豐, 張萬(wàn)磊, 等. 秦皇島海域微微型藻華期間葉綠素a分級(jí)研究[J]. 生態(tài)科學(xué), 2012, 31(4): 357-363.

    Zhang Yong, Zhang Yongfeng, Zhang Wanlei, et al. Size fraction of chlorophylladuring and after brown tide in Qinhuangdao coastal waters[J]. Ecological Science, 2012, 31(4): 357-363.

    [8] 熊竹, 丁世敏, 肖紅艷, 等. 水體中葉綠素a濃度遙感反演的算法研究進(jìn)展[J]. 中國(guó)西部科技, 2015, 14(10): 16-20.

    Xiong Zhu, Ding Shimin, Xiao Hongyan, et al. Progress of algorithms research on remote sensing retrieval of chlorophyll-aconcentration in water body[J]. Science and Technology of West China, 2015, 14(10): 16-20.

    [9] Magnuson A, Harding Jr L W, Mallonee M E, et al. Bio-optical model for Chesapeake Bay and the Middle Atlantic Bight[J]. Estuarine, Coastal and Shelf Science, 2004, 61(3): 403-424.

    [10] D’Sa E J, Hu C, Muller-Karger F E, et al. Estimation of colored dissolved organic matter and salinity fields in case 2 waters using SeaWiFS: Examples from Florida Bay and Florida Shelf[J]. Journal of Earth System Science, 2002, 111(3): 197-207.

    [11] 矯京均, 辛?xí)灾? 余珊珊, 等. HJ-1衛(wèi)星數(shù)據(jù)估算地表能量平衡[J]. 遙感學(xué)報(bào), 2014, 18(5): 1048-1058.

    Jiao Jingjun, Xin Xiaozhou, Yu Shanshan, et al. Estimation of surface energy balance from HJ-1 satellite data[J]. Journal of Remote Sensing, 2014, 18(5): 1048-1058.

    [12] 黃薇. 膠州灣葉綠素a濃度遙感產(chǎn)品檢驗(yàn)與反演算法研究[D]. 青島: 中國(guó)海洋大學(xué), 2012.

    Huang Wei. MERIS chlorophyll-aconcentration product validation and researching on the inversion algorithm of chlorophyll-aconcentration for the Jiaozhou Bay[D]. Qingdao: Ocean University of China, 2012.

    [13] 楊錦坤, 陳楚群. 珠江口二類水體水色三要素的優(yōu)化反演[J]. 熱帶海洋學(xué)報(bào), 2007, 26(5): 15-20.

    Yang Jingkun, Chen Chuqun. An optimal algorithm for retrieval of chlorophyll, suspended sediments and gelbstoff of case Ⅱ waters in Zhujiang River estuary[J]. Journal of Tropical Oceanography, 2007, 26(5): 15-20.

    [14] Keiner L E, Yan Xiaoyan. A Neural network model for estimating sea surface chlorophyll and sediments from thematic mapper imagery[J]. Remote Sensing of Environment, 1998, 66(2): 153-165.

    [15] Tanaka A, Kishino M, Doerffer R, et al. Development of a neural network algorithm for retrieving concentrations of chlorophyll, suspended matter and yellow substance from radiance data of the ocean color and temperature scanner[J]. Journal of Oceanography, 2004, 60(3): 519-530.

    [16] 陳莉瓊, 田禮喬, 邱鳳, 等. HJ-1A/B衛(wèi)星CCD影像的武漢市東湖水色三要素遙感研究[J]. 武漢大學(xué)學(xué)報(bào)·信息科學(xué)版, 2011, 36(11): 1280-1283, 1335.

    Chen Liqiong, Tian Liqiao, Qiu Feng, et al. Water color constituents remote sensing in Wuhan Donghu Lake using HJ-1A/B CCD imagery[J]. Geomatics and Information Science of Wuhan University, 2011, 36(11): 1280-1283, 1335.

    [17] 張彥喆. 渤海海域葉綠素濃度反演方法研究[D]. 天津: 天津科技大學(xué), 2010.

    Zhang Yanzhe. Study on inversion method of chlorophyll concentration in the Bohai Sea[D]. Tianjin: Tianjin University of Science and Technology, 2010.

    [18] 吳倩, 林蕾, 王學(xué)軍, 等. 福海葉綠素含量的人工神經(jīng)網(wǎng)絡(luò)反演模型[J]. 地理與地理信息科學(xué), 2004, 20(4): 27-30.

    Wu Qian, Lin Lei, Wang Xuejun, et al. The neural network model on chlorophyll concentration in the Fuhai Lake[J]. Geography and Geo-Information Science, 2004, 20(4): 27-30.

    [19] 楊燕明, 劉貞文, 陳本清, 等. 用偏最小二乘法反演二類水體的水色要素[J]. 遙感學(xué)報(bào), 2005, 9(2): 123-130.

    Yang Yanming, Liu Zhenwen, Chen Benqing, et al. Retrieval of oceanic color constituents from case Ⅱ water reflectance by partial least squares regression[J]. Journal of Remote Sensing, 2005, 9(2): 123-130.

    [20] 楊曦光. 黃海葉綠素及初級(jí)生產(chǎn)力的遙感估算[D]. 青島: 中國(guó)科學(xué)院海洋研究所, 2013.

    Yang Xiguang. Chlorophyll and primary productivity estimation based on remote sensing in Yellow Sea[D]. Qingdao: Institute of Oceanology, Chinese Academy of Sciences, 2013.

    [21] Yapa K K A S. Seasonal variability of sea surface chlorophyll-aof waters around Sri Lanka[J]. Journal of Earth System Science, 2000, 109(4): 427-432.

    [22] Goes J I, Oaku T S H, Liang D L. A method for estimating sea surface nitrate concentrations from remotely sensed SST and chlorophylla-a case study for the north Pacific Ocean using OCTS/ADEOS data[J]. IEEE Transactions on Geoscience and Remote Sensing, 1999, 37(3): 1633-1644.

    [23] Nagamani P V, Chauhan P, Dwivedi R M. Estimation of chlorophyll-aconcentration using an artificial neural network (ANN)-based algorithm with oceansat-I OCM data[J]. Journal of the Indian Society of Remote Sensing, 2007, 35(3): 201-207.

    [24] Sancak S, Besiktepe S T, Yilmaz A, et al. Evaluation of SeaWiFS chlorophyll-ain the Black and Mediterranean Seas[J]. International Journal of Remote Sensing, 2005, 26(10): 2045-2060.

    [25] Dogliotti A I, Schloss I R, Almandoz G O, et al. Evaluation of SeaWiFS and MODIS Chlorophyll-aProducts in the Argentinean Patagonian Continental Shelf (38°S-55°S)[J]. International Journal of Remote Sensing, 2009, 30(1): 251-273.

    [26] Vaiciūtè D, Olenina I, Kavolytè T, et al. Validation of MERIS chlorophyllaproducts in the Lithuanian Baltic Sea case 2 coastal waters[C]// 2010 IEEE/OES US/EU BALTIC International Symposium. Riga: Institute of Electrical and Electronics Engineers, 2010.

    [27] Carder K L, Chen F R, Cannizzaro J P, et al. Performance of the MODIS semi-analytical ocean color algorithm for chlorophyll-a[J]. Advances in Space Research, 2004,33(2):1152-1159

    [28] 丁曉英, 陳曉翔. 基于環(huán)境因素的沿岸水域葉綠素遙感探測(cè)研究[J]. 遙感學(xué)報(bào), 2005, 9(4): 446-451.

    Ding Xiaoying, Chen Xiaoxiang. The study of chlorophyll detection in coastal waters based on environmental factors[J]. Journal of Remote Sensing, 2005, 9(4): 446-451.

    [29] 沈春燕, 陳楚群, 詹海剛. 人工神經(jīng)網(wǎng)絡(luò)反演珠江口海域葉綠素濃度[J]. 熱帶海洋學(xué)報(bào), 2005, 24(6): 38-43.

    Shen Chunyan, Chen Chuqun, Zhan Haigang. Inverse of chlorophyll concentration in Zhujiang River estuary using artifical neural network[J]. Journal of Tropical Oceanography, 2005, 24(6): 38-43.

    [30] 于堃, 陸殿梅, 熊黑鋼. 近7年來(lái)渤海海區(qū)冬季表層海水葉綠素濃度的遙感反演及其變化分析[J]. 遙感信息, 2009(6): 55-62.

    Yu Kun, Lu Dianmei, Xiong Heigang. Analysis and retrieving of surface layer chlorophyll concentration of Bohai Sea in past 7 years[J]. Remote Sensing Application, 2009(6): 55-62.

    [31] 萬(wàn)幼川, 黃家潔, 劉良明. 基于MODIS數(shù)據(jù)的二類水體葉綠素反演算法理論與實(shí)現(xiàn)[J]. 武漢大學(xué)學(xué)報(bào)·信息科學(xué)版, 2007, 32(7): 572-575.

    Wan Youchuan, Huang Jiajie, Liu Liangming. Theory and implementation of chlorophyll concentration retrieval algorithm for case 2 water with Modis[J]. Geomatics and Information Science of Wuhan University, 2007, 32(7): 572-575.

    [32] 沙慧敏, 李小恕, 楊文波, 等. 用MODIS遙感數(shù)據(jù)反演東海海表溫度、葉綠素a濃度年際變化的研究[J]. 大連水產(chǎn)學(xué)院學(xué)報(bào), 2009, 24(2): 151-156.

    Sha Huimin, Li Xiaoshu, Yang Wenbo, et al. Annual variation in sea surface temperature and chlorophyll-aconcentration retrieved by MODIS in East China Sea[J]. Journal of Dalian Fisheries University, 2009, 24(2): 151-156.

    [33] 李新星, 張亭祿, 田林, 等. 多衛(wèi)星傳感器南海葉綠素a濃度遙感數(shù)據(jù)融合[J]. 遙感學(xué)報(bào), 2015, 19(4): 680-689.

    Li Xinxing, Zhang Tinglu, Tian Lin, et al. Merging chlorophyll-adata from multiple ocean color sensors in South China Sea[J]. Journal of Remote Sensing, 2015, 19(4): 680-689.

    [34] 于堃, 陸殿梅. 二類水體的MODIS數(shù)據(jù)大氣校正及在渤海表層水體葉綠素濃度監(jiān)測(cè)反演中的應(yīng)用研究[J]. 海洋通報(bào), 2009, 28(5): 21-27.

    Yu Kun, Lu Dianmei. Atmospheric correction of modis data for case-ii water and research of monitoring and retrieving chlorophyll concentration in the Bohai Sea[J]. Marine Science Bulletin, 2009, 28(5): 21-27.

    [35] 張麗華, 戴學(xué)芳, 包玉海, 等. 基于TM影像的烏梁素海葉綠素a濃度反演[J]. 環(huán)境工程, 2015, 33(6): 133-138.

    Zhang Lihua, Dai Xuefang, Bao Yuhai, et al. Inversion of chlorophyll-aconcentration based on tm remote sensing image in Wuliangsuhai Lake[J]. Environmental Monitoring & Assessment, 2015, 33(6): 133-138.

    [36] 張宏亮. FLH算法反演渤海葉綠素濃度的能力評(píng)價(jià)研究[D]. 呼和浩特: 內(nèi)蒙古大學(xué), 2007.

    Zhang Hongliang. Evaluation research on chlorophyll concentration retrieval by FLH algorithm in Bohai Sea[D]. Huhhot: Inner Mongolia University, 2007.

    [37] 鄭小慎, 敖翔. 利用遙感數(shù)據(jù)反演渤海海域葉綠素的濃度[J]. 天津科技大學(xué)學(xué)報(bào), 2009, 24(2): 35-37, 69.

    Zheng Xiaoshen, Ao Xiang. Inversion of chlorophyll concentrations in Bohai Sea using remote sensing data[J]. Journal of Tianjin University of Science & Technology, 2009, 24(2): 35-37, 69.

    [38] 鄭國(guó)強(qiáng), 史同廣, 孫林, 等. 基于Hyperion數(shù)據(jù)的南四湖葉綠素濃度反演研究[J]. 地理與地理信息科學(xué), 2008, 24(1): 31-34.

    Zheng Guoqiang, Shi Tongguang, Sun Lin, et al. Regression of chlorophyll content based Hyperion data in Nansi Lake[J]. Geography and Geo-Information Science, 2008, 24(1): 31-34.

    [39] 聞建光, 肖青, 楊一鵬, 等. 基于Hyperion數(shù)據(jù)的太湖水體葉綠素a濃度遙感估算[J]. 湖泊科學(xué), 2006, 18(4): 327-336.

    Wen Jianguang, Xiao Qing, Yang Yipeng, et al. Remote sensing estimation of aquatic chlorophyll-aconcentration based on Hyperion data in Lake Taihu[J]. Journal of Lake Sciences, 2006, 18(4): 327-336.

    [40] 錢莉. 渤海海域赤潮監(jiān)測(cè)中葉綠素濃度反演算法的研究[D]. 天津: 天津科技大學(xué), 2011.

    Qian Li. Research on inversion algorithms of chlorophyll concentration about red tides monitoring in Bohai Sea[D]. Tianjin: Tianjin University of Science and Technology, 2011.

    [41] 郭宇龍. 面向內(nèi)陸水色遙感的多源遙感影像融合方法研究[D]. 南京: 南京師范大學(xué), 2015.

    Guo Yulong. Research of multisource image fusion algorithm for inland water color remote sensing[D]. Nanjing: Nanjing Normal University, 2015.

    [42] 鄒亞榮. 渤海葉綠素a時(shí)空分布特征分析[J]. 遙感信息, 2004, 1(3): 30-31.

    Zou Yarong. Analysis spatial distribution characteristics for chlain Bohai Sea[J]. Remote Sensing Information, 2004, 1(3): 30-31.

    [43] He Xianqiang, Bai Yan, Pan Delu, et al. The atmospheric correction algorithm for HY-1B/COCTS[C]// Geoinformatics 2008 and Joint Conference on GIS and Built. Washington: SPIE, 2008.

    [44] 林珊, 張玉超, 錢新, 等. 廣州流溪河水庫(kù)葉綠素a遙感反演研究[J]. 環(huán)境保護(hù)科學(xué), 2010, 36(3): 43-46, 104.

    Lin Shan, Zhang Yuchao, Qian Xin, et al. Chlorophyllaconcentration retrieval of Liuxihe Reservoir in Guangzhou using remote sensing data[J]. Environmental Protection Science, 2010, 36(3): 43-46, 104.

    [45] 王貴臣, 劉田田. 基于線性回歸分析的東平湖葉綠素含量監(jiān)測(cè)[J]. 首都師范大學(xué)學(xué)報(bào)(自然科學(xué)版), 2015, 36(2): 72-77.

    Wang Guichen, Liu Tiantian. Chlorophyllaconcentration monitoring of Dongping Lake based on linear regression analysis[J]. Journal of Capital Normal University (Natural Sciences Edition), 2015, 36(2): 72-77.

    [46] 雷坤, 鄭丙輝, 王橋. 基于中巴地球資源1號(hào)衛(wèi)星的太湖表層水體水質(zhì)遙感[J]. 環(huán)境科學(xué)學(xué)報(bào), 2004, 24(3): 376-380.

    Lei Kun, Zheng Binghui, Wang Qiao. Monitoring the surface water quality of Taihu Lake based on the data of CBERS-1[J]. Acta Scientiae Circumstantiae, 2004, 24(3): 376-380.

    [47] 朱利, 李云梅, 趙少華, 等. 基于GF-1號(hào)衛(wèi)星WFV數(shù)據(jù)的太湖水質(zhì)遙感監(jiān)測(cè)[J]. 國(guó)土資源遙感, 2015, 27(1): 113-120.

    Zhu Li, Li Yunmei, Zhao Shaohua, et al. Remote sensing monitoring of Taihu Lake water quality by using GF-1 satellite WFV data[J]. Remote Sensing for Land & Resources, 2015, 27(1): 113-120.

    [48] 郭宇龍, 李云梅, 朱利, 等. 基于HJ1A-CCD數(shù)據(jù)的高光譜影像重構(gòu)研究[J]. 環(huán)境科學(xué), 2013, 34(1): 69-76.

    Guo Yulong, Li Yunmei, Zhu Li, et al. Research of hyperspectral reconstruction based on HJ1A-CCD data[J]. Environmental Science, 2013, 34(1): 69-76.

    [49] 王琦, 孟偉, 馬云峰, 等. 基于HJ-1衛(wèi)星的大伙房水庫(kù)葉綠素a濃度反演模型研究[J]. 安全與環(huán)境學(xué)報(bào), 2013, 13(4): 137-141.

    Wang Qi, Meng Wei, Ma Yunfeng, et al. On the inversion models for chlorophyll-aconcentration based on the HJ-1 satellite images of Dahuofang Reservoir, Liaoning[J]. Journal of Safety and Environment, 2013, 13(4): 137-141.

    [50] 安如, 劉影影, 曲春梅, 等. NDCI法Ⅱ類水體葉綠素a濃度高光譜遙感數(shù)據(jù)估算[J]. 湖泊科學(xué), 2013, 25(3): 437-444.

    An Ru, Liu Yingying, Qu Chunmei, et al. Estimation of chlorophyll-aconcentration of case Ⅱ waters from hyperspectral remote sensing data in NDCI method[J]. Journal of Lake Sciences, 2013, 25(3): 437-444.

    [51] 劉王兵, 于之鋒, 周斌, 等. 杭州灣HJ CCD影像懸浮泥沙遙感定量反演[J]. 遙感學(xué)報(bào), 2013, 17(4): 905-918.

    Liu Wangbing, Yu Zhifeng, Zhou Bin, et al. Assessment of suspended sediment concentration at the Hangzhou Bay using HJ CCD imagery[J]. Journal of Remote Sensing, 2013, 17(4): 905-918.

    [52] 吳川, 張玉龍, 張克榮, 等. 丹江口水庫(kù)水質(zhì)的遙感監(jiān)測(cè)[J]. 南水北調(diào)與水利科技, 2013, 11(6): 75-80.

    Wu Chuan, Zhang Yulong, Zhang Kerong, et al. Remote sensing monitoring of water quality in the Danjiangkou Reservoir[J]. South-to-North Water Transfers and Water Science & Technology, 2013, 11(6): 75-80.

    [53] 杜挺, 王旭紅, 梁海艷, 等. 湖泊葉綠素a反演的大氣校正模型比較研究[J]. 云南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2014, 36(2): 241-248.

    Du Ting, Wang Xuhong, Liang Haiyan, et al. Comparative study on atmospheric correction methods in lakes chlorophyll-aretrieval[J]. Journal of Yunnan University, 2014, 36(2): 241-248.

    [54] 曹曉峰. 基于HJ-1A/1B影像的滇池水質(zhì)遙感監(jiān)測(cè)研究[D]. 西安: 西安科技大學(xué), 2012.

    Cao Xiaofeng. The study of water quality monitoring in Dianchi Lake based on HJ-1A/1B remote sensing data[D]. Xi’an: Xi’an University of Science and Technology, 2012.

    [55] 周正, 何連, 劉良明. 基于HJ-1A/B CCD數(shù)據(jù)的東湖葉綠素a濃度反演可行性研究[J]. 測(cè)繪通報(bào), 2011(3): 11-14.

    Zhou Zheng, He Lian, Liu Liangming. Feasibility study for retrieval of chlorophyll-aconcentration in Donghu Lake using HJ-1 A/B CCD data[J]. Bulletin of Surveying and Mapping, 2011(3): 11-14.

    [56] 楊榮, 楊昆, 洪亮, 等. 基于HJ-1A衛(wèi)星CCD1數(shù)據(jù)的滇池水質(zhì)參數(shù)反演模型研究[J]. 節(jié)水灌溉, 2015(3): 42-45.

    Yang Rong, Yang Kun, Hong Liang, et al. A study of inversion modeling of water quality parameters in Dianchi using HJ-1A CCD1 data[J]. Water Saving Irrigation, 2015(3): 42-45.

    [57] 朱玉霞, 張杰, 劉銳, 等. 基于HJ-1衛(wèi)星的丹江口庫(kù)區(qū)水質(zhì)遙感監(jiān)測(cè)研究[J]. 環(huán)境科技, 2014, 27(5): 52-58.

    Zhu Yuxia, Zhang Jie, Liu Rui, et al. Study on remote sensing monitoring of water quality for Danjiangkou Reservoir by HJ-1 satellite data[J]. Environmental Science and Technology, 2014, 27(5): 52-58.

    [58] 劉朝相, 宮兆寧, 趙文吉, 等. 基于SVM模型的媯水河葉綠素a濃度的遙感反演[J]. 遙感技術(shù)與應(yīng)用, 2014, 29(3): 419-427.

    Liu Chaoxiang, Gong Zhaoning, Zhao Wenji, et al. Remote sensing retrieval of chlorophyll-aconcentration in Beijing Guishuihe River using support vector machine model[J]. Remote Sensing Technology and Application, 2014, 29(3): 419-427.

    [59] 潘梅娥, 楊昆, 洪亮. 基于環(huán)境一號(hào)衛(wèi)星影像的內(nèi)陸水體葉綠素a濃度遙感定量反演模型研究[J]. 科學(xué)技術(shù)與工程, 2013, 13(15): 4228-4233.

    Pan Mei’e, Yang Kun, Hong Liang. Remote sensing quantitative retrieval models research of inland water chlorophyll-aconcentration based on HJ-1 satellite images[J]. Science Technology and Engineering, 2013, 13(15): 4228-4233.

    [60] 王珊珊, 李云梅, 王永波, 等. 太湖水體葉綠素濃度反演模型適宜性分析[J]. 湖泊科學(xué), 2015, 27(1): 150-162.

    Wang Shanshan, Li Yunmei, Wang Yongbo, et al. Suitability of the retrieval models for estimating chlorophyll concentration in Lake Taihu[J]. Journal of Lake Sciences, 2015, 27(1): 150-162.

    [61] 鄭高強(qiáng), 陳蕓芝, 汪小欽, 等. 基于HJ-1衛(wèi)星CCD數(shù)據(jù)的廈門海域葉綠素a濃度反演[J]. 遙感技術(shù)與應(yīng)用, 2015, 30(2): 235-241, 284.

    Zheng Gaoqiang, Chen Yunzhi, Wang Xiaoqin, et al. The chlorophyll-aconcentration extraction of Xiamen Sea based on the HJ-1 satellite image[J]. Remote Sensing Technology and Application, 2015, 30(2): 235-241, 284.

    [62] 王皓. 基于HJ-1星數(shù)據(jù)的大洋河河口水域水質(zhì)參數(shù)反演研究[D]. 大連: 遼寧師范大學(xué), 2012.

    Wang Hao. Study of water quality index retrievals in Dayang River estuary based on HJ-1 satellite remote-sensed imagery[D]. Dalian: Liaoning Normal University, 2012.

    [63] 孟慶輝, 王林, 毛志華, 等. 基于HJ CCD的杭州灣海域水質(zhì)評(píng)價(jià)[J]. 廣西科學(xué), 2015, 22(3): 322-328, 336.

    Meng Qinghui, Wang Lin, Mao Zhihua, et al. Preliminary study on the application of HJ CCD imagery to water quality assessment for Hangzhou Bay[J]. Guangxi Sciences, 2015, 22(3): 322-328, 336.

    [64] 秦平, 沈鉞, 牟冰, 等. 基于進(jìn)化建模方法的HJ-1 CCD黃海懸浮物和葉綠素a濃度遙感反演模型研究[J]. 海洋學(xué)報(bào), 2014, 36(11): 142-149.

    Qin Ping, Shen Yue, Mu Bing, et al. Retrieval models of totalsuspended matter and chlorophyllaconcentration in Yellow Sea based on HJ-1 CCD data and evolutionary modeling method[J]. Haiyang Xuebao, 2014, 36(11): 142-149.

    [65] 楊婷, 張慧, 王橋, 等. 基于HJ-1A衛(wèi)星超光譜數(shù)據(jù)的太湖葉綠素a濃度及懸浮物濃度反演[J]. 環(huán)境科學(xué), 2011, 32(11): 3207-3214.

    Yang Ting, Zhang Hui, Wang Qiao, et al. Retrieving for chlorophyll-aconcentration and suspended substance concentration based on HJ-1A HIS image[J]. Environmental Science, 2011, 32(11): 3207-3214.

    [66] 叢丕福, 牛錚, 蒙繼華, 等. 1998~2003年衛(wèi)星反演的中國(guó)陸架海葉綠素a濃度變化分析[J]. 海洋環(huán)境科學(xué), 2006, 25(1): 30-33.

    Cong Pifu, Niu Zheng, Meng Jihua, et al. Variability of chlorophyllaretrieved from satellite in Chinese shelf sea from 1998 to 2003[J]. Marine Environmental Science, 2006, 25(1): 30-33.

    [67] 韓曉慶, 蘇藝, 李靜, 等. 海岸帶地區(qū)SPOT衛(wèi)星影像大氣校正方法比較及精度驗(yàn)證[J]. 地理研究, 2012, 31(11): 2007-2016.

    Han Xiaoqing, Su Yi, Li Jing, et al. Atmospheric correction and verification of the SPOT remote sensing image in coastal zones[J]. Geographical Research, 2012, 31(11): 2007-2016.

    [68] 鄧書斌. ENVI遙感圖像處理方法[M]. 北京: 科學(xué)出版社, 2010.

    Deng Shubin. ENVI Remote Sensing Image Processing Method[M]. Beijing: Science Press, 2010.

    [69] 王玉玨, 劉哲, 張永, 等. 2010-2011年膠州灣葉綠素a與環(huán)境因子的時(shí)空變化特征[J]. 海洋學(xué)報(bào), 2015, 37(4): 103-116.

    Wang Yujue, Liu Zhe, Zhang Yong, et al. Temporal and spatial variations of chlorophyllaand environmental factors in Jiaozhou Bay in 2010-2011[J]. Haiyang Xuebao, 2015, 37(4): 103-116.

    Inversion of chlorophyllaconcentration in offshore Ⅱ waters using HJ satellite data——Example in the north of the Luanhe Delta

    Luo Jianmei1,2,3, Huo Yongwei4, Han Xiaoqing4

    (1.CollegeofLandResourcesandRural-UrbanPlanning,HebeiGEOUniversity,Shijiazhuang050031,China; 2.KeyLaboratoryofAgriculturalWaterResources,CenterforAgriculturalResourcesResearch,InstituteofGeneticsandDevelopmentalBiology,ChineseAcademyofSciences,Shijiazhuang050022,China; 3.UniversityofChineseAcademyofSciences,Beijing100049,China;4.HebeiUtilizationandPlanningInstituteofLandResources,Shijiazhuang050051,China)

    In order to verify the suitability of environmental satellite images being used to conduct concentration inversion of chlorophylla(Chla) in nearshore Ⅱ waters, and to obtain the inversion algorithm in the north nearshore water areas of Luanhe River, several inversion models of Chlaconcentration were established based on the HJ-1A CCD2 images and synchronous measurement data. By comparing the simulated results with the validation data, the most prosperity model to simulate Chlaconcentration for study area was selected. The results show that HJ-1A CCD2 images are effective enough to inverse the Chlaconcentration of nearshore Ⅱ waters under the complex hydrological situation. We therefore built the reciprocal model by B3/B1 which was the optimal band combination for the study area. Inversion results show that the Chlaconcentration gradually decreases from the nearshore to the inner sea, and the same is true from the north to the south, thus forming a high nutrient area to the sea outside the estuary. The maximum value of 25.0 μg/L occurs in the south of Jinshanzui, the minimum value of 4.6 μg/L in offshore areas occurs at about 13.2 km east of the Luanhe Delta, while the minimum value of 11.0 μg/L in nearshore areas occurs nearby the Luanhe Delta. The methods and processes analysis in this study can be used as references for the similar researches. In addition, this study can provide decision-making references for the local government to manage marine environment and formulate related policies, and the important supporting data for the sustainable development in this area.

    HJ satellite;the north of the Luanhe Delta;chlorophyllaconcentration;model;inversion

    10.3969/j.issn.0253-4193.2017.04.012

    2016-06-02;

    2016-08-06。

    河北省科技計(jì)劃項(xiàng)目(15273302D)。

    羅建美(1978—),女,山西省大同市人,副教授,主要從事資源環(huán)境與城鄉(xiāng)規(guī)劃管理研究。E-mail:jm3003@126.com

    *通信作者:韓曉慶(1981—),男,工程師,主要從事基于遙感的環(huán)境演變研究。E-mail:xqhan312@126.com

    X87

    A

    0253-4193(2017)04-0117-13

    羅建美,霍永偉,韓曉慶. 基于HJ衛(wèi)星的近岸Ⅱ類水體葉綠素a濃度定量遙感反演研究——以灤河口北部海域?yàn)槔齕J].海洋學(xué)報(bào),2017,39(4):117—129,

    Luo Jianmei, Huo Yongwei, Han Xiaoqing. Inversion of chlorophyllaconcentration in offshore Ⅱ waters using HJ satellite data——Example in the north of the Luanhe Delta[J]. Haiyang Xuebao,2017,39(4):117—129, doi:10.3969/j.issn.0253-4193.2017.04.012

    猜你喜歡
    波段葉綠素反演
    春日暖陽(yáng)
    反演對(duì)稱變換在解決平面幾何問(wèn)題中的應(yīng)用
    提取葉綠素
    桃樹(shù)葉綠素含量與SPAD值呈極顯著正相關(guān)
    基于低頻軟約束的疊前AVA稀疏層反演
    基于自適應(yīng)遺傳算法的CSAMT一維反演
    葉綠素家族概述
    M87的多波段輻射過(guò)程及其能譜擬合
    日常維護(hù)對(duì)L 波段雷達(dá)的重要性
    西藏科技(2015年4期)2015-09-26 12:12:58
    由松針制取三種葉綠素鈉鹽及其穩(wěn)定性的研究

    黑丝袜美女国产一区| 91精品国产国语对白视频| 少妇熟女欧美另类| 在线播放无遮挡| 欧美xxⅹ黑人| 日日撸夜夜添| 丰满少妇做爰视频| 熟妇人妻不卡中文字幕| tube8黄色片| 高清在线视频一区二区三区| 精品亚洲乱码少妇综合久久| 成人综合一区亚洲| 蜜桃在线观看..| 中文字幕久久专区| 在线播放无遮挡| 国产精品一区二区三区四区免费观看| 大陆偷拍与自拍| 欧美精品一区二区大全| 美女大奶头黄色视频| 中文欧美无线码| 精品一区二区免费观看| 97超视频在线观看视频| 免费不卡的大黄色大毛片视频在线观看| 一本—道久久a久久精品蜜桃钙片| 国国产精品蜜臀av免费| 国产在视频线精品| 久久精品国产鲁丝片午夜精品| 欧美老熟妇乱子伦牲交| 成人亚洲精品一区在线观看| 日韩,欧美,国产一区二区三区| 精品久久久噜噜| 国产成人freesex在线| 久久人妻熟女aⅴ| 亚洲精品亚洲一区二区| 欧美精品一区二区免费开放| 18禁裸乳无遮挡动漫免费视频| 国产精品久久久久久av不卡| 欧美精品高潮呻吟av久久| 精品国产一区二区三区久久久樱花| 国产极品粉嫩免费观看在线 | av国产精品久久久久影院| 妹子高潮喷水视频| 99re6热这里在线精品视频| 最近中文字幕2019免费版| 久久这里有精品视频免费| av卡一久久| 亚洲精品日本国产第一区| 亚洲精品456在线播放app| 久久久a久久爽久久v久久| 80岁老熟妇乱子伦牲交| 成人免费观看视频高清| 在线观看av片永久免费下载| 看非洲黑人一级黄片| 亚洲av福利一区| 国产黄片视频在线免费观看| av不卡在线播放| 国产无遮挡羞羞视频在线观看| 欧美精品高潮呻吟av久久| 欧美激情国产日韩精品一区| 欧美三级亚洲精品| 日本av免费视频播放| av播播在线观看一区| 日韩不卡一区二区三区视频在线| 亚洲精品乱久久久久久| 亚洲国产精品一区三区| 日本av手机在线免费观看| 国产精品伦人一区二区| 亚洲欧美日韩东京热| 日本色播在线视频| 大片免费播放器 马上看| 国产深夜福利视频在线观看| 亚洲av日韩在线播放| 国产男女内射视频| 51国产日韩欧美| 日韩电影二区| 男女啪啪激烈高潮av片| 国产精品一区二区在线不卡| a级毛色黄片| 欧美+日韩+精品| av视频免费观看在线观看| 国产精品一区二区在线观看99| 看非洲黑人一级黄片| av卡一久久| 中文在线观看免费www的网站| 久久人妻熟女aⅴ| 国产精品国产三级国产av玫瑰| 日韩欧美精品免费久久| 国产欧美日韩精品一区二区| 中文欧美无线码| 一区二区三区四区激情视频| 中文字幕免费在线视频6| 性高湖久久久久久久久免费观看| 亚洲av福利一区| 欧美 亚洲 国产 日韩一| 男人添女人高潮全过程视频| 国产一区二区三区av在线| 国产黄片视频在线免费观看| 99热这里只有是精品50| 免费人妻精品一区二区三区视频| 日本av手机在线免费观看| 精品99又大又爽又粗少妇毛片| 日韩熟女老妇一区二区性免费视频| 99久久精品热视频| 欧美xxxx性猛交bbbb| a级片在线免费高清观看视频| 欧美精品一区二区免费开放| 美女福利国产在线| 成人综合一区亚洲| 亚洲精品日韩av片在线观看| 国产在视频线精品| 人人妻人人看人人澡| 亚洲精品,欧美精品| 人人妻人人添人人爽欧美一区卜| 老熟女久久久| 国产av码专区亚洲av| 国产成人精品无人区| 又大又黄又爽视频免费| 久久久久久久精品精品| av在线app专区| 永久免费av网站大全| av一本久久久久| 国产无遮挡羞羞视频在线观看| 女人久久www免费人成看片| 亚洲激情五月婷婷啪啪| 亚洲av福利一区| 中文乱码字字幕精品一区二区三区| 国产精品伦人一区二区| 夜夜看夜夜爽夜夜摸| 精品午夜福利在线看| 亚洲人与动物交配视频| 岛国毛片在线播放| 色哟哟·www| 久久99热6这里只有精品| 免费大片黄手机在线观看| 国产成人91sexporn| 日韩欧美精品免费久久| 亚洲成人一二三区av| 久久精品国产a三级三级三级| 中文字幕久久专区| 久久鲁丝午夜福利片| 看十八女毛片水多多多| 99视频精品全部免费 在线| 成年美女黄网站色视频大全免费 | 在线观看av片永久免费下载| 日韩大片免费观看网站| 亚洲色图综合在线观看| 麻豆乱淫一区二区| 人妻 亚洲 视频| 中文字幕人妻丝袜制服| 成人国产麻豆网| 视频中文字幕在线观看| 国产美女午夜福利| 在线观看一区二区三区激情| h日本视频在线播放| av免费观看日本| 欧美另类一区| 欧美日韩国产mv在线观看视频| 人妻少妇偷人精品九色| 亚洲精品国产av蜜桃| 色5月婷婷丁香| 亚洲av国产av综合av卡| 中文精品一卡2卡3卡4更新| 欧美日韩国产mv在线观看视频| 新久久久久国产一级毛片| 亚洲av欧美aⅴ国产| 综合色丁香网| av在线老鸭窝| av免费观看日本| 高清欧美精品videossex| 欧美变态另类bdsm刘玥| 亚洲av男天堂| 国产精品久久久久久久久免| 色94色欧美一区二区| 亚洲国产av新网站| 人体艺术视频欧美日本| 亚洲无线观看免费| 国产精品成人在线| 久久人人爽av亚洲精品天堂| 妹子高潮喷水视频| 成人亚洲精品一区在线观看| 国产精品国产三级国产av玫瑰| 亚洲四区av| videos熟女内射| 两个人免费观看高清视频 | 少妇人妻久久综合中文| 女人精品久久久久毛片| 热re99久久国产66热| 自线自在国产av| 女性生殖器流出的白浆| av福利片在线| 嫩草影院新地址| 韩国高清视频一区二区三区| 精品一品国产午夜福利视频| 日本黄大片高清| 中文字幕亚洲精品专区| 人妻人人澡人人爽人人| 亚洲av不卡在线观看| 亚洲在久久综合| 亚州av有码| 七月丁香在线播放| av卡一久久| 亚洲精品456在线播放app| 亚洲精品一二三| 国产一区二区三区综合在线观看 | 你懂的网址亚洲精品在线观看| 我的女老师完整版在线观看| 日本黄色片子视频| 亚洲欧美精品自产自拍| 天堂8中文在线网| 内地一区二区视频在线| 午夜日本视频在线| 午夜福利视频精品| 久热这里只有精品99| 少妇人妻一区二区三区视频| 亚洲四区av| 国产伦理片在线播放av一区| 少妇人妻 视频| 婷婷色综合www| 夜夜爽夜夜爽视频| 婷婷色综合大香蕉| 狂野欧美白嫩少妇大欣赏| 成年人免费黄色播放视频 | 亚洲av综合色区一区| 国产精品一区二区在线观看99| 久久久国产一区二区| 欧美3d第一页| 国产av精品麻豆| 久久久精品94久久精品| 国产日韩欧美亚洲二区| 欧美日韩视频精品一区| 亚洲欧美一区二区三区国产| 欧美区成人在线视频| 91久久精品国产一区二区三区| 精品99又大又爽又粗少妇毛片| 国产免费一级a男人的天堂| 我要看日韩黄色一级片| 国产乱人偷精品视频| 亚洲精品自拍成人| 日本欧美国产在线视频| 日韩熟女老妇一区二区性免费视频| 久久人人爽人人爽人人片va| 男人狂女人下面高潮的视频| 欧美精品一区二区大全| 亚洲婷婷狠狠爱综合网| 欧美+日韩+精品| 国产伦精品一区二区三区四那| 少妇人妻一区二区三区视频| 日韩 亚洲 欧美在线| 三级经典国产精品| 亚洲国产精品一区三区| 久久精品熟女亚洲av麻豆精品| 亚洲精品日韩av片在线观看| 欧美xxxx性猛交bbbb| 色94色欧美一区二区| 国国产精品蜜臀av免费| 91成人精品电影| 国产精品成人在线| 久久久久精品久久久久真实原创| 国产乱来视频区| 99精国产麻豆久久婷婷| 久久精品久久精品一区二区三区| 国产伦精品一区二区三区四那| 日韩精品有码人妻一区| 丰满饥渴人妻一区二区三| 你懂的网址亚洲精品在线观看| 亚洲情色 制服丝袜| 涩涩av久久男人的天堂| 日韩中字成人| a级毛色黄片| 久久人妻熟女aⅴ| 少妇裸体淫交视频免费看高清| 超碰97精品在线观看| 免费观看在线日韩| 国产成人a∨麻豆精品| 国产av一区二区精品久久| 精品国产一区二区久久| 妹子高潮喷水视频| 久久免费观看电影| 热99国产精品久久久久久7| 日本黄色日本黄色录像| 99九九线精品视频在线观看视频| 99久久精品一区二区三区| 少妇人妻 视频| 人妻系列 视频| 夜夜看夜夜爽夜夜摸| 校园人妻丝袜中文字幕| 国产又色又爽无遮挡免| 大片电影免费在线观看免费| 精品亚洲成a人片在线观看| 欧美日韩在线观看h| 国产片特级美女逼逼视频| 免费大片黄手机在线观看| 2018国产大陆天天弄谢| 免费黄频网站在线观看国产| 精品国产一区二区久久| 狠狠精品人妻久久久久久综合| 日韩av不卡免费在线播放| 好男人视频免费观看在线| 国产伦理片在线播放av一区| 亚洲精品国产成人久久av| av在线老鸭窝| 国产成人91sexporn| 亚洲av综合色区一区| 天天操日日干夜夜撸| 嫩草影院入口| 精品酒店卫生间| 久久ye,这里只有精品| 国产一区有黄有色的免费视频| 亚洲美女黄色视频免费看| 国产精品嫩草影院av在线观看| 精品酒店卫生间| 偷拍熟女少妇极品色| 欧美97在线视频| 国内揄拍国产精品人妻在线| 香蕉精品网在线| 男人舔奶头视频| 欧美 亚洲 国产 日韩一| 精品酒店卫生间| 国产精品女同一区二区软件| 91精品一卡2卡3卡4卡| 中文精品一卡2卡3卡4更新| 蜜桃久久精品国产亚洲av| 一区二区三区精品91| av女优亚洲男人天堂| 色94色欧美一区二区| 久久久久视频综合| 极品教师在线视频| 久久99精品国语久久久| 69精品国产乱码久久久| 日韩一本色道免费dvd| 在线天堂最新版资源| av视频免费观看在线观看| 我的老师免费观看完整版| 久久人人爽人人爽人人片va| 国产中年淑女户外野战色| 在线观看一区二区三区激情| 亚洲av欧美aⅴ国产| 久久影院123| 美女视频免费永久观看网站| 人人妻人人澡人人爽人人夜夜| 菩萨蛮人人尽说江南好唐韦庄| 国产日韩一区二区三区精品不卡 | 草草在线视频免费看| 日本av手机在线免费观看| 亚洲av免费高清在线观看| 在线观看免费视频网站a站| 中文欧美无线码| 黄色毛片三级朝国网站 | 啦啦啦在线观看免费高清www| 热re99久久精品国产66热6| 夫妻性生交免费视频一级片| 精品久久久久久久久av| 天堂俺去俺来也www色官网| 成人免费观看视频高清| 亚州av有码| 香蕉精品网在线| 人体艺术视频欧美日本| a级毛片免费高清观看在线播放| 中文字幕免费在线视频6| 欧美区成人在线视频| 一级毛片久久久久久久久女| videos熟女内射| 国产午夜精品久久久久久一区二区三区| 99久久综合免费| 午夜视频国产福利| 欧美bdsm另类| 简卡轻食公司| 国产黄片美女视频| 色哟哟·www| 黑人巨大精品欧美一区二区蜜桃 | 精品99又大又爽又粗少妇毛片| 日韩av不卡免费在线播放| av在线播放精品| 亚洲av成人精品一区久久| 少妇人妻一区二区三区视频| 老司机影院成人| 99久久精品一区二区三区| 少妇丰满av| 久久久久国产精品人妻一区二区| 精品少妇黑人巨大在线播放| 成人国产av品久久久| 亚洲精品一二三| 男女啪啪激烈高潮av片| 一区在线观看完整版| 午夜老司机福利剧场| 人妻 亚洲 视频| 欧美日韩视频精品一区| 国产精品99久久99久久久不卡 | 免费看日本二区| 99国产精品免费福利视频| 亚洲美女黄色视频免费看| 亚洲成人手机| 日韩一区二区视频免费看| 日本av免费视频播放| 日日摸夜夜添夜夜爱| 国产高清三级在线| 日韩一本色道免费dvd| 18禁在线播放成人免费| 国产日韩欧美视频二区| 性高湖久久久久久久久免费观看| 国产精品人妻久久久影院| 丰满饥渴人妻一区二区三| 国产精品一区二区性色av| 色吧在线观看| 国产精品久久久久久久久免| 日本黄大片高清| 久久6这里有精品| 中文字幕久久专区| 寂寞人妻少妇视频99o| 国产日韩一区二区三区精品不卡 | 亚洲精品国产av成人精品| 黄色配什么色好看| 一级二级三级毛片免费看| 日韩电影二区| 国产伦精品一区二区三区视频9| 曰老女人黄片| 看免费成人av毛片| 一边亲一边摸免费视频| 成人二区视频| 成年美女黄网站色视频大全免费 | 亚洲国产av新网站| 国产视频内射| 大话2 男鬼变身卡| av在线播放精品| 免费看日本二区| 国产白丝娇喘喷水9色精品| 草草在线视频免费看| 我要看黄色一级片免费的| 久久久精品免费免费高清| 亚洲精品视频女| 久久韩国三级中文字幕| 亚洲丝袜综合中文字幕| 2022亚洲国产成人精品| 五月天丁香电影| 91午夜精品亚洲一区二区三区| 国产精品国产三级国产专区5o| 亚洲熟女精品中文字幕| 久久精品国产鲁丝片午夜精品| 免费不卡的大黄色大毛片视频在线观看| 国产高清三级在线| 有码 亚洲区| a级一级毛片免费在线观看| 精品一区在线观看国产| 在线观看免费视频网站a站| 亚洲av.av天堂| 成人美女网站在线观看视频| 午夜久久久在线观看| 精品99又大又爽又粗少妇毛片| 成人毛片a级毛片在线播放| 国产精品福利在线免费观看| 国产精品秋霞免费鲁丝片| 91精品国产九色| 男女免费视频国产| 久久国内精品自在自线图片| 国产高清不卡午夜福利| 亚洲色图综合在线观看| 18+在线观看网站| 国产男人的电影天堂91| 涩涩av久久男人的天堂| 久久国产乱子免费精品| 少妇熟女欧美另类| 男的添女的下面高潮视频| 国模一区二区三区四区视频| av福利片在线观看| 精品一区二区三区视频在线| 国产精品成人在线| 亚洲不卡免费看| 中文精品一卡2卡3卡4更新| 熟女av电影| 啦啦啦视频在线资源免费观看| 99视频精品全部免费 在线| 黄色视频在线播放观看不卡| 婷婷色综合大香蕉| 99热这里只有是精品50| 偷拍熟女少妇极品色| 亚洲av二区三区四区| 男人和女人高潮做爰伦理| 成人毛片60女人毛片免费| av不卡在线播放| 亚洲第一区二区三区不卡| 日韩中文字幕视频在线看片| 美女福利国产在线| 久久av网站| xxx大片免费视频| 高清在线视频一区二区三区| 九九久久精品国产亚洲av麻豆| 国产欧美日韩一区二区三区在线 | 国产男女超爽视频在线观看| 一区二区三区四区激情视频| 午夜免费男女啪啪视频观看| 大码成人一级视频| 日本91视频免费播放| 少妇被粗大猛烈的视频| 国产精品久久久久久精品古装| 黄片无遮挡物在线观看| 国产中年淑女户外野战色| 男女边摸边吃奶| 亚洲欧美中文字幕日韩二区| 久久99精品国语久久久| 精品人妻一区二区三区麻豆| 久久ye,这里只有精品| 日韩熟女老妇一区二区性免费视频| 五月天丁香电影| 如日韩欧美国产精品一区二区三区 | 欧美少妇被猛烈插入视频| 伊人久久精品亚洲午夜| 日韩亚洲欧美综合| 国产熟女午夜一区二区三区 | 三级国产精品欧美在线观看| 一级av片app| 免费不卡的大黄色大毛片视频在线观看| 国产精品人妻久久久久久| 嫩草影院入口| 午夜日本视频在线| 亚洲精品一区蜜桃| 国产高清国产精品国产三级| 大陆偷拍与自拍| 久久人人爽人人片av| 亚洲av成人精品一二三区| 午夜免费观看性视频| 91久久精品国产一区二区成人| 国产女主播在线喷水免费视频网站| 极品教师在线视频| 午夜激情久久久久久久| 黄色配什么色好看| 下体分泌物呈黄色| 日韩伦理黄色片| 黄色欧美视频在线观看| 日本爱情动作片www.在线观看| 亚洲欧美精品自产自拍| 欧美成人精品欧美一级黄| 又粗又硬又长又爽又黄的视频| 欧美日本中文国产一区发布| 老司机亚洲免费影院| 欧美精品一区二区免费开放| 又黄又爽又刺激的免费视频.| 免费大片18禁| 制服丝袜香蕉在线| 亚洲av不卡在线观看| 免费久久久久久久精品成人欧美视频 | √禁漫天堂资源中文www| 少妇高潮的动态图| 欧美xxxx性猛交bbbb| 日韩制服骚丝袜av| 国产精品不卡视频一区二区| 夜夜骑夜夜射夜夜干| 自拍偷自拍亚洲精品老妇| 精品亚洲成国产av| 在线观看美女被高潮喷水网站| 在线天堂最新版资源| 日韩一区二区三区影片| 一本色道久久久久久精品综合| 在线播放无遮挡| 亚洲成人手机| 国产亚洲5aaaaa淫片| videos熟女内射| 少妇裸体淫交视频免费看高清| 制服丝袜香蕉在线| 中文字幕制服av| 一级a做视频免费观看| 日韩熟女老妇一区二区性免费视频| 久久久久国产精品人妻一区二区| 亚洲av男天堂| 欧美人与善性xxx| 国产老妇伦熟女老妇高清| 国产一区二区在线观看日韩| 两个人免费观看高清视频 | 永久网站在线| 三上悠亚av全集在线观看 | 日本黄色日本黄色录像| 国产 精品1| 精品一区二区免费观看| 亚洲精华国产精华液的使用体验| 春色校园在线视频观看| 亚洲天堂av无毛| 美女中出高潮动态图| 毛片一级片免费看久久久久| 男人舔奶头视频| 国产高清三级在线| 久久6这里有精品| 日韩av在线免费看完整版不卡| 高清在线视频一区二区三区| 综合色丁香网| 22中文网久久字幕| 亚洲av在线观看美女高潮| 成人影院久久| 免费观看av网站的网址| 欧美 日韩 精品 国产| 国产亚洲午夜精品一区二区久久| av天堂中文字幕网| 在线观看www视频免费| 国产亚洲91精品色在线| 亚洲经典国产精华液单| 欧美丝袜亚洲另类| 精品久久久久久久久av| 美女脱内裤让男人舔精品视频| 妹子高潮喷水视频| 国产亚洲午夜精品一区二区久久| 人体艺术视频欧美日本| 特大巨黑吊av在线直播| 国产亚洲最大av| 久久人妻熟女aⅴ| 国产爽快片一区二区三区| 久久ye,这里只有精品| 国产成人精品久久久久久| 亚洲美女搞黄在线观看| 久久影院123| 永久网站在线| 丝瓜视频免费看黄片| 免费av中文字幕在线| 熟女电影av网| 中文字幕亚洲精品专区| 亚洲国产精品一区三区| 中文天堂在线官网| 久久久国产欧美日韩av|