亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        遙感技術(shù)在我國(guó)冬小麥產(chǎn)量估算中的應(yīng)用研究進(jìn)展

        2023-03-22 05:23:30涂巧針劉釔廷吳保升
        南方農(nóng)業(yè) 2023年1期
        關(guān)鍵詞:估產(chǎn)植被指數(shù)遙感技術(shù)

        涂巧針,李 旭,劉釔廷,吳保升

        (塔里木大學(xué)信息學(xué)院,新疆阿拉爾 843017)

        冬小麥?zhǔn)俏覈?guó)的重要經(jīng)濟(jì)作物之一,其播種面積約為全國(guó)糧食總播種面積的1/5[1]。收割之前,適時(shí)精確的估算糧食作物生產(chǎn)水平既有利于我國(guó)國(guó)家糧食戰(zhàn)略的調(diào)整和制訂,也有利于進(jìn)一步完善糧食作物管理,合理調(diào)整種植業(yè)布局,深入發(fā)掘其生產(chǎn)潛力,為我國(guó)進(jìn)一步調(diào)整糧油儲(chǔ)備結(jié)構(gòu)提供科技支撐。實(shí)現(xiàn)我國(guó)糧食安全目標(biāo),將對(duì)我國(guó)糧食市場(chǎng)的宏觀調(diào)節(jié)、糧食出口貿(mào)易發(fā)展產(chǎn)生重大意義[2]。

        由于冬小麥分布區(qū)廣泛,北方地區(qū)復(fù)雜,對(duì)播種面積、生產(chǎn)和統(tǒng)計(jì)資料的獲取一般要利用常規(guī)的估產(chǎn)手段,包括統(tǒng)計(jì)手段和傳統(tǒng)的土壤監(jiān)測(cè)手段,其速度慢、工程量大、成本高,無(wú)法滿足對(duì)冬小麥的監(jiān)測(cè)需求。遙感資料具備覆蓋面大、檢測(cè)周期短、資料豐富、現(xiàn)勢(shì)性好、代價(jià)低等優(yōu)點(diǎn),可以迅速和精確地對(duì)冬小麥進(jìn)行估產(chǎn)。遙感估產(chǎn)是一種使用衛(wèi)星傳感器記錄作物表面信號(hào),并利用資料采集與數(shù)據(jù)分析所獲取的農(nóng)作物表面光譜特征,確定作物種類,監(jiān)控作物的生長(zhǎng)過(guò)程,從而構(gòu)建將作物光譜數(shù)據(jù)與農(nóng)業(yè)生產(chǎn)過(guò)程相互聯(lián)系的新科技手段[3]。

        1 農(nóng)作物遙感估產(chǎn)研究概況

        美國(guó)首先開(kāi)展的大面積作物遙感技術(shù)估產(chǎn),效果也最為顯著。美國(guó)早在1974年就開(kāi)始實(shí)施“大規(guī)模耕地估產(chǎn)實(shí)驗(yàn)” 即LACIE 規(guī)劃(Large Area Crop Inventory and Experiment)和“利用空間遙感技術(shù)開(kāi)展農(nóng)產(chǎn)品和資源利用研究” 即AGRISARS 規(guī)劃(Agriculture and Resources Inventory Surveys through Aerospace Remote Sensing)[4]。這是遙感技術(shù)在農(nóng)業(yè)上應(yīng)用的成功典型案例,極大的推動(dòng)了農(nóng)業(yè)遙感技術(shù)的發(fā)展。此后,意大利、俄羅斯、英國(guó)、日本和其他多個(gè)國(guó)家也相繼開(kāi)展農(nóng)業(yè)遙感技術(shù)相關(guān)的研究工作。2002 年,美國(guó)農(nóng)業(yè)部與NASA 再次合作,首次將MODIS 數(shù)據(jù)應(yīng)用到了農(nóng)業(yè)遙感研究中并采用Landsat數(shù)據(jù)對(duì)多個(gè)州進(jìn)行了農(nóng)作物監(jiān)測(cè)實(shí)驗(yàn),到2010年實(shí)現(xiàn)了基于遙感的長(zhǎng)期農(nóng)業(yè)監(jiān)測(cè)。目前,經(jīng)過(guò)多年的實(shí)驗(yàn)研究,國(guó)外農(nóng)業(yè)遙感技術(shù)已經(jīng)逐步走向成熟。

        我國(guó)的農(nóng)業(yè)遙感研究起步較晚。1981 年,我國(guó)使用遙感技術(shù)成功預(yù)測(cè)了冬小麥的實(shí)際生產(chǎn)。1983 年,我國(guó)研發(fā)了估算京津冀地區(qū)冬小麥產(chǎn)量的技術(shù)方法,在此后的3年里,GIS系統(tǒng)已經(jīng)被廣泛應(yīng)用于黃淮海地區(qū)冬小麥的產(chǎn)量估算研究中。通過(guò)開(kāi)發(fā),我國(guó)遙感技術(shù)在農(nóng)作物長(zhǎng)勢(shì)監(jiān)測(cè)與產(chǎn)量估測(cè)領(lǐng)域的應(yīng)用逐步走向完善。為更好推動(dòng)我國(guó)農(nóng)業(yè)遙感科技的開(kāi)發(fā),中國(guó)科學(xué)院制定了“九五”重大科研項(xiàng)目,結(jié)合“3S”(RS,GPS,GIS)技術(shù)推動(dòng)我國(guó)的農(nóng)業(yè)遙感向著周期更短、精度更高、成本更低的方向深入發(fā)展。2008 年,中科院和中華糧網(wǎng)合作將農(nóng)業(yè)遙感獲得的產(chǎn)量成果等推向了服務(wù)應(yīng)用,使該項(xiàng)技術(shù)得到了更為廣泛的應(yīng)用[5]。

        2 冬小麥估產(chǎn)應(yīng)用的遙感資料

        冬小麥估產(chǎn)中應(yīng)用的遙感資料主要有資源衛(wèi)星資料、氣象衛(wèi)星資料、航空遙感和地面遙感資料等。1)資源衛(wèi)星資料。Landsat 系列數(shù)據(jù)憑借優(yōu)勢(shì),從20 世紀(jì)70年代開(kāi)始就成為農(nóng)業(yè)遙感領(lǐng)域應(yīng)用最多的數(shù)據(jù)源之一,極大地推動(dòng)了農(nóng)業(yè)遙感監(jiān)測(cè)技術(shù)。除此之外,CBERS 衛(wèi)星也是常見(jiàn)的資源衛(wèi)星,但由于衛(wèi)星重訪周期長(zhǎng)、作物生育期可用數(shù)據(jù)有限等原因,限制了這一類衛(wèi)星數(shù)據(jù)的用途。2)氣象衛(wèi)星資料。氣象衛(wèi)星重訪周期快,因此氣象衛(wèi)星多基于作物的季相節(jié)律信息形成時(shí)間序列數(shù)據(jù),被應(yīng)用于大規(guī)模的作物遙感監(jiān)測(cè)。AVHRR 是最早的氣象衛(wèi)星,1987 年美國(guó)基于該數(shù)據(jù)進(jìn)行農(nóng)作物估產(chǎn)試驗(yàn)。之后,美國(guó)發(fā)射中分辨率成像儀(MODIS)用于觀測(cè)全球植被,相比之前的氣象衛(wèi)星,MODIS 數(shù)據(jù)擁有了更高的空間、光譜分辨率,從而大大增強(qiáng)了對(duì)地觀測(cè)能力,為大面積作物遙感識(shí)別提供了重要的數(shù)據(jù)支撐。但由于氣象衛(wèi)星空間分辨率較低,導(dǎo)致這類數(shù)據(jù)的總體提取精度不高。3)航空遙感和地面遙感資料[6]。

        各種遙感數(shù)據(jù)既有優(yōu)點(diǎn),也有不足之處。千懷遂按照各類資料源的清晰度(空氣清晰度和時(shí)間清晰度)和空中云層,將我國(guó)小麥估產(chǎn)的資料源類型劃分為陸地衛(wèi)星型、陸地衛(wèi)星為主型、氣象衛(wèi)星為主型和氣象衛(wèi)星型[7]。根據(jù)遙感技術(shù)數(shù)據(jù)的實(shí)際使用效果考慮,AVHRR 數(shù)據(jù)主要用來(lái)觀察冬小麥的發(fā)育狀況,利用相應(yīng)的綠度指標(biāo)建立產(chǎn)量預(yù)測(cè)模型。由于高分辨率的優(yōu)勢(shì),雖然TM 數(shù)據(jù)已廣泛用作大冬小麥面積數(shù)據(jù)的獲取手段,但因?yàn)樯a(chǎn)成本較高,目前還難以實(shí)現(xiàn)基于TM數(shù)據(jù)的大面積小麥監(jiān)測(cè)與生產(chǎn)預(yù)測(cè)。

        3 冬小麥光譜特性

        根據(jù)中國(guó)西北干燥、半干旱地區(qū)的玉米生長(zhǎng)發(fā)育情況,通過(guò)田間樣點(diǎn)取樣和室內(nèi)研究的方法,研究了不同生產(chǎn)階段冬小麥生理生化數(shù)據(jù)和冠層光譜反射率[8-9]。通過(guò)研究,冬小麥在冠層原始光譜反射率趨勢(shì)與不同生產(chǎn)階段基本相同,均在可見(jiàn)光波段反射率很低,分別在550 nm 和400 nm,670 nm 波段形成了1 個(gè)反射峰和2 個(gè)吸收谷,原因主要是基于葉片葉綠素對(duì)陸光的強(qiáng)烈反射,對(duì)藍(lán)色和紅光都具有強(qiáng)烈的吸收功能。在670 nm 后反射率急劇增大,在700~900 nm 的近紅外波束中產(chǎn)生了1 個(gè)反射率大于0.34 的高反射峰值,主要在起身期到抽穗期之間,冬小麥冠層覆蓋度變化較大,在不同層位的葉內(nèi)對(duì)太陽(yáng)光產(chǎn)生了多次反射,主要是由于葉片的細(xì)胞內(nèi)部多次反射、散射引起的。在1 300~2 500 nm 的短波紅外范圍內(nèi),冠層原始反射率逐漸降低。但是,由于施氮量的提高,在700~1 000 nm之間的近紅外波束,在冬小麥冠層葉片密度增加,冠層結(jié)構(gòu)變得更加復(fù)雜,使得近紅外部分的反射總量增加,光譜反射率增大,而在可見(jiàn)光350~670 nm區(qū)域內(nèi)差異不顯著。

        4 冬小麥遙感估產(chǎn)過(guò)程

        對(duì)冬小麥的遙感技術(shù)估算過(guò)程,一般可以分成以下8個(gè)步驟[10]。

        1)遙感信息的獲取與處理。遙感信息源的選用,首先要考察能否達(dá)到工藝要求,同時(shí)還要兼顧經(jīng)濟(jì)性。

        2)遙感估產(chǎn)區(qū)劃。冬小麥生長(zhǎng)發(fā)育條件動(dòng)態(tài)觀測(cè)和估算的遙感技術(shù)的一個(gè)重大應(yīng)用,根據(jù)需要對(duì)自然環(huán)境條件、社會(huì)環(huán)境條件與冬小麥生長(zhǎng)發(fā)育條件相同的地區(qū)進(jìn)行劃分,優(yōu)化了冬小麥生長(zhǎng)發(fā)育條件的觀測(cè)與估測(cè)模型的建立。

        3)地面采集點(diǎn)布設(shè)及觀測(cè)。遙感估產(chǎn)的數(shù)據(jù)大部分來(lái)自遙感數(shù)據(jù),但僅僅依靠遙感數(shù)據(jù)獲取精確的冬小麥播種規(guī)模與數(shù)量是不夠的。必須從土壤中收集適當(dāng)?shù)臉悠?,用作遙感信息的補(bǔ)充與驗(yàn)證。

        4)建立數(shù)據(jù)庫(kù)。在遙感估產(chǎn)中,建立數(shù)據(jù)庫(kù)系統(tǒng)是一個(gè)相當(dāng)關(guān)鍵的基礎(chǔ)工作。背景數(shù)據(jù)庫(kù)主要具有2 種功能:①農(nóng)業(yè)遙感信息的劃分提供背景信息,以增強(qiáng)農(nóng)業(yè)劃分的正確性。②在難以獲得農(nóng)業(yè)遙感信息的狀況下,它可以支持模擬數(shù)據(jù)分析技術(shù),從歷史數(shù)據(jù)和現(xiàn)實(shí)調(diào)查樣本中獲取的農(nóng)作物數(shù)據(jù)加以處理,從而得出當(dāng)年的農(nóng)作物實(shí)際種植面積和產(chǎn)量數(shù)據(jù)。

        5)冬小麥種植面積的提取。冬小麥播種面積的提取,是估算全國(guó)糧食作物生產(chǎn)能力的重要基礎(chǔ)。TM數(shù)據(jù)適用于對(duì)GIS 所支持的NOAA 數(shù)據(jù)的像素分解以及作物種子的提取數(shù)據(jù)的計(jì)算機(jī)自動(dòng)分類。

        6)冬小麥長(zhǎng)勢(shì)及災(zāi)害監(jiān)測(cè)。檢測(cè)的主要手段,是通過(guò)植被指數(shù)的變動(dòng)以及與歷史數(shù)據(jù)的對(duì)比,檢測(cè)各個(gè)生長(zhǎng)季的植被指數(shù)。能夠準(zhǔn)確掌握各個(gè)生長(zhǎng)期冬小麥的發(fā)育狀況,并能夠根據(jù)發(fā)生狀態(tài)預(yù)測(cè)冬小麥的發(fā)育變化趨勢(shì)。

        7)建立冬小麥遙感估產(chǎn)模型。建立估產(chǎn)模型是冬小麥產(chǎn)量估算的關(guān)鍵。

        8)估算冬小麥總產(chǎn)并對(duì)其精度進(jìn)行評(píng)估。通過(guò)遙感估算綜合技術(shù)估計(jì)冬小麥,由于"精度"可以代表了估算成果的可靠性,為保證結(jié)果的準(zhǔn)確度水平,應(yīng)該把各個(gè)環(huán)節(jié)錯(cuò)誤的概率減至最小化。

        5 冬小麥遙感估產(chǎn)的方法

        在農(nóng)業(yè)實(shí)踐中,植被指數(shù)一般作為判斷植物生長(zhǎng)發(fā)育情況的重要指標(biāo)。一個(gè)最基本的方法就是把植被指數(shù)特性與作物生產(chǎn)特性聯(lián)系起來(lái),方法最好使用回歸方程,而模型試驗(yàn)指標(biāo)則可以采用模型擬合系數(shù)。

        目前,遙感產(chǎn)量估算的方式主要包括統(tǒng)計(jì)建模(經(jīng)驗(yàn)?zāi)P停┓?、機(jī)理模型法、半經(jīng)驗(yàn)(半機(jī)理)模擬方法等。1)統(tǒng)計(jì)建模(經(jīng)驗(yàn)?zāi)P停┓ㄍㄟ^(guò)選擇能夠反映冬小麥生長(zhǎng)發(fā)育階段的遙感技術(shù)特征的數(shù)據(jù),如遙感植地系數(shù)、葉面積系數(shù)、植被的凈初級(jí)生產(chǎn)力等,從而構(gòu)建出與冬小麥生產(chǎn)有關(guān)的數(shù)學(xué)模型,由于具有簡(jiǎn)單、統(tǒng)計(jì)簡(jiǎn)便的優(yōu)點(diǎn),因此比較普遍地運(yùn)用在冬小麥遙感估產(chǎn)。2)機(jī)理模型法根據(jù)生物相似性原理,研究冬小麥生長(zhǎng)與發(fā)育的物理機(jī)制、物理環(huán)境與自然因素等,以物理學(xué)規(guī)律的數(shù)學(xué)模型研究冬小麥生長(zhǎng)與發(fā)育、產(chǎn)量形成的原理,并在某種假設(shè)情況下,形成對(duì)冬小麥估產(chǎn)的數(shù)值模擬模型,最常見(jiàn)的作物生長(zhǎng)模式為DSSAT、WOFOST、CERES-Wheat 等[11-14]。機(jī)理模型的應(yīng)用領(lǐng)域廣泛,但參數(shù)修正難度較大,在特定狀態(tài)下有些參數(shù)不能修正。3)半經(jīng)驗(yàn)(半機(jī)理)模擬方法。半經(jīng)驗(yàn)(半機(jī)理)模型相較于機(jī)理模型有所簡(jiǎn)化,但是模型輸入中需要的2 個(gè)關(guān)鍵變量(光能利用效率、作物收獲指數(shù)),都難以在時(shí)空分布中進(jìn)行準(zhǔn)確地定量模擬。

        統(tǒng)計(jì)模型法因?yàn)椴僮骱?jiǎn)便,已成為迅速開(kāi)展冬小麥遙感估產(chǎn)的重要常規(guī)方式,并獲得了廣泛應(yīng)用。它主要通過(guò)建立關(guān)于植被指數(shù)和冬小麥產(chǎn)量的線性或非線性回歸方程,然后經(jīng)過(guò)對(duì)比分析它們的相互關(guān)系,并選擇影響最大的植被指數(shù)和冬小麥產(chǎn)量的回歸方程,成為冬小麥生產(chǎn)預(yù)測(cè)模型??紤]到不同植物植被指數(shù)的共同作用,產(chǎn)量估算的準(zhǔn)確性相對(duì)較低。近年來(lái),機(jī)器學(xué)習(xí)算法被應(yīng)用于建立作物性狀的遙感估測(cè)模式。其中,BP神經(jīng)網(wǎng)絡(luò)、支持向量機(jī)和隨機(jī)森林算法能基于單一時(shí)相進(jìn)行遙感分類并達(dá)到較高的準(zhǔn)確率,有效解決了衛(wèi)星影像有限的問(wèn)題,同時(shí)打破了傳統(tǒng)方法中僅依靠波段反射率來(lái)實(shí)現(xiàn)信息提取的局限性,實(shí)現(xiàn)了冬小麥遙感識(shí)別最佳特征的優(yōu)化。

        6 冬小麥產(chǎn)量估算模型

        國(guó)內(nèi)外研究者在作物估產(chǎn)模型上已有深入的研究,針對(duì)冬小麥種植的區(qū)域環(huán)境差異研發(fā)出了不同的估產(chǎn)模型,具體可以分為直接和間接2 種類型(見(jiàn)表1)。江東等給出了基于人工神經(jīng)網(wǎng)絡(luò)的農(nóng)作物遙感估產(chǎn)模型[15]。朱再春等采用信息擴(kuò)散方法,結(jié)合農(nóng)業(yè)關(guān)鍵時(shí)期的遙感數(shù)據(jù)構(gòu)建的冬小麥估產(chǎn)模式,更好地模擬了冬小麥在遙感估產(chǎn)中歸一化植被指數(shù)與產(chǎn)量之間的非線性關(guān)系[16]。程千等使用偏小二乘回歸法、支持向量機(jī)回歸法和隨機(jī)森林回歸3 個(gè)機(jī)器學(xué)習(xí)算法在估算冬小麥產(chǎn)量中,總結(jié)出RFR 模式有著良好的穩(wěn)定性[17]。陶惠林等采用多元線性回歸的MLR 模型、隨機(jī)森林RF 模型在植物不同生長(zhǎng)期的12 種光譜參數(shù)的產(chǎn)量估測(cè),得出MLR 模型預(yù)測(cè)準(zhǔn)確度更高,如挑旗期、開(kāi)花期、灌漿期等,與產(chǎn)量相關(guān)最好的分別為TCARI/OSAVI、SR、TCARI/OSAVI。就生長(zhǎng)期而言,效果最佳的是灌漿期[18]。林少喆等通過(guò)對(duì)比以植物含氮率、植物含水率、水氮耦合為中間參量構(gòu)建估產(chǎn)模型,發(fā)現(xiàn)采用水氮耦合的高光譜估產(chǎn)模型準(zhǔn)確率最高,模型最為穩(wěn)定,在灌漿期時(shí)效果亦為最好[19]。而張松等比較直接與間接模型的精度,采用LAI、SPAD,地上干生物量作為中間量進(jìn)行估產(chǎn),對(duì)比表明采用LAI 構(gòu)建的間接模型效果較好,而估產(chǎn)的最佳時(shí)間則是開(kāi)花期[20]。

        表1 冬小麥產(chǎn)量估測(cè)模型表

        7 討論與展望

        研究者們對(duì)農(nóng)作物估產(chǎn)的研究成果不盡相同,今后要分生育期地考慮冬小麥產(chǎn)量預(yù)測(cè),以更精確地掌握產(chǎn)量預(yù)測(cè)生育階段。另外,有些方面需要更深層次地研究,但由于冬小麥種類、生長(zhǎng)期、生產(chǎn)地域的差異等原因,估產(chǎn)監(jiān)測(cè)模型仍無(wú)法大范圍推進(jìn),因此必須盡量完善模型適用性,并研究提升模型精度的途徑,分生育時(shí)期進(jìn)行生長(zhǎng)指數(shù)的監(jiān)測(cè)和估產(chǎn)模型的構(gòu)建。在反演小麥產(chǎn)量方法中,農(nóng)學(xué)參數(shù)也可借助一些生理生化指標(biāo)來(lái)進(jìn)行估產(chǎn)。

        遙感技術(shù)在現(xiàn)代農(nóng)業(yè)中的應(yīng)用已獲得初步成功,而更進(jìn)一步深層次的應(yīng)用仍擁有著巨大的發(fā)展前景。因此,需要根據(jù)我國(guó)國(guó)內(nèi)現(xiàn)階段和未來(lái)現(xiàn)代農(nóng)業(yè)技術(shù)的發(fā)展趨勢(shì),充分地開(kāi)展農(nóng)業(yè)長(zhǎng)遠(yuǎn)規(guī)劃研究和加強(qiáng)現(xiàn)代農(nóng)產(chǎn)品應(yīng)用同傳統(tǒng)遙感技術(shù)之間的聯(lián)系,從而利用現(xiàn)代遙感技術(shù)的巨大空間和農(nóng)業(yè)科技資源優(yōu)勢(shì),對(duì)中國(guó)今后現(xiàn)代農(nóng)業(yè)的發(fā)展將起到很大的促進(jìn)作用。

        猜你喜歡
        估產(chǎn)植被指數(shù)遙感技術(shù)
        基于無(wú)人機(jī)多光譜遙感數(shù)據(jù)的煙草植被指數(shù)估產(chǎn)模型研究
        遙感技術(shù)在大豆種植情況監(jiān)測(cè)中的應(yīng)用
        AMSR_2微波植被指數(shù)在黃河流域的適用性對(duì)比與分析
        河南省冬小麥產(chǎn)量遙感監(jiān)測(cè)精度比較研究
        遙感技術(shù)在林業(yè)中的應(yīng)用
        基于地級(jí)市的區(qū)域水稻遙感估產(chǎn)與空間化研究
        主要植被指數(shù)在生態(tài)環(huán)評(píng)中的作用
        西藏科技(2015年1期)2015-09-26 12:09:29
        遙感技術(shù)在農(nóng)村土地承包經(jīng)營(yíng)權(quán)確權(quán)登記中的應(yīng)用
        河北遙感(2015年2期)2015-07-18 11:11:14
        基于MODIS數(shù)據(jù)的植被指數(shù)與植被覆蓋度關(guān)系研究
        基于SAR技術(shù)的高原山區(qū)煙草估產(chǎn)模型
        久久久久人妻一区精品| 黄色国产一区在线观看| 亚洲国产成人久久精品美女av| 日本伦理视频一区二区| 日本人妻伦理在线播放| 99久久精品免费看国产| 好屌草这里只有精品| а中文在线天堂| 国产成人精品亚洲午夜| 无码熟妇人妻av在线c0930| 白色白色在线视频播放平台| 日本按摩偷拍在线观看| 久久久99精品成人片| 国产成+人+综合+亚洲欧美丁香花| 97免费人妻在线视频| 精品无码国产一二三区麻豆| 精品国模人妻视频网站| 91色老久久偷偷精品蜜臀懂色 | 男女做爰高清免费视频网站| 婷婷综合缴情亚洲| 日韩精品国产自在欧美| 国产丝袜美腿诱惑在线观看| 久久精品不卡一区二区三区| 久久成人国产精品免费软件| 亚洲精品久久久久久| 无码伊人久久大蕉中文无码| 中文字幕一区二区三区乱码人妻| 吃奶呻吟打开双腿做受视频 | 亚洲成av人片天堂网无码| 色偷偷av亚洲男人的天堂| 亚洲av美女在线播放啊| 亚洲精品日本久久久中文字幕 | 国产婷婷色综合av蜜臀av| 九九久久精品国产| 国产欧美日韩专区毛茸茸| 日本视频一区二区这里只有精品| aaa日本高清在线播放免费观看| 国产成人精品日本亚洲| 国产人澡人澡澡澡人碰视频| 最好的99精品色视频大全在线| 久久人妻av一区二区软件|