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

    基于高光譜數(shù)據(jù)的楊樹葉片干物質含量的估算*

    2017-02-10 03:11:41程志慶張勁松李巖泉
    中國農業(yè)氣象 2017年1期
    關鍵詞:冠層楊樹波段

    程志慶,張勁松**,孟 平,李巖泉,鄭 寧

    ?

    基于高光譜數(shù)據(jù)的楊樹葉片干物質含量的估算*

    程志慶1,2,3,張勁松1,2,3**,孟 平1,李巖泉1,鄭 寧1

    (1.中國林業(yè)科學研究院林業(yè)研究所,北京 100091;2.南京林業(yè)大學南方現(xiàn)代林業(yè)協(xié)同創(chuàng)新中心,南京 210037;3.國家林業(yè)局林木培育重點實驗室,北京 100091)

    將楊樹葉片實測的理化數(shù)據(jù)和土壤背景光譜數(shù)據(jù)作為PROSPECT和PROSAIL模型的輸入?yún)?shù),輸出楊樹葉片高光譜數(shù)據(jù)模擬值,通過實測獲得葉片尺度和冠層尺度干物質含量、等效水厚度以及高光譜數(shù)據(jù),利用統(tǒng)計方法,分別對兩種尺度楊樹葉片干物質含量進行分析。結果表明:基于歸一化指數(shù)計算方法,楊樹葉片尺度和冠層尺度的最佳估算干物質含量的干物質含量歸一化指數(shù)(NDMI)波段組合分別為1685,1704nm和1551,2143nm,使用偏最小二乘法分別對篩選波段組成的NDMI(1685,1704)指數(shù)和NDMI(1551,2143)指數(shù)構建葉片尺度干物質含量和冠層尺度干物質含量的估算模型,葉片尺度干物質含量估算模型精度為R2=0.663,RMSE=0.001g·cm-2,冠層尺度精度為R2=0.91,RMSE=16.7g·m-2??梢?,高光譜技術對楊樹葉片干物質含量的估算具有較高的精度,可為楊樹葉片干物質含量的快速、無損估算提供參考依據(jù)。

    高光譜技術;楊樹;干物質含量;估算模型

    葉片干物質含量是葉片的重要形狀參數(shù)之一,能夠反映葉片對光合有效輻射的捕獲能力,而且與抗逆性、葉片壽命等具有很大的相關性[1-2]。楊樹因其生長迅速、成林早、木材用途廣、防護效益高等優(yōu)勢,是中國人工林的主要樹種。因此,快速、準確測定楊樹葉片干物質含量,對于監(jiān)測楊樹生長狀態(tài),加強防護林的經營管理,提升楊樹人工林的生態(tài)功能具有重要意義。

    干物質含量在植物生態(tài)學研究中有3種表示:葉片干重和鮮重比、比葉面積以及單位葉面積的干物質含量。由于后兩種葉片的干物質含量表示形式與葉片結構參數(shù)以及植物的生長和生存對策有緊密聯(lián)系,能反映植物對不同生境的適應性,使其成為植物生態(tài)學研究中的首選指標,也是輻射傳輸模型的重要參數(shù)[3-6]。目前,葉片干物質含量的測量主要是利用稱重測量葉面積的方法進行計算[7]。盡管該方法較簡便,但要獲取大面積干物質含量時,該方法耗時耗力,破壞性大。因此,需要一種新的技術對其進行快速的測定。而隨著遙感技術的發(fā)展,特別是當前的高光譜遙感技術,為葉片干物質含量的快速估算提供了巨大的潛力。Ria?o等[8]利用光譜數(shù)據(jù)分別對鮮葉的干物質含量與干葉的干物質含量進行估算,發(fā)現(xiàn)干葉的干物質含量估算較鮮葉估算精度高。Maire[9]利用指數(shù)的方法估算葉片尺度與冠層尺度干物質含量,發(fā)現(xiàn)該指數(shù)對葉片尺度干物質含量估算精度不高。Broge等[10]研究發(fā)現(xiàn)干物質含量與Landsat TM數(shù)據(jù)中的紅光、紅外以及中紅外波段具有顯著的相關性,并認為利用經驗模型可以在光譜分辨率的Landsat TM數(shù)據(jù)中估算冠層的干物質含量。從現(xiàn)有研究可見,至今還未見報道一種植被指數(shù)或者干物質反演模型能夠高精度地估算任何植物干物質含量。

    為了對楊樹葉片干物質含量進行快速、準確的估算,本研究通過分析已有葉片干物質含量植被指數(shù)對葉片干物質含量的估算精度,根據(jù)楊樹干物質含量與不同波段之間的相關性,提取楊樹干物質含量特征波段,并建立新植被指數(shù)。在建立葉片尺度的干物質含量高光譜估算模型的基礎上,進一步構建葉片冠層尺度干物質含量的高光譜估算模型,獲得楊樹干物質含量估算精度最高的植被指數(shù)與反演模型,以期為楊樹葉片干物質含量的準確、快速估算提供理論支持。

    1 材料與方法

    1.1 樣地概況

    本研究分別設置楊樹盆栽實驗與楊樹人工林野外試驗。

    盆栽實驗在河北農業(yè)大學科技園(115°24'31.31"E,38°48'22.31"N)進行。實驗材料采用盆栽一年生,形態(tài)相近的107速生楊樹苗。為了保證結果的全面性與精確性,獲得不同生長狀態(tài)的楊樹,共設6個處理,分別是極度干旱、不施氮處理,干旱、不施氮處理,正常水分管理、不施氮處理(即為空白對照);極度干旱、施氮處理,干旱、施氮處理,正常水分、施氮處理。水分處理條件為:最大田間持水量的70%為正常、55%為干旱和40%為極度干旱;氮素處理條件為:每株每千克土壤施氮量為0.167g。

    楊樹人工林位于衡水市武邑縣(38°49′N,115°53′E)。樹高10m,株行距為2m×3m,面積約30000m2,林下少有植被。在林地中間設立一座25m觀測鐵塔,在鐵塔上于楊樹冠層底部、內部以及頂部1m處各修建觀測平臺。

    1.2 實測數(shù)據(jù)

    1.2.1 光譜數(shù)據(jù)測定

    (1)葉片光譜使用Field Spec Pro光譜儀(光譜波段范圍為350?2500nm,美國產)自帶的植被探頭與葉片夾式光譜探測器在野外直接測定。 楊樹葉片數(shù)據(jù)采樣時間為2014年6月10?12日、 7月8?10日、8月10?13日和9月10?12日的每日10:00?14:00。為減少誤差,測量前用標準白板進行校正,同時將葉片表層的浮塵用干凈干燥的試紙輕拂去。每個葉片觀測記錄6次采集數(shù)據(jù)(測定時避開葉脈),取其平均值作為該處的光譜反射率。利用植被探頭與葉片夾,采用光譜儀自帶光源,保證葉片平整且被探測面積相同,有效消除了背景反射、葉片表面彎曲導致的光譜波動的影響,保證了研究結果的精確性,共獲得304組葉片尺度高光譜實測數(shù)據(jù)。

    (2)冠層光譜。利用25m觀測塔上位于樹頂上方2m處的觀測平臺,在天氣晴朗無云、無風的2014年10月5日10:00?14:00,光纖視場角為25°的Field Spec Pro光譜儀(美國產),從平臺上垂直朝下測量楊樹冠層反射光譜。并利用LI-2200測定該視場角范圍內(0.05m2)的冠層葉面積指數(shù)(LAI,圖1)。獲取冠層光譜的同時,利用高光譜儀在楊樹林區(qū)的空地處選擇與測定處一致的土壤進行土壤反射率測定。共獲得9組冠層尺度高光譜實測數(shù)據(jù)。

    1.2.2 葉片干物質含量、等效水厚度的測定

    葉片采集。將測定光譜后的楊樹葉片,摘取后立刻放入自封袋,并冷藏;冠層葉片等效水厚度采樣,采取光譜測定視場角范圍內楊樹頂層、中層、底層三部位葉片各20片,放入自封袋,并冷藏。

    葉片等效水厚度測定。首先利用Li-3000測定采集鮮葉的葉面積,并稱重,然后立刻放入105℃烘箱中15min殺青,再將溫度調至90℃烘干至恒重。

    葉片等效水厚度EWTleaf的計算

    葉片尺度干物質含量的計算

    冠層尺度干物質含量的計算

    式中,LM和CM分別表示葉片干物質含量(g·cm?2)和冠層尺度干物質含量(g·m?2),Dw為葉片干重(g),F(xiàn)W為葉片鮮葉重量(g),A為葉片面積(cm2),LAI為葉面積指數(shù)(m2·m?2)。

    1.3 輻射傳輸模型參數(shù)

    植物葉片光譜分析時需要大量光譜數(shù)據(jù),而且在野外進行葉片光譜數(shù)據(jù)采集時,其它物質光譜的干擾造成葉片光譜有用信息的提取與分析更為困難。因此,可以根據(jù)植被的參數(shù)變化范圍利用輻射傳輸模型獲取模擬的光譜數(shù)據(jù)。

    PROSPECT模型[6]是一個基于“平板模型”發(fā)展而來的葉片輻射傳輸模型,模型中的輸入主要參數(shù)包括葉肉結構參數(shù)(N)、葉片葉綠素含量、葉片等效水厚度(EWT)以及葉片干物質含量(LM),能夠模擬葉片400?2500nm的半球反射和透過率。

    由Verhoef[11]在Suit模型的基礎上發(fā)展起來的SAIL模型是使用最為廣泛的冠層二向性反射物理模型之一,該模型描述了水平均勻植被冠層中直射和上行下行散射光通量的輻射傳輸過程,通過輸入葉片的反射率和透過率、葉面積指數(shù)、葉傾角參數(shù)、土壤反射率、太陽輻射的散射分量以及光線入射和觀測的幾何參數(shù),可以實現(xiàn)對植被冠層的反射光譜的模擬。

    PROSPECT模型能夠輸出葉片的反射率和透過率,而SAIL模型中又需要這些作為其參數(shù)進行冠層光譜的擬合,因此,研究者將兩種模型進行耦合形成PROSAIL模型,該模型能夠在其它冠層參數(shù)已知的條件下對植被的化學參數(shù)進行估算。PROSAIL模型最后所需的模型參數(shù)為葉片理化學參數(shù)(葉肉結構參數(shù)、葉綠素含量、等效水厚度、葉片干物質含量),冠層結構參數(shù)(葉面積指數(shù)、葉傾角分布類型)和其它參數(shù)(熱點效應參數(shù)、太陽天頂角、觀測天頂角、相對方位角、土壤反射率)。

    本文主要考慮為了獲取楊樹較為相近的光譜數(shù)據(jù),利用實測楊樹葉片生理參數(shù)數(shù)據(jù)范圍和土壤反射率作為模型輸入?yún)?shù),PROSAIL模型輸入的葉片理化參數(shù)與PROSPECT模型相同,具體模型參數(shù)設置見表1。

    2 結果與分析

    2.1 基于敏感波段的歸一化干物質含量指數(shù)構建

    2.1.1 干物質含量對光譜反射率的影響

    根據(jù)楊樹葉片等效水厚度的實測值將等效水厚度設為0.06cm,葉片干物質含量依據(jù)表1中范圍變化,由PROSPECT模型分別模擬的葉片反射光譜曲線見圖2a;葉片干物質含量和等效水厚度同時按表1中范圍變化,由PROSPECT模型模擬的葉片反射光譜曲線見圖2b。

    由圖2a可見,葉片干物質含量影響的波段主要集中于近紅外波段(750?1350nm)以及短波紅外的兩個波段范圍(1500?1870nm和2000?2500nm),在可見光范圍(400?700nm)的反射率光譜差異較小。圖2b為干物質含量逐漸增大而等效水厚度逐漸減小的反射光譜曲線圖,當葉片中等效水厚度變化與干物質含量相當時,由于受到葉片等效水厚度的影響,對干物質含量敏感的短波紅外波段的敏感性明顯減弱。因此,認為對干物質含量敏感的波段范圍為近紅外波段以及部分短波紅外波段(2000?2250nm)。

    表1 模型參數(shù)

    注:圖a中等效水厚度為0.06cm;圖b中干物質含量和不同等效水厚度均按照表1中變化

    Note: In fig.a, EWT is 0.06cm;In fig.b, EWT and LM change according to table 1

    2.1.2 葉片干物質含量的特征波段提取

    采用實測數(shù)據(jù)與PROSPECT模型擬合數(shù)據(jù),分析葉片干物質特征波段,擬合數(shù)據(jù)采用2.1.1節(jié)中擬合的光譜數(shù)據(jù),實測數(shù)據(jù)為153組實測楊樹葉片干物質含量。利用實測楊樹光譜數(shù)據(jù)和PROSPECT模型擬合的光譜數(shù)據(jù)分別與葉片干物質含量作相關性分析,結果如圖3。

    由圖3可見,在750?1350nm近紅外波段兩種數(shù)據(jù)與干物質含量的相關性均較大,但卻呈現(xiàn)出相反的相關方向,圖2a在該波段區(qū)間內的光譜反射率隨LM增大而減少,根據(jù)干物質的光學特征分析,該波段區(qū)間應呈現(xiàn)負相關,但實測結果卻顯示相反的趨勢。從圖2b的不同干物質含量和等效水厚度光譜圖以及圖4干物質含量與等效水厚度的關系圖可知,在楊樹葉片干物質含量較小、等效水厚度較大的條件下,750?1350nm近紅外波段的反射率隨干物質的增大而增大。從圖3可見,所獲得的楊樹葉片的等效水厚度較大,干物質含量較小。由于干物質含量在所有的特征波段中都受到其它因素的影響,因此,在干物質含量進行特征波段選取時,應根據(jù)實際光譜信息進行分析提取。同時從圖3中可見,實測光譜波段與干物質含量的相關系數(shù)均較低。因此,為獲取更好的估算波段,對兩波段組合的歸一化干物質含量指數(shù)(NDMI)對干物質含量的估算能力進行比較。NDMI的計算式為

    式中,R(i)與R(j)分別表示為第i波段的反射率和第j波段的反射率。

    由圖5可見,決定系數(shù)較大的區(qū)間(R2>0.62)主要分布在1250?2250nm,黑色圓圈內表示決定系數(shù)最大(R2=0.696)的區(qū)域,即以(1441nm,2309nm)坐標點為中心的附近,可見以上隨機選取兩波段的NDMI對葉片干物質含量具有很高的敏感性。但是,大氣水分在1400nm附近產生較強的吸收作用,地面遙感由于距離近,大氣水分對其影響較小,可以忽略,而衛(wèi)星遙感獲取地表植被光譜時,需要透過大氣,所以在利用衛(wèi)星遙感數(shù)據(jù)時應避免選取非大氣窗口波段[12]。因此,本研究選取非1400nm附近的最大系數(shù)的波段組合,即圖5中紅色圈位置,其坐標點為(1685nm,1704nm),其決定系數(shù)為0.661。從圖3中可以看出,這兩個波段內無論是擬合數(shù)據(jù)還是實測數(shù)據(jù),其與干物質均具有較大的負相關性,同時1704nm是最大的負相關波段,選取的這兩個波段,在大氣中的透射率為90%,能夠避免大氣對其它波段的吸收而產生干擾,故利用NDMI進行衛(wèi)星遙感估算楊樹葉片干物質含量時可使用(1685nm,1704nm)波段組合。而為獲得能在LandSAT 8等寬波段光譜數(shù)據(jù)中應用的干物質含量植被指數(shù),選取波段范圍應為LandSAT 8所包含的波段。由此,本研究選取最大決定系數(shù)的波段組合為(1551nm,2143nm)。

    2.2 葉片尺度干物質含量光譜估算模型的構建

    為了實現(xiàn)葉片尺度NDMI指數(shù)的干物質含量高精度估算,在304組實測數(shù)據(jù)中以153組數(shù)據(jù)進行建模,另外151組數(shù)據(jù)驗證模型估算精度。

    利用NDMI指數(shù)進行干物質含量LM估算的關系式為

    式中,K和b為估算模型的線性參數(shù),由選定NDMI波段與實測干物質含量求出。

    由表2可見,NDMI(1441,2309)和NDMI(1685,1704)干物質含量估算模型的殘差相同,為0.0001,兩者的擬合決定系數(shù)均為極顯著,但是前者擬合決定系數(shù)比后者大0.036??梢?,在葉片尺度的光譜估算中NDMI(1441,2309)具有一定優(yōu)勢。

    表2 葉片尺度優(yōu)選的3種波段組合下干物質含量的NDMI估算模型

    注:*、**分別表示相關系數(shù)(R)通過0.05、0.01水平的顯著性檢驗。下同。

    Note:*is P<0.05,**is P<0.01.The same as below.

    從圖6可以看出,3種兩波段組合的NDMI估算模型的估算精度,雖然NDMI(1441,2309)模型估算結果擬合度略高于NDMI(1685,1704),但后者擬合結果更接近1:1線,加之兩個估算結果的RMSE相近。因此,考慮到大氣水汽吸收時,NDMI(1685,1704)應為最佳選擇指數(shù)。NDMI(1551,2143)模型估算精度雖然較前者低,但其均方根誤差及斜率均較可觀,關鍵是可以應用于Landsat 8等衛(wèi)星寬波段數(shù)據(jù)的NDMI構建,所以NDMI(1551,2143)指數(shù)可以作為衛(wèi)星光譜數(shù)據(jù)估算干物質含量的最佳指數(shù)。

    2.3 冠層尺度干物質含量光譜模型的構建

    由式(3)可見,冠層尺度的干物質含量除受葉片尺度干物質含量的影響外,冠層葉面積指數(shù)也會對其產生影響。因此,利用PROSAIL模型擬合的冠層光譜數(shù)據(jù)對冠層葉面積和干物質含量的影響進行分析。

    由圖7可見,3種NDMI指數(shù)由于受LAI的影響,在相同葉面積指數(shù)條件下,這3種指數(shù)與冠層尺度干物質含量均具有顯著相關性,其關系表達式與葉片尺度下雖然有差異,但其表達式形式相同。由于不同葉面積指數(shù)的影響表達式的系數(shù)產生了變化,為此,設想利用LAI對上述模型的系數(shù)進行修正,以分析LAI與模型表達式中的系數(shù)關系(圖8)。

    在鉆進機構抵達限幅機構預定位置時,對接鎖合組件能夠實現(xiàn)限幅機構與鉆進機構的鎖定,從而使得限幅機構與鉆進機構進行隨動,保證后續(xù)的鉆取采樣作業(yè)順利進行。

    從圖8可以看出,3種NDMI估算CM模型系數(shù)a、b、c與LAI的擬合決定系數(shù)R2均大于0.99(P<0.001),因此,可以利用LAI替換模型系數(shù),得到利用LAI改進的3種NDMI估算冠層尺度干物質含量的模型。

    利用9組楊樹冠層實測數(shù)據(jù)對冠層尺度干物質含量模型進行驗證,結果見圖9。由圖可見,利用LAI修正后的估算模型NDMI(1441,2309)模型估算精度最高,而NDMI(1551,2143)估算值與實測值的擬合決定系數(shù)比NDMI(1685,1704)模型的高出0.2,且NDMI(1551,2143)模型估算的結果與實測值更接近1:1的擬合線。因此,結合2.2節(jié)中的NDMI(1551,2143)與NDMI(1551,2143)指數(shù)分析可知,利用地面高光譜數(shù)據(jù)進行楊樹冠層尺度干物質含量估算時,NDMI(1441,2309)構建的模型最優(yōu),而利用衛(wèi)星遙感數(shù)據(jù)對冠層尺度的干物質含量進行估算時,NDMI(1551,2143)構建的模型為最佳選擇。

    3 結論與討論

    3.1 結論

    通過對楊樹高光譜數(shù)據(jù)的實測和利用PROSPECT與PROSAIL模型分別進行葉片尺度和冠層尺度的干物質含量分析,利用歸一化指數(shù)的計算方法,篩選出楊樹葉片尺度和冠層尺度的最佳NDMI指數(shù)波段組合,根據(jù)所選對應葉片尺度和冠層尺度的高光譜特征波段組合的NDMI指數(shù),構建了楊樹葉片尺度及冠層尺度干物質含量模型。對兩種不同尺度的干物質含量模型的檢驗與驗證表明,所構建的估算模型具有較高精度。

    3.2 討論

    (1)由于模型模擬光譜數(shù)據(jù)只需葉片的部分生理生化參數(shù),排除了非葉片本身所引起的光譜變化,因此,在進行特征波段提取時具有很大的優(yōu)勢。但是由于野外環(huán)境復雜多變,若簡單地利用模型模擬光譜數(shù)據(jù)進行特征波段的提取會讓所選特征波段的適用性降低[13-14],因此,利用模型模擬光譜數(shù)據(jù)的同時應結合實測數(shù)據(jù),并利用實測數(shù)據(jù)進行驗證,最終才能提取適用的光譜特征波段。本研究利用輻射傳輸模型模擬的光譜數(shù)據(jù)與實測數(shù)據(jù)相結合的辦法,最終提取了分別對楊樹葉片尺度和冠層尺度葉片干物質含量準確估算的特征波段。

    (2)本試驗在研究冠層尺度干物質含量的估算中,通過對葉片尺度獲取的植被指數(shù)中引入葉面積,以提高冠層尺度干物質含量的估算精度,該法獲得了較高精度的冠層尺度干物質含量的估算模型。程曉娟等[15]在研究冠層尺度等效水厚度時指出,利用植被指數(shù)對冠層尺度等效水厚度進行估算,由于直接利用冠層等效水厚度選取的植被指數(shù),沒有從等效水厚度自身考慮,缺乏理論基礎,從而影響模型的準確性。野外中冠層尺度干物質含量的高光譜不僅受葉面積的影響,測定的冠層尺度高光譜數(shù)據(jù)也會受到其它因素的影響[16-17]。因此,若要獲得高精度、高通用性的冠層尺度干物質含量估算的模型,尚需將更多的影響因素在建模時引入估算模型中。

    [1]于龍鳳,李福恒,安福全,等.番茄不同節(jié)位葉片形態(tài)特征與干物質含量的關系[J].東北農業(yè)大學學報,2009,40(6):58-62.

    Yu L F,Li F H,An F Q,et al.Relation between leaf morphology features and dry matter content in different site of node of tomatoes[J].Journal of Northeast Acultural University,2009, 40(6):58-62.(in Chinese)

    [2]Romero A,Aguado I,Yebra M.Estimation of dry matter content in leaves using normalized indexes and PROSPECT model inversion[J].International Journal of Remote Sensing, 2012,33(2):396-414.

    [3]李玉霖,崔健垣,蘇永中.不同沙丘生境主要植物比葉面積和葉干物質含量的比較[J].生態(tài)學報,2005,25(2):304-311.

    Li Y L,Cui J H,Su Y Z.Specific leaf area and leaf dry matter content of some plants in different dune habitats[J].Act Ecologica Sinica,2005,25(2):304-311. (in Chinese)

    [4]潘瑞熾,王小青,李娘輝.植物生理[M].北京:高等教育出版社,2004:1-309.

    Pan R C,Wang X Q,Li N H.Plant physiology[M]. Beijing:Higher Education Press,2004:1-309.(in Chinese)

    [5]肖艷芳,周德民,趙文吉.輻射傳輸模型多尺度反演植被理化參數(shù)研究進展[J].生態(tài)學報,2013,33(11): 3291-3297.

    Xiao Y F,Zhou D M,Zhao W J.Review of inversing biophysical and biochemical vegetation parameters in various spatial scales using radiative transfer models[J]. Acta Ecologica Sinica,2013,33(11): 3291-3297. (in Chinese)

    [6]Feret J,Fran?ois C,Asner G P,et al.PROSPECT-4 and 5:advances in the leaf optical properties model separating photosynthetic pigments[J].Remote Sensing of Environment,2008,112(6): 3030- 3043.

    [7]張林,羅天祥,鄧坤枚,等.云南松比葉面積和葉干物質含量隨冠層高度的垂直變化規(guī)律[J].北京林業(yè)大學學報,2008, 30(1):40-44.

    Zhang L,Luo T X,Deng K M,et al.Vertical variations in specific leaf area and leaf dry matter content with canopy height in pinus yunnanensis[J].Journal of Beijing Forestry University,2008,30(1):40-44.(in Chinese)

    [8]Ria?o D,Vaughan P,Chuvieco E,et al.Estimation of fuel moisture content by inversion of radiative transfer models to simulate equivalent water thickness and dry matter content: analysis at leaf and canopy level[J].Geoscience and Remote Sensing, 2005,43(4):819-826.

    [9]Maire G,Fran?ois C,Soudani K,et al.Calibration and validation of hyperspectral indices for the estimation of broadleaved forest leaf chlorophyll content,leaf mass per area,leaf area index and leaf canopy biomass[J].Remote Sensing of Environment,2008, 112(10):3846-3864.

    [10]Broge N H,Mortensen J V.Deriving green crop area index and canopy chlorophyll density of winter wheat from spectral reflectance data[J].Remote Sensing of Environment, 2002,81(1):45-57.

    [11]Verhoef W.Light scattering by leaf layers with application to canopy reflectance modeling:the SAIL Model[J].Remote Sensing of Environment, 1984,16(2):125-141.

    [12]張佳華,郭文娟,姚鳳梅.植被水分遙感監(jiān)測模型的研究[J].應用基礎與工程科學學報,2007,15(1): 45-53.

    Zhang J H,Guo W J,Yao F M.The study on vegetation water content estimating model based on remote sensing technique [J].Journal of Basic Science and Engineering,2007,15(1): 45-53.(in Chinese)

    [13]趙永光,馬靈玲,李傳榮,等.基于冠層輻射傳輸模型的地表反射率光譜重建方法[J].光譜學與光譜分析,2015,35(7): 1763-1769.

    Zhao Y G,Ma L L,Li C R,et al.A method to reconstruct surface reflectance spectrum from multispectral image based on canopy radiation transfer[J]. Spectroscopy and Spectral Analysis, 2015,35(7): 1763-1769.(in Chinese)

    [14]馬文勇,王訓明.基于高光譜分析的草地葉綠素含量估算研究進展[J].地理科學進展,2016,35(1): 25-34.

    Ma W Y,Wang X M.Progress on grassland chlorophyll content estimation by hyperspectral analysis[J]. Progress in Geography, 2016,35(1):25-34. (in Chinese)

    [15]程曉娟,楊貴軍,徐新剛,等.新植被水分指數(shù)的冬小麥冠層水分遙感估算[J].光譜學與光譜分析,2014, 34(12): 3391-3396.

    Cheng X J,Yang G J,Xu X G,et al.Estimating canopy water content in wheat based on new vegetation water index[J]. Spectroscopy and Spectral Analysis,2014,34 (12):3391-3396.(in Chinese)

    [16]佃袁勇.高光譜數(shù)據(jù)反演植被信息的研究[D].武漢:武漢大學,2011.

    Dian Y Y.Inversion vegetation information base on hyper- spectral remote sensing data[D]. Wuhan: Wuhan University, 2011.(in Chinese)

    [17]Ali A M,Darvishzadeh R,Skidmore A K,et al. Estimating leaf functional traits by inversion of PROSPECT:assessing leaf dry matter content and specific leaf area in mixed mountainous forest[J]. International Journal of Applied Earth Observation, 2015,45(A):66-76.

    Estimating Dry Matter Content of Poplar Leaf by Hyperspectral Data

    CHENG Zhi-qing1,2,3, ZHANG Jin-song1,2,3, MENG Ping1, LI Yan-quan1, ZHENG Ning1

    (1.Research Institute of Forestry, Chinese Academy of Forestry, Beijing 100091,China; 2.Collaborative Innovation Center of Sustaintable Forestry in Southern China, Nanjing Forestry University, Nanjing 210037; 3. Key Laboratory of Tree Breeding and Cultivation, State Forestry Administration, Beijing 100091)

    In order to explore estimating the dry matter content of poplar leaves accurately, rapidly and accurately and nondestructive with hyperspectral data. Simulated values of poplar leave hyperspectral data which were leaf scales and canopy scales were obtained by PROSPECT and PROSAIL models which input parameters were measured by the physical and chemical data and soil background spectral data. At the same time, the dry matter content, equivalent water thickness and hyperspectral data were obtained from the measured leaf scales and canopy scales. Finally, the statistical methods were used to analyze the dry matter content in two scales of poplar leaves. The results showed that, the best band combination of the normalized difference vegetation index of dry matter was respectively band combination between 1685nm and 1704nm and band combination between 1551nm and 2143nm, respectively in leaf scales and canopy scales. Dry matter content models of leaf and canopy scales were built respectively by NDMI(1685,1704)and NDMI(1551,2143).The model validation accuracy of the leaf scale dry matter content was R2=0.663, RMSE=0.001g·cm-2, and the model validation accuracy of the canopy scale dry matter content was R2=0.91, RMSE=16.7g·m-2. Obviously, the hyperspectral technique has a high precision for estimating the dry matter content of poplar leaves, which can provide reference for the rapid and nondestructive estimation of dry matter content of poplar leaves.

    Technology of hyperspectrum;Poplar;Dry matter content;Estimation model

    10.3969/j.issn.1000-6362.2017.01.006

    2015?12?16

    國家科技支撐計劃課題(2015BAD07B05);林業(yè)公益性行業(yè)科研專項(201204105);國家自然科學基金項目(41105076);國家重大科學研究計劃(2012CB956202)

    程志慶(1986?),博士,講師,主要從事林業(yè)高光譜模型及氣象方面研究。E-mail:chengzhiqing1@126.com

    **通訊作者。E-mail:zhangjs@caf.ac.cn

    猜你喜歡
    冠層楊樹波段
    春日暖陽
    基于低空遙感的果樹冠層信息提取方法研究
    基于激光雷達的樹形靶標冠層葉面積探測模型研究
    安徽省淮南森林冠層輻射傳輸過程的特征
    大氣科學(2021年1期)2021-04-16 07:34:18
    跟蹤導練(五)
    口算本失蹤記
    楊樹山漫畫作品欣賞
    喜劇世界(2017年10期)2017-06-01 12:39:29
    施氮水平對冬小麥冠層氨揮發(fā)的影響
    M87的多波段輻射過程及其能譜擬合
    日常維護對L 波段雷達的重要性
    西藏科技(2015年4期)2015-09-26 12:12:58
    美女福利国产在线| 亚洲国产中文字幕在线视频| 天堂动漫精品| 久久久久久人人人人人| av电影中文网址| 大香蕉久久成人网| 欧美午夜高清在线| 一级黄色大片毛片| 亚洲欧洲精品一区二区精品久久久| av超薄肉色丝袜交足视频| 国产精品久久久人人做人人爽| 亚洲精品一卡2卡三卡4卡5卡| 黑人猛操日本美女一级片| 精品国产乱码久久久久久男人| 女人高潮潮喷娇喘18禁视频| 午夜福利影视在线免费观看| 国产精品亚洲av一区麻豆| 久久人妻熟女aⅴ| 纯流量卡能插随身wifi吗| 久久国产亚洲av麻豆专区| 两性夫妻黄色片| 日本黄色日本黄色录像| 日韩精品青青久久久久久| 一级黄色大片毛片| 少妇的丰满在线观看| a级毛片在线看网站| 色精品久久人妻99蜜桃| 在线观看午夜福利视频| 国产伦人伦偷精品视频| 少妇裸体淫交视频免费看高清 | 波多野结衣高清无吗| 日韩国内少妇激情av| www.熟女人妻精品国产| 久9热在线精品视频| 黄色视频,在线免费观看| 亚洲成人久久性| 亚洲五月色婷婷综合| 亚洲av片天天在线观看| e午夜精品久久久久久久| 国内毛片毛片毛片毛片毛片| 一个人免费在线观看的高清视频| 亚洲男人天堂网一区| 伊人久久大香线蕉亚洲五| 国产精品一区二区精品视频观看| 欧美乱色亚洲激情| 久久99一区二区三区| 欧美乱色亚洲激情| 熟女少妇亚洲综合色aaa.| 久久久国产成人免费| 91国产中文字幕| www.www免费av| 69av精品久久久久久| 91成年电影在线观看| 天堂俺去俺来也www色官网| 怎么达到女性高潮| 国产精品 欧美亚洲| 成人黄色视频免费在线看| 丝袜人妻中文字幕| 高清黄色对白视频在线免费看| 俄罗斯特黄特色一大片| 国产精品九九99| 999久久久精品免费观看国产| 搡老熟女国产l中国老女人| 丰满迷人的少妇在线观看| 青草久久国产| 久久天躁狠狠躁夜夜2o2o| 日韩av在线大香蕉| 久久亚洲精品不卡| 99国产极品粉嫩在线观看| 亚洲国产看品久久| 啦啦啦免费观看视频1| 亚洲成人国产一区在线观看| 国产精品二区激情视频| 88av欧美| 12—13女人毛片做爰片一| 香蕉久久夜色| √禁漫天堂资源中文www| av福利片在线| 欧美成人性av电影在线观看| 99热只有精品国产| 欧美+亚洲+日韩+国产| 欧美不卡视频在线免费观看 | 午夜影院日韩av| 午夜福利影视在线免费观看| 麻豆av在线久日| 两性夫妻黄色片| 人成视频在线观看免费观看| 99国产综合亚洲精品| 久久中文看片网| 免费看a级黄色片| 岛国在线观看网站| 免费在线观看日本一区| 操美女的视频在线观看| 色精品久久人妻99蜜桃| 五月开心婷婷网| 国产主播在线观看一区二区| 不卡一级毛片| 高清欧美精品videossex| 国产片内射在线| 国产麻豆69| 亚洲国产欧美一区二区综合| 日韩欧美三级三区| 99国产极品粉嫩在线观看| 少妇 在线观看| 好男人电影高清在线观看| 美女国产高潮福利片在线看| 人妻丰满熟妇av一区二区三区| 国产亚洲欧美98| 国产成人精品久久二区二区91| 成年人免费黄色播放视频| 美女高潮喷水抽搐中文字幕| 不卡一级毛片| 久久 成人 亚洲| 18美女黄网站色大片免费观看| 国产又色又爽无遮挡免费看| 精品福利观看| 亚洲欧美激情综合另类| 午夜福利,免费看| 亚洲国产欧美网| 国产黄色免费在线视频| 91字幕亚洲| 久久久久久久久免费视频了| 亚洲精品中文字幕一二三四区| 每晚都被弄得嗷嗷叫到高潮| 黑人操中国人逼视频| 精品卡一卡二卡四卡免费| 国产人伦9x9x在线观看| 久久久久久久久中文| 久久久久国产精品人妻aⅴ院| 亚洲精品成人av观看孕妇| 黄色 视频免费看| 日韩欧美国产一区二区入口| 欧美日韩黄片免| 亚洲熟女毛片儿| 激情视频va一区二区三区| 97碰自拍视频| 身体一侧抽搐| 色哟哟哟哟哟哟| 亚洲少妇的诱惑av| 嫩草影视91久久| 一进一出抽搐gif免费好疼 | 99久久精品国产亚洲精品| 一级毛片精品| 法律面前人人平等表现在哪些方面| 久久这里只有精品19| 后天国语完整版免费观看| 夜夜躁狠狠躁天天躁| 欧美在线一区亚洲| 黄色视频不卡| 99国产精品免费福利视频| 国产蜜桃级精品一区二区三区| 男女下面进入的视频免费午夜 | 成人国产一区最新在线观看| 自线自在国产av| 国产亚洲精品第一综合不卡| 国产亚洲精品第一综合不卡| 精品一区二区三区视频在线观看免费 | 久久久久久久午夜电影 | 久久精品国产99精品国产亚洲性色 | 老熟妇仑乱视频hdxx| 亚洲va日本ⅴa欧美va伊人久久| 手机成人av网站| 亚洲熟女毛片儿| 国产精品影院久久| av在线播放免费不卡| 十八禁网站免费在线| 亚洲成av片中文字幕在线观看| 国产一区在线观看成人免费| 久久久久国产一级毛片高清牌| 男人的好看免费观看在线视频 | 少妇被粗大的猛进出69影院| 国产精品久久久av美女十八| 涩涩av久久男人的天堂| 欧美一区二区精品小视频在线| 日韩中文字幕欧美一区二区| 亚洲欧美精品综合一区二区三区| 又黄又爽又免费观看的视频| 欧美性长视频在线观看| netflix在线观看网站| 麻豆一二三区av精品| 高潮久久久久久久久久久不卡| 高清在线国产一区| 电影成人av| 国产伦人伦偷精品视频| 国产国语露脸激情在线看| 色婷婷av一区二区三区视频| 超碰成人久久| 亚洲欧美激情在线| 亚洲av五月六月丁香网| 女人爽到高潮嗷嗷叫在线视频| 午夜日韩欧美国产| 18+在线观看网站| 小蜜桃在线观看免费完整版高清| 99久久精品一区二区三区| 午夜精品一区二区三区免费看| 男女那种视频在线观看| 中文在线观看免费www的网站| 久久久久精品国产欧美久久久| 日韩欧美在线二视频| 国产成人啪精品午夜网站| 蜜桃久久精品国产亚洲av| 亚洲av不卡在线观看| 日韩欧美在线二视频| 天堂av国产一区二区熟女人妻| 精品久久久久久久人妻蜜臀av| 日日摸夜夜添夜夜添av毛片 | 内地一区二区视频在线| 麻豆一二三区av精品| 国语自产精品视频在线第100页| 欧美一区二区精品小视频在线| 中文亚洲av片在线观看爽| 1000部很黄的大片| 老熟妇仑乱视频hdxx| 在线观看av片永久免费下载| 91久久精品电影网| 91久久精品电影网| 亚洲精品在线观看二区| 亚洲欧美日韩无卡精品| eeuss影院久久| 久久精品国产清高在天天线| 亚洲精品色激情综合| 久久久久久国产a免费观看| 精品久久久久久久久av| 亚洲av熟女| 国产成人av教育| 观看免费一级毛片| 在线十欧美十亚洲十日本专区| 欧美三级亚洲精品| 悠悠久久av| 少妇人妻一区二区三区视频| 亚洲熟妇熟女久久| 黄片小视频在线播放| 国产黄色小视频在线观看| 在线免费观看不下载黄p国产 | 18禁在线播放成人免费| 国产伦人伦偷精品视频| 欧美在线一区亚洲| 亚洲成人久久爱视频| 国产视频内射| 级片在线观看| 欧美日韩国产亚洲二区| 成人毛片a级毛片在线播放| 午夜a级毛片| 国产精品自产拍在线观看55亚洲| 一个人免费在线观看电影| 天天躁日日操中文字幕| 亚洲五月婷婷丁香| 可以在线观看的亚洲视频| 国产视频内射| 在线播放国产精品三级| 脱女人内裤的视频| 国产伦在线观看视频一区| 久久国产精品影院| 精品一区二区三区人妻视频| 国产真实伦视频高清在线观看 | 亚洲电影在线观看av| 亚洲欧美日韩无卡精品| 国产一区二区亚洲精品在线观看| www.www免费av| 亚洲人成伊人成综合网2020| 免费在线观看亚洲国产| 久久久久精品国产欧美久久久| 真人做人爱边吃奶动态| 内射极品少妇av片p| 国产精品综合久久久久久久免费| 亚洲片人在线观看| 国产国拍精品亚洲av在线观看| 岛国在线免费视频观看| 国产麻豆成人av免费视频| 亚洲av成人精品一区久久| 国产黄片美女视频| 亚洲欧美日韩无卡精品| 我要看日韩黄色一级片| 免费人成视频x8x8入口观看| 国产精华一区二区三区| 婷婷丁香在线五月| 在线播放国产精品三级| 国产综合懂色| 赤兔流量卡办理| 人妻制服诱惑在线中文字幕| 一区二区三区激情视频| 日韩免费av在线播放| 男女下面进入的视频免费午夜| 亚洲内射少妇av| 禁无遮挡网站| 久久伊人香网站| 99国产精品一区二区三区| 婷婷精品国产亚洲av在线| 国产成人福利小说| 美女cb高潮喷水在线观看| 国产美女午夜福利| 欧美国产日韩亚洲一区| 男插女下体视频免费在线播放| 在线天堂最新版资源| 中文字幕人成人乱码亚洲影| 超碰av人人做人人爽久久| 国产精品亚洲av一区麻豆| 欧美日本视频| 在线a可以看的网站| 97人妻精品一区二区三区麻豆| 亚洲欧美日韩卡通动漫| 欧美成人免费av一区二区三区| 白带黄色成豆腐渣| 亚洲美女视频黄频| 国产精品爽爽va在线观看网站| 国产亚洲精品久久久com| 蜜桃久久精品国产亚洲av| 亚洲经典国产精华液单 | 岛国在线免费视频观看| 色av中文字幕| 婷婷亚洲欧美| 国产真实伦视频高清在线观看 | 日本与韩国留学比较| 久久久久久久久大av| 欧美一区二区亚洲| 国产亚洲精品久久久com| 老司机午夜福利在线观看视频| 99热这里只有是精品50| 欧美激情国产日韩精品一区| 欧美三级亚洲精品| 91av网一区二区| 国产成人福利小说| 老鸭窝网址在线观看| 在线看三级毛片| 精品一区二区三区av网在线观看| 日韩精品中文字幕看吧| 亚洲 国产 在线| 高清日韩中文字幕在线| 欧美一区二区亚洲| 国产真实乱freesex| 蜜桃久久精品国产亚洲av| 男女视频在线观看网站免费| 成人无遮挡网站| 久久99热6这里只有精品| 老熟妇乱子伦视频在线观看| 黄片小视频在线播放| 非洲黑人性xxxx精品又粗又长| 亚洲天堂国产精品一区在线| 一级黄色大片毛片| 免费av毛片视频| 日本 欧美在线| 国产v大片淫在线免费观看| 国产精品一区二区三区四区免费观看 | 国产三级中文精品| 中文字幕熟女人妻在线| 一级黄片播放器| 午夜精品在线福利| 精品久久久久久久久久久久久| 亚洲一区二区三区色噜噜| 18禁在线播放成人免费| 久久久久精品国产欧美久久久| 12—13女人毛片做爰片一| 亚洲av第一区精品v没综合| av在线老鸭窝| 亚洲 国产 在线| 如何舔出高潮| 给我免费播放毛片高清在线观看| 亚洲av电影不卡..在线观看| 少妇熟女aⅴ在线视频| 久久精品国产亚洲av香蕉五月| 日本 欧美在线| 老司机午夜福利在线观看视频| 国产伦一二天堂av在线观看| 波多野结衣高清无吗| 天堂动漫精品| 男人的好看免费观看在线视频| www.熟女人妻精品国产| 欧美不卡视频在线免费观看| 97人妻精品一区二区三区麻豆| 18禁在线播放成人免费| 少妇人妻精品综合一区二区 | 精品久久久久久久久亚洲 | 成年女人永久免费观看视频| 熟女人妻精品中文字幕| 九色国产91popny在线| 老熟妇乱子伦视频在线观看| 国产91精品成人一区二区三区| 好男人在线观看高清免费视频| 午夜视频国产福利| 亚洲美女黄片视频| 性插视频无遮挡在线免费观看| 99久久成人亚洲精品观看| 欧美日韩黄片免| 高清日韩中文字幕在线| 亚洲精华国产精华精| 校园春色视频在线观看| 久久天躁狠狠躁夜夜2o2o| 男女做爰动态图高潮gif福利片| 日韩欧美免费精品| 在线看三级毛片| 午夜日韩欧美国产| 成人国产一区最新在线观看| 免费大片18禁| 免费看日本二区| 日韩精品中文字幕看吧| avwww免费| 91久久精品国产一区二区成人| 国产综合懂色| 精品人妻偷拍中文字幕| 听说在线观看完整版免费高清| 夜夜夜夜夜久久久久| 成人高潮视频无遮挡免费网站| 一进一出抽搐gif免费好疼| 亚洲av免费在线观看| 亚洲成av人片免费观看| 国产色爽女视频免费观看| 亚洲中文字幕日韩| 黄片小视频在线播放| 午夜福利视频1000在线观看| 黄色丝袜av网址大全| 女同久久另类99精品国产91| 色综合亚洲欧美另类图片| 18禁裸乳无遮挡免费网站照片| 国产精品一区二区免费欧美| 1000部很黄的大片| 久久久久久大精品| 国产精品亚洲av一区麻豆| 精品福利观看| 波多野结衣高清作品| 国内精品久久久久久久电影| 日本熟妇午夜| 免费av毛片视频| 国产精品嫩草影院av在线观看 | 国产精品av视频在线免费观看| 毛片一级片免费看久久久久 | 丰满人妻一区二区三区视频av| 禁无遮挡网站| 小蜜桃在线观看免费完整版高清| 村上凉子中文字幕在线| 成人午夜高清在线视频| 亚洲欧美精品综合久久99| 久久国产精品影院| 精品久久久久久久久av| 一进一出抽搐动态| 熟妇人妻久久中文字幕3abv| 草草在线视频免费看| 欧美成狂野欧美在线观看| 又黄又爽又刺激的免费视频.| 三级国产精品欧美在线观看| 亚洲精品粉嫩美女一区| 国产中年淑女户外野战色| 国产伦一二天堂av在线观看| av国产免费在线观看| 性插视频无遮挡在线免费观看| 日韩欧美国产在线观看| 欧美日韩黄片免| 成人特级av手机在线观看| 最新在线观看一区二区三区| 少妇熟女aⅴ在线视频| 男女做爰动态图高潮gif福利片| 俺也久久电影网| 国产精品野战在线观看| 2021天堂中文幕一二区在线观| 男女之事视频高清在线观看| 亚洲在线自拍视频| 久久精品国产自在天天线| 午夜福利视频1000在线观看| 99热精品在线国产| 亚洲熟妇中文字幕五十中出| 99热6这里只有精品| 天天一区二区日本电影三级| 亚洲av.av天堂| 三级毛片av免费| 内射极品少妇av片p| 在线免费观看的www视频| 日韩 亚洲 欧美在线| 精品福利观看| 欧美一区二区国产精品久久精品| 熟女人妻精品中文字幕| 亚洲一区二区三区不卡视频| 精品国产三级普通话版| 亚洲成av人片免费观看| 午夜日韩欧美国产| www日本黄色视频网| 午夜老司机福利剧场| 99久久精品一区二区三区| 美女大奶头视频| 久久久久九九精品影院| 中文字幕高清在线视频| www日本黄色视频网| 男插女下体视频免费在线播放| a级一级毛片免费在线观看| 成人三级黄色视频| 高清毛片免费观看视频网站| 国产伦精品一区二区三区四那| 亚洲国产精品999在线| 国产精品久久视频播放| 怎么达到女性高潮| 午夜影院日韩av| 最后的刺客免费高清国语| 国产精品野战在线观看| 亚洲精品成人久久久久久| 国产亚洲av嫩草精品影院| 小蜜桃在线观看免费完整版高清| eeuss影院久久| 久久久色成人| h日本视频在线播放| 亚洲性夜色夜夜综合| 亚洲不卡免费看| 嫩草影视91久久| 久久6这里有精品| 在线观看舔阴道视频| 丰满人妻一区二区三区视频av| 欧美3d第一页| 午夜福利高清视频| 午夜老司机福利剧场| 少妇熟女aⅴ在线视频| a级毛片a级免费在线| 九色国产91popny在线| 熟妇人妻久久中文字幕3abv| 91麻豆精品激情在线观看国产| 男人狂女人下面高潮的视频| av在线蜜桃| 免费搜索国产男女视频| 亚洲午夜理论影院| 国产精品一及| 性欧美人与动物交配| 成年女人看的毛片在线观看| 午夜亚洲福利在线播放| 精品一区二区三区av网在线观看| 国产真实伦视频高清在线观看 | 欧美日本亚洲视频在线播放| 高清毛片免费观看视频网站| 91狼人影院| 亚州av有码| 如何舔出高潮| 午夜福利高清视频| 亚洲精品影视一区二区三区av| 国产激情偷乱视频一区二区| 在线观看一区二区三区| 能在线免费观看的黄片| 婷婷精品国产亚洲av| 亚洲狠狠婷婷综合久久图片| 国产 一区 欧美 日韩| 国产午夜精品久久久久久一区二区三区 | 亚洲精品亚洲一区二区| 亚洲精品久久国产高清桃花| 成人永久免费在线观看视频| 亚洲熟妇中文字幕五十中出| 亚洲人成网站在线播| 波多野结衣高清无吗| 男插女下体视频免费在线播放| 亚洲va日本ⅴa欧美va伊人久久| 免费高清视频大片| 成熟少妇高潮喷水视频| 欧美极品一区二区三区四区| 国产高清激情床上av| 午夜精品久久久久久毛片777| h日本视频在线播放| 极品教师在线免费播放| 十八禁国产超污无遮挡网站| 午夜福利欧美成人| 精品无人区乱码1区二区| 夜夜看夜夜爽夜夜摸| 性色av乱码一区二区三区2| 99精品在免费线老司机午夜| www日本黄色视频网| 免费在线观看成人毛片| 色哟哟哟哟哟哟| 在线观看午夜福利视频| 亚洲国产精品合色在线| 乱人视频在线观看| 欧美一级a爱片免费观看看| 久久99热这里只有精品18| 女人被狂操c到高潮| 91午夜精品亚洲一区二区三区 | 亚洲美女视频黄频| 亚洲精品456在线播放app | 精品不卡国产一区二区三区| 狠狠狠狠99中文字幕| 9191精品国产免费久久| 亚洲熟妇熟女久久| 午夜久久久久精精品| 久久精品夜夜夜夜夜久久蜜豆| 人妻久久中文字幕网| 中文字幕人成人乱码亚洲影| 色视频www国产| 村上凉子中文字幕在线| 欧美黄色淫秽网站| 午夜免费成人在线视频| 国产成人aa在线观看| 12—13女人毛片做爰片一| 十八禁人妻一区二区| 亚洲av二区三区四区| 狠狠狠狠99中文字幕| 一本综合久久免费| 久久久久久久午夜电影| 国产欧美日韩精品一区二区| 天堂动漫精品| 国内精品久久久久精免费| 美女黄网站色视频| 黄色丝袜av网址大全| 日本 av在线| 最新在线观看一区二区三区| 日本五十路高清| 国产视频内射| 自拍偷自拍亚洲精品老妇| 日本撒尿小便嘘嘘汇集6| 国产三级中文精品| 亚洲人与动物交配视频| 亚洲专区中文字幕在线| 成人av在线播放网站| 国产高潮美女av| 亚洲成av人片免费观看| 欧美日韩黄片免| av女优亚洲男人天堂| 国产午夜福利久久久久久| 亚洲成av人片免费观看| 999久久久精品免费观看国产| 亚洲成av人片免费观看| 欧美日韩黄片免| 麻豆成人午夜福利视频| 国内少妇人妻偷人精品xxx网站| 亚洲精品影视一区二区三区av| 一个人看视频在线观看www免费|