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

    一種基于陣風鋒時空變化特征的大風識別算法研究

    2023-10-31 09:49:36吳奕霄王珂清周學東
    雷達科學與技術 2023年5期
    關鍵詞:陣風反射率分值

    吳奕霄,王珂清,程 婷,項 瑛,周學東,徐 芬

    (1.江蘇省氣候中心,江蘇南京 210041;2.江蘇省人工影響天氣中心,江蘇南京 210041;3.江蘇省氣象信息中心,江蘇南京 210041;4.南京氣象科技創(chuàng)新研究院,江蘇南京 210041)

    0 引言

    陣風鋒過境伴有氣壓升高、風速突增和溫度陡降等現(xiàn)象,陣風鋒及其母雷暴經(jīng)常會產(chǎn)生地面大風,引發(fā)人員傷亡或者經(jīng)濟損失,因此陣風鋒研究對防災、減災有重要意義[1]。提前識別和預警陣風鋒有利于采取有效的措施避免危險事故發(fā)生,進而降低和減少經(jīng)濟損失和人員傷亡。雷達目標檢測能夠判斷回波信號中目標存在與否,并提取目標位置信息[2-3]。因此,利用雷達進行易引發(fā)大風的陣風鋒自動檢測技術研究引起了廣泛關注。模板匹配技術是在設定一個目標模板后,與具有相似目標的圖像進行比對,以獲取目標模板在圖像中的位置,以此達到圖像識別的目的,其優(yōu)點是匹配的算法更加簡單,在圖像變化不大的情況下,識別率理想[4-5]。

    在天氣雷達反射率因子圖像上陣風鋒表現(xiàn)為一條高度較低(低仰角層觀測)、形態(tài)呈弧線結構和具有移動特點這3 個特征的窄帶回波。在雷達徑向速度數(shù)據(jù)場上呈現(xiàn)速度輻合特征。Uyeda等[6]首先提出陣風鋒自動探測算法(Automatic Detection of Gust Fronts,ADGF),該算法只使用沿雷達徑向的速度輻合,將徑向切變超過預設閾值的段組合在一起,用曲線擬合出每個段上的最大切變值,這條曲線就代表探測到的陣風鋒。后來提出了改進的陣風鋒探測算法(Gust Front Detection Algorithm,GFDA)[7-9],利用兩個低仰角垂直關聯(lián)減少了虛假預警,該算法被用于美國TDWR 雷達業(yè)務系統(tǒng)。美國強風暴實驗室和麻省理工學院林肯實驗室聯(lián)合開發(fā)了改進的陣風鋒自動探測算法[10-11],包含了陣風鋒雷達回波細線探測,通過尋找局地最大反射率因子與陣風鋒反射率因子和寬度一致性探測陣風鋒回波細線。Delanoy等[12]開發(fā)了機器智能識別陣風鋒算法(Machine Intelligent Gust Front Algorithm,MIGFA),該算法首次使用模板函數(shù)計算雷達回波圖和速度標準偏差得分,并設計了相應的質(zhì)量控制算法抑制虛假信號,最后在組合得分圖像上提取陣風鋒。隨著需求的提升和技術的發(fā)展,一些新的或改進的陣風鋒識別技術也得到全面發(fā)展,徐芬等[13]通過統(tǒng)計江蘇3種窄帶回波分布特征和徑向波形特征,在此基礎上設計了基于回波平坦性計算和陣風鋒徑向波形尋找算法的陣風鋒自動識別算法。冷亮等[14]基于陣風鋒的3 個特征,提出了基于數(shù)字圖像處理技術中的數(shù)學形態(tài)學算法的陣風鋒自動識別算法,通過建立雷達觀測范圍內(nèi)的二維特征量場,對其處理后得到一條反映陣風鋒骨架特征的曲線。徐芬等[15]改進和調(diào)整了MIGFA 算法,設計了較高仰角(1.5°/2.4°)反射率的陣風鋒細線函數(shù)模板,替代原空間差分反射率函數(shù)模板,新的算法對于CINRAD/SA 雷達探測到的江蘇沿江地區(qū)陣風鋒回波的識別更有針對性。HWANG 等[16]使用模糊邏輯理論開發(fā)了新的神經(jīng)模糊陣風鋒檢測算法(Neuro-Fuzzy Gust Front Detection Algorithm,NFGDA),該算法通過使用神經(jīng)網(wǎng)絡的訓練進行優(yōu)化,對于陣風鋒的識別率高達92%。Yuan等[17]基于統(tǒng)計特征分析和圖像處理技術,提出一種基于本地二進制雙模板(Local Binary with Dual-Template,LBDT)算法改進的雷達數(shù)據(jù)自動檢測陣風鋒的方法。該算法對反射率值和陣風鋒寬度的變化不敏感,在自動識別陣風鋒過程中檢測概率高、誤報率低。徐月飛等[18]基于Faster RCNN 算法和InceptionV2 模型,提出了改進的深度卷積網(wǎng)絡陣風鋒識別算法(d-CNNGFDA),利用窄帶弱回波的繞雷達中心旋轉不變性的特性,只對垂直走向的窄帶弱回波進行標記,擴充了數(shù)據(jù)量用以提升復雜回波中陣風鋒引起的窄帶弱回波的自動識別能力。

    上述這些方法均需選取陣風鋒樣本符合各自算法設定條件上才可以較好的識別,各有使用局限性。本文試圖融合易引發(fā)大風的陣風鋒通用特征:雷達回波圖像所在高度較低、圖像特征呈弧線結構和具有移動這3 個特征,設計能表征這3個特征的函數(shù)模板,再采用風暴回波特征抑制技術和其他弧線回波特征抑制技術以降低檢測虛警率,最終實現(xiàn)更具普適性的一套陣風鋒檢測技術。

    1 識別方法研究

    1.1 算法流程概述

    根據(jù)觀測事實以及學者們總結[12]的易引發(fā)大風的陣風鋒雷達反射率因子圖像具有以下特征:1)陣風鋒雷達反射率因子圖像呈細長線狀,其數(shù)值一般小于25 dBz;2)陣風鋒圖像高度較低;3)呈現(xiàn)移動的特點,如圖1方框內(nèi)弧線狀圖像所示。根據(jù)這3個時空變化特征,采用相應預處理方法提高雷達反射率因子圖像質(zhì)量后,借鑒模板匹配方法[12]得到3種不同的得分值結果。

    圖1 陣風鋒天氣雷達反射率因子圖像(方框內(nèi)弧線狀圖像)

    融合3種時空變化特征匹配結果后,采用區(qū)域增長法提取陣風鋒弧線圖像結構,達到自動檢測陣風鋒的目的?;陉囷L鋒時空變化特征的大風識別算法流程如圖2所示,由該圖可見通過對包含有不同特征的雷達回波圖像進行相應的模板匹配后,融合多組得分值,經(jīng)過兩次有效提取,最終實現(xiàn)對易引發(fā)大風的陣風鋒雷達反射率因子圖像的提取。

    在利用函數(shù)模板提取陣風鋒圖像特征前,為提高圖像提取成功率,對數(shù)據(jù)質(zhì)量進行如下預處理:將0.5°和1.5°仰角的雷達反射率因子圖像采用9 點平滑濾除孤立點、奇異點和回波圖像中的空點。將這兩個仰角的數(shù)據(jù)從極坐標下轉換到直角坐標下,空間分辨率為1 km×1 km。后續(xù)基于函數(shù)模板和得分函數(shù)的計算均基于直角坐標系下的雷達回波圖像。

    1.2 函數(shù)模板設計

    根據(jù)陣風鋒圖像在不同仰角、仰角間圖像特征差異、同一仰角前后時刻圖像差異等設計相應函數(shù)模板提取能夠反映陣風鋒的圖像特征。根據(jù)陣風鋒在空間變化上的特征設計了陣風鋒線狀函數(shù)模板;根據(jù)陣風鋒在空間上線狀差異特征分別設計了風暴邊緣函數(shù)模板和放射狀回波圖像函數(shù)模板;根據(jù)陣風鋒在時間變化上的特征設計了陣風鋒差分反射率因子圖像函數(shù)模板。以上幾種函數(shù)模板和得分函數(shù)分別設計如下:由于陣風鋒天氣雷達反射率因子圖像特征呈細長弧線結構,根據(jù)這一特征設計函數(shù)模板和得分函數(shù)如圖2所示,利用相應的得分函數(shù)可得到高的陣風鋒圖像得分值,抑制其他圖像得分值。該設計思路是依據(jù)盡管陣風鋒圖像長度和彎曲弧度不等,但將其分為若干段,每一段均近似為直線,因此圖3(a)中的模板被設計為直線段主要基于假設陣風鋒圖像的一個段為近似直線。黑色和灰色格點代表參與計算的點,灰色代表模板中陣風鋒回波圖像上的點,對應圖3(b)中的得分函數(shù)1。因為陣風鋒回波一般在25 dBz 以下,所以給25 dBz 以下回波一個較高的得分值,特別是9 dBz 附近[10],是陣風鋒最常出現(xiàn)的反射率因子值。給大于35 dBz 的回波一個緩慢下降的得分值,因為陣風鋒在某些情況下會與風暴邊緣結合在一起,而陣風鋒一般不會出現(xiàn)負值,所以給負的反射率因子值一個負得分值。黑色的點代表陣風鋒回波兩側的點,對應圖3(b)中的得分函數(shù)0,給得分函數(shù)0 低反射率因子一個高的分值,而給高反射率因子一個低得分值。需要說明的是,函數(shù)模板和得分函數(shù)的長寬是可調(diào)的,因為各地發(fā)生的陣風鋒長度、寬度和反射率因子值可能不盡相同,筆者結合江蘇地區(qū)陣風鋒反射率因子值和相關文獻確定的得分函數(shù)值僅供參考。

    因為陣風鋒回波圖像出現(xiàn)的高度較低,一般為3 km 以下,可利用高低仰角層差分反射率因子來識別陣風鋒,這樣做的好處是排除較多其他回波干擾,特別是貫穿了高低層的線狀對流風暴圖像。將上下兩個仰角的反射率因子圖像經(jīng)上述數(shù)據(jù)預處理方法質(zhì)控后,上下仰角同一格點數(shù)據(jù)相減得到差分反射率因子圖像(圖4(a))。如圖4(a)所示,圖中大部分回波已被消除,陣風鋒回波圖像基本被保留。根據(jù)陣風鋒差分反射率因子圖像特征,設計圖5中的函數(shù)模板和得分函數(shù)對差分反射率因子圖像進行特征提取。與圖3 陣風鋒線狀函數(shù)模板和得分函數(shù)對比可見,采用的函數(shù)模板不變,得分函數(shù)根據(jù)陣風鋒差分反射率圖像特征作了一定調(diào)整。同理,因為易引發(fā)大風的陣風鋒具有移動的特征,所以可以利用前后時刻差分反射率因子圖像特征來識別陣風鋒。將前后連續(xù)兩個時刻的反射率因子圖像經(jīng)上述數(shù)據(jù)預處理方法質(zhì)控后,前后時刻同一格點數(shù)據(jù)相減得到陣風鋒差分反射率因子圖像(圖4(b)),仍舊采用圖5的函數(shù)模板和得分函數(shù)對前后時刻陣風鋒差分反射率圖像特征進行提取。

    圖4 陣風鋒天氣雷達反射率因子圖像(方框內(nèi)弧線狀圖像)

    圖5 陣風鋒差分反射率因子圖像函數(shù)模板

    圖4(b)中的前后時刻差分反射率因子排除了較多的回波,卻保留了風暴邊緣圖像,該風暴邊緣呈弧線狀,很容易被誤識別為陣風鋒圖像。本文設計風暴函數(shù)模板和相應得分函數(shù)用以削弱線狀風暴邊緣對識別結果的影響。如圖6(a)所示,函數(shù)模板被設計為一個近似的圓形,圖中的每個格點代表1 km 距離,該函數(shù)能夠抑制線狀回波圖像得分值,突出呈片狀的風暴和層云圖像得分值。圖6(b)為風暴邊緣得分函數(shù)圖。

    圖6 風暴邊緣函數(shù)模板

    另外,在圖1 和圖4 的圖像中都受到一種放射狀回波的干擾。如圖4(a)第一象限中能看到一個沿雷達掃描方向的線狀圖像(圖中箭頭所示),這是由于放射狀回波干擾所致。為了排除放射狀回波干擾,設計如圖7所示的函數(shù)模板和得分函數(shù),圖7(a)中每個格點代表1 km,得分函數(shù)根據(jù)放射狀回波圖像特征設計。函數(shù)1 代表在放射狀回波圖像上,函數(shù)0代表在放射狀回波圖像兩側。因為放射狀回波干擾圖像一般只會出現(xiàn)在雷達掃描方向,所以該模板僅在極坐標反射率圖像上沿雷達掃描方向進行計算。再分別利用圖3、圖5 函數(shù)模板和得分函數(shù)處理每幅雷達反射率圖像得到的得分值結果中減去使用該模板計算得到的結果,即可排除放射狀回波干擾。

    圖7 放射狀回波圖像函數(shù)模板

    為了讓得分值結果便于提取,本文還設計了如圖8所示的函數(shù)模板和得分函數(shù)對融合的得分值結果進行平滑處理。如圖8(a)所示,以黑色圓點為中心點,灰色方格代表陣風鋒得分值上的點,黑色方格代表陣風鋒得分值兩側的點,模板中心點兩側設計了較多的灰色格點,目的是讓陣風鋒得分值更加連續(xù)。

    圖8 平滑圖像函數(shù)模板

    由于陣風鋒圖像的朝向是不確定的,在某個格點上,上述圖3、圖5、圖8 中的函數(shù)模板以黑色小圓點為中心,讓其以10°的間隔旋轉180°,得到不同角度下模板匹配的得分值,最后在這個格點上只保留最高得分值。然后以1 個格點的步長移過所有格點,并在所有格點上重復以上步驟。旋轉公式如下:

    式中,θ代表旋轉角度,x'和y'代表旋轉后新坐標,x和y代表原坐標。

    經(jīng)過上述多個函數(shù)模板和得分函數(shù)處理后,對多組得分值進行融合處理,并對得分值結果作二值化處理,得分值大于0.5 賦值為1,小于0.5 則賦值為0。為提取屬于陣風鋒的圖像,本文借鑒區(qū)域增長法[11]對二值化圖像進行處理,得到的候選陣風鋒圖像如果長度小于6 km,則排除這個候選陣風鋒。然后再采用一個有條件的區(qū)域增長法對初次提取的陣風鋒進行處理:如果一個格點滿足與原候選陣風鋒包含的點相鄰,得分值不低于0.2且該點得分值角度與原候選陣風鋒平均角度相差不超過30°,則認為該點是候選陣風鋒中的一個點,并入候選陣風鋒中。如果候選陣風鋒長度超過20 km,則認為是陣風鋒,否則排除。

    2 識別方法檢測

    以2012年5月16日鹽城新一代多普勒天氣雷達探測到的陣風鋒圖像為例(圖1方框內(nèi))。如圖1(a)中0.5°仰角雷達站東南方向的陣風鋒圖像被風暴邊緣掩蓋,圖1(b)中1.5°仰角雷達站東南方向的陣風鋒圖像較明顯,基本未受影響,但是一般陣風鋒圖像出現(xiàn)高度較低,抬高仰角至1.5°時,隨著離開雷達站距離增加,探測高度也在增加,在雷達站東偏北約25°方向第二個距離圈處已探測不到陣風鋒圖像??梢?.5°雷達資料探測高度較低,在近距離處陣風鋒圖像易被地物回波所掩蓋;1.5°仰角探測時,受地物回波影響較小,卻由于遠距離處探測高度較高而不能探測到陣風鋒。由此可以看出,單仰角的雷達資料并不一定能探測到完整的陣風鋒圖像,這也是文中利用了兩個仰角的圖像進行陣風鋒自動識別的原因。

    采用本文圖3 設計的函數(shù)模板和得分函數(shù)分別對0.5°仰角和1.5°仰角雷達圖像進行處理后,用灰度值表達計算得到的得分值如圖9(a)和圖9(b)所示,圖像中越高的灰度值代表得分值越高。將圖9 與圖1 實況圖對比可以看出,陣風鋒回波處明顯得到了較高的得分值,說明模板函數(shù)匹配技術能夠在陣風鋒識別中發(fā)揮重要作用。圖9(a)和圖9(b)中部分非陣風鋒圖像也得到了較高的得分值,這也是陣風鋒自動識別中最困難的問題所在——相似回波較多且難以區(qū)分,所以文中基于多種陣風鋒特征來對其進行自動識別。圖9(c)是由圖9(a)和圖9(b)綜合后的組合得分值。綜合兩個仰角的探測優(yōu)勢,將得分值圖像組合起來,圖9(c)中得到較完整、連續(xù)的陣風鋒得分值,但同時也組合了兩個得分值圖像中的干擾回波圖像得分值。

    圖9 陣風鋒線狀函數(shù)模板處理的0.5°和1.5°仰角雷達圖像

    再利用圖5的函數(shù)模板和得分函數(shù)分別對圖4(a)和圖4(b)進行得分計算,得到圖10兩個得分值圖像,從圖10 中可以看出陣風鋒圖像處得分值明顯較高,且具有較好的連續(xù)性,與圖4 中陣風鋒圖像所在位置和形狀有較高的一致性。

    圖10 陣風鋒差分反射率因子圖像函數(shù)模板計算的兩個得分值圖像

    根據(jù)圖2 的算法流程,利用圖6 函數(shù)模板和得分函數(shù)處理前后時刻差分反射率圖像(圖4(b))得到風暴邊緣圖像得分值(圖略),用圖10(b)得分值結果減去風暴邊緣圖像得分值,降低風暴邊緣圖像對識別結果的影響。再分別對圖9(c)、圖10(a)和圖10(b)得分值減去放射狀回波圖像得分值結果后,采用加權平均的方式融合3種得分值得到綜合得分值結果,如圖11(a)所示。由該圖可見,雷達站南偏東約20°方向的放射狀回波圖像已得到較好的抑制。利用圖8 平滑圖像函數(shù)模板和得分函數(shù)處理圖11(a),得到圖11(b)最終的得分值結果。由圖11(b)可見,本文設計的檢測技術成功突出了陣風鋒得分值,而抑制其他得分值。在實際中可根據(jù)算法在本地運行情況調(diào)整權重,進一步提高識別成功率。最后,利用區(qū)域增長法處理圖11(b),得到結果如圖11(c),對比實況圖1 可以看出,算法能夠較好地識別出陣風鋒,且識別結果與原陣風鋒位置和形狀比較一致。

    圖11 多次處理后得分值及自動識別結果

    3 結束語

    本文根據(jù)我國易引發(fā)大風的陣風鋒雷達回波圖像高度較低、呈弧線結構和具有移動特點的時空變化特征,在3 種不同反射率因子圖像上,采用模板匹配技術分別對其進行處理,得到3種不同得分值結果。對3 種不同得分值結果進行加權平均后,采用能夠使得分值平滑細化的函數(shù)模板技術進行處理,最后采用兩次不同的區(qū)域增長法提取陣風鋒,并排除虛假預警。在得到3種得分值結果的過程中,還采用了風暴邊緣圖像抑制模板技術和放射狀回波圖像模板抑制技術。結果表明,算法能夠有效識別陣風鋒,識別結果與原陣風鋒位置和形狀比較一致。

    作為會對生活和生產(chǎn)帶來嚴重影響的天氣因素,陣風鋒常年帶來損失不容忽視。本文所設計的陣風鋒模板匹配技術為能準確地監(jiān)測和預警陣風鋒導致的災害性天氣提供多一種可能。目前本文中涉及的算法函數(shù)模板以及得分函數(shù)中的閾值集均來自江淮地區(qū)多年陣風鋒圖像特征統(tǒng)計結果,適用于該地區(qū)的陣風鋒特征識別,若要在其他地區(qū)進行推廣,還需要結合其他地區(qū)陣風鋒圖像特征進行相應閾值的調(diào)試;另外本文設計的算法僅利用了雷達反射率因子圖像,陣風鋒在速度圖像上還具有輻合特征,但是雷達僅能探測到徑向速度,且算法對徑向速度圖像質(zhì)量有較高要求,導致在目前陣風鋒圖像特征識別中雷達徑向速度圖像數(shù)據(jù)難以利用。后期作者將嘗試結合天氣雷達徑向速度圖像特征以及地面自動站氣象要素變化特征,以進一步提升基于陣風鋒的大風識別算法。

    猜你喜歡
    陣風反射率分值
    影響Mini LED板油墨層反射率的因素
    近岸水體異源遙感反射率產(chǎn)品的融合方法研究
    海洋通報(2022年4期)2022-10-10 07:40:26
    一起來看看交通違法記分分值有什么變化
    工會博覽(2022年8期)2022-06-30 12:19:30
    具有顏色恒常性的光譜反射率重建
    陣風戰(zhàn)斗機
    軍事文摘(2021年19期)2021-10-10 13:29:36
    法國陣風戰(zhàn)斗機
    軍事文摘(2021年17期)2021-09-24 01:11:46
    陣風勁吹
    航空世界(2018年12期)2018-07-16 08:34:50
    化學腐蝕硅表面結構反射率影響因素的研究*
    電子器件(2017年2期)2017-04-25 08:58:37
    夾在兩陣風之間
    宿遷城鎮(zhèn)居民醫(yī)保按病種分值結算初探
    亚洲国产看品久久| 中文字幕精品免费在线观看视频| 欧美黄色片欧美黄色片| 国产男靠女视频免费网站| 欧美日本视频| 久久久久久久久久久久大奶| 日本五十路高清| 国产亚洲精品av在线| 最近最新中文字幕大全电影3 | 一级毛片精品| 人妻丰满熟妇av一区二区三区| 男人操女人黄网站| 国产主播在线观看一区二区| 9热在线视频观看99| 国产精品野战在线观看| 国内久久婷婷六月综合欲色啪| 免费不卡黄色视频| 亚洲人成电影免费在线| 视频在线观看一区二区三区| 久久久水蜜桃国产精品网| 日韩精品免费视频一区二区三区| 好男人电影高清在线观看| 亚洲av片天天在线观看| 淫妇啪啪啪对白视频| 91大片在线观看| 久99久视频精品免费| 日韩av在线大香蕉| 久久久久久久精品吃奶| 亚洲国产精品成人综合色| 亚洲五月婷婷丁香| 精品一区二区三区四区五区乱码| 欧美另类亚洲清纯唯美| 狂野欧美激情性xxxx| 亚洲avbb在线观看| 可以在线观看毛片的网站| 久久精品亚洲精品国产色婷小说| 又黄又爽又免费观看的视频| 老汉色∧v一级毛片| 精品电影一区二区在线| 身体一侧抽搐| 身体一侧抽搐| 好男人电影高清在线观看| 亚洲精品在线观看二区| 久久人人97超碰香蕉20202| 久久婷婷成人综合色麻豆| 亚洲精品中文字幕在线视频| 动漫黄色视频在线观看| 亚洲电影在线观看av| 亚洲精品在线观看二区| 男人舔女人下体高潮全视频| 国产精品九九99| 啪啪无遮挡十八禁网站| 在线观看免费午夜福利视频| 操美女的视频在线观看| 制服人妻中文乱码| 国产人伦9x9x在线观看| 日日摸夜夜添夜夜添小说| 国产成+人综合+亚洲专区| 九色亚洲精品在线播放| 天天躁夜夜躁狠狠躁躁| 亚洲专区中文字幕在线| 日韩成人在线观看一区二区三区| 黄色女人牲交| 一二三四社区在线视频社区8| 国产一区二区激情短视频| 免费观看人在逋| 国产高清有码在线观看视频 | 精品国产乱子伦一区二区三区| 精品卡一卡二卡四卡免费| 美女国产高潮福利片在线看| 亚洲成人免费电影在线观看| 神马国产精品三级电影在线观看 | 热re99久久国产66热| 成人国产综合亚洲| 人人妻人人爽人人添夜夜欢视频| 真人一进一出gif抽搐免费| 人妻丰满熟妇av一区二区三区| 999精品在线视频| 亚洲精品粉嫩美女一区| 久久久精品国产亚洲av高清涩受| 天天一区二区日本电影三级 | 伦理电影免费视频| 亚洲成人久久性| 色综合站精品国产| 欧美黄色淫秽网站| 亚洲专区国产一区二区| 黄片小视频在线播放| 日韩大码丰满熟妇| 丁香六月欧美| 日韩成人在线观看一区二区三区| 亚洲精华国产精华精| 国产片内射在线| 黑人巨大精品欧美一区二区mp4| 99热只有精品国产| 可以在线观看的亚洲视频| 啦啦啦韩国在线观看视频| 纯流量卡能插随身wifi吗| 色综合站精品国产| 欧美色欧美亚洲另类二区 | 中国美女看黄片| 久久午夜综合久久蜜桃| 午夜视频精品福利| 国产伦人伦偷精品视频| 久久影院123| 久久久久久免费高清国产稀缺| 欧美黑人欧美精品刺激| 欧美色欧美亚洲另类二区 | 两性夫妻黄色片| 怎么达到女性高潮| 午夜影院日韩av| 大型av网站在线播放| 国产av一区在线观看免费| 一边摸一边做爽爽视频免费| 国产一区二区激情短视频| 99久久精品国产亚洲精品| 国产精品乱码一区二三区的特点 | 精品国内亚洲2022精品成人| 日韩欧美国产一区二区入口| 久久久久久久精品吃奶| 我的亚洲天堂| 成人av一区二区三区在线看| 男人舔女人下体高潮全视频| 日韩欧美一区二区三区在线观看| 欧美日本中文国产一区发布| av超薄肉色丝袜交足视频| aaaaa片日本免费| 看免费av毛片| 成人国产综合亚洲| 自线自在国产av| 在线观看免费视频网站a站| 国产一卡二卡三卡精品| 免费在线观看完整版高清| 久久久精品欧美日韩精品| 一夜夜www| 高清毛片免费观看视频网站| 别揉我奶头~嗯~啊~动态视频| av福利片在线| 欧美成人一区二区免费高清观看 | 天天添夜夜摸| 高清毛片免费观看视频网站| 两性夫妻黄色片| 99riav亚洲国产免费| 777久久人妻少妇嫩草av网站| 日韩精品中文字幕看吧| 麻豆成人av在线观看| 妹子高潮喷水视频| 999精品在线视频| 午夜福利高清视频| 男人舔女人的私密视频| 久久人人精品亚洲av| 99热只有精品国产| 久久中文字幕人妻熟女| 亚洲一码二码三码区别大吗| 人成视频在线观看免费观看| 精品免费久久久久久久清纯| 一边摸一边抽搐一进一小说| 在线观看日韩欧美| 国内久久婷婷六月综合欲色啪| 日韩大尺度精品在线看网址 | 久久精品aⅴ一区二区三区四区| 午夜a级毛片| 亚洲久久久国产精品| 国产av精品麻豆| 亚洲国产欧美一区二区综合| av福利片在线| 日本三级黄在线观看| 免费在线观看视频国产中文字幕亚洲| 国产99白浆流出| 亚洲自拍偷在线| 欧美激情极品国产一区二区三区| 色播亚洲综合网| 国产区一区二久久| 人人澡人人妻人| 窝窝影院91人妻| 好看av亚洲va欧美ⅴa在| 亚洲一区中文字幕在线| 91字幕亚洲| 国产主播在线观看一区二区| 亚洲熟妇中文字幕五十中出| 国产精品日韩av在线免费观看 | 日韩欧美在线二视频| 国产精品一区二区免费欧美| 欧美最黄视频在线播放免费| 老熟妇乱子伦视频在线观看| 此物有八面人人有两片| 国产午夜精品久久久久久| 亚洲五月天丁香| 男人操女人黄网站| 国产成人欧美在线观看| 大型av网站在线播放| 久久久精品欧美日韩精品| 国产精品 国内视频| 国产精华一区二区三区| aaaaa片日本免费| 亚洲第一欧美日韩一区二区三区| 久久伊人香网站| 精品国内亚洲2022精品成人| 日韩国内少妇激情av| 变态另类成人亚洲欧美熟女 | 日日干狠狠操夜夜爽| 在线观看免费午夜福利视频| 亚洲情色 制服丝袜| 亚洲午夜理论影院| 精品国产超薄肉色丝袜足j| 亚洲片人在线观看| 国产成年人精品一区二区| 亚洲性夜色夜夜综合| 电影成人av| 精品国产一区二区久久| 中文字幕最新亚洲高清| 欧美乱色亚洲激情| 亚洲精品在线美女| 色在线成人网| АⅤ资源中文在线天堂| 亚洲av熟女| 三级毛片av免费| 国产片内射在线| 国产精品1区2区在线观看.| 久久国产乱子伦精品免费另类| 视频区欧美日本亚洲| 香蕉国产在线看| e午夜精品久久久久久久| 成人18禁在线播放| 男男h啪啪无遮挡| 国产精品亚洲美女久久久| 在线观看www视频免费| 国产三级在线视频| 久久中文看片网| 夜夜爽天天搞| 久久久久久亚洲精品国产蜜桃av| 看黄色毛片网站| 动漫黄色视频在线观看| 18禁国产床啪视频网站| 如日韩欧美国产精品一区二区三区| 日本在线视频免费播放| 久久精品亚洲熟妇少妇任你| 日本五十路高清| 天堂√8在线中文| 日韩av在线大香蕉| 国产亚洲欧美98| 黄色视频,在线免费观看| 久久人妻福利社区极品人妻图片| 欧美日韩福利视频一区二区| 亚洲情色 制服丝袜| 一级a爱视频在线免费观看| 巨乳人妻的诱惑在线观看| av电影中文网址| 亚洲精品一卡2卡三卡4卡5卡| 成人三级黄色视频| 99在线人妻在线中文字幕| 电影成人av| 精品午夜福利视频在线观看一区| 欧美一级a爱片免费观看看 | 亚洲,欧美精品.| 午夜亚洲福利在线播放| 操出白浆在线播放| 午夜久久久久精精品| av电影中文网址| 无遮挡黄片免费观看| 久久久久久国产a免费观看| 香蕉国产在线看| 法律面前人人平等表现在哪些方面| 久久久久国产精品人妻aⅴ院| 伊人久久大香线蕉亚洲五| 在线天堂中文资源库| 中文字幕人妻熟女乱码| 免费一级毛片在线播放高清视频 | 欧美乱妇无乱码| 亚洲精品国产色婷婷电影| 少妇熟女aⅴ在线视频| 久久欧美精品欧美久久欧美| 国产私拍福利视频在线观看| 亚洲av熟女| www.精华液| 亚洲精品一卡2卡三卡4卡5卡| 国产精品久久久人人做人人爽| 久久人妻福利社区极品人妻图片| 在线观看免费日韩欧美大片| 国产成人av激情在线播放| 老司机深夜福利视频在线观看| 老熟妇仑乱视频hdxx| 国产午夜精品久久久久久| 很黄的视频免费| 超碰成人久久| 两个人免费观看高清视频| 国产成人精品无人区| 丰满人妻熟妇乱又伦精品不卡| 久久久久久免费高清国产稀缺| 亚洲av成人av| 久久人妻熟女aⅴ| 性少妇av在线| 国产99久久九九免费精品| 亚洲人成网站在线播放欧美日韩| 亚洲片人在线观看| 久久精品91无色码中文字幕| 丰满人妻熟妇乱又伦精品不卡| 动漫黄色视频在线观看| 亚洲av日韩精品久久久久久密| 亚洲成国产人片在线观看| 中文字幕精品免费在线观看视频| av欧美777| 一区二区三区精品91| 精品第一国产精品| 亚洲精品在线美女| 久久精品人人爽人人爽视色| 久久久国产成人精品二区| 一级a爱视频在线免费观看| 色综合亚洲欧美另类图片| 国产av又大| 中文字幕高清在线视频| 在线视频色国产色| 91av网站免费观看| 久久久久国产一级毛片高清牌| 一级a爱视频在线免费观看| АⅤ资源中文在线天堂| 亚洲成a人片在线一区二区| 国产伦一二天堂av在线观看| 国产真人三级小视频在线观看| 嫁个100分男人电影在线观看| 国产精品野战在线观看| 免费在线观看日本一区| 精品一区二区三区视频在线观看免费| 岛国视频午夜一区免费看| 午夜精品久久久久久毛片777| 咕卡用的链子| 天天躁狠狠躁夜夜躁狠狠躁| 日本 欧美在线| 欧美 亚洲 国产 日韩一| 女人高潮潮喷娇喘18禁视频| av视频在线观看入口| 成熟少妇高潮喷水视频| 麻豆成人av在线观看| 国产精品,欧美在线| 国产精品久久视频播放| 国产高清videossex| 手机成人av网站| 久久精品影院6| 一级作爱视频免费观看| 午夜福利影视在线免费观看| 国产精品综合久久久久久久免费 | 51午夜福利影视在线观看| 成人国语在线视频| 变态另类成人亚洲欧美熟女 | 中文字幕最新亚洲高清| 此物有八面人人有两片| 亚洲国产精品成人综合色| 久久热在线av| 亚洲一区二区三区色噜噜| 在线天堂中文资源库| 久久草成人影院| 亚洲中文字幕日韩| 啦啦啦 在线观看视频| 久久草成人影院| 亚洲欧美精品综合一区二区三区| АⅤ资源中文在线天堂| 欧美色欧美亚洲另类二区 | 露出奶头的视频| 欧美中文综合在线视频| 夜夜看夜夜爽夜夜摸| 亚洲av日韩精品久久久久久密| 国产单亲对白刺激| 一区二区三区激情视频| 久久久久久大精品| 成人亚洲精品av一区二区| 欧美亚洲日本最大视频资源| 免费高清在线观看日韩| 日韩欧美在线二视频| 久久久久久人人人人人| 少妇裸体淫交视频免费看高清 | 中文字幕另类日韩欧美亚洲嫩草| 美国免费a级毛片| 黑丝袜美女国产一区| 亚洲国产日韩欧美精品在线观看 | 国产欧美日韩综合在线一区二区| 此物有八面人人有两片| 免费不卡黄色视频| 97人妻精品一区二区三区麻豆 | 9色porny在线观看| АⅤ资源中文在线天堂| 成人欧美大片| 美女高潮到喷水免费观看| 久久 成人 亚洲| 美女高潮到喷水免费观看| 久久人人爽av亚洲精品天堂| 丁香欧美五月| 久久九九热精品免费| 亚洲午夜理论影院| 久久婷婷人人爽人人干人人爱 | 色哟哟哟哟哟哟| 国产一区二区三区在线臀色熟女| √禁漫天堂资源中文www| 不卡一级毛片| 国产av精品麻豆| 9色porny在线观看| 人人妻人人澡欧美一区二区 | 成人永久免费在线观看视频| 国产精品免费一区二区三区在线| 免费不卡黄色视频| 亚洲欧美激情在线| 国产成人精品在线电影| 国产av一区在线观看免费| 精品电影一区二区在线| 国产精品 国内视频| 成人欧美大片| 免费少妇av软件| av福利片在线| 国产免费男女视频| 国产av一区二区精品久久| 中文字幕人妻熟女乱码| 中文亚洲av片在线观看爽| 欧美日韩亚洲综合一区二区三区_| 级片在线观看| 韩国精品一区二区三区| 一级,二级,三级黄色视频| 中文字幕最新亚洲高清| 欧美日本亚洲视频在线播放| 久久性视频一级片| 天天一区二区日本电影三级 | 熟女少妇亚洲综合色aaa.| 麻豆久久精品国产亚洲av| 免费观看人在逋| 校园春色视频在线观看| 男人舔女人的私密视频| 国内毛片毛片毛片毛片毛片| 日本三级黄在线观看| 午夜视频精品福利| 国产亚洲精品久久久久5区| 国产人伦9x9x在线观看| www.自偷自拍.com| 一级,二级,三级黄色视频| 成人欧美大片| 91成人精品电影| 免费看十八禁软件| 亚洲专区字幕在线| 日本精品一区二区三区蜜桃| 国产激情久久老熟女| 妹子高潮喷水视频| 免费女性裸体啪啪无遮挡网站| 国产日韩一区二区三区精品不卡| 美女免费视频网站| 一进一出抽搐动态| 亚洲 欧美 日韩 在线 免费| 亚洲一码二码三码区别大吗| 久久 成人 亚洲| 国产一区在线观看成人免费| 两人在一起打扑克的视频| 不卡一级毛片| 在线观看午夜福利视频| 老司机深夜福利视频在线观看| 国产亚洲欧美在线一区二区| 久久影院123| 美女午夜性视频免费| 免费在线观看亚洲国产| 日本黄色视频三级网站网址| 色播亚洲综合网| 国产免费男女视频| 淫秽高清视频在线观看| 我的亚洲天堂| 成在线人永久免费视频| 日日爽夜夜爽网站| 日韩视频一区二区在线观看| 成年版毛片免费区| 国产色视频综合| 黑人巨大精品欧美一区二区蜜桃| 精品日产1卡2卡| 国产精品,欧美在线| 久久 成人 亚洲| 性色av乱码一区二区三区2| 少妇 在线观看| 好男人在线观看高清免费视频 | 91精品三级在线观看| 精品高清国产在线一区| 亚洲aⅴ乱码一区二区在线播放 | 一级黄色大片毛片| 国产精品美女特级片免费视频播放器 | 19禁男女啪啪无遮挡网站| 好看av亚洲va欧美ⅴa在| 一个人观看的视频www高清免费观看 | 日韩欧美三级三区| 国产亚洲精品综合一区在线观看 | 日韩视频一区二区在线观看| 国语自产精品视频在线第100页| 制服诱惑二区| 精品久久久久久,| 正在播放国产对白刺激| 成人亚洲精品一区在线观看| 曰老女人黄片| 别揉我奶头~嗯~啊~动态视频| 无遮挡黄片免费观看| 级片在线观看| or卡值多少钱| 97人妻天天添夜夜摸| 麻豆国产av国片精品| 精品久久久久久,| 午夜久久久久精精品| 视频区欧美日本亚洲| 亚洲一区中文字幕在线| 久久香蕉激情| 亚洲精品一区av在线观看| 久久伊人香网站| 久久香蕉国产精品| 色老头精品视频在线观看| 国产av又大| 久久国产乱子伦精品免费另类| 正在播放国产对白刺激| 嫩草影院精品99| 在线观看66精品国产| 国产麻豆成人av免费视频| 一a级毛片在线观看| 中文字幕久久专区| 免费观看人在逋| 日韩中文字幕欧美一区二区| av片东京热男人的天堂| 久久中文看片网| 91在线观看av| 中国美女看黄片| www.自偷自拍.com| 高清在线国产一区| 嫩草影院精品99| 麻豆成人av在线观看| 国产成人一区二区三区免费视频网站| 两个人免费观看高清视频| 99精品在免费线老司机午夜| 99国产精品一区二区蜜桃av| 少妇粗大呻吟视频| 熟女少妇亚洲综合色aaa.| 啦啦啦 在线观看视频| 少妇熟女aⅴ在线视频| 亚洲色图 男人天堂 中文字幕| 天天躁狠狠躁夜夜躁狠狠躁| 50天的宝宝边吃奶边哭怎么回事| 亚洲av第一区精品v没综合| 精品一区二区三区视频在线观看免费| 久久久国产成人精品二区| 亚洲男人天堂网一区| 欧美乱色亚洲激情| 午夜精品在线福利| 日本欧美视频一区| 窝窝影院91人妻| 夜夜爽天天搞| 日韩精品青青久久久久久| 欧美国产日韩亚洲一区| 久久精品人人爽人人爽视色| 麻豆成人av在线观看| 日本vs欧美在线观看视频| 国产91精品成人一区二区三区| 精品久久久精品久久久| 国产色视频综合| 一级a爱片免费观看的视频| 国产在线精品亚洲第一网站| 中文字幕另类日韩欧美亚洲嫩草| 身体一侧抽搐| 美女扒开内裤让男人捅视频| 久久精品人人爽人人爽视色| 黑人巨大精品欧美一区二区mp4| 99久久精品国产亚洲精品| 国产私拍福利视频在线观看| 久久精品亚洲精品国产色婷小说| 亚洲色图 男人天堂 中文字幕| 美女午夜性视频免费| 丁香六月欧美| 18美女黄网站色大片免费观看| 亚洲精品在线观看二区| 久久热在线av| 日日干狠狠操夜夜爽| 纯流量卡能插随身wifi吗| 久久伊人香网站| 大陆偷拍与自拍| 夜夜看夜夜爽夜夜摸| 亚洲免费av在线视频| av欧美777| 亚洲情色 制服丝袜| 国产亚洲欧美在线一区二区| 搡老熟女国产l中国老女人| 国产伦人伦偷精品视频| 亚洲一区中文字幕在线| 中亚洲国语对白在线视频| 欧美成人午夜精品| 给我免费播放毛片高清在线观看| 妹子高潮喷水视频| 一本久久中文字幕| 99久久久亚洲精品蜜臀av| 国产麻豆69| 国产精品久久久av美女十八| 色尼玛亚洲综合影院| 日韩欧美三级三区| 国产精品野战在线观看| 日本免费a在线| 少妇熟女aⅴ在线视频| 日韩高清综合在线| 麻豆久久精品国产亚洲av| 两人在一起打扑克的视频| 99re在线观看精品视频| 亚洲精品国产一区二区精华液| 黄色女人牲交| 精品久久久久久久久久免费视频| 黄色视频不卡| 色播在线永久视频| 91成年电影在线观看| 在线观看舔阴道视频| 精品免费久久久久久久清纯| 国产免费男女视频| 一个人观看的视频www高清免费观看 | 黄片小视频在线播放| 国产人伦9x9x在线观看| 天堂动漫精品| 黑人操中国人逼视频| 国内毛片毛片毛片毛片毛片| 中文字幕精品免费在线观看视频| 黄频高清免费视频| 岛国在线观看网站| 视频在线观看一区二区三区| 国产一级毛片七仙女欲春2 | 97超级碰碰碰精品色视频在线观看| avwww免费|