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

    基于蒸騰速率與光譜特征的松材線蟲病害預(yù)測(cè)

    2015-02-21 17:40:53張曉麗王書涵黃雅婷
    關(guān)鍵詞:針葉蒸騰速率松材

    劉 寧,張曉麗,王書涵,黃雅婷

    (北京林業(yè)大學(xué) 林學(xué)院/省部共建森林培育與保護(hù)教育部重點(diǎn)實(shí)驗(yàn)室,北京 100083)

    基于蒸騰速率與光譜特征的松材線蟲病害預(yù)測(cè)

    劉 寧,張曉麗,王書涵,黃雅婷

    (北京林業(yè)大學(xué) 林學(xué)院/省部共建森林培育與保護(hù)教育部重點(diǎn)實(shí)驗(yàn)室,北京 100083)

    【目的】 明確馬尾松感染松材線蟲病初期蒸騰速率與光譜特征的變化規(guī)律及相互關(guān)系,以解決病害潛伏期表征不易鑒別的問題,為在遙感監(jiān)測(cè)中特征參數(shù)的選取提供理論參考。【方法】 以福建省泰寧縣馬尾松林為例,分別對(duì)感染松材線蟲病后4個(gè)病害階段(健康期及病害初期、中期、末期)馬尾松的蒸騰速率進(jìn)行測(cè)定,分析其在不同病害階段的變化規(guī)律;再對(duì)健康期及病害初期馬尾松的光譜數(shù)據(jù)進(jìn)行采集,以這2個(gè)階段的蒸騰速率為基準(zhǔn),分析其蒸騰速率與原始光譜及光譜特征參數(shù)(紅邊位置、紅邊斜率、紅谷位置、紅谷反射率、紅邊比值植被指數(shù)和紅邊歸一化植被指數(shù))的關(guān)系?!窘Y(jié)果】 馬尾松蒸騰速率隨著病害的加重而逐漸減弱;健康期及病害初期馬尾松蒸騰速率與原始光譜有一定的相關(guān)性,相關(guān)關(guān)系因波長(zhǎng)的不同而不同,最大值在753 nm 處,相關(guān)系數(shù)達(dá)到0.86;蒸騰速率與紅邊參數(shù)有較高相關(guān)性,擬合結(jié)果R2均大于0.7,其中與紅邊比值植被指數(shù)相關(guān)性最高,R2達(dá)到0.851 5?!窘Y(jié)論】 在外界因素一致的情況下,馬尾松蒸騰速率對(duì)預(yù)測(cè)病害的發(fā)生有一定的指示作用,近紅外波段及紅邊參數(shù)與蒸騰速率有較高相關(guān)性,可以利用特征光譜對(duì)松材線蟲病進(jìn)行早期預(yù)測(cè)。

    馬尾松;光譜特征;蒸騰速率;松材線蟲?。患t邊參數(shù);病害預(yù)測(cè)

    松材線蟲病是松類樹種的一種毀滅性病害,被稱為松樹的“癌癥”和“艾滋病”,已成為我國森林病蟲害之首,是我國目前發(fā)生最嚴(yán)重、危害最大的外來入侵森林病害,其傳播媒介主要為松墨天牛。感染了松材線蟲病的松樹,從發(fā)病到死亡只要2~3個(gè)月的時(shí)間,從出現(xiàn)癥狀到死亡只要40 d左右,從單株發(fā)病到整片松樹毀滅只需3~5年。在我國,松材線蟲病害的發(fā)生形勢(shì)十分嚴(yán)峻,其發(fā)病周期短、危害性極強(qiáng),給防治工作帶來了極大困難[1]。在松材線蟲病害監(jiān)測(cè)方面,雖然有部分學(xué)者已運(yùn)用遙感技術(shù)對(duì)病害區(qū)進(jìn)行了大尺度的監(jiān)測(cè)[2-11],但大多數(shù)研究局限在病害已出現(xiàn)外部癥狀的前提下,由于松材線蟲病害周期極短,因此預(yù)測(cè)預(yù)報(bào)效果難以得到保證。松材線蟲病害的發(fā)生會(huì)引起植物生理指標(biāo)的變化,如樹脂分泌減少、蒸騰速率降低、乙烯含量增加、葉綠素及水分含量變化等[12]。部分學(xué)者已對(duì)葉綠素及水分含量變化的特點(diǎn)做了相應(yīng)研究[13-19],但葉綠素及水分變化與諸多因素有關(guān),且通常與外部特征變化同時(shí)發(fā)生,因此能否用其對(duì)病蟲害進(jìn)行早期預(yù)測(cè)仍有待進(jìn)一步研究。1977年,日本學(xué)者真宮靖治與田村弘忠對(duì)感染了松材線蟲病的赤松蒸騰速率的日變化進(jìn)行了研究,結(jié)果表明,蒸騰速率起初每日下降,然后有一個(gè)劇烈下降過程,最后趨緩以至可以忽略不計(jì);而出現(xiàn)明顯枯萎和針葉變黃的時(shí)間,比蒸騰速率開始下降要延遲15~30 d[12]。說明蒸騰速率的變化先于外部特征,這一結(jié)論有助于在潛伏期內(nèi)診斷病害。馬菁等[20]利用Yaxin-1102便攜式光合蒸騰儀對(duì)不同病害時(shí)期馬尾松的蒸騰速率進(jìn)行了測(cè)定,結(jié)果顯示,隨著病害的加重,蒸騰速率在減弱。還有部分學(xué)者研究表明,病害的發(fā)生改變了植被的光譜特征,如王震等[21]運(yùn)用ASD FieldSpec HH采集了不同病害時(shí)期馬尾松的光譜特征,發(fā)現(xiàn)馬尾松不同受害類型的光譜特征變化明顯,而且呈現(xiàn)出較好的變化規(guī)律。由于植物感染松材線蟲病后其蒸騰速率變化先于外部特征,如果能運(yùn)用這一特點(diǎn)進(jìn)行病害預(yù)測(cè),對(duì)松材線蟲病害的防治將有極大幫助。為此,本研究以馬尾松蒸騰速率作為判斷松材線蟲病害嚴(yán)重程度的標(biāo)準(zhǔn),分析其與病害馬尾松光譜數(shù)據(jù)之間的關(guān)系,尋找出對(duì)病害變化敏感的波段范圍,建立生理指標(biāo)變化與光譜變化的聯(lián)系,以期為利用光譜特征對(duì)病害進(jìn)行早期預(yù)測(cè)提供理論參考。

    1 研究區(qū)概況

    泰寧縣隸屬福建省三明市,位于福建省西北部,武夷山脈中段的杉嶺支脈東南側(cè),居兩省三地市交界處。地理坐標(biāo)為北緯26°33′~27°07′,東經(jīng)116°53′~117°24′,其北靠邵武,東連將樂,南鄰明溪,西接建寧。年降雨量1 725 mm,年均氣溫17 ℃,總面積1 533.8 km2(圖1)。泰寧縣是閩西北重要林業(yè)縣之一,全縣林地面積11.47萬hm2,其中馬尾松(含混交林)面積7.4萬hm2,約占全縣林地面積的65.9%。自2007年發(fā)生松材線蟲病以來,全縣累計(jì)發(fā)生面積1 440.53 hm2,枯死木17.6萬株,主要分布在杉城(含國有林場(chǎng))、朱口和下渠3個(gè)鄉(xiāng)鎮(zhèn)[22]。

    2 研究方法

    2.1 松材線蟲病病害階段的劃分

    在多數(shù)情況下,植株染病后,于當(dāng)年秋季即表現(xiàn)為全株枯死。根據(jù)馬尾松的生長(zhǎng)狀況、外部癥狀、內(nèi)部生理指標(biāo)的變化情況及研究區(qū)松材線蟲病害現(xiàn)狀,將染病馬尾松病害程度大致分為4個(gè)階段:(1)健康期(Health period,HP)。外觀健康,樹脂分泌正常。(2)病害初期(Early period of disease,EPD)。植株少量針葉失去光澤、略變黃,樹脂分泌減少。(3)病害中期(Middle period of disease,MPD)。植株多數(shù)針葉變黃,植株開始萎蔫。(4)病害末期(Terminal period of disease,TPD)。整個(gè)樹冠部針葉呈紅褐色,植株表面多處出現(xiàn)天牛取食坑道,全株枯死[23]。

    2.2 不同病害階段馬尾松蒸騰速率的測(cè)定

    蒸騰作用是指植物體內(nèi)的水分通過體表向外散失的過程,蒸騰作用的強(qiáng)弱通常用蒸騰速率(植物在一定時(shí)間內(nèi)單位葉面積蒸騰的水量)來衡量。影響植被蒸騰速率的因素有很多,包括光照、空氣相對(duì)濕度、溫度、病蟲害及內(nèi)部的氣孔頻度、氣孔大小等[24]。為了有效地反映病蟲害對(duì)蒸騰速率的影響,避免其他因素的干擾,筆者運(yùn)用美國LI-COR公司生產(chǎn)的第3代氣體交換測(cè)量系統(tǒng)LI-6400XT進(jìn)行測(cè)定。此儀器可以測(cè)定葉片的光合速率、胞間CO2濃度、氣孔導(dǎo)度和蒸騰速率,同時(shí)具有較好地控制其他環(huán)境因子干擾的能力[25-26]。選擇LI-6400XT選配的LED紅藍(lán)光源,將光照強(qiáng)度設(shè)定為1 000 μmol/(m2·s),對(duì)其進(jìn)行恒定控制;選擇遠(yuǎn)離道路、立地條件相似、人為干擾較少的成片馬尾松林區(qū),加裝2.5 L自制緩沖瓶,減少人員走動(dòng),保證測(cè)試過程中氣流的穩(wěn)定性;選擇晴朗、天氣狀況相似的午間10:00-14:00進(jìn)行測(cè)量,減小溫度及濕度變化的影響。本試驗(yàn)選取了60株馬尾松進(jìn)行測(cè)量,處于健康及病害初期、中期、末期的馬尾松各15株,每個(gè)階段的樣本依次編號(hào)。在當(dāng)?shù)毓と说膮f(xié)助下,按照東南西北4個(gè)方向逐一剪取植株外側(cè)短枝,在每個(gè)方位的短枝中選取3束針葉,平鋪在2 cm×3 cm的標(biāo)準(zhǔn)葉室上,分別進(jìn)行離體測(cè)定(4個(gè)方位共計(jì)4次),將其平均值作為該病害階段所測(cè)樣本的蒸騰速率值,如健康時(shí)期1號(hào)樣本的蒸騰速率值,標(biāo)記為HP_1。針葉的選取原則為:在健康時(shí)期保證針葉全部健康;病害初期選取有少量針葉變黃的短枝中表征為綠色的針葉,便于后期的預(yù)測(cè)分析;病害中期選取表征變黃的針葉;病害末期選取表征變紅褐色的針葉。同時(shí),為減小時(shí)間對(duì)測(cè)定結(jié)果的影響,首先完成各病害時(shí)期馬尾松針葉1號(hào)樣本的測(cè)定,再按照相同的病害選取順序?qū)Ω鲿r(shí)期2號(hào)樣本進(jìn)行測(cè)定,依次類推,重復(fù)15次,完成對(duì)所有病害時(shí)期樣本的測(cè)定。在測(cè)定時(shí),每更換一個(gè)樣本,儀器需要重新匹配一次,以保證測(cè)定的可靠性。

    2.3 光譜數(shù)據(jù)的采集

    采用美國ASD公司生產(chǎn)的手持野外光譜輻射儀,該設(shè)備可在325~1 075 nm波長(zhǎng)連續(xù)測(cè)定,采樣間隔為1.6 nm,光譜分辨率為3.5 nm,視場(chǎng)角25°,共512個(gè)波段。測(cè)定光譜時(shí)按照上述蒸騰速率的測(cè)定順序接替進(jìn)行,將采集的東南西北4個(gè)方位共12束針葉放置在預(yù)先備好的黑布上,防止地表反射的干擾,每次測(cè)定新樣本時(shí),進(jìn)行10次光譜采集,取平均值作為該樣本的反射光譜值。以白板為標(biāo)準(zhǔn)進(jìn)行校正。由于儀器穩(wěn)定性的限制,僅在400~900 nm波長(zhǎng)進(jìn)行光譜分析。

    2.4 馬尾松蒸騰速率與光譜數(shù)據(jù)的相關(guān)性分析

    運(yùn)用的軟件包括ViewSpec Proversion 4.02、SPSS 18、Excel 2010、ENVI 4.8、MATLAB R2012b。分析方法如下:(1)選取不同病害階段馬尾松,分析其受害情況與蒸騰速率的關(guān)系;(2)運(yùn)用MATLAB R2012b軟件,分析健康和病害初期馬尾松蒸騰速率與相應(yīng)階段原始光譜的相關(guān)性,尋找預(yù)測(cè)病蟲害發(fā)生的最佳波長(zhǎng)范圍;(3)運(yùn)用ENVI 4.8+IDL對(duì)健康和病害初期馬尾松光譜進(jìn)行一階導(dǎo)數(shù)變換,分析其蒸騰速率與相應(yīng)階段光譜特征參數(shù)的相關(guān)性,尋找出預(yù)測(cè)病蟲害發(fā)生的最佳光譜特征參數(shù)。

    3 結(jié)果與分析

    3.1 不同病害階段馬尾松的蒸騰速率

    病害馬尾松蒸騰速率的變化規(guī)律見圖2。

    從圖2可以看出,隨著病害的加重,馬尾松蒸騰速率整體呈減弱的趨勢(shì)。從表1可知,蒸騰速率的標(biāo)準(zhǔn)差和標(biāo)準(zhǔn)誤均在合理范圍,說明實(shí)測(cè)的樣本數(shù)據(jù)具有一定的代表性,可以評(píng)估整體的變化規(guī)律。因此,在生長(zhǎng)條件相似、環(huán)境因素大致相同的情況下,馬尾松蒸騰速率可以作為預(yù)測(cè)松材線蟲病害發(fā)生的有效因子,其蒸騰速率越快,馬尾松越健康,反之病害越嚴(yán)重。

    3.2 馬尾松蒸騰速率與原始光譜的相關(guān)性

    圖3為健康期和病害初期馬尾松的光譜特征曲線。由圖3可知,健康期和病害初期的馬尾松在550 nm波長(zhǎng)(綠峰)處光譜差別不大,而在680 nm處及大于740 nm的范圍光譜差別明顯,進(jìn)一步說明試驗(yàn)選取的健康期及病害初期樣本在表征顏色方面差別不大,從而保證了運(yùn)用蒸騰速率預(yù)測(cè)的優(yōu)勢(shì)。

    選取試驗(yàn)中健康期及病害初期馬尾松樣本數(shù)據(jù)(共30株),將其蒸騰速率與其在400~900 nm波長(zhǎng)的光譜值做相關(guān)性分析,得到相關(guān)系數(shù)隨波長(zhǎng)變化的曲線,結(jié)果如圖4所示。

    圖3 健康期和病害初期馬尾松的原始光譜曲線
    Fig.3 Raw spectrum curve ofPinusmassonianaat HP and EPD

    由圖4可知,蒸騰速率與原始光譜反射率的最大相關(guān)系數(shù)為0.86,位于753 nm處,呈正相關(guān)。從總體上看,光譜反射值在400~521 nm及565~711 nm與蒸騰速率呈負(fù)相關(guān),且在406~507 nm、600~704 nm處達(dá)到了極顯著相關(guān)水平;在522~564 nm及712~900 nm呈正相關(guān),且在721~896 nm達(dá)到極顯著相關(guān)水平;在522~564 nm相關(guān)系數(shù)較低,沒有達(dá)到顯著水平。綜上所述,在馬尾松外部特征沒有明顯變化的前提下,隨著蒸騰速率的改變,首先變化的是藍(lán)光(425~480 nm)、紅光(600~700 nm)及近紅外(700~900 nm)這3個(gè)光譜范圍,在藍(lán)光與紅光范圍內(nèi),蒸騰速率越快,植物對(duì)此處光譜的吸收率越高,相反,蒸騰速率越慢,吸收率越低,反射率就會(huì)越高;而在近紅外范圍內(nèi),隨著蒸騰速率的加快,光譜反射率升高,相反,蒸騰速率越慢,光譜反射率越低。馬尾松原始光譜反射率與蒸騰速率的相關(guān)系數(shù)較大值主要集中在近紅外范圍內(nèi),說明近紅外波段是預(yù)測(cè)松材線蟲病害發(fā)生的最佳波長(zhǎng)范圍。

    3.3 馬尾松蒸騰速率與光譜特征參數(shù)的相關(guān)性分析

    鑒于馬尾松蒸騰速率與近紅外波段有較高的相關(guān)性,進(jìn)一步分析其與此范圍光譜特征參數(shù)之間的關(guān)系,這對(duì)松材線蟲病害的預(yù)測(cè)有著積極的作用,這些參數(shù)包括:

    (1)紅邊位置(Red edge position, REP):指綠色植物在680~740 nm處反射率增高最快點(diǎn)的位置,也是一階導(dǎo)數(shù)光譜在該區(qū)間內(nèi)拐點(diǎn)的位置。運(yùn)用ENVI 4.8+IDL調(diào)用計(jì)算一階導(dǎo)數(shù)函數(shù),獲取紅邊信息,其名稱及代碼如式(1)。(2)紅邊斜率(Red edge slope,RES):指在紅邊范圍內(nèi)一階導(dǎo)數(shù)最大峰值[14]。(3)紅谷位置(Red band position,RBP):指600~700 nm內(nèi)反射率最低值的波長(zhǎng)。(4)紅谷反射率(Reflectance of red band, RRB):指600~700 nm內(nèi)反射率最低值。(5)紅邊歸一化植被指數(shù)(Red edge normalized difference vegetation index,RENDVI):借鑒NDVI在紅邊范圍內(nèi)構(gòu)建一個(gè)指數(shù)(如式(2),式中Rx代表在x波長(zhǎng)處的反射率,下同),可以指示植被的生長(zhǎng)狀況;(6)紅邊比值植被指數(shù)(Red edge ratio vegetation index,RERVI):借鑒植被指數(shù)RVI構(gòu)建紅邊比值植被指數(shù)(如式(3))。

    pf_1st_derivative.pro:

    function pf_1st_derivative,x,y,bbl,bbl_list,

    _extra=_extra

    ptr=where(bbl_list eq 1,count)

    result=fltarr(n_elements(y))

    if(count ge 3) then $

    result(ptr)=deriv(x[ptr],y[ptr])

    return,result

    end

    (1)

    RENDVI=(R740-R680)/(R740+R680)。

    (2)

    RERVI=R740/R680。

    (3)

    圖5為馬尾松健康期及病害初期蒸騰速率與6個(gè)光譜特征參數(shù)之間的相關(guān)曲線。

    以線性函數(shù)關(guān)系模擬蒸騰速率與光譜特征參數(shù)之間的關(guān)系,結(jié)果(圖5)顯示,除紅谷反射率與蒸騰速率呈負(fù)相關(guān)外,其余光譜特征參數(shù)均與蒸騰速率呈正相關(guān),R2介于0.411 4~0.851 5。從圖5可以看出,隨著蒸騰速率的降低,紅邊位置及紅谷位置向短波方向移動(dòng),紅邊斜率降低,紅谷反射率增高,RERVI及RENDVI的值有所降低。由圖5還可知,蒸騰速率與紅邊特征參數(shù)有較高的相關(guān)性,其中相關(guān)性最高的是RERVI,R2達(dá)到了0.851 5,說明隨著蒸騰速率的變化,紅邊范圍內(nèi)光譜特征變化明顯,從而得出紅邊參數(shù),尤其是紅邊范圍內(nèi)2個(gè)植被指數(shù)是預(yù)測(cè)病蟲害發(fā)生的有效參數(shù)。

    4 結(jié)論與討論

    本研究分析了不同病害階段馬尾松蒸騰速率的變化情況,再以健康期及病害初期馬尾松蒸騰速率為基準(zhǔn),分析其與相應(yīng)時(shí)期光譜特征的相關(guān)性,尋找相關(guān)性較高的光譜波段和光譜特征參數(shù)作為預(yù)測(cè)病害的因子。其研究結(jié)果表明:

    1)松材線蟲病對(duì)馬尾松蒸騰速率有一定影響,隨著病害的加重,蒸騰速率逐漸減慢;蒸騰速率可以作為馬尾松染病前期病害診斷的依據(jù),有效解決了病害潛伏期內(nèi)從外觀不易鑒別的弊端。

    2)染病馬尾松蒸騰速率與其原始光譜有一定的相關(guān)性,在406~507 nm、600~704 nm、721~896 nm相關(guān)性較高,達(dá)到了極顯著水平,其中在波長(zhǎng)為753 nm時(shí)相關(guān)性最高,而此波長(zhǎng)處于近紅外光譜范圍內(nèi),說明在利用光譜特征預(yù)測(cè)松材線蟲病時(shí),近紅外波段應(yīng)作為重點(diǎn)考慮的波段。

    3)染病馬尾松蒸騰速率與6個(gè)光譜特征參數(shù)存在一定的相關(guān)性,其中與紅邊參數(shù)的相關(guān)性較高,R2均在0.7以上,與紅邊比值植被指數(shù)相關(guān)性最高,R2更是達(dá)到了0.851 5,說明光譜紅邊特征、特別是紅邊比值植被指數(shù)對(duì)松材線蟲病害有較好的指示作用。

    4)利用馬尾松蒸騰速率與其光譜特征的相關(guān)性來診斷病害,能有效克服由于潛伏期等原因?qū)е碌耐庥^特征不明顯的問題,為遙感監(jiān)測(cè)中特征參數(shù)的選取提供了參考。

    由于儀器、環(huán)境、人為等因素對(duì)試驗(yàn)結(jié)果有一定影響,如何更為合理地排除其他因素的干擾變得較為重要,這方面仍有待進(jìn)一步研究。在遙感應(yīng)用中,由于機(jī)載、星載等傳感器在監(jiān)測(cè)尺度上與地測(cè)光譜存在差異,本研究在“林木葉組分”范疇內(nèi)分析的蒸騰速率與光譜特征的相關(guān)性,僅為遙感應(yīng)用中光譜特征參數(shù)的選取提供理論參考,在后續(xù)的應(yīng)用中仍有待進(jìn)一步加以驗(yàn)證。

    [1] 黃明祥,龔建華,張健欽.松材線蟲病害遙感監(jiān)測(cè)與傳播模擬研究 [M].北京:中國環(huán)境科學(xué)出版社,2012:144.

    Huang M X,Gong J H,Zhang J Q.Study on pine wilt disease remote sensing monitoring and diffusion simulation [M].Beijing:China Environmental Science Press,2012:144.(in Chinese)

    [2] White J C,Wulder M A,Darin B,et al.Detection of red attack stage mountain pine beetle infestation with high spatial resolution satellite imagery [J].Remote Sensing Environment,2005,96:340-351.

    [3] Qin J,Burks T F,Ritenour M A,et al.Detection of citrus canker using hyperspectral reflectance imaging with spectral information divergence [J].Journal of Food Engineering,2009,93(2):183-191.

    [4] 張?jiān)龉?基于過程模型的松材線蟲災(zāi)害的遙感早期識(shí)別 [D].北京:北京林業(yè)大學(xué),2011:53.

    Zhang Z G.The early identification of remote sensing aboutBuraphelenchusxylophilusbased on process model [D].Beijing:Beijing Forestry University,2011:53.(in Chinese)

    [5] 陳福海.馬尾松毛蟲害遙感監(jiān)測(cè)及其影響因子分析 [D].福州:福建農(nóng)林大學(xué),2012:70.

    Chen F H.Dendrolimuspunctantusdamage remote sensing monitoring and influence factors analysis [D].Fuzhou:Fujian Agriculture and Forestry University,2012:70.(in Chinese)

    [6] 亓興蘭.SPOT-5遙感影像馬尾松毛蟲害信息提取技術(shù)研究 [D].福州:福建農(nóng)林大學(xué),2011:122.

    Qi X L.Study on information extraction technology ofDendrolimuspunctatusdamage based on SPOT-5 remote sensing images [D].Fuzhou:Fujian Agriculture and Forestry University,2011:122.(in Chinese)

    [7] 云麗麗,高國平,欒慶書,等.遼西油松毛蟲蟲害的TM影像監(jiān)測(cè)技術(shù) [J].安徽農(nóng)業(yè)科學(xué),2010,38(1):201-204.

    Yun L L,Gao G P,Luan Q S,et al.Monitoring technique for Chinese pine caterpillars by TM imagery in West Liaoning Province [J].Journal of Anhui Agriculture Science,2010,38(1):201-204.(in Chinese)

    [8] Franklin S E,Raske A G.Satellite remote sensing of spruce budworm forest defoliation in western New found land [J].Canadian Journal of Remote Sensing,1994,20(1):30-48.

    [9] 張 田.福建省三明市松材線蟲病遙感監(jiān)測(cè)預(yù)測(cè)研究 [D].北京:北京林業(yè)大學(xué),2010:63.

    Zhang T.Remote sensing monitoring and predicting ofBursaphelenchusxylophilusin Sanming of Fujian Province [D].Beijing:Beijing Forestry University,2010:63.(in Chinese)

    [10] Shafri H Z M,Hamdan N.Hyperspectral imagery for mapping disease infection in oil palm plantation using vegetation indices and red edge techniques [J].American Journal of Applied Sciences,2009,6(6):1031-1035.

    [11] Radeloff V C,Mladenoff D J,Boyce M S.Detecting jack pine budworm defoliation using spectral mixture analysis,separating effects from determinants [J].Remote Sensing Environment,1999,69:156-169.

    [12] 楊寶君,潘宏陽,湯 堅(jiān),等.松材線蟲病 [M].北京:中國林業(yè)出版社,2003:263.

    Yang B J,Pan H Y,Tang J,et al.Bursaphelenchusxylophilus[M].Beijing:China Forestry Publishing House,2003:263.(in Chinese)

    [13] 劉秀英,熊建利,臧 卓,等.馬尾松葉綠素含量與高光譜數(shù)據(jù)相關(guān)性分析 [J].廣東農(nóng)業(yè)科學(xué),2012(10):35-37,50.

    Liu X Y,Xiong J L,Zang Z,et al.Study on correlation between chlorophyll content with hyperspectral dataPinusmassonianaLamb [J].Guangdong Agricultural Sciences,2012(10):35-37,50.(in Chinese)

    [14] 杜華強(qiáng),葛宏立,范文義,等.馬尾松針葉光譜特征與其葉綠素含量間關(guān)系研究 [J].光譜學(xué)與光譜分析,2009(11):3033-3037.

    Du H Q,Ge H L,Fan W Y,et al.Study on relationships between total chlorophyll with hyperspectral features for leaves ofPinusmassonianaforest [J].Spectroscopy and Spectral Analysis,2009(11):3033-3037.(in Chinese)

    [15] 伍 南,劉君昂,閆瑞坤,等.利用光譜特征參數(shù)估算病害脅迫下杉木葉綠素含量 [J].植物保護(hù),2012(4):72-76,88.

    Wu N,Liu J A,Yan R K,et al.Using spectral feature parameters to estimate the chlorophyll content of Chinese fir under disease stress [J].Plant Protection,2012(4):72-76,88.(in Chinese)

    [16] 姚付啟,張振華,楊潤(rùn)亞,等.基于紅邊參數(shù)的植被葉綠素含量高光譜估算模型 [J].農(nóng)業(yè)工程學(xué)報(bào),2009(S2):123-129.

    Yao F Q,Zhang Z H,Yang R Y,et al.Hyperspectral models for estimating vegetation chlorophyll content based on red edge parameter [J].Transactions of CSAE,2009(S2):123-129.(in Chinese)

    [17] 梁守真,施 平,馬萬棟,等. 植被葉片光譜及紅邊特征與葉片生化組分關(guān)系的分析 [J].中國生態(tài)農(nóng)業(yè)學(xué)報(bào),2010(4):804-809.

    Liang S Z,Shi P,Ma W D,et al.Relational analysis of spectra and red-edge characteristics of plant leaf and leaf biochemical constituent [J].Chinese Journal of Eco-Agriculture,2010(4): 804-809.(in Chinese)

    [18] 王 潔,徐端松,馬躍良,等.植被含水量的遙感反演方法及研究進(jìn)展 [J].遙感信息,2008(1):100-105,90.

    Wang J,Xu D S,Ma Y L,et al.Methods and research developments for retrival of vegetable water content by remote sensing [J].Remote Sensing Information,2008(1):100-105,90.(in Chinese)

    [19] 于君明,王世新,周 藝,等.植被水分遙感監(jiān)測(cè)研究綜述 [J].遙感信息,2008(2):97-102.

    Yu J M,Wang S X,Zhou Y,et al.Review on monitoring models of vegetation water using remote sensing [J].Remote Sensing Information,2008(2):97-102.(in Chinese)

    [20] 馬 菁,張曉麗,劉 維,等.松材線蟲病對(duì)馬尾松蒸騰速率和光譜特征的影響 [J].東北林業(yè)大學(xué)學(xué)報(bào),2013(11):115-117.

    Ma J,Zhang X L,Liu W,et al.Influence ofBursaphelenchusxylophiluson transpiration rate and spectral characteristics ofPinusmassoniana[J].Journal of Northeast Forestry University,2013(11):115-117.(in Chinese)

    [21] 王 震,張曉麗,安樹杰.松材線蟲病危害的馬尾松林木光譜特征分析 [J].遙感技術(shù)與應(yīng)用,2007(3):367-370.

    Wang Z,Zhang X L,An S J.Spectral characteristics analysis ofPinusmassonianasuffered byBursaphelenchusxylophilus[J].Remote Sensing Technology and Application,2007(3):367-370.(in Chinese)

    [22] 詹遠(yuǎn)順.泰寧縣松材線蟲病發(fā)生現(xiàn)狀與防控對(duì)策 [J].生物災(zāi)害科學(xué),2014(1):91-94.

    Zhan Y S.Status quo ofBuraphelenchusxylophilus(Sterner et Buhrer) nickle occurrence in Taining and its prevention and control measures [J].Biological Disaster Science,2014(1):91-94.(in Chinese)

    [23] 徐華潮,駱有慶,張廷廷,等.松材線蟲自然侵染后松樹不同感病階段針葉光譜特征變化 [J].光譜學(xué)與光譜分析,2011(5):1352-1356.

    Xu H C,Luo Y Q,Zhang T T,et al.Changes of reflectance spectra of pine needles in different stage after being infected by pine wood nematode [J].Spectroscopy and Spectral Analysis,2011(5):1352-1356.(in Chinese)

    [24] 楊世杰.植物生物學(xué) [M].北京:高等教育出版社,2010:365.

    Yang S J.Plant biology [M].Beijing:Higher Education Press,2010:365.(in Chinese)

    [25] 劉 燦,李 宏,李丕軍.不同土壤水分對(duì)九里香光合特性的影響 [J].東北林業(yè)大學(xué)學(xué)報(bào),2009(9):61-67.

    Liu C,Li H,Li P J.Effect of photosynthetic characteristics ofMurrayaexoticaunder different soil water conditions [J].Journal of Northeast Forestry University,2009(9):61-67.(in Chinese)

    [26] 張向峰,王玉杰,李云霞,等.模擬酸雨對(duì)盆栽馬尾松葉片光合特性的影響 [J].水土保持學(xué)報(bào),2013(1):197-202.

    Zhang X F,Wang Y J,Li Y X,et al.Effect of simulated acid rain on photosynthetic characteristics of pottedPinusmassonian[J].Journal of Soil and Water Conservation,2013(1):197-202.(in Chinese)

    Prediction ofBursaphelenchusxylophilusbased on transpiration rate and spectral characteristics

    LIU Ning,ZHANG Xiao-li,WANG Shu-han,HUANG Ya-ting

    (CollegeofForestry/KeyLaboratoryforSilvicultureandConservationofMinistryofEducation,BeijingForestryUniversity,Beijing100083,China)

    【Objective】 This study investigated the changes in transpiration rate and spectral characteristics ofPinusmassonianaat the initial stage ofBursaphelenchusxylophilusinfection and their relationship to improve the identification ofBursaphelenchusxylophilusduring incubation and the selection of feature parameters in remote sensing monitoring.【Method】 Transpiration rates at four different disease periods (HP,EPD,MPD,and TPD) ofP.massonianainfected byB.xylophiluswere measured to investigate the change laws in Taining,Fujian.Then transpiration rates at HP and EPD were analyzed with collected spectral data and spectral characteristics (REP,RES,RBP,RRB,RERVI,and RENDVI).【Result】 The transpiration rate ofP.massonianagradually weakened with the deterioration of disease.There was certain relationship between the transpiration rates and raw spectral data at health and early disease stages.The correlation coefficients were dependent on wave length (the maximal value was 0.86 at 753 nm).Transpiration rate and red edge parameters had high relationship withR2of larger than 0.7,especially for RERVI withR2of 0.851 5.【Conclusion】 The transpiration rate ofP.massonianacould be an indicator of disease.Transpiration rate had high correlation with infrared band and red edge parameters.Thus spectral characteristics can be used for early forecast ofB.xylophilus.

    Pinusmassoniana;spectral characteristics;transpiration rate;Bursaphelenchusxylophilus;red edge parameters;disease forecast

    2014-08-15

    國家自然科學(xué)基金項(xiàng)目(30871964)

    劉 寧(1987-),男,吉林白城人,在讀碩士,主要從事3S技術(shù)在資源環(huán)境中的應(yīng)用研究。E-mail:103807711@qq.com

    張曉麗(1967-),女,河北張家口人,教授,博士生導(dǎo)師,主要從事林業(yè)遙感和地理信息系統(tǒng)研究。 E-mail:zhang-xl@263.net

    時(shí)間:2015-05-11 15:03

    10.13207/j.cnki.jnwafu.2015.06.011

    S431.9

    A

    1671-9387(2015)06-0129-07

    網(wǎng)絡(luò)出版地址:http://www.cnki.net/kcms/detail/61.1390.S.20150511.1503.011.html

    猜你喜歡
    針葉蒸騰速率松材
    土壤水分狀況對(duì)獼猴桃葉片蒸騰速率影響研究
    更 正
    黃土丘陵區(qū)山杏人工林蒸騰速率與環(huán)境因子的關(guān)系
    《松材線蟲病生態(tài)災(zāi)害督辦追責(zé)辦法》出臺(tái)
    綠色中國(2019年13期)2019-11-26 07:11:02
    風(fēng)雨海棠
    松材線蟲病的發(fā)生防治及對(duì)策——以重慶市為例
    我國松材線蟲病的主要防治技術(shù)探析
    今年1季度德國針葉材出口量增長(zhǎng)7%
    加拿大針葉鋸材繼續(xù)以美國南部為中心
    吉安市松材線蟲病發(fā)生規(guī)律初探
    国产中年淑女户外野战色| 国产视频一区二区在线看| 国产三级中文精品| 国产亚洲91精品色在线| 免费人成视频x8x8入口观看| 欧美高清性xxxxhd video| 中文资源天堂在线| 精品一区二区三区人妻视频| 国产蜜桃级精品一区二区三区| 在线观看美女被高潮喷水网站| 午夜激情福利司机影院| 久久久久免费精品人妻一区二区| xxxwww97欧美| 成人精品一区二区免费| 老熟妇仑乱视频hdxx| 精品乱码久久久久久99久播| 国产亚洲精品久久久com| 在线播放国产精品三级| 国产亚洲精品久久久久久毛片| 美女大奶头视频| 综合色av麻豆| 两人在一起打扑克的视频| 午夜视频国产福利| 成人国产一区最新在线观看| 极品教师在线免费播放| 日本撒尿小便嘘嘘汇集6| 国产激情偷乱视频一区二区| 又粗又爽又猛毛片免费看| 国产精品一区二区三区四区久久| 国产精品久久视频播放| 淫妇啪啪啪对白视频| 国产亚洲精品久久久久久毛片| 婷婷六月久久综合丁香| 日韩高清综合在线| 亚洲av成人精品一区久久| av天堂中文字幕网| www.色视频.com| 性插视频无遮挡在线免费观看| 女的被弄到高潮叫床怎么办 | 亚洲熟妇中文字幕五十中出| 乱系列少妇在线播放| 精品久久久久久久久久免费视频| 免费在线观看影片大全网站| 亚洲国产欧洲综合997久久,| 18禁黄网站禁片午夜丰满| 免费在线观看日本一区| 亚洲av中文av极速乱 | 一个人看视频在线观看www免费| 欧美日本亚洲视频在线播放| 日日干狠狠操夜夜爽| 啦啦啦啦在线视频资源| 国产午夜精品论理片| 人人妻,人人澡人人爽秒播| 色综合站精品国产| 国产69精品久久久久777片| 亚洲精华国产精华液的使用体验 | 国产熟女欧美一区二区| 国产亚洲91精品色在线| 真人一进一出gif抽搐免费| 最近最新中文字幕大全电影3| 亚洲人成网站在线播放欧美日韩| 深夜精品福利| 亚洲国产色片| 成人欧美大片| 人人妻人人澡欧美一区二区| 岛国在线免费视频观看| 精品久久久久久久久久久久久| 成年版毛片免费区| 97超视频在线观看视频| 国产美女午夜福利| 免费一级毛片在线播放高清视频| 午夜福利在线在线| 色播亚洲综合网| 精品久久久久久成人av| 亚洲国产高清在线一区二区三| 精品乱码久久久久久99久播| 日韩亚洲欧美综合| 亚洲国产日韩欧美精品在线观看| 午夜福利欧美成人| 高清毛片免费观看视频网站| 欧美最黄视频在线播放免费| 国产av在哪里看| 最近视频中文字幕2019在线8| 国产视频内射| 成年版毛片免费区| 久久精品影院6| 亚洲国产色片| 不卡视频在线观看欧美| 免费av毛片视频| 日韩欧美精品免费久久| 可以在线观看毛片的网站| 久久久久久九九精品二区国产| 色哟哟·www| 在线免费观看不下载黄p国产 | 性插视频无遮挡在线免费观看| 直男gayav资源| 变态另类成人亚洲欧美熟女| 99在线人妻在线中文字幕| 可以在线观看的亚洲视频| 国产伦精品一区二区三区四那| 日韩人妻高清精品专区| 一个人免费在线观看电影| 极品教师在线视频| 国产视频一区二区在线看| 亚洲自偷自拍三级| 国内精品久久久久久久电影| 伊人久久精品亚洲午夜| 久久亚洲精品不卡| 精品欧美国产一区二区三| 日韩一区二区视频免费看| 成人国产一区最新在线观看| 亚洲国产精品成人综合色| 99久久成人亚洲精品观看| 亚洲成人免费电影在线观看| 国国产精品蜜臀av免费| 极品教师在线视频| 亚洲中文字幕一区二区三区有码在线看| 麻豆国产97在线/欧美| 非洲黑人性xxxx精品又粗又长| 国产午夜精品论理片| 午夜福利欧美成人| 美女cb高潮喷水在线观看| 人妻丰满熟妇av一区二区三区| 国产高清视频在线播放一区| 超碰av人人做人人爽久久| 老司机福利观看| av天堂中文字幕网| 亚洲精品国产成人久久av| 欧美丝袜亚洲另类 | 桃红色精品国产亚洲av| 联通29元200g的流量卡| 亚洲成人精品中文字幕电影| 少妇猛男粗大的猛烈进出视频 | 婷婷丁香在线五月| 亚洲中文字幕一区二区三区有码在线看| 国产精品日韩av在线免费观看| 日本黄色片子视频| 又紧又爽又黄一区二区| 久久久精品欧美日韩精品| 免费看日本二区| 一个人看的www免费观看视频| 国产精品免费一区二区三区在线| 在线看三级毛片| 欧美日韩瑟瑟在线播放| 十八禁网站免费在线| 日韩欧美免费精品| av黄色大香蕉| 国产精品综合久久久久久久免费| av视频在线观看入口| 99久久成人亚洲精品观看| 12—13女人毛片做爰片一| 国产一区二区三区在线臀色熟女| 亚洲男人的天堂狠狠| 国产熟女欧美一区二区| 国产国拍精品亚洲av在线观看| 免费人成视频x8x8入口观看| 看片在线看免费视频| 亚洲人成网站在线播放欧美日韩| 亚洲熟妇中文字幕五十中出| 亚洲av电影不卡..在线观看| 91久久精品电影网| 一个人观看的视频www高清免费观看| 国产大屁股一区二区在线视频| 老熟妇仑乱视频hdxx| 少妇裸体淫交视频免费看高清| 别揉我奶头 嗯啊视频| 九九在线视频观看精品| 免费看美女性在线毛片视频| 又黄又爽又刺激的免费视频.| 亚洲 国产 在线| 一个人免费在线观看电影| 午夜久久久久精精品| 成人毛片a级毛片在线播放| 嫩草影院新地址| 又爽又黄a免费视频| 国产高潮美女av| 亚洲av一区综合| 亚洲不卡免费看| 亚洲一区二区三区色噜噜| 美女xxoo啪啪120秒动态图| 直男gayav资源| 日韩人妻高清精品专区| 亚洲五月天丁香| 亚洲三级黄色毛片| 日韩精品有码人妻一区| 国产成人aa在线观看| 一个人看视频在线观看www免费| 国产精品免费一区二区三区在线| 国产精品久久久久久久久免| 国产爱豆传媒在线观看| 一区二区三区四区激情视频 | 亚洲人成网站在线播| 1000部很黄的大片| 久久亚洲真实| 小蜜桃在线观看免费完整版高清| 久久久久久久久大av| 国产麻豆成人av免费视频| 人妻制服诱惑在线中文字幕| 免费av毛片视频| 热99re8久久精品国产| 成人国产综合亚洲| 国产av不卡久久| 亚洲精品456在线播放app | 精品无人区乱码1区二区| 99久久九九国产精品国产免费| 欧美黑人欧美精品刺激| 麻豆av噜噜一区二区三区| 99视频精品全部免费 在线| 久久久久久大精品| 真实男女啪啪啪动态图| 黄片wwwwww| 1024手机看黄色片| 国语自产精品视频在线第100页| 少妇的逼好多水| 欧美成人免费av一区二区三区| 1000部很黄的大片| 亚洲aⅴ乱码一区二区在线播放| 成人美女网站在线观看视频| 久久精品国产亚洲av香蕉五月| 18禁黄网站禁片午夜丰满| 看十八女毛片水多多多| 一个人免费在线观看电影| 人妻丰满熟妇av一区二区三区| 五月伊人婷婷丁香| bbb黄色大片| 国语自产精品视频在线第100页| 久久久久久久午夜电影| 老女人水多毛片| 国产午夜精品久久久久久一区二区三区 | 亚洲国产欧洲综合997久久,| 尤物成人国产欧美一区二区三区| eeuss影院久久| 午夜精品久久久久久毛片777| 国内精品一区二区在线观看| 欧美一区二区亚洲| 国产伦精品一区二区三区四那| 国产伦精品一区二区三区视频9| 亚洲人成伊人成综合网2020| 色哟哟哟哟哟哟| 91在线精品国自产拍蜜月| 成人综合一区亚洲| 他把我摸到了高潮在线观看| 欧美极品一区二区三区四区| 午夜影院日韩av| 日韩亚洲欧美综合| 中文亚洲av片在线观看爽| 校园人妻丝袜中文字幕| 精品国内亚洲2022精品成人| 久久精品夜夜夜夜夜久久蜜豆| 美女 人体艺术 gogo| 亚洲人成网站在线播放欧美日韩| 久久亚洲真实| 欧美激情久久久久久爽电影| 国产精品精品国产色婷婷| 久久国产精品人妻蜜桃| 国产免费av片在线观看野外av| 97超级碰碰碰精品色视频在线观看| 国产熟女欧美一区二区| 亚洲五月天丁香| 国产精品亚洲美女久久久| 九色成人免费人妻av| 亚洲不卡免费看| 久久人人精品亚洲av| 99在线视频只有这里精品首页| 老熟妇乱子伦视频在线观看| 岛国在线免费视频观看| 欧美人与善性xxx| 99在线人妻在线中文字幕| 久久精品国产鲁丝片午夜精品 | 国产av一区在线观看免费| 免费看a级黄色片| 久久午夜亚洲精品久久| 精品福利观看| 欧美成人性av电影在线观看| 国产成人一区二区在线| 国国产精品蜜臀av免费| 99精品在免费线老司机午夜| 国产亚洲91精品色在线| 午夜久久久久精精品| 一级a爱片免费观看的视频| 三级毛片av免费| 夜夜看夜夜爽夜夜摸| 免费人成视频x8x8入口观看| 色5月婷婷丁香| 哪里可以看免费的av片| 久久久久久久亚洲中文字幕| 国产不卡一卡二| 热99re8久久精品国产| 亚洲一区二区三区色噜噜| 国产精品久久视频播放| 国产精品98久久久久久宅男小说| 亚洲av熟女| 在线播放无遮挡| 国产成人福利小说| 99精品在免费线老司机午夜| 99在线人妻在线中文字幕| 国产又黄又爽又无遮挡在线| www.色视频.com| 成熟少妇高潮喷水视频| 黄片wwwwww| 国产69精品久久久久777片| 亚洲专区中文字幕在线| 麻豆成人午夜福利视频| 亚洲中文字幕一区二区三区有码在线看| 亚洲五月天丁香| 少妇熟女aⅴ在线视频| 舔av片在线| 欧美最黄视频在线播放免费| av视频在线观看入口| 国产亚洲精品av在线| 国产一区二区激情短视频| 国产aⅴ精品一区二区三区波| 在线观看一区二区三区| 女的被弄到高潮叫床怎么办 | 尤物成人国产欧美一区二区三区| 久久精品久久久久久噜噜老黄 | 色噜噜av男人的天堂激情| 国模一区二区三区四区视频| 久久久久久久久久成人| 免费高清视频大片| 精品日产1卡2卡| 免费在线观看影片大全网站| 99精品久久久久人妻精品| 午夜福利在线观看免费完整高清在 | 婷婷精品国产亚洲av在线| 国产真实乱freesex| 国产午夜福利久久久久久| 在线观看午夜福利视频| 在线a可以看的网站| 一个人观看的视频www高清免费观看| 老司机深夜福利视频在线观看| 我要看日韩黄色一级片| av在线亚洲专区| 欧美一区二区国产精品久久精品| 99久久精品一区二区三区| 亚洲自偷自拍三级| 人妻丰满熟妇av一区二区三区| 日韩欧美三级三区| 欧美一级a爱片免费观看看| 亚洲综合色惰| 亚洲天堂国产精品一区在线| 香蕉av资源在线| 亚洲不卡免费看| avwww免费| 91精品国产九色| 国产精品福利在线免费观看| 成年免费大片在线观看| 精品久久久久久,| 欧美xxxx性猛交bbbb| 悠悠久久av| 国产亚洲av嫩草精品影院| 欧美日韩综合久久久久久 | 国产探花极品一区二区| 最近中文字幕高清免费大全6 | 丰满乱子伦码专区| 少妇的逼水好多| 国产精品三级大全| 看黄色毛片网站| 女生性感内裤真人,穿戴方法视频| 中文字幕久久专区| 亚洲内射少妇av| 国产精品电影一区二区三区| 欧美日韩中文字幕国产精品一区二区三区| 午夜日韩欧美国产| or卡值多少钱| 欧美日韩国产亚洲二区| 少妇裸体淫交视频免费看高清| 国产午夜精品论理片| 欧美一区二区亚洲| 看免费成人av毛片| 欧美成人a在线观看| 一个人看的www免费观看视频| 啪啪无遮挡十八禁网站| 女同久久另类99精品国产91| 午夜免费激情av| 日本黄色视频三级网站网址| 日韩亚洲欧美综合| 在线观看av片永久免费下载| xxxwww97欧美| 无遮挡黄片免费观看| 美女高潮的动态| 自拍偷自拍亚洲精品老妇| 桃色一区二区三区在线观看| 91av网一区二区| 91麻豆av在线| 国产免费av片在线观看野外av| 黄色丝袜av网址大全| 精品一区二区免费观看| 99热这里只有是精品50| 有码 亚洲区| a级一级毛片免费在线观看| 亚洲国产欧美人成| 国产精品人妻久久久久久| 狂野欧美激情性xxxx在线观看| 中文字幕av在线有码专区| 丰满乱子伦码专区| 久久久久久久久久成人| 少妇猛男粗大的猛烈进出视频 | 两性午夜刺激爽爽歪歪视频在线观看| 女人十人毛片免费观看3o分钟| 少妇熟女aⅴ在线视频| 欧美在线一区亚洲| 久久香蕉精品热| 免费观看的影片在线观看| 在线a可以看的网站| 嫁个100分男人电影在线观看| 国产精品1区2区在线观看.| 在线观看午夜福利视频| av在线蜜桃| av女优亚洲男人天堂| 一进一出抽搐gif免费好疼| 国产成人aa在线观看| 午夜免费成人在线视频| 老司机福利观看| 国产伦精品一区二区三区视频9| 天堂影院成人在线观看| 变态另类丝袜制服| 欧美日本视频| 岛国在线免费视频观看| 久久香蕉精品热| 1024手机看黄色片| 午夜免费成人在线视频| 男人狂女人下面高潮的视频| 久久精品夜夜夜夜夜久久蜜豆| 尤物成人国产欧美一区二区三区| 亚洲av免费在线观看| 99精品在免费线老司机午夜| 国产精品精品国产色婷婷| 天天一区二区日本电影三级| 亚洲精品久久国产高清桃花| 国产一级毛片七仙女欲春2| 麻豆成人午夜福利视频| 日本与韩国留学比较| 免费看日本二区| 在线观看一区二区三区| 又黄又爽又刺激的免费视频.| 国产av不卡久久| 午夜日韩欧美国产| 欧美激情在线99| 国产精品亚洲一级av第二区| 麻豆国产av国片精品| 波多野结衣高清无吗| 成人二区视频| 99热网站在线观看| 可以在线观看毛片的网站| 桃红色精品国产亚洲av| 欧美三级亚洲精品| 精品无人区乱码1区二区| 国产三级在线视频| 精品久久久久久成人av| 天堂动漫精品| 午夜日韩欧美国产| 91久久精品国产一区二区成人| 亚洲人成网站在线播放欧美日韩| 可以在线观看的亚洲视频| 欧美日本亚洲视频在线播放| 精品99又大又爽又粗少妇毛片 | 成人精品一区二区免费| 狠狠狠狠99中文字幕| 小说图片视频综合网站| 一进一出抽搐动态| 亚洲自拍偷在线| 老司机福利观看| 久久99热这里只有精品18| 九色国产91popny在线| 天堂动漫精品| 国产成人一区二区在线| 国模一区二区三区四区视频| 两个人的视频大全免费| 成人二区视频| 一进一出好大好爽视频| 国产在线男女| 国产高清有码在线观看视频| 国产亚洲精品久久久久久毛片| 午夜免费激情av| 日韩一本色道免费dvd| 国产69精品久久久久777片| 草草在线视频免费看| 日韩大尺度精品在线看网址| 亚洲av一区综合| 变态另类成人亚洲欧美熟女| 可以在线观看的亚洲视频| 美女大奶头视频| 中文在线观看免费www的网站| 成年女人永久免费观看视频| 两性午夜刺激爽爽歪歪视频在线观看| 嫁个100分男人电影在线观看| 欧美+亚洲+日韩+国产| 极品教师在线视频| 美女高潮喷水抽搐中文字幕| 国产亚洲精品av在线| 欧美黑人巨大hd| 我的老师免费观看完整版| 日本撒尿小便嘘嘘汇集6| 村上凉子中文字幕在线| 亚洲第一电影网av| 精品久久久久久久久av| 亚洲无线观看免费| 午夜福利在线观看吧| 97超视频在线观看视频| 小蜜桃在线观看免费完整版高清| 欧美在线一区亚洲| 日韩欧美 国产精品| 99热这里只有是精品在线观看| 男女做爰动态图高潮gif福利片| 免费在线观看日本一区| 国产高清有码在线观看视频| 日本黄大片高清| 日本撒尿小便嘘嘘汇集6| 久久99热这里只有精品18| 日韩一区二区视频免费看| 亚洲七黄色美女视频| 日本-黄色视频高清免费观看| 国产视频内射| 我的女老师完整版在线观看| 韩国av一区二区三区四区| 一个人观看的视频www高清免费观看| 两个人的视频大全免费| 特大巨黑吊av在线直播| 久久人人爽人人爽人人片va| 熟妇人妻久久中文字幕3abv| 免费观看的影片在线观看| 日本与韩国留学比较| 老师上课跳d突然被开到最大视频| or卡值多少钱| 欧美区成人在线视频| 欧美xxxx性猛交bbbb| 国产主播在线观看一区二区| 亚洲av中文字字幕乱码综合| 欧美成人a在线观看| 日韩人妻高清精品专区| 亚洲精品乱码久久久v下载方式| 男女边吃奶边做爰视频| 国产乱人伦免费视频| 精品人妻熟女av久视频| 色噜噜av男人的天堂激情| 啦啦啦韩国在线观看视频| 一区二区三区四区激情视频 | 网址你懂的国产日韩在线| 国产成人aa在线观看| 精品国产三级普通话版| 成人精品一区二区免费| videossex国产| 亚洲精品在线观看二区| 在线播放国产精品三级| 九九热线精品视视频播放| av福利片在线观看| 午夜激情欧美在线| 国产精品久久久久久精品电影| 99精品在免费线老司机午夜| 日韩av在线大香蕉| 男女啪啪激烈高潮av片| 日本三级黄在线观看| 琪琪午夜伦伦电影理论片6080| 亚洲av二区三区四区| 18禁黄网站禁片午夜丰满| 草草在线视频免费看| 18+在线观看网站| 久久午夜福利片| 免费看a级黄色片| 又黄又爽又刺激的免费视频.| 国产伦人伦偷精品视频| 国产一级毛片七仙女欲春2| 婷婷丁香在线五月| 色哟哟·www| av福利片在线观看| 夜夜夜夜夜久久久久| 国产成年人精品一区二区| 特级一级黄色大片| 精品一区二区三区视频在线| 少妇丰满av| 中文字幕高清在线视频| 日韩欧美精品免费久久| 国产精品无大码| 国产探花在线观看一区二区| 久久亚洲真实| av在线观看视频网站免费| 99热这里只有是精品50| 欧美激情在线99| 国产三级中文精品| 尾随美女入室| 性色avwww在线观看| 两人在一起打扑克的视频| 日本熟妇午夜| 免费搜索国产男女视频| 一区二区三区激情视频| aaaaa片日本免费| 少妇被粗大猛烈的视频| 亚洲avbb在线观看| 久久精品综合一区二区三区| 免费看光身美女| 久9热在线精品视频| 欧美bdsm另类| 成人欧美大片| 69人妻影院| 国产女主播在线喷水免费视频网站 | 午夜免费激情av| av中文乱码字幕在线| 亚洲不卡免费看| 国产一区二区亚洲精品在线观看| 国产精品1区2区在线观看.| 51国产日韩欧美| 久久久久久伊人网av| 精品久久久久久久人妻蜜臀av| 18禁裸乳无遮挡免费网站照片| 无遮挡黄片免费观看| 免费观看精品视频网站| 免费搜索国产男女视频| 日本与韩国留学比较| 亚洲人成网站高清观看| 亚洲av第一区精品v没综合| ponron亚洲| av天堂在线播放| 禁无遮挡网站| 69av精品久久久久久|