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

    基于哨兵二號數(shù)據(jù)的礦區(qū)葉面降塵信息提取及塵源識別

    2023-10-19 02:48:10帥爽張志呂新彪陳思馬梓程謝翠容
    遙感學(xué)報 2023年9期
    關(guān)鍵詞:降塵葉面反射率

    帥爽,張志,呂新彪,陳思,馬梓程,謝翠容

    1.中國地質(zhì)大學(xué)(武漢)地質(zhì)調(diào)查研究院,武漢 430074;2.中國地質(zhì)大學(xué)(武漢)地球物理與空間信息學(xué)院,武漢 430074;3.湖北省國土測繪院,武漢 430010

    1 引言

    灰塵是氣溶膠系統(tǒng)的主要成分(Kaufman 等,2001),會影響空氣質(zhì)量、大氣化學(xué)成分和氣候狀況(Fairlie 等,2007)。葉面降塵指大氣中的浮塵經(jīng)重力沉降后,在植物葉片表面所形成的積塵(彭杰 等,2013)。植被通過葉面降塵現(xiàn)象固定空氣中的灰塵,達(dá)到凈化空氣的效果(Zhu 等,2020),但過量的葉面降塵會阻礙植物生長,降低生產(chǎn)力,改變植物的群落結(jié)構(gòu),還會導(dǎo)致葉片壞死(Kayet 等,2019;Vellak 等,2014)。葉面降塵的塵源可分為自然源和人為源兩大類(Tegen 等,2004),自然源包括沙塵暴、火山爆發(fā)、海洋飛沫等,人為源包括工程建設(shè)、交通運輸、礦產(chǎn)資源開發(fā)等(Ma等,2020)。其中礦山粉塵富集大量重金屬的顆粒物(Allen等,2001),對礦山周邊居民健康和植被生態(tài)造成嚴(yán)重威脅。因此在礦山及周邊區(qū)域監(jiān)測葉面降塵狀況能有效反應(yīng)礦山粉塵污染狀況。本研究提出了一種利用衛(wèi)星遙感數(shù)據(jù)快速獲取礦山周邊區(qū)域葉面降塵范圍、強(qiáng)度,同時區(qū)分自然降塵和礦山降塵的方法。

    葉面降塵的光譜響應(yīng)特征是利用光學(xué)遙感數(shù)據(jù)監(jiān)測葉面降塵量的基礎(chǔ)。前人對不同區(qū)域、不同塵源、不同植被類型的葉面降塵光譜特征進(jìn)行了定量研究。Chudnovsky 和Ben-Dor(2008)利用可見光—短波紅外波段反射率數(shù)據(jù)研究室內(nèi)降塵量和降塵礦物成分,發(fā)現(xiàn)可見光區(qū)間光譜指標(biāo)可用于估算樣品中的礦物和有機(jī)物相對含量。彭杰等(2013,2015)對南疆受沙塵暴影響的榆樹葉片和香梨樹葉片的葉面降塵光譜特征進(jìn)行了研究,認(rèn)為葉面降塵可提高綠色植物葉片可見光波段的反射率,降低近紅外波段的反射率,對三邊位置沒有影響。Yan 等(2015)和Zhu 等(2019)等通過分析北京市大葉黃楊的葉面降塵光譜特征,研究城市的空氣粉塵污染狀況。Yan等(2015)發(fā)現(xiàn)在近紅外波段(780—1300 nm),降塵葉面光譜的反射率明顯低于無塵葉面。Zhu等(2019)發(fā)現(xiàn)對大葉黃楊的葉面降塵量最敏感的光譜區(qū)間為695—1400 nm,降塵葉面光譜紅邊位置明顯藍(lán)移,且“藍(lán)移”距離隨降塵量增加而增大。在此基礎(chǔ)上,Zhu 等(2020)分析了不同葉面降塵成分(煤塵,水泥塵和純土塵)的光譜響應(yīng)特征差異,發(fā)現(xiàn)在760—1400 nm,隨著葉面降塵量增加,煤塵,水泥塵和純土塵葉面反射率的減小幅度依次降低。Kayet 等(2019)、Ma 等(2020)、Ma 等(2017)、Zhao 等(2020)研究鐵礦、煤礦礦區(qū)林地、草地植被葉面降塵的光譜響應(yīng)特征,獲得了相似的結(jié)果,即葉面降塵導(dǎo)致植被可見光區(qū)間反射率升高,近紅外區(qū)間反射率降低。Ma 等(2020)通過模擬礦山粉塵的點源、線源擴(kuò)散特征,利用三維離散各向異性輻射傳輸(DART)模型模擬了植被冠層的滯塵響應(yīng)光譜,同時分析了像元尺度上,植被覆蓋度對葉面滯塵量估算的影響。這些研究定量分析了葉面、冠層、像元不同尺度上降塵量變化的植被光譜響應(yīng)特征,但少有研究分析降塵的礦物成分差異,及其引起的葉面降塵光譜響應(yīng)特征差異。尤其對于礦山區(qū)域,不同礦物組成的礦山粉塵滯留于葉面上可能形成不同的光譜特征。

    在開展葉面降塵的光譜響應(yīng)特征研究的基礎(chǔ)上,前人利用多光譜和高光譜遙感數(shù)據(jù)進(jìn)行了葉面降塵量定量反演。Jing等(2018)利用支持向量機(jī)(SVM)、分類回歸樹(CART)和隨機(jī)森林(RF)等機(jī)器學(xué)習(xí)算法構(gòu)建了葉面沉降量與光譜參數(shù)的關(guān)系模型,結(jié)果表明RF 的預(yù)測模型精度更高(R2=0.879)。Zhao等(2020)利用實測光譜與實測葉面降塵量數(shù)據(jù)建立葉面降塵量反演模型,并應(yīng)用于機(jī)載高光譜數(shù)據(jù)完成區(qū)域尺度葉面降塵量制圖。Ma 等(2017)和Zhao 等(2020)針對礦區(qū)周邊的粉塵污染問題,利用實測光譜與實測葉面降塵量數(shù)據(jù)建立葉面降塵量反演模型,并應(yīng)用于Hyperion 和OLI 數(shù)據(jù),獲取了不同礦山地物周邊的葉面降塵量分布。王皓飛等(2016)利用OLI的綠波段和NDVI 構(gòu)建了降塵量反演模型。Yan 等(2015)基于實測數(shù)據(jù),利用BP 神經(jīng)網(wǎng)絡(luò)方法構(gòu)建降塵量預(yù)測模型,并應(yīng)用于MODIS 數(shù)據(jù),獲取了北京市葉面降塵量狀況。這些研究大多在葉面尺度建立預(yù)測模型,直接應(yīng)用于航空航天遙感數(shù)據(jù),未考慮像元尺度下混合像元效應(yīng)、植被類型、冠層結(jié)構(gòu)和植被覆蓋度差異對預(yù)測模型適應(yīng)性的影響(Sims 和Gamon,2002;Stagakis 等,2010)。另外,基于實測數(shù)據(jù)構(gòu)建經(jīng)驗預(yù)測模型進(jìn)行葉面降塵量遙感監(jiān)測的方法依賴野外采樣測試工作,難以實現(xiàn)大區(qū)域葉面降塵的快速監(jiān)測。同時,構(gòu)建的預(yù)測模型也難以直接應(yīng)用于其他區(qū)域。

    前人還通過構(gòu)建新植被指數(shù)、評估已有植被指數(shù)對葉面降塵的敏感程度來獲取葉面降塵范圍和強(qiáng)度。彭杰等(2013)構(gòu)建了葉面降塵差值指數(shù)FDCDI(R852-R673)、葉面降塵比值指數(shù)FDCRI(R852/R673)以及葉面降塵歸一化指數(shù)FDCNI((R852-R673)/(R852+R673))。Zhang 等(2014)發(fā)現(xiàn)結(jié)構(gòu)不敏感色素指數(shù)SIPI((R800-R445)/(R800-R680))、植被水分指數(shù)PWI((R900/R970)×(R860+R1240)/(R860-R1240))與葉面煤炭粉塵量的相關(guān)性較強(qiáng)。Su 等(2019)認(rèn)為歸一化差異物候指數(shù)(NDPI)(Wang 等,2017)更適宜用于構(gòu)建葉面降塵量回歸模型。然而,這些植被指數(shù)大多是針對高光譜數(shù)據(jù)設(shè)計的,無法應(yīng)用于OLI、哨兵二號等多光譜衛(wèi)星遙感數(shù)據(jù)。

    特征向量主成分選擇方法FPCS(Featureoriented Principal Components Selection)由Crósta 和Moore(1989)提出,廣泛應(yīng)用于鐵氧化物、羥基礦物、碳酸根礦物等遙感礦化蝕變“弱信息”提取(Liu等,2011,2018;Jain等,2018;Moradpour等,2022;Wambo 等,2020)。傳統(tǒng)的主分成分析(PCA)方法使用遙感數(shù)據(jù)的全部波段進(jìn)行主成分變換,礦化蝕變等“弱信息”往往包含于高順位的主分量中,而這些高順位主分量圖像中含有大量噪聲,對目標(biāo)信息提取造成了較大干擾(Crósta和Moore,1989)。FPCS 方法是根據(jù)目標(biāo)地物波譜特征,針對性選擇若干個診斷性波譜特征對應(yīng)的波段進(jìn)行主成分分析(?zyava?,2016),這樣目標(biāo)地物信息就更有可能包含于噪聲信息較少的低順位主分量中(Crósta 和Moore,1989)。葉面降塵信息對于正常植被來說,也是一種“弱”的“異常信息”。此次研究中首次嘗試?yán)肍PCS 方法提取葉面降塵信息。

    此次研究旨在分析評價哨兵二號遙感數(shù)據(jù)對不同來源葉面降塵的光譜響應(yīng)敏感性,同時提出一種快速、可靠的區(qū)域性葉面降塵信息提取方法,為礦區(qū)粉塵污染監(jiān)測提供一種新思路。

    2 研究區(qū)及數(shù)據(jù)源

    2.1 研究區(qū)概況

    研究區(qū)位于內(nèi)蒙古自治區(qū)新巴爾虎右旗阿敦楚魯蘇木的甲烏拉—查干鉛鋅銀礦區(qū)(圖1(a)),礦區(qū)大規(guī)模建設(shè)始于2001 年,2007 年礦山開采系統(tǒng)逐漸完善,投生產(chǎn)運營的礦井逐漸增多(高德榮 等,2016)。礦區(qū)廢石堆分布在礦井周邊,無防塵措施,尾礦庫壩頂多修筑防塵網(wǎng),礦山道路聯(lián)通各礦山地物,道路兩側(cè)無防塵設(shè)施。礦區(qū)植被類型主要為草地。區(qū)內(nèi)葉面降塵空間分布與風(fēng)要素關(guān)系密切(Rodríguez 等,2009),由國家氣象信息中心網(wǎng)站(http://www.nmic.cn/data[2021-01-05])獲取新巴爾虎右旗氣象觀測站近20 年風(fēng)向和風(fēng)速數(shù)據(jù),分析得知區(qū)內(nèi)春秋冬季干旱,多西北風(fēng),每年2—5 月和9—10 月累年月優(yōu)勢風(fēng)向為西北風(fēng),累年月平均風(fēng)速常年在3 m/s 左右,其中4—5 月份達(dá)到4 m/s 以上。區(qū)內(nèi)葉面降塵的塵源主要有兩種,一是來自鉛鋅礦礦山開發(fā)中礦石、廢石的轉(zhuǎn)運和堆放,在晴天風(fēng)速大于2.5 m/s,尾礦庫和礦石堆即可發(fā)生風(fēng)蝕揚塵(?or?evi? 等,2004;Etyemezian,2003)。鉛鋅礦揚塵的主要化學(xué)成分為SiO2,CaO,Al2O3,F(xiàn)e2O3等,顏色呈淡藍(lán)色(宋國策和張志,2020);二是區(qū)內(nèi)沙地和干涸的湖泊,2004 年全國沙漠化普查顯示,呼倫貝爾沙化土地總面積131萬ha,新巴爾虎右旗又是呼倫貝爾市沙地分布的核心區(qū)(趙慧穎,2007),沙地內(nèi)松散物質(zhì)在風(fēng)力作用下傳播形成葉面降塵。Google Earth 影像顯示礦山尾礦庫(圖1(b))和干涸湖泊(圖1(c))優(yōu)勢風(fēng)向的下風(fēng)向草地葉面滯塵現(xiàn)象十分明顯。

    圖1 研究區(qū)、訓(xùn)練區(qū)位置與礦山、干涸湖泊揚塵影像特征(干涸湖泊揚塵區(qū)指,在枯水期,湖泊干涸后,湖底細(xì)粒沉積物露出地表在風(fēng)力作用下形成的揚塵區(qū)域)Fig.1 The location of the study area and training areas,and the image of dust from mines and dry lake

    2.2 遙感數(shù)據(jù)及處理

    哨兵二號遙感數(shù)據(jù)由歐州航天局(https://scihub.copernicus.eu/[2020-11-20])提供。由于區(qū)內(nèi)草地生長期為5—9 月,此次研究中收集了研究區(qū)2019 年9 月18 日采集的L1C 級哨兵二號A 星遙感數(shù)據(jù)和2019 年8 月22 日采集的高分一號遙感數(shù)據(jù),數(shù)據(jù)云覆蓋量均小于5%。哨兵二號A 星(S2A)波段設(shè)置及參數(shù)如表1。利用Sen2cor 軟件對S2A 數(shù)據(jù)進(jìn)行輻射定標(biāo)和大氣校正,排除水汽(B9)和卷云(B10)波段,獲得11 波段(B1-8A,B11-12)空間分辨率為10 m 的地表反射率數(shù)據(jù)。同時,對高分一號遙感數(shù)據(jù)進(jìn)行大氣校正、幾何校正、圖像融合,獲得4波段(藍(lán)、綠、紅、近紅外)空間分辨率為2 m的正射影像。

    表1 哨兵2號A星(S2A)波段設(shè)置及參數(shù)Table 1 Band settings and parameters of Sentinel-2A

    3 研究方法

    此次工作在選擇合適訓(xùn)練區(qū),在分析葉面降塵植被與正常植被的光譜特征差異以及不同塵源的葉面降塵光譜特征差異基礎(chǔ)上,利用FPCS 方法提取葉面降塵范圍和強(qiáng)度;構(gòu)建塵源光譜指數(shù)DSI(Dust-source Spectrum Index),區(qū)分提取的葉面降塵信息的塵源類別;結(jié)合利用高分遙感數(shù)據(jù)提取礦山地物信息,分析葉面降塵信息與礦山地物的空間相關(guān)性;針對不同塵源葉面降塵訓(xùn)練區(qū),利用MNF(Minimum Noise Fraction)、PPI(Pixel Purity Index)等方法,提取、對比不同塵源純凈像元光譜差異,并分析不同塵源的礦物組成差異。研究方法和流程如圖2。

    圖2 研究方法與技術(shù)流程Fig.2 Methods and technical flow of this research

    3.1 光譜特征分析

    結(jié)合野外踏勘和哨兵二號遙感影像特征,分別選取礦山揚塵區(qū)、干涸湖泊揚塵區(qū)和遠(yuǎn)離塵源的正常植被區(qū)3 種不同葉面降塵狀況的訓(xùn)練區(qū)(圖1(a)),其中礦山揚塵和干涸湖泊揚塵訓(xùn)練區(qū)分別選取在礦山廢石堆和干涸湖泊的下風(fēng)向揚塵影像特征明顯的區(qū)域,正常植被訓(xùn)練區(qū)選取在遠(yuǎn)離礦山和自然塵源,植被光譜特征典型的區(qū)域。計算訓(xùn)練區(qū)內(nèi)像元的反射率均值光譜(圖3),對比S2A 反射率數(shù)據(jù)對不同葉面降塵狀況的響應(yīng)特征。結(jié)果顯示,在可見光波段(B1—B4),葉面降塵像元(礦山、干涸湖泊揚塵區(qū))反射率大于正常植被像元,而對于近紅外波段(B5—B8A),葉面降塵像元反射率小于正常植被。葉面降塵像元綠波段(B3,559.8 nm 附近)至紅波段(B4,664.6 nm 附近)斜率大于正常草地像元。紅邊特征方面,葉面降塵像元紅邊斜率小于正常植被,紅邊位置也稍向短波方向移動。礦山揚塵區(qū)像元均值光譜與干涸湖泊揚塵區(qū)像元均值光譜相比,可見光波段和短波紅外(B11、B12)波段反射率值更大,近紅外波段兩者反射率接近。另外,礦山揚塵區(qū)像元在8A 波段(864.7 nm 附近)顯示出微弱吸收特征。

    圖3 訓(xùn)練區(qū)像元均值反射率光譜曲線Fig.3 Mean reflectance of training areas pixels

    前人研究表明,葉面降塵量與距塵源的距離成反比(Ma 等,2020)。同時選取影像上葉面降塵特征明顯的典型廢石堆揚塵區(qū)(圖1(e))和干涸湖泊揚塵區(qū)(圖1(d)),按照距離塵源由近及遠(yuǎn)、等間距(40 m)選取光譜采樣點,對比各光譜采樣點反射率變化情況,分析對于不同塵源,S2A 反射率數(shù)據(jù)對葉面降塵量變化的響應(yīng)特征,如圖4。結(jié)果顯示隨著遠(yuǎn)離塵源,礦山廢石堆揚塵區(qū)與干涸湖泊揚塵區(qū)像元光譜特征的變化顯示出不同規(guī)律。對于廢石堆揚塵區(qū),隨著遠(yuǎn)離塵源,像元反射率曲線形態(tài)變化不明顯,各波段反射率值整體降低,紅邊位置逐漸向長波方向移動。另外,大多數(shù)采樣像元在8A 波段(864.7 nm附近)顯示出吸收特征。對于干涸湖泊揚塵區(qū),隨著遠(yuǎn)離塵源,像元反射率曲線形態(tài)有一定變化,紅波段(B4,664.6 nm 附近)處吸收特征逐漸顯現(xiàn),在可見光波段,反射率逐漸降低,而在近紅外波段,反射率逐漸升高,紅邊位置逐漸向長波方向移動。

    圖4 距塵源不同距離采樣點光譜曲線對比Fig.4 Comparison of sampling points reflectance curves at different distances to the dust source

    3.2 特征向量主成分選擇方法(FPCS)

    在分析S2A 反射率對葉面降塵信息的光譜響應(yīng)特征的基礎(chǔ)上,利用特征向量主成分選擇(FPCS)方法提取研究區(qū)葉面降塵范圍和強(qiáng)度。本文將葉面降塵信息視作對于正常植被的是一種“弱”的“異常信息”,按照FPCS 方法,選擇若干個區(qū)分降塵葉面與正常葉面診斷性波譜對應(yīng)的波段進(jìn)行主成分分析,獲取特征向量載荷因子矩陣(吳志春 等,2018)。特征向量載荷因子表達(dá)了各主分量與原始波段的相關(guān)性,載荷因子的正、負(fù)符號指示主分量與對應(yīng)波段的正、負(fù)相關(guān)性,載荷因子數(shù)值的大小則表示主分量與對應(yīng)波段相關(guān)性大?。▍侵敬?等,2018)。根據(jù)葉面降塵像元光譜吸收、反射的特征波段位置,選擇與目標(biāo)地物光譜吸收、反射特征最相關(guān)的主分量進(jìn)行信息提取,并通過設(shè)置適當(dāng)閾值獲取研究區(qū)葉面降塵的分布和強(qiáng)度。

    3.3 葉面降塵的塵源光譜指數(shù)

    在提取研究區(qū)葉面降塵范圍的基礎(chǔ)上,增強(qiáng)自然降塵像元與礦山降塵像元的圖像反射率差異,以達(dá)到識別并區(qū)分葉面降塵像元的塵源的目的,此次研究利用礦山塵源葉面降塵像元與自然塵源葉面降塵像元的光譜差異構(gòu)建光譜指數(shù),以區(qū)分礦山塵源葉面降塵與自然塵源葉面降塵的范圍。訓(xùn)練區(qū)像元均值光譜(圖3)和采樣點像元光譜(圖4)分析結(jié)果顯示:相比自然塵源葉面降塵像元,礦山塵源像元在864.7 nm(B8A)附近呈現(xiàn)特殊的吸收特征。依據(jù)此特征,選取S2A 數(shù)據(jù)B8、B8A和B11波段構(gòu)建了適宜本研究區(qū)的葉面降塵的塵源光譜指數(shù)DSI(Dust-source Spectrum Index):

    式中,R代表B8、B8A和B11波段反射率,由于相對自然塵源像元,礦山塵源像元反射率在B8A 波段顯示吸收特征,B8 和B11 波段為該吸收谷兩側(cè)的反射峰,所以R(B8)與R(B11)相加,再除以R(B8A)可以增大礦山塵源像元的值,以及礦山塵源與自然塵源像元的差異。

    3.4 礦山地物信息提取與空間相關(guān)性分析

    利用高分一號遙感數(shù)據(jù)對研究區(qū)礦山地物進(jìn)行目視解譯,重點解譯尾礦庫、廢石堆、礦山道路等礦山塵源的分布狀況。在利用哨兵二號數(shù)據(jù)和FPCS 方法提取研究區(qū)葉面降塵范圍和強(qiáng)度,利用塵源光譜指數(shù)(DSI)區(qū)分自然塵源和礦山塵源的基礎(chǔ)上,按照葉面降塵強(qiáng)度類型對研究區(qū)像元進(jìn)行分類賦值(高降塵強(qiáng)度、中降塵強(qiáng)度、低降塵強(qiáng)度和無降塵像元分別賦值為3,2,1,0),并計算各類礦山塵源外圍等距離范圍內(nèi)(0—100 m,100—200 m,200—300 m,…,900—1000 m)像元的降塵強(qiáng)度分類賦值的均值變化情況,分析提取的葉面降塵信息與礦山地物分布的空間相關(guān)性,同時分析各類型礦山地物周邊揚塵擴(kuò)散的空間特征。

    3.5 塵源純凈像元提取與光譜特征匹配

    為了探討礦山塵源與自然塵源葉面降塵光譜特征差異的原因,提取并對比研究區(qū)S2A 數(shù)據(jù)上典型廢石堆揚塵區(qū)(圖1(e))和干涸湖泊揚塵區(qū)(圖1(d))塵源純凈像元光譜,分析研究區(qū)礦山塵源與自然塵源的光譜與礦物成分差異,具體步驟如下:

    (1)正交化處理。利用最小噪聲分離(MNF)(Green 等,1988)變換方法對研究區(qū)哨兵二號影像數(shù)據(jù)進(jìn)行正交化處理,去除各波段間冗余信息。

    (2)光譜端元提取。選擇MNF變換后的前4 個主分量,使用純凈像元指數(shù)PPI(Pixel Purity Index)方法對圖像中的像元進(jìn)行迭代運算,找出圖像中相對純凈的像元。利用ENVI 5.3 平臺的N維可視化功能(The n-D Visualizer),結(jié)合影像特征,分別獲取礦山塵源和干涸湖泊塵源的端元光譜。

    4 實驗結(jié)果

    4.1 葉面降塵信息提取

    此次葉面降塵信息提取的思路是,利用FPCS方法,增強(qiáng)葉面降塵像元與正常植被像元的差異。首先計算研究區(qū)歸一化植被指數(shù),并通過試錯法,確定0.32 的閾值排除植被以外的其他地物。由于葉面降塵植被光譜與正常植被光譜的差異主要體現(xiàn)在可見光波段和近紅外波段,選擇S2A 數(shù)據(jù)B1(442.7 nm附近)、B2(492.4 nm附近)、B3(559.8 nm附近)、B4(664.6 nm 附近)、B8(832.8 nm 附近)和B8A(864.7 nm附近)波段進(jìn)行主成分分析,獲得了特征向量載荷因子矩陣(表2)。

    表2 特征向量載荷因子矩陣Table 2 Load factor matrix of eigenvector

    葉面降塵植被像元與正常植被像元反射率特征差異主要表現(xiàn)為:可見光波段,葉面降塵像元反射率大于正常植被像元;近紅外波段,葉面降塵像元反射率小于正常植被像元。PC2與可見光波段特征向量載荷因子均為正,與近紅外波段特征向量載荷因子均為負(fù),表明PC2與可見光波段(B1、B2、B3 和B4)呈正相關(guān),與近紅外波段(B8 和B8A)呈負(fù)相關(guān)。同時,PC2對應(yīng)B1、B2、B8 和B8A 的特征向量載荷因子值較高,表明PC2與這些波段相關(guān)性較高。所以在PC2圖像上,葉面降塵植被像元應(yīng)表現(xiàn)為高值,而正常植被像元應(yīng)表現(xiàn)為低值。選擇PC2圖像設(shè)置適當(dāng)閾值提取葉面降塵范圍、劃分降塵強(qiáng)度。閾值的選取使用礦物蝕變信息提取中劃分蝕變異常等級的方法(張玉君 等,2003),利用均值(X)和標(biāo)準(zhǔn)差(σ)的線性組合(X+kσ)確定葉面降塵下限和劃分降塵強(qiáng)度等級。此次研究中以影像上植被葉面降塵特征明顯的區(qū)域為訓(xùn)練區(qū)(圖1(d)、(e)區(qū)域),使用試錯法確定系數(shù)k的取值,最終確定葉面降塵下限(低降塵區(qū)下限)為X+1.5σ,中降塵區(qū)下限為X+2.5σ,高降塵區(qū)下限為X+3.5σ,葉面降塵提取效果如圖5(b)。

    圖5 不同方法葉面降塵信息提取效果對比Fig.5 Comparison of foliar dust information extraction results with different methods

    為對比本文方法的效果,利用(Ma等,2017)等構(gòu)建的葉面降塵量反演模型(-286.5·NDVI+244.9,R2=0.999),同樣利用試錯法確定葉面降塵下限和劃分降塵強(qiáng)度等級,最終確定葉面降塵下限(低降塵區(qū)下限)為X+1.5σ,中降塵區(qū)下限為X+2σ,高降塵區(qū)下限為X+3σ,葉面降塵提取效果如圖5(a)。

    4.2 塵源的識別

    在提取研究區(qū)葉面降塵范圍的基礎(chǔ)上,利用葉面降塵的塵源光譜指數(shù)(DSI)區(qū)分礦山塵源葉面降塵與自然塵源葉面降塵的范圍。利用葉面降塵提取結(jié)果對研究區(qū)S2A 反射率數(shù)據(jù)進(jìn)行掩膜處理后計算DSI 指數(shù),以X+1.5σ為閾值,獲取了研究區(qū)礦山塵源葉面降塵與自然塵源葉面降塵的分布(圖6)。提取的礦山塵源葉面降塵像元主要分布在廢石堆(圖6(a))、尾礦庫、礦山道路等礦山地物周邊,而自然塵源葉面像元主要分布于干涸湖泊(圖6(b))、沙地等自然塵源周邊。

    圖6 葉面降塵塵源分類結(jié)果Fig.6 Classification results of foliar dust source

    4.3 葉面降塵與礦山地物的空間相關(guān)性

    研究區(qū)內(nèi)礦山地物(礦山道路、尾礦庫、廢石堆等)與自然塵源、礦山塵源葉面降塵分布及強(qiáng)度如圖7(a)。礦山地物周邊0—100 m,100—200 m,200—300 m,…,900—1000 m 等距離范圍內(nèi)的像元自然塵源、礦山塵源降塵強(qiáng)度分類賦值均值變化情況如圖7(b)。結(jié)果顯示提取的礦山塵源葉面降塵像元與礦山地物的空間相關(guān)性更強(qiáng),礦山地物周邊700 m 范圍內(nèi),各等距離范圍內(nèi)礦山塵源像元強(qiáng)度賦值均值均大于自然塵源分類像元,并隨遠(yuǎn)離礦山塵源規(guī)律性減小,表明提取的礦山塵源葉面降塵像元在靠近礦山地物方向呈逐漸集中分布趨勢,而提取的自然塵源葉面降塵像元在礦山地物周邊各等距離范圍內(nèi)分布均勻且強(qiáng)度較弱,與礦山地物的空間相關(guān)性較弱。這驗證了本文葉面降塵的塵源區(qū)分方法的可靠性。另外,礦山塵源分類像元降塵強(qiáng)度賦值的均值隨距離礦山地物距離變化的趨勢還顯示出礦山揚塵的擴(kuò)散特征。圖7(b)中顯示,距離大于800 m,礦山塵源像元強(qiáng)度賦值的均值極?。ㄐ∮?.05)并處于穩(wěn)定狀態(tài),表明研究區(qū)礦山塵源揚塵擴(kuò)散的平均最遠(yuǎn)距離為700—800 m。這與前人在類似區(qū)域的實地調(diào)查結(jié)果類似,高原等(2016)等發(fā)現(xiàn)白音諾爾鉛鋅礦(平均風(fēng)速4.5 m/s)區(qū)尾礦庫周邊揚塵風(fēng)積物隨距離變化規(guī)律性變化的區(qū)間為0—900 m。

    圖7 葉面降塵信息與礦山地物的空間相關(guān)性Fig.7 The spatial correlation between foliar dust information and mine objects

    同樣計算了廢石堆、礦山道路、尾礦庫周邊等距離范圍內(nèi)的提取礦山塵源降塵像元降塵強(qiáng)度分類賦值均值變化情況(圖7(c))。結(jié)果顯示提取的礦山塵源葉面降塵像元與廢石堆和礦山道路的空間相關(guān)性更強(qiáng),與尾礦庫的空間相關(guān)性較弱,表明研究區(qū)內(nèi)礦山塵源主要為廢石堆和礦山道路。同時,廢石堆周邊各等距離范圍內(nèi)像元降塵強(qiáng)度賦值均值均大于礦山道路,表明廢石堆周邊葉面降塵分布更集中、強(qiáng)度更高。另外,廢石堆的揚塵擴(kuò)散距離大于礦山道路。圖7(c)顯示距離礦山道路大于700 m,礦山塵源像元降塵強(qiáng)度賦值均值極?。ㄐ∮?.05)并處于穩(wěn)定狀態(tài),而距離廢石堆700—1000 m,礦山塵源像元均值仍處于變化中,且明顯大于同距離等級的礦山道路的像元降塵強(qiáng)度賦值均值。礦山塵源的擴(kuò)散特征還與是否有保護(hù)措施有關(guān),高分影像上顯示研究區(qū)幾處尾礦庫邊緣均修建了防護(hù)網(wǎng),有效抑制了揚塵擴(kuò)散。

    5 討論

    5.1 塵源礦物成分的光譜差異分析

    影像上典型礦山廢石堆揚塵區(qū)(圖1(e))和干涸湖泊揚塵區(qū)(圖1(d))的塵源和揚塵擴(kuò)散影像特征明顯。利用MNF、PPI等方法提取了塵源近似純凈像元,兩區(qū)域塵源純凈像元分布如圖8(a)(b),可見典型尾礦滯塵區(qū)塵源純凈像元主要分布于廢石堆中心(圖8(a)),湖積物滯塵區(qū)塵源純凈像元主要分布于干涸湖泊東南側(cè)的揚塵中心(圖8(b))。

    圖8 礦山塵源與干涸湖泊塵源純凈像元分布及反射率均值曲線對比Fig.8 Pure pixels distribution of mine and dry lake dust sources and comparison of mean reflectance curves of pure pixels between mine and lake sediment dust source

    分別提取礦山塵源和干涸湖泊塵源純凈像元的均值光譜(圖8(c)),可見不同塵源光譜特征差異主要體現(xiàn)在可見光至近紅外波段。塵源光譜差異主要受塵源礦物成分影響。礦山塵源純凈像元均值光譜在864.7 nm(B8A)附近顯示明顯吸收特征,同時在可見光區(qū)間表現(xiàn)出強(qiáng)反射特征,這與礦山塵源的降塵草地像元均值光譜表現(xiàn)出相似的吸收特征。前人研究表明864.7 nm 附近的反射率吸收特征與三價鐵的氧化物相關(guān)(張玉君 等,2002;Hunt,1977)。區(qū)域礦床地質(zhì)資料顯示甲烏拉—查干鉛鋅銀礦區(qū)鉛鋅銀礦化的各個階段均有黃鐵礦伴生(高德榮 等,2016),而純凈、無風(fēng)化的黃鐵礦(FeS2)在可見光至短波紅外波段無明顯吸收特征,864.7 nm附近顯示出吸收特征,表明黃鐵礦(FeS2)可能部分氧化形成赤鐵礦(Fe2O3)。前人對鉛鋅礦揚塵的礦物成分分析結(jié)果(鉛鋅礦揚塵的主要化學(xué)成分為SiO2,CaO,Al2O3,F(xiàn)e2O3等)也證實了這一點(?or?evi? 等,2004)。吉磊等(1994)測試了工作區(qū)以東呼倫湖表層沉積物的礦物組成,主要為蒙脫石、綠泥石和伊利石等粘土礦物,證明了該區(qū)域湖積物包含蒙脫石等粘土礦物。蒙脫石等含鋁羥基粘土礦物的反射率特征吸收譜段為2200 nm 附近,由于哨兵二號遙感數(shù)據(jù)在短波紅外波段設(shè)置波段較少,沒能顯示出相應(yīng)的吸收特征。下一步工作中將嘗試?yán)肁STER等在短波紅外區(qū)間有更多波段設(shè)置的數(shù)據(jù)研究干涸湖泊塵源的光譜特征。

    5.2 哨兵二號反射率對葉面降塵的敏感性

    此次研究確認(rèn)了哨兵二號反射率數(shù)據(jù)對于植被葉面降塵的光譜響應(yīng)特征,和區(qū)分葉面降塵像元與非降塵像元、塵源類型的潛力。通過選取訓(xùn)練區(qū),對比了葉面降塵植被與正常植被的光譜差異,發(fā)現(xiàn)可見光波段,降塵植被像元反射率大于正常植被,而對于近紅外波段,降塵植被像元反射率小于正常植被,這與基于地面實測光譜(Ma等,2017;彭杰 等,2015;Su 等,2019;Xu 和Yu,2013;Yan 等,2015;Zhang 等,2014;Zhu等,2020)和機(jī)載高光譜數(shù)據(jù)(Zhao 等,2020)的研究結(jié)論一致??梢姽獠ǘ畏瓷渎实纳吲c葉面降塵對葉綠素光合作用的抑制作用相關(guān),近紅外波段反射率的降低與葉面降塵導(dǎo)致植物葉片組織結(jié)構(gòu)的變化相關(guān)(如柵欄組織細(xì)胞收縮、變短,甚至彎曲)。研究中還發(fā)現(xiàn),降塵植被像元紅邊斜率較正常草地像元減小,紅邊位置藍(lán)移,這與Zhu等(2019,2020)的結(jié)論一致,而與王濤等(2012)、彭杰等(2013)認(rèn)為葉面降塵對植被三邊位置沒有明顯影響的觀點不一致,這可能與葉面降塵量累計對葉片組織結(jié)構(gòu)的影響程度差異相關(guān)。

    還確認(rèn)了哨兵二號數(shù)據(jù)礦山塵源與自然塵源葉面降塵光譜響應(yīng)特征差異,發(fā)現(xiàn)礦山塵源周邊的葉面降塵像元大多在864.7 nm 附近顯示出吸收特征(圖3、圖4);利用MNF 和PPI方法提取的礦山塵源純凈像元也在864.7 nm 附近顯示吸收特征(圖8),經(jīng)過礦床地質(zhì)資料確認(rèn),研究區(qū)礦山塵源中可能含有一定量黃鐵礦氧化形成的鐵氧化物,進(jìn)而使礦山廢石堆光譜、礦山粉塵與草地的混合光譜顯示出864.7 nm 附近吸收特征。逐漸遠(yuǎn)離自然塵源和礦山塵源的等間距采樣點反射率顯現(xiàn)出不同的光譜響應(yīng)特征。兩種塵源的采樣點反射率在可見光區(qū)間逐漸降低,紅藍(lán)波段的植被吸收特征逐漸顯現(xiàn),紅邊位置向長波方向移動,自然塵源采樣點近紅外波段反射率逐漸升高(圖4),這與前人基于地面實測光譜的研究結(jié)果一致(Kayet等,2019;Ma 等,2017;Zhu 等,2019)。而礦山塵源的采樣點反射率隨著遠(yuǎn)離塵源逐漸降低,這與前人研究結(jié)果不相符,可能與礦山粉塵獨特的礦物成分相關(guān),需要在下一步工作,通過采樣測試,確定研究區(qū)礦山粉塵具體礦物成分,明確隨著葉面降塵量減少,礦山塵源降塵葉面的反射率變化特征及其成因。

    5.3 葉面降塵信息提取與塵源類型區(qū)分方法效果

    本文方法將葉面降塵視為對于正常植被的一種“異?!毙畔?,基于葉面降塵像元與正常植被像元光譜特征差異,應(yīng)用特征向量主成分選擇法(FPCS)提取了葉面降塵信息,并劃分降塵強(qiáng)度;提出了針對此次研究區(qū)的塵源區(qū)分光譜指數(shù)(DSI),實現(xiàn)了對提取的葉面降塵信息塵源的判別。與前人提出的基于NDVI 的降塵量線性反演模型相比,F(xiàn)PCS 的提取結(jié)果更符合葉面降塵分布的客觀規(guī)律(提取的降塵像元主要集中在廢石堆、礦山道路、尾礦庫等礦山塵源和湖積物、沙地等自然塵源周邊)。而NDVI 反演模型結(jié)果隨機(jī)性更強(qiáng),而且在部分葉面降塵影像特征明顯的區(qū)域提取效果較差,這表現(xiàn)出經(jīng)驗?zāi)P蛯τ诓煌瑓^(qū)域的適用性較差。此次研究中未對利用塵源光譜指數(shù)(DSI)區(qū)別自然塵源與礦山塵源的結(jié)果進(jìn)行實地驗證,但典型區(qū)域(圖6(b)、圖(c))的降塵像元的塵源識別效果證實了該指數(shù)對于區(qū)分葉面降塵塵源的可行性。同時葉面降塵提取信息與礦山地物的空間相關(guān)性分析結(jié)果也驗證了本文方法的可行性,提取的礦山塵源葉面降塵像元與礦山地物的空間相關(guān)性更強(qiáng),礦山塵源葉面降塵像元在靠近礦山地物方向呈逐漸集中分布趨勢,而提取的自然塵源葉面降塵像元在礦山地物周邊呈均勻分布趨勢(圖7),與礦山地物的空間相關(guān)性弱。綜上,本文方法可用于區(qū)域性葉面降塵信息快速提取與塵源類型區(qū)分,不依賴地面采樣測試,獲取區(qū)域尺度半定量的葉面降塵分布和強(qiáng)度信息,可嘗試用于其他礦區(qū)葉面降塵信息提取。需要注意的是,此次提出的DSI 是依據(jù)自然塵源和礦山塵源光譜差異建立的,而研究區(qū)礦山塵源像元864.7 nm附近的特殊吸收特征與鐵氧化物相關(guān)。所以,應(yīng)用于其他類型金屬礦山時,需要結(jié)合礦山塵源的礦物成分具體分析。同時,降水、大風(fēng)等天氣現(xiàn)象均會對礦區(qū)植被葉面降塵范圍和強(qiáng)度造成影響,所以在影像時相選取上需要考慮數(shù)據(jù)采集前的天氣條件。另外研究中是利用DNVI 結(jié)合閾值分割提取植被范圍的基礎(chǔ)上,提取葉面降塵像元。由于哨兵二號數(shù)據(jù)空間分辨率的限制,當(dāng)植被覆蓋度較低,像元內(nèi)植被與裸土地像元混合可能造成植被光譜被掩蓋(Ma等,2020),進(jìn)而無法爭取提取部分植被覆蓋度較低的像元。同時,植被像元葉面降塵很嚴(yán)重時,可能完全掩蓋植被光譜特征,導(dǎo)致此類葉面降塵像元無法被正確識別。下一步工作中將定量分析植被覆蓋度對葉面降塵信息提取的影響,同時探索結(jié)合圖像紋理和光譜信息更準(zhǔn)確提取植被覆蓋度低或葉面降塵極嚴(yán)重的植被像元。

    6 結(jié)論

    基于哨兵二號遙感數(shù)據(jù),本文提出了一種基于光譜特征分析的葉面降塵信息快速提取和塵源類型識別的方法,成功獲取了研究區(qū)葉面降塵范圍、強(qiáng)度,和塵源類別信息。主要取得了以下結(jié)論:

    (1)驗證了哨兵二號反射率數(shù)據(jù)對區(qū)分葉面降塵強(qiáng)度和塵源類別的敏感性。與正常植被像元相比,葉面降塵像元可見光波段(B1—B4)反射率升高、近紅外波段(B7—B8A)反射率降低、紅邊位置向短波方向移動的特征;礦山塵源葉面降塵像元在864.7 nm(B8A)附近顯示出診斷性吸收特征,顯示出黃鐵礦氧化為赤鐵礦的光譜特征。

    (2)利用FPCS 方法,提取了研究區(qū)葉面降塵范圍和強(qiáng)度,提出了一種塵源光譜指數(shù)(DSI:(R(B8)+R(B11))/R(B8A),區(qū)分了葉面降塵像元的塵源類型(自然塵源與礦山塵源)。

    (3)與提取的自然塵源葉面降塵像元相比,礦山塵源葉面降塵像元與礦山地物的空間相關(guān)性更強(qiáng),驗證了本文方法的可行性。各類型礦山地物中,廢石堆、礦山道路與礦山塵源葉面降塵信息的相關(guān)性更強(qiáng),是研究區(qū)內(nèi)的主要礦山塵源。

    下一步工作中,將通過采樣測試,定量評估研究區(qū)礦山塵源葉面降塵的具體礦物成分和光譜響應(yīng)特征。同時開展其他礦產(chǎn)類型礦區(qū)的葉面降塵遙感信息提取方法研究。

    猜你喜歡
    降塵葉面反射率
    噴施葉面阻隔劑對抑制小麥吸收鎘的效果初報
    影響Mini LED板油墨層反射率的因素
    近岸水體異源遙感反射率產(chǎn)品的融合方法研究
    海洋通報(2022年4期)2022-10-10 07:40:26
    具有顏色恒常性的光譜反射率重建
    連云港市大氣降塵時空分布特征
    葉面施鈣對冬棗品質(zhì)的影響
    嶺石葉面鐵肥 缺鐵黃化
    化學(xué)腐蝕硅表面結(jié)構(gòu)反射率影響因素的研究*
    電子器件(2017年2期)2017-04-25 08:58:37
    南京市大氣降塵固碳微生物群落多樣性研究
    大白菜葉面噴施氮肥效果試驗
    上海蔬菜(2015年2期)2015-12-26 05:03:40
    亚洲精品在线美女| 日韩av不卡免费在线播放| 欧美成人精品欧美一级黄| 国产片特级美女逼逼视频| 国产成人免费无遮挡视频| 国产亚洲最大av| 爱豆传媒免费全集在线观看| 国产熟女午夜一区二区三区| 国产精品一国产av| 丝袜人妻中文字幕| 午夜福利免费观看在线| 少妇被粗大的猛进出69影院| 两个人免费观看高清视频| 80岁老熟妇乱子伦牲交| 久久精品aⅴ一区二区三区四区| 国产精品国产av在线观看| 妹子高潮喷水视频| 亚洲精品在线美女| 天天躁夜夜躁狠狠躁躁| 少妇人妻精品综合一区二区| 亚洲av福利一区| 午夜av观看不卡| 51午夜福利影视在线观看| 天天添夜夜摸| 国产免费又黄又爽又色| 免费黄网站久久成人精品| 蜜桃国产av成人99| 亚洲欧美一区二区三区久久| 亚洲av福利一区| 七月丁香在线播放| 国产男人的电影天堂91| 成人国产av品久久久| www.av在线官网国产| 亚洲中文av在线| 国产熟女欧美一区二区| 99精国产麻豆久久婷婷| 啦啦啦 在线观看视频| 赤兔流量卡办理| 亚洲人成网站在线观看播放| 亚洲精品一二三| 午夜福利一区二区在线看| 午夜福利,免费看| 尾随美女入室| 人人妻人人澡人人爽人人夜夜| 国产免费现黄频在线看| 成人毛片60女人毛片免费| 又大又爽又粗| 亚洲精品视频女| 欧美日韩综合久久久久久| 亚洲欧美一区二区三区黑人| 日韩av在线免费看完整版不卡| 国产黄色免费在线视频| 亚洲av国产av综合av卡| 在线观看免费日韩欧美大片| 亚洲精品aⅴ在线观看| av在线观看视频网站免费| 黄片无遮挡物在线观看| 欧美日韩视频精品一区| 99热国产这里只有精品6| 最新在线观看一区二区三区 | 啦啦啦中文免费视频观看日本| 国产精品亚洲av一区麻豆 | 尾随美女入室| 国产欧美亚洲国产| 天堂中文最新版在线下载| 国产成人精品久久二区二区91 | 视频在线观看一区二区三区| 少妇人妻 视频| 免费不卡黄色视频| 国产成人91sexporn| 一级片'在线观看视频| 黄色怎么调成土黄色| 亚洲国产欧美在线一区| 777米奇影视久久| 黑人猛操日本美女一级片| 啦啦啦中文免费视频观看日本| 亚洲欧美清纯卡通| 免费看av在线观看网站| 国产亚洲av片在线观看秒播厂| 少妇 在线观看| 777久久人妻少妇嫩草av网站| 99热全是精品| 天天躁夜夜躁狠狠躁躁| 91成人精品电影| 99热网站在线观看| 午夜福利,免费看| 亚洲av在线观看美女高潮| 老司机亚洲免费影院| 夫妻午夜视频| 各种免费的搞黄视频| 欧美日韩成人在线一区二区| 叶爱在线成人免费视频播放| av福利片在线| 国产一区二区在线观看av| 一区二区三区四区激情视频| 久热爱精品视频在线9| av在线老鸭窝| 欧美变态另类bdsm刘玥| 欧美激情高清一区二区三区 | 婷婷成人精品国产| 亚洲精品一二三| 纵有疾风起免费观看全集完整版| 少妇人妻久久综合中文| 视频区图区小说| 99热网站在线观看| 亚洲精品久久午夜乱码| a级毛片黄视频| 国产99久久九九免费精品| 欧美精品亚洲一区二区| 亚洲五月色婷婷综合| 狠狠精品人妻久久久久久综合| 一本一本久久a久久精品综合妖精| 久久狼人影院| 精品国产一区二区三区久久久樱花| 成人手机av| 欧美日韩一级在线毛片| 亚洲欧美清纯卡通| 男人爽女人下面视频在线观看| 亚洲国产精品成人久久小说| 丰满饥渴人妻一区二区三| 激情五月婷婷亚洲| 黄片无遮挡物在线观看| 大片免费播放器 马上看| 99久久精品国产亚洲精品| 青春草亚洲视频在线观看| 久久97久久精品| 一边摸一边做爽爽视频免费| www日本在线高清视频| 三上悠亚av全集在线观看| 热99久久久久精品小说推荐| 狠狠婷婷综合久久久久久88av| 久久久久精品久久久久真实原创| 久久影院123| 国产片特级美女逼逼视频| av在线app专区| 中文字幕亚洲精品专区| 精品一区二区三卡| 国产一级毛片在线| 麻豆精品久久久久久蜜桃| 久久久欧美国产精品| 你懂的网址亚洲精品在线观看| 亚洲熟女精品中文字幕| 亚洲综合色网址| 亚洲精品av麻豆狂野| 一级a爱视频在线免费观看| 欧美亚洲日本最大视频资源| 91成人精品电影| 尾随美女入室| 大话2 男鬼变身卡| 精品一区二区三区四区五区乱码 | 亚洲精品自拍成人| 亚洲精品国产一区二区精华液| 欧美日韩视频高清一区二区三区二| 中文字幕另类日韩欧美亚洲嫩草| 十八禁人妻一区二区| 亚洲精品乱久久久久久| 日日摸夜夜添夜夜爱| av电影中文网址| 纯流量卡能插随身wifi吗| 午夜福利乱码中文字幕| 国产精品av久久久久免费| 欧美精品一区二区免费开放| 亚洲av福利一区| 亚洲国产最新在线播放| 蜜桃国产av成人99| 青春草视频在线免费观看| 99re6热这里在线精品视频| 亚洲第一区二区三区不卡| 国产精品无大码| 最新在线观看一区二区三区 | 国产高清不卡午夜福利| 亚洲色图 男人天堂 中文字幕| 亚洲一级一片aⅴ在线观看| 久久久久精品国产欧美久久久 | 亚洲精品第二区| 人人妻,人人澡人人爽秒播 | 久久精品久久久久久久性| 国产免费又黄又爽又色| 精品酒店卫生间| 亚洲av国产av综合av卡| 精品一区二区免费观看| 91老司机精品| 色婷婷久久久亚洲欧美| 国产精品一二三区在线看| 美女午夜性视频免费| 欧美最新免费一区二区三区| 如日韩欧美国产精品一区二区三区| 美女扒开内裤让男人捅视频| 丁香六月欧美| 欧美国产精品va在线观看不卡| 亚洲av成人不卡在线观看播放网 | 午夜激情久久久久久久| 精品久久久精品久久久| 国产精品蜜桃在线观看| 亚洲视频免费观看视频| 一级毛片黄色毛片免费观看视频| 国产 一区精品| 午夜久久久在线观看| 日韩免费高清中文字幕av| 如何舔出高潮| 久热这里只有精品99| 国产成人精品久久久久久| 视频区图区小说| 精品少妇久久久久久888优播| a 毛片基地| av在线app专区| 久久久久久久久久久免费av| 国产 一区精品| 久久热在线av| 亚洲欧美一区二区三区国产| 乱人伦中国视频| 亚洲成av片中文字幕在线观看| 电影成人av| 久久久国产一区二区| 一级毛片我不卡| 免费黄网站久久成人精品| 岛国毛片在线播放| 伊人亚洲综合成人网| 国产精品免费大片| 亚洲自偷自拍图片 自拍| 精品亚洲乱码少妇综合久久| 99re6热这里在线精品视频| 18禁裸乳无遮挡动漫免费视频| 少妇猛男粗大的猛烈进出视频| 欧美日韩亚洲高清精品| 中文字幕人妻丝袜一区二区 | 777久久人妻少妇嫩草av网站| 青春草视频在线免费观看| 欧美在线一区亚洲| 亚洲熟女精品中文字幕| www.精华液| 国产精品嫩草影院av在线观看| 狠狠精品人妻久久久久久综合| 日韩,欧美,国产一区二区三区| 亚洲美女黄色视频免费看| 啦啦啦在线观看免费高清www| 亚洲成人手机| 黄色 视频免费看| 欧美人与善性xxx| 成人手机av| 人人澡人人妻人| 亚洲五月色婷婷综合| 少妇被粗大的猛进出69影院| 色精品久久人妻99蜜桃| 亚洲国产中文字幕在线视频| 新久久久久国产一级毛片| 亚洲精品av麻豆狂野| 在现免费观看毛片| 久久久久久人人人人人| 亚洲欧美成人精品一区二区| 你懂的网址亚洲精品在线观看| 日韩不卡一区二区三区视频在线| 成人漫画全彩无遮挡| 日韩免费高清中文字幕av| 亚洲专区中文字幕在线 | 欧美日韩成人在线一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 王馨瑶露胸无遮挡在线观看| 国产精品久久久久久久久免| 久久毛片免费看一区二区三区| 国产福利在线免费观看视频| 精品午夜福利在线看| 一边摸一边抽搐一进一出视频| 搡老岳熟女国产| 美女脱内裤让男人舔精品视频| 黄片无遮挡物在线观看| 免费观看性生交大片5| 久久精品久久久久久久性| 国产精品.久久久| 又大又爽又粗| 成人影院久久| 亚洲精品在线美女| 精品人妻熟女毛片av久久网站| 波多野结衣一区麻豆| 亚洲,一卡二卡三卡| 亚洲专区中文字幕在线 | 天堂8中文在线网| 久久久久久免费高清国产稀缺| 免费少妇av软件| 亚洲精品日韩在线中文字幕| 亚洲av电影在线观看一区二区三区| 国产欧美日韩一区二区三区在线| 久久精品国产亚洲av涩爱| 国产熟女午夜一区二区三区| 美女福利国产在线| 亚洲图色成人| 国产免费福利视频在线观看| 久久久精品94久久精品| 国产精品国产三级国产专区5o| 人妻一区二区av| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜福利视频精品| 成年av动漫网址| 国产伦人伦偷精品视频| 最近最新中文字幕免费大全7| 我要看黄色一级片免费的| 一区福利在线观看| 国产黄色免费在线视频| 欧美另类一区| av不卡在线播放| 精品一品国产午夜福利视频| 日韩av不卡免费在线播放| 一级毛片电影观看| 青春草国产在线视频| 最近中文字幕2019免费版| 日本猛色少妇xxxxx猛交久久| 日韩中文字幕欧美一区二区 | 精品少妇黑人巨大在线播放| 香蕉国产在线看| 操美女的视频在线观看| 老司机亚洲免费影院| 两个人免费观看高清视频| 国产精品无大码| 久久 成人 亚洲| 涩涩av久久男人的天堂| 大话2 男鬼变身卡| 亚洲av电影在线观看一区二区三区| 亚洲精品日本国产第一区| 亚洲欧美色中文字幕在线| 日韩一本色道免费dvd| av天堂久久9| 婷婷色av中文字幕| 人妻一区二区av| 久久久精品国产亚洲av高清涩受| 亚洲四区av| 91精品三级在线观看| 操美女的视频在线观看| 水蜜桃什么品种好| 欧美老熟妇乱子伦牲交| 欧美日韩一级在线毛片| av不卡在线播放| 宅男免费午夜| 亚洲av成人精品一二三区| 老司机影院成人| 无遮挡黄片免费观看| 日本91视频免费播放| 丝袜人妻中文字幕| 国产成人免费无遮挡视频| 免费看不卡的av| 各种免费的搞黄视频| 久久久精品区二区三区| 人体艺术视频欧美日本| av在线播放精品| 国产精品久久久久久精品古装| 80岁老熟妇乱子伦牲交| 少妇人妻久久综合中文| 又黄又粗又硬又大视频| 免费少妇av软件| 51午夜福利影视在线观看| 最黄视频免费看| 国产片特级美女逼逼视频| 80岁老熟妇乱子伦牲交| 99精国产麻豆久久婷婷| 国产精品免费大片| 最近手机中文字幕大全| 欧美黑人欧美精品刺激| 成人国语在线视频| av不卡在线播放| 日韩精品免费视频一区二区三区| 亚洲一码二码三码区别大吗| av片东京热男人的天堂| 亚洲欧美日韩另类电影网站| 男人添女人高潮全过程视频| 久久久久久久久久久久大奶| 9色porny在线观看| 最近手机中文字幕大全| 丝袜喷水一区| 成人亚洲精品一区在线观看| 成年女人毛片免费观看观看9 | 99国产精品免费福利视频| 精品福利永久在线观看| 国产成人欧美在线观看 | 精品第一国产精品| 免费看av在线观看网站| 成人免费观看视频高清| 少妇被粗大猛烈的视频| 国产黄色视频一区二区在线观看| 国产成人午夜福利电影在线观看| 欧美日韩福利视频一区二区| 免费日韩欧美在线观看| 亚洲精品一二三| 日本wwww免费看| 亚洲四区av| 尾随美女入室| 精品一区二区三区四区五区乱码 | 亚洲欧美一区二区三区久久| av在线app专区| 丝袜人妻中文字幕| 最近手机中文字幕大全| 美女午夜性视频免费| 免费在线观看黄色视频的| 亚洲三区欧美一区| 久久人妻熟女aⅴ| 日日摸夜夜添夜夜爱| 国产精品二区激情视频| 欧美精品一区二区免费开放| 制服人妻中文乱码| 丰满饥渴人妻一区二区三| 欧美老熟妇乱子伦牲交| 成人黄色视频免费在线看| 亚洲 欧美一区二区三区| 女人被躁到高潮嗷嗷叫费观| 亚洲av欧美aⅴ国产| 少妇的丰满在线观看| 蜜桃国产av成人99| 免费高清在线观看视频在线观看| 色婷婷久久久亚洲欧美| 久久久欧美国产精品| 黑人猛操日本美女一级片| 美女中出高潮动态图| 可以免费在线观看a视频的电影网站 | 美女大奶头黄色视频| 精品少妇久久久久久888优播| 国产乱人偷精品视频| 欧美日韩福利视频一区二区| 免费看av在线观看网站| 十八禁人妻一区二区| 狂野欧美激情性bbbbbb| 成人毛片60女人毛片免费| 中文字幕精品免费在线观看视频| 日韩熟女老妇一区二区性免费视频| 高清不卡的av网站| 亚洲第一青青草原| 最新在线观看一区二区三区 | 国产日韩欧美视频二区| 国产福利在线免费观看视频| 久久青草综合色| 亚洲欧美色中文字幕在线| 赤兔流量卡办理| 一本色道久久久久久精品综合| 男人舔女人的私密视频| 国产精品久久久久久人妻精品电影 | 中文乱码字字幕精品一区二区三区| 欧美精品高潮呻吟av久久| 亚洲欧美精品自产自拍| 国产亚洲午夜精品一区二区久久| 欧美日韩精品网址| 看免费成人av毛片| 国产亚洲最大av| 精品人妻一区二区三区麻豆| 色精品久久人妻99蜜桃| 老司机影院成人| 十八禁网站网址无遮挡| 免费久久久久久久精品成人欧美视频| 午夜福利视频在线观看免费| 中文字幕制服av| 男人舔女人的私密视频| 成人黄色视频免费在线看| 一边摸一边做爽爽视频免费| 久久99一区二区三区| 一级毛片我不卡| 午夜福利在线免费观看网站| 国产高清国产精品国产三级| 巨乳人妻的诱惑在线观看| 十分钟在线观看高清视频www| 亚洲av成人不卡在线观看播放网 | 国产免费视频播放在线视频| 欧美日韩av久久| 最近中文字幕2019免费版| 观看美女的网站| 欧美老熟妇乱子伦牲交| 欧美日韩亚洲国产一区二区在线观看 | 高清欧美精品videossex| 女性生殖器流出的白浆| 久久人人爽av亚洲精品天堂| 国产午夜精品一二区理论片| 精品午夜福利在线看| 国产亚洲av片在线观看秒播厂| 少妇人妻久久综合中文| 99久久人妻综合| 日本wwww免费看| 韩国精品一区二区三区| 国产精品人妻久久久影院| 久久久欧美国产精品| 啦啦啦在线免费观看视频4| 香蕉丝袜av| 久热爱精品视频在线9| 国产 一区精品| 女人高潮潮喷娇喘18禁视频| www.熟女人妻精品国产| 成人免费观看视频高清| 曰老女人黄片| 亚洲美女黄色视频免费看| 亚洲伊人色综图| 国产激情久久老熟女| 亚洲在久久综合| 一二三四中文在线观看免费高清| 国产福利在线免费观看视频| 亚洲精品中文字幕在线视频| 操美女的视频在线观看| 国产精品香港三级国产av潘金莲 | 久久午夜综合久久蜜桃| 中文字幕高清在线视频| 日韩精品免费视频一区二区三区| 又粗又硬又长又爽又黄的视频| 搡老岳熟女国产| 国产激情久久老熟女| 国产午夜精品一二区理论片| 欧美精品av麻豆av| 亚洲,欧美精品.| av又黄又爽大尺度在线免费看| 亚洲第一av免费看| 在线观看三级黄色| 久久国产精品男人的天堂亚洲| 国产亚洲欧美精品永久| 国产精品香港三级国产av潘金莲 | 精品午夜福利在线看| 久久国产精品男人的天堂亚洲| 黑人猛操日本美女一级片| 人妻 亚洲 视频| 午夜免费男女啪啪视频观看| 哪个播放器可以免费观看大片| 国产乱来视频区| 两性夫妻黄色片| 丝袜在线中文字幕| 两性夫妻黄色片| 国产精品一二三区在线看| 久久精品亚洲av国产电影网| 亚洲精品成人av观看孕妇| 婷婷色麻豆天堂久久| 黄色视频在线播放观看不卡| 免费高清在线观看视频在线观看| 飞空精品影院首页| 欧美日韩成人在线一区二区| 女的被弄到高潮叫床怎么办| 亚洲国产av影院在线观看| 国产精品久久久av美女十八| 啦啦啦在线观看免费高清www| 国产日韩欧美在线精品| 国产成人欧美在线观看 | 水蜜桃什么品种好| 一区二区日韩欧美中文字幕| 在线观看三级黄色| 老熟女久久久| 一级,二级,三级黄色视频| 国产精品国产三级国产专区5o| 街头女战士在线观看网站| 丝袜在线中文字幕| 精品国产乱码久久久久久小说| 国产一级毛片在线| 欧美人与性动交α欧美精品济南到| 亚洲精品国产av蜜桃| 操出白浆在线播放| 悠悠久久av| 视频区图区小说| 天堂俺去俺来也www色官网| 亚洲美女黄色视频免费看| 晚上一个人看的免费电影| 欧美日韩一级在线毛片| 极品少妇高潮喷水抽搐| 国产精品免费大片| a级片在线免费高清观看视频| 精品少妇久久久久久888优播| 看十八女毛片水多多多| 在现免费观看毛片| 大香蕉久久成人网| 九九爱精品视频在线观看| 18禁裸乳无遮挡动漫免费视频| 久热这里只有精品99| xxxhd国产人妻xxx| 最近中文字幕2019免费版| 欧美日韩精品网址| 哪个播放器可以免费观看大片| 男女午夜视频在线观看| 久久久国产一区二区| 欧美黄色片欧美黄色片| 丝袜美足系列| 叶爱在线成人免费视频播放| 最近2019中文字幕mv第一页| 国产精品一国产av| 男女午夜视频在线观看| 男女国产视频网站| 在线天堂中文资源库| 超碰成人久久| 国产精品偷伦视频观看了| 亚洲伊人色综图| 亚洲国产中文字幕在线视频| h视频一区二区三区| 精品第一国产精品| 色94色欧美一区二区| 日本欧美国产在线视频| 国产精品成人在线| 伊人久久国产一区二区| 亚洲一级一片aⅴ在线观看| 精品国产一区二区三区四区第35| 中文字幕最新亚洲高清| 精品酒店卫生间| 欧美国产精品va在线观看不卡| 成年美女黄网站色视频大全免费| 国产片特级美女逼逼视频| 只有这里有精品99| 天天操日日干夜夜撸| 日韩一区二区三区影片| 日韩伦理黄色片| 午夜影院在线不卡| 成人影院久久| 亚洲中文av在线| 日韩免费高清中文字幕av| 久久免费观看电影| 在线亚洲精品国产二区图片欧美| 少妇人妻 视频| kizo精华| 色婷婷久久久亚洲欧美| 午夜老司机福利片| 欧美黑人欧美精品刺激| 丝袜脚勾引网站| 女人精品久久久久毛片| 在线天堂中文资源库| 秋霞在线观看毛片| 在线观看三级黄色| 免费高清在线观看日韩| 蜜桃国产av成人99| 激情五月婷婷亚洲| 制服丝袜香蕉在线| 国产在线一区二区三区精|