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

    基于無人機遙感的黃海綠潮擱淺生物量估算

    2021-11-18 09:30:00尚偉濤高志強姜曉鵬田信鵬郭少方
    海洋科學 2021年10期
    關(guān)鍵詞:乳山潮間帶文登

    尚偉濤, 高志強, 姜曉鵬, 田信鵬, 郭少方

    基于無人機遙感的黃海綠潮擱淺生物量估算

    尚偉濤1, 2, 3, 高志強1, 2, 姜曉鵬1, 2, 3, 田信鵬1, 2, 郭少方4

    (1. 中國科學院 海岸帶環(huán)境過程與生態(tài)修復(fù)重點實驗室, 中國科學院煙臺海岸帶研究所, 山東 煙臺 264003; 2. 山東省海岸帶環(huán)境過程重點實驗室, 中國科學院 煙臺海岸帶研究所, 山東 煙臺 264003; 3. 中國科學院大學, 北京 100049; 4. 煙臺市科技情報研究所, 山東 煙臺 264003)

    滸苔在近岸擱淺后會破壞海岸景觀, 干擾水上運動, 給濱海旅游業(yè)造成嚴重影響。本文使用無人機搭載的多光譜和可見光傳感器對山東半島的海陽、乳山和文登的三個海灘擱淺的滸苔進行航拍監(jiān)測, 并結(jié)合地物光譜測量數(shù)據(jù), 分別選擇歸一化植被指數(shù)(NDVI)、差值植被指數(shù)(DVI)和虛擬基線高度浮藻指數(shù)(VB-FAH)對海灘擱淺滸苔與岸邊植被及非植被(海水、沙灘)進行識別評估, 并分別估算了三個研究區(qū)擱淺滸苔的生物量。研究結(jié)果表明: NDVI可以識別植被和非植被, 但無法區(qū)分潮間帶上部和潮間帶下部分布的滸苔; DVI和VB-FAH對植被和非植被的區(qū)分度不高, 但對不同分布的擱淺滸苔具有一定的區(qū)分度, 其中, DVI對潮間帶上部和潮間帶下部分布滸苔的識別能力優(yōu)于VB-FAH。因此, 通過對岸邊植被進行腌膜, 利用DVI構(gòu)建海灘擱淺滸苔生物量估算模型, 實現(xiàn)了海灘擱淺滸苔生物量的估算。海陽、乳山和文登三個海灘擱淺滸苔的生物量分別為1 468 t、745 t和5 034 t, 本文提出的方法可以為擱淺滸苔的清理和資源合理分配提供技術(shù)支持。

    無人機遙感; 綠潮; 擱淺生物量; 多光譜

    自2007年以來, 我國黃海海域每年都會暴發(fā)以滸苔為優(yōu)勢種的綠潮災(zāi)害, 已經(jīng)持續(xù)了13 a[1]。滸苔本身是一種無毒無害的大型綠藻, 但是大量綠藻在近岸堆積、腐爛后會對沿海的自然景觀造成嚴重影響, 尤其是對青島、煙臺、威海的濱海旅游業(yè)產(chǎn)生重大影響, 當?shù)卣磕甓家ㄙM大量的人力財力來清理堆積在海灘上的滸苔, 據(jù)估算, 2008年用于治理綠潮災(zāi)害的總費約為20億元[2]。

    衛(wèi)星遙感技術(shù)由于觀測尺度大、時間分辨率高、覆蓋范圍廣等特點, 在綠潮災(zāi)害的監(jiān)測中發(fā)揮了重要作用[3-4]。已有學者通過使用MODIS、GOCI、Landsat、GF-1 WFV、HJ-1A/1B衛(wèi)星遙感數(shù)據(jù)對綠潮的起源、時空變化過程進行了監(jiān)測研究[5-8]。但對于濱海景區(qū)這樣的小范圍綠潮重點防控區(qū)域, 衛(wèi)星遙感數(shù)據(jù)的時間和空間分辨率已經(jīng)不能滿足實時監(jiān)測和打撈的需求。

    無人機遙感技術(shù)因其體積小、靈活機動、運行成本低、實時傳輸?shù)忍攸c, 在精準農(nóng)業(yè)、森林資源調(diào)查、濕地植被監(jiān)測等方面得到了廣泛的應(yīng)用[9-11]。在綠潮研究方面, 已有學者使用無人機航拍輔助船只進行漂浮綠潮的打撈[12], 使用無人機結(jié)合衛(wèi)星遙感數(shù)據(jù)估算了黃海綠潮初始生物量[13], 也有學者使用不同的植被指數(shù)算法對岸灘滸苔和岸邊植被進行了區(qū)分研究, 在此基礎(chǔ)之上估算了岸灘滸苔的生物量[14],但僅研究了一個岸灘的實例, 并未對多個海灘的滸苔分布和生物量情況進行對比研究。

    本文通過無人機獲取的多光譜和可見光數(shù)據(jù)對山東半島的海陽、乳山、文登的三個海灘的擱淺滸苔進行了航拍監(jiān)測, 并結(jié)合現(xiàn)場獲取的地物光譜數(shù)據(jù), 選擇了歸一化植被指數(shù)(NDVI)、差值植被指數(shù)(DVI)和虛擬基線高度浮藻指數(shù)(VB-FAH)三種指數(shù)算法對海灘擱淺滸苔與岸邊植被、海水、沙灘的區(qū)分能力進行了對比研究, 在此基礎(chǔ)之上, 結(jié)合現(xiàn)場采樣數(shù)據(jù), 對三個海灘的滸苔的生物量進行了估算, 為更好地清理擱淺滸苔和合理分配資源提供技術(shù)支持。

    1 材料與方法

    1.1 研究區(qū)概況

    本研究選取了山東半島三個典型的海灘, 分別位于海陽市愛琴海海灘(121°9′48″E, 36°40′12″N)、乳山市銀灘海灘(121°46′40″E, 36°53′47″N)和文登市南海公園海灘(121°53′47″E, 36°55′48″N)(圖1), 這三個海灘都是重要的濱海旅游景區(qū), 自2007年綠潮災(zāi)害暴發(fā)以來, 這三個濱海旅游景區(qū)的海灘每年都會有大量的滸苔擱淺, 如果不及時清理, 堆積的滸苔會腐爛發(fā)臭, 污染周邊環(huán)境, 對濱海景區(qū)的旅游業(yè)造成影響。當?shù)卣磕甓家ㄙM大量的人力和物理來清理堆積在海灘上的滸苔。

    1.2 無人機數(shù)據(jù)獲取與處理

    2019年7月4日—6日使用大疆M600 Pro無人機分別對海陽市愛琴海海灘、乳山市銀灘海灘和文登市南海公園海灘進行正射航拍, 無人機飛行高度為200 m, 選擇晴朗或少云的天氣進行飛行作業(yè)。無人機平臺搭載一個Parrot Sequoia 多光譜相機, 它由兩個傳感器組成: 多光譜傳感器, 包括五個同步波段-四個多光譜波段(Green、Red、 Red edge、NIR)和一個RGB相機, 還有一個安裝在無人機頂部的光照傳感器。在飛行過程中, 光照傳感器獲取和記錄當前的光照環(huán)境數(shù)據(jù), 用于后續(xù)的數(shù)據(jù)處理, 多光譜相機的詳細參數(shù)見表1。多光譜相機拍攝獲取的多光譜數(shù)據(jù)和RGB數(shù)據(jù)帶有GPS模塊獲取的定位信息, 可以直接使用Pix4Dmapper軟件對數(shù)據(jù)進行處理, 生成正射影像(包括多光譜影像和RGB影像), 然后使用ENVI 5.3軟件和ArcMap 10.1軟件對數(shù)據(jù)進行分析和制圖。

    圖1 研究區(qū)概況

    表1 Parrot Sequoia多光譜相機主要參數(shù)

    “—”: 表格內(nèi)容為空

    1.3 光譜測量

    為了研究滸苔和周邊地物的光譜特征, 使用USB4000(Ocean Optics Inc.)地物光譜儀對海灘擱淺的滸苔、岸邊植被、海水、沙灘的光譜進行測量, 地物光譜儀波長范圍為345.30~1 024.12 nm, 光譜分辨率為0.2 nm。測量時, 選擇晴朗的天氣, 測量人員穿黑色衣服, 地物光譜儀探頭距離被測量物體50 cm, 首先測量標準參考板的輻射亮度值, 然后再測量地物的輻射亮度值, 每種地物重復(fù)測量3次取平均值, 反射率計算公式為:

    式中,()表示被測量地物的反射率,()表示被測量地物的輻射亮度值,L()表示標準參考板的輻射亮度值,()表示標準參考板的反射率, 本次測量中使用的標準參考板的反射率是0.25。

    1.4 海灘擱淺滸苔提取方法

    為了研究不同的指數(shù)方法對海灘擱淺滸苔的提取精度, 本文選取了歸一化植被指數(shù)(NDVI)[15]、差值植被指數(shù)(DVI)[16]和虛擬基線高度浮藻指數(shù)(VB- FAH)[7]三個指數(shù)對海灘擱淺滸苔進行了提取研究, 計算公式如下:

    NDVI=(NIR–R)/(NIR+R), (2)

    DVI=NIR–R, (3)

    公式中,G、R和NIR分別表示無人機多光譜影像的綠波段、紅波段和近紅外波段。G、R、NIR是這三個波段對應(yīng)的中心波長, 分別為550 nm、660 nm和790 nm。

    為了比較NDVI、DVI和VB-FAH這三種指數(shù)算法對海灘擱淺滸苔的識別能力, 借鑒模式識別領(lǐng)域用來區(qū)分兩類樣本的“類間距”概念[17]來判斷不同算法對滸苔和三種地物(岸邊植被、海水和沙灘)的區(qū)分能力:

    1.5 海灘擱淺滸苔生物量估算方法

    1.5.1 海灘擱淺滸苔厚度和生物量數(shù)據(jù)采集

    由于潮汐作用的影響, 海灘擱淺滸苔在潮間帶上部和潮間帶下部的分布有所差異, 為了更好地研究擱淺滸苔的分布情況, 將研究區(qū)分為潮間帶上部采樣區(qū)和潮間帶下部采樣區(qū), 每個采樣區(qū)平行間隔100 m左右設(shè)置一個采樣點, 海陽、乳山和文登三個研究區(qū)分別設(shè)置了16個、10個和18個采樣點(由于乳山研究區(qū)海灘擱淺滸苔分布不均勻, 所以根據(jù)實際情況只設(shè)置了10個采樣點)。對于每個采樣點選取1 m×1 m范圍進行采樣, 首先使用精度為0.1 cm的量尺測量正方形樣點的四個角和中心的滸苔的厚度, 取平均值作為該采樣點的平均厚度。然后使用精度為0.01 kg的電子秤對樣點范圍內(nèi)的滸苔重量進行測量, 測量3次取平均值, 再計算單位體積滸苔的重量(單位為kg/m3)作為該采樣點滸苔的密度。

    圖2 研究區(qū)采樣點

    1.5.2 海灘擱淺滸苔生物量計算方法

    為了計算海灘擱淺滸苔的生物量, 首先使用DVI指數(shù)算法提取每個研究區(qū)擱淺滸苔在潮間帶上部和潮間帶下部的分布范圍, 然后按照擱淺滸苔在無人機RGB影像上的分布紋理和采樣點的分布位置將潮間帶上部和潮間帶下部提取的擱淺滸苔的分布范圍劃分為若干個子研究區(qū), 盡量使每個子研究區(qū)都有采樣點分布, 然后使用每個采樣點的滸苔分布密度和厚度數(shù)據(jù)作為采樣點所在子研究區(qū)的滸苔平均密度和平均厚度。使用公式(6)分別計算每個子研究區(qū)的擱淺滸苔生物量, 然后統(tǒng)計計算研究區(qū)滸苔的總生物量。

    其中,表示研究區(qū)滸苔的總生物量,A表示子研究區(qū)海灘擱淺滸苔的分布面積,D表示子研究區(qū)海灘擱淺滸苔的密度,T表示子研究區(qū)海灘擱淺滸苔的的厚度。

    2 結(jié)果

    2.1 海灘不同地物光譜特征

    使用地物光譜儀對海灘滸苔、岸邊植被、海水和沙灘進行地物光譜測量并繪制了不同地物的光譜曲線, 如圖3a所示。海灘滸苔和岸邊植被的光譜曲線很相似, 在可見光波段, 有一個小的反射峰, 位置在550 nm處(綠波段), 而在450 nm(藍波段)和670 nm(紅波段)附近則是兩個吸收帶, 這是由于葉綠素對藍光和紅光有較強的吸收作用, 對綠光有較強的反射作用引起的。在近紅外波段(700~800 nm), 由于反射率快速增長, 形成了植被特有的紅邊特征, 這是植被區(qū)別于其他地物的主要特征。海水在藍波段有較強的反射率, 而在其他波段則有很強的吸收率, 尤其是近紅外波段, 吸收率更強, 使得近紅外遙感圖上的海水呈現(xiàn)黑色。沙灘的光譜曲線從藍波段到近紅外波段反射率呈現(xiàn)逐漸增高的趨勢, 沒有出現(xiàn)明顯的反射峰特征。

    圖3 不同地物光譜曲線

    通過對無人機搭載的多光譜相機獲取的影像中不同地物的反射率進行統(tǒng)計, 如圖3b所示, 該多光譜相機包括了綠波段(550 nm)、紅波段(660 nm)、紅邊波段(735 nm)和近紅外波段(790 nm)四個波段, 滸苔由于所含水分的不同, 在這四個波段上的反射率都比岸邊植被的高, 在紅波段, 滸苔的反射率為0.14, 岸邊植被的反射率為0.08, 在近紅外波段, 滸苔的反射率為0.5, 岸邊植被的反射率為0.3。海水在紅邊波段和近紅外波段的反射率都比滸苔的低, 而在紅波段的反射率高于滸苔的反射率。沙灘在四個波段的反射率差異不大, 都在0.2以下, 在紅波段, 沙灘的反射率高于滸苔, 而在紅邊波段和近紅外波段, 滸苔的反射率則高于沙灘。

    通過對以上四種地物的光譜曲線的分析, 滸苔和植被的光譜特征比較明顯, 有一個反射峰和紅邊特征, 而海水和沙灘則沒有這樣的光譜特征, 利用這一光譜特征的差異, 通過波段之間的差值運算, 選擇差值比較大的植被指數(shù)(如NDVI、DVI、VB-FAH)就可以將這四種地物進行區(qū)分, 以提高這四種地物的提取精度, 為后續(xù)的研究提供精確的數(shù)據(jù)。

    2.2 不同指數(shù)對比研究

    圖4顯示了海陽、乳山、文登三個研究區(qū)的RGB、NDVI、DVI和VB-FAH四種圖像的對比圖。其中NDVI指數(shù)的取值范圍分別是–0.62至0.85、–0.54至0.88、–0.81至0.93, 乳山研究區(qū)NDVI的取值范圍大于海陽和乳山研究區(qū), DVI指數(shù)的取值范圍分別是–0.33至0.50、–0.45至0.70、–1.04至1.31, 文登研究區(qū)DVI的取值范圍大于海陽和乳山研究區(qū), VB-FAH指數(shù)的取值范圍分別是–0.25至0.47、–0.19至0.73、–0.59至1.31, 文登研究區(qū)VB-FAH的取值范圍大于海陽和乳山研究區(qū)。滸苔在三種指數(shù)圖像上都有很高的值, 在NDVI指數(shù)圖像上, 岸邊植被也有很高的值, 和滸苔的區(qū)分度不高, 而在DVI和VB-FAH圖像上, 只有很小一部分植被和滸苔的值接近。

    圖4 海陽、乳山和文登三個研究區(qū)的RGB、NDVI、DVI和VB-FAH四種圖像的對比圖

    為了進一步研究這三種指數(shù)的差異性, 從NDVI、DVI、VB-FAH三種指數(shù)圖像上分別選取100個四種地物(岸邊植被、海水、滸苔、沙灘)的樣本, 統(tǒng)計四種地物的直方圖分布情況, 結(jié)果如圖5所示。從圖中可以看到, NDVI指數(shù)對滸苔和岸邊植被的區(qū)分度不高, 但是滸苔與海水和沙灘有明顯的區(qū)分界限。DVI和VB-FAH指數(shù)對滸苔和岸邊植被有一定的區(qū)分度, 但是有一小部分滸苔和岸邊植被的取值相同, 很難通過單一閾值將二者進行區(qū)分。岸邊植被和沙灘在DVI圖像上的區(qū)分度很高, 而VB-FAH圖像上岸邊植被和沙灘有一部分的取值相同, 也很難通過單一閾值將二者進行區(qū)分。

    滸苔在海灘擱淺后, 滸苔所含的水分會隨著時間推移而逐漸減少, 先擱淺的滸苔比后擱淺的滸苔所含的水分少, 由于所含水分的不同, 在不同的指數(shù)圖像上, 不同時間段擱淺的滸苔會有所差異, 在DVI和VB-FAH指數(shù)圖像上, 先擱淺的滸苔比后擱淺的滸苔具有更高的值, 而在NDVI指數(shù)圖像上, 這一差異并不明顯。

    從圖4和圖5三種指數(shù)的比較中可以看到, NDVI可以識別植被和非植被, 但無法區(qū)分潮間帶上部和潮間帶下部分布的滸苔; DVI和VB-FAH對植被和非植被的區(qū)分度不高, 但對不同分布的擱淺滸苔具有一定的區(qū)分度。

    2.3 滸苔與岸邊植被、海水、沙灘樣本之間的距離比較

    為了進一步研究NDVI、DVI、VB-FAH這三個指數(shù)中滸苔與岸邊植被、海水、沙灘樣本之間的區(qū)分度, 分別計算了滸苔與岸邊植被、滸苔與海水、滸苔與沙灘樣本之間的距離, 計算結(jié)果如表2所示。在三個指數(shù)中, 滸苔和岸邊植被樣本之間的距離最大的是DVI指數(shù), 為0.46, NDVI指數(shù)的距離值最小, 為0.06。滸苔和海水樣本之間的距離最大的是NDVI指數(shù), 為2.42; 最小的是VB-FAH指數(shù), 為1.56。滸苔和沙灘樣本之間的距離最大的是NDVI指數(shù), 為2.26, 最小的是VB-FAH指數(shù), 為1.29。

    表2 滸苔與岸邊植被、海水、沙灘樣本之間的距離比較結(jié)果

    注:gv: 滸苔與岸邊植被樣本之間的距離;gw: 滸苔與海水樣本之間的距離;gs: 滸苔與沙灘樣本之間的距離

    從以上的結(jié)果中可以看出, NDVI指數(shù)可以很好地把滸苔與海水和沙灘區(qū)分開來, 但是由于岸邊植被和滸苔有很相似的光譜曲線, 很難將滸苔和岸邊植被區(qū)分開來。在區(qū)分潮間帶上部和潮間帶下部分布的滸苔時, DVI優(yōu)于VB-FAH。

    由于這三個指數(shù)都不能很好的將滸苔和岸邊植被區(qū)分開來, 在實際的滸苔提取過程中, 為了減少岸邊植被對滸苔提取精度的影響, 在滸苔和岸邊植被區(qū)分明顯的海灘過度地帶對岸邊植被做一個掩膜, 以消除岸邊植被對滸苔提取精度的影響。

    為了區(qū)分潮間帶上部和潮間帶下部滸苔的分布范圍, 綜合考慮這三個指數(shù)的特點, 選擇DVI指數(shù)作為滸苔提取的指數(shù)算法來分別提取潮間帶上部和潮間帶下部滸苔的分布面積。

    2.4 海灘擱淺生物量估算

    使用DVI指數(shù)按照不同的閾值分別提取了海陽、乳山、文登三個研究區(qū)的潮間帶上部和潮間帶下部的滸苔的分布面積, 結(jié)果如表3所示。海陽潮上帶和潮下帶滸苔的分布面積分別為41 785.27 m2和26 643.97 m2。乳山潮間帶上部和潮間帶下部滸苔的分布面積分別為48 148.76 m2和5 529.75 m2。文登潮間帶上部和潮間帶下部滸苔的分布面積分別為46 732.13 m2和40 492.97 m2。在三個研究區(qū),潮間帶上部滸苔的厚度都高于潮間帶下部滸苔的厚度。滸苔在三個研究區(qū)的密度分布也有所不同, 整體表現(xiàn)為潮間帶上部的滸苔密度大于潮間帶下部的滸苔密度, 滸苔密度最大的是乳山潮間帶上部, 而滸苔密度最小的是文登潮下帶。

    通過對每個研究區(qū)劃分的子研究區(qū)擱淺滸苔生物量的統(tǒng)計得到三個研究區(qū)擱淺滸苔的總生物量。海陽、乳山和文登這三個研究區(qū)的滸苔總生物量分別是1 468 t、745 t和5 034 t, 其中文登海灘的滸苔生物量最大, 而乳山海灘的滸苔生物量最小。

    3 討論

    3.1 三種指數(shù)提取海灘擱淺滸苔結(jié)果比較

    NDVI、DVI、VB-FAH這三個指數(shù)都是提取漂浮滸苔常用的指數(shù)算法, 在提取海灘擱淺滸苔方面的研究還未有報道。在過去的研究中, 這三個指數(shù)主要用于衛(wèi)星遙感數(shù)據(jù), 而在搭載多光譜傳感器的無人機遙感方面的應(yīng)用較少。已有的研究結(jié)果表明[7], VB-FAH提取漂浮滸苔的效果優(yōu)于NDVI, 而與DVI的提取效果相當, 而本文通過對這三個指數(shù)的比較發(fā)現(xiàn), DVI在區(qū)分潮間帶上部和潮間帶下部海灘擱淺滸苔方面優(yōu)于NDVI和VB-FAH, 但在區(qū)分植被和非植被方面, NDVI則優(yōu)于DVI和VB-FAH, 這可能是由于漂浮滸苔的背景環(huán)境是海水而擱淺滸苔的背景環(huán)境是海灘造成的。

    表3 三個研究區(qū)海灘擱淺滸苔生物量估算結(jié)果

    3.2 海灘擱淺生物量結(jié)果分析

    通過對海陽、乳山、文登三個典型的海灘的擱淺滸苔生物量進行計算, 滸苔總生物量分別為1 468 t、745 t和5 034 t。這三個海灘擱淺滸苔的生物量差距比較明顯, 文登海灘的擱淺滸苔生物量最多, 分別是海陽和乳山擱淺滸苔生物量的3.4倍和6.8倍, 如果不及時清理, 在很短的時間之內(nèi)大量滸苔就會腐爛, 發(fā)出惡臭的氣味, 嚴重影響濱海景觀, 給當?shù)氐穆糜螛I(yè)造成沖擊。當?shù)卣磕甓紩M織大量的人力和機械設(shè)備來清除海灘上的擱淺滸苔, 如果以載重為10 t的卡車來運輸這些滸苔, 清理完海陽、乳山、文登這三個研究區(qū)的海灘擱淺滸苔需要的車次分別為147、75和504。為了更高效地清理海灘上的擱淺滸苔, 需要根據(jù)滸苔的總生物量合理分配人力和機械設(shè)備。本文提出的海灘擱淺生物量估算方法, 可以為滸苔的清理工作提供數(shù)據(jù)支持, 方便決策部門合理配置資源。

    4 結(jié)論

    本文通過無人機獲取的多光譜和可見光數(shù)據(jù)對山東半島海陽、乳山、文登的三個海灘的擱淺滸苔進行了航拍監(jiān)測, 并結(jié)合現(xiàn)場獲取的地物光譜數(shù)據(jù), 選擇了NDVI、DVI和VB-FAH三種指數(shù)算法對海灘擱淺滸苔與岸邊植被、海水、沙灘的區(qū)分能力進行了對比研究, 在此基礎(chǔ)之上, 結(jié)合現(xiàn)場采樣數(shù)據(jù), 對三個海灘的滸苔的生物量進行了估算, 主要結(jié)論如下:

    (1) 海灘擱淺滸苔和植被的光譜特征比較相似, 在可見光波段有一個反射峰, 在近紅外波段有紅邊特征, 使用這一特性, 選擇差值植被指數(shù)(如NDVI、DVI、VB-FAH)就可以將海灘擱淺滸苔和非植被地物進行區(qū)分。

    (2) NDVI可以識別植被和非植被, 但無法區(qū)分潮間帶上部和潮間帶下部分布的滸苔; DVI和VB- FAH對植被和非植被的區(qū)分度不高, 但對不同分布的擱淺滸苔具有一定的區(qū)分度, 其中, DVI對潮間帶上部和潮間帶下部分布滸苔的識別能力優(yōu)于VB-FAH。

    (3) 海陽、乳山、文登三個研究區(qū)海灘擱淺滸苔的分布特征主要是潮間帶上部擱淺滸苔的厚度大于潮下帶擱淺滸苔的厚度, 潮間帶上部擱淺滸苔的密度大于潮下帶擱淺滸苔的密度, 滸苔密度最大的是乳山潮間帶上部, 而滸苔密度最小的是文登潮間帶下部。

    (4) 海陽、乳山和文登這三個研究區(qū)的滸苔總生物量分別是1 468 t、745 t和5 034 t, 其中文登海灘擱淺滸苔生物量最大, 而乳山海灘擱淺滸苔的生物量最小。

    使用無人機遙感技術(shù)結(jié)合實地采樣數(shù)據(jù), 可以快速地估算海灘擱淺滸苔的生物量, 為擱淺滸苔的清理和資源合理分配提供技術(shù)支持。

    [1] 于仁成, 孫松, 顏天, 等. 黃海綠潮研究: 回顧與展望[J]. 海洋與湖沼, 2018, 49(5): 942-949.

    YU Rencheng, SUN Song, YAN Tian, et al. Progresses and perspectives on green tide studies in the Yellow Sea[J]. Oceanologia et Limnologia Sinica, 2018, 49(5): 942-949.

    [2] YE Naihao, ZHANG Xiaowen, MAO Yuze, et al. ‘Green tides’ are overwhelming the coastline of our blue planet: taking the world’s largest example[J]. Ecological Research, 2011, 26(3): 477-485.

    [3] HU Chuanmin, HE Mingxia. Origin and offshore extent of floating algae in Olympic sailing area[J]. Eos, Tran-sactions American Geophysical Union, 2008, 89(33): 302-303.

    [4] 蔣興偉, 劉建強, 鄒斌, 等. 滸苔災(zāi)害衛(wèi)星遙感應(yīng)急監(jiān)視監(jiān)測系統(tǒng)及其應(yīng)用[J]. 海洋學報, 2009, 31(1): 52-64.

    JIANG Xingwei, LIU Jianqiang, ZOU Bin, et al. The satellite remote sensing system used in emergency response monitoring for Entermorpha prolifera disaster and its application[J]. Acta Oceanologica Sinica, 2009, 31(1): 52-64.

    [5] HU Chuanmin, LI Daqiu, CHEN Changsheng, et al. On the recurrent Ulva prolifera blooms in the Yellow Sea and East China Sea[J]. Journal of Geophysical Research-Oceans, 2010, 115: C05017.

    [6] SON Y B, MIN J-E, RYU J-H. Detecting massive green algae (Ulva prolifera) blooms in the Yellow Sea and East China Sea using Geostationary Ocean Color Ima-ger (GOCI) data[J]. Ocean Science Journal, 2012, 47(3): 359-375.

    [7] XING Qianguo, HU Chuanmin. Mapping macroalgal blooms in the Yellow Sea and East China Sea using HJ-1 and Landsat data: Application of a virtual baseline reflectance height technique[J]. Remote Sensing of Environment, 2016, 178: 113-126.

    [8] ZHANG Jianheng, SHI Jinting, GAO Song, et al. Annual patterns of macroalgal blooms in the Yellow Sea during 2007-2017[J]. Plos One, 2019, 14(1): e0210460.

    [9] DENG Lei, MAO Zhihui, LI Xiaojuan, et al. UAV- based multispectral remote sensing for precision agriculture: A comparison between different cameras[J]. Isprs Journal of Photogrammetry and Remote Sensing, 2018, 146: 124-136.

    [10] LIN Jiayuan, WANG Meimei, MA Mingguo, et al. Abo-veground tree biomass estimation of sparse subal-pine coniferous forest with UAV oblique photography[J]. Remote Sensing, 2018, 10(11): 1849.

    [11] ZHOU Zaiming, YANG Yanming, CHEN Benqing. Estimating spartina alterniflora fractional vegetation cover and aboveground biomass in a coastal wetland using SPOT6 satellite and UAV data[J]. Aquatic Botany, 2018, 144: 38-45.

    [12] XU Fuxiang, GAO Zhiqiang, SHANG Weitao, et al. Validation of MODIS-based monitoring for a green tide in the Yellow Sea with the aid of unmanned aerial vehicle[J]. Journal of Applied Remote Sensing, 2017, 11(1): 012007.

    [13] XU Fuxiang, GAO Zhiqiang, JIANG Xiaopeng, et al. A UAV and S2A data-based estimation of the initial biomass of green algae in the South Yellow Sea[J]. Marine Pollution Bulletin, 2018, 128: 408-414.

    [14] 李冬雪, 高志強, 尚偉濤, 等. 基于無人機多光譜影像的海濱景區(qū)滸苔信息提取研究[J]. 海洋環(huán)境科學, 2020, 39(3): 438-446.

    LI Dongxue, GAO Zhiqiang, SHANG Weitao, et al. Information extraction of Ulva Prolifera from coastal landscape using UAV multispectral remote sensing images[J]. Marine Environmental Science, 2020, 39(3): 438-446.

    [15] ROUSEJW, HAAS R H, SCHELL J A, et al. Monitoring Vegetation Systems in the Great Plains with ERTS[C]//Proceedings of the Third ERTS Symposium, Washington DC, 10-14 December 1973, 309-317.

    [16] RICHARDSON A J, EVERITT J H. Using spectral vegetation indices to estimate rangeland productivity[J]. Geocarto International, 1992, 7(1): 63-69.

    [17] WEBB A. 統(tǒng)計模式識別(第二版)[M]. 北京: 電子工業(yè)出版社, 2004.

    WEBB A. Statistical pattern recognition (Second edition)[M]. Beijing: Electronic Industry Press, 2004.

    Estimation of green tide stranded biomass in the Yellow Sea based on unmanned aerial vehicle remote sensing

    SHANG Wei-tao1, 2, 3, GAO Zhi-qiang1, 2, JIANG Xiao-peng1, 2, 3, TIAN Xin-peng1, 2, GUO Shao-fang4

    (1.CAS Key Laboratory of Coastal Environmental Processes and Ecological Remediation, Yantai Institute of Coastal Zone Research, Chinese Academy of Sciences, Yantai 264003, China; 2. Shandong Key Laboratory of Coastal Environmental Processes, Yantai Institute of Coastal Zone Research, Chinese Academy of Sciences, Yantai 264003, China; 3. University of Chinese Academy of Sciences, Beijing 100049, China; 4. Yantai Institute of S&T Information, Yantai 264003, China)

    UAV remote sensing; green tide; stranded biomass; multispectral

    stranded near a shore damages the coastal landscape, interferes with water sports, and seriously affects coastal tourism. In this study, the multispectral and visible light sensors carried using an unmanned aerial vehicle (UAV) were used to monitor the strandedon the beaches of Haiyang, Rushan, and Wendeng in the Shandong Peninsula. The normalized difference vegetation index (NDVI), difference vegetation index (DVI), and virtual baseline height floating algae index (VB-FAH) were selected to identify and evaluate the beach-strandedand shore vegetation and nonvegetation (sea water and sand), respectively, and the biomass of the strandedin three research areas was estimated. The results showed that the NDVI can identify vegetation and nonvegetation but it cannot identify thedistributed in the upper and lower tidal zones. DVI and VB-FAH cannot distinguish between vegetation and nonvegetation but they can identifydistributed in the upper and lower tidal zones. In these zones, DVI was superior to the VB-FAH in identifying the distribution of. Therefore, the biomass estimation model of the beach-strandedwas established by masking shore vegetation; and using DVI, the biomass of the beach-strandedwas estimated. The biomass of the strandedwas 1468, 745, and 5034 t in the beaches of Haiyang, Rushan, and Wendeng. The method proposed in this study can provide technical support for the cleaning up of the strandedand rational use of resources.

    Jul. 8, 2020

    [NSFC fund project, No. 41876107; NSFC-Shandong joint fund project, No. U1706219; Key Deployment Project of Center for Ocean Mega-Science, Chinese Academy of Sciences, No. COMS2019J02; Key Research Program of Frontier Science, Chinese Academy of Sciences, No. ZDBS-LY-7010; CAS Key Laboratory of Marine Ecology and Environmental Sciences, Institute of Oceanology, Chinese Academy of Sciences No. KLMEES202005; The National Key R&D Program of China, No. 2019YFD0900705]

    P76

    A

    1000-3096(2021)10-0011-09

    10.11759/hykx20200708001

    2020-07-08;

    2020-08-19

    國家自然科學基金項目(41876107); 山東省聯(lián)合基金項目(U1706219); 中國科學院海洋大科學研究中心重點部署項目(COMS2019J02); 中國科學院前沿科學重點研究計劃(ZDBS-LY-7010); 中國科學院海洋生態(tài)與環(huán)境科學重點實驗室(中國科學院海洋研究所)開放基金資助(KLMEES202005); 國家重點研發(fā)計劃“藍色糧倉科技創(chuàng)新”項目(2019YFD0900705)

    尚偉濤(1987—), 男, 甘肅慶陽人, 博士研究生, 主要從事無人機遙感海岸帶災(zāi)害監(jiān)測及應(yīng)對研究, E-mail: wtshang@yic.ac.cn; 高志強(1966—), 通信作者, 研究員, 博士生導(dǎo)師, E-mail: zqgao@yic.ac.cn

    (本文編輯: 康亦兼)

    猜你喜歡
    乳山潮間帶文登
    走進潮間帶
    文化經(jīng)濟融合視域下乳山民間文學的傳承與發(fā)展研究
    赴文登,靜享“陌上花開”
    金橋(2022年5期)2022-08-24 01:41:28
    山東乳山:草莓小鎮(zhèn)的“莓”好生活
    基于海陸融合DEM的大尺度潮間帶提取方法
    乳山長牡蠣(Crassostrea gigas)的抗性基因表達和生存環(huán)境的季節(jié)差異
    海洋通報(2020年5期)2021-01-14 09:27:02
    文登西洋參“紅”了!
    金橋(2020年11期)2020-12-14 07:52:54
    福建閩江口潮間帶大型底棲動物次級生產(chǎn)力時空特征
    海洋通報(2020年3期)2020-10-20 09:26:36
    “參”得人心的文登西洋參
    金橋(2020年7期)2020-08-13 03:07:02
    文登:給我三天,還你千年
    金橋(2018年7期)2018-09-25 02:28:24
    69精品国产乱码久久久| 亚洲一区二区三区不卡视频| 国产成人av激情在线播放| 亚洲成国产人片在线观看| 香蕉国产在线看| 一区在线观看完整版| 国产色视频综合| 国产成人影院久久av| av网站免费在线观看视频| 免费高清在线观看日韩| 亚洲精品粉嫩美女一区| 宅男免费午夜| 久热这里只有精品99| 伊人久久大香线蕉亚洲五| 高清在线国产一区| 91麻豆av在线| 国产男女内射视频| 欧美久久黑人一区二区| 婷婷丁香在线五月| 在线播放国产精品三级| 午夜久久久在线观看| aaaaa片日本免费| 亚洲性夜色夜夜综合| 免费看a级黄色片| 午夜91福利影院| 久久午夜综合久久蜜桃| 无人区码免费观看不卡| 如日韩欧美国产精品一区二区三区| 如日韩欧美国产精品一区二区三区| 欧美日韩成人在线一区二区| 高清毛片免费观看视频网站 | 亚洲 欧美一区二区三区| 黑人欧美特级aaaaaa片| 亚洲欧美激情综合另类| 午夜老司机福利片| 女人被狂操c到高潮| 国产精品偷伦视频观看了| 99riav亚洲国产免费| 在线看a的网站| 国产精品一区二区在线观看99| 宅男免费午夜| 色在线成人网| 激情在线观看视频在线高清 | 男人舔女人的私密视频| 精品一区二区三区视频在线观看免费 | 国产精品免费视频内射| 中文亚洲av片在线观看爽 | 69av精品久久久久久| 免费女性裸体啪啪无遮挡网站| 天天操日日干夜夜撸| 国产精品 国内视频| 黄片大片在线免费观看| 别揉我奶头~嗯~啊~动态视频| 亚洲人成伊人成综合网2020| 色婷婷av一区二区三区视频| 青草久久国产| 国产一区二区三区综合在线观看| 自线自在国产av| 精品卡一卡二卡四卡免费| 国产精品免费大片| 后天国语完整版免费观看| 国产成人精品无人区| 丁香六月欧美| 伦理电影免费视频| 日韩有码中文字幕| 久久天躁狠狠躁夜夜2o2o| 午夜福利,免费看| 午夜福利欧美成人| 老熟妇仑乱视频hdxx| 怎么达到女性高潮| x7x7x7水蜜桃| 日本撒尿小便嘘嘘汇集6| 老司机午夜十八禁免费视频| 男人操女人黄网站| 99精品久久久久人妻精品| av网站在线播放免费| 国产单亲对白刺激| 少妇的丰满在线观看| 大片电影免费在线观看免费| 最新在线观看一区二区三区| 精品国产一区二区久久| 久久九九热精品免费| 高清毛片免费观看视频网站 | 啦啦啦在线免费观看视频4| 国产单亲对白刺激| 亚洲欧美激情综合另类| 国产极品粉嫩免费观看在线| 18禁黄网站禁片午夜丰满| 69精品国产乱码久久久| 99香蕉大伊视频| 国产精品成人在线| 精品熟女少妇八av免费久了| 高清毛片免费观看视频网站 | 黄色 视频免费看| 少妇 在线观看| 国产成人免费无遮挡视频| 自拍欧美九色日韩亚洲蝌蚪91| 高清毛片免费观看视频网站 | 日本五十路高清| 在线视频色国产色| 一区福利在线观看| 亚洲av成人av| 狠狠狠狠99中文字幕| 天天操日日干夜夜撸| 日韩人妻精品一区2区三区| 五月开心婷婷网| 国产欧美日韩一区二区三| 国产精品久久久av美女十八| 十八禁高潮呻吟视频| 三上悠亚av全集在线观看| 岛国在线观看网站| av电影中文网址| 窝窝影院91人妻| 成年女人毛片免费观看观看9 | av国产精品久久久久影院| 亚洲情色 制服丝袜| 欧美另类亚洲清纯唯美| 69精品国产乱码久久久| 国产精华一区二区三区| 在线观看免费日韩欧美大片| 国产熟女午夜一区二区三区| 国产精品久久视频播放| 亚洲人成77777在线视频| 18禁国产床啪视频网站| 精品乱码久久久久久99久播| 在线免费观看的www视频| 日韩欧美一区二区三区在线观看 | 大陆偷拍与自拍| 极品教师在线免费播放| 国产亚洲精品久久久久5区| 国产av精品麻豆| 精品国产一区二区三区久久久樱花| 成在线人永久免费视频| 国产一区二区三区在线臀色熟女 | 90打野战视频偷拍视频| 国产av精品麻豆| 国产伦人伦偷精品视频| 老鸭窝网址在线观看| 青草久久国产| 老司机午夜十八禁免费视频| 欧美一级毛片孕妇| 欧美中文综合在线视频| av天堂久久9| 亚洲在线自拍视频| 丁香六月欧美| 好看av亚洲va欧美ⅴa在| 国产成+人综合+亚洲专区| a级毛片黄视频| 久久精品亚洲精品国产色婷小说| 一本大道久久a久久精品| 搡老熟女国产l中国老女人| 亚洲一区二区三区不卡视频| 别揉我奶头~嗯~啊~动态视频| 久久久国产精品麻豆| 国产精品久久久人人做人人爽| 精品亚洲成国产av| 欧美日韩亚洲高清精品| 午夜视频精品福利| 纯流量卡能插随身wifi吗| 国产无遮挡羞羞视频在线观看| ponron亚洲| 丝袜在线中文字幕| 久久香蕉精品热| 午夜两性在线视频| 亚洲精品在线美女| 久久香蕉精品热| 亚洲成人免费电影在线观看| 91九色精品人成在线观看| 99精品在免费线老司机午夜| 9色porny在线观看| 黑人欧美特级aaaaaa片| 成人三级做爰电影| 国产成人av教育| 精品无人区乱码1区二区| 色婷婷av一区二区三区视频| 午夜老司机福利片| 色综合欧美亚洲国产小说| 热99re8久久精品国产| 在线观看午夜福利视频| 亚洲精品自拍成人| 男人的好看免费观看在线视频 | 国产精品99久久99久久久不卡| 中文亚洲av片在线观看爽 | a在线观看视频网站| 国产精品成人在线| 欧美日韩一级在线毛片| 国产成人av教育| 久久久国产成人精品二区 | 国产亚洲欧美精品永久| 国产片内射在线| av国产精品久久久久影院| 欧美乱色亚洲激情| 美女国产高潮福利片在线看| 亚洲欧美精品综合一区二区三区| 手机成人av网站| 激情在线观看视频在线高清 | 69av精品久久久久久| 久久九九热精品免费| 黄色片一级片一级黄色片| 亚洲成a人片在线一区二区| 亚洲综合色网址| 日韩欧美在线二视频 | 国产免费av片在线观看野外av| 亚洲男人天堂网一区| 天天添夜夜摸| av超薄肉色丝袜交足视频| 免费看十八禁软件| 69精品国产乱码久久久| 夜夜夜夜夜久久久久| 日本撒尿小便嘘嘘汇集6| 成人国语在线视频| 久久人妻av系列| 啦啦啦视频在线资源免费观看| 黑人操中国人逼视频| 久久国产精品影院| 自线自在国产av| 波多野结衣av一区二区av| 国产男靠女视频免费网站| 日本五十路高清| 亚洲国产看品久久| av视频免费观看在线观看| 亚洲熟妇中文字幕五十中出 | 久久久国产精品麻豆| 高清黄色对白视频在线免费看| 成年人黄色毛片网站| 日本一区二区免费在线视频| 女人精品久久久久毛片| 高清在线国产一区| 少妇猛男粗大的猛烈进出视频| 在线观看舔阴道视频| 午夜激情av网站| 18禁观看日本| 免费黄频网站在线观看国产| av欧美777| 黄色视频不卡| 高清黄色对白视频在线免费看| a在线观看视频网站| 很黄的视频免费| 成人18禁在线播放| 久久国产亚洲av麻豆专区| 欧美日韩瑟瑟在线播放| 不卡一级毛片| tube8黄色片| 中文字幕人妻丝袜一区二区| 亚洲黑人精品在线| 日本vs欧美在线观看视频| 真人做人爱边吃奶动态| 在线观看免费视频日本深夜| 乱人伦中国视频| 人妻一区二区av| 中亚洲国语对白在线视频| 女人被狂操c到高潮| 久久久久视频综合| 久久人妻av系列| 老司机亚洲免费影院| 悠悠久久av| 亚洲成人国产一区在线观看| 91av网站免费观看| 一级作爱视频免费观看| 国产男女内射视频| 久久精品国产综合久久久| 亚洲精品在线美女| 他把我摸到了高潮在线观看| 男女下面插进去视频免费观看| 人人妻人人澡人人看| 制服人妻中文乱码| 亚洲av片天天在线观看| 国产精品免费视频内射| 丝袜在线中文字幕| 国产精品久久久av美女十八| 午夜免费观看网址| 国产精品九九99| 91av网站免费观看| 一边摸一边抽搐一进一小说 | 国产乱人伦免费视频| 久久久精品区二区三区| 久久久久久久国产电影| 久久久水蜜桃国产精品网| 国产男女超爽视频在线观看| 亚洲av欧美aⅴ国产| 女警被强在线播放| 色尼玛亚洲综合影院| 最近最新中文字幕大全电影3 | 成人国语在线视频| 国产男女内射视频| 搡老熟女国产l中国老女人| 精品国产亚洲在线| av线在线观看网站| 一二三四社区在线视频社区8| 一级,二级,三级黄色视频| 色婷婷av一区二区三区视频| 黄色丝袜av网址大全| 岛国在线观看网站| 黑人巨大精品欧美一区二区蜜桃| 久久亚洲精品不卡| 国产精品一区二区在线观看99| 18禁裸乳无遮挡免费网站照片 | 国产一卡二卡三卡精品| 一个人免费在线观看的高清视频| 最近最新中文字幕大全免费视频| 国产色视频综合| 老司机午夜福利在线观看视频| 国产色视频综合| 成人三级做爰电影| 日日摸夜夜添夜夜添小说| 免费在线观看视频国产中文字幕亚洲| 女人高潮潮喷娇喘18禁视频| 两个人免费观看高清视频| 国产欧美日韩精品亚洲av| 亚洲性夜色夜夜综合| 国产精品成人在线| 国产精品免费大片| 亚洲一区二区三区不卡视频| 国产在线一区二区三区精| 大香蕉久久网| 精品熟女少妇八av免费久了| 一区二区三区国产精品乱码| 久久婷婷成人综合色麻豆| 精品人妻在线不人妻| www.999成人在线观看| 欧美精品人与动牲交sv欧美| 国产麻豆69| 国产成人啪精品午夜网站| 免费观看精品视频网站| 国产精品亚洲一级av第二区| tube8黄色片| 精品国产乱码久久久久久男人| a级毛片在线看网站| 精品熟女少妇八av免费久了| 老司机午夜福利在线观看视频| 日日爽夜夜爽网站| 在线十欧美十亚洲十日本专区| 涩涩av久久男人的天堂| 在线天堂中文资源库| 黄色视频,在线免费观看| 在线观看免费日韩欧美大片| 欧美一级毛片孕妇| 两性午夜刺激爽爽歪歪视频在线观看 | 极品少妇高潮喷水抽搐| 深夜精品福利| 亚洲色图 男人天堂 中文字幕| 人人澡人人妻人| 亚洲中文字幕日韩| 日韩 欧美 亚洲 中文字幕| 宅男免费午夜| 电影成人av| 午夜久久久在线观看| 女性被躁到高潮视频| 国产精品久久久av美女十八| 午夜亚洲福利在线播放| 久久这里只有精品19| 久久国产亚洲av麻豆专区| 91字幕亚洲| 大陆偷拍与自拍| 自拍欧美九色日韩亚洲蝌蚪91| 啦啦啦视频在线资源免费观看| 老熟妇仑乱视频hdxx| 国产成人精品在线电影| 国产精品秋霞免费鲁丝片| 咕卡用的链子| 欧美av亚洲av综合av国产av| av免费在线观看网站| 精品人妻1区二区| 一夜夜www| www.精华液| 亚洲精品在线美女| 国产亚洲欧美在线一区二区| 性少妇av在线| 一边摸一边抽搐一进一出视频| 久久精品国产亚洲av香蕉五月 | 国产在视频线精品| 丁香欧美五月| 一级片免费观看大全| 人人妻人人澡人人看| 日日夜夜操网爽| 每晚都被弄得嗷嗷叫到高潮| 一二三四社区在线视频社区8| 高清毛片免费观看视频网站 | 麻豆av在线久日| 高清在线国产一区| 国产精品1区2区在线观看. | 久久天躁狠狠躁夜夜2o2o| 久久精品国产99精品国产亚洲性色 | 夜夜夜夜夜久久久久| 窝窝影院91人妻| 久久精品国产a三级三级三级| 免费高清在线观看日韩| 777米奇影视久久| 成人特级黄色片久久久久久久| av天堂在线播放| av网站免费在线观看视频| 一级毛片精品| 午夜免费观看网址| 亚洲成国产人片在线观看| 午夜激情av网站| 亚洲成av片中文字幕在线观看| 久久久国产成人精品二区 | 久热爱精品视频在线9| 欧美老熟妇乱子伦牲交| av片东京热男人的天堂| 欧美成人午夜精品| 一夜夜www| 国产精品电影一区二区三区 | 在线观看免费日韩欧美大片| 制服人妻中文乱码| 久久久国产欧美日韩av| 亚洲人成电影免费在线| 国产精品秋霞免费鲁丝片| 国产aⅴ精品一区二区三区波| 国产成人欧美| 日韩熟女老妇一区二区性免费视频| 欧美丝袜亚洲另类 | 91老司机精品| 天天躁狠狠躁夜夜躁狠狠躁| 波多野结衣一区麻豆| 欧美精品亚洲一区二区| 视频区欧美日本亚洲| 久久久久久久久免费视频了| 一级a爱片免费观看的视频| 两性夫妻黄色片| 97人妻天天添夜夜摸| 日韩欧美免费精品| 日日摸夜夜添夜夜添小说| 国产一卡二卡三卡精品| 自拍欧美九色日韩亚洲蝌蚪91| 欧美+亚洲+日韩+国产| 桃红色精品国产亚洲av| 精品福利观看| 性色av乱码一区二区三区2| 亚洲中文日韩欧美视频| 久久中文看片网| 久久草成人影院| 欧美在线黄色| 国产日韩一区二区三区精品不卡| 欧美日本中文国产一区发布| 国产亚洲精品一区二区www | 男人舔女人的私密视频| 首页视频小说图片口味搜索| 午夜影院日韩av| 精品午夜福利视频在线观看一区| 午夜免费观看网址| 看片在线看免费视频| 亚洲熟妇熟女久久| 天天躁日日躁夜夜躁夜夜| cao死你这个sao货| 不卡一级毛片| 国产97色在线日韩免费| 少妇被粗大的猛进出69影院| 黄色怎么调成土黄色| 午夜福利在线免费观看网站| 老熟妇乱子伦视频在线观看| 国产黄色免费在线视频| 一进一出抽搐gif免费好疼 | 老熟妇乱子伦视频在线观看| 久久天躁狠狠躁夜夜2o2o| 中文亚洲av片在线观看爽 | 亚洲欧美色中文字幕在线| 国内久久婷婷六月综合欲色啪| www.999成人在线观看| 99精品久久久久人妻精品| 日本黄色视频三级网站网址 | 午夜福利乱码中文字幕| 亚洲精品国产色婷婷电影| www.999成人在线观看| 黑人操中国人逼视频| 亚洲七黄色美女视频| 国产精品一区二区精品视频观看| 午夜精品久久久久久毛片777| av线在线观看网站| 欧洲精品卡2卡3卡4卡5卡区| 天堂√8在线中文| 黑人操中国人逼视频| 免费人成视频x8x8入口观看| 免费在线观看影片大全网站| av片东京热男人的天堂| 丰满的人妻完整版| 精品高清国产在线一区| 欧美激情高清一区二区三区| 国产亚洲精品久久久久久毛片 | 亚洲第一av免费看| 欧美大码av| 亚洲精品粉嫩美女一区| 欧美日韩瑟瑟在线播放| 欧美色视频一区免费| 99国产极品粉嫩在线观看| 欧美中文综合在线视频| 亚洲熟妇中文字幕五十中出 | 亚洲第一av免费看| 国产成人av激情在线播放| 国产精品一区二区在线观看99| 亚洲午夜精品一区,二区,三区| 桃红色精品国产亚洲av| 涩涩av久久男人的天堂| 成人影院久久| 亚洲av日韩精品久久久久久密| 欧美日韩瑟瑟在线播放| 麻豆国产av国片精品| 新久久久久国产一级毛片| 国内久久婷婷六月综合欲色啪| 久久香蕉国产精品| 热99国产精品久久久久久7| 大香蕉久久成人网| 美女国产高潮福利片在线看| 制服诱惑二区| 搡老岳熟女国产| 最新美女视频免费是黄的| 亚洲人成77777在线视频| 老司机深夜福利视频在线观看| 亚洲人成电影观看| 涩涩av久久男人的天堂| 久久影院123| 亚洲性夜色夜夜综合| 少妇 在线观看| 热99re8久久精品国产| 国产主播在线观看一区二区| 亚洲精品av麻豆狂野| 黄频高清免费视频| 天堂俺去俺来也www色官网| 美女 人体艺术 gogo| 精品人妻在线不人妻| 免费一级毛片在线播放高清视频 | 99在线人妻在线中文字幕 | 欧美日韩视频精品一区| 在线观看一区二区三区激情| xxx96com| 99国产综合亚洲精品| 一级黄色大片毛片| av片东京热男人的天堂| 亚洲,欧美精品.| 搡老熟女国产l中国老女人| 国产成人影院久久av| 男女午夜视频在线观看| 久久人妻av系列| 欧美激情高清一区二区三区| 妹子高潮喷水视频| 成年女人毛片免费观看观看9 | 国产片内射在线| 日韩视频一区二区在线观看| 久久国产精品影院| 夜夜夜夜夜久久久久| 免费在线观看完整版高清| 午夜福利影视在线免费观看| 青草久久国产| 欧美在线黄色| 国产成人免费无遮挡视频| 欧美在线一区亚洲| 久久人人97超碰香蕉20202| 午夜久久久在线观看| 久久久国产精品麻豆| 国产精品免费视频内射| av片东京热男人的天堂| 美女高潮喷水抽搐中文字幕| 欧美精品啪啪一区二区三区| 日韩有码中文字幕| 国产精品国产av在线观看| 精品久久久精品久久久| 如日韩欧美国产精品一区二区三区| 亚洲少妇的诱惑av| 久久人人爽av亚洲精品天堂| 国内毛片毛片毛片毛片毛片| 身体一侧抽搐| 国产精品久久电影中文字幕 | 国产成人av激情在线播放| а√天堂www在线а√下载 | 777米奇影视久久| 在线天堂中文资源库| 成年人午夜在线观看视频| 一级作爱视频免费观看| av欧美777| 亚洲午夜理论影院| 国产黄色免费在线视频| 欧美av亚洲av综合av国产av| 国产99久久九九免费精品| 男女下面插进去视频免费观看| 999久久久国产精品视频| 熟女少妇亚洲综合色aaa.| 成人亚洲精品一区在线观看| 亚洲欧美激情在线| 777米奇影视久久| 欧美日韩中文字幕国产精品一区二区三区 | 人人妻人人澡人人爽人人夜夜| 亚洲欧美激情在线| 777米奇影视久久| 日韩熟女老妇一区二区性免费视频| 国产精品98久久久久久宅男小说| 亚洲成a人片在线一区二区| 在线永久观看黄色视频| 91九色精品人成在线观看| 国产av又大| 日韩 欧美 亚洲 中文字幕| 男男h啪啪无遮挡| 这个男人来自地球电影免费观看| 涩涩av久久男人的天堂| 久久青草综合色| 最近最新免费中文字幕在线| 无遮挡黄片免费观看| 丁香六月欧美| 视频在线观看一区二区三区| 午夜影院日韩av| 国精品久久久久久国模美| 最近最新免费中文字幕在线| 波多野结衣av一区二区av| 日日摸夜夜添夜夜添小说| 最新的欧美精品一区二区| 欧美在线黄色| 黑人操中国人逼视频| 成人黄色视频免费在线看| 国产免费av片在线观看野外av| 亚洲欧美一区二区三区久久| 中文欧美无线码| 久久精品aⅴ一区二区三区四区| 国产亚洲精品一区二区www | 国产亚洲一区二区精品| 免费在线观看黄色视频的| 老汉色∧v一级毛片|