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

    不同病蟲害危害程度下云南松針葉光譜模擬及其生化參數(shù)敏感性分析

    2016-06-15 16:40:58林起楠黃華國俞琳鋒
    光譜學與光譜分析 2016年8期
    關(guān)鍵詞:云南松針葉白化

    林起楠, 黃華國, 陳 玲, 俞琳鋒, 黃 侃

    北京林業(yè)大學省部共建森林培育與保護教育部重點實驗室, 北京 100083

    不同病蟲害危害程度下云南松針葉光譜模擬及其生化參數(shù)敏感性分析

    林起楠, 黃華國*, 陳 玲, 俞琳鋒, 黃 侃

    北京林業(yè)大學省部共建森林培育與保護教育部重點實驗室, 北京 100083

    明確植被葉片反射率對生化參數(shù)敏感程度是遙感定量反演生化參數(shù)含量的前提。 本文以不同健康狀態(tài)下的云南松針葉光譜為例, 基于查找表法(LUT)對針葉光譜模型LIBERTY的葉綠素和水分吸收系數(shù)進行了標定。 提出了一種新的光譜合成方法和枯黃指數(shù)YI, 解決LIBERTY模型不能有效擬合嚴重受害程度的云南松針葉反射光譜曲線的問題。 并應(yīng)用EFAST全局敏感性分析方法, 利用標定后的LIBERTY模型定量分析了不同受害程度的針葉生化參數(shù)敏感性變化。 結(jié)果表明: (1) 標定后的LIBERTY模型能有效擬合健康針葉光譜(R2=0.999, RMSE<0.01)、 輕度(R2=0.991, RMSE<0.02)和中度受害針葉光譜(R2=0.992, RMSE<0.03); 但是嚴重受害針葉光譜模擬能力很差(R2=0.803, RMSE>0.1)。 (2) 新光譜合成方法能有效擬合嚴重受害光譜(R2=0.991, RMSE<0.03); 枯黃指數(shù)(YI)能夠定量分析受害程度。 (3) 隨著受害程度加深, 針葉反射率對葉綠素、 水分參數(shù)敏感度降低, 而基吸收、 白化吸收等參數(shù)敏感度增大, 且存在敏感波段。 如在505~565和705~850 nm范圍白化吸收對針葉反射率作用明顯; 對于嚴重受害針葉, 除了葉綠素之外白化吸收對可見光區(qū)針葉反射率影響也非常顯著; 通過引入白化吸收先驗知識, 發(fā)現(xiàn)如果能進一步確定白化吸收情況, 能夠提高嚴重受害針葉葉綠素反演的精度。

    LIBERTY; 針葉反射率; 敏感度分析; 云南松; EFAST

    引 言

    云南松(Pinusyunnanensis)在我國主要分布于西南地區(qū)。 近年來, 在云南地區(qū)發(fā)生了嚴重干旱, 進而引起松毛蟲和切梢小蠹等害蟲大量爆發(fā), 導致云南松人工林遭受大面積的破壞。 通過遙感的手段對云南松病蟲害進行監(jiān)測預(yù)警, 對預(yù)防保護和災(zāi)后評估有重要意義。 光譜分析是森林病蟲害遙感監(jiān)測的常用手段[1]。 其主要理論假設(shè)是受害森林冠層的結(jié)構(gòu)、 生理和生化參數(shù)發(fā)生變化, 該變化與光譜存在良好的相關(guān)性。 盡管云南地區(qū)的病蟲害遙感應(yīng)用研究不少, 但是鮮有從機理上闡述病蟲害危害與針葉光譜的定量關(guān)系, 并進行病害程度及其生化參數(shù)敏感性分析的研究。

    針葉光譜模型LIBERTY(leaf incorporating biochemistry exhibiting reflectance and transmittance yields)[2], 能夠從機理上模擬各類理化參數(shù)在不同含量下的針葉光譜反射和透射特征, 在生化參數(shù)反演、 光譜特征分析等方面研究運用廣泛[3]。 Barton等利用LIBERTY模型分析了不同黃化病斑塊比率對葉片的反射率和冠層光譜特征的影響[4]; 王雷等基于標定的LIBERTY模型分析不同水分條件下針葉光譜特征, 為落葉松病蟲害監(jiān)測提供了理論基礎(chǔ)[5]; Coops等通過測定不同受害程度的加拿大輻射松(Pinusradiata)葉綠素含量結(jié)合LIBERTY模型構(gòu)建光譜指數(shù)實現(xiàn)對葉綠素含量估算[6]。 尚未發(fā)現(xiàn)利用LIBERTY模型對云南松的研究案例, 及枯黃針葉光譜模擬相關(guān)研究。

    明確不同受害程度下針葉反射率對理化參數(shù)的敏感性, 是精確反演植被針葉理化參數(shù)的前提。 在LIBERTY模型理化參數(shù)敏感性分析上, Dawson等以LIBERTY模型分析了干、 濕針葉的木質(zhì)素和纖維素含量變化, 對針葉光譜的影響[2]; 顏春燕將LIBERTY模型參數(shù)歸一化后以不同參考點定量計算理化參數(shù)對針葉反射率的平均敏感度。 Li和Wang采用sobol法定量分析了LIBERTY模型的理化參數(shù)對針葉可見光區(qū)反射率的貢獻度[7]。 但在對LIBERTY模型敏感性分析研究中, 都忽略了針葉結(jié)構(gòu)參數(shù)與生化參數(shù)對針葉反射率作用的相互關(guān)系。

    為此, 本文試圖解決以上三個問題。 首先, 基于實測的針葉光譜數(shù)據(jù), 解決LIBERTY模型在云南松樹種上的標定問題。 然后, 利用標定后的模型對健康、 輕度、 中度、 嚴重受害的云南松針葉反射光譜曲線進行模擬, 重點解決嚴重受害針葉光譜模擬問題。 最后, 結(jié)合EFAST(extended Fourier amptitude sensitivity test)方法分析不同受害程度下針葉生化參數(shù)的敏感性分析。

    1 實驗部分

    1.1 光譜數(shù)據(jù)

    2014年7月中旬至9月中旬, 在云南省玉溪市紅塔區(qū)、 江川縣和墨江縣聯(lián)珠鎮(zhèn)、 魚塘鎮(zhèn)地區(qū)云南松人工林進行健康和受害云南松的針葉采樣, 并同步測量光譜。

    實驗針對受云南松毛蟲和小蠹等侵害的云南松進行研究, 按針葉受害的程度分為健康、 輕度受害、 中度受害和重度受害4種類型。 從樹冠的四個垂直方向采集針葉, 將采集到的針葉用塑封袋保存好, 立即帶回實驗室進行葉綠素測量和光譜測量。 采用AvaSpec-EDU-VIS(350~850 nm, 分辨率0.29 mn)和AvaSpec-NIR1.7(900~1 750 nm, 分辨率2.0~5.0 nm)兩個光譜儀互補聯(lián)合觀測, 可獲得350~1 750 nm的光譜數(shù)據(jù)。 實驗觀測時將同一類型的針葉多層緊密平鋪于觀測的卡槽容器中, 每層平均40個針葉共5層, 將光纖探頭放置在觀測支架中進行觀測, 每個針葉樣本重復(fù)5次取平均值作為樣本最終光譜值。 葉綠素含量測定采用分光光度計方法測定, 將針葉剪碎后用石英砂研磨均勻, 用無水乙醇于黑暗室溫內(nèi)浸提, 測定649和665 nm處的吸光度, 以無水乙醇作為空白對照組。 測定的不同類型云南松針葉葉綠素含量范圍為: 健康針葉約為288.35~337.30 mg·m-2; 輕度受害針葉約為144.03~230.06 mg·m-2; 中度受害針葉約為115.44~128.64 mg·m-2; 嚴重受害針葉約為29.80~57.25 mg·m-2。

    對實測云南松針葉反射光譜, 每種受害類型分別取10個樣本。 為了與LIBERTY模型模擬光譜分辨率一致將實測光譜分辨率重采樣為5 nm, 分析的光譜波長范圍為400~1 700 nm。

    1.2 LIBERTY針葉反射模型及其標定方法

    LIBERTY模型是基于針葉特征發(fā)展而來的輻射傳輸模型[2]。 該模型假設(shè)針葉細胞是球形顆粒的漫射體并且遵從朗伯余弦定理, 由若干層排列緊密的細胞的反射率和透過率迭代過程求得單葉的反射率和透過率。 模型能夠模擬針葉單葉或針葉簇可見光到近紅外波段范圍(400~2 500 nm, 間隔5 nm)內(nèi)的光譜特性。 模型有9個輸入?yún)?shù)(詳見表1), 包括3個結(jié)構(gòu)參數(shù)(細胞直徑、 細胞間隙、 針葉厚度)和6個生化參數(shù)(基吸收、 白化吸收、 葉綠素、 水分、 木質(zhì)素和纖維素、 氮)。 除此之外還有5個缺省的生化吸收光譜系數(shù): 葉綠素吸收系數(shù)fchl、 水分吸收系數(shù)fw、 木質(zhì)素吸收系數(shù)flc、 白化吸收系數(shù)fa、 氮元素吸收系數(shù)fn。

    表1 LIBERTY模型9個輸入?yún)?shù)、 單位、 取值范圍和查找表設(shè)置

    基于實測的反射光譜數(shù)據(jù)利用LIBERTY模型進行逆向反演, 對模型輸入?yún)?shù)進行標定, 來彌補部分參數(shù)未測量或難于測量的缺陷。 通過建立查找表(LUT)(相關(guān)設(shè)定見表1)的方法來進行光譜反演, 為減小計算量舍去不合理組合共模擬10 000條針葉反射率光譜。 為減少LIBERTY模型標定過程中輸入?yún)?shù)個數(shù), 對不同參數(shù)分別在不同波段范圍內(nèi)進行標定[8]。 例如D和xu兩個細胞結(jié)構(gòu)參數(shù)其主要作用于700~1 000 nm波段范圍之間[9]; 葉綠素的光譜吸收范圍在400~800 nm的可見光區(qū); 而水分的光譜吸收范圍在大于1 000 nm的近紅外波段[3]。 以實測光譜值和模型模擬光譜值的RMSE[式(1)]值優(yōu)化各參數(shù)值, 當RMSE值最小時即為參數(shù)的最優(yōu)標定值。

    (1)

    式中:Rmod為模型反射率模擬值,Rmes為反射率實測值

    LIBERTY模型默認的吸收系數(shù)是基于北美短葉松(Pinusbanksiana)獲得的[2], 對新的樹種云南松而言并不一定完全適用, 因此有必要對缺少的吸收系數(shù)進行重新標定[7, 10]。 葉綠素和水分是影響植物葉片可見光到近紅外光區(qū)光譜反射率的主要因素, 因此葉綠素吸收系數(shù)和水分吸收系數(shù)標定是重點考慮部分。 利用實測的反射光譜數(shù)據(jù)和LIBINV反演程序計算總吸收系數(shù)(fglobal)值, 再根據(jù)總的吸收系數(shù)和5個獨立的生化吸收系數(shù)之間的關(guān)系[式(2)], 將計算得到參數(shù)的標定值帶入式(2)中, 來實現(xiàn)對葉綠素和水分吸收系數(shù)的標定。 為確保標定系數(shù)的獨立性, 以實測光譜70%數(shù)據(jù)進行標定, 剩余30%數(shù)據(jù)用來驗證標定結(jié)果。

    (2)

    式中,cBaseline,calbino為基吸收和白化吸收;cchl,cw,clc,cn為葉綠素, 水分, 木質(zhì)素和纖維素中氮含量;fchl,fw,fa,flc,fn為5個吸收系數(shù)。

    最后, 用重新標定后的葉綠素、 水分系數(shù)的LIBERTY模型開展云南松針葉反射光譜曲線模擬。 通過對不同受害程度云南松針葉反射光譜模擬的大量實驗, 發(fā)現(xiàn)LIBERTY模型只能夠模擬出有“峰”、 “谷”特征的光譜曲線, 對嚴重受害針葉反射光譜模擬的誤差較大, 需要考慮新的模型改進嚴重受害類型的光譜模擬結(jié)果。

    為此, 提出一種新的加權(quán)模擬方法。 從ENVI光譜庫中篩選與實測嚴重受害反射光譜最接近的一條光譜, 經(jīng)計算Chris Elvidge的干植被波譜庫中灌木蒿(Sagebrush)葉片最接近(R2=0.913, RMSE=0.059), 并將該高光譜數(shù)據(jù)作為先驗知識, 通過線性加權(quán)擬合方法來改進LIBERTY模擬反射率值, 改進關(guān)系表達式如式(3)所示

    (3)

    式中:fsageb和fmod為加權(quán)系數(shù);Rsageb,Rmod和RM分別為灌木蒿、 LIBERTY和最終光譜模擬反射率值。fsageb能反映病蟲害針葉的受害嚴重程度, 值越接近于1, 越是枯黃。 因此, 將其命名為枯黃指數(shù)(yellow index, YI)。

    需要說明的是, 灌木蒿的光譜只是反映嚴重受害光譜的先驗知識, 和灌木蒿這個物種無本質(zhì)聯(lián)系, 也可以用其他更合適的光譜代替。

    1.3 EFAST方法

    EFAST方法是Saltelli等[11]將sobol方法與傅里葉振幅靈敏度檢驗法(FAST)相結(jié)合的全局敏感性分析方法, 其特點是穩(wěn)定、 樣本量少、 計算效率高。 該方法從方差分析角度出發(fā), 認為模型輸出結(jié)果方差是由輸入?yún)?shù)和參數(shù)之間相互作用引起, 以參數(shù)對模型輸出總方差的貢獻率來反映參數(shù)的敏感性。 輸入?yún)?shù)的敏感性可由主要敏感度和總敏感度表示, 其中主要敏感度為該參數(shù)單獨作用的敏感性, 總敏感度還包括參數(shù)間相互作用的敏感性。

    EFAST法主要是應(yīng)用一個搜索函數(shù)在參數(shù)的多維空間內(nèi)搜索, 對模型的輸入?yún)?shù)引入一個共同獨立參數(shù)函數(shù)并分別對各參數(shù)定義一個整數(shù)頻率, 使模型成為獨立參數(shù)的周期函數(shù)[12]。 然后, 對模型的輸出結(jié)果進行傅里葉分析, 以傅里葉振幅大小來表示模型對參數(shù)的敏感度, 振幅越大敏感度越高。

    假定一個模型Y=f(X),X=(x1,x2, …,xn), 模型總方差與輸入?yún)?shù)方差關(guān)系如式(4)

    (4)

    其中,V(Y)為模型總方差,Vi表示參數(shù)xi的方差,Vij, …,V(1, 2, …, n)表示各參數(shù)間相互作用的方差。

    參數(shù)xi的主要敏感度(一階敏感性指數(shù)), 由參數(shù)xi的方差對總方差V(Y)的貢獻率表示為式(5)

    (5)

    同理,xi參數(shù)與其他參數(shù)相互作用二階、 三階和高階敏感性指數(shù)表達如式(6)

    (6)

    由式(5)和式(6)可知, 模型輸入?yún)?shù)xi總的敏感度為

    (7)

    運用simlab2.2工具軟件中的EFAST敏感性分析模塊進行LIBERTY模型的敏感性分析。 LIBERTY模型輸入?yún)?shù)按均勻分布采樣, 利用simlab2.2軟件以EFAST采樣方法隨機生成5 000組輸入樣本。 將輸入樣本導入LIBERTY模型程序輸出5 000條針葉光譜, 將5 000個樣本輸出結(jié)果輸入simlab2.2軟件執(zhí)行Monte Carlo模擬計算各參數(shù)的敏感性。

    2 結(jié)果與討論

    2.1 云南松實測光譜曲線

    圖1為實測的不同蟲害受害程度的云南松針葉光譜反射率平均值。 可以看出, 在綠光區(qū)(500~560 nm), 不同受害程度的云南松光譜反射率大小出現(xiàn)遞減的規(guī)律, 即健康>輕度>中度>重度; 而在紅光區(qū)(620~760 nm), 不同受害程度的針葉反射率大小沒有明顯的規(guī)律性, 呈現(xiàn)出重度>健康>輕度>中度; 在760~900, 1 100~1 300, 1 400~1 600 nm波段范圍內(nèi), 隨著針葉受害加重, 針葉的反射率大小呈現(xiàn)出上升的趨勢, 即重度>中度>輕度>健康。 引起針葉光譜反射率發(fā)生變化的原因, 主要因不同受害程度云南松針葉其葉綠素和水分的含量不同導致。

    圖1 不同受害程度云南松針葉光譜曲線

    2.2 LIBERTY模型標定結(jié)果

    對葉綠素和水分吸收系數(shù)重新標定的結(jié)果, 經(jīng)驗證精度為95.4%, 標定后葉綠素、 水分吸收系數(shù)曲線變化明顯(圖2)。 基于LIBERTY模型默認的和標定后吸收系數(shù)分別擬合實測的不同受害程度的云南松針葉光譜, 發(fā)現(xiàn)默認的吸收系數(shù)的擬合誤差較大, 標定后的吸收系數(shù)擬合誤差減小明顯(圖3)。 其中, 健康類型針葉反射光譜曲線模擬精度最好(R2=0.999, RMSE<0.01), 輕度受害類型模擬結(jié)果(R2=0.991, RMSE<0.02)和中度受害類型模擬結(jié)果(R2=0.992, RMSE<0.03)次之, 嚴重受害類型譜模擬較差(R2=0.803, RMSE>0.1)。 經(jīng)過新的加權(quán)擬合方法, 嚴重受害針葉光譜的模擬精度顯著提升(R2=0.991, RMSE<0.03), 枯黃指數(shù)YI均值為0.664。 利用新的光譜擬合方法對非嚴重受害的針葉光譜進行擬合, 精度與不加權(quán)時未發(fā)生顯著變化, 效果仍然很好(圖3), YI指數(shù)均值與嚴重程度呈正向相關(guān): 健康0.018<輕度0.041<中度0.182<重度0.664。 由于本研究基于同一標定參數(shù), 在四個不同受害類型上, 都獲得了良好的擬合效果, 因此認為標定結(jié)果是合理的。

    圖2 LIBERTY吸收系數(shù)的標定值和缺省值

    圖3 LIBERTY標定前后不同受害程度云南松針葉實測的光譜曲線和模擬的光譜曲線

    2.3 LIBERTY模型針葉理化參數(shù)敏感性分析

    圖4顯示了在400~1 700 nm區(qū)間針葉反射率對模型各參數(shù)的一階敏感指數(shù)和參數(shù)間相互作用的結(jié)果。 在可見光區(qū), 葉綠素是影響針葉反射率的主要因素, 葉綠素對針葉反射率變化貢獻度約為78%; 在540 nm附近葉綠素一階敏感度值形成一個“低谷”, 之后敏感度迅速回升, 在660 nm達到峰值。 細胞直徑、 細胞間隙、 針葉厚度3個結(jié)構(gòu)參數(shù)在400~1 700 nm整個波段范圍內(nèi)對針葉反射率都有影響, 其中細胞直徑和細胞間隙參數(shù)貢獻度較大, 對700~1 000 nm對針葉反射率變化起主導作用, 該結(jié)果也驗證了Alasdair等[10]的結(jié)論。 水分敏感波段范圍是在近紅外和短波紅外光區(qū), 在1 400 nm附近敏感度達到最高值80%; 白化吸收敏感波段位于400~900 nm, 基吸收作用波段為700~1 400 nm, 木質(zhì)素纖維素和氮元素的影響波段大于1 400 nm。 總體而言, 基吸收、 白化吸收、 木質(zhì)素和纖維素、 氮對400~1 700 nm的針葉反射率影響較低<20%, 葉綠素和水分敏感波段明顯而且對針葉反射率的影響具有絕對主導性。

    為進一步分析模型結(jié)構(gòu)參數(shù)與生化參數(shù)對針葉反射率影響之間的相互作用。 首先探討在不同葉綠素和水分含量水平下, 理化參數(shù)之間對針葉反射率作用是否存在制約關(guān)系。 葉綠素和水分含量大小范圍設(shè)置為: 低(葉綠素: 0~50, 水分: 0~50), 中(葉綠素: 50~300, 水分: 50~200)、 高(葉綠素: 300~600, 水分: 200~500)三個水平, 計算各水平下參數(shù)的敏感度。

    圖4 理化參數(shù)對針葉反射率的敏感性

    由圖5(a)—(c)可知, 隨著葉綠素和水分含量升高, 細胞直徑和細胞間隙對針葉反射率的主導作用越來越明顯, 結(jié)構(gòu)參數(shù)對針葉反射率的貢獻度抑制了其他生化參數(shù)的作用, 且葉綠素、 水分含量越高抑制作用越明顯。 結(jié)果證明了Moorthy等[13]的試驗結(jié)論, 即模型在沒有細胞直徑和細胞間隙參數(shù)的先驗知識前提下, 無法取得精確的生化參數(shù)反演結(jié)果。 因此, 在不同受害程度針葉生化參數(shù)敏感性分析中, 以云南松針葉細胞直徑和細胞間隙的標定結(jié)果作為先驗知識來進行相關(guān)計算分析。

    2.4 不同受害程度針葉生化參數(shù)敏感性

    為了分析云南松在不同受害程度條件下對針葉生化參數(shù)敏感性的影響, 以葉綠素、 水分含量大小范圍來表征云南松針葉各受害類型[6]。 新的光譜擬合方法是基于LIBERTY模型模擬既定的結(jié)果加權(quán)得來, 而新增加權(quán)部分是固定光譜曲線, 敏感性不變, 因此敏感性分析為LIBERTY模型敏感度乘以(1-YI)。 不同受害類型的葉綠素和水分含量取值范圍參照葉綠素測量值及其與水分的標定結(jié)果和相關(guān)文獻[4, 6]設(shè)定。 設(shè)定情況如下: 嚴重受害(葉綠素: 0~50、 水分: 0~50、 YI=0.664), 中度受害(葉綠素: 50~150、 水分: 50~100、 YI=0.182), 輕度受害(葉綠素: 150~300 、 水分: 100~200、 YI=0.041), 健康(葉綠素: 300~600、 水分: 200~500、 YI=0.018)。 對于同齡的云南松針葉, 相同受害程度的針葉的結(jié)構(gòu)參數(shù)大小相對一致, 大小以標定值設(shè)定。 模擬對應(yīng)狀態(tài)下的針葉反射率并計算反射率對參數(shù)的敏感度。

    圖5 不同葉綠素、 水分含量大小的理化參數(shù)敏感性

    Fig.5 Sensitivity of leaf reflectance to biochemical parameters in three level of chlorophyll and water contents for LIBERTY

    圖6(a)—(d)分別是云南松四種不同受害程度針葉反射率對6個生化參數(shù)的一階敏感性及其參數(shù)間相互作用的影響。 隨著受害程度增強, 葉綠素和水分對針葉反射率的主導性減弱, 其他四個生化參數(shù)對針葉反射率的影響逐漸增大。

    在400~1 700 nm, 氮對不同受害類型的針葉反射率的貢獻較低<8%且遠小于其他生化參數(shù)的作用。 木質(zhì)素和纖維素與基吸收隨著受害程度加深對針葉反射率的影響增大, 但是在近紅外和短波紅外, 水分的影響還是占據(jù)主導, 對針葉水分反演的干擾作用小。 受害的針葉[圖6(b)—(d)], 在535和730 nm左右, 白化吸收有兩個很窄的高敏感區(qū); 尤其對嚴重受害針葉, 各參數(shù)敏感波段明顯, 白化吸收與葉綠素在可見光區(qū)對針葉反射率的貢獻相當。 因此, 在沒有白化吸收參數(shù)的先驗知識情況下, 對反演嚴重受害的針葉葉綠素含量將造成很大的誤差。

    圖6 不同受害程度下針葉生化參數(shù)對光譜的敏感性

    為此, 借鑒肖艷芳[14]的方法引入白化吸收參數(shù)先驗知識。 將白化吸收設(shè)置為: 0~1, 1~3, 3~4分別代表白化吸收低、 中、 高三種水平。 葉綠素對針葉反射率的作用主要在可見光區(qū), 因此計算基吸收、 白化吸收、 葉綠素在可見光區(qū)對針葉反射率的敏感度, 而其他參數(shù)在該波段范圍的敏感度很低, 因此不計算其敏感度。 如圖7(a)—(c)所示, 在一定范圍內(nèi), 白化吸收對針葉反射率的影響隨著白化吸收水平的增大而增大; 但增大到一定程度后白化吸收對針葉反射率的作用趨于飽和, 之后針葉反射率對白化吸收變化的敏感度降低。 在三種白化吸收水平下, 針葉反射率對嚴重受害針葉的葉綠素敏感度顯著提高, 主導作用明顯。 這表明, 在已知嚴重受害前提下, 如果能進一步確定其白化吸收水平, 就能夠得到比較準確的葉綠素反演結(jié)果。

    圖7 嚴重受害針葉在不同白化吸收水平下參數(shù)的敏感性

    Fig.7 Sensitivity of leaf reflectance to chemical parameters in three levels of albino absorption for seriously stressed needles

    3 結(jié) 論

    利用LIBERTY模型對不同受害程度的云南松針葉光譜進行模擬, 結(jié)果表明: 健康的針葉光譜模擬的精度最高, 輕度、 中度受害類型次之, 嚴重受害類型精度較低。 通過加權(quán)擬合的方法來改進LIBERTY模擬光譜值, 得到較好的光譜模擬結(jié)果, 并獲得能反映針葉受害程度的枯黃指數(shù)(YI)。

    討論了模型的結(jié)構(gòu)參數(shù)和生化參數(shù)對針葉反射率影響之間的相互作用, 發(fā)現(xiàn)兩者之間存在明顯的制約關(guān)系。 在不同葉綠素和水分含量水平下, 細胞直徑和細胞間隙兩個結(jié)構(gòu)參數(shù)對針葉反射率的貢獻度極易掩蓋其他參數(shù)的貢獻度。 利用LIBERTY模型進行生化參數(shù)反演時, 需要針葉細胞直徑和細胞間隙兩個結(jié)構(gòu)參數(shù)的先驗知識, 才能夠獲得較為精確的反演值。

    基于LIBERTY模型利用EFAST方法定量分析了生化參數(shù)對不同受害程度針葉反射率的敏感性。 結(jié)果表明: (1)總體來說, 葉綠素和水分是影響針葉反射率的主要因素。 (2)在不同的受害程度下, 各種生化參數(shù)對針葉反射率的敏感性具有明顯的差異。 對健康的針葉, 葉綠素和水分的貢獻度極易掩蓋其他生化參數(shù)對針葉反射率的作用, 而隨著受害程度加深其主導作用降低; 對于嚴重受害的針葉, 各參數(shù)都有明顯的主導敏感性波段, 尤其是白化吸收對可見光反射率影響與葉綠素參數(shù)相當。 (3)對于嚴重受害針葉, 如果明確白化吸收水平則能夠提高模型葉綠素反演的精度。 同時如何利用不同受害程度的針葉生化參數(shù)敏感性結(jié)果, 針對不同受害程度的針葉生化參數(shù)反演設(shè)定合理的分步反演策略也需要進一步探討。

    本研究是在針葉尺度上, 利用LIBERTY模型定量分析不同受害程度針葉生化參數(shù)的敏感度和敏感波段變化, 為病蟲害遙感監(jiān)測提供理論基礎(chǔ)。 下一步將耦合冠層模型研究在冠層尺度上分析病蟲害問題。

    [1] ZHANG Heng, PAN Jie, LIAO Zhen-feng(張 衡, 潘 潔, 廖振峰). Forest Pest and Disease(中國森林病蟲), 2014, 33(2): 26.

    [2] Dawson T P, Curran P J, Plummer S E. Remote Sensing of Environment, 1998, 65(1): 50.

    [3] Zarco-Tejada P J, Miller J R, Harron J, et al. Remote Sensing of Environment, 2004, 89(2): 189.

    [4] Barton C V. Tree Physiol., 2001, 21(12-13): 789.

    [5] Wang Lei, Luo Youqing, Huang Huaguo, et al. Forestry Studies in China, 2008, 1(11): 28.

    [6] Coops N C, Stone C. Australian Journal of Botany, 2005, 53(5): 417.

    [7] Li P, Wang Q. International Journal of Remote Sensing, 2013, 34(7): 2402.

    [8] Moorthy I, Miller J R, Noland T L. Remote Sensing of Environment, 2008, 112(6): 2824.

    [9] Alasdair Mac Arthur, Tim Malthus. International Journal of Remote Sensing, 2012, 33(16): 5214.

    [10] Feret J, Francois C, Asner G P, et al. Remote Sensing of Environment, 2008, 112(6): 3030.

    [11] Saltelli A, Tarantola S, Chan K. Technometrics, 1999, 41(1): 39.

    [12] LI Yan, HUANG Chun-lin, LU Ling, et al(李 艷, 黃春林, 盧 玲, 等). Remote Sensing Technology and Application(遙感技術(shù)與應(yīng)用), 2014, (5): 719.

    [13] Moorthy I, Miller J R, Zarco-Tejada P J, et al. IEEE International Geoscience & International Geoscience & Remote Sensing Symposium, 2004, (3): 1676.

    [14] XIAO Yan-fang, ZHOU De-min, GONG Hui-li, et al(肖艷芳, 周德民, 宮輝力, 等). Journal of Remote Sensing(遙感學報), 2015, (3): 368.

    *Corresponding author

    Simulation of Needle Reflectance Spectrum and Sensitivity Analysis of Biochemical Parameters ofPinusYunnanensisin Different Healthy Status

    LIN Qi-nan, HUANG Hua-guo*, CHEN Ling, YU Lin-feng, HUANG Kan

    The Key Laboratory for Silviculture and Conservation of Ministry of Education, Beijing Forestry University, Beijing 100083, China

    The sensitivity of biochemical effects on leaf reflectance is vital for retieving biochemical parameters with remote sensing. In this study, the chlorophyll and water absorption coefficients of the commonly used model LIBERTY (leaf incorporation biochemistry exhibiting reflectance and transmittance yields) were calibrated using field measured needle spectral reflectance curves based on a look up table (LUT) method. A novel spectra reflectance fitting method were presented by involving a new index (named as yellow index, YI), which could obviously improve the fitting accuracy of Pinus yunnanensis reflection spectrum at highly-stressed status. As a global sensitivity analysis method, the EFAST (extended Fourier amplitued sensitivity test) was implemented to quantitatively assess the sensitivity of biochemical parameters on needle reflectance. Results show that: (1) the reflectauce spectrum of healthy needles (R2=0.999, RMSE<0.01), slightly stressed needles (R2=0.991, RMSE<0.02) and moderately stressed needles (R2=0.992, RMSE<0.03) are simulated fairly well by calibrated LIBERTY model which has less potential in fitting the reflectance spectrum of seriously stressed needles (R2=0.803, RMSE>0.1). (2) the reflectance spectrum of seriously stressed needles can be successfully simulated by our proposed spectrum reflectance fitting method (R2=0.991, RMSE<0.03), because YI can quantitatively describe different degrees of stress, and (3) the sensitivity of leaf reflectance tochlorophyllandwaterparameters decreases with the degree of stress; while the sensitivity to other biochemical parameters is increasing, which includ baseline absorption, albino absorption, Lignin and Cellulose content, and nitrogen content, increases with the stress degree. Needle reflectance spectrum also have sensitivive bands for these parameters. For example, the albino absorption have a significant effect on needle reflectance in 505~565 and 705~850 nm). In addition,Albinoabsorptionandchlorophyllalso have significant effects on needle reflectance in visible region for seriously stressed needles, which indicates that the prior knowledge of thealbinoabsorptionlevel can help obtain the valid inversion result ofchlorophyllcontent.

    LIBERTY; Needle reflectance; Sensitivity analysis;Pinusyunnanensis; EFAST

    Jan. 24, 2016; accepted Apr. 21, 2016)

    2016-01-24,

    2016-04-21

    林業(yè)公益性行業(yè)科研專項(201404401)和國家自然科學基金項目(41571332)資助

    林起楠, 1991年生, 北京林業(yè)大學省部共建森林培育與保護教育部重點實驗室碩士研究生 e-mail: 1056446928@qq.com *通訊聯(lián)系人 e-mail: huaguo_huang@bjfu.edu.cn

    TP79, S763

    A

    10.3964/j.issn.1000-0593(2016)08-2538-08

    猜你喜歡
    云南松針葉白化
    更 正
    白化黃喉擬水龜人工培育研究①
    云南松正交膠合木膠合工藝探索?
    最嚴重白化
    風雨海棠
    云南松天然更新幼苗和幼樹期的生長動態(tài)分析
    云南松及其變種栽培與適應(yīng)性研究進展
    廣西細葉云南松森林資源動態(tài)分析
    今年1季度德國針葉材出口量增長7%
    加拿大針葉鋸材繼續(xù)以美國南部為中心
    别揉我奶头~嗯~啊~动态视频| 色在线成人网| 成年人免费黄色播放视频| 久久亚洲真实| 久久久久视频综合| 咕卡用的链子| 国产一区二区三区视频了| 一级a爱片免费观看的视频| 亚洲熟妇熟女久久| 在线观看日韩欧美| 日日爽夜夜爽网站| 精品电影一区二区在线| 成人亚洲精品一区在线观看| 三级毛片av免费| 50天的宝宝边吃奶边哭怎么回事| 999精品在线视频| 热99久久久久精品小说推荐| 亚洲五月天丁香| 欧美一级毛片孕妇| 色老头精品视频在线观看| 午夜影院日韩av| 久久精品人人爽人人爽视色| 午夜精品久久久久久毛片777| √禁漫天堂资源中文www| 国产av又大| 国产一区在线观看成人免费| 免费黄频网站在线观看国产| 男女床上黄色一级片免费看| 国产色视频综合| 久久午夜综合久久蜜桃| 欧美日韩黄片免| 国产一卡二卡三卡精品| 制服诱惑二区| 美女福利国产在线| 亚洲少妇的诱惑av| 90打野战视频偷拍视频| 黄色片一级片一级黄色片| 国产精品电影一区二区三区 | 亚洲精品国产精品久久久不卡| 国产色视频综合| 久久人妻av系列| 午夜两性在线视频| 日韩精品免费视频一区二区三区| 在线观看免费视频日本深夜| 如日韩欧美国产精品一区二区三区| 亚洲七黄色美女视频| 午夜福利在线免费观看网站| 午夜福利影视在线免费观看| 美女福利国产在线| 老司机深夜福利视频在线观看| 建设人人有责人人尽责人人享有的| 悠悠久久av| 涩涩av久久男人的天堂| 黄色女人牲交| 一个人免费在线观看的高清视频| 捣出白浆h1v1| 国产男靠女视频免费网站| 亚洲第一av免费看| 热re99久久精品国产66热6| 少妇粗大呻吟视频| av天堂在线播放| 悠悠久久av| 午夜福利视频在线观看免费| 欧美黑人精品巨大| 91麻豆精品激情在线观看国产 | 免费在线观看视频国产中文字幕亚洲| 黑人欧美特级aaaaaa片| 99久久99久久久精品蜜桃| 亚洲精品一卡2卡三卡4卡5卡| ponron亚洲| 精品国产美女av久久久久小说| 亚洲国产精品一区二区三区在线| 亚洲第一欧美日韩一区二区三区| 777久久人妻少妇嫩草av网站| 一a级毛片在线观看| 狠狠狠狠99中文字幕| 精品久久蜜臀av无| 国产亚洲av高清不卡| 国产免费av片在线观看野外av| 日韩人妻精品一区2区三区| svipshipincom国产片| 侵犯人妻中文字幕一二三四区| 国产欧美日韩一区二区三| 亚洲av日韩在线播放| 午夜福利在线观看吧| 成人国产一区最新在线观看| 免费在线观看影片大全网站| 老司机靠b影院| 国产精品1区2区在线观看. | 免费观看精品视频网站| 日韩有码中文字幕| 亚洲aⅴ乱码一区二区在线播放 | 久久香蕉精品热| 精品一区二区三区视频在线观看免费 | 人人妻人人澡人人看| 久久久国产一区二区| 日日爽夜夜爽网站| 亚洲欧美精品综合一区二区三区| 老熟女久久久| 午夜福利一区二区在线看| 91精品三级在线观看| 国产精品自产拍在线观看55亚洲 | 欧美日韩精品网址| 两个人免费观看高清视频| 免费看十八禁软件| 我的亚洲天堂| 男人舔女人的私密视频| 天天操日日干夜夜撸| av一本久久久久| 999久久久精品免费观看国产| 亚洲av成人av| 亚洲精品乱久久久久久| 麻豆乱淫一区二区| 国产一区二区三区在线臀色熟女 | 一本大道久久a久久精品| 精品国产一区二区三区四区第35| 免费看十八禁软件| 国产成人精品久久二区二区免费| 日韩精品免费视频一区二区三区| 别揉我奶头~嗯~啊~动态视频| 丝袜美足系列| 国产又色又爽无遮挡免费看| 中文欧美无线码| 黄色毛片三级朝国网站| 久久国产精品男人的天堂亚洲| av超薄肉色丝袜交足视频| 国产精品一区二区在线不卡| 午夜视频精品福利| 亚洲一码二码三码区别大吗| 国产高清videossex| 动漫黄色视频在线观看| 精品国产一区二区三区久久久樱花| 久久久国产成人精品二区 | 亚洲片人在线观看| 99精品久久久久人妻精品| 在线观看舔阴道视频| 国产男靠女视频免费网站| 岛国毛片在线播放| 麻豆av在线久日| 中文字幕色久视频| 看免费av毛片| 在线观看舔阴道视频| 自线自在国产av| 国产免费av片在线观看野外av| 精品久久久久久久毛片微露脸| 在线观看午夜福利视频| 欧美一级毛片孕妇| 亚洲中文日韩欧美视频| 国产精品久久久久久精品古装| 欧美在线黄色| 露出奶头的视频| 久久九九热精品免费| 亚洲avbb在线观看| 香蕉久久夜色| 黄色 视频免费看| 1024香蕉在线观看| 国产伦人伦偷精品视频| 国产精品一区二区在线不卡| 亚洲国产中文字幕在线视频| 黄色毛片三级朝国网站| 亚洲一区二区三区不卡视频| 一进一出抽搐gif免费好疼 | 91成人精品电影| 久久久国产精品麻豆| 男女免费视频国产| av超薄肉色丝袜交足视频| 国产精品秋霞免费鲁丝片| 操美女的视频在线观看| 日日摸夜夜添夜夜添小说| 亚洲精品中文字幕在线视频| 在线免费观看的www视频| 亚洲av日韩在线播放| 制服人妻中文乱码| 国产又爽黄色视频| 18禁黄网站禁片午夜丰满| 老司机亚洲免费影院| 人妻一区二区av| 精品无人区乱码1区二区| 国产精品乱码一区二三区的特点 | 日韩 欧美 亚洲 中文字幕| 丝袜人妻中文字幕| 久久精品亚洲精品国产色婷小说| www.999成人在线观看| 国产精品亚洲av一区麻豆| 亚洲专区国产一区二区| 中文字幕人妻丝袜制服| xxxhd国产人妻xxx| 国产真人三级小视频在线观看| a在线观看视频网站| 亚洲国产欧美网| 午夜福利一区二区在线看| 女警被强在线播放| 国产精品国产高清国产av | 好看av亚洲va欧美ⅴa在| 人妻 亚洲 视频| 视频区图区小说| 一边摸一边抽搐一进一出视频| 亚洲精品美女久久久久99蜜臀| 在线av久久热| 久久久久久人人人人人| 亚洲av熟女| xxxhd国产人妻xxx| 热99久久久久精品小说推荐| 亚洲人成伊人成综合网2020| 欧美亚洲日本最大视频资源| 麻豆av在线久日| 看黄色毛片网站| 人妻 亚洲 视频| 午夜老司机福利片| 999久久久精品免费观看国产| 中文字幕另类日韩欧美亚洲嫩草| 久久久久久久久免费视频了| 村上凉子中文字幕在线| 99精品久久久久人妻精品| 欧美不卡视频在线免费观看 | 亚洲人成伊人成综合网2020| 国产视频一区二区在线看| 一个人免费在线观看的高清视频| 人妻丰满熟妇av一区二区三区 | 午夜老司机福利片| 成人特级黄色片久久久久久久| 国产男女超爽视频在线观看| 看免费av毛片| 欧美午夜高清在线| 色综合欧美亚洲国产小说| 国产成人欧美| 亚洲一区二区三区欧美精品| 欧美精品啪啪一区二区三区| 亚洲性夜色夜夜综合| 欧美中文综合在线视频| videos熟女内射| 91国产中文字幕| 中文字幕精品免费在线观看视频| 亚洲熟妇中文字幕五十中出 | 国产乱人伦免费视频| 亚洲精品粉嫩美女一区| 成人18禁高潮啪啪吃奶动态图| 久久国产乱子伦精品免费另类| 在线观看免费日韩欧美大片| 久热爱精品视频在线9| 国产精品一区二区在线观看99| 欧美大码av| 日韩成人在线观看一区二区三区| 女警被强在线播放| 国产精品久久久久久精品古装| 在线看a的网站| 一级毛片高清免费大全| 日本黄色日本黄色录像| av天堂久久9| 丁香欧美五月| 99久久精品国产亚洲精品| 亚洲中文字幕日韩| 最近最新中文字幕大全免费视频| 国产在视频线精品| 亚洲精品国产精品久久久不卡| 成人18禁在线播放| 国产又爽黄色视频| 女人被狂操c到高潮| 亚洲精品国产区一区二| 国产黄色免费在线视频| 成年人午夜在线观看视频| 久久久国产精品麻豆| 最新美女视频免费是黄的| 国产xxxxx性猛交| bbb黄色大片| 天堂俺去俺来也www色官网| 性色av乱码一区二区三区2| 窝窝影院91人妻| av视频免费观看在线观看| 国产精品亚洲av一区麻豆| 纯流量卡能插随身wifi吗| 亚洲精华国产精华精| 好看av亚洲va欧美ⅴa在| 亚洲精品美女久久av网站| 日韩欧美在线二视频 | 国产高清国产精品国产三级| 久久久精品免费免费高清| 亚洲免费av在线视频| 亚洲国产欧美日韩在线播放| 高清毛片免费观看视频网站 | 日本黄色视频三级网站网址 | 女人被躁到高潮嗷嗷叫费观| 中国美女看黄片| 老司机深夜福利视频在线观看| 女人久久www免费人成看片| 国产精品欧美亚洲77777| 久久国产乱子伦精品免费另类| 国产精品美女特级片免费视频播放器 | 欧洲精品卡2卡3卡4卡5卡区| 一级a爱片免费观看的视频| 精品无人区乱码1区二区| 欧美日韩黄片免| 1024视频免费在线观看| 欧美精品亚洲一区二区| 国产成人精品久久二区二区91| 午夜91福利影院| 日韩三级视频一区二区三区| 我的亚洲天堂| 十八禁网站免费在线| 搡老乐熟女国产| 在线观看免费视频网站a站| 国产精品亚洲av一区麻豆| 亚洲 国产 在线| 国产蜜桃级精品一区二区三区 | 欧美黑人精品巨大| 国产精品免费大片| 亚洲国产看品久久| 在线av久久热| 亚洲欧美色中文字幕在线| 无人区码免费观看不卡| 欧美亚洲日本最大视频资源| 又大又爽又粗| 变态另类成人亚洲欧美熟女 | 国产av精品麻豆| 欧美最黄视频在线播放免费 | 久久中文字幕一级| 咕卡用的链子| 国产有黄有色有爽视频| 好看av亚洲va欧美ⅴa在| 五月开心婷婷网| 精品少妇久久久久久888优播| 日韩欧美一区二区三区在线观看 | av中文乱码字幕在线| 国产精品成人在线| 国产一区二区三区视频了| 十八禁人妻一区二区| 久久精品亚洲精品国产色婷小说| 国产精品av久久久久免费| 午夜日韩欧美国产| 精品国产一区二区三区久久久樱花| 少妇猛男粗大的猛烈进出视频| 欧美黑人欧美精品刺激| 50天的宝宝边吃奶边哭怎么回事| 老汉色av国产亚洲站长工具| 免费观看精品视频网站| 久久香蕉国产精品| 中文字幕高清在线视频| 91成年电影在线观看| 国产1区2区3区精品| 久久久久精品人妻al黑| 免费观看精品视频网站| 国产三级黄色录像| 亚洲av成人一区二区三| 免费高清在线观看日韩| 久久精品亚洲精品国产色婷小说| 久久人妻熟女aⅴ| 亚洲专区字幕在线| 最近最新中文字幕大全免费视频| 国产精品国产高清国产av | 国产精品久久久久久精品古装| 大香蕉久久成人网| 黄色 视频免费看| 国产欧美日韩一区二区精品| 纯流量卡能插随身wifi吗| 亚洲一卡2卡3卡4卡5卡精品中文| 18禁观看日本| 精品久久久久久电影网| 伊人久久大香线蕉亚洲五| 在线观看免费午夜福利视频| 欧美日韩亚洲高清精品| 国产精品秋霞免费鲁丝片| 黑人巨大精品欧美一区二区蜜桃| 极品人妻少妇av视频| 人人妻,人人澡人人爽秒播| 午夜老司机福利片| 国产亚洲精品第一综合不卡| 狂野欧美激情性xxxx| 国产精品久久久av美女十八| av网站在线播放免费| 99re在线观看精品视频| av网站在线播放免费| 美女午夜性视频免费| av网站在线播放免费| 国产精品香港三级国产av潘金莲| 亚洲精品国产一区二区精华液| 国产xxxxx性猛交| 久久久国产一区二区| avwww免费| 国产精品 欧美亚洲| 久热这里只有精品99| svipshipincom国产片| 丰满的人妻完整版| 久久人人爽av亚洲精品天堂| 天堂动漫精品| 一级毛片高清免费大全| 亚洲九九香蕉| a级毛片黄视频| 男女下面插进去视频免费观看| 欧美不卡视频在线免费观看 | 久久国产精品人妻蜜桃| 99国产精品免费福利视频| 人妻丰满熟妇av一区二区三区 | 久久人妻av系列| 99精品欧美一区二区三区四区| 亚洲欧美日韩高清在线视频| 久久久久久久午夜电影 | 午夜免费成人在线视频| 欧美黄色淫秽网站| 亚洲精品久久成人aⅴ小说| 色尼玛亚洲综合影院| 国产精品久久久人人做人人爽| 亚洲精品国产色婷婷电影| 热99久久久久精品小说推荐| 高潮久久久久久久久久久不卡| 亚洲熟妇熟女久久| 欧美日韩国产mv在线观看视频| 国产av精品麻豆| 国产精品久久久久成人av| 国产成人精品无人区| 黄色成人免费大全| 久久九九热精品免费| 久久久久久人人人人人| 久久久久久亚洲精品国产蜜桃av| avwww免费| 中文字幕最新亚洲高清| 捣出白浆h1v1| 咕卡用的链子| 精品福利永久在线观看| 久久久精品区二区三区| 国产熟女午夜一区二区三区| 十八禁高潮呻吟视频| 久久草成人影院| 亚洲片人在线观看| 亚洲在线自拍视频| 午夜福利在线免费观看网站| 欧美精品啪啪一区二区三区| 18禁裸乳无遮挡免费网站照片 | 超色免费av| 少妇粗大呻吟视频| 最新在线观看一区二区三区| 久久精品熟女亚洲av麻豆精品| svipshipincom国产片| 黄片小视频在线播放| 黄色怎么调成土黄色| 欧美激情极品国产一区二区三区| 9191精品国产免费久久| 交换朋友夫妻互换小说| 欧美日韩黄片免| 午夜老司机福利片| 精品一区二区三卡| 亚洲欧洲精品一区二区精品久久久| 老熟妇乱子伦视频在线观看| 日韩大码丰满熟妇| 12—13女人毛片做爰片一| 精品国产一区二区久久| 午夜福利欧美成人| 最近最新中文字幕大全电影3 | 国产在线精品亚洲第一网站| 亚洲av日韩精品久久久久久密| 亚洲专区字幕在线| 免费看十八禁软件| 天天躁夜夜躁狠狠躁躁| 欧美成人免费av一区二区三区 | 国产亚洲欧美精品永久| 亚洲,欧美精品.| 国产1区2区3区精品| 欧美日韩乱码在线| 亚洲欧美激情在线| 人人妻人人澡人人看| tube8黄色片| 黄色丝袜av网址大全| 久久久水蜜桃国产精品网| 91在线观看av| 色精品久久人妻99蜜桃| ponron亚洲| 国产单亲对白刺激| 欧美激情 高清一区二区三区| 日日摸夜夜添夜夜添小说| 91av网站免费观看| 在线观看免费视频日本深夜| 人人妻人人澡人人看| 国产又爽黄色视频| 亚洲在线自拍视频| 欧美日韩亚洲综合一区二区三区_| 老熟妇乱子伦视频在线观看| 麻豆国产av国片精品| 亚洲伊人色综图| 精品午夜福利视频在线观看一区| 人人妻人人爽人人添夜夜欢视频| 亚洲欧美一区二区三区久久| 老鸭窝网址在线观看| 不卡一级毛片| e午夜精品久久久久久久| 精品一区二区三卡| 亚洲精品久久午夜乱码| 国产亚洲精品久久久久5区| 麻豆成人av在线观看| 99热国产这里只有精品6| 天天影视国产精品| 最近最新中文字幕大全电影3 | 高清视频免费观看一区二区| 欧美乱码精品一区二区三区| a级毛片在线看网站| 老熟妇仑乱视频hdxx| 最新在线观看一区二区三区| 国产一卡二卡三卡精品| 黑人操中国人逼视频| 老鸭窝网址在线观看| 每晚都被弄得嗷嗷叫到高潮| 久久久久国产一级毛片高清牌| 亚洲精品一二三| 99国产精品99久久久久| 久久久久久久久久久久大奶| 18禁观看日本| √禁漫天堂资源中文www| 国产欧美亚洲国产| 亚洲一区中文字幕在线| 中文字幕人妻熟女乱码| 一a级毛片在线观看| 精品少妇久久久久久888优播| 少妇裸体淫交视频免费看高清 | 日韩免费av在线播放| 自拍欧美九色日韩亚洲蝌蚪91| 美女高潮喷水抽搐中文字幕| 久久精品亚洲精品国产色婷小说| 亚洲五月婷婷丁香| 午夜影院日韩av| 国产精品一区二区在线不卡| 99精国产麻豆久久婷婷| 一边摸一边做爽爽视频免费| 国产成人免费观看mmmm| 亚洲,欧美精品.| 女人爽到高潮嗷嗷叫在线视频| 亚洲五月天丁香| 韩国精品一区二区三区| 高清在线国产一区| 国产精品自产拍在线观看55亚洲 | 国产成人一区二区三区免费视频网站| 成年动漫av网址| 中文字幕av电影在线播放| 69av精品久久久久久| 欧美日韩亚洲高清精品| 国产精品免费大片| 视频区欧美日本亚洲| 手机成人av网站| av福利片在线| 久热这里只有精品99| 99国产精品99久久久久| 亚洲精品自拍成人| 新久久久久国产一级毛片| 91大片在线观看| 精品国产一区二区三区久久久樱花| 欧美成人午夜精品| 国产亚洲欧美98| 人妻一区二区av| 国产精品久久久久成人av| 制服诱惑二区| 国产不卡一卡二| 亚洲 欧美一区二区三区| 交换朋友夫妻互换小说| 嫩草影视91久久| 国产高清videossex| 国产av精品麻豆| 动漫黄色视频在线观看| 精品久久蜜臀av无| 亚洲精品国产精品久久久不卡| 亚洲色图av天堂| 99热国产这里只有精品6| 人妻 亚洲 视频| 黑人操中国人逼视频| 男人的好看免费观看在线视频 | 亚洲aⅴ乱码一区二区在线播放 | 中文字幕人妻丝袜制服| 亚洲精品国产区一区二| 高清av免费在线| 五月开心婷婷网| 精品少妇一区二区三区视频日本电影| 两人在一起打扑克的视频| 人人妻人人爽人人添夜夜欢视频| 色婷婷久久久亚洲欧美| 女警被强在线播放| 久久精品国产清高在天天线| 91九色精品人成在线观看| 两性夫妻黄色片| 久久国产精品男人的天堂亚洲| 一级毛片精品| 1024视频免费在线观看| 日韩成人在线观看一区二区三区| 脱女人内裤的视频| 亚洲av欧美aⅴ国产| 成人国产一区最新在线观看| 久久午夜综合久久蜜桃| 亚洲av第一区精品v没综合| 三上悠亚av全集在线观看| 黄色丝袜av网址大全| 91av网站免费观看| 一进一出抽搐gif免费好疼 | 动漫黄色视频在线观看| av网站在线播放免费| 国产成人精品在线电影| 在线观看www视频免费| 欧美日韩亚洲国产一区二区在线观看 | 女警被强在线播放| 国产91精品成人一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久久国产一级毛片高清牌| 国产午夜精品久久久久久| 亚洲五月色婷婷综合| 久久精品国产99精品国产亚洲性色 | 国产免费男女视频| 在线观看日韩欧美| 91麻豆精品激情在线观看国产 | 在线观看免费视频网站a站| 99久久99久久久精品蜜桃| 香蕉国产在线看| 人妻 亚洲 视频| 久久精品亚洲精品国产色婷小说| 亚洲人成电影免费在线| 18禁国产床啪视频网站| 欧美精品啪啪一区二区三区| 精品一区二区三区四区五区乱码| 亚洲久久久国产精品| 一二三四在线观看免费中文在| 亚洲成国产人片在线观看| 五月开心婷婷网|