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

    尺度效應(yīng)對多光譜遙感監(jiān)測土壤水分影響研究

    2023-03-04 06:28:08羅亞東郭宇宏孫智鵬金煜龍陳俊英余衛(wèi)華
    節(jié)水灌溉 2023年2期
    關(guān)鍵詞:特征模型

    羅亞東,許 齊,郭宇宏,2,孫智鵬,金煜龍,陳俊英,2,余衛(wèi)華,2

    (1.西北農(nóng)林科技大學(xué)水利與建筑工程學(xué)院,陜西 楊凌 712100;2.西北農(nóng)林科技大學(xué) 旱區(qū)農(nóng)業(yè)水土工程教育部重點實驗室,陜西 楊凌 712100)

    0 引言

    土壤含水率是分析冬小麥生育情況的主要依據(jù),對于理解和解決農(nóng)業(yè)生產(chǎn)、生態(tài)規(guī)劃以及水資源管理中的科學(xué)與實際問題至關(guān)重要[1]。傳統(tǒng)的土壤含水率監(jiān)測方法操作復(fù)雜、且數(shù)據(jù)滯后。遙感技術(shù)因其快速、簡便以及無損的特點,已廣泛應(yīng)用于農(nóng)田農(nóng)情的監(jiān)測[2-7]?;谶b感的墑情監(jiān)測是農(nóng)村水利科技工作的重點領(lǐng)域[8]。但是在無人機多光譜遙感實踐中,存在尺度效應(yīng)影響監(jiān)測精度的問題。

    近年來,遙感監(jiān)測已廣泛應(yīng)用于土壤含水率的研究,楊珺博[9]等人利用無人機遙感監(jiān)測數(shù)據(jù)反演土壤含水率,發(fā)現(xiàn)反演精度可達80%以上,且反演效果最好的土壤深度為10~20 cm。Cheng Minghan[10]等利用無人機多模態(tài)數(shù)據(jù)對玉米田間含水率進行4種模型建立與反演,得出多光譜數(shù)據(jù)在不同模型下均提供最精確的土壤水分估算。

    尺度效應(yīng)是指在無人機遙感觀測中,隨著遙感分辨率的變化,遙感反演得到的數(shù)據(jù)與實際表現(xiàn)不一致的現(xiàn)象,作為遙感科學(xué)的核心問題,近年來受到了學(xué)者的廣泛關(guān)注[11,12]。歐盛華[13]等利用機載航空WIDAS 數(shù)據(jù),發(fā)現(xiàn)尺度效應(yīng)是空間異質(zhì)性和反演函數(shù)非線性的綜合體現(xiàn),空間異質(zhì)性和反演函數(shù)非線性程度的增加均會造成尺度效應(yīng)的增大。徐凱健[14]通過對多光譜影像的森林樹種識別及其空間尺度響應(yīng)的研究。結(jié)果表明,影像空間分辨率對區(qū)域樹種識別結(jié)果具有顯著影響。Zhang[15]等利用偏振遙感器,根據(jù)地物目標(biāo)的空間異質(zhì)性和成像分辨率帶來的限制,發(fā)現(xiàn)了尺度效應(yīng)對地物偏振圖像的影響。Tang[16]等對SEBAL 模型的空間尺度效應(yīng)進行探究,得出高低分辨率對觀測熱通量的影響??梢?,尺度效應(yīng)對于無人機遙感觀測有一定影響,但目前尺度效應(yīng)對于無人機遙感監(jiān)測土壤水分的研究還鮮有報道,在實踐中缺乏最佳的監(jiān)測尺度。在提取作物冠層信息過程中,土壤、陰影背景均會對冠層信息準(zhǔn)確性造成影響[17,18]。楊帥[19]等基于RGRI指數(shù)法獲取的玉米冠層溫度與實測冠層溫度的相關(guān)性,得出不剔除土壤背景效果最差。王佳兒[20]等通過剔除土壤背景反演玉米田間含水率,發(fā)現(xiàn)數(shù)據(jù)的相關(guān)性和回歸模型的精度始終比未剔除土壤背景數(shù)據(jù)的相關(guān)性和模型精度差。紋理特征在冠層信息提取中中是一種相當(dāng)重要的特征,它反映了圖像或物體本身的屬性[21]。王思宇[22]利用紋理特征為基礎(chǔ),基于多光譜影像并構(gòu)建玉米冠層FAPAR 估算模型,發(fā)現(xiàn)其精度較高。目前仍缺乏通過去土壤去陰影掩膜方法提取作物冠層信息。

    本文通過無人機搭載六波段多光譜相機在不同高度下獲取冬小麥紋理特征,剔除土壤背景、陰影背景,在此基礎(chǔ)上研究土壤含水率和植被紋理特征之間的關(guān)系。探究在不同分辨率下無人機遙感數(shù)據(jù)與土壤含水率的相關(guān)性并分析尺度效應(yīng)對于監(jiān)測數(shù)據(jù)的影響。最后得出尺度效應(yīng)對無人機多光譜遙感監(jiān)測土壤含水率的影響,獲取無人機監(jiān)測土壤水分的最佳飛行高度。

    1 材料與方法

    1.1 研究區(qū)域概況

    研究區(qū)域位于陜西省楊陵區(qū)陜西楊凌西北農(nóng)林科技大學(xué)旱區(qū)農(nóng)業(yè)水土工程教育部重點實驗室的灌溉試驗站(34°17'N,108°04'E,海拔521 m)。該地區(qū)為典型的溫帶大陸性季風(fēng)季候,夏季炎熱,冬季寒冷,蒸發(fā)強度較大。年平均氣溫在12 ℃左右,多年平均降水量為630 mm,年平均蒸發(fā)量為884 mm,土壤干容重為1.38 g/cm3。

    1.2 試驗材料與設(shè)計

    本試驗冬小麥選用小偃22 號,于2020年10月20日進行人工播種,共播種7.68 kg種子,并且施用20 kg的復(fù)合肥料作為底肥。同時在其生長過程中嚴(yán)格控制水分條件,試驗田灌水上限分別采用田間持水量的50%、65%、80%和95%,設(shè)置4 個水分處理,3 個重復(fù),共12 個小區(qū),以獲得不同田持下的土壤含水率數(shù)據(jù),使數(shù)據(jù)更有代表性。每個小區(qū)的面積為4 m×4 m=16 m2,每個小區(qū)間隔150 cm,相鄰兩個小區(qū)之間設(shè)置保護行。小區(qū)內(nèi)種植的冬小麥行距為25 cm,每一行播種40 g種子。為減少外界環(huán)境的干擾,每個小區(qū)都配有擋雨棚,圖1為試驗小區(qū)布置。

    圖1 實驗小區(qū)布置Fig.1 Layout of experimental plot

    1.3 數(shù)據(jù)采集

    1.3.1 土壤水分獲取

    分別在3 個生育期在在每個實驗小區(qū)取3 個采樣點A、B、C,在采樣點處用土鉆采集10、20、30、40、60 cm 深度的土樣,每個土樣30 g 左右,土樣采出后立即裝入鋁盒中。將采集完成的鋁盒放入干燥箱,在105 ℃干燥后稱重并計算土壤含水率。烘干法測定土壤質(zhì)量含水率的計算式為:

    式中:W1為鋁盒加濕土壤質(zhì)量;W2為鋁盒加干土質(zhì)量;W3為空鋁盒質(zhì)量。

    1.3.2 無人機平臺與冠層數(shù)據(jù)獲取

    在試驗田內(nèi)選取晴朗天氣下的土壤樣本,以配合無人機多光譜遙感圖像的獲取。采集時間為冬小麥分蘗期(2020年12月24日)、拔節(jié)期(2021年3月23日)、抽穗期(2021年4月19日),每日11∶00-14∶00。

    本試驗采用由深圳市大疆創(chuàng)新科技有限公司生產(chǎn)的經(jīng)緯Matrice 600 六旋翼無人機,搭載Micro-MCA 多光譜相機進行多光譜遙感圖像采集。Micro-MCA 相機共有6個鏡頭,分別對應(yīng)6 個波段,波長分別為900 nm(近紅外)、800 nm(近紅外)、720 nm(紅邊)、680 nm(紅光)、550 nm(綠光)、490 nm(藍光)。無人機飛行高度設(shè)置為19、37、55、74 和92 m,對應(yīng)的分辨率分別為10、20、30、40 和50 mm,鏡頭垂直向下采集光譜信息。拍攝之前在實驗區(qū)附近放置白板和黑板,用于研究區(qū)內(nèi)進行標(biāo)定,從而獲取冬小麥光譜反射率。

    1.4 多光譜圖像處理

    使用多光譜相機完成圖像采集后,使用PixelWrench2[23]軟件對分別對分蘗期、拔節(jié)期、抽穗期無人機多光譜遙感原始圖像進行預(yù)處理,將預(yù)處理后遙感圖像導(dǎo)入Pix4Dmapper 軟件,進行幾何校正、輻射定標(biāo)、試驗區(qū)拼接處理,得到試驗區(qū)多光譜遙感拼接圖像。

    將拼接處理后的圖像導(dǎo)入ENVI5.3(64 bit)軟件中,剔除土壤和陰影背景。使用ENVI5.3(64 bit)軟件中Co-occurrence measure 工具,獲取多光譜圖像6 個波段基于二階概率統(tǒng)計濾波的均值、方差、協(xié)同性、對比度、相異性、信息熵、二階矩、相關(guān)性等八種紋理特征。即從每個多光譜圖像中提取48組紋理特征。

    1.5 灰色關(guān)聯(lián)法優(yōu)選紋理特征

    不同冠層光譜反射率和紋理特征之間差異較大,存在較強的共線性問題[23]。從理論上來[24]說,預(yù)測指標(biāo)與參考變量之間的相關(guān)性越強,其預(yù)測結(jié)果越可靠[25]。采用全變量建模時,會導(dǎo)致模型預(yù)測不準(zhǔn)確。

    基于灰色系統(tǒng)理論的灰色關(guān)聯(lián)分析是一種利用比較序列和參考序列之間的相似度來判別各變量相關(guān)性強弱的方法,其模型準(zhǔn)確度高、數(shù)據(jù)流失性小[26,27]。

    本文使用灰色關(guān)聯(lián)法研究步驟如下:

    構(gòu)建參考序列:

    構(gòu)建比較序列:

    計算關(guān)聯(lián)系數(shù):

    計算關(guān)聯(lián)度:

    式中:y表示土壤含水率;xa、xb表示一個生育期下的紋理特征;ξi代表關(guān)聯(lián)系數(shù);si代表關(guān)聯(lián)度;n表示3 個生育期下12塊小區(qū)各自對應(yīng)的土壤含水率,n=36;k為6 個波段下的的紋理特征總數(shù),k=48;n為樣本組數(shù),n=36;i表示3 組掩膜處理,i=3。

    1.6 反演模型構(gòu)建

    本研究每組樣本經(jīng)灰色關(guān)聯(lián)法優(yōu)選后為8個輸入變量,采用多元線性回歸(MLR)、前饋神經(jīng)網(wǎng)絡(luò)(BPNN)和隨機森林(RF)3種回歸方法在5 個高度和3 種掩膜方法下分別建立土壤水分反演模型,共45個土壤水分模型。將36組樣本按2∶1的比例隨機取樣作為建模集和驗證集。即選取24 個樣本用于建模,其余12個用于驗證。

    MLR 利用線性模型來擬合多個自變量和因變量的關(guān)系,從而確定多元線性回歸模型的參數(shù),回歸至原假設(shè)方程中,通過回歸方程來預(yù)測因變量的趨勢[28-30]。

    BPNN 作為一種根據(jù)誤差逆?zhèn)鞑ニ惴ㄓ?xùn)練的多次前饋網(wǎng)絡(luò),它通過對人類神經(jīng)元的功能進行模擬,可以儲存及學(xué)習(xí)大量的輸入數(shù)據(jù)和輸出數(shù)據(jù),且不需對變量的映射關(guān)系進行描述,利用輸入和輸出數(shù)據(jù)建模,其對非線性系統(tǒng)具有很強的模擬能力[31,32]。本文搭建4 層BP 神經(jīng)網(wǎng)絡(luò),采用兩個隱含層,傳遞函數(shù)分別為雙曲正切S型函數(shù)和線性函數(shù),節(jié)點數(shù)分別為8和2。

    RF 是一種基于分類樹的算法,通過對大量分類樹的匯總提高了模型的預(yù)測精度,是取代神經(jīng)網(wǎng)絡(luò)等傳統(tǒng)機器學(xué)習(xí)方法的新的模型。當(dāng)林中的樹數(shù)變大時,林的泛化誤差會收斂到極限??梢泽w現(xiàn)變量間的交互作用,具有很高的預(yù)測準(zhǔn)確率,對異常值和噪聲具有很好的容忍度,且不容易出現(xiàn)過擬合[33,34]。

    構(gòu)建反演模型,比較3種掩膜方法下各反演模型的擬合優(yōu)度R2和均方根誤差RMSE,分析尺度效應(yīng)對于無人機遙感監(jiān)測土壤水分的影響。R2越接近1,RMSE越接近0,則反演模型的精度和準(zhǔn)確性越好[23]。

    2 結(jié)果與分析

    2.1 土壤含水率

    通過取土烘干法獲取的不同生育期下每個小區(qū)土壤含水率梯度值,結(jié)果如圖2所示。

    從圖2可以看出,不同生育期的土壤含水率有所不同,各小區(qū)分別在分蘗期、拔節(jié)期和抽穗期的土壤含水率依次降低,土壤含水率的覆蓋范圍較大且具有代表性。

    圖2 各生育期各小區(qū)土壤含水率Fig.2 SWC of each plot at each fertility stage

    2.2 灰色關(guān)聯(lián)法處理結(jié)果

    獲得3組不同掩膜方法下的紋理特征,每組數(shù)據(jù)統(tǒng)計各自48 個紋理特征。本文通過灰色關(guān)聯(lián)法結(jié)果顯示紋理特征與土壤含水率的相關(guān)度,進而對各分辨率與掩膜方法進行初步評價與分析。

    使用灰色關(guān)聯(lián)法計算紋理特征與土壤含水率相關(guān)度,結(jié)果如圖3所示。

    圖3 灰色關(guān)聯(lián)法處理結(jié)果Fig.3 Results of the GRA method

    無掩膜(NM)處理中,10 mm 和30 mm 下紋理特征與土壤含水率相關(guān)性最顯著,20 mm和40 mm下紋理特征與土壤含水率相關(guān)性總體最佳,關(guān)聯(lián)度分布在0.7~0.8。去土壤掩膜(SRM)處理中,紋理特征與土壤含水率的相關(guān)性變化不大。去土壤去陰影掩膜(SSRM)處理中,總體相關(guān)性較差,分辨率為10 mm時紋理特征與土壤含水率相關(guān)性最顯著,分辨率為40 mm 時次之。

    從圖3可以看出,在10 mm 下,紋理特征與土壤含水率呈現(xiàn)最高相關(guān)性,其次則為30 mm、20 mm、40 mm、50 mm。這是由于飛行高度提高,分辨率隨之降低,混合像元在總像元中的比例上升,對冬小麥紋理特征造成影響,紋理特征與土壤含水率的相關(guān)性下降。

    在掩膜分析結(jié)果中,無掩膜與去土壤掩膜的結(jié)果差異較小,且優(yōu)于去土壤去陰影掩膜。本文認(rèn)為在進行去土壤掩膜處理時,混合像元容易分類,而陰影和冠層的差異較??;而在進行去土壤去陰影處理時,混合像元難以分類,導(dǎo)致過多地保留陰影或刪除必要的冠層像元導(dǎo)致紋理特征不能準(zhǔn)確反映冠層情況。

    2.3 三種掩膜方法與不同分辨率相關(guān)性分析

    2.3.1 無掩膜處理下反演結(jié)果

    圖4可以看出,多元線性回歸模型在分辨率為10 mm 時擬合優(yōu)度(R2)最高,其次是40 mm。BP 神經(jīng)網(wǎng)絡(luò)模型在分辨率為10 mm 時R2最高,30 mm 次之。隨機森林模型在分辨率為10 mm時R2較高,其次是30 mm。

    圖4 無掩膜處理反演結(jié)果Fig.4 Regression results without mask

    在無掩膜處理中,3 種回歸方法的最佳反演結(jié)果均在10 mm 時取得較好結(jié)果,擬合優(yōu)度在0.7 以上,其中BP 神經(jīng)網(wǎng)絡(luò)的擬合優(yōu)度為0.86,RMSE為0.011。值得注意的是,BP 神經(jīng)網(wǎng)絡(luò)和隨機森林兩種方法在30 mm 分辨率時R2超過了10 mm,因此在10 mm和30 mm分辨率下,冠層情況在紋理特征中得到準(zhǔn)確表達。

    同時,3組回歸模型中擬合優(yōu)度總體上呈現(xiàn)出先降低再升高再降低的趨勢,在分辨率降低到30 mm 至40 mm 時R2出現(xiàn)升高趨勢。在此掩膜處理下,多元線性回歸模型與BP 神經(jīng)網(wǎng)絡(luò)的擬合優(yōu)度低于隨機森林模型。

    2.3.2 去土壤掩膜處理下反演結(jié)果

    將篩選后去土壤掩膜處理的36 組濾波信息分別使用多元線性回歸、BP 神經(jīng)網(wǎng)絡(luò)、隨機森林進行反演處理,結(jié)果如圖5所示。

    由圖5得出,多元線性回歸模型在分辨率為40 mm時R2最高,其次是30 mm。BP 神經(jīng)網(wǎng)絡(luò)模型在分辨率為30 mm 時R2最高,其次是40 mm。隨機森林模型在分辨率為10 mm,40 mm次之。

    圖5 去土壤掩膜處理反演結(jié)果Fig.5 Regression results with soil background removed by masking

    在去土壤掩膜處理中,3種回歸方法的最佳反演結(jié)果不盡相同,總體分布顯示在該掩膜處理下,40 mm分辨率下冠層情況在紋理特征中得到準(zhǔn)確表達。

    3組回歸模型中擬合優(yōu)度總體上呈現(xiàn)出先降低再升高再降低的趨勢,除10 mm 分辨率的結(jié)果較好之外,30 mm 與40 mm分辨率下的反演結(jié)果擬合優(yōu)度較高。在去土壤掩膜處理之下,3種模型的擬合優(yōu)度均在0.6~0.7之間。

    2.3.3 去土壤去陰影掩膜處理下反演結(jié)果

    將篩選后去土壤去陰影掩膜處理的36 組濾波信息分別使用多元線性回歸、BP 神經(jīng)網(wǎng)絡(luò)、隨機森林進行反演處理,結(jié)果如圖6所示。

    圖6 去土壤去陰影掩膜處理反演結(jié)果Fig.6 Regression results with soil and shadow background removed by masking

    圖6可以看出,在去土壤去陰影掩膜處理中,多元線性回歸模型在分辨率為10 mm 時R2最高,其次是20 mm。BP 神經(jīng)網(wǎng)絡(luò)模型在分辨率為10 mm時R2最高,其次是40 mm。隨機森林模型在分辨率為10 mm,40 mm次之。

    在此種掩膜處理中,3 種回歸方法的擬合結(jié)果顯示,在10 mm與40 mm分辨率下取得相對較高的反演精度,見圖7。

    圖7 各掩膜方法下擬合優(yōu)度對比Fig.7 Comparison of the goodness of fit under different masking method

    3組回歸模型中擬合優(yōu)度總體上呈現(xiàn)出逐漸降低趨勢,值得注意的是,在去土壤去陰影掩膜處理之下,除10 mm 的擬合優(yōu)度高于0.7,其余分辨率在各自模型下的擬合優(yōu)度均在0.5左右。

    2.4 討論與分析

    分析反演結(jié)果可以得出在分辨率為10 mm和40 mm時,模型反演效果較好。本文認(rèn)為飛行高度較低時,遙感圖像分辨率較高,能較為精準(zhǔn)地表現(xiàn)冠層狀況,同樣使得冠層情況也能在紋理特征中得到準(zhǔn)確表達;分辨率為40 mm 時,反演效果僅次于最高分辨率10 mm,本文認(rèn)為該現(xiàn)象由于是混合像元在總像元中占比降低,當(dāng)分辨率由10 mm下降至30 mm時,混合像元在在總像元中的比重增加,在用閾值進行背景分割時,會保留部分混合像元中的背景,或?qū)⒉糠只旌舷裨械淖魑锾蕹?,?dǎo)致紋理特征不能準(zhǔn)確表達冠層信息。當(dāng)分辨率下降至40 mm 時,分辨率與葉面寬度相近,此時混合像元占比反而下降,從而獲得較為準(zhǔn)確的反演精度。由于遙感圖像中存在土壤、陰影和田間雜草等背景,在用紋理特征進行建模時作為無關(guān)變量干擾了模型精度。導(dǎo)致回歸模型的魯棒性較差。當(dāng)高度增加,分辨率降低至與作物葉面尺度相近時,混合像元中背景為優(yōu)勢類別的像元被完全剔除,冠層為優(yōu)勢類別的像元被保留,并且較少的背景像元與冠層像元混合后對紋理特征的影響降低,能一定程度上減小背景對回歸模型的干擾。

    最后,在各掩膜處理下,3種建模方法的反演結(jié)果不盡相同。在無掩膜處理下,多元線性回歸的擬合精度最差;在去土壤掩膜處理下,隨機森林的擬合精度最差;而在去土壤陰影掩膜處理中,BP 神經(jīng)網(wǎng)絡(luò)的擬合精度最差。本文認(rèn)為原因如下:多元線性回歸對數(shù)據(jù)有比較嚴(yán)格的要求如:①隨機誤差項是一個期望值為0 的隨機變量;②解釋變量的所有觀測值,隨機誤差項有相同的方差;③隨機誤差項彼此不相關(guān)等。在無掩膜處理的紋理特征中,干擾值很多隨機誤差不可控,所以這些條件在無掩膜的情況下很難滿足。隨機森林在回歸分析中不能連順序輸出,也不能預(yù)測超出訓(xùn)練集范圍的數(shù)據(jù),導(dǎo)致在出現(xiàn)特定噪聲的數(shù)據(jù)集中建模時容易出現(xiàn)過擬合,在去土壤掩膜處理的紋理特征數(shù)據(jù)中,出現(xiàn)了陰影的噪聲數(shù)據(jù),對隨機森林回歸造成了干擾,并且在用支持向量機進行監(jiān)督分類時,部分土壤并沒有被掩膜處理,也對隨機森林回歸造成了干擾。而BP 神經(jīng)網(wǎng)絡(luò)存在局部極小化的問題,采取局部搜索的方法的容易陷入局部最小值, 并且BP神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)自由度很大,目前并沒有完整的理論指導(dǎo),因此本文采用的隱含層結(jié)構(gòu)可能較適合于有完整信息的數(shù)據(jù)集,而對噪聲的容忍能力較高,去土壤去陰影掩膜刪去過多有效信息,所以在去土壤去陰影掩膜處理的紋理特征中效果較差。

    3 結(jié)論

    (1)本文通過對不同分辨率下獲得的紋理特征進行土壤含水率的反演,發(fā)現(xiàn)在拍攝分辨率為10、40 mm 時,反演模型的擬合優(yōu)度較好。但分辨率在40 mm 時對不同反演方法表現(xiàn)更穩(wěn)定。在實際應(yīng)用中,過高分辨率下工作效率低,不能大范圍進行遙感監(jiān)測,為工作帶來不便。故本文認(rèn)為無人機遙感監(jiān)測土壤水分的最佳分辨率為40 mm。

    (2)30 mm-NM-RF 模型為改組數(shù)據(jù)反演冬小麥土壤含水率的最佳反演模型,但在30 mm 分辨率下其他模型反演效果不佳。

    (3)對于3 種反演模型,表現(xiàn)最好的是NM-RF 模型,擬合優(yōu)度較高且在各分辨率下表現(xiàn)穩(wěn)定。綜合各掩膜處理方法表現(xiàn)最好反演模型是BPNN 模型,在大部分情況下表現(xiàn)穩(wěn)定有較好的反演效果,為機器學(xué)習(xí)在遙感監(jiān)測的應(yīng)用驗證了可行性。

    (4)無掩膜處理和去土壤掩膜處理的結(jié)果優(yōu)于去土壤去陰影掩膜,且去土壤去陰影掩膜處理下分辨率較低時,反演效果較差。

    猜你喜歡
    特征模型
    一半模型
    抓住特征巧觀察
    重要模型『一線三等角』
    新型冠狀病毒及其流行病學(xué)特征認(rèn)識
    重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
    如何表達“特征”
    不忠誠的四個特征
    抓住特征巧觀察
    3D打印中的模型分割與打包
    FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
    精品亚洲乱码少妇综合久久| 久久国产精品男人的天堂亚洲| 另类精品久久| 人人妻,人人澡人人爽秒播 | 十分钟在线观看高清视频www| 黄片小视频在线播放| 男女边摸边吃奶| 亚洲精品美女久久久久99蜜臀 | 人人妻人人添人人爽欧美一区卜| 欧美亚洲 丝袜 人妻 在线| 亚洲国产欧美一区二区综合| 高清黄色对白视频在线免费看| videosex国产| 精品一区二区免费观看| av电影中文网址| 18禁观看日本| 久久鲁丝午夜福利片| 在线观看免费日韩欧美大片| 两个人看的免费小视频| 老汉色∧v一级毛片| 1024香蕉在线观看| 欧美精品一区二区免费开放| 欧美另类一区| 777久久人妻少妇嫩草av网站| 精品亚洲成国产av| 精品国产露脸久久av麻豆| 亚洲欧洲国产日韩| 色播在线永久视频| 999精品在线视频| 亚洲第一青青草原| 如日韩欧美国产精品一区二区三区| 精品一区在线观看国产| 免费黄色在线免费观看| 亚洲,欧美精品.| 亚洲精品日韩在线中文字幕| 亚洲情色 制服丝袜| 精品少妇一区二区三区视频日本电影 | 水蜜桃什么品种好| 久久精品熟女亚洲av麻豆精品| 天堂中文最新版在线下载| 亚洲av男天堂| 亚洲国产最新在线播放| 国产xxxxx性猛交| 一级毛片黄色毛片免费观看视频| 亚洲综合色网址| 中文字幕精品免费在线观看视频| 亚洲国产看品久久| 满18在线观看网站| 色综合欧美亚洲国产小说| 老司机亚洲免费影院| 亚洲国产毛片av蜜桃av| 男女边摸边吃奶| www.熟女人妻精品国产| 国产精品久久久久久精品古装| 亚洲精华国产精华液的使用体验| 亚洲av中文av极速乱| 大片电影免费在线观看免费| 免费黄色在线免费观看| 高清视频免费观看一区二区| 99热国产这里只有精品6| 国产免费视频播放在线视频| 国产日韩欧美亚洲二区| 制服丝袜香蕉在线| 亚洲综合色网址| 丝袜脚勾引网站| 久久狼人影院| 国产成人精品久久二区二区91 | 日日啪夜夜爽| 国产高清不卡午夜福利| 自拍欧美九色日韩亚洲蝌蚪91| 无限看片的www在线观看| 精品亚洲成国产av| 亚洲精品自拍成人| 一二三四中文在线观看免费高清| 国产亚洲av片在线观看秒播厂| 最近中文字幕2019免费版| 免费不卡黄色视频| 9热在线视频观看99| 国产精品香港三级国产av潘金莲 | 国产成人啪精品午夜网站| 久久久久精品性色| 亚洲国产精品999| 黄片无遮挡物在线观看| 国产1区2区3区精品| 午夜精品国产一区二区电影| 亚洲欧美精品自产自拍| 久久久国产一区二区| a 毛片基地| 热99久久久久精品小说推荐| 如日韩欧美国产精品一区二区三区| 久久青草综合色| 纯流量卡能插随身wifi吗| 欧美激情极品国产一区二区三区| 桃花免费在线播放| 下体分泌物呈黄色| 亚洲,欧美,日韩| 大香蕉久久成人网| 久久99热这里只频精品6学生| 欧美 日韩 精品 国产| 欧美激情极品国产一区二区三区| 中文字幕色久视频| 亚洲五月色婷婷综合| 熟女av电影| 如日韩欧美国产精品一区二区三区| 中文字幕av电影在线播放| 久久人人97超碰香蕉20202| 久久久久久久久久久免费av| 色综合欧美亚洲国产小说| 一级a爱视频在线免费观看| 久久毛片免费看一区二区三区| 汤姆久久久久久久影院中文字幕| 在线观看免费高清a一片| 人人妻,人人澡人人爽秒播 | 精品国产乱码久久久久久小说| 精品少妇内射三级| 欧美变态另类bdsm刘玥| 99香蕉大伊视频| 中文字幕高清在线视频| 日韩电影二区| 男人操女人黄网站| 精品久久久精品久久久| 久热这里只有精品99| 又大又爽又粗| 久久韩国三级中文字幕| 男人添女人高潮全过程视频| 街头女战士在线观看网站| 最近中文字幕高清免费大全6| 久久久久久久国产电影| 永久免费av网站大全| 激情五月婷婷亚洲| 亚洲伊人色综图| 免费女性裸体啪啪无遮挡网站| 夫妻午夜视频| 韩国av在线不卡| 亚洲精品美女久久久久99蜜臀 | 丰满少妇做爰视频| 久久天堂一区二区三区四区| 亚洲美女视频黄频| 少妇人妻精品综合一区二区| 久久鲁丝午夜福利片| 女人精品久久久久毛片| 亚洲一区中文字幕在线| 老汉色av国产亚洲站长工具| 成年女人毛片免费观看观看9 | 丝袜在线中文字幕| 精品国产露脸久久av麻豆| 欧美中文综合在线视频| 飞空精品影院首页| 人成视频在线观看免费观看| 91精品国产国语对白视频| videos熟女内射| 永久免费av网站大全| 中文字幕最新亚洲高清| 精品亚洲乱码少妇综合久久| 大香蕉久久成人网| 免费日韩欧美在线观看| 久久久欧美国产精品| 一区在线观看完整版| 亚洲国产欧美网| 亚洲第一av免费看| 日本色播在线视频| 国产老妇伦熟女老妇高清| 欧美激情 高清一区二区三区| 在线观看免费午夜福利视频| 不卡av一区二区三区| 午夜福利一区二区在线看| 日本色播在线视频| 丝袜脚勾引网站| 久久鲁丝午夜福利片| 黄片无遮挡物在线观看| 性色av一级| 国产精品99久久99久久久不卡 | 国精品久久久久久国模美| a 毛片基地| 秋霞伦理黄片| 99香蕉大伊视频| 国产激情久久老熟女| 热99国产精品久久久久久7| 美女午夜性视频免费| 男女之事视频高清在线观看 | 9色porny在线观看| 伊人久久国产一区二区| 最近的中文字幕免费完整| 亚洲av成人不卡在线观看播放网 | 日韩视频在线欧美| 亚洲综合精品二区| 精品国产一区二区久久| 亚洲精品久久久久久婷婷小说| 老鸭窝网址在线观看| 精品福利永久在线观看| 国产精品 国内视频| a级毛片在线看网站| 亚洲人成电影观看| 视频区图区小说| 最新在线观看一区二区三区 | 最新的欧美精品一区二区| 校园人妻丝袜中文字幕| 欧美在线黄色| 少妇被粗大猛烈的视频| 亚洲欧美日韩另类电影网站| 亚洲成人av在线免费| 国产成人啪精品午夜网站| 国产精品国产三级国产专区5o| av有码第一页| 制服人妻中文乱码| 欧美在线黄色| 精品一区二区三区四区五区乱码 | 一级黄片播放器| 国产欧美亚洲国产| 国产精品欧美亚洲77777| 国产成人欧美| 伊人亚洲综合成人网| 校园人妻丝袜中文字幕| 秋霞伦理黄片| 成人国语在线视频| 欧美成人精品欧美一级黄| 丝袜美腿诱惑在线| 久久久久久久精品精品| 精品人妻在线不人妻| 欧美人与性动交α欧美精品济南到| 午夜91福利影院| 一区二区日韩欧美中文字幕| 精品一区二区三区四区五区乱码 | 晚上一个人看的免费电影| 香蕉丝袜av| av天堂久久9| 老熟女久久久| 久久性视频一级片| 久久久久网色| 波野结衣二区三区在线| 我的亚洲天堂| 日韩一区二区视频免费看| 免费观看a级毛片全部| 又黄又粗又硬又大视频| 日韩,欧美,国产一区二区三区| 久久久久视频综合| 国产欧美日韩一区二区三区在线| 亚洲精品自拍成人| 不卡av一区二区三区| 欧美人与善性xxx| 国产精品亚洲av一区麻豆 | 午夜老司机福利片| 国产亚洲一区二区精品| 在线看a的网站| xxxhd国产人妻xxx| 赤兔流量卡办理| 午夜老司机福利片| 黄色 视频免费看| 亚洲一码二码三码区别大吗| 欧美乱码精品一区二区三区| 看免费成人av毛片| 久久久久精品国产欧美久久久 | 精品卡一卡二卡四卡免费| 免费黄网站久久成人精品| 少妇精品久久久久久久| 亚洲国产av影院在线观看| 国产福利在线免费观看视频| 一级片免费观看大全| 青草久久国产| 成年女人毛片免费观看观看9 | 欧美精品人与动牲交sv欧美| 少妇被粗大猛烈的视频| 亚洲一级一片aⅴ在线观看| 久久青草综合色| 亚洲第一av免费看| 精品久久久久久电影网| 亚洲精品中文字幕在线视频| 色精品久久人妻99蜜桃| 不卡av一区二区三区| 国产一区二区激情短视频 | 国产一区二区在线观看av| 久热爱精品视频在线9| 免费高清在线观看日韩| 中文字幕av电影在线播放| 午夜免费观看性视频| 精品一区二区三区四区五区乱码 | 伊人久久国产一区二区| 操出白浆在线播放| 亚洲精品成人av观看孕妇| 美女脱内裤让男人舔精品视频| 国产一区二区 视频在线| 国产免费又黄又爽又色| 成人国语在线视频| 十八禁高潮呻吟视频| 成人手机av| 老司机影院成人| 欧美最新免费一区二区三区| 免费观看性生交大片5| 日韩一区二区视频免费看| 女人精品久久久久毛片| 九九爱精品视频在线观看| 亚洲精品久久久久久婷婷小说| 不卡av一区二区三区| 成人免费观看视频高清| 久久久久人妻精品一区果冻| 考比视频在线观看| 成人三级做爰电影| av在线老鸭窝| 日本vs欧美在线观看视频| 性色av一级| 亚洲精品aⅴ在线观看| 日韩欧美精品免费久久| 日韩熟女老妇一区二区性免费视频| 亚洲av电影在线观看一区二区三区| 丝袜喷水一区| 性少妇av在线| www.自偷自拍.com| 中文字幕人妻熟女乱码| 老汉色∧v一级毛片| 久久久亚洲精品成人影院| 国语对白做爰xxxⅹ性视频网站| 十八禁网站网址无遮挡| 欧美久久黑人一区二区| 国产伦理片在线播放av一区| www.熟女人妻精品国产| 精品一品国产午夜福利视频| 七月丁香在线播放| 亚洲伊人色综图| 在线观看免费午夜福利视频| 老汉色∧v一级毛片| 国产精品国产三级专区第一集| 久久久久精品性色| 中文天堂在线官网| 亚洲精品日本国产第一区| 90打野战视频偷拍视频| 久久 成人 亚洲| 欧美人与善性xxx| 国产视频首页在线观看| 80岁老熟妇乱子伦牲交| 免费久久久久久久精品成人欧美视频| 97精品久久久久久久久久精品| 久久久久久免费高清国产稀缺| 啦啦啦在线观看免费高清www| 少妇猛男粗大的猛烈进出视频| 性色av一级| 亚洲成人免费av在线播放| 自线自在国产av| 丝瓜视频免费看黄片| 久久狼人影院| 五月开心婷婷网| 别揉我奶头~嗯~啊~动态视频 | 丁香六月天网| 亚洲七黄色美女视频| 欧美日韩亚洲国产一区二区在线观看 | 亚洲色图综合在线观看| 晚上一个人看的免费电影| 精品酒店卫生间| 青草久久国产| 热re99久久精品国产66热6| 色精品久久人妻99蜜桃| 一级黄片播放器| 国产精品 国内视频| av女优亚洲男人天堂| 国产免费福利视频在线观看| 99久国产av精品国产电影| 国产成人精品在线电影| 精品少妇内射三级| 老司机影院毛片| 国产成人免费无遮挡视频| 日本黄色日本黄色录像| 久久国产亚洲av麻豆专区| 久久精品国产亚洲av涩爱| 青春草亚洲视频在线观看| 国产又色又爽无遮挡免| 久久ye,这里只有精品| h视频一区二区三区| xxx大片免费视频| 精品亚洲成国产av| 精品午夜福利在线看| 久久久国产欧美日韩av| 一本色道久久久久久精品综合| 亚洲伊人色综图| 精品久久久久久电影网| 亚洲成人免费av在线播放| 亚洲免费av在线视频| 人人妻,人人澡人人爽秒播 | 日本猛色少妇xxxxx猛交久久| 搡老岳熟女国产| 亚洲成国产人片在线观看| 十八禁高潮呻吟视频| 99久久精品国产亚洲精品| 97在线人人人人妻| 日韩中文字幕视频在线看片| 国产又爽黄色视频| 亚洲综合色网址| 国产免费又黄又爽又色| 国产亚洲午夜精品一区二区久久| 亚洲视频免费观看视频| 国产日韩一区二区三区精品不卡| 国产成人一区二区在线| 啦啦啦在线免费观看视频4| 岛国毛片在线播放| 中文字幕高清在线视频| 美女脱内裤让男人舔精品视频| 亚洲综合精品二区| 男女免费视频国产| 丝袜脚勾引网站| 亚洲精品aⅴ在线观看| 黄色视频在线播放观看不卡| 男女边摸边吃奶| 久热爱精品视频在线9| 亚洲精品久久久久久婷婷小说| 韩国高清视频一区二区三区| 欧美日韩福利视频一区二区| 美女视频免费永久观看网站| 亚洲av欧美aⅴ国产| 国产片内射在线| 国产一区二区在线观看av| 国产精品一区二区精品视频观看| 十八禁网站网址无遮挡| 啦啦啦在线观看免费高清www| 婷婷成人精品国产| 中文字幕另类日韩欧美亚洲嫩草| 亚洲av电影在线观看一区二区三区| 欧美av亚洲av综合av国产av | 色婷婷av一区二区三区视频| 天堂俺去俺来也www色官网| 男女免费视频国产| 亚洲国产欧美网| 亚洲精品美女久久久久99蜜臀 | 丰满乱子伦码专区| 乱人伦中国视频| 欧美精品一区二区免费开放| 黄色 视频免费看| 精品国产乱码久久久久久小说| 亚洲av在线观看美女高潮| 亚洲av男天堂| 欧美黑人精品巨大| 亚洲情色 制服丝袜| 日韩,欧美,国产一区二区三区| 99久久综合免费| 欧美日韩一区二区视频在线观看视频在线| 在线 av 中文字幕| 亚洲av国产av综合av卡| 999精品在线视频| 国产精品秋霞免费鲁丝片| 在线观看人妻少妇| 欧美国产精品一级二级三级| 一本一本久久a久久精品综合妖精| 日本av手机在线免费观看| 午夜福利免费观看在线| 国产老妇伦熟女老妇高清| 一级片免费观看大全| 在线观看www视频免费| 日本一区二区免费在线视频| 国产亚洲精品第一综合不卡| 2018国产大陆天天弄谢| 亚洲,欧美精品.| 一二三四中文在线观看免费高清| 欧美日韩亚洲综合一区二区三区_| 国产精品熟女久久久久浪| 母亲3免费完整高清在线观看| 亚洲欧美激情在线| 少妇人妻 视频| 日本欧美视频一区| 中文字幕最新亚洲高清| 黄色视频不卡| 无限看片的www在线观看| 午夜福利网站1000一区二区三区| 18在线观看网站| 下体分泌物呈黄色| 如日韩欧美国产精品一区二区三区| 久久久久久人妻| 18禁裸乳无遮挡动漫免费视频| 午夜免费观看性视频| 观看美女的网站| 老司机影院毛片| av女优亚洲男人天堂| 久久精品亚洲av国产电影网| 伊人久久大香线蕉亚洲五| 国产精品亚洲av一区麻豆 | 日本爱情动作片www.在线观看| 国产免费现黄频在线看| 在现免费观看毛片| 两性夫妻黄色片| 亚洲国产精品国产精品| 一级片'在线观看视频| 国产又色又爽无遮挡免| 91精品伊人久久大香线蕉| 国产免费福利视频在线观看| 精品国产超薄肉色丝袜足j| 国产一区二区在线观看av| 国产精品香港三级国产av潘金莲 | 999精品在线视频| 操美女的视频在线观看| 成人午夜精彩视频在线观看| 丁香六月欧美| 纵有疾风起免费观看全集完整版| 国产乱来视频区| 国产欧美日韩综合在线一区二区| 天天躁狠狠躁夜夜躁狠狠躁| 成人18禁高潮啪啪吃奶动态图| 女性生殖器流出的白浆| 成年人午夜在线观看视频| 两性夫妻黄色片| 新久久久久国产一级毛片| www.熟女人妻精品国产| 18禁国产床啪视频网站| 国产成人免费无遮挡视频| 国产av码专区亚洲av| av电影中文网址| 久久免费观看电影| 纯流量卡能插随身wifi吗| 在线观看三级黄色| 亚洲精品美女久久久久99蜜臀 | 亚洲国产毛片av蜜桃av| 亚洲av电影在线观看一区二区三区| 亚洲av电影在线进入| 亚洲一区中文字幕在线| 如日韩欧美国产精品一区二区三区| 精品亚洲成a人片在线观看| 久久久久精品国产欧美久久久 | 日本av手机在线免费观看| av国产久精品久网站免费入址| 看十八女毛片水多多多| 午夜免费鲁丝| 亚洲精品久久成人aⅴ小说| 两性夫妻黄色片| 欧美国产精品va在线观看不卡| 国产av码专区亚洲av| 国产在视频线精品| 欧美日韩亚洲综合一区二区三区_| av国产久精品久网站免费入址| 欧美日韩综合久久久久久| 亚洲欧美清纯卡通| 国产成人精品久久久久久| 国产老妇伦熟女老妇高清| 成人毛片60女人毛片免费| 在现免费观看毛片| 免费观看av网站的网址| 欧美人与性动交α欧美软件| 满18在线观看网站| 亚洲精品日本国产第一区| 国产有黄有色有爽视频| 国产精品久久久久成人av| 国产成人精品久久二区二区91 | 婷婷色av中文字幕| 香蕉丝袜av| 国产免费一区二区三区四区乱码| 久久女婷五月综合色啪小说| 熟女少妇亚洲综合色aaa.| videos熟女内射| 亚洲精品乱久久久久久| 日韩一区二区视频免费看| 国产亚洲最大av| 亚洲欧美成人精品一区二区| av天堂久久9| 国产淫语在线视频| 亚洲综合色网址| 精品一品国产午夜福利视频| 免费av中文字幕在线| 考比视频在线观看| 丰满迷人的少妇在线观看| 国产伦理片在线播放av一区| 日本av手机在线免费观看| 女人久久www免费人成看片| 亚洲成人手机| 国产高清不卡午夜福利| 亚洲欧美一区二区三区久久| 叶爱在线成人免费视频播放| 色婷婷久久久亚洲欧美| 交换朋友夫妻互换小说| 国产精品国产三级国产专区5o| 成人黄色视频免费在线看| 国产欧美日韩一区二区三区在线| 日韩制服骚丝袜av| 精品亚洲成a人片在线观看| 中文精品一卡2卡3卡4更新| 久久精品熟女亚洲av麻豆精品| 国产片内射在线| 人体艺术视频欧美日本| 免费看av在线观看网站| 久久久国产一区二区| 久久久国产精品麻豆| 麻豆乱淫一区二区| 十分钟在线观看高清视频www| 在线精品无人区一区二区三| 久久精品亚洲熟妇少妇任你| 别揉我奶头~嗯~啊~动态视频 | 伊人久久国产一区二区| 国产精品女同一区二区软件| 在线看a的网站| 久久久久久人妻| 制服人妻中文乱码| 亚洲欧美一区二区三区黑人| 午夜精品国产一区二区电影| 亚洲精品国产色婷婷电影| 中文字幕人妻丝袜制服| 中文字幕亚洲精品专区| 男女边吃奶边做爰视频| 欧美日韩亚洲综合一区二区三区_| 亚洲精品国产区一区二| 亚洲国产欧美在线一区| 久久久久久久久久久免费av| 黄网站色视频无遮挡免费观看| 777米奇影视久久| 久久精品人人爽人人爽视色| 久久人人爽人人片av| 久久久久国产一级毛片高清牌| 亚洲国产日韩一区二区| 久热爱精品视频在线9| 国产成人欧美在线观看 | 在线观看免费午夜福利视频| 国产精品国产av在线观看| 亚洲自偷自拍图片 自拍| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲人成网站在线观看播放| 欧美人与性动交α欧美软件| 操出白浆在线播放| 欧美亚洲 丝袜 人妻 在线| 天天躁狠狠躁夜夜躁狠狠躁| 国产成人免费无遮挡视频| 十八禁人妻一区二区| 国产日韩欧美在线精品| 嫩草影院入口|