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

    基于回歸分析的玉米冠層葉綠素含量高光譜反演分析

    2021-09-08 07:45:02趙占輝張叢志張佳寶吳運金張宏敏魯春陽
    中國農(nóng)學(xué)通報 2021年20期
    關(guān)鍵詞:診斷模型冠層反射率

    趙占輝,張叢志,張佳寶,吳運金,張宏敏,魯春陽

    (1河南城建學(xué)院測繪與城市空間信息學(xué)院,河南平頂山 467036;2封丘農(nóng)田生態(tài)系統(tǒng)國家試驗站/土壤與農(nóng)業(yè)可持續(xù)發(fā)展國家重點實驗室/中國科學(xué)院南京土壤研究所,南京 210008;3生態(tài)環(huán)境部南京環(huán)境科學(xué)研究所,南京 210042)

    0 引言

    葉綠素是高等植物和其他所有能進(jìn)行光合作用的生物體含有的一類綠色色素,其含量高低能夠反映植物光合速率的強弱、氮利用效率、生理脅迫和生長狀況等信息[1-2]。當(dāng)前農(nóng)作物葉綠素含量測定主要采用分光光度法、活體葉綠素儀法和光譜法[3]。分光光度法通過化學(xué)方法提取葉片中的葉綠素,以此測定葉片葉綠素含量,測定結(jié)果穩(wěn)定可靠,但提取速度較慢、易損壞植株葉面,容易影響農(nóng)作物生長。根據(jù)工作原理不同,活體葉綠素儀法可分為透射型活體葉綠素儀法和反射型活體葉綠素儀法,通過測定紅光(650 nm)與紅外光(940 nm)被葉片吸收比率,測量植物的葉綠素相對含量,該方法是一種速測法,測定結(jié)果穩(wěn)定可靠,不會對植被葉片造成破壞[4]。光譜法主要光譜信息技術(shù)進(jìn)行建模分析,建立植物葉片葉綠素含量化學(xué)值與光譜數(shù)據(jù)之間的定量關(guān)系,間接獲取植物葉片葉綠素含量信息。

    由于傳統(tǒng)的化學(xué)提取方法具有非常大的機械破壞性且不可恢復(fù),其工作量巨大,無法進(jìn)行大規(guī)模實時監(jiān)測,活體葉綠素儀法雖能大大提高測量工作效率,但依舊受樣本數(shù)量限制,無法進(jìn)行大面積的植物葉綠素含量變化實時監(jiān)測。相反,隨著遙感技術(shù)快速發(fā)展,基于遙感數(shù)據(jù)的光譜分析法依靠少量樣本就能夠快速、無損的獲取大面積農(nóng)作物葉綠素含量,對快速監(jiān)測作物光合作用強弱、生長發(fā)育狀況等具有重要意義,關(guān)于農(nóng)作物葉綠素含量與光譜關(guān)系的研究也已經(jīng)成為近年來農(nóng)業(yè)領(lǐng)域研究的熱點。

    無人機搭載的高光譜儀具有光譜分辨率高、波段多、波段連續(xù)性強等優(yōu)點,且能夠快速獲取納米級的光譜信息,利用該技術(shù)能直接對植被冠層微弱光譜進(jìn)行定量分析,在農(nóng)作物長勢監(jiān)測方面表現(xiàn)出強大的優(yōu)勢[2,5-6]。因此,根據(jù)作物冠層的光譜反射率與葉綠素含量建立反演模型,快速獲取空間尺度作物葉綠素含量分布特征,對評估作物光合作用強弱、生長發(fā)育狀況具有重要意義。早在20世紀(jì)Hinzman等[7]學(xué)者就已經(jīng)提出利用作物冠層反射光譜監(jiān)測作物葉綠素含量具有可行性,隨后,國內(nèi)外學(xué)者[8-10]先后在實驗室利用高光譜數(shù)據(jù)對植物葉片葉綠素含量進(jìn)行了估算,基本實現(xiàn)了作物葉片無損監(jiān)測與葉綠素分布可視化。梁爽等[5]、鄭濤等[9]也通過篩選敏感波段,建立了葉綠素診斷模型,驗證了作物葉片葉綠素含量與光譜特征具有密切關(guān)系。肖武等[11]研究還發(fā)現(xiàn)玉米冠層葉綠素含量對紅邊波段的響應(yīng)同樣敏感,以此為基礎(chǔ)建立葉綠素含量遙感反演模型,并應(yīng)用于土地?fù)p毀監(jiān)測與評價,估測精度達(dá)到了83.4%。

    目前,由于所獲取的影像數(shù)據(jù)波段數(shù)量較少、光譜信息有限,而且光譜診斷波段及參數(shù)的范圍較窄,導(dǎo)致作物葉綠素含量高光譜反演模型估測精度普遍低于85%,難以滿足當(dāng)前數(shù)字農(nóng)業(yè)的需求,有待進(jìn)一步深入研究更高精度的作物冠層高光譜估測模型。為此,本研究以中國科學(xué)院封丘農(nóng)業(yè)生態(tài)實驗站長期定位氮肥梯度試驗田為研究對象,研究了玉米冠層光譜反射率與作物葉綠素含量之間的關(guān)系,建立了基于高光譜反射率的玉米冠層葉綠素含量反演模型,并對模型精度、可靠性進(jìn)行了驗證分析,以期為推進(jìn)高光譜遙感技術(shù)在農(nóng)業(yè)中的應(yīng)用提供數(shù)據(jù)支撐。

    1 材料與方法

    1.1 研究區(qū)概況

    試驗區(qū)位于河南省封丘縣(35°01'N,114°24'E),在中國科學(xué)院封丘農(nóng)田生態(tài)系統(tǒng)國家試驗站內(nèi)進(jìn)行,該試驗站地處黃淮海平原,屬于半干旱半濕潤季風(fēng)氣候,多年平均降雨量為615 mm,且60%~90%分布在6—10月,年平均氣溫為13.9℃。定位試驗自2008 年一直延續(xù)至今,夏玉米品種為‘鄭單958’,共設(shè)7個氮肥水平:0、30、70、110、150、190、240 kg(N)hm2/a,磷肥按氮肥施入量的40%施入,不施鉀肥,玉米季磷肥作為基肥一次施入,氮肥肥分別在苗期、灌漿期各施入40%、60%。每個處理設(shè)置3次重復(fù),共計21個小區(qū),每個小區(qū)面積約28 m2,試驗當(dāng)天在每個小區(qū)采集3 個樣本,共計63個樣本,其中前2次重復(fù)小區(qū)中42個樣本數(shù)據(jù)設(shè)為數(shù)據(jù)集Ⅰ,用于模型校正,第3次重復(fù)小區(qū)中21個樣本數(shù)據(jù)設(shè)為數(shù)據(jù)集Ⅱ,用于模型驗證(如圖1)。

    圖1 基于高光譜數(shù)據(jù)合成的實驗區(qū)真彩色影像圖

    1.2 數(shù)據(jù)獲取與處理方法

    無人機搭載的傳感器為Resonon PIKA L 高光譜成像儀,PIKA L 采用非凝視成像原理,拍攝方式為線掃式,光譜范圍為400~1000 nm,波段數(shù)200,光譜分辨率3 nm。無人機飛行高80 m,設(shè)定航速1 m/s。傳感器鏡頭焦距17 mm,視場角17.6°,鏡頭垂直向下,地面分辨率5 cm,航拍面積0.25 km2覆蓋整個研究區(qū)。

    本研究以小型系統(tǒng)性試驗田為研究靶區(qū),獲取的高光譜影像的預(yù)處理流程如下:(1)由于原始影像數(shù)據(jù)為信號DN值,無單位,需要通過輻射定標(biāo)文件將原始影像轉(zhuǎn)換為輻射率圖像,本研究利用Spectronon Pro軟件將高光譜影像輻射亮度轉(zhuǎn)換成反射率數(shù)據(jù);(2)在Spectronon Pro軟件中根據(jù)飛行高度、視場角等參數(shù)信息對高光譜影像進(jìn)行粗校正;(3)基于差分GPS(Trimble GeoXH,美國天寶公司)獲得的地面靶標(biāo)點坐標(biāo)數(shù)據(jù),利用ArcGIS 10.3軟件對高光譜影像進(jìn)行手動精幾何校正;(4)最后,在Spectronon Pro軟件中對獲得的影像進(jìn)行拼接,利用ENVI 5.1軟件裁切生成研究區(qū)高光譜影像。

    根據(jù)葉綠素地面測量對應(yīng)的樣點位置構(gòu)建感興趣區(qū)域(region of interest,ROI),采集該樣點玉米葉片的光譜反射率,得到各樣點的反射率光譜數(shù)據(jù)。

    1.3 玉米葉片SPAD值測量

    為同步、快速、無損獲得玉米葉片葉綠素含量,采用JN-4N型葉綠素測定儀,測得葉片SPAD值。JN-4N型手持葉綠素儀2 個LED 光源發(fā)射2 種光,一種是紅光(峰波長650 nm),一種是紅外線(940 nm),2 種光穿透葉片,打到接收器上,光信號轉(zhuǎn)換成模擬信號,模擬信號被放大器放大,由模擬/數(shù)字轉(zhuǎn)換器轉(zhuǎn)換成數(shù)字信號,數(shù)字信號被微處理器處理后計算出SPAD 值。本研究在每個小區(qū)選取3個樣點,63個樣點,同步測得玉米葉片SPAD值。每個樣點的玉米選擇3片葉子測量,每片葉子不同部位測3個SPAD值,取平均值作為該片葉子的SPAD值,并記錄樣點位置。

    1.4 反演模型構(gòu)建與檢驗

    過去研究發(fā)現(xiàn),反射光譜在可見光及中遠(yuǎn)紅外區(qū)域是診斷作物葉片葉綠素含量的敏感波段,葉綠素作物葉片的重要賦色成分,其含量高低對反射光譜產(chǎn)生一定影響[9,12-13],研究區(qū)不同葉綠素含量的光譜反射曲線形態(tài)特征基本一致(圖2)。在可見光波段400~650 nm 范圍內(nèi)光譜曲線呈先上升后下降趨勢,650~760 nm范圍呈先下降后快速上升趨勢,具有明顯的反射峰和吸收谷;在近紅波段(760~900 nm)范圍趨于上升趨勢,且不同葉綠素含量的光譜反射率較為容易區(qū)分。由于原始光譜數(shù)據(jù)是離散的數(shù)據(jù),本研究采用如下公式對光譜儀采集原始光譜數(shù)據(jù)進(jìn)行一階微分計算[14]。選取夏玉米葉片光譜“三邊”特征的參數(shù)、反射峰與吸收帶位置、植被光譜指數(shù)等17個前人構(gòu)建的高光譜特征參數(shù)(表1);以光譜特征參數(shù)為自變量利用回歸法建立葉綠素含量的估算模型[12,15]。將試驗獲得的數(shù)據(jù)集Ⅰ用于模型校正,數(shù)據(jù)集Ⅱ用于模型驗證。通過對比各模型的建模決定系數(shù)R2、極值(Max、Min)、極差(Range)、變異系數(shù)(Coefficient of Variation)、平均絕對誤差率(Mean Absolute Percentage Error,MAPE)和殘差均方根(Root Mean Squared Residual,RMSR)篩選出最優(yōu)模型用以研究區(qū)春小麥的葉綠素含量的反演。R2用以估測模型的穩(wěn)定程度,越接近于1 說明模型的穩(wěn)定性越好,最大值(Max)、最小值(Min)、極差(Range)和變異系數(shù)(Coefficient of Variation,C.V.)能反映不同診斷模型對空間變異性的表達(dá)能力,平均絕對誤差率(Mean Absolute Percentage Error,MAPE)和殘差均方根(Root Mean Squared Residual)則從另一角度表征診斷模型的估測結(jié)果與“真值”的擬合程度,其計算公式如式(1)~(2)所示。

    表1 反射光譜特征參數(shù)

    式中,xi為采樣調(diào)查法土壤有機碳密度測定值,xi為土壤有機碳密度估測值,n為樣本數(shù)量。

    本研究利用ENVI 5.1 軟件進(jìn)行光譜數(shù)據(jù)讀取與波段運算,在SPSS 22.0 軟件中進(jìn)行模型構(gòu)建與檢驗,采用ArcGIS 10.3軟件制圖。

    2 結(jié)果與分析

    2.1 不同SPAD值的玉米冠層葉片反射光譜特征

    從高光譜影像中提取SPAD測樣點位置的平均反射光譜特征曲線,得到玉米葉片反射光譜對7 種不同SPAD 值水平的響應(yīng),并在玉米地臨近地塊獲取裸地的平均反射光譜特征曲線(圖2)??傮w上,不同施氮水平下葉片光譜曲線形狀相似,可見光區(qū)域(400~700 nm)植物冠層的光譜反射率均較低,其光譜曲線存在一個反射峰和一個吸收谷,即550 nm左右區(qū)域的綠光與690 nm左右區(qū)域的紅光,該波長范圍內(nèi)玉米葉片SPAD高值區(qū)的光譜反射率低于SPAD低值區(qū)。與可見光區(qū)域相比,近紅外區(qū)域光譜反射率大幅上升,特別是在690~750 nm之間的光譜反射率急劇上升,形成植被光譜的最重要特征——紅邊。不同SPAD值的玉米冠層葉片在750~900 nm之間的光譜反射率差異顯著,高量施氮處理冠層葉片光譜反射率顯著高于低量施氮處理。裸地光譜曲線隨波長增加均勻增加,符合該地物的典型特征,說明本文獲取的高光譜數(shù)據(jù)具有一定的可靠性。

    圖2 玉米與裸地光譜特征曲線

    2.2 玉米冠層葉片SPAD值診斷模型構(gòu)建及精度分析

    本研究分別將數(shù)據(jù)集Ⅰ中42個訓(xùn)練樣本SPAD(R)值進(jìn)行了Log(R)、Log(1/R)、1/R、SQRT(R)、SQRT(1/R)5種數(shù)學(xué)變換,將17個光譜參數(shù)根據(jù)其特性劃分為4組,分別為反射曲線“三邊”斜率及其位置、葉綠素的反射峰與吸收帶位置、反射曲線“三邊”面積、光譜指數(shù)與葉綠素含量相關(guān)的植被指數(shù)。為了簡化模型參數(shù),在建模前對各參數(shù)進(jìn)行篩選,篩選出與SPAD值或其變換值關(guān)系緊密的的光譜參數(shù)作為診斷模型備選參數(shù),具體為:根據(jù)相關(guān)系數(shù)大小,篩選出各組參數(shù)中與SPAD值變換前后相關(guān)性最大且達(dá)到顯著水平的備選參數(shù),分別為λr、Rg、SDb和NDVI。備選光譜參數(shù)λr、Rg、SDb、NDVI與SPAD變換前后值相關(guān)系數(shù)變化范圍分別為:-0.89~0.79、-0.91~0.90、-0.84~0.84、-0.88~0.89,由Pearson相關(guān)矩陣圖(圖3)可以看出,變換前后的SPAD值變換后與備選光譜參數(shù)相關(guān)性系數(shù)均大于同組其他參數(shù)。因此,本研究中將變換前后的APAD 值作為因變量、備選光譜參數(shù)作為自變量,做進(jìn)一步診斷模型建模分析。

    圖3 玉米冠層SPAD值變換前后與光譜參數(shù)Pearson相關(guān)性分析矩陣圖

    將變換前后的SPAD值與光譜參數(shù)回歸分析結(jié)果顯示(表2),備選光譜參數(shù)與變換前后的SPAD值擬合程度均達(dá)到了顯著性水平,視為該分析結(jié)果具有參考價值。Rg與變換前后的SPAD值回歸方程的R2最大,變化范圍為0.79~0.83,光譜參數(shù)SDb與NDVI的回歸方程次之,變化范圍分別為0.69~0.71、0.74~0.80,λr的R2最小,變化范圍為0.15~0.19。分析結(jié)果還顯示,與SPAD值相比,數(shù)學(xué)變換后的SPAD值與光譜參數(shù)回歸分析的R2值呈下降趨勢,MAPE值呈顯著上升趨勢,RMSR值無明顯變化。由此可見,通過一元回歸分析所建立的診斷模型預(yù)測能力隨參數(shù)變化而呈現(xiàn)顯著差異,最高可解釋因變量83%,最低不足20%,對SPAD進(jìn)行數(shù)學(xué)變換并不能提升診斷模型預(yù)測能力,反而增加了預(yù)測結(jié)果誤差。因此,未進(jìn)行數(shù)學(xué)變換的SPAD值與備選光譜參數(shù)所建立的一元回歸方程可認(rèn)為是本研究中的精度較高的一類預(yù)測模型,本文將該預(yù)測模型統(tǒng)稱為“一元診斷模型”。

    表2 玉米冠層SPAD值變換前后與光譜參數(shù)單變量線性回歸分析

    為了進(jìn)一步提升診斷模型預(yù)測能力,基于上述分析結(jié)果,本文將SPAD值作為因變量、4 個備選光譜參數(shù)作為自變量,采用逐步回歸方法做多變量回歸分析,以回歸分析方程的R2值最接近1 為原則,最終得出基于多變量的SPAD值最優(yōu)預(yù)測模型,結(jié)果如表3 所示,診斷模型的R2值為0.93,MAPE、RMSR值分別為0.13、0.26??梢?,與“一元診斷模型”相比,該模型的預(yù)測精度得到顯著提升,預(yù)測結(jié)果誤差也未出現(xiàn)明顯增加。因此,本文將基于光譜參數(shù)λr、Rg、SDb、NDVI的多元回歸方程作為最優(yōu)預(yù)測模型,并將該預(yù)測模型定義為“多元診斷模型”。

    2.3 玉米葉片SPAD值高光譜影像反演及精度

    為了檢驗診斷模型預(yù)測能力,本文利用數(shù)據(jù)集Ⅱ中的21個樣本數(shù)據(jù)驗證“一元診斷模型”和“多元診斷模型”的精度?;凇耙辉\斷模型”SPAD預(yù)測值與實測值的回歸分析結(jié)果如圖4 所示,擬合方程R2分別為0.6183、0.7345、0.2867、0.8945。盡管基于光譜參數(shù)NDVI的“一元診斷模型”的預(yù)測值與實測值達(dá)到較高的決定系數(shù),但存在偏離1:1 等值線問題,從此說明“一元診斷模型”存在預(yù)測結(jié)果整體偏移現(xiàn)象,其預(yù)測極值能力依舊較弱。同理,圖5 為基于“多元診斷模型”預(yù)測結(jié)果與SPAD實測值的回歸分析結(jié)果,由此可以看出,該模型所預(yù)測的結(jié)果與實測值決定系數(shù)較高,R2為0.9285,擬合直線接近于1:1 等值線,由此證明該診斷模型對極值的預(yù)測能力較強。因此,基于光譜參數(shù)λr、Rg、SDb、NDVI的“多元診斷模型”為本研究所獲取的最優(yōu)SPAD值預(yù)測模型,估測結(jié)果準(zhǔn)確性較高。

    圖4 基于“一元診斷模型”SPAD預(yù)測值與實測值的關(guān)系

    圖5 基于“多元診斷模型”SPAD預(yù)測值與實測值的關(guān)系

    表3“多元診斷模型”精度分析

    綜合上述分析結(jié)果,利用ENVI+IDL 軟件平臺,使用“多元診斷模型”(表3)對試驗區(qū)高光譜影像進(jìn)行填圖(圖6)。結(jié)果顯示,不同施肥處理SPAD值差異十分顯著,能清晰反映不同地塊SPAD值差異,影像上高值區(qū)呈條帶狀分布。這是由于玉米正值拔節(jié)期,部分地塊玉米葉片不能完全覆蓋地面,從而導(dǎo)致玉米行間存在裸露土壤,診斷模型僅對玉米冠層葉片SPAD值估測有效,玉米間空隙區(qū)域因而SPAD估測值較低,故而玉米行間顏色多為紅、黃色圖斑,表示該區(qū)域SPAD較低。在高量施氮處理(F7)地塊,由于氮肥充足,玉米枝葉比較繁茂,完全覆蓋行間土壤,因而影像中條帶狀特征圖斑消失。玉米試驗田臨近裸地區(qū)域多為紅、黃色圖斑,僅存在少量藍(lán)色圖斑,這是由于裸地區(qū)域存在少量雜草,其葉片中的葉綠素引起預(yù)測模型響應(yīng),該結(jié)果也間接證明了該多元診斷模型對植物葉片葉綠素具有一定的敏感性。

    圖6 基于多元診斷模型玉米冠層葉片SPAD值反演圖

    3 結(jié)論

    本研究以中國科學(xué)院封丘農(nóng)業(yè)生態(tài)實驗站長期定位試驗的夏玉米為研究靶區(qū),利用42 個玉米冠層SPAD值為訓(xùn)練樣本,采用回歸分析方法,對高光譜反射率數(shù)據(jù)進(jìn)行反演分析。得出以下結(jié)論。

    (1)對選定的17 個光譜參數(shù)涵蓋了400~900 nm波段范圍內(nèi)光譜反射率曲線特征、曲線面積及植被指數(shù)等信息,采用相關(guān)性分析方法,從同一類型光譜參數(shù)中遴選與數(shù)學(xué)變換前后SPAD值關(guān)系最緊密的參數(shù)作為建模參數(shù)。光譜參數(shù)λr、Rg、SDb、NDVI分別與SPAD原始值具有較高的決定系數(shù),分別為0.19、0.83、0.71、0.80,所構(gòu)建的“多元診斷模型”在預(yù)測精度、結(jié)果穩(wěn)定性等方面均得到明顯提升,在高光譜或多光譜遙感監(jiān)測方面具有一定的應(yīng)用潛力。

    (2)本文采用的參數(shù)篩選與回歸分析建模方法獲取了精度較高的玉米冠層葉綠素含量高光譜預(yù)測模型,并通過獨立樣本對預(yù)測模型精度進(jìn)行驗證分析,其預(yù)測精度與可靠性較前人研究具有一定的提升。本研究論證了回歸分析建模方法在高光譜遙感在作物冠層葉綠素含量估測分析的可行性,“多元診斷模型”的決定系數(shù)高達(dá)0.93,同時揭示了光譜參數(shù)λr、Rg、SDb、NDVI能夠較全面的反演作物冠層葉綠素空間分布特征。

    4 討論

    玉米拔節(jié)期是整個生育期的關(guān)鍵時期,該時期植株莖葉生長快速快,光合作用強度大,葉綠素作為重要光合色素,其含量高低可直接決定植株光合效率并進(jìn)一步影響光合產(chǎn)量[16]。過去學(xué)者研究認(rèn)為,通過監(jiān)測植株葉綠素含量變化能診斷出該時期作物生長健康狀況,已經(jīng)成為制定農(nóng)田管理措施的重要參考依據(jù)[17]。研究還表明,葉綠素含量高低對可見光波段的光譜反射率具有重要影響,其顯著差異主要表現(xiàn)在500~600 nm波段上,葉綠素能吸收大部分的紅光和紫光但反射綠光[18]。基于葉綠素對特定波長光譜吸收和反射等特性的認(rèn)識,國內(nèi)外已經(jīng)開展了大量關(guān)于作物葉綠素光譜學(xué)診斷研究,如Coste 等[19]利用650 nm 和940 nm 處的透射光的相對比值來計算葉綠素含量指標(biāo)。湯旭光等[20]基于可見、近紅外反射光譜反射率估算了作物葉片葉綠素含量,并建立了反演模型。本研究借助長期定位梯度施肥試驗平臺,構(gòu)建出夏玉米冠層葉片SPAD值梯度變化特征,玉米拔節(jié)期冠層光譜反射率變化表現(xiàn)為:在可見光區(qū)域呈現(xiàn)先增加后降低,最后快速上升,在400~725 nm 波段范圍內(nèi)SPAD值越高對光譜吸收能力越強,呈低反射率特征,相反,725~900 nm近紅外波段SPAD值越高光譜反射率越高。這些光譜特征與植物葉片葉綠素的特征密切相關(guān),葉綠素是一種能吸收光能含脂的色素物質(zhì),將吸收的紅光和藍(lán)紫光能量轉(zhuǎn)換成化學(xué)能,而不能吸收綠光光能,因而SPAD值越高光譜反射率越低,在550 nm 左右波段呈現(xiàn)反射峰[21]。由于色素和纖維素在近紅外波段的吸收較低,光線在葉片細(xì)胞結(jié)構(gòu)經(jīng)過多次散射致使接近50%的光線被反射,當(dāng)植株營養(yǎng)狀況較好時葉片內(nèi)葉綠素含量高、細(xì)胞層較多、葉肉與細(xì)胞間空隙較厚能進(jìn)一步增加光譜反射率,最終表現(xiàn)為SPAD值越高,725~900 nm波段光譜反射率越高[22]。根據(jù)上述光譜特征能清晰區(qū)分不同SPAD值對應(yīng)光譜特征曲線(圖2),為進(jìn)一步構(gòu)建玉米冠層葉片SPAD值診斷模型提供了可靠的數(shù)據(jù)基礎(chǔ)。

    以光譜特征曲線為依據(jù),本研究將“三邊”特征參數(shù)與植被指數(shù)作為構(gòu)建診斷模型的候選參數(shù),采用微分光譜的方法,對原始光譜做一階微分,用來表征“三邊”極值、位置、面積等特征信息。微分光譜具有與SPAD值相關(guān)性高、獨立于土壤背景及減弱或消除測量誤差干擾等優(yōu)點,如茹克亞·薩吾提等[23]對原始光譜的數(shù)據(jù)進(jìn)行微分預(yù)處理后,通過0.01水平顯著性檢驗的波段數(shù)量明顯增加,并且指數(shù)微分處理能夠增加葉綠素估算精度。袁杰等[24]研究也認(rèn)為一階微分光譜數(shù)值可以較好的反演作物生長信息和葉綠素密度等信息。本研究以長期定位試驗田為研究靶區(qū),為簡化模型參數(shù),有效消除同類型參數(shù)同時參數(shù)模型構(gòu)建所帶來的數(shù)據(jù)冗余,將獲取的微分光譜數(shù)值和植被指數(shù)等17個光譜參數(shù)根據(jù)其特性劃分為四類。采用相關(guān)性分析方法,根據(jù)光譜參數(shù)與葉綠素值相關(guān)系數(shù)大小遴選的4個構(gòu)建診斷模型的光譜參數(shù),首先將單個光譜參數(shù)與數(shù)學(xué)變換前后的SPAD值構(gòu)建一元回歸預(yù)測方程,與已有研究結(jié)果不同[23],本研究發(fā)現(xiàn)數(shù)學(xué)變換后的SPAD值與光譜參數(shù)間的擬合度并未得到明顯提高。因此,將SPAD原始數(shù)值與λr、Rg、SDb、NDVI的回歸方程作為優(yōu)選的“一元診斷模型”,并進(jìn)一步構(gòu)建“多元診斷模型”。通過對比決定系數(shù)R2大小,“一元診斷模型”預(yù)測能力低于“多元診斷模型”,驗證樣本對這兩種模型驗證結(jié)果更進(jìn)一步印證了“多元診斷模型”的可靠性,其預(yù)測結(jié)果與實測結(jié)果決定系數(shù)R2達(dá)到0.9285,擬合直線斜率為1.0132。與前人研究結(jié)果相比,本研究中構(gòu)建的“多元診斷模型”預(yù)測精度得到明顯提升,因此,本研究推薦基于光譜參數(shù)λr、Rg、SDb、NDVI的“多元診斷模型”為SPAD值最優(yōu)高光譜診斷模型。

    本研究利用高光譜400~900 nm 波段診斷玉米冠層SPAD值變化,其結(jié)果與前人研究具有一致性,但該診斷模型決定系數(shù)有一定提高,光譜參數(shù)數(shù)量更少,模型結(jié)構(gòu)更加簡單[12-13]?!岸嘣\斷模型”包含的光譜參數(shù)綜合考慮了光譜反射曲線特征、曲線面積及植被指數(shù)等信息,能夠完整反映整個光譜波段范圍(400~900 nm)內(nèi)高光譜信息,并有效規(guī)避了因同類型光譜參數(shù)攜帶類似信息所造成的數(shù)據(jù)冗余,提高模型診斷能力的同時簡化了模型參數(shù),在高光譜或多光譜遙感監(jiān)測中具有一定的應(yīng)用潛力。然而,由于驗證樣品數(shù)量有限,本研究所建立的SPAD高光譜診斷模型可靠性及普適性仍需進(jìn)一步驗證。

    猜你喜歡
    診斷模型冠層反射率
    影響Mini LED板油墨層反射率的因素
    近岸水體異源遙感反射率產(chǎn)品的融合方法研究
    海洋通報(2022年4期)2022-10-10 07:40:26
    基于低空遙感的果樹冠層信息提取方法研究
    具有顏色恒常性的光譜反射率重建
    基于激光雷達(dá)的樹形靶標(biāo)冠層葉面積探測模型研究
    安徽省淮南森林冠層輻射傳輸過程的特征
    CD4細(xì)胞計數(shù)聯(lián)合IGRA預(yù)測AIDS/Ⅲ型TB影像診斷模型分析
    甘肅科技(2020年20期)2020-04-13 00:30:56
    一種電網(wǎng)通信網(wǎng)絡(luò)的故障檢測方法及系統(tǒng)
    化學(xué)腐蝕硅表面結(jié)構(gòu)反射率影響因素的研究*
    電子器件(2017年2期)2017-04-25 08:58:37
    施氮水平對冬小麥冠層氨揮發(fā)的影響
    免费观看性生交大片5| 亚洲三级黄色毛片| 男女无遮挡免费网站观看| 精品视频人人做人人爽| 国产色婷婷99| 日韩视频在线欧美| a级片在线免费高清观看视频| 国产欧美亚洲国产| 欧美+日韩+精品| 国产成人一区二区在线| 精品少妇黑人巨大在线播放| 亚洲,欧美,日韩| 欧美最新免费一区二区三区| 人人妻人人澡人人爽人人夜夜| 国产男女超爽视频在线观看| 99视频精品全部免费 在线| 大话2 男鬼变身卡| 亚洲av欧美aⅴ国产| 男女免费视频国产| 亚洲成人av在线免费| 永久免费av网站大全| 国产精品无大码| 精品久久蜜臀av无| 搡老乐熟女国产| 亚洲,一卡二卡三卡| 久久青草综合色| 咕卡用的链子| 下体分泌物呈黄色| 国产精品久久久久久久久免| 新久久久久国产一级毛片| 久久久精品免费免费高清| 亚洲精品美女久久av网站| 亚洲,欧美精品.| 伦精品一区二区三区| 久久久久久久久久久免费av| 熟女电影av网| 午夜日本视频在线| 久久女婷五月综合色啪小说| 高清黄色对白视频在线免费看| 91国产中文字幕| www.色视频.com| √禁漫天堂资源中文www| 亚洲色图综合在线观看| 亚洲婷婷狠狠爱综合网| 久久久久精品性色| 国产精品蜜桃在线观看| 欧美精品高潮呻吟av久久| 日韩在线高清观看一区二区三区| 九草在线视频观看| 人体艺术视频欧美日本| 视频在线观看一区二区三区| 五月伊人婷婷丁香| 中文精品一卡2卡3卡4更新| 久久人人爽av亚洲精品天堂| 我要看黄色一级片免费的| 久久毛片免费看一区二区三区| 精品酒店卫生间| 你懂的网址亚洲精品在线观看| 久久久亚洲精品成人影院| 亚洲欧美成人综合另类久久久| 亚洲,一卡二卡三卡| 热re99久久精品国产66热6| 日韩成人伦理影院| 蜜臀久久99精品久久宅男| 国产精品秋霞免费鲁丝片| 日本猛色少妇xxxxx猛交久久| 国产精品三级大全| 亚洲 欧美一区二区三区| 好男人视频免费观看在线| 亚洲国产精品成人久久小说| 日韩一区二区三区影片| 日日啪夜夜爽| 在线观看免费视频网站a站| 99re6热这里在线精品视频| 99精国产麻豆久久婷婷| 国产永久视频网站| 26uuu在线亚洲综合色| 一级片'在线观看视频| 热re99久久国产66热| 久久精品国产亚洲av天美| 欧美最新免费一区二区三区| 高清av免费在线| 精品一区二区三区视频在线| 免费在线观看完整版高清| 91国产中文字幕| 日韩成人av中文字幕在线观看| 在线观看www视频免费| 高清视频免费观看一区二区| 国产欧美日韩综合在线一区二区| 大片免费播放器 马上看| 韩国精品一区二区三区 | 在线观看免费高清a一片| 九九爱精品视频在线观看| 亚洲人成77777在线视频| 99热网站在线观看| 国产精品久久久久久精品电影小说| 秋霞伦理黄片| av在线观看视频网站免费| 亚洲情色 制服丝袜| 国产精品久久久久久精品古装| 久久久久久久久久成人| 亚洲综合色惰| a 毛片基地| 国产精品人妻久久久久久| 免费高清在线观看视频在线观看| 99久久综合免费| 男女午夜视频在线观看 | 亚洲欧美日韩另类电影网站| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | freevideosex欧美| av卡一久久| 亚洲少妇的诱惑av| 91午夜精品亚洲一区二区三区| 亚洲av电影在线观看一区二区三区| 免费看av在线观看网站| 九九在线视频观看精品| xxx大片免费视频| 欧美精品亚洲一区二区| 国产在线视频一区二区| 人人妻人人澡人人看| 最近手机中文字幕大全| a级片在线免费高清观看视频| 欧美日韩视频精品一区| 国产1区2区3区精品| 亚洲少妇的诱惑av| 国产永久视频网站| 毛片一级片免费看久久久久| 日韩精品免费视频一区二区三区 | 欧美人与性动交α欧美软件 | 看免费av毛片| 亚洲欧美色中文字幕在线| 永久免费av网站大全| 国产精品一区二区在线不卡| 日韩制服丝袜自拍偷拍| 22中文网久久字幕| 国产成人免费无遮挡视频| 丝袜喷水一区| 街头女战士在线观看网站| 啦啦啦视频在线资源免费观看| 男女午夜视频在线观看 | 国产无遮挡羞羞视频在线观看| 亚洲中文av在线| 国产熟女午夜一区二区三区| 久久影院123| 欧美 亚洲 国产 日韩一| 亚洲三级黄色毛片| 五月伊人婷婷丁香| 下体分泌物呈黄色| 久久精品人人爽人人爽视色| 一边亲一边摸免费视频| 香蕉精品网在线| 99精国产麻豆久久婷婷| 亚洲情色 制服丝袜| 内地一区二区视频在线| 日韩人妻精品一区2区三区| 在线观看三级黄色| 精品卡一卡二卡四卡免费| 纯流量卡能插随身wifi吗| 国产av码专区亚洲av| 欧美国产精品一级二级三级| 高清毛片免费看| 国产在线免费精品| 国产成人精品无人区| 日本与韩国留学比较| 成人国产av品久久久| 久久久精品区二区三区| 日本91视频免费播放| 性高湖久久久久久久久免费观看| 超色免费av| 2018国产大陆天天弄谢| av电影中文网址| 久久久久久人人人人人| 蜜桃在线观看..| 成人毛片60女人毛片免费| 亚洲国产av新网站| 亚洲国产精品一区三区| 高清毛片免费看| 边亲边吃奶的免费视频| 天堂8中文在线网| 欧美成人午夜免费资源| 99热6这里只有精品| 日本黄色日本黄色录像| 少妇的逼好多水| 啦啦啦啦在线视频资源| kizo精华| 亚洲精品视频女| 日韩中文字幕视频在线看片| 99久久综合免费| 久热久热在线精品观看| 中国美白少妇内射xxxbb| 交换朋友夫妻互换小说| 日韩中文字幕视频在线看片| 久久99热6这里只有精品| 国产精品一区二区在线观看99| 26uuu在线亚洲综合色| 久久久久久久精品精品| 亚洲欧美日韩另类电影网站| 少妇被粗大的猛进出69影院 | 日韩不卡一区二区三区视频在线| 国产综合精华液| 亚洲欧美日韩另类电影网站| 亚洲精品久久成人aⅴ小说| 成人毛片60女人毛片免费| 国产综合精华液| 国产精品麻豆人妻色哟哟久久| av在线观看视频网站免费| 精品卡一卡二卡四卡免费| 乱码一卡2卡4卡精品| 国产高清不卡午夜福利| 日本与韩国留学比较| 国产男女内射视频| 高清欧美精品videossex| 日本91视频免费播放| 亚洲成人av在线免费| 中国三级夫妇交换| 欧美成人精品欧美一级黄| 大话2 男鬼变身卡| 久久午夜福利片| 国产av一区二区精品久久| 亚洲欧洲国产日韩| 啦啦啦啦在线视频资源| 国产永久视频网站| 亚洲四区av| 精品少妇黑人巨大在线播放| 亚洲内射少妇av| 18禁观看日本| 日本与韩国留学比较| 亚洲四区av| 在线天堂中文资源库| 国产亚洲午夜精品一区二区久久| 免费黄频网站在线观看国产| 亚洲美女黄色视频免费看| av不卡在线播放| 精品人妻熟女毛片av久久网站| 大码成人一级视频| 亚洲av在线观看美女高潮| 丝袜人妻中文字幕| 亚洲欧美中文字幕日韩二区| 狠狠婷婷综合久久久久久88av| 伦理电影免费视频| 日韩制服丝袜自拍偷拍| 国产精品偷伦视频观看了| 一本久久精品| 国产在线视频一区二区| 亚洲av电影在线观看一区二区三区| 老女人水多毛片| 欧美日韩视频精品一区| 韩国av在线不卡| 精品国产露脸久久av麻豆| 黄片播放在线免费| 国产精品国产av在线观看| 日本vs欧美在线观看视频| 热re99久久精品国产66热6| 制服丝袜香蕉在线| 99九九在线精品视频| 日韩成人伦理影院| 黑人猛操日本美女一级片| av黄色大香蕉| 亚洲美女视频黄频| 国产亚洲一区二区精品| 亚洲精品视频女| 春色校园在线视频观看| 欧美最新免费一区二区三区| 国产成人免费无遮挡视频| 赤兔流量卡办理| 国产精品国产三级国产av玫瑰| 国产成人欧美| 免费人妻精品一区二区三区视频| 久久久久国产精品人妻一区二区| 亚洲av日韩在线播放| 亚洲色图综合在线观看| 99re6热这里在线精品视频| 免费久久久久久久精品成人欧美视频 | 妹子高潮喷水视频| 亚洲精品视频女| 亚洲av免费高清在线观看| 国产精品.久久久| av天堂久久9| 男人添女人高潮全过程视频| 亚洲成国产人片在线观看| 女的被弄到高潮叫床怎么办| 制服诱惑二区| a级毛色黄片| 国产免费一级a男人的天堂| 少妇猛男粗大的猛烈进出视频| 宅男免费午夜| 欧美国产精品一级二级三级| 国产淫语在线视频| 国产一区亚洲一区在线观看| 国产亚洲av片在线观看秒播厂| 国产熟女午夜一区二区三区| 国产精品久久久久久精品古装| 一级片免费观看大全| 精品久久久精品久久久| 国产男女超爽视频在线观看| 少妇猛男粗大的猛烈进出视频| 亚洲,欧美精品.| 色哟哟·www| 久久久a久久爽久久v久久| 国产一区二区在线观看日韩| 国产成人欧美| 男人操女人黄网站| 国产亚洲一区二区精品| av播播在线观看一区| 午夜激情av网站| 日本-黄色视频高清免费观看| 国产1区2区3区精品| 午夜福利,免费看| 免费av中文字幕在线| 欧美精品高潮呻吟av久久| 精品一区二区免费观看| 久久久久国产精品人妻一区二区| 成人免费观看视频高清| 亚洲成国产人片在线观看| 青春草国产在线视频| 国产av精品麻豆| 成人影院久久| 欧美 日韩 精品 国产| 免费观看av网站的网址| 日韩欧美精品免费久久| 国产一区二区激情短视频 | 在线天堂最新版资源| 2018国产大陆天天弄谢| 日本欧美视频一区| 777米奇影视久久| 国产亚洲精品第一综合不卡 | 男人添女人高潮全过程视频| 国产男人的电影天堂91| 男女无遮挡免费网站观看| 亚洲一级一片aⅴ在线观看| 韩国精品一区二区三区 | 久久久久久久亚洲中文字幕| 各种免费的搞黄视频| 久久精品国产a三级三级三级| av免费观看日本| 久久99热6这里只有精品| 男女啪啪激烈高潮av片| 51国产日韩欧美| 欧美精品人与动牲交sv欧美| 久久影院123| 精品人妻在线不人妻| 少妇人妻久久综合中文| 亚洲欧美一区二区三区国产| 午夜福利,免费看| 美女国产视频在线观看| 亚洲图色成人| 免费播放大片免费观看视频在线观看| 韩国精品一区二区三区 | 国产色爽女视频免费观看| 精品国产一区二区久久| 久久人人97超碰香蕉20202| av福利片在线| 国产爽快片一区二区三区| 欧美日韩精品成人综合77777| 热re99久久精品国产66热6| 母亲3免费完整高清在线观看 | 人妻系列 视频| 久久久亚洲精品成人影院| 亚洲av免费高清在线观看| 全区人妻精品视频| 老司机亚洲免费影院| 午夜福利在线观看免费完整高清在| 午夜91福利影院| 波野结衣二区三区在线| 丝袜美足系列| 免费高清在线观看日韩| 中文字幕精品免费在线观看视频 | 久久精品夜色国产| 欧美亚洲 丝袜 人妻 在线| 男人爽女人下面视频在线观看| 国产一区二区三区综合在线观看 | 熟女av电影| av福利片在线| 免费黄色在线免费观看| av有码第一页| 国产国拍精品亚洲av在线观看| 久久精品国产鲁丝片午夜精品| 黄片无遮挡物在线观看| 夜夜骑夜夜射夜夜干| 深夜精品福利| 极品人妻少妇av视频| 国产亚洲一区二区精品| av.在线天堂| 日本爱情动作片www.在线观看| 国产亚洲精品第一综合不卡 | 99九九在线精品视频| 午夜免费男女啪啪视频观看| www.色视频.com| 亚洲成人一二三区av| 在线看a的网站| 少妇猛男粗大的猛烈进出视频| av视频免费观看在线观看| 哪个播放器可以免费观看大片| 99久久精品国产国产毛片| 少妇的逼水好多| 日本vs欧美在线观看视频| 国产乱来视频区| 国产av精品麻豆| 国精品久久久久久国模美| 久久久欧美国产精品| av又黄又爽大尺度在线免费看| 巨乳人妻的诱惑在线观看| 五月玫瑰六月丁香| 久久久欧美国产精品| 精品第一国产精品| 成年人免费黄色播放视频| 国产成人免费无遮挡视频| 亚洲经典国产精华液单| 99久久精品国产国产毛片| 亚洲欧美一区二区三区黑人 | 亚洲精品国产av蜜桃| 一区二区三区精品91| 亚洲欧美清纯卡通| 国产探花极品一区二区| 丰满乱子伦码专区| 丝瓜视频免费看黄片| 免费女性裸体啪啪无遮挡网站| 国产精品女同一区二区软件| 男女午夜视频在线观看 | 18禁观看日本| 最近最新中文字幕免费大全7| 国产色婷婷99| av天堂久久9| 亚洲久久久国产精品| 一二三四中文在线观看免费高清| 欧美日韩综合久久久久久| 国内精品宾馆在线| 日韩av不卡免费在线播放| 少妇 在线观看| 中文欧美无线码| 中文字幕制服av| 亚洲伊人久久精品综合| 国产又爽黄色视频| 欧美国产精品一级二级三级| 国产探花极品一区二区| 国产免费一区二区三区四区乱码| 午夜福利网站1000一区二区三区| 人妻系列 视频| 亚洲天堂av无毛| 免费黄网站久久成人精品| 一级毛片 在线播放| 国产色婷婷99| 国产精品女同一区二区软件| 久久精品人人爽人人爽视色| www日本在线高清视频| 一级片免费观看大全| 黄网站色视频无遮挡免费观看| 日本午夜av视频| 亚洲图色成人| 亚洲内射少妇av| 日韩成人伦理影院| 99久久综合免费| 欧美xxⅹ黑人| 人妻系列 视频| 中文字幕人妻熟女乱码| 日本91视频免费播放| 色哟哟·www| 人妻一区二区av| 久久97久久精品| 波野结衣二区三区在线| a级毛片黄视频| 欧美精品一区二区免费开放| 韩国av在线不卡| 婷婷色综合大香蕉| 满18在线观看网站| 一级毛片我不卡| 午夜福利网站1000一区二区三区| 建设人人有责人人尽责人人享有的| 欧美人与性动交α欧美精品济南到 | 久久青草综合色| 久久久久久久久久人人人人人人| 久久久久精品性色| 一区二区三区精品91| 欧美 日韩 精品 国产| 日韩精品有码人妻一区| 国产精品久久久久久久久免| 亚洲美女黄色视频免费看| 久久精品久久久久久久性| 久久精品国产自在天天线| 久久精品国产亚洲av天美| 777米奇影视久久| 亚洲,欧美精品.| 国产免费视频播放在线视频| 国产成人精品无人区| 丝袜美足系列| 狠狠婷婷综合久久久久久88av| 亚洲国产精品成人久久小说| 如日韩欧美国产精品一区二区三区| 亚洲精品一区蜜桃| 午夜久久久在线观看| 国产精品国产av在线观看| 久久97久久精品| 中文字幕另类日韩欧美亚洲嫩草| 99热全是精品| 街头女战士在线观看网站| 99久久人妻综合| 啦啦啦在线观看免费高清www| 熟女av电影| 热re99久久国产66热| 成人毛片60女人毛片免费| 在线观看www视频免费| 精品一区二区免费观看| 午夜免费男女啪啪视频观看| 男男h啪啪无遮挡| 黑丝袜美女国产一区| 日韩 亚洲 欧美在线| 国产 一区精品| 一个人免费看片子| 国产日韩一区二区三区精品不卡| 亚洲经典国产精华液单| 黄色视频在线播放观看不卡| 最黄视频免费看| 亚洲国产毛片av蜜桃av| 久久影院123| 国产精品一区二区在线观看99| 国产精品一区www在线观看| av免费观看日本| 日本欧美视频一区| 亚洲 欧美一区二区三区| 51午夜福利影视在线观看| 久久久久久免费高清国产稀缺| 两人在一起打扑克的视频| 国产免费av片在线观看野外av| 91精品国产国语对白视频| 最近最新中文字幕大全免费视频| 久久久久久免费高清国产稀缺| 夜夜夜夜夜久久久久| 国产不卡av网站在线观看| 亚洲中文日韩欧美视频| 亚洲精华国产精华精| 丝瓜视频免费看黄片| 欧美中文综合在线视频| 国产精品.久久久| 欧美不卡视频在线免费观看 | 国产1区2区3区精品| 男女下面插进去视频免费观看| 国产一区二区三区视频了| 国产单亲对白刺激| 99re6热这里在线精品视频| 又黄又粗又硬又大视频| 成人手机av| 国产一区有黄有色的免费视频| 夜夜躁狠狠躁天天躁| 性色av乱码一区二区三区2| 超碰97精品在线观看| 中出人妻视频一区二区| 久久久久精品人妻al黑| 欧美激情久久久久久爽电影 | 又大又爽又粗| 国产日韩欧美亚洲二区| 久久热在线av| 51午夜福利影视在线观看| 老汉色∧v一级毛片| 老司机午夜福利在线观看视频| 深夜精品福利| 色婷婷久久久亚洲欧美| 99精品久久久久人妻精品| 国产精品香港三级国产av潘金莲| 国产1区2区3区精品| 日韩人妻精品一区2区三区| 日韩欧美国产一区二区入口| 国产人伦9x9x在线观看| 国产男女超爽视频在线观看| 亚洲三区欧美一区| 国产精品亚洲av一区麻豆| 一二三四社区在线视频社区8| 欧美精品啪啪一区二区三区| 99热国产这里只有精品6| 精品视频人人做人人爽| 窝窝影院91人妻| 中亚洲国语对白在线视频| 婷婷精品国产亚洲av在线 | 五月开心婷婷网| 亚洲精品一卡2卡三卡4卡5卡| 日本vs欧美在线观看视频| 精品乱码久久久久久99久播| 国产精品乱码一区二三区的特点 | 亚洲性夜色夜夜综合| 亚洲va日本ⅴa欧美va伊人久久| 一级片免费观看大全| 国产深夜福利视频在线观看| 国产日韩一区二区三区精品不卡| 少妇裸体淫交视频免费看高清 | 一级a爱片免费观看的视频| 亚洲精品久久午夜乱码| 久久亚洲真实| 日本五十路高清| 人人澡人人妻人| 男女下面插进去视频免费观看| 丝瓜视频免费看黄片| 久久午夜综合久久蜜桃| 一夜夜www| 精品电影一区二区在线| 亚洲欧洲精品一区二区精品久久久| 亚洲九九香蕉| 身体一侧抽搐| 法律面前人人平等表现在哪些方面| 成在线人永久免费视频| 亚洲精品国产区一区二| 亚洲精品国产一区二区精华液| 免费高清在线观看日韩| 在线观看www视频免费| 国产男女内射视频| 久久国产精品大桥未久av| 国产99白浆流出| 国产精品久久电影中文字幕 | av国产精品久久久久影院| 大型av网站在线播放| 午夜免费成人在线视频| xxx96com| 免费人成视频x8x8入口观看| 精品亚洲成a人片在线观看| 脱女人内裤的视频| 中国美女看黄片| 欧美在线黄色|