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

    松材線蟲病變色松樹遙感監(jiān)測研究進(jìn)展

    2020-07-20 07:10:50李存軍蔣麗雅胡海棠
    林業(yè)科學(xué)研究 2020年3期
    關(guān)鍵詞:松材線蟲病變色

    陶 歡,李存軍*,程 成,蔣麗雅,胡海棠

    (1.北京農(nóng)業(yè)信息技術(shù)研究中心,北京 100097;2.安徽省林業(yè)有害生物防治檢疫局,安徽合肥 230031)

    松材線蟲?。╬ine wilt disease, PWD)被稱為松樹的“癌癥”,松樹感染此病后,從發(fā)病到死亡只要2~3 月,最快的40 d 左右即可表現(xiàn)出枯死狀(以下簡稱“變色松樹”),3~5 a 間便可造成大面積毀林的惡性災(zāi)害。松材線蟲源于北美,傳入日本后于1905 年在長崎被發(fā)現(xiàn),進(jìn)而逐漸擴散到東亞的中國、韓國和歐洲的葡萄牙等國家和地區(qū),但直到1971 年才被證實其病原是松材線蟲( Bursaphelenchus xylophilus, 或 Pine wood nematode),傳播的媒介是甲蟲[1]。目前國際上關(guān)于松材線蟲病的致病機理仍存在爭議,達(dá)成的共識是松材線蟲病是由多種病原體、寄主松樹、昆蟲媒介和氣候條件共同作用的結(jié)果[2]。松樹感染松材線蟲病后最顯著的外部特征變化是松針葉的變色和松脂分泌的減少,究其原因是松樹內(nèi)部光合生理參數(shù)(如光合作用、葉綠素含量、含水量)的變化[3]。這些變化在遙感圖像上的響應(yīng)產(chǎn)生光譜值的變化,為變色松樹的遙感監(jiān)測提供了可能。

    近年來,我國高分系列衛(wèi)星的發(fā)射,為我國林業(yè)有害生物的監(jiān)測提供了很好的數(shù)據(jù)源,如可監(jiān)測大尺度的高分一號衛(wèi)星數(shù)據(jù)[4],可監(jiān)測單株樹級別的高分二號衛(wèi)星數(shù)據(jù),可對樹高和水分含量進(jìn)行監(jiān)測的高分三號雷達(dá)數(shù)據(jù),可用于森林病蟲害定量監(jiān)測的高分五號衛(wèi)星數(shù)據(jù)等。無人機遙感作為一種快速、實時、可靠、無破壞性且便于大面積覆蓋的監(jiān)測技術(shù),在第3 次全國林業(yè)有害生物普查工作中也已有應(yīng)用[5]。此外,計算機處理能力的不斷提高,尤其是機器視覺領(lǐng)域的圖像處理算法和硬件的發(fā)展,為衛(wèi)星和無人機遙感影像的應(yīng)用創(chuàng)造了很好的條件。松樹枯萎的原因較多,如干旱脅迫、山松甲蟲、松材線蟲病、衰老等。關(guān)于山松甲蟲導(dǎo)致的松樹枯萎遙感衛(wèi)星和無人機遙感監(jiān)測研究在北美地區(qū)開展了大量系統(tǒng)的研究[6-7],而關(guān)于松材線蟲病的致病機理和防控措施在日本也有大量相關(guān)研究[8],但有關(guān)松材線蟲病變色松樹的遙感監(jiān)測研究尚較分散。

    本文對有關(guān)松材線蟲病變色松樹遙感監(jiān)測研究主題中涉及地面調(diào)查、衛(wèi)星遙感、無人機遙感和新出現(xiàn)的遙感解譯方法等方面的研究現(xiàn)狀與存在的問題進(jìn)行了綜合評述,思考一些未來可能的研究方向,最后提出一套適合于我國山區(qū)復(fù)雜環(huán)境星-機-地多源數(shù)據(jù)協(xié)同監(jiān)測變色松樹的方案。

    1 松材線蟲病病感樹的生理生態(tài)特征和冠層光譜變化

    1.1 寄主松樹病感階段的生理生態(tài)特征變化

    松樹遭受松材線蟲病后,不同樹種間表現(xiàn)出一定的差異,但主要有兩種表現(xiàn)形式:一是樹木外部形態(tài)發(fā)生變化,如針葉褪色,枯萎;二是樹木內(nèi)部光合生理發(fā)生變化,如葉綠素含量和水分含量降低,光合作用和蒸騰作用下降。這些外部形態(tài)或內(nèi)部光合生理變化導(dǎo)致森林光譜反射和輻射特征的變化,在遙感影像上表現(xiàn)為光譜值的變化。這種內(nèi)外部變化和光譜值變化之間的關(guān)系構(gòu)成了遙感影像進(jìn)行變色松樹監(jiān)測的理論基礎(chǔ)。由光合生理參數(shù)的變化最終將表現(xiàn)為寄主松樹外部特征的變化[9],相關(guān)變化整理結(jié)果如表1 所示。

    病死疫木的病害發(fā)展過程分4 個階段,不同階段采用的遙感監(jiān)測方法會不同。

    第1 階段外觀正常,樹脂分泌減少,蒸騰作用下降,在嫩枝上可見天??惺硺淦さ暮圹E。這一階段松樹仍表現(xiàn)為綠色,健康松樹和變色松樹的光譜信號差異微弱。為了檢測這種變化,單個影像像素內(nèi)的地物數(shù)量必須最小化,相對差異必須最大化,這需要一個具有高空間分辨率和高光譜分辨率的傳感器。無人機機載高光譜的方式為這種早期監(jiān)測提供了可能。目前常采用紅邊數(shù)據(jù)或雷達(dá)數(shù)據(jù)的“水分含量”指示值因子實現(xiàn)疫木的早期監(jiān)測[10-11],但目前的監(jiān)測精度不高,最精確的方式為防護(hù)工人地面調(diào)查取樣然后實驗室分析。Roberts et al[12]的研究表明無人機的飛行時間對于松材線蟲病的早期監(jiān)測是非常關(guān)鍵的,越早監(jiān)測到變色松樹,其監(jiān)測的誤差就越大,減少這類誤差可采用多時相方法[14]。此外,云量、干旱脅迫、積雪等環(huán)境因素可能會進(jìn)一步削弱這一階段的遙感監(jiān)測能力。

    第2 階段針葉開始變色,樹脂分泌停止,除見天牛補充營養(yǎng)痕跡外,還可發(fā)現(xiàn)產(chǎn)卵刻槽及其它甲蟲侵害的痕跡。這一階段的松樹上會出現(xiàn)紅葉松針和綠葉松針混合的現(xiàn)象,且紅葉松針較少。若要實現(xiàn)精確的監(jiān)測,需要較高空間分辨率和較高光譜分辨率的遙感影像。

    第3 階段大部分針葉變?yōu)辄S褐色,萎蔫,可見到天牛及其它甲蟲的蛀屑。在某些條件下,利用像素的光譜混合可以檢測到這一階段的變色松樹和健康松樹之間的光譜差異[13]。如果受到感染的變色松樹冠層直徑較大,且被感染的變色松樹較為集中,則可以提高變色松樹的監(jiān)測能力。因為變色樹冠層直徑越大,傳感器所需的空間分辨率和光譜分辨率越小。

    第4 階段針葉全部變?yōu)辄S褐色至紅褐色,病樹整株干枯死亡,但針葉不脫落。對于第3、4 階段的變色松樹,現(xiàn)有的衛(wèi)星遙感監(jiān)測和無人機遙感監(jiān)測均能實現(xiàn)較高精度的識別。

    表1 松樹接種松材線蟲后的特征變化Table 1 Variation of characteristics pine trees before and after the inoculation of Bursaphelenchus Xylophilus

    1.2 松材線蟲病變色松樹冠層光譜

    不同的植物具有不同的反射光譜特征,且同種植物也會因植物種類、生長階段、葉綠素含量、葉片含水量、健康狀況而異[15]。因此,采用高光譜技術(shù)間接診斷松材線蟲病的方法可為大區(qū)域遙感監(jiān)測提供可靠依據(jù)[16-17]。目前基于高光譜的松材線蟲病監(jiān)測目的主要是為了實現(xiàn)對松材線蟲病的早期監(jiān)測,如Ju et al[18]采用松樹冠層反射光譜數(shù)據(jù)來區(qū)分染病松樹和健康松樹、確定其敏感波段或波段組合并分析葉綠素含量和光譜指數(shù)在不同病感階段的關(guān)系以推測松樹的染病階段。Kim et al[19]通過分析松材線蟲病感染的松樹光譜變化確定早期監(jiān)測的最優(yōu)光譜指數(shù)。黃明祥等[20]連續(xù)觀測了健康和發(fā)病馬尾松針葉的時序高光譜數(shù)據(jù),分析松材線蟲病害的光譜時序特征、最佳觀測時期和最佳光譜指數(shù)。徐華潮等[21]采用高光譜數(shù)據(jù)分析了不同病感階段松材線蟲病害的敏感光譜波段和葉綠素含量變化。

    這些研究都證實了紅邊波段(680~780 nm)的位置和斜率變化是描述植物健康狀況的重要指示參數(shù),可被用于松樹健康狀況的探測。從圖1a-c 可以看出[22],紅邊斜率越大,松樹(馬尾松、黑松和黃山松)的生長狀態(tài)越好,紅邊位置會出現(xiàn)“紅移”;反之,則紅邊位置會相應(yīng)“藍(lán)移”。此外,光譜的一階、二階和高階導(dǎo)數(shù)具有消除背景噪聲,分辨重疊光譜,提取光譜參數(shù)的特點,如提取吸收峰位置,紅邊位置等,可用于研究植被受病蟲害的侵害導(dǎo)致的生化組分變化。

    圖1 不同松樹不同病感階段冠層光譜反射率Fig.1 The canopy spectral reflectance in different stages of affection and different pinus L.

    2 基于衛(wèi)星遙感的變色松樹監(jiān)測

    2.1 衛(wèi)星遙感數(shù)據(jù)源

    目前衛(wèi)星遙感技術(shù)仍然是森林病蟲害監(jiān)測的重要手段之一。光學(xué)衛(wèi)星遙感由于具有多波段和多時相的特點,對于變色松樹的監(jiān)測、定位、疫情評估均有優(yōu)勢。Wulder et al[23]以一個管理者的角度,提出了區(qū)域、景觀、局部3 種尺度山松甲蟲災(zāi)害變色松樹遙感監(jiān)測的時空需求,并對遙感監(jiān)測手段的優(yōu)缺點進(jìn)行了綜合論述。歐洲合資研究中心(JRC)在2014 年和2015 年進(jìn)行了一項先導(dǎo)性研究[24],建立一套基于衛(wèi)星遙感的方法監(jiān)測受松材線蟲病感染的樹木。JCR 收集了該時間段不同傳感器來源和不同分辨率的衛(wèi)片和航片數(shù)據(jù)。目的是評價這些影像是否能夠 1)區(qū)分健康冠層和輕度枯死的冠層;2)區(qū)分不同級別的枯死。結(jié)果表明通過對這些影像進(jìn)行合適的處理,尤其是一些高空間分辨率和高光譜分辨率的航片足以識別松樹冠層的枯死。目前大多數(shù)的基于遙感進(jìn)行變色松樹監(jiān)測的對象均是山松甲蟲,表2 中總結(jié)了山松甲蟲衛(wèi)星遙感監(jiān)測的數(shù)據(jù)源、采用的方法和識別精度等,以期對我國的松材線蟲病遙感監(jiān)測能夠提供一些借鑒意義。

    從表2 中可以看出,目前,常用于森林病蟲害的衛(wèi)星影像主要有Landsat TM/ETM+、QuickBird、IKONOS、Spot、WorldView 等遙感數(shù)據(jù)。對于粗空間分辨率衛(wèi)星影像的應(yīng)用方面,大量的研究證實了Landsat TM 或ETM+影像能夠?qū)ι剿杉紫x侵害的林分尺度的松樹實現(xiàn)70%~85%的監(jiān)測精度[15,30-32],但如果對于發(fā)病的晚期受災(zāi)面積比較大時采用Landsat 可以實現(xiàn)超過85%的精度[32]。Landsat 影像的空間分辨率很低,由此產(chǎn)生的混合像元問題給變色松樹的監(jiān)測帶來挑戰(zhàn),但是多時相的Landsat影像彌補了空間分辨率上的不足。這些研究中Landsat 的多光譜數(shù)據(jù)經(jīng)過纓帽變換(TC)得到明度、綠度和濕度后尋找與變色松樹數(shù)量之間的相關(guān)關(guān)系。Sharma and Murtha[35]論證了TC 變換后的參數(shù)在變色松樹和健康松樹之間存在顯著的差異。除了單景Landsat 影像的應(yīng)用外,多時相的影像經(jīng)過TCT 變換后可以用于評價山松甲蟲的損害程度。Franklin et al[13]使用自動檢測算法對30 m Landsat TM 影像的變色松樹識別,精度達(dá)到73.3%。Skakun et al.[29]使用多時相Landsat ETM+影像,采用TC 變換和EWDI 增強濕度差異指數(shù)對紅色變色松樹的解譯精度達(dá)到67%~78%。

    表2 變色松樹衛(wèi)星遙感監(jiān)測數(shù)據(jù)源Table 2 Data sources of satellite remote sensing imagery for red attacked pine trees

    對于高空間分辨率的衛(wèi)星應(yīng)用方面,White et al.[32]使用一種無監(jiān)督聚類方法(ISODATA)在4 m多光譜IKONOS 圖像上對損害程度為中低水平的山松甲蟲監(jiān)測,并與航空攝影影像的解譯結(jié)果進(jìn)行對比。結(jié)果表明對于松樹損害程度低的區(qū)域該方法實現(xiàn)70.1%的變色松樹識別精度,對于中度損害區(qū)域可實現(xiàn)92.5%的識別精度。Hicke et al.[33]應(yīng)用單時相或多時相QuickBird 遙感影像,認(rèn)為綠色指數(shù)(GVI)、紅綠指數(shù)(RGI) 及其差值,可以探測山松甲蟲危害的林木,分析松甲蟲的種群變化情況和評估治理成效等。

    對于高光譜衛(wèi)星應(yīng)用方面White et al.[36]使用EO-1 高光譜數(shù)據(jù)計算6 個水分指數(shù),并比較其與紅葉侵襲像素之間的相關(guān)性。結(jié)果表明含有短波紅外(SWIR)和近紅外(NIR)區(qū)域的高光譜水分指數(shù)與變色松樹疫情級別顯著相關(guān)。

    2.2 基于衛(wèi)星多光譜遙感數(shù)據(jù)的變色松樹監(jiān)測方法

    就分類方法而言,隨著研究尺度的不同,基于重訪周期、空間分辨率和光譜分辨率的考慮采用的數(shù)據(jù)源會不同,因此選用的分類方法也會不同。目前變色松樹衛(wèi)星遙感監(jiān)測方法主要包括基于像素的影像分類方法、面向?qū)ο蠓椒?、深度學(xué)習(xí)方法、植被指數(shù)法等。

    (1)基于像素的影像分類方法

    早期的衛(wèi)星影像分類主要針對的是低分辨率的影像,多采用基于像素的分類方法,包括K-Means、支持向量機法(SVM)、最大似然方法(MLC)、BP 神經(jīng)網(wǎng)絡(luò)、ISODATA 聚類等。Hicke et al.[33]利用Quick Bird 高分影像的4 波段圖像和MLC 分類方法成功提取了大面積山松甲蟲侵襲致死的白皮松狀況,并將提取制圖結(jié)果分為非植被、草地、健康松樹和感染松樹4 類。Dennison et al.[26]采用MLC方法和GeoEye-l 高分影像對美國黑松森林區(qū)域的病害松樹(灰色或紅色)和健康松樹(綠色)進(jìn)行分離。Poona et al.[37]在Quick Bird 數(shù)據(jù)上提取多種植被指數(shù)和穗帽變換成分后,將神經(jīng)網(wǎng)絡(luò)方法和貝葉斯分類器結(jié)合應(yīng)用于南非地區(qū)松樹脂潰瘍病感染樹木的提取。White et al.[32]利用ISODATA 方法和IKONOS 高分圖像對山松甲蟲的紅色攻擊區(qū)域進(jìn)行了監(jiān)測。由于松材線蟲病引起的變色松樹是一個變化的過程,難以利用單景高分衛(wèi)星影像實現(xiàn)高的分類精度,Wulder et al.[14]認(rèn)為可以利用多時相的高分?jǐn)?shù)據(jù)以獲得更高的分類精度。

    (2)面向?qū)ο蠓椒?/p>

    基于像素的分類方法是以單個像元為基本分類單元,通常使用的是影像的像元光譜信息,而不考慮空間信息,并且大多已經(jīng)發(fā)展成熟并打包到很多商業(yè)軟件中。隨著高分辨率影像的可獲得性提高,傳統(tǒng)的基于像素的分類方法對高分影像的對象提取結(jié)果噪音較大,并且缺乏提取對象的語義特征,難以實現(xiàn)對象級別的提取。面向?qū)ο蟮姆椒ǖ某霈F(xiàn)很好地彌補了這類分類方法的缺陷,目前己經(jīng)在高分影像上的應(yīng)用上取得較好的效果[38]。面向?qū)ο蟮姆诸惙椒ㄌ幚淼氖且幌盗邢裨募希ɑ蚍Q對象)。該類方法在分類前的第一步是對像元級的高分影像進(jìn)行分割,分割后的單元視為一個個對象。在森林病害感染樹木提取中,一些學(xué)者嘗試采用面向?qū)ο蠓椒ㄌ岣叻诸惖木取uo et al.[39]利用面向?qū)ο蠓椒ńY(jié)合專家知識分類法,成功地探測了櫟猝死病導(dǎo)致的病死樹區(qū)域。Johnson et al.[40]在監(jiān)測日本松墨天牛和小蠹蟲引起的樹木枯萎狀況時,采用多尺度的面向?qū)ο髨D像分類方法從HIS 融合后的高分衛(wèi)星影像上區(qū)分出變色樹和健康樹。

    (3)深度學(xué)習(xí)方法

    變色松樹的識別,受到噪音、光照、季節(jié)和很多其他因素的影響,給高分辨率的遙感分類帶來技術(shù)難度。傳統(tǒng)的基于像素的遙感影像分類方法需要人工選擇一些特征圖層,基于這些特征圖層形成劃分各類別的特征空間,然后再設(shè)計一個好的分類器實現(xiàn)不同類別的分類(如SVMs,K-NN)。但這些方法只能提取到原始影像數(shù)據(jù)的淺層特征(如光譜特征和紋理空間特征等),難以滿足山區(qū)復(fù)雜地形條件下地物分類的精度要求。即使使用面向?qū)ο蟮姆椒?,在有些混合像元較為嚴(yán)重的變色松樹和稀疏植被裸土交錯帶的分類結(jié)果也不是很理想。影響面向?qū)ο蠓椒ńY(jié)果精度的主要原因是分類特征的選擇,由于人為選擇分類特征的不具代表性,導(dǎo)致面向?qū)ο蠓诸惖慕Y(jié)果不理想的研究很多。近年來發(fā)展起來的深度學(xué)習(xí)方法[41-42]為解決這類問題提供了很好的解決方案。深度學(xué)習(xí)方法具有自學(xué)習(xí)功能,能夠?qū)W習(xí)到一些高級抽象的空間特征或光譜特征,其具有更強的“泛化”能力[41],因此能夠提高復(fù)雜環(huán)境下的地物分類精度。

    目前深度學(xué)習(xí)方法應(yīng)用于遙感領(lǐng)域主要包括影像的預(yù)處理[43-44]、基于像素的分類、目標(biāo)識別和場景理解4 個方面。而對于變色松樹的識別屬于目標(biāo)識別領(lǐng)域。隨著衛(wèi)星影像的空間分辨率的提高,遙感影像的分析開始從單個像素轉(zhuǎn)向地理對象,甚至語義級別的地理對象(如房屋,松樹,交通工具等)。基于高分衛(wèi)星影像的目標(biāo)識別要比自然照片中的目標(biāo)識別更具挑戰(zhàn)。一方面是因為自然照片中的目標(biāo)對象尺寸更大,而光學(xué)遙感影像上的目標(biāo)對象通常只有幾個至幾十個像素不等的小對象特征。Ren et al.[45]嘗試對Faster R-CNN 模型進(jìn)行修正以適應(yīng)遙感圖像小對象的較高精度識別。另一方面,大尺度的光學(xué)遙感影像對象檢測受復(fù)雜背景地物的干擾較大,難以實現(xiàn)高精度的識別結(jié)果。此外,遙感領(lǐng)域很難獲取特定對象的大量樣本,使得深度學(xué)習(xí)網(wǎng)絡(luò)提取的對象特征“泛化”能力不強[46]。Yang et al.[47]借鑒人類認(rèn)知過程中的對同一對象的多角度認(rèn)識過程實現(xiàn)遙感影像上對象特征的提取。

    (4)其他輔助方法

    一些特征變化方法可將原始特征通過特定函數(shù)變換到新的特征空間,目前常用的方法有: Gra-Sdimidt 變換、主成分分析(PCA)、纓帽變換(TC)、光譜混合分析、變化向量分析[48-49]、小波變換、最小噪聲分析變化等。衛(wèi)星遙感數(shù)據(jù)還能夠很好地與其他測繪數(shù)據(jù)、地形數(shù)據(jù)、氣象數(shù)據(jù)和森林林分調(diào)查數(shù)據(jù)進(jìn)行集成應(yīng)用。為提高高分辨率衛(wèi)星影像監(jiān)測松材線蟲病變色松樹的精度,一方面可以加入一些輔助數(shù)據(jù),如地形要素、樹種分布或一些地面輔助調(diào)查數(shù)據(jù)等;另外采用衛(wèi)星遙感數(shù)據(jù)監(jiān)測松材線蟲病變色松樹的最大優(yōu)勢在于其能夠獲取疫區(qū)的多時相數(shù)據(jù)。Skakun et al.[29]通過多時相的ETM+數(shù)據(jù)和經(jīng)TC 變換導(dǎo)出的濕度,計算出濕度差指數(shù)EWDI,利用EWDI 進(jìn)行全葉分析來識別變色松樹。

    3 基于無人機的變色松樹遙感監(jiān)測

    3.1 無人機數(shù)據(jù)的獲取

    無人機(Unmanned Aerial Vehicle, UAV)是一種由動力驅(qū)動、機上無人駕駛、依靠空氣提供升力、可重復(fù)使用航空飛行器的簡稱。盡管大量研究均表明多光譜和多時相衛(wèi)星遙感影像能夠用于探測病蟲害引起的變色松樹監(jiān)測,然而受其空間分辨率、大氣影響和重訪周期的限制,對于特定時間窗口單株樹級別的變色松樹監(jiān)測難以準(zhǔn)確定位。無人機遙感是一種低空遙感監(jiān)測手段,在數(shù)據(jù)獲取過程中受大氣的影響較小、靈活性、時空分辨率均可根據(jù)需求調(diào)整[50]。隨著無人機平衡架、慣性測量單元(IMU)、全球?qū)Ш叫l(wèi)星系統(tǒng)(GNSS)和其他零部件技術(shù)的發(fā)展,無人機的定位和圖像質(zhì)量得到顯著的提升[51]。

    近年來無人機技術(shù)的發(fā)展,特別是一些關(guān)鍵技術(shù)的突破,為自然資源的遙感監(jiān)測提供了新的搭載平臺[52]。除無人機硬件設(shè)施的提高外,山區(qū)復(fù)雜環(huán)境下無人機數(shù)據(jù)的獲取在航線規(guī)劃時還需考慮相機參數(shù)、航向疊片率、旁向疊片率、空域要求等因素。起伏不大的地形區(qū)域,采取一般航線規(guī)劃技術(shù)按平行帶狀進(jìn)行航空攝影拍照即可。對地勢起伏較大的局部地區(qū),學(xué)者[53]提出從可見光相機硬件和航線規(guī)劃兩方面改進(jìn)以保證獲取的航片質(zhì)量:一個是采取通過無人機搭載多個數(shù)字相機以克服地形起伏給航片帶來的畸變;另一個是根據(jù)山地的等高線來規(guī)劃航拍的路線,按照等高線分割成小的區(qū)域依次實施航測或者采取增加重疊度和疊片率的辦法保證圖像質(zhì)量。

    3.2 基于無人機可見光遙感變色松樹監(jiān)測

    我國松材線蟲病發(fā)病或病死的疫木是空間分散的,很難采用衛(wèi)星遙感進(jìn)行監(jiān)測。而近年來發(fā)展起來的無人機遙感技術(shù)為我國變色松樹的快速定位和總量估計帶來了解決辦法。然而,目前在分散的單株變色松樹無人機影像提取中仍然停留在目視判讀的水平[53-54],完全依靠人工目視判讀變色松樹的工作方案效率低、主觀性強[55-56]?,F(xiàn)有的支持向量機[57-58]、面向?qū)ο骩59]、神經(jīng)網(wǎng)絡(luò)[60-61]等遙感數(shù)據(jù)分類提取算法雖然一定程度上提高人工目視判讀的效率,但就時間復(fù)雜度而言難以滿足GB 甚至TB 數(shù)量級的超高空間分辨率大區(qū)域無人機影像計算需求[62]。

    Cheng et al.[42]將光學(xué)遙感中地理對象的探測方法分為5 大類別:模板匹配對象檢測方法[63]、基于知識的對象檢測方法[64]、基于OBIA 對象分析的對象檢測方法[65]和基于機器學(xué)習(xí)的對象檢測方法[66-67]。卷積神經(jīng)網(wǎng)絡(luò)方法(CNNs)是機器學(xué)習(xí)方法的一種(如R-CNN 或Faster RCNN 模型的改進(jìn)),目前基于CNNs 的地理對象探測中雖然識別精度非常高,但其探測對象主要為機場[68]、汽車[47]、建筑物[69]等清晰可見的目標(biāo),針對山區(qū)復(fù)雜地形條件下變色松樹的單株木探測的精度如何還未曾有相關(guān)研究報道。

    與一些基于高分衛(wèi)星遙感的單株樹樹冠監(jiān)測的研究[70]不同,無人機高分辨率影像在分辨率提高的同時也給變色松樹提取工作帶來了一定的困難,比如裸地陰影明顯、樹冠變化過多、變色松樹和健康松樹之間邊界復(fù)雜等問題。在研究變色松樹無人機遙感監(jiān)測的精度評價時,還需要選擇由松樹與其他不同地面地物組合的不同情景下的多個樣方實驗,評價各情景下的不確定性情況以增強研究方法的魯棒性。如表3 中的常見樣方典型地物組成成份。

    3.3 基于無人機高光譜遙感的松材線蟲病疫情評價

    機載高光譜的出現(xiàn)為高分影像上“同物異譜”或“異物同譜”的情況提供了解決辦法,如區(qū)分變色松樹和其他紅色闊葉樹、稀疏植被裸土等。Kim et al.[19]基于多時相的高光譜1 m 空間分辨率航拍數(shù)據(jù),采用NDVI 和VIgreen 識別砍伐的和未砍伐的變色松樹。但由于機載高光譜的空間分辨率很高,且在數(shù)據(jù)采集時受局部視場條件和光照條件以及冠層結(jié)構(gòu)的差異,同一棵樹的冠層亮度值差異很明顯,因此變色松樹像元位置的選擇對光譜曲線的影響很大。N?si et al.[71]使用UAV 遙感制圖山松甲蟲的損害程度,發(fā)現(xiàn)使用機載高光譜影像對森林病蟲害的探測要優(yōu)于多光譜影像,主要是因為其具有詳細(xì)的連續(xù)的地物光譜信息。Leckie et al.[72]比較了整株樹樹冠,向陽部分的樹冠、陰影部分的樹冠和樹頂最亮部分4 個位置,結(jié)果發(fā)現(xiàn)整個樹冠和向陽部分的樹冠對于變色松樹的屬性反演結(jié)果最好。Coops et al.[73]在機載高光譜單株樹的選擇時采用“樹暈去中心”的方法獲得了很好的變色松樹分級結(jié)果圖。

    表3 典型地物的RGB 影像Table 3 RGB images of typical categories

    此外,高光譜遙感變色松樹監(jiān)測的另一個“核心問題”是如何選擇好的特征波段來刻畫地物光譜與空間影像特征。Fassnacht et al.[74]的變色松樹高光譜識別研究中采用遺傳算法(GA)篩選特征波段,發(fā)現(xiàn)特征波段的選擇對變色松樹分類的精度很重要。最新的研究中Zhang et al.[75]采用了3 種特征波段的選擇方法,實現(xiàn)了松樹純林環(huán)境下分段擬合變色松樹95.23%的監(jiān)測精度。然而,這種手工選擇影像特征的方法很難有效區(qū)分變色松樹和其他地物,因為人為設(shè)計的特征過多地依賴于專家知識和影像源數(shù)據(jù)本身。目前研究中所采用的特征波段選擇均是針對某種特定情景下適合的,缺乏足夠的“泛化”能力。近年來深度學(xué)習(xí)在高光譜的應(yīng)用中得到了大量地應(yīng)用[41,76],借助Hinton[77]提出的深度學(xué)習(xí)方法,遙感影像地物自動識別效率取得了較大的提高。深度學(xué)習(xí)通過多層處理機制揭示隱藏在數(shù)據(jù)中的非線性特征,能夠從大量的訓(xùn)練數(shù)據(jù)集中得到“自動學(xué)習(xí)特征”,使得遙感識別特征的選擇從“手工選擇”向“自動選擇”的轉(zhuǎn)變。

    4 松材線蟲病變色松樹遙感監(jiān)測存在的問題及展望

    4.1 監(jiān)測手段單一的問題

    松材線蟲病疫區(qū)變色松樹的監(jiān)測手段包括地面調(diào)查、無人機遙感調(diào)查、衛(wèi)星遙感調(diào)查,具體選擇哪種方法取決于調(diào)查的目的。不同的監(jiān)測手段即選擇不同的數(shù)據(jù)源,Lee et al.[78]探討了多源(IKONOS,QuickBird,航拍照片等)多時相高空間分辨率遙感數(shù)據(jù)對于松材線蟲病導(dǎo)致的變色松樹識別的潛力,發(fā)現(xiàn)不同遙感平臺的數(shù)據(jù)各有其優(yōu)勢,如衛(wèi)星遙感的多時相優(yōu)勢、無人機遙感的高空間分辨率和高靈活性優(yōu)勢、地面調(diào)查高可靠性優(yōu)勢等。未來,多種監(jiān)測手段的協(xié)同技術(shù)將成為制約松材線蟲病遙感監(jiān)測的關(guān)鍵瓶頸。

    第3 次全國林業(yè)有害生物普查工作中是以地面調(diào)查為主要手段。鑒于松材線蟲病具有“致病力強、發(fā)病時間短、傳播速度快、治理難度大”的特點,目前最有效的防治辦法是清除發(fā)病或病死疫木[79]。及時發(fā)現(xiàn)和準(zhǔn)確定位變色松樹對于病死疫木的清理及該病的防控具有重要價值。但是地面調(diào)查方式的成本高、周期長,不利于大面積變色松樹分布信息的獲取,且難以滿足大區(qū)域的應(yīng)急性需要。而且發(fā)病或病死疫木在山區(qū)是空間分散的,僅靠人工地面逐株排查不但工作量巨大,時間周期長,容易遺漏,難以全面地掌握疫情動態(tài)的問題。

    在遙感監(jiān)測方面,與加拿大采用衛(wèi)星遙感監(jiān)測大面積地勢平坦地區(qū)的山松甲蟲變色松樹不同[30],我國的松材線蟲病發(fā)病地區(qū)主要分布在中國東南部的山地區(qū)域,這些地區(qū)地勢復(fù)雜,而且光學(xué)遙感影像受云層的影響非常嚴(yán)重[80];關(guān)于無人機遙感監(jiān)測,目前大多數(shù)變色松樹無人機遙感影像的解譯工作仍然是靠人工目視判讀,這種方式對判讀人員要求很高,具有主觀性強、多人協(xié)作難以統(tǒng)一標(biāo)準(zhǔn)且難以滿足大區(qū)域快速判讀要求的特點[55]。這也一定程度上延遲了森防人員對疫木的清理,加劇了疫情擴散。找到一種適合于超高空間分辨率、大數(shù)據(jù)量的無人機影像自動處理算法對于變色松樹的快速識別、消除人工判讀誤差顯得尤為重要。

    4.2 變色松樹遙感監(jiān)測中的定位、識別和計數(shù)問題

    松材線蟲病變色松樹的遙感監(jiān)測包括識別、分級和計數(shù)3 個階段:識別階段是從遙感影像上快速定位出變色松樹的位置,并對定位的每株變色松樹是否為變色松樹進(jìn)行判定;分級階段是對變色松樹的病感階段(或稱枯死程度)進(jìn)行定量評估;計數(shù)階段是對所有定位的變色松樹按照識別的病感階段進(jìn)行分類統(tǒng)計計數(shù)。關(guān)于松材線蟲病監(jiān)測,國家層面早在2009 年已發(fā)布了技術(shù)規(guī)程《松材線蟲普查監(jiān)測技術(shù)規(guī)程》[81],該規(guī)程有利于指導(dǎo)全國松材線蟲病普查和監(jiān)測的工作。關(guān)于松材線蟲病遙感監(jiān)測,2016 年安徽省質(zhì)量技術(shù)監(jiān)督局發(fā)布了安徽省地方標(biāo)準(zhǔn)《基于無人機平臺的松材線蟲病枯死松樹監(jiān)測技術(shù)規(guī)程》[22],2018 年中國林學(xué)會預(yù)發(fā)布了行業(yè)標(biāo)準(zhǔn)《無人機遙感監(jiān)測異常變色木操作規(guī)程》,規(guī)程中規(guī)定了應(yīng)用無人機遙感監(jiān)測和識別異常變色松樹的航攝系統(tǒng)要求、航攝計劃與設(shè)計方案、攝影質(zhì)量要求、信息處理與地面驗證方法、監(jiān)測報告和檔案管理等操作規(guī)程,為松材線蟲病所致異常變色松樹的早期發(fā)現(xiàn)、疫情核查與防治效果評價提供了支持。利用光學(xué)遙感識別、分級和計數(shù)變色松樹目前亟需解決以下幾個科學(xué)問題:

    (1)定位方面包括遙感圖像拼接不佳導(dǎo)致的局部區(qū)域圖像扭曲,受地形、坐標(biāo)系轉(zhuǎn)換、風(fēng)速對無人機姿態(tài)等的影響導(dǎo)致的遙感影像幾何偏移問題;

    (2)識別方面包括變色樹漏判和誤判問題。漏判的情況如變色樹冠層被遮擋、變色樹冠層直徑太小、圖像拼接不佳導(dǎo)致的部分地區(qū)變色樹難以判別;誤判情況如把別的地物如黃色灌木、裸土、枯死草地或地面已伐未封袋死樹等判為變色樹。

    (3)計數(shù)方面包括漏判導(dǎo)致的少計數(shù)和誤判導(dǎo)致的多計數(shù)。無人機影像上統(tǒng)計的結(jié)果總是低估了變色松樹的實際總量,究其原因:一是多株變色松樹在無人機遙感影像上邊界模糊,因此容易將多株相鄰的變色松樹判為1 株,造成漏判;另一方面,處于下層的變色松樹被上層的健康松樹遮擋,而無法在遙感影像上體現(xiàn),從而導(dǎo)致下層的變色松樹漏判。

    4.3 多源遙感數(shù)據(jù)協(xié)同監(jiān)測松材線蟲病變色松樹的展望

    為結(jié)合多種調(diào)查手段的優(yōu)勢,彌補其不足,未來的研究中可采用“星-機-地”多源數(shù)據(jù)協(xié)同監(jiān)測變色松樹(如圖2)。首先采用地面調(diào)查統(tǒng)計數(shù)據(jù)和非成像高光譜數(shù)據(jù),精細(xì)DEM 數(shù)據(jù)作為衛(wèi)星遙感影像校正和疫情分區(qū)分級的基礎(chǔ),然后通過疫情分級專題圖得到的疫情嚴(yán)重區(qū)邊界作為指導(dǎo)無人機影像作業(yè)的范圍,基于獲取的無人機影像采用適合于超高分辨率、大數(shù)據(jù)量的無人機影像高效處理的深度學(xué)習(xí)算法對單株疫木進(jìn)行定位和計數(shù),最后導(dǎo)出單株疫木的經(jīng)緯度坐標(biāo)指導(dǎo)疫情嚴(yán)重區(qū)的地面疫木砍伐工作。

    (1)基于地面調(diào)查數(shù)據(jù)和衛(wèi)星影像提取疫情發(fā)生前的松樹邊界

    基于衛(wèi)星遙感手段的大區(qū)域松樹的空間分布信息是松材線蟲病防治的第一步。然而,當(dāng)前基于光學(xué)遙感提取山區(qū)松林仍受山區(qū)陰影的嚴(yán)重影響,這給利用光學(xué)遙感技術(shù)提取山區(qū)松樹分布帶來困難。地形校正方法是一種減弱山區(qū)地形陰影對遙感光譜值影響的常用方法,可提高光學(xué)遙感影像的數(shù)據(jù)質(zhì)量。可利用資源三號立體像對提取的高分辨率DEM 來計算不同地形校正模型的地形校正參數(shù),對校正后的多時相高分辨率光學(xué)遙感影像進(jìn)行分類提取,從而準(zhǔn)確地獲取疫情發(fā)生前山區(qū)松樹空間分布信息。

    (2)基于地面非成像光譜數(shù)據(jù)和衛(wèi)星影像識別疫情嚴(yán)重區(qū)邊界

    受災(zāi)范圍、面積和受災(zāi)程度的動態(tài)監(jiān)測和損失評估是松樹疫情監(jiān)測的一項重要工作,尤其是大尺度松材線蟲病的防治對于森林防護(hù)部門來說是一大挑戰(zhàn)。由于松材線蟲病的傳播速度快,疫情嚴(yán)重地區(qū)是需要最先采取措施的區(qū)域,因此區(qū)域疫情的分級工作尤為重要??刹捎玫孛婧褪覂?nèi)高光譜定量測量,測定不同程度枯死松針葉和健康松針葉區(qū)分的敏感波段,構(gòu)建松材線蟲病疫情分級的光譜指數(shù),采用多時相高分辨率衛(wèi)星遙感數(shù)據(jù)實現(xiàn)疫情的分區(qū)分級研究。對于不同疫級的區(qū)域采用不同的防護(hù)措施,如對疫情嚴(yán)重區(qū)域采取無人機進(jìn)行監(jiān)測以識別需要砍伐的單株疫木位置。

    圖2 多源遙感協(xié)同的松材線蟲病變色松樹監(jiān)測框架Fig.2 Framework of PWD monitoring by multi-source data collaboration

    (3)基于地面數(shù)據(jù)和無人機影像數(shù)據(jù)定位單株疫木并總量計數(shù)統(tǒng)計松材線蟲病的防治在目前來說最有效的辦法是清除發(fā)病和病死疫木。而病死疫木在松材線蟲病發(fā)病的早期階段是分散的,森林防護(hù)人員在山地區(qū)域?qū)ふ乙吣敬嬖诠ぷ髁烤薮?、且容易遺漏的問題。結(jié)合深度學(xué)習(xí)目標(biāo)檢測的方法和無人機影像實現(xiàn)對單株疫木的智能自動定位與計數(shù),可為林業(yè)疫木砍伐工程實施提供基礎(chǔ)數(shù)據(jù)支撐。

    猜你喜歡
    松材線蟲病變色
    肉牛消化道線蟲病和呼吸道線蟲病的流行病學(xué)、臨床癥狀和防治
    家畜類圓線蟲病的分析、診斷和治療
    《松材線蟲病生態(tài)災(zāi)害督辦追責(zé)辦法》出臺
    綠色中國(2019年13期)2019-11-26 07:11:02
    變色的真相
    變色花
    童話世界(2018年23期)2018-08-23 09:13:16
    松材線蟲病的發(fā)生防治及對策——以重慶市為例
    我國松材線蟲病的主要防治技術(shù)探析
    為什么削皮后的蘋果易變色
    獼猴桃根結(jié)線蟲病的發(fā)生與防治
    吉安市松材線蟲病發(fā)生規(guī)律初探
    成人鲁丝片一二三区免费| 免费搜索国产男女视频| 美女cb高潮喷水在线观看| 最近的中文字幕免费完整| 精品久久久久久久久久免费视频| 国产精品一区二区在线观看99 | 热99re8久久精品国产| 国产高清三级在线| 亚洲国产精品合色在线| www日本黄色视频网| 九色成人免费人妻av| 狠狠狠狠99中文字幕| 能在线免费观看的黄片| 国产日韩欧美在线精品| 成年版毛片免费区| 精品国产三级普通话版| 国产高清有码在线观看视频| 久久草成人影院| 18+在线观看网站| 99久久久亚洲精品蜜臀av| 国产精品,欧美在线| 亚洲五月天丁香| 亚洲av电影不卡..在线观看| 看十八女毛片水多多多| 国产蜜桃级精品一区二区三区| 校园人妻丝袜中文字幕| 国产91av在线免费观看| 国产白丝娇喘喷水9色精品| 亚洲精品成人久久久久久| 美女xxoo啪啪120秒动态图| 麻豆乱淫一区二区| 成人永久免费在线观看视频| av免费观看日本| 免费观看人在逋| 麻豆av噜噜一区二区三区| 国产午夜精品一二区理论片| 久久欧美精品欧美久久欧美| 一级毛片电影观看 | 国产一级毛片在线| 黄色一级大片看看| 久久精品国产亚洲av香蕉五月| 综合色av麻豆| 欧美丝袜亚洲另类| 久久鲁丝午夜福利片| 91aial.com中文字幕在线观看| 激情 狠狠 欧美| 国产人妻一区二区三区在| 国产日韩欧美在线精品| 欧美又色又爽又黄视频| 乱人视频在线观看| 男人狂女人下面高潮的视频| av天堂中文字幕网| 国产精品国产三级国产av玫瑰| 日韩欧美精品免费久久| 真实男女啪啪啪动态图| 国产一级毛片七仙女欲春2| 97超视频在线观看视频| 国产人妻一区二区三区在| 国产欧美日韩精品一区二区| av天堂在线播放| 亚洲国产精品sss在线观看| 久久精品综合一区二区三区| 久99久视频精品免费| 美女被艹到高潮喷水动态| 1000部很黄的大片| 成人三级黄色视频| 日韩三级伦理在线观看| 国产精品国产三级国产av玫瑰| 免费不卡的大黄色大毛片视频在线观看 | 亚洲熟妇中文字幕五十中出| 亚洲欧洲日产国产| 校园人妻丝袜中文字幕| 性插视频无遮挡在线免费观看| 久久午夜亚洲精品久久| 乱系列少妇在线播放| 国产麻豆成人av免费视频| 国产视频内射| 国内久久婷婷六月综合欲色啪| 欧美变态另类bdsm刘玥| 大型黄色视频在线免费观看| 国产精品日韩av在线免费观看| 日韩亚洲欧美综合| 亚洲丝袜综合中文字幕| 免费看a级黄色片| 亚洲18禁久久av| 国产一区二区亚洲精品在线观看| 久久久久久伊人网av| 国国产精品蜜臀av免费| 欧美一区二区亚洲| 欧美zozozo另类| av免费在线看不卡| 欧美一区二区亚洲| 国产黄色视频一区二区在线观看 | 午夜福利在线观看吧| 51国产日韩欧美| 国产成人精品一,二区 | 久久人人精品亚洲av| 免费av毛片视频| 午夜亚洲福利在线播放| 国产中年淑女户外野战色| 中文在线观看免费www的网站| 日韩欧美在线乱码| 亚洲人成网站在线播放欧美日韩| 边亲边吃奶的免费视频| 国产 一区精品| 亚洲一级一片aⅴ在线观看| 国产精品精品国产色婷婷| 全区人妻精品视频| 美女脱内裤让男人舔精品视频 | 国产精品乱码一区二三区的特点| 中国美白少妇内射xxxbb| 亚洲国产欧美人成| 99久久久亚洲精品蜜臀av| 亚洲中文字幕日韩| av免费观看日本| 热99在线观看视频| 日日干狠狠操夜夜爽| 一本久久中文字幕| a级毛色黄片| 观看免费一级毛片| 少妇高潮的动态图| 午夜福利在线在线| 黑人高潮一二区| 午夜亚洲福利在线播放| 午夜免费激情av| 99热全是精品| 天堂网av新在线| 免费观看精品视频网站| 精品少妇黑人巨大在线播放 | 国产综合懂色| 国产探花极品一区二区| 欧美区成人在线视频| 午夜福利高清视频| 亚洲成人中文字幕在线播放| 99热只有精品国产| 亚洲精品久久国产高清桃花| 大又大粗又爽又黄少妇毛片口| 男人舔奶头视频| 天美传媒精品一区二区| 日韩欧美国产在线观看| 午夜福利在线在线| 高清毛片免费观看视频网站| 国产探花在线观看一区二区| 久久99精品国语久久久| 超碰av人人做人人爽久久| 国产精品久久电影中文字幕| 日本免费一区二区三区高清不卡| av视频在线观看入口| 十八禁国产超污无遮挡网站| 中文欧美无线码| 好男人视频免费观看在线| 久久久久久久亚洲中文字幕| 97人妻精品一区二区三区麻豆| videossex国产| 精品国内亚洲2022精品成人| 久久亚洲精品不卡| 久久久久性生活片| 免费av不卡在线播放| 久久久久久久久大av| 成人av在线播放网站| 91精品一卡2卡3卡4卡| 日韩精品有码人妻一区| 国产中年淑女户外野战色| 色综合色国产| 亚洲国产欧美在线一区| 特大巨黑吊av在线直播| 麻豆久久精品国产亚洲av| 舔av片在线| 欧美另类亚洲清纯唯美| 亚洲精品日韩在线中文字幕 | 寂寞人妻少妇视频99o| 激情 狠狠 欧美| 日本三级黄在线观看| 亚洲中文字幕一区二区三区有码在线看| 亚洲在久久综合| 久久精品国产自在天天线| 久久午夜亚洲精品久久| 大香蕉久久网| 日韩强制内射视频| 免费看av在线观看网站| 欧美人与善性xxx| 国产精品久久久久久精品电影| 在线免费十八禁| 日本熟妇午夜| 久久久精品94久久精品| 中文字幕人妻熟人妻熟丝袜美| 日韩欧美 国产精品| 久久久久久大精品| 国产 一区 欧美 日韩| 日韩一本色道免费dvd| 欧美zozozo另类| 亚洲三级黄色毛片| 搡女人真爽免费视频火全软件| 国产精品嫩草影院av在线观看| 国产伦在线观看视频一区| 我的女老师完整版在线观看| 免费看日本二区| 国产极品精品免费视频能看的| 最近视频中文字幕2019在线8| 色尼玛亚洲综合影院| 男的添女的下面高潮视频| 观看美女的网站| 久久久久久久亚洲中文字幕| 精品国产三级普通话版| 国产麻豆成人av免费视频| 伦理电影大哥的女人| 亚洲18禁久久av| 又爽又黄无遮挡网站| 夜夜看夜夜爽夜夜摸| 成人漫画全彩无遮挡| 午夜福利成人在线免费观看| 国产精品综合久久久久久久免费| 一本久久精品| 尾随美女入室| 国产男人的电影天堂91| 亚洲av电影不卡..在线观看| 国产乱人视频| 国产亚洲精品久久久com| 尾随美女入室| 亚洲自偷自拍三级| 午夜福利成人在线免费观看| 日本与韩国留学比较| 国产成人影院久久av| 欧美日本视频| 色哟哟·www| 国产精品久久电影中文字幕| 精品人妻视频免费看| 日本欧美国产在线视频| av福利片在线观看| 亚洲av一区综合| 国产亚洲av片在线观看秒播厂 | 日韩成人伦理影院| 99在线人妻在线中文字幕| 极品教师在线视频| 久久久午夜欧美精品| 一本久久精品| 亚洲av二区三区四区| 丰满乱子伦码专区| 日韩视频在线欧美| 亚洲综合色惰| 亚洲人成网站在线播| 亚洲精品自拍成人| 国产片特级美女逼逼视频| 久久久久久久久久久免费av| 中文字幕熟女人妻在线| 大又大粗又爽又黄少妇毛片口| 午夜a级毛片| 美女 人体艺术 gogo| 成人亚洲欧美一区二区av| 国产精品精品国产色婷婷| 国产成人a区在线观看| 精华霜和精华液先用哪个| 可以在线观看的亚洲视频| 欧美最新免费一区二区三区| а√天堂www在线а√下载| av专区在线播放| 久久久久久伊人网av| 国产av不卡久久| 亚洲精品456在线播放app| 欧美精品一区二区大全| 亚洲精品国产av成人精品| 人妻制服诱惑在线中文字幕| 久久综合国产亚洲精品| 少妇的逼水好多| 91麻豆精品激情在线观看国产| 性欧美人与动物交配| 日韩欧美 国产精品| 好男人视频免费观看在线| 伦精品一区二区三区| 国产av不卡久久| 国产乱人偷精品视频| 黄片无遮挡物在线观看| 国产麻豆成人av免费视频| 熟女电影av网| 97热精品久久久久久| 亚洲在久久综合| 国产视频首页在线观看| 亚洲丝袜综合中文字幕| 国产av不卡久久| 级片在线观看| 亚洲综合色惰| 尤物成人国产欧美一区二区三区| 久久这里只有精品中国| 精品久久久久久成人av| 看非洲黑人一级黄片| 亚洲成a人片在线一区二区| 国产精品一区二区性色av| 欧美xxxx黑人xx丫x性爽| 少妇熟女aⅴ在线视频| 狂野欧美激情性xxxx在线观看| 国产色爽女视频免费观看| 看非洲黑人一级黄片| 桃色一区二区三区在线观看| 神马国产精品三级电影在线观看| 日韩 亚洲 欧美在线| 欧美三级亚洲精品| 亚洲精品国产av成人精品| 免费无遮挡裸体视频| 观看美女的网站| 日韩欧美在线乱码| 狂野欧美激情性xxxx在线观看| 日韩大尺度精品在线看网址| 亚洲av免费在线观看| 日本三级黄在线观看| 免费观看a级毛片全部| 午夜激情欧美在线| 亚洲精品日韩在线中文字幕 | 成人三级黄色视频| 国产精品伦人一区二区| 亚洲国产精品sss在线观看| 中文字幕精品亚洲无线码一区| 欧美三级亚洲精品| 国产亚洲精品久久久com| 变态另类成人亚洲欧美熟女| 日日摸夜夜添夜夜爱| 日韩成人伦理影院| 美女xxoo啪啪120秒动态图| 小蜜桃在线观看免费完整版高清| 国产一区二区在线av高清观看| 丝袜喷水一区| 亚洲18禁久久av| 国产综合懂色| 日韩三级伦理在线观看| 精品久久久噜噜| 国产精品,欧美在线| 国内久久婷婷六月综合欲色啪| 国产不卡一卡二| 国产一级毛片在线| 日韩中字成人| 人妻久久中文字幕网| 精品人妻视频免费看| 桃色一区二区三区在线观看| 日韩欧美精品v在线| 国产中年淑女户外野战色| 久久久久久九九精品二区国产| 久久午夜亚洲精品久久| 边亲边吃奶的免费视频| 日韩 亚洲 欧美在线| 亚洲成人久久性| 老司机影院成人| 不卡一级毛片| 简卡轻食公司| 日本免费一区二区三区高清不卡| 免费观看精品视频网站| 麻豆国产97在线/欧美| 午夜精品一区二区三区免费看| 国产精品一及| 男女那种视频在线观看| 97超视频在线观看视频| 天堂网av新在线| 精品午夜福利在线看| 老女人水多毛片| 91麻豆精品激情在线观看国产| 高清在线视频一区二区三区 | 国产成人福利小说| 啦啦啦韩国在线观看视频| 精品午夜福利在线看| 精品久久久久久久久久免费视频| 久久久久久久久大av| 久久精品国产鲁丝片午夜精品| 热99在线观看视频| 日韩欧美国产在线观看| 亚洲自偷自拍三级| 美女被艹到高潮喷水动态| 精品久久久久久久久av| 国产午夜精品一二区理论片| 美女xxoo啪啪120秒动态图| 给我免费播放毛片高清在线观看| 两个人的视频大全免费| 三级毛片av免费| 国产精品久久久久久亚洲av鲁大| 夜夜夜夜夜久久久久| 久久99精品国语久久久| 国产精品久久久久久久久免| 亚洲美女视频黄频| 国产又黄又爽又无遮挡在线| 国产亚洲91精品色在线| 国产午夜精品一二区理论片| 一级毛片久久久久久久久女| 黄色配什么色好看| 久久午夜福利片| 在线a可以看的网站| 99riav亚洲国产免费| 在线观看免费视频日本深夜| 国产精华一区二区三区| ponron亚洲| 精品人妻熟女av久视频| 2021天堂中文幕一二区在线观| 日本撒尿小便嘘嘘汇集6| 亚洲真实伦在线观看| 麻豆乱淫一区二区| 亚洲精品成人久久久久久| 国产91av在线免费观看| 色吧在线观看| 国产探花在线观看一区二区| or卡值多少钱| 天天一区二区日本电影三级| 一个人观看的视频www高清免费观看| 岛国毛片在线播放| 18禁裸乳无遮挡免费网站照片| 欧美日韩一区二区视频在线观看视频在线 | 如何舔出高潮| 日韩欧美一区二区三区在线观看| 久久久a久久爽久久v久久| 五月玫瑰六月丁香| 看黄色毛片网站| 国产精品爽爽va在线观看网站| 久久久久免费精品人妻一区二区| 日日干狠狠操夜夜爽| 成年版毛片免费区| 禁无遮挡网站| 男的添女的下面高潮视频| 午夜激情福利司机影院| 哪个播放器可以免费观看大片| 午夜福利成人在线免费观看| 我的女老师完整版在线观看| av专区在线播放| 中国国产av一级| 不卡视频在线观看欧美| 国产成人a∨麻豆精品| 美女脱内裤让男人舔精品视频 | 国产精品1区2区在线观看.| 99riav亚洲国产免费| 精品一区二区三区人妻视频| 免费在线观看成人毛片| 亚洲成人中文字幕在线播放| 国产在视频线在精品| 国产亚洲91精品色在线| 欧美日本视频| 国产成人a∨麻豆精品| 高清午夜精品一区二区三区 | 国产一区二区三区av在线 | 国产黄a三级三级三级人| 亚洲人成网站在线观看播放| 欧美日韩乱码在线| 免费av不卡在线播放| 亚洲最大成人手机在线| 在线观看av片永久免费下载| 亚洲精品456在线播放app| 亚洲欧美成人综合另类久久久 | 亚洲电影在线观看av| 麻豆成人午夜福利视频| 亚洲欧美精品专区久久| 99久久成人亚洲精品观看| 国产一区亚洲一区在线观看| 赤兔流量卡办理| 日本av手机在线免费观看| 国产精品1区2区在线观看.| 毛片女人毛片| 又粗又爽又猛毛片免费看| 亚洲一区高清亚洲精品| 搡女人真爽免费视频火全软件| 又黄又爽又刺激的免费视频.| 成年女人永久免费观看视频| 亚洲在久久综合| 好男人视频免费观看在线| 亚洲一区二区三区色噜噜| 秋霞在线观看毛片| 欧美日韩在线观看h| 久久精品国产亚洲av天美| 免费看a级黄色片| 三级毛片av免费| 国产真实伦视频高清在线观看| 午夜福利成人在线免费观看| 亚洲自偷自拍三级| 中国美白少妇内射xxxbb| 欧美在线一区亚洲| 精品熟女少妇av免费看| 能在线免费观看的黄片| 精华霜和精华液先用哪个| 亚洲av电影不卡..在线观看| 精品午夜福利在线看| 亚洲成av人片在线播放无| 久久6这里有精品| 在线观看免费视频日本深夜| 久久精品夜色国产| 成人国产麻豆网| 黄色欧美视频在线观看| 亚洲av一区综合| 一级黄色大片毛片| 亚洲欧美中文字幕日韩二区| 五月伊人婷婷丁香| 精品久久久噜噜| 一级二级三级毛片免费看| 国产精品一区二区三区四区免费观看| 一本一本综合久久| 免费看a级黄色片| 日本在线视频免费播放| 国产一区二区在线av高清观看| 一个人看视频在线观看www免费| 日韩欧美精品免费久久| 久久久午夜欧美精品| 精品99又大又爽又粗少妇毛片| 黑人高潮一二区| 尾随美女入室| 插阴视频在线观看视频| 国产精品久久久久久亚洲av鲁大| 亚洲国产精品久久男人天堂| 久久久久久国产a免费观看| 人妻制服诱惑在线中文字幕| 在现免费观看毛片| 黄片无遮挡物在线观看| 观看免费一级毛片| av在线天堂中文字幕| 三级男女做爰猛烈吃奶摸视频| 久久久国产成人免费| 日韩欧美精品v在线| 老女人水多毛片| 嫩草影院精品99| 久久久久久久久久久免费av| 国产精品美女特级片免费视频播放器| eeuss影院久久| 日韩大尺度精品在线看网址| 国产黄片视频在线免费观看| 国内精品美女久久久久久| 国产69精品久久久久777片| 久久婷婷人人爽人人干人人爱| 精品国内亚洲2022精品成人| 亚洲欧美日韩卡通动漫| 两个人的视频大全免费| 中文字幕制服av| 精品午夜福利在线看| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲精品国产av成人精品| 人人妻人人澡欧美一区二区| 精品久久久久久久久久久久久| 一级毛片电影观看 | 精品一区二区三区人妻视频| 蜜桃亚洲精品一区二区三区| 国产成人精品婷婷| 亚洲在久久综合| 在线观看午夜福利视频| 一级毛片电影观看 | 欧美日韩国产亚洲二区| 亚洲va在线va天堂va国产| 蜜臀久久99精品久久宅男| av黄色大香蕉| 深爱激情五月婷婷| 日本-黄色视频高清免费观看| 干丝袜人妻中文字幕| 99久久中文字幕三级久久日本| 日韩成人伦理影院| 免费观看的影片在线观看| 亚洲欧美日韩无卡精品| 国产色婷婷99| 女人十人毛片免费观看3o分钟| 91久久精品电影网| 欧美成人精品欧美一级黄| 国产高清视频在线观看网站| 看免费成人av毛片| 熟女人妻精品中文字幕| 亚洲欧美成人综合另类久久久 | 亚洲真实伦在线观看| 久久精品国产鲁丝片午夜精品| 97超视频在线观看视频| 国产伦精品一区二区三区四那| 嫩草影院新地址| 亚洲成人久久性| 亚洲一区高清亚洲精品| 精品国产三级普通话版| 99热这里只有是精品50| 色综合站精品国产| 婷婷精品国产亚洲av| 久久人人爽人人片av| 国产私拍福利视频在线观看| 搞女人的毛片| 亚洲国产精品成人久久小说 | 美女被艹到高潮喷水动态| 男女视频在线观看网站免费| 亚洲精品久久久久久婷婷小说 | 尾随美女入室| 男人和女人高潮做爰伦理| 免费观看a级毛片全部| 成人欧美大片| 久久午夜福利片| 国产三级在线视频| 精品久久久久久久久久免费视频| 成人午夜高清在线视频| 亚洲精品乱码久久久久久按摩| 高清在线视频一区二区三区 | 亚洲欧美清纯卡通| 欧美日韩在线观看h| 国产精品电影一区二区三区| 六月丁香七月| 人人妻人人澡欧美一区二区| 午夜精品一区二区三区免费看| 欧美在线一区亚洲| 久久久午夜欧美精品| 精品久久久久久久久久免费视频| 亚洲aⅴ乱码一区二区在线播放| 精品人妻偷拍中文字幕| 亚洲第一电影网av| 青青草视频在线视频观看| 国产女主播在线喷水免费视频网站 | 日本免费a在线| 91久久精品电影网| 又爽又黄a免费视频| 日本免费a在线| 精品人妻偷拍中文字幕| 国产亚洲av片在线观看秒播厂 | 草草在线视频免费看| 在线国产一区二区在线| 99久国产av精品国产电影| 3wmmmm亚洲av在线观看| 免费av观看视频| 99久国产av精品国产电影| avwww免费| av在线亚洲专区| 成人av在线播放网站| 日日撸夜夜添| 日韩,欧美,国产一区二区三区 | 91狼人影院| 亚洲国产色片| 久久久久久久久大av| 久久久久国产网址|