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

    基于衰減消去蜻蜓算法的小麥粉蛋白質(zhì)近紅外特征波長優(yōu)選

    2022-08-02 03:11:02吳彩娥熊智新
    食品科學(xué) 2022年14期
    關(guān)鍵詞:小麥粉波長光譜

    陳 勇,吳彩娥,熊智新*

    (南京林業(yè)大學(xué)輕工與食品學(xué)院,江蘇 南京 210037)

    小麥粉中蛋白質(zhì)質(zhì)量分?jǐn)?shù)約為7%~15%,決定著小麥粉的營養(yǎng)品質(zhì)和加工品質(zhì),所以小麥粉蛋白含量的檢測具有重要意義。小麥蛋白含量測定有化學(xué)法和物理法,前者雖然具有較高的準(zhǔn)確度,但需要對實(shí)驗(yàn)材料進(jìn)行繁雜的化學(xué)處理,比較耗時(shí)費(fèi)力,同時(shí)實(shí)驗(yàn)使用的化學(xué)藥品會(huì)造成一定的環(huán)境污染。與傳統(tǒng)的化學(xué)方法相比,近紅外光譜分析是一種綠色、無損快速檢測技術(shù),已經(jīng)廣泛應(yīng)用于農(nóng)業(yè)、食品、醫(yī)藥、化工等行業(yè)。目前,國內(nèi)外學(xué)者利用近紅外光譜法對小麥粉中蛋白質(zhì)、灰分和水分等含量進(jìn)行了檢測分析。然而,近紅外光譜有吸收信號弱、譜峰重疊嚴(yán)重以及易受外界環(huán)境干擾等缺點(diǎn),并且一條光譜往往包含著數(shù)量眾多的波長點(diǎn),這給建立高質(zhì)量的預(yù)測模型帶來了很大的挑戰(zhàn)。通過合適的波長選擇方法可以在眾多近紅外光譜波長中篩選出具有特征信息的波長,并對數(shù)據(jù)進(jìn)行有效地降維,降低建立預(yù)測模型的復(fù)雜程度,在一定程度上可提高預(yù)測模型的準(zhǔn)確性。

    常見的波長選擇方法有無信息變量消除法、競爭自適應(yīng)加權(quán)采樣法、連續(xù)投影算法以及群體智能優(yōu)化算法。群體智能優(yōu)化算法具有強(qiáng)大的全局搜索能力,使其在特征變量篩選方面潛力巨大,其中遺傳算法、粒子群算法、灰狼算法等在波長選擇上已有很多成功研究案例。其中,蜻蜓算法(dragonfly algorithm,DA)是Mirjalili等在2016年通過對自然界蜻蜓行為進(jìn)行觀察、總結(jié)和抽象后,提出的一種新智能群體優(yōu)化算法。之后,分別采用連續(xù)DA、二進(jìn)制DA(binary DA,BDA)及多目標(biāo)DA,通過對幾類典型函數(shù)進(jìn)行優(yōu)化驗(yàn)證了算法的有效性。2019年,Chen Yuanyuan等首次把BDA應(yīng)用于近紅外特征波長篩選,并比較了單群BDA(single-BDA)、多群BDA(multi-BDA)、基于集成學(xué)習(xí)BDA(ensemble learning-based BDA,ELB-BDA),表明后兩種算法可以明顯提高波長選擇的穩(wěn)定性以及分析模型的泛化能力。然而,使用BDA進(jìn)行波長篩選最終得到的波長位置和波長數(shù)量具有強(qiáng)隨機(jī)性。盡管使用multi-BDA和ELB-BDA在一定程度上降低了波長數(shù)量并進(jìn)一步篩選出特征波長,提高了算法的穩(wěn)定性,但是這兩種策略操作復(fù)雜,需要多次進(jìn)行single-BDA運(yùn)算,導(dǎo)致運(yùn)行時(shí)間遠(yuǎn)多于single-BDA,不利于快速建模。并且由于multi-BDA和ELB-BDA的內(nèi)核還是single-BDA算法,選出的波長分布隨機(jī)性仍然較大。為了加快變量選擇速度,Yun Yonghuan等基于達(dá)爾文自然進(jìn)化論中簡單有效的“適者生存”原則,在變量組合集群分析的變量選擇策略中采用指數(shù)衰減函數(shù)不斷縮小搜尋空間以加快算法速度。借鑒這一思想,在single-BDA算法中引入指數(shù)衰減函數(shù),使算法能在迭代過程中由快趨緩地剔除無用變量,有利于保留重要變量,但不利的是由于衰減曲線末端趨向平緩可能造成無效計(jì)算,或末端迭代刪除變量數(shù)大于1而漏掉某一組或多組變量組合中可能存在的最優(yōu)變量。為此,本研究嘗試用線性衰減函數(shù)代替末期的指數(shù)衰減函數(shù),保證算法末期每迭代一次剔除一個(gè)變量。Single-BDA算法的改進(jìn)因此分為兩個(gè)階段:第1階段使用指數(shù)衰減函數(shù)對變量進(jìn)行快速挑選;第2階段采用線性衰減函數(shù)對變量進(jìn)行精細(xì)挑選。本研究把改進(jìn)算法命名為衰減消去BDA(attenuation elimination-BDA,AEBDA),以期能從小麥粉近紅外光譜中盡快挑選出數(shù)量少且穩(wěn)定的特征波長,并建立精度較高的蛋白質(zhì)近紅外分析模型。

    1 材料與方法

    1.1 材料

    市購不同品牌以及不同批次的160 個(gè)小麥粉的樣品,置于保鮮袋中于室溫(20~23 ℃)保存。

    1.2 儀器與設(shè)備

    Micro NIR Pro 1700便攜式近紅外光譜儀 美國Viavi Solutions公司;D200杜馬斯定氮儀 海能未來科技有限公司。

    1.3 方法

    1.3.1 小麥粉蛋白質(zhì)含量的測定

    參照GB 5009.5ü2016《食品中蛋白質(zhì)的測定》中的燃燒法進(jìn)行測定。

    1.3.2 近紅外光譜采集

    對小麥粉樣品在室溫(20~23 ℃)環(huán)境下不做前處理進(jìn)行光譜掃描。光譜儀機(jī)身采用金屬試管架夾持固定,探頭向下垂直對準(zhǔn)深1 cm的圓盤樣品池,樣品池頂部與探頭底部相距1 cm。采集光譜時(shí),小麥粉樣品鋪平樣品池,按120°的角度間隔采集,得到3 條不同檢測點(diǎn)的光譜,取平均作為該樣品的最終采集光譜。儀器的波長范圍為908~1 650 nm,光譜分辨率為6.25 nm,有125 個(gè)波長通道,并用與儀器配套的測控軟件MicroNIRTM Pro v2.5采集和儲(chǔ)存信號。

    1.3.3 預(yù)處理方法

    常見的預(yù)處理方法有移動(dòng)平滑法(moving average filter,MAF)、卷積平滑法(Savitzky-Golay filter,SGF)、一階導(dǎo)數(shù)(1st D)、標(biāo)準(zhǔn)正態(tài)變量變換(standard normal variate transformation,SNV)等。使用這4 種預(yù)處理方法的3 種組合(MAF與SNV、SGF與SNV、1st D與SNV,預(yù)處理窗口寬度為5)以求達(dá)到最好的建模預(yù)測效果。

    1.3.4 建模與模型評估方法

    建模方法采用偏最小二乘回歸(partial least square regression,PLSR)算法,交互驗(yàn)證采用留一法;single-BDA和AE-BDA算法中采用的適應(yīng)度函數(shù)為PLSR回歸建模中交互驗(yàn)證標(biāo)準(zhǔn)偏差(root mean square error of cross validation,RMSECV)。

    使用RMSECV、預(yù)測標(biāo)準(zhǔn)偏差(root mean square errors of prediction,RMSEP)、決定系數(shù)()對所建立的校正模型進(jìn)行評估,計(jì)算方法如下:

    式(1)~(3)中:y為第個(gè)樣品化學(xué)值的測定值;y為預(yù)測過程中第個(gè)樣品的預(yù)測值;為所用樣品集樣品數(shù);為化學(xué)值測定值的均值。

    在建模預(yù)測過程中,越接近1,表示模型的回歸或者預(yù)測結(jié)果越好;如果為負(fù)值,表示擬合效果極差。RMSECV與RMSEP的值越小表示所建立的模型的穩(wěn)定性與精確度越好。

    1.3.5 波長變量選擇方法

    1.3.5.1 Single-BDA波長選擇

    蜻蜓通常通過5 種主要策略改變它們的位置:避撞、結(jié)隊(duì)、聚集、覓食、躲避。應(yīng)用到波長選擇方法中,食物即最優(yōu)解,敵人即極差解,通過不斷改變位置即迭代更新運(yùn)算直至到達(dá)食物的位置即求得最優(yōu)解。據(jù)此建立的位置更新計(jì)算方法如表1所示。

    表1 蜻蜓5 種主要位置更新策略的數(shù)學(xué)建模Table 1 Mathematical modeling of five major position shift strategies of dragonfly

    采用single-BDA算法進(jìn)行波長選擇時(shí),位置向量X的每個(gè)元素值只能是0或1,因此蜻蜓位置更新并非如表1所示的直接在原有位置X后加上更新速度值Δ,而是只能在0和1之間切換。這就需要對Δ進(jìn)行連續(xù)域到離散域的轉(zhuǎn)換。最簡單、最有效的方法是采用傳遞函數(shù)(式(4)),以Δ作為函數(shù)輸入,返回一個(gè)[0,1]之間的數(shù)值,表示位置變化的概率。

    計(jì)算所有個(gè)體位置的變化概率后,應(yīng)用式(5)更新蜻蜓在空間中的搜索位置:

    式(5)中:為0到1之間的隨機(jī)數(shù),負(fù)號表示邏輯取反運(yùn)算。

    1.3.5.2 AE-BDA波長選擇

    AE-BDA算法是經(jīng)single-BDA算法結(jié)合指數(shù)和線性衰減函數(shù)改進(jìn)而來,其在波長篩選過程中主要分為前后兩個(gè)部分:快速挑選和精細(xì)挑選,兩部分流程如圖1所示。設(shè)定最優(yōu)波長組合比例,AE-BDA算法在迭代運(yùn)行過程中,統(tǒng)計(jì)每次迭代single-BDA計(jì)算完畢產(chǎn)生的所有個(gè)體中前h個(gè)最優(yōu)RMSECV個(gè)體中波長出現(xiàn)的次數(shù)。出現(xiàn)的頻率越高,表示波長越重要。然后,按指數(shù)或線性衰減函數(shù)所計(jì)算的當(dāng)前AE-BDA迭代應(yīng)保留的變量數(shù),確定需要剔除的無用波長。AE-BDA算法采用均等抽樣法(equal sampling method,ESM),產(chǎn)生算法迭代中single-BDA第一代初始二進(jìn)制矩陣,保證每個(gè)波長具有相同的被選中機(jī)會(huì)。

    圖1 AE-BDA流程圖Fig.1 Flow chart of AE-BDA

    波長快速挑選部分:算法開始,在全波段中無用波長數(shù)所占比例高,隨著迭代的進(jìn)行其所占比例逐步降低,根據(jù)指數(shù)衰減函數(shù)曲線特點(diǎn),波長剔除速度由快到緩。在這一步中,設(shè)置種群數(shù)中個(gè)體數(shù)為,原始變量數(shù)為,快速篩選部分篩選變化率,算法運(yùn)行結(jié)束保留的波長個(gè)數(shù)為p(1為快速挑選部分迭代次數(shù)),如果當(dāng)前迭代變量數(shù)hrp,則快速挑選部分結(jié)束;通過不斷的迭代更新,波長數(shù)量呈指數(shù)級減少,達(dá)到所設(shè)置保留的波長個(gè)數(shù)p為止。每次迭代所保留的波長數(shù)為:

    式(6)中:為當(dāng)前迭代次數(shù);為波長篩選率;p為第次迭代后所剩的理論波長數(shù)。

    精細(xì)挑選部分:根據(jù)線性衰減函數(shù)曲線特點(diǎn),保證在算法末期每代可剔除一個(gè)變量。設(shè)置最終波長數(shù)p(2為精細(xì)挑選部分迭代次數(shù)),篩選過程中每次迭代并剔除出現(xiàn)次數(shù)最少的1 個(gè)波長,利用當(dāng)次迭代得到的波長組合建立PLSR模型,記錄校正集的RMSECV。每次迭代所保留的波長數(shù)為:

    式(7)中:為當(dāng)前迭代次數(shù)。

    當(dāng)波長數(shù)減少到最小設(shè)定值p時(shí),迭代終止。選取全部RMSECV記錄中最小的波長組合作為AE-BDA波長優(yōu)選結(jié)果。

    圖2中,曲線1為AE-BDA算法波長數(shù)量的衰減趨勢圖(以初始變量數(shù)=500 個(gè),p=30 個(gè),最終波長數(shù)p=10 個(gè)為例)。在快速挑選部分,受指數(shù)函數(shù)性質(zhì)影響,波長數(shù)量經(jīng)過1次迭代快速降低到p個(gè),然后再由精細(xì)挑選進(jìn)一步壓縮波長數(shù)量至最小波長數(shù)p個(gè);曲線2為指數(shù)衰減函數(shù)變量衰減曲線,曲線1相比較于曲線2在算法末端(pp階段)逐個(gè)波長剔除,不會(huì)導(dǎo)致某一波長點(diǎn)數(shù)被跳過。圖2的柱狀圖為每一次迭代剔除的無用變量數(shù)量,在快速挑選部分,隨著迭代次數(shù)的增加每次迭代剔除的無用變量數(shù)也呈指數(shù)減少;在精細(xì)挑選部分,每一代消去1 個(gè)無用變量,保證無用變量被精細(xì)均勻剔除。

    圖2 AE-BDA算法波長數(shù)衰減圖Fig.2 Decreasing number of variables in AE-BDA

    1.4 數(shù)據(jù)處理

    使用MATLAB 2016a軟件以及實(shí)驗(yàn)室自主開發(fā)的NIRSA 5.9.4化學(xué)計(jì)量學(xué)軟件,采用主成分分析和馬氏距離相結(jié)合的方法檢測異常樣本,主成分分析結(jié)合Kennard-Stone的方法對校正集和預(yù)測集進(jìn)行劃分。

    2 結(jié)果與分析

    2.1 小麥粉近紅外光譜圖

    圖3為小麥粉近紅外光譜部分原始光譜,共有160 條小麥粉近紅外光譜,每條光譜含125 個(gè)波長通道。

    圖3 小麥粉近紅外光譜部分原始光譜Fig.3 Selected original near-infrared spectra of wheat flour

    2.2 小麥粉蛋白質(zhì)含量

    由表2可知,所選樣本的蛋白質(zhì)含量基本覆蓋小麥粉中蛋白質(zhì)質(zhì)量分?jǐn)?shù)(7%~15%)且分布較為均勻,表明實(shí)驗(yàn)所選樣本具有代表性。

    表2 小麥粉蛋白質(zhì)含量統(tǒng)計(jì)表Table 2 Statistics of protein content of wheat flour

    2.3 樣本異常值剔除與樣本集劃分

    樣品、采集環(huán)境和儀器在一定程度上會(huì)產(chǎn)生異常樣本數(shù)據(jù),嚴(yán)重影響所建模型的穩(wěn)定性與預(yù)測能力,所以在建模前必須將異常樣本數(shù)據(jù)從集合中剔除。3 種預(yù)處理組合各自檢測出的異常樣本集的并集(共8 個(gè)樣本)作為總異常樣本集予以剔除。最后,劃分剩余的152 個(gè)樣本為校正集和預(yù)測集;共得到98 個(gè)校正集樣本,54 個(gè)預(yù)測集樣本。如表3所示,校正集與預(yù)測集的樣本分布較寬,具有良好的代表性。

    表3 校正集與預(yù)測集小麥粉蛋白質(zhì)含量統(tǒng)計(jì)表Table 3 Statistics of protein content of wheat flour in calibration and prediction sets

    2.4 預(yù)處理方法對全譜段PLSR建模及預(yù)測的影響

    由表4可知,MAF與SNV預(yù)處理組合的PLSR模型除

    表4 不同預(yù)處理方法全譜PLSR模型性能Table 4 Performance of full-band PLSR models developed using different preprocessing methods

    值略低于SGF與SNV預(yù)處理組合的PLSR模型外,各指標(biāo)均優(yōu)于其他建模方法。這說明采用MAF和SNV組合的預(yù)處理方法建立PLSR全譜校正模型可以提高小麥粉近紅外全譜段模型穩(wěn)定性和預(yù)測精度。因此,后續(xù)波長選擇算法研究中光譜都采用MAF與SNV組合進(jìn)行預(yù)處理。

    2.5 2 種算法對波長選擇建模及預(yù)測的影響

    分別用single-BDA和AE-BDA 2 種方法進(jìn)行波長篩選后建模,16 次建模及預(yù)測實(shí)驗(yàn)結(jié)果的比較列于表5。實(shí)驗(yàn)S1~S8波長選擇方法為single-BDA,實(shí)驗(yàn)A1~A8波長選擇方法為AE-BDA。其中,AE-BDA算法設(shè)置種群數(shù)為200 個(gè),衰減率為0.9,p設(shè)置為20 個(gè),p設(shè)置為10 個(gè),每次迭代中內(nèi)核single-BDA迭代次數(shù)為10 次。經(jīng)過程序運(yùn)行后AE-BDA共迭代28 次(快速挑選迭代18 次,精細(xì)挑選迭代10 次),其內(nèi)核single-BDA共迭代280 次。為了便于改進(jìn)前后方法對比,single-BDA設(shè)置種群數(shù)量為200 個(gè),迭代次數(shù)為280 次,篩選率為20%。

    由表4和表5可知,無論是single-BDA還是AE-BDA算法,挑選后的波長進(jìn)行建模及預(yù)測效果均明顯優(yōu)于全譜段建模。使用single-BDA進(jìn)行波長選擇時(shí),實(shí)驗(yàn)S1~S8所挑選出的平均波長數(shù)為28.38 個(gè),占原始波長數(shù)(125 個(gè))的22.7%;而使用AE-BDA方法進(jìn)行波長選擇時(shí),實(shí)驗(yàn)A1~A8所挑選出平均波長數(shù)為15.8 個(gè),為原始波長數(shù)的12.6%,相較于single-BDA方法篩選出建模波長數(shù)量更少。同時(shí)可以看到,AE-BDA方法挑選出的波長建模和預(yù)測評價(jià)指標(biāo)相應(yīng)都較接近,且其預(yù)測效果總體上好于single-BDA方法。盡管single-BDA方法建立的模型評價(jià)指標(biāo)及RMSECV總體上比AE-BDA好,但預(yù)測效果卻相對略差,這說明single-BDA方法選擇的波長過多,建立的模型可能存在過擬合現(xiàn)象,導(dǎo)致預(yù)測結(jié)果變差。

    表5 Single-BDA與AE-BDA建模及預(yù)測性能比較Table 5 Comparison of modeling and prediction performance between single-BDA and AE-BDA

    由圖4和圖5可知,盡管兩種方法都是采用隨機(jī)算法更新決策變量,但AE-BDA與single-BDA相比所挑選的波長點(diǎn)少,分布范圍更加集中,且AE-BDA中有3 個(gè)特征波長(1 385.07、1 508.95 nm和1 589.48 nm)在8 次實(shí)驗(yàn)中均被挑選出來,而single-BDA方法挑選出頻次最多的2 個(gè)波長(908.1、1 502.76 nm)只有7 次。這說明AE-BDA挑選波長的隨機(jī)性要弱于single-BDA方法,結(jié)果更加穩(wěn)定。

    圖4 Single-BDA(A)和AE-BDA(B)實(shí)驗(yàn)入選波長位置Fig.4 Selected wavelength locations in single-BDA (A) and AE-BDA (B) experiments

    圖5 Single-BDA(A)和AE-BDA(B)實(shí)驗(yàn)中波長出現(xiàn)頻次圖Fig.5 Frequency of appearance of wavelengths in single-BDA (A) and AE-BDA (B) experiments

    AE-BDA相較于single-BDA方法計(jì)算效率更高。由圖6可知,種群數(shù)相同(200 個(gè))的情況下,隨著迭代次數(shù)的增加,AE-BDA參與迭代計(jì)算的總波長數(shù)呈指數(shù)下降,意味計(jì)算量也呈指數(shù)級降低,計(jì)算過程逐步加速;而single-BDA算法在迭代過程中,進(jìn)入每一次迭代的總波長數(shù)是不變的,計(jì)算量不變,運(yùn)算速度保持相對恒定,運(yùn)算時(shí)間較長。

    圖6 Single-BDA和AE-BDA在迭代過程中初始波長數(shù)變化曲線Fig.6 Change in initial number of wavelength in single-BDA and AE-BDA as a function of number of iterations

    雖然近紅外光譜帶重疊嚴(yán)重且很難清晰地指出各譜峰對應(yīng)的具體基團(tuán)信息,但根據(jù)近紅外主要譜帶歸屬仍可對所選波長的合理性做出必要分析。本研究僅討論8 次AE-BDA實(shí)驗(yàn)中每次都被挑選出的3 個(gè)波長點(diǎn)。蛋白質(zhì)分子中含有氨基,NüH鍵伸縮振動(dòng)的一級倍頻吸收帶位于1 510 nm波長附近,1 508.95nm波長處為-氨基酸中的氨基的吸收峰;甲基和亞甲基中CüH鍵的伸縮振動(dòng)的一級倍頻和變形振動(dòng)的基頻吸收帶位于1 360~1 395nm波長處,1 385.07 nm波長處為CüH吸收峰;波長1 589.48 nm處位于OüH鍵的一級倍頻吸收譜區(qū),由于模型中其標(biāo)準(zhǔn)回歸系數(shù)為負(fù)且較大,可以判斷為主要由非蛋白質(zhì)分子中OüH鍵吸收引起(波長點(diǎn)位置見圖7)。這可以作為對小麥粉中淀粉、脂肪等強(qiáng)背景組分的修正。通過進(jìn)一步對照分析可以看出,所選其他波長點(diǎn)大部分也位于以上譜帶區(qū)域內(nèi)或者附近。

    圖7 出現(xiàn)8 次的波長點(diǎn)在全譜段的位置Fig.7 Locations of wavelengths that appeared eight times in the full-band spectrum

    3 結(jié) 論

    針對常規(guī)的single-BDA在近紅外光譜波長篩選過程中出現(xiàn)的問題,引入指數(shù)衰減函數(shù)和線性衰減函數(shù)對其進(jìn)行改進(jìn),形成AE-BDA,并結(jié)合PLSR建立了小麥粉蛋白質(zhì)近紅外分析模型。通過對AE-BDA所挑選出的波長進(jìn)行特征分析表明,這些波長較好地分布于蛋白質(zhì)分子的主要官能團(tuán)吸收區(qū)域。蛋白質(zhì)分子中重要的NüH鍵的吸收波長(1 508.95 nm)每次都能被挑選出來;其他反映蛋白質(zhì)分子的甲基、亞甲基以及次甲基的CüH鍵吸收波長也因出現(xiàn)頻次高而作為挑選出的特征波長用于建立小麥粉蛋白質(zhì)近紅外分析模型,提高了模型的可解釋性。

    改進(jìn)后的AE-BDA運(yùn)算速度快,所選波長穩(wěn)定,且一定程度上克服了過擬合,所建模型精度高。這為小麥粉蛋白質(zhì)近紅外光譜提供了一種可靠、有效的波長挑選方法,對于其他分析對象的近紅外光譜的變量選擇也有重要的借鑒意義。在后續(xù)研究中可以優(yōu)化算法中各因子的取值范圍,探討其對AE-BDA的影響,進(jìn)一步提高算法的運(yùn)行效率和自動(dòng)化水平。

    猜你喜歡
    小麥粉波長光譜
    HPLC-PDA雙波長法同時(shí)測定四季草片中沒食子酸和槲皮苷的含量
    小麥粉加工及選購
    糧食加工(2022年5期)2022-12-28 22:34:59
    白俄羅斯延長小麥粉和蕎麥出口禁令
    基于三維Saab變換的高光譜圖像壓縮方法
    SPE-GC法在小麥粉中典型農(nóng)藥殘留物測定中的應(yīng)用
    高光譜遙感成像技術(shù)的發(fā)展與展望
    雙波長激光治療慢性牙周炎的療效觀察
    日本研發(fā)出可完全覆蓋可見光波長的LED光源
    中國照明(2016年4期)2016-05-17 06:16:15
    星載近紅外高光譜CO2遙感進(jìn)展
    便攜式多用途光波波長測量儀
    欧美一区二区国产精品久久精品| 人人妻,人人澡人人爽秒播| 亚洲,欧美精品.| av在线观看视频网站免费| 小说图片视频综合网站| .国产精品久久| 欧美绝顶高潮抽搐喷水| 日韩免费av在线播放| 亚洲va日本ⅴa欧美va伊人久久| 三级男女做爰猛烈吃奶摸视频| 91午夜精品亚洲一区二区三区 | 搡老熟女国产l中国老女人| 久久久久久久久久成人| 欧美精品啪啪一区二区三区| 老熟妇乱子伦视频在线观看| 一进一出抽搐gif免费好疼| 亚洲精品一区av在线观看| 久久精品夜夜夜夜夜久久蜜豆| 国产精品一区二区三区四区久久| 国产精品精品国产色婷婷| 色精品久久人妻99蜜桃| 麻豆国产av国片精品| 久久草成人影院| 欧美高清性xxxxhd video| 怎么达到女性高潮| 国产熟女xx| 小说图片视频综合网站| 此物有八面人人有两片| 日韩成人在线观看一区二区三区| 亚洲av美国av| 中文字幕熟女人妻在线| 十八禁人妻一区二区| 免费在线观看成人毛片| 久久伊人香网站| 国产激情偷乱视频一区二区| 搞女人的毛片| 欧美日韩黄片免| 欧美在线黄色| 国产av不卡久久| 久久久久久久精品吃奶| 麻豆成人午夜福利视频| 精品日产1卡2卡| 精品国内亚洲2022精品成人| 欧美精品啪啪一区二区三区| 人妻夜夜爽99麻豆av| 亚洲成人久久性| 国产一区二区在线观看日韩| 男女床上黄色一级片免费看| 久久久精品欧美日韩精品| 精品久久久久久久末码| 麻豆国产av国片精品| 蜜桃久久精品国产亚洲av| 欧美潮喷喷水| 国产免费av片在线观看野外av| 五月玫瑰六月丁香| 波多野结衣高清作品| 少妇的逼好多水| а√天堂www在线а√下载| 精品欧美国产一区二区三| 在线看三级毛片| 国产成人欧美在线观看| 欧美黄色淫秽网站| 性欧美人与动物交配| 精品99又大又爽又粗少妇毛片 | 亚州av有码| 精品午夜福利在线看| 制服丝袜大香蕉在线| 国产精品日韩av在线免费观看| 亚洲精品成人久久久久久| 别揉我奶头~嗯~啊~动态视频| 午夜老司机福利剧场| 久久午夜亚洲精品久久| 国产视频内射| 欧美黄色淫秽网站| 国产精品一区二区免费欧美| 欧美激情国产日韩精品一区| 精品国产亚洲在线| 村上凉子中文字幕在线| 他把我摸到了高潮在线观看| 亚洲av.av天堂| 国产高清三级在线| 亚洲avbb在线观看| 成熟少妇高潮喷水视频| 日本在线视频免费播放| 国产真实乱freesex| 久久99热6这里只有精品| 极品教师在线视频| 中出人妻视频一区二区| 精品午夜福利在线看| 在线免费观看的www视频| 美女黄网站色视频| 嫩草影院精品99| 亚洲人成伊人成综合网2020| 伊人久久精品亚洲午夜| 精品国产三级普通话版| 亚洲成av人片免费观看| 亚洲人成网站高清观看| 国产精品久久久久久亚洲av鲁大| 成人午夜高清在线视频| 久久人妻av系列| 国产极品精品免费视频能看的| 男女床上黄色一级片免费看| 久久99热6这里只有精品| 欧美乱色亚洲激情| 国产蜜桃级精品一区二区三区| 一区二区三区免费毛片| 成人欧美大片| 午夜福利欧美成人| 亚洲最大成人手机在线| 日韩欧美在线乱码| 男女下面进入的视频免费午夜| 99久久无色码亚洲精品果冻| 高潮久久久久久久久久久不卡| 一区二区三区高清视频在线| 在线观看美女被高潮喷水网站 | 国产精品一区二区三区四区免费观看 | 中出人妻视频一区二区| 久久香蕉精品热| 欧美黄色片欧美黄色片| 美女 人体艺术 gogo| 91九色精品人成在线观看| 日韩欧美精品v在线| 成人午夜高清在线视频| 成人一区二区视频在线观看| 好男人电影高清在线观看| 噜噜噜噜噜久久久久久91| 日本黄色片子视频| 三级男女做爰猛烈吃奶摸视频| 亚洲18禁久久av| 国产精品美女特级片免费视频播放器| 亚洲欧美日韩东京热| 五月伊人婷婷丁香| 日本免费a在线| 欧美日韩乱码在线| 青草久久国产| 国产精品免费一区二区三区在线| 免费看a级黄色片| 国产真实伦视频高清在线观看 | 特大巨黑吊av在线直播| 久久午夜亚洲精品久久| 夜夜夜夜夜久久久久| 国产又黄又爽又无遮挡在线| 99久久精品一区二区三区| 韩国av一区二区三区四区| 国产精品一区二区三区四区免费观看 | 精品免费久久久久久久清纯| 日韩 亚洲 欧美在线| www.熟女人妻精品国产| 在线播放无遮挡| 99热精品在线国产| 日韩有码中文字幕| 亚洲国产精品久久男人天堂| 校园春色视频在线观看| 亚洲人成网站在线播| 啦啦啦观看免费观看视频高清| 亚洲,欧美,日韩| 国产精品一区二区性色av| 91av网一区二区| 免费电影在线观看免费观看| 哪里可以看免费的av片| 国产精品精品国产色婷婷| 麻豆国产av国片精品| 国产av一区在线观看免费| 小说图片视频综合网站| 香蕉av资源在线| 69人妻影院| 国产在视频线在精品| 国产免费一级a男人的天堂| 乱人视频在线观看| 夜夜躁狠狠躁天天躁| 午夜影院日韩av| 中亚洲国语对白在线视频| 国产色婷婷99| 久久精品久久久久久噜噜老黄 | 国产精品三级大全| 国产中年淑女户外野战色| 91狼人影院| av欧美777| 成年人黄色毛片网站| 日韩欧美在线乱码| 精品一区二区免费观看| 亚洲成a人片在线一区二区| 白带黄色成豆腐渣| 综合色av麻豆| 热99re8久久精品国产| 亚洲最大成人中文| 亚洲无线在线观看| 免费黄网站久久成人精品 | 国产精品伦人一区二区| 1000部很黄的大片| 久久天躁狠狠躁夜夜2o2o| 看片在线看免费视频| 婷婷色综合大香蕉| 啪啪无遮挡十八禁网站| 99精品在免费线老司机午夜| 国产高清视频在线播放一区| 欧美激情在线99| 久久久精品大字幕| 精品国产三级普通话版| 国产成人av教育| 婷婷精品国产亚洲av| 亚洲内射少妇av| 首页视频小说图片口味搜索| 91九色精品人成在线观看| 国产野战对白在线观看| 女同久久另类99精品国产91| 可以在线观看毛片的网站| 我的老师免费观看完整版| 亚洲美女黄片视频| 宅男免费午夜| 简卡轻食公司| 99热精品在线国产| 免费在线观看影片大全网站| 亚洲av免费在线观看| 亚洲五月婷婷丁香| 色综合婷婷激情| 亚洲精品在线观看二区| 男女之事视频高清在线观看| 精品99又大又爽又粗少妇毛片 | 天堂av国产一区二区熟女人妻| 久久久久久九九精品二区国产| 日本精品一区二区三区蜜桃| 国产亚洲精品久久久com| 最近最新中文字幕大全电影3| 亚洲av二区三区四区| 精品一区二区三区人妻视频| 黄色丝袜av网址大全| 又爽又黄无遮挡网站| 波多野结衣巨乳人妻| 亚洲最大成人中文| 婷婷精品国产亚洲av在线| 久久精品综合一区二区三区| 1024手机看黄色片| 国产久久久一区二区三区| 国产黄片美女视频| 国产伦人伦偷精品视频| 免费一级毛片在线播放高清视频| 香蕉av资源在线| 男人狂女人下面高潮的视频| 又粗又爽又猛毛片免费看| 亚洲五月天丁香| 亚洲成av人片在线播放无| 国产精品免费一区二区三区在线| 亚洲精品乱码久久久v下载方式| 成人国产一区最新在线观看| 午夜老司机福利剧场| 欧美乱色亚洲激情| 亚洲精品粉嫩美女一区| 久久久久精品国产欧美久久久| 亚洲美女黄片视频| 很黄的视频免费| 夜夜看夜夜爽夜夜摸| 嫩草影院新地址| a级一级毛片免费在线观看| 黄色女人牲交| 亚洲欧美精品综合久久99| 99久久久亚洲精品蜜臀av| 国产精品久久久久久人妻精品电影| 成人欧美大片| 长腿黑丝高跟| 天堂网av新在线| 少妇丰满av| 国产成人av教育| 亚洲 欧美 日韩 在线 免费| 亚洲精品影视一区二区三区av| 日韩欧美在线乱码| 国产成年人精品一区二区| 蜜桃久久精品国产亚洲av| av中文乱码字幕在线| 国产亚洲av嫩草精品影院| 激情在线观看视频在线高清| 观看美女的网站| 亚洲av中文字字幕乱码综合| 少妇的逼好多水| 日韩欧美国产在线观看| 精品一区二区免费观看| 日韩欧美在线二视频| 日韩成人在线观看一区二区三区| 夜夜夜夜夜久久久久| 成人毛片a级毛片在线播放| 最后的刺客免费高清国语| 在线十欧美十亚洲十日本专区| 女生性感内裤真人,穿戴方法视频| 日韩欧美国产在线观看| 精品一区二区免费观看| 91麻豆av在线| 91av网一区二区| 免费在线观看亚洲国产| 精品国产亚洲在线| 国产精品精品国产色婷婷| 97超级碰碰碰精品色视频在线观看| 国产精品不卡视频一区二区 | 午夜亚洲福利在线播放| 欧美高清性xxxxhd video| 成年版毛片免费区| 深夜a级毛片| 综合色av麻豆| 一级黄片播放器| 免费av观看视频| 美女大奶头视频| 日本免费a在线| 亚洲av电影不卡..在线观看| 成人国产综合亚洲| 黄色女人牲交| 青草久久国产| 在线国产一区二区在线| 日韩欧美国产一区二区入口| 一个人免费在线观看电影| 美女高潮的动态| 亚洲最大成人中文| 亚洲精品色激情综合| 国产美女午夜福利| 婷婷亚洲欧美| 国产高清视频在线观看网站| 日本黄色片子视频| 国产亚洲精品久久久久久毛片| 国产亚洲精品久久久久久毛片| 亚洲欧美日韩高清在线视频| .国产精品久久| 国产精品日韩av在线免费观看| 在线观看午夜福利视频| 精品久久久久久久人妻蜜臀av| 深夜a级毛片| 午夜福利视频1000在线观看| 搡老熟女国产l中国老女人| 一区福利在线观看| 亚洲精品久久国产高清桃花| 美女被艹到高潮喷水动态| 日韩中字成人| 美女 人体艺术 gogo| 国产精品1区2区在线观看.| 91麻豆精品激情在线观看国产| 国产人妻一区二区三区在| 91狼人影院| 嫩草影视91久久| 免费av不卡在线播放| 国产精品一及| 熟女人妻精品中文字幕| 国产乱人伦免费视频| 午夜福利免费观看在线| 日韩欧美国产在线观看| 男女下面进入的视频免费午夜| 在现免费观看毛片| 国产久久久一区二区三区| 欧美日韩黄片免| 深爱激情五月婷婷| 亚洲天堂国产精品一区在线| 大型黄色视频在线免费观看| 俄罗斯特黄特色一大片| 日本在线视频免费播放| 欧美日韩瑟瑟在线播放| 宅男免费午夜| 国产一区二区在线观看日韩| 亚洲精品影视一区二区三区av| 国产高清三级在线| 亚洲真实伦在线观看| 亚洲欧美日韩卡通动漫| 亚洲三级黄色毛片| 99国产极品粉嫩在线观看| 欧美成人一区二区免费高清观看| 午夜日韩欧美国产| 人人妻,人人澡人人爽秒播| 91久久精品国产一区二区成人| 观看免费一级毛片| 国内少妇人妻偷人精品xxx网站| 亚洲av第一区精品v没综合| 天堂网av新在线| 男女下面进入的视频免费午夜| 无人区码免费观看不卡| 午夜两性在线视频| 我要搜黄色片| 真实男女啪啪啪动态图| 亚洲国产精品999在线| 成人性生交大片免费视频hd| 激情在线观看视频在线高清| 欧美性猛交黑人性爽| 久久伊人香网站| 精品久久久久久,| 久久九九热精品免费| 日本熟妇午夜| 久久精品综合一区二区三区| 99久久成人亚洲精品观看| 日韩精品青青久久久久久| 精品久久久久久久末码| 亚洲成av人片在线播放无| 亚洲av成人av| x7x7x7水蜜桃| 成人特级黄色片久久久久久久| 性色av乱码一区二区三区2| a级一级毛片免费在线观看| 国产蜜桃级精品一区二区三区| 午夜老司机福利剧场| 18美女黄网站色大片免费观看| 国产大屁股一区二区在线视频| 亚洲中文字幕日韩| 久久精品影院6| 99在线视频只有这里精品首页| a级一级毛片免费在线观看| 国产精品一及| 日韩亚洲欧美综合| 黄色配什么色好看| 国产精品,欧美在线| 色综合欧美亚洲国产小说| 99视频精品全部免费 在线| 69av精品久久久久久| 国产精品电影一区二区三区| 男女下面进入的视频免费午夜| 中文字幕精品亚洲无线码一区| 亚洲国产日韩欧美精品在线观看| 日本五十路高清| 99国产精品一区二区蜜桃av| 国产伦精品一区二区三区四那| 美女黄网站色视频| 91麻豆av在线| 国产精品亚洲一级av第二区| 99国产综合亚洲精品| 久久久精品大字幕| www.色视频.com| 国产精品伦人一区二区| 欧美又色又爽又黄视频| 亚洲精华国产精华精| 国产免费男女视频| 免费av毛片视频| 国模一区二区三区四区视频| 国产亚洲精品综合一区在线观看| 国产精品久久久久久精品电影| 熟女人妻精品中文字幕| 欧美激情国产日韩精品一区| 亚州av有码| 丰满人妻熟妇乱又伦精品不卡| 亚洲乱码一区二区免费版| 热99re8久久精品国产| 中文亚洲av片在线观看爽| 深爱激情五月婷婷| 真实男女啪啪啪动态图| 久久亚洲精品不卡| 成年女人毛片免费观看观看9| 久久热精品热| 成人无遮挡网站| 国产精品98久久久久久宅男小说| 亚洲五月天丁香| 激情在线观看视频在线高清| 日韩欧美三级三区| 一进一出抽搐动态| 午夜福利在线观看免费完整高清在 | 亚洲乱码一区二区免费版| 久久久久性生活片| 欧美+日韩+精品| 亚洲内射少妇av| 天堂av国产一区二区熟女人妻| 日韩欧美三级三区| 国产久久久一区二区三区| 婷婷精品国产亚洲av在线| 青草久久国产| 国产乱人伦免费视频| 国产黄色小视频在线观看| 国产精品,欧美在线| 国内毛片毛片毛片毛片毛片| 免费av观看视频| 噜噜噜噜噜久久久久久91| 亚洲成av人片在线播放无| 乱人视频在线观看| 中文字幕高清在线视频| 国产国拍精品亚洲av在线观看| 色尼玛亚洲综合影院| 午夜福利欧美成人| 国产精品久久久久久久久免 | 一个人免费在线观看的高清视频| 国产高清三级在线| 久久婷婷人人爽人人干人人爱| 日本撒尿小便嘘嘘汇集6| 成人毛片a级毛片在线播放| 一边摸一边抽搐一进一小说| 在线a可以看的网站| 禁无遮挡网站| 99久久久亚洲精品蜜臀av| 一本一本综合久久| 国产极品精品免费视频能看的| 国产精品久久视频播放| 窝窝影院91人妻| 久久久久久久午夜电影| 国产欧美日韩精品一区二区| 国产精品免费一区二区三区在线| 在线播放国产精品三级| 九九久久精品国产亚洲av麻豆| 亚洲不卡免费看| 精品久久久久久久久亚洲 | 欧美成人免费av一区二区三区| 日本与韩国留学比较| 国产69精品久久久久777片| 日本五十路高清| 能在线免费观看的黄片| 熟女电影av网| 免费观看人在逋| 精品国内亚洲2022精品成人| 男女床上黄色一级片免费看| 久久久国产成人精品二区| 婷婷精品国产亚洲av在线| www.www免费av| 两个人的视频大全免费| 性插视频无遮挡在线免费观看| 日本成人三级电影网站| 亚洲自偷自拍三级| 欧美丝袜亚洲另类 | 搞女人的毛片| 听说在线观看完整版免费高清| 亚洲国产日韩欧美精品在线观看| 亚洲狠狠婷婷综合久久图片| 欧美一级a爱片免费观看看| 在线免费观看的www视频| 精品久久久久久久久久久久久| 久久亚洲精品不卡| 国产精品亚洲av一区麻豆| 国产精品影院久久| 熟女人妻精品中文字幕| 精品一区二区三区人妻视频| 亚洲人成电影免费在线| 97碰自拍视频| 极品教师在线视频| 久久精品久久久久久噜噜老黄 | 很黄的视频免费| 亚洲乱码一区二区免费版| 最新在线观看一区二区三区| 亚洲美女黄片视频| 中文字幕免费在线视频6| 国产主播在线观看一区二区| 国产综合懂色| 免费高清视频大片| 9191精品国产免费久久| 一级毛片久久久久久久久女| 激情在线观看视频在线高清| 一边摸一边抽搐一进一小说| 人妻夜夜爽99麻豆av| 中文字幕精品亚洲无线码一区| 色在线成人网| 内地一区二区视频在线| 97超视频在线观看视频| 国内精品久久久久精免费| 琪琪午夜伦伦电影理论片6080| 国产激情偷乱视频一区二区| 露出奶头的视频| 亚洲无线观看免费| 嫩草影院新地址| 特级一级黄色大片| av视频在线观看入口| 小说图片视频综合网站| 久9热在线精品视频| av天堂在线播放| 亚洲欧美激情综合另类| 搡老岳熟女国产| 亚洲美女搞黄在线观看 | 99精品在免费线老司机午夜| 90打野战视频偷拍视频| 国内精品一区二区在线观看| 亚洲黑人精品在线| 亚洲色图av天堂| 男女视频在线观看网站免费| 日本成人三级电影网站| 久久九九热精品免费| 日韩人妻高清精品专区| 欧美xxxx性猛交bbbb| 精品人妻1区二区| 极品教师在线视频| 香蕉av资源在线| 免费人成在线观看视频色| 久久久久久大精品| 欧美激情国产日韩精品一区| 欧美成人性av电影在线观看| 亚洲专区国产一区二区| 天堂√8在线中文| 成人一区二区视频在线观看| 成年女人永久免费观看视频| 午夜免费男女啪啪视频观看 | 久久久久久久久久成人| 免费观看人在逋| 日韩精品中文字幕看吧| 99久久无色码亚洲精品果冻| 亚洲国产精品合色在线| 91在线精品国自产拍蜜月| 网址你懂的国产日韩在线| 国产精华一区二区三区| av在线蜜桃| 成人性生交大片免费视频hd| 成人国产一区最新在线观看| 久久草成人影院| 亚洲av美国av| 日日摸夜夜添夜夜添av毛片 | 国产成人影院久久av| 中文字幕av成人在线电影| 两人在一起打扑克的视频| 欧美bdsm另类| 乱码一卡2卡4卡精品| 久久久久久国产a免费观看| 国产熟女xx| 一个人观看的视频www高清免费观看| 99视频精品全部免费 在线| 日韩 亚洲 欧美在线| 午夜福利18| 中文字幕免费在线视频6| 亚洲国产欧洲综合997久久,| 国产精品人妻久久久久久| 啪啪无遮挡十八禁网站| 动漫黄色视频在线观看| 很黄的视频免费| 一区福利在线观看| 99热这里只有是精品50| a级毛片免费高清观看在线播放| 久久久久免费精品人妻一区二区| 在线观看美女被高潮喷水网站 | 日韩欧美三级三区| 99精品在免费线老司机午夜| 深夜精品福利| 脱女人内裤的视频| 国产伦人伦偷精品视频| 亚洲 欧美 日韩 在线 免费|