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

    日冕物質(zhì)拋射檢測(cè)研究進(jìn)展?

    2020-01-16 08:02:08尚振宏強(qiáng)振平黎敬濤冼祥貴楊志鵬
    天文研究與技術(shù) 2020年1期
    關(guān)鍵詞:日冕分類器灰度

    郭 敏,尚振宏,強(qiáng)振平,黎敬濤,冼祥貴,楊志鵬

    (1.昆明理工大學(xué)信息工程與自動(dòng)化學(xué)院,云南昆明 650500;2.西南林業(yè)大學(xué)大數(shù)據(jù)與智能工程學(xué)院,云南昆明 650224)

    日冕物質(zhì)拋射(Coronal Mass Ejection,CME)是等離子體從太陽(yáng)的日冕層噴射到行星際空間的一種動(dòng)態(tài)事件。它在20世紀(jì)70年代早期的天空實(shí)驗(yàn)室(Skylab)任務(wù)中被首次觀測(cè)到,是源自太陽(yáng)的最大和最具活力的爆發(fā)現(xiàn)象,并且可以通過(guò)白光日冕儀在延伸的日冕中觀察到[1]。由于日冕物質(zhì)拋射是一種巨大的磁化等離子體噴發(fā),因此,它能以高達(dá)每秒數(shù)千千米的速度傳播[2]。鑒于這種爆發(fā)活動(dòng)引發(fā)的太陽(yáng)高能粒子事件和地磁暴會(huì)影響航空安全、衛(wèi)星運(yùn)行、通信系統(tǒng)、電力設(shè)施等,因此有必要對(duì)日冕物質(zhì)拋射進(jìn)行檢測(cè),為空間天氣預(yù)報(bào)提供預(yù)警信號(hào),從而采取相應(yīng)的措施以避免不必要的損失[3]。此外,有關(guān)日冕物質(zhì)拋射的統(tǒng)計(jì)信息對(duì)于更好地了解其性質(zhì)也非常重要。日冕物質(zhì)拋射的概念一直比較模糊,這種定義的不明確性導(dǎo)致不同檢測(cè)方法得到的日冕物質(zhì)拋射特征和數(shù)量有所不同。文[4]將日冕物質(zhì)拋射定義為發(fā)生在幾分鐘到幾小時(shí)的時(shí)間尺度上,在日冕視場(chǎng)中出現(xiàn)的新的、離散的、明亮的白光特征。由于日冕物質(zhì)拋射現(xiàn)象與其他太陽(yáng)活動(dòng)有相似的特征,為了對(duì)日冕物質(zhì)拋射的認(rèn)識(shí)更加準(zhǔn)確,文[5]將日冕物質(zhì)拋射的定義修改為:在日冕視場(chǎng)中出現(xiàn)的一種新的,離散的,明亮的白光特征,并具有徑向向外的速度。圖1是空間環(huán)境預(yù)報(bào)中心的SOHO/LASCO C2衛(wèi)星觀測(cè)到的一次日冕物質(zhì)拋射的爆發(fā)過(guò)程。

    隨著SOHO,Wind和STEREO等衛(wèi)星的發(fā)射,日冕物質(zhì)拋射的檢測(cè)和識(shí)別成為可能。其中1995年發(fā)射的SOHO衛(wèi)星搭載的LASCO日冕儀可以觀測(cè)到日冕物質(zhì)拋射的發(fā)生過(guò)程,得到日冕物質(zhì)拋射的影像資料,通過(guò)這些資料可以檢測(cè)日冕物質(zhì)拋射的運(yùn)動(dòng)狀態(tài),進(jìn)而得到日冕物質(zhì)拋射的描述參數(shù)(位置角、角寬度、速度等)。目前,國(guó)內(nèi)外關(guān)于日冕物質(zhì)拋射的檢測(cè)研究已經(jīng)取得很大進(jìn)展,形成了許多日冕物質(zhì)拋射目錄,為以后的進(jìn)一步研究提供了基準(zhǔn)和參考。

    為了方便讀者詳細(xì)了解目前日冕物質(zhì)拋射檢測(cè)的研究進(jìn)展,并為相應(yīng)研究提供借鑒,本文梳理了國(guó)內(nèi)外關(guān)于日冕物質(zhì)拋射檢測(cè)的一些方法,給出了具有代表性方法的具體識(shí)別過(guò)程。同時(shí)為了便于讀者理解,本文將目前日冕物質(zhì)拋射檢測(cè)方法分成了兩大類:基于手工識(shí)別的方法和基于自動(dòng)檢測(cè)的方法,同時(shí)又根據(jù)所選取的特征和采用的技術(shù)將自動(dòng)檢測(cè)方法分為4類:基于灰度特征、基于紋理特征、基于光流法和基于學(xué)習(xí)的方法。

    1 日冕物質(zhì)拋射的特征

    日冕物質(zhì)拋射檢測(cè)實(shí)質(zhì)上是利用人工標(biāo)記或計(jì)算機(jī)自動(dòng)識(shí)別的方法檢測(cè)日冕物質(zhì)拋射事件,從而得到日冕物質(zhì)拋射的特征。常用的日冕物質(zhì)拋射的特征包括速度、張角、發(fā)生規(guī)律等。

    1.1 日冕物質(zhì)拋射的速度

    日冕物質(zhì)拋射是一種運(yùn)動(dòng)的過(guò)程,其速度是一個(gè)重要的物理特性。由于目前觀測(cè)到的日冕物質(zhì)拋射通常是在天空平面的投影,因此一般情況下,日冕物質(zhì)拋射的速度是指日冕物質(zhì)拋射最快前沿的運(yùn)動(dòng)速度[6],它由日冕物質(zhì)拋射前沿的“時(shí)間-高度”數(shù)據(jù)經(jīng)過(guò)線性擬合得到。圖2是協(xié)調(diào)數(shù)據(jù)分析研討會(huì)中心(Coordinated Data Analysis Workshop,CDAW①CDAW: https://cdaw.gsfc.nasa.gov/CME_list/)目錄中2017年7月2日13時(shí)25分42秒發(fā)生的一個(gè)日冕物質(zhì)拋射的時(shí)間-高度圖。日冕物質(zhì)拋射的速度范圍為10~3 500 km/s,一般情況下,速度>1 000 km/s為快速日冕物質(zhì)拋射,速度<100 km/s為緩慢日冕物質(zhì)拋射[7]。

    1.2 日冕物質(zhì)拋射的張角

    日冕物質(zhì)拋射的大小一般用張角(角寬度)衡量,張角是日冕物質(zhì)拋射區(qū)域兩個(gè)邊緣所構(gòu)成的角度,在數(shù)值上等于兩邊緣的位置角之差。圖3是CACTus軟件檢測(cè)到的一個(gè)日冕物質(zhì)拋射,其中a,b兩處為某時(shí)刻發(fā)生的日冕物質(zhì)拋射區(qū)域的位置角,張角則為圖中白色實(shí)線中間的角度。日冕物質(zhì)拋射的張角從幾度到360°,360°的日冕物質(zhì)拋射稱為暈狀日冕物質(zhì)拋射,一般情況下,寬度>120°為寬日冕物質(zhì)拋射,寬度<30°為窄日冕物質(zhì)拋射[7]。

    1.3 日冕物質(zhì)拋射與太陽(yáng)黑子數(shù)的關(guān)系

    日冕物質(zhì)拋射的發(fā)生規(guī)律與太陽(yáng)黑子活動(dòng)周有關(guān),在太陽(yáng)活動(dòng)極小年時(shí)為0.5個(gè)/天,在太陽(yáng)活動(dòng)極大年時(shí)為2~6個(gè)/天。圖4為2008年到2017年間協(xié)調(diào)數(shù)據(jù)分析研討會(huì)目錄記錄的日冕物質(zhì)拋射的數(shù)量,圖5中的數(shù)據(jù)來(lái)源于太陽(yáng)黑子指數(shù)和長(zhǎng)期太陽(yáng)觀測(cè)(SILSO②SILSO: http://sidc.oma.be/silso/home)中近13年太陽(yáng)黑子的數(shù)量,可以看到,在這兩幅圖的重疊時(shí)間段(2008年~2017年),日冕物質(zhì)拋射的發(fā)生頻率和黑子的基本一致,因此,可以說(shuō)明日冕物質(zhì)拋射的發(fā)生規(guī)律與太陽(yáng)黑子活動(dòng)周有關(guān)。

    2 日冕物質(zhì)拋射檢測(cè)方法

    十多年來(lái),研究者們利用多種方法檢測(cè)到日冕物質(zhì)拋射事件及其特征,并且制成了日冕物質(zhì)拋射目錄供科研人員學(xué)習(xí)研究。目前日冕物質(zhì)拋射的檢測(cè)方法主要分為兩大類:

    (1)基于手工方法標(biāo)識(shí)日冕物質(zhì)拋射:主要有兩個(gè)著名的目錄:CDAW目錄和海軍研究實(shí)驗(yàn)室(Naval Research Laboratory,NRL)目錄,這種方法主要靠觀測(cè)者每天手動(dòng)記錄數(shù)據(jù)進(jìn)行日冕物質(zhì)拋射的編目,比較耗時(shí);

    (2)基于自動(dòng)方法檢測(cè)日冕物質(zhì)拋射:由于日冕物質(zhì)拋射是一種動(dòng)態(tài)變化的過(guò)程,通常具有明亮的強(qiáng)度特征和復(fù)雜的紋理結(jié)構(gòu)。依據(jù)日冕物質(zhì)拋射檢測(cè)選用的特征及技術(shù)手段將現(xiàn)有的自動(dòng)檢測(cè)方法分為4類進(jìn)行概述和分析,表1是對(duì)日冕物質(zhì)拋射檢測(cè)方法的總結(jié)。

    圖2 CDAW目錄中2017年7月2日13時(shí)25分42秒發(fā)生的一個(gè)日冕物質(zhì)拋射的時(shí)間-高度圖Fig.2 The time-height plot of a CME occurring at 13:25:42 on July 2,2017 in the CDAW catalog

    圖3 日冕物質(zhì)拋射的角寬度Fig.3 The angular width of a CME

    圖4 2008年到2017年間CDAW目錄記錄的日冕物質(zhì)拋射的數(shù)量Fig.4 Number of CMEs recorded between 2008 and 2017 in the CDAW catalog

    圖5 近13年太陽(yáng)黑子數(shù)量的變化情況Fig.5 Evolution of sunspot numbers in the past 13 years

    2.1 手工識(shí)別目錄與自動(dòng)檢測(cè)目錄的比較

    日冕物質(zhì)拋射與耀斑、太陽(yáng)高能粒子、地磁風(fēng)暴等許多現(xiàn)象有關(guān),因此檢測(cè)日冕物質(zhì)拋射現(xiàn)象并編制事件目錄對(duì)于國(guó)內(nèi)外研究工作者理解這些現(xiàn)象并開展相關(guān)工作非常重要。日冕物質(zhì)拋射目錄是將檢測(cè)到的日冕物質(zhì)拋射事件特征以數(shù)據(jù)和圖表的形式記錄下來(lái),相關(guān)參數(shù)特征主要包括速度、角寬度、位置角等。根據(jù)上述介紹的日冕物質(zhì)拋射檢測(cè)方法可以得到兩類日冕物質(zhì)拋射目錄:手工識(shí)別目錄和自動(dòng)檢測(cè)目錄,通過(guò)以CDAW目錄和CACTus目錄為代表進(jìn)行特征(速度v和角寬度W)的比較[7]得到這兩類目錄的特點(diǎn):

    (1)識(shí)別狹窄的日冕物質(zhì)拋射(W<30°):研究表明,在同一時(shí)期CACTus檢測(cè)到的窄日冕物質(zhì)拋射遠(yuǎn)遠(yuǎn)多于CDAW,主要原因在于CDAW目錄在太陽(yáng)活動(dòng)極大年遺漏了許多窄日冕物質(zhì)拋射。

    (2)識(shí)別快速的日冕物質(zhì)拋射(v>1000km/s):在同一時(shí)期雖然CACTus檢測(cè)到的快速日冕物質(zhì)拋射多于CDAW,但是在檢測(cè)到的32%的真正日冕物質(zhì)拋射中只有6%是快速日冕物質(zhì)拋射。

    因此可以得出結(jié)論:手工目錄和自動(dòng)檢測(cè)目錄各有優(yōu)缺點(diǎn),由于手工目錄是基于人工標(biāo)記的,因此不可避免地遺漏許多狹窄的日冕物質(zhì)拋射;而自動(dòng)檢測(cè)目錄雖然比手工目錄識(shí)別的快速日冕物質(zhì)拋射多,但是準(zhǔn)確率不高。

    2.2 手工方法標(biāo)識(shí)日冕物質(zhì)拋射

    1971年發(fā)現(xiàn)日冕物質(zhì)拋射現(xiàn)象以來(lái),已經(jīng)被多種星載儀器觀測(cè)到。自1995年SOHO航天器發(fā)射,日冕物質(zhì)拋射觀測(cè)主要采用LASCO日冕儀。LASCO CME的識(shí)別和編目是一項(xiàng)重要任務(wù),為進(jìn)一步的科學(xué)研究提供基礎(chǔ)知識(shí)。

    NRL目錄由LASCO觀測(cè)員編制,他們查看LASCO日冕圖像的序列并記錄每天發(fā)生的事件。這是一份初步目錄,提供有關(guān)日冕物質(zhì)拋射時(shí)間和大致的位置信息。

    CDAW目錄是由美國(guó)宇航局、美國(guó)天主教大學(xué)以及海軍研究實(shí)驗(yàn)室合作,在CDAW數(shù)據(jù)中心生成和維護(hù)。標(biāo)記的每個(gè)日冕物質(zhì)拋射由發(fā)生的日期和時(shí)刻確定,該目錄提供了日冕物質(zhì)拋射屬性的一些測(cè)量值,包括速度和角寬度等,同時(shí)還提供了確定日冕物質(zhì)拋射速度的“時(shí)間-高度”圖。除了得到的數(shù)據(jù)和圖表之外,該目錄還包含由日冕儀得到的日冕物質(zhì)拋射圖像序列,可供用戶下載查看。

    由于NRL目錄和CDAW目錄都是觀測(cè)者手工標(biāo)識(shí)的,即使對(duì)于同一觀測(cè)者,識(shí)別能力也不是恒定的。因此,所提供的事件和測(cè)量的參數(shù)受人類主觀因素影響較大,同時(shí)編目過(guò)程也非常耗時(shí)。

    表1 日冕物質(zhì)拋射檢測(cè)方法總結(jié)Table 1 Summary of CME detection methods

    2.3 自動(dòng)方法檢測(cè)日冕物質(zhì)拋射

    傳統(tǒng)的日冕物質(zhì)拋射檢測(cè)方法基于人工觀測(cè),這種方法效率低,容易受到個(gè)體主觀因素的影響,隨著近年來(lái)自動(dòng)識(shí)別技術(shù)的迅速發(fā)展,出現(xiàn)了很多自動(dòng)檢測(cè)日冕物質(zhì)拋射的算法,為日冕物質(zhì)拋射的檢測(cè)研究提供了新方法。由于日冕物質(zhì)拋射具有不同的形態(tài)、尺度和特征,通常表現(xiàn)為一個(gè)明亮、紋理復(fù)雜的增強(qiáng)結(jié)構(gòu),尾隨著一個(gè)亮度不足的暗區(qū)域[3],可以將現(xiàn)有的自動(dòng)檢測(cè)方法分為4類。

    2.3.1 基于灰度特征的方法

    灰度特征是基于像素點(diǎn)表示圖像亮度的特征,與背景灰度特征不同,日冕物質(zhì)拋射具有更明亮的結(jié)構(gòu),通常表現(xiàn)為明亮的白光特征,因此通過(guò)提取灰度特征檢測(cè)日冕物質(zhì)拋射是一種很直觀的思路。目前,基于灰度特征檢測(cè)日冕物質(zhì)拋射的方法研究較為成熟,國(guó)內(nèi)外已經(jīng)有不少人利用該方法進(jìn)行日冕物質(zhì)拋射檢測(cè)并且取得了良好的效果。計(jì)算機(jī)輔助跟蹤軟件包(CACTus)和太陽(yáng)能爆發(fā)事件檢測(cè)系統(tǒng)(SEEDS)是其中具有代表性的兩種檢測(cè)方法。

    CACTus方法由文[5,8]提出,該方法通過(guò)濾波去噪、極坐標(biāo)變換、幀間差分以及整合LASCO C2和C3圖像完成預(yù)處理。對(duì)不同時(shí)刻t的差分圖像[θ,r]進(jìn)行整合,得到[t,θ,r]數(shù)據(jù)立方體,該立方體沿指定角度θ的切片[t,r]中的一條脊線是一個(gè)日冕物質(zhì)拋射沿指定角度θ上的噴射表現(xiàn)。為了能清楚地看到脊線,將切片圖顏色進(jìn)行反轉(zhuǎn)得到傾斜的暗黑脊線,如圖6。由于霍夫變換[26]是一種圖像處理技術(shù),它可以從圖像中檢測(cè)出線段,因此霍夫變換被用來(lái)檢測(cè)[t,r]切片中的脊線。同時(shí)霍夫變換檢測(cè)到的脊線是沿一個(gè)角度切片的日冕物質(zhì)拋射,為了得到完整的日冕物質(zhì)拋射,將圖像投影到[v,θ,t]數(shù)據(jù)立方體中,速度v可通過(guò)脊線段長(zhǎng)度計(jì)算得到,此時(shí)日冕物質(zhì)拋射識(shí)別便轉(zhuǎn)換為[v,θ,t]立方體中點(diǎn)的聚類問(wèn)題。由于CACTus認(rèn)為同一個(gè)日冕物質(zhì)拋射沿不同θ方向的徑向速度相似,v方向信息對(duì)聚類貢獻(xiàn)較小,為了降低計(jì)算復(fù)雜度,對(duì)[v,θ,t]沿v方向積分,使得聚類轉(zhuǎn)換到二維空間[θ,t]中進(jìn)行,從而在日冕物質(zhì)拋射概覽圖中得到聚類的位置,該位置為日冕物質(zhì)拋射發(fā)生的起始時(shí)間和角度,垂直方向上聚類的長(zhǎng)度表示日冕物質(zhì)拋射噴發(fā)持續(xù)的時(shí)間。

    圖6 時(shí)間-高度圖中檢測(cè)脊線,圖像來(lái)源于BERGHMANS等(2004)[8]Fig.6 Detected ridges in a time-height diagram,image from BERGHMANS et al.(2004)[8]

    CACTus軟件第1個(gè)利用自動(dòng)方法實(shí)現(xiàn)日冕物質(zhì)拋射檢測(cè),相比于人工標(biāo)識(shí),該方法檢測(cè)速度快,檢測(cè)準(zhǔn)確率高,可以檢測(cè)到CDAW目錄中約75%的日冕物質(zhì)拋射,并且能夠檢測(cè)到CDAW目錄中沒有檢測(cè)到的弱日冕物質(zhì)拋射。目前,CACTus軟件編制了一個(gè)從1996年至今的在線目錄③CACTus: http://sidc.oma.be/cactus/。

    SEEDS方法由文[9-10]提出,檢測(cè)過(guò)程使用基于區(qū)域增長(zhǎng)算法的圖像分割技術(shù)。與CACTus類似,該方法通過(guò)歸一化輸入圖像、濾波去噪、極坐標(biāo)變換完成預(yù)處理。然后將極坐標(biāo)變換得到的[θ,r]圖像序列沿θ方向進(jìn)行投影得到信號(hào)強(qiáng)度曲線,如圖7(a)。接著通過(guò)閾值處理得到日冕物質(zhì)拋射的中心和最亮的部分,即核心角度,同時(shí)對(duì)核心角度利用區(qū)域增長(zhǎng)算法得到整個(gè)日冕物質(zhì)拋射的角寬度。日冕物質(zhì)拋射前沿的確定是通過(guò)將圖像沿r方向進(jìn)行投影得到信號(hào)強(qiáng)度曲線,如圖7(b),曲線峰值對(duì)應(yīng)于最大高度處,曲線峰值的一半在前沿和后繼處。通過(guò)對(duì)比連續(xù)兩幀圖像,如果后張圖像的前沿和后繼以及兩者之差高于前張圖像時(shí),表示日冕物質(zhì)拋射是擴(kuò)張的,同時(shí)若前沿溢出視場(chǎng)或者日冕物質(zhì)拋射太暗識(shí)別不出時(shí),停止追蹤,尋找下一個(gè)日冕物質(zhì)拋射。

    在檢測(cè)準(zhǔn)確率上,與CACTus類似,SEEDS能檢測(cè)到CDAW目錄約75%的日冕物質(zhì)拋射。但在檢測(cè)到的日冕物質(zhì)拋射總數(shù)量上,通過(guò)CACTus和SEEDS方法檢測(cè)得到的日冕物質(zhì)拋射數(shù)量是CDAW目錄的兩倍多。同時(shí)SEEDS方法還可以得到日冕物質(zhì)拋射的運(yùn)動(dòng)前沿輪廓,這可以在SEEDS目錄④SEEDS: http://spaceweather.gmu.edu/seeds/看到。

    圖7 (a)沿θ軸的一維投影強(qiáng)度分布圖;(b)沿r軸的一維投影強(qiáng)度分布圖。圖像來(lái)源于OLMEDO等(2008)[9]Fig.7 (a)The intensity profile of the 1D projection along the θ axis; (b)The intensity profile of the 1D projection along the r axis,image from OLMEDO et al.(2008)[9]

    除了上述詳細(xì)介紹的兩種方法之外,還有許多基于灰度特征檢測(cè)日冕物質(zhì)拋射的方法。例如文[11]提出了一種基于支持向量機(jī)(Support Vector Machine,SVM)的日冕物質(zhì)拋射檢測(cè)方法,為了將日冕物質(zhì)拋射與其他結(jié)構(gòu)區(qū)分開來(lái),該算法通過(guò)預(yù)處理產(chǎn)生差分圖像和分割圖像,然后再次利用閾值方法分割預(yù)處理得到的差分圖像和分割圖像,最終的分割結(jié)果是來(lái)自兩個(gè)分割結(jié)果的總和。接著將所有分割區(qū)域視為日冕物質(zhì)拋射候選區(qū)域,通過(guò)設(shè)定閾值將候選區(qū)域分為日冕物質(zhì)拋射區(qū)域和背景,最后在日冕物質(zhì)拋射區(qū)域利用支持向量機(jī)分類器區(qū)分強(qiáng)日冕物質(zhì)拋射和弱日冕物質(zhì)拋射。與其他方法相比,該方法的最大優(yōu)點(diǎn)是能夠?qū)?qiáng)日冕物質(zhì)拋射和大多數(shù)弱日冕物質(zhì)拋射進(jìn)行準(zhǔn)確而快速的檢測(cè)。

    與上述3種方法將LASCO圖像轉(zhuǎn)換為差分處理的極坐標(biāo)圖不同,文[12-13]提出的ARTEMIS方法將LASCO C2圖轉(zhuǎn)換為天氣圖。如圖8,天氣圖的橫坐標(biāo)為時(shí)間,表示日冕物質(zhì)拋射的流出時(shí)間,縱坐標(biāo)為緯度,表示日冕物質(zhì)拋射的角度寬度。該算法利用中值濾波器去除噪聲,然后對(duì)天氣圖進(jìn)行閾值處理返回一個(gè)二維掩碼,該掩碼定義了圖像的感興趣區(qū)域。最后通過(guò)引入日冕物質(zhì)拋射特性的先驗(yàn)知識(shí)正確識(shí)別日冕物質(zhì)拋射。算法最終得到一個(gè)ARTEMIS目錄,列出了每個(gè)檢測(cè)到的日冕物質(zhì)拋射事件以及它們的主要參數(shù)。與CDAW目錄比較可以發(fā)現(xiàn),該方法的優(yōu)點(diǎn)是能夠檢測(cè)到一些小的、微弱的日冕物質(zhì)拋射,與CACTus,SEEDS一樣,檢測(cè)到的日冕物質(zhì)拋射事件多于CDAW目錄。

    圖8 將日冕物質(zhì)拋射圖轉(zhuǎn)換為天氣圖,圖像來(lái)源于BOURSIER等(2009)[13]Fig.8 Converting a CME image to a synoptic map,image from BOURSIER et al.(2009)[13]

    基于灰度特征的4種方法先對(duì)LASCO圖像預(yù)處理得到極坐標(biāo)圖或天氣圖,再進(jìn)行閾值處理和形態(tài)學(xué)處理提取日冕物質(zhì)拋射特征,最后通過(guò)對(duì)日冕物質(zhì)拋射特征的識(shí)別或跟蹤計(jì)算出日冕物質(zhì)拋射的各種屬性。可以發(fā)現(xiàn)這些方法有一個(gè)共同點(diǎn),就是都采用了亮度增強(qiáng)的方式突出日冕圖像中的感興趣區(qū)域即日冕物質(zhì)拋射區(qū)域。同時(shí)上述幾種方法在同一時(shí)期的LASCO圖像上檢測(cè)出的日冕物質(zhì)拋射數(shù)量都比CDAW目錄多,但也都不能完全檢測(cè)出CDAW中100%的日冕物質(zhì)拋射,這與自動(dòng)檢測(cè)方法使用的圖像預(yù)處理技術(shù)和檢測(cè)規(guī)則、閾值的選取以及日冕物質(zhì)拋射的不明確定義都有關(guān)系。

    2.3.2 基于紋理特征的方法

    日冕物質(zhì)拋射不僅具有明顯的亮度結(jié)構(gòu),而且還有復(fù)雜的紋理結(jié)構(gòu)。對(duì)于一些比較暗但具有明顯紋理結(jié)構(gòu)的日冕物質(zhì)拋射區(qū)域來(lái)說(shuō),利用紋理特征檢測(cè)是一個(gè)不錯(cuò)的選擇。紋理特征的提取一般需要設(shè)定一定大小、包含多個(gè)像素點(diǎn)的窗口區(qū)域,然后從中取得紋理特征,尤其在檢索具有明顯粗細(xì)、疏密等圖像時(shí)有良好的表現(xiàn)。

    文[14]使用小波分析方法對(duì)日冕物質(zhì)拋射進(jìn)行分類,該算法是一種基于頻域分析的方法。通過(guò)在LASCO C2圖像上放置固定像素的窗口(實(shí)驗(yàn)數(shù)據(jù)表明,窗口可以放置在圖像的任意位置),使用標(biāo)準(zhǔn)快速傅里葉變換算法在窗口上執(zhí)行卷積,同時(shí)使用墨西哥帽作為基函數(shù)獲得窗口區(qū)域的小波譜,得到的小波譜圖橫坐標(biāo)為尺度,縱坐標(biāo)為通量。窗口區(qū)域表示的曲線都有一個(gè)拐點(diǎn),該點(diǎn)區(qū)分了通量隨尺度變化的速度,通過(guò)計(jì)算每條曲線的二階導(dǎo)數(shù)獲得拐點(diǎn)的尺度(ac)和通量(fc)。如圖9,在ac-fc圖中可以看到小波譜明顯被分為兩組:低通量、小空間尺度組(Homogeneous Group)和高通量、大空間尺度組(Collimated Group),其中box1和box2表示圖像上的兩個(gè)窗口。將該方法應(yīng)用到更大的日冕物質(zhì)拋射圖像樣本中,通過(guò)在ac-fc圖中繪制拐點(diǎn)的尺度和通量值,將這些點(diǎn)進(jìn)行線性最小二乘擬合,就可以區(qū)分日冕物質(zhì)拋射圖像屬于哪一組。該方法是一種圖像濾波的分析方法,利用頻域的高通和低通特征分類日冕物質(zhì)拋射圖像,可以對(duì)由不同衛(wèi)星獲得的日冕物質(zhì)拋射進(jìn)行直接分類。

    圖9 日冕物質(zhì)拋射圖像分類結(jié)果,圖像來(lái)源于GONZ?LEZ-GóMEZ等(2010)[14]Fig.9 Classification results of CMEs images,image from GONZ?LEZ-GóMEZ et al.(2010)[14]

    日冕物質(zhì)拋射圖像有多尺度特性,多尺度圖像處理技術(shù)對(duì)增強(qiáng)日冕物質(zhì)拋射前沿的可見性和抑制噪聲有很大作用。由于線段更能反應(yīng)圖像信息,而小波更適合識(shí)別點(diǎn)狀特征,如噪聲或背景恒星,不適合檢測(cè)日冕物質(zhì)拋射特征的線結(jié)構(gòu),為此文[15]研究了高階多尺度技術(shù),如脊波和曲線波。與小波變換不同,脊波變換首先進(jìn)行Randon變換,用線參數(shù)取代點(diǎn)參數(shù),之后進(jìn)行小波變換。該變換能夠有效地表示直線奇異性特征,但是不能很好地表示圖像邊緣曲線。而曲線波變換彌補(bǔ)了脊波變換的不足,該方法能夠最優(yōu)地表示圖像的曲線奇異性特征。小波、脊波和曲線波類似,都采用基函數(shù)與信號(hào)的內(nèi)積實(shí)現(xiàn)信號(hào)的稀疏表示,不同的是脊波和曲線波變換的去噪效果好,同時(shí)能夠較好地表達(dá)圖像邊緣信息。分別利用小波變換和曲線波變換過(guò)濾原始圖像,對(duì)比實(shí)驗(yàn)結(jié)果表明,與小波相比,曲線波變換能夠獲得較好的去噪效果,有效地增強(qiáng)了日冕物質(zhì)拋射圖像的前沿結(jié)構(gòu)。

    圖像紋理可以表示為像素灰度值的空間變化的函數(shù),紋理特征的提取除了利用信號(hào)處理領(lǐng)域的線性多尺度變換方法增強(qiáng)圖像的邊緣信息來(lái)檢測(cè)日冕物質(zhì)拋射前沿結(jié)構(gòu)之外,另一種有效的方法是利用灰度共生矩陣[27]捕獲并表征不同區(qū)域的紋理信息。文[16-17]提出了一種基于灰度共生矩陣的非參數(shù)監(jiān)督的日冕物質(zhì)拋射分割方法[16-17],稱之為CORSET算法(CORonal Segmentation Technique),該方法是在有監(jiān)督的區(qū)域競(jìng)爭(zhēng)模型分割方法上的改進(jìn)。區(qū)域競(jìng)爭(zhēng)的目標(biāo)是將圖像分割成多個(gè)區(qū)域,使得每個(gè)區(qū)域中的點(diǎn)具有相似的圖像特征。而將圖像每個(gè)像素的灰度值用作區(qū)域競(jìng)爭(zhēng)方法中的特征向量時(shí),有監(jiān)督的區(qū)域競(jìng)爭(zhēng)模型分割日冕物質(zhì)拋射存在兩個(gè)問(wèn)題:(1)日冕物質(zhì)拋射區(qū)域和背景的灰度直方圖有重疊;(2)日冕物質(zhì)拋射事件和背景的灰度直方圖不遵循正態(tài)分布,即統(tǒng)計(jì)模型未知。這兩個(gè)問(wèn)題將會(huì)導(dǎo)致分類錯(cuò)誤。針對(duì)第1個(gè)問(wèn)題,為了解決直方圖重疊,新算法利用灰度共生矩陣描述日冕物質(zhì)拋射的紋理信息;針對(duì)第2個(gè)問(wèn)題,使用非參數(shù)模型中的卡方統(tǒng)計(jì)檢驗(yàn),它可以評(píng)估觀察到的事件是否遵循某個(gè)特定分布。通過(guò)修改區(qū)域競(jìng)爭(zhēng)運(yùn)動(dòng)方程,將卡方檢驗(yàn)和紋理信息引入方程,同時(shí)在圖像分割中,使用快速水平集算法實(shí)現(xiàn)分割曲線的演變。該分割算法首先圍繞C2遮擋盤形成環(huán)狀區(qū)域,如圖10(a),如果存在日冕物質(zhì)拋射,曲線演變?yōu)槿彰嵛镔|(zhì)拋射輪廓,如圖10(b)~(d);如果不存在日冕物質(zhì)拋射,輪廓就會(huì)消失,如圖11。為了跟蹤日冕物質(zhì)拋射,使用當(dāng)前圖像的分割結(jié)果作為下一圖像的初始輪廓,并根據(jù)修改的區(qū)域競(jìng)爭(zhēng)運(yùn)動(dòng)方程演變輪廓,如果輪廓消失意味著某個(gè)日冕物質(zhì)拋射事件結(jié)束。該算法使用卡方檢驗(yàn)和灰度共生矩陣,以便正確捕獲在差分圖像中觀察到的日冕物質(zhì)拋射紋理信息,可以檢測(cè)和跟蹤具有不同形狀和強(qiáng)度的日冕物質(zhì)拋射,同時(shí)得到了一個(gè)類似CACTus的日冕物質(zhì)拋射概覽圖。

    圖10 日冕物質(zhì)拋射出現(xiàn)時(shí)的輪廓演變,圖像來(lái)源于GOUSSIES等(2010)[17]Fig.10 Contour evolution with a CME,image from GOUSSIES et al.(2010)[17]

    圖11 沒有日冕物質(zhì)拋射時(shí)的輪廓演變,圖像來(lái)源于GOUSSIES等(2010)[17]Fig.11 Contour evolution without CMEs,image from GOUSSIES et al.(2010)[17]

    雖然CORSET算法提供了日冕物質(zhì)拋射事件的邊界信息,但是沒有對(duì)日冕物質(zhì)拋射參數(shù)進(jìn)行明確的定量評(píng)估,為此,在文[17]的基礎(chǔ)上,文[18]通過(guò)添加幾個(gè)新的功能,即自動(dòng)計(jì)算不同的形態(tài)學(xué)和運(yùn)動(dòng)學(xué)參數(shù)擴(kuò)展CORSET,增強(qiáng)了算法的功能,并將獲得的參數(shù)(中心位置角、角寬度和速度)與現(xiàn)有的手工目錄和自動(dòng)檢測(cè)目錄進(jìn)行比較。

    由于日冕物質(zhì)拋射是一種不斷運(yùn)動(dòng)的狀態(tài),因此可以通過(guò)將日冕物質(zhì)拋射結(jié)構(gòu)與背景區(qū)分開來(lái)進(jìn)行檢測(cè),目前大多數(shù)研究采用基于差分的方法檢測(cè)圖像的運(yùn)動(dòng)區(qū)域,但是這種數(shù)值差分會(huì)將噪聲增強(qiáng)到與目標(biāo)相當(dāng)?shù)乃剑m然利用中值濾波器會(huì)抑制噪聲,但也會(huì)平滑小型日冕物質(zhì)拋射特征,同時(shí)差分方法還會(huì)引入時(shí)空串?dāng)_的問(wèn)題。為此,文[19]提出了一種反卷積方法,該算法利用歸一化徑向漸變?yōu)V波器[28]將日冕物質(zhì)拋射圖像分離成靜態(tài)的背景結(jié)構(gòu)和動(dòng)態(tài)的日冕物質(zhì)拋射運(yùn)動(dòng)結(jié)構(gòu),實(shí)驗(yàn)結(jié)果表明,該方法能夠很好地檢測(cè)出弱日冕物質(zhì)拋射。

    日冕物質(zhì)拋射區(qū)域與背景以及其他太陽(yáng)結(jié)構(gòu)如冕流等在形態(tài)、紋理等特征上有明顯差異,在視覺上,冕流的頻譜分布相對(duì)均勻,而日冕物質(zhì)拋射區(qū)域的頻譜發(fā)生頻繁突變,為此,文[20]提出了一種基于頻譜突變分析的日冕物質(zhì)拋射識(shí)別方法。該算法與上述多尺度變換方法一樣,也是一種基于頻域分析的方法,算法利用傅里葉變換得到預(yù)處理后的日冕圖像的頻譜圖,分離出非突變和突變信息后,通過(guò)傅里葉反變換得到突變結(jié)構(gòu)在原圖中的對(duì)應(yīng)信息,確定日冕物質(zhì)拋射的初步分離,最后利用局部穩(wěn)定極值區(qū)域檢測(cè)方法確定日冕物質(zhì)拋射區(qū)域的輪廓。由于日冕物質(zhì)拋射區(qū)域輪廓不規(guī)則,因此提出基于區(qū)域協(xié)方差的外接橢圓方法得到角寬度、速度等參數(shù)。由于該方法分析的是圖像的頻譜,因此在識(shí)別多重和晦暗日冕物質(zhì)拋射時(shí)具有良好的效果。

    上述日冕物質(zhì)拋射檢測(cè)方法是建立在時(shí)域、頻域以及多尺度分析的基礎(chǔ)上,通過(guò)信號(hào)處理中的小波、脊波和曲波等線性多尺度變換方法,研究圖像的頻譜信息以及利用統(tǒng)計(jì)中的灰度共生矩陣描述紋理特征,該類方法在檢測(cè)弱日冕物質(zhì)拋射以及暗日冕物質(zhì)拋射時(shí)表現(xiàn)良好。

    2.3.3 基于光流法

    光流是測(cè)量數(shù)字圖像運(yùn)動(dòng)的一種強(qiáng)大的圖像處理工具,它包含圖像中目標(biāo)運(yùn)動(dòng)的信息。光流算法可以從連續(xù)圖像序列中估計(jì)每個(gè)像素的速度矢量,從而形成該圖像的運(yùn)動(dòng)場(chǎng),可以用來(lái)確定目標(biāo)的運(yùn)動(dòng)情況。

    文[21]提出了一種基于光流法的日冕物質(zhì)拋射檢測(cè)和跟蹤算法,對(duì)光流約束方程添加平滑正則化項(xiàng)——馬爾科夫隨機(jī)場(chǎng)。然后將光流的估計(jì)表示為全局優(yōu)化問(wèn)題,為了找到全局最小值,采用多重網(wǎng)格松弛法[29]進(jìn)行求解。將光流估計(jì)方法應(yīng)用于LASCO C2圖像,觀察到1999年到2004年間大約330個(gè)日冕物質(zhì)拋射并進(jìn)行編目,最終對(duì)10個(gè)具有明顯結(jié)構(gòu)的日冕物質(zhì)拋射進(jìn)行分析,推導(dǎo)出它們的速度場(chǎng),結(jié)果如圖12,上面一行表示輸入圖像,下面一行表示光流結(jié)果,陰影顏色表示速度大小,箭頭方向表示在8×8網(wǎng)格上采樣的速度方向。得到的速度測(cè)量結(jié)果能夠可視化日冕物質(zhì)拋射等離子體的演化。該算法速度非???,可以應(yīng)用于所有日冕物質(zhì)拋射圖像,同時(shí)基于光流算法檢測(cè)到的圖像可以很容易地看到日冕物質(zhì)拋射的擴(kuò)展。

    圖12 光流法檢測(cè)日冕物質(zhì)拋射的結(jié)果,圖像來(lái)源于COLANINNO等(2006)[21]Fig.12 CMEs detection results from an optical flow method,image from COLANINNO et al.(2006)[21]

    與文[21]應(yīng)用全局優(yōu)化的光流算法不同,文[22]提出了一種利用梯度估計(jì)的局部參數(shù)化光流方法檢測(cè)日冕物質(zhì)拋射,通過(guò)Lucas-Kanade算法計(jì)算相鄰兩幀圖像在某一時(shí)段每個(gè)像素點(diǎn)位置的移動(dòng),即速度矢量。而該光流算法有一定的約束條件,僅適應(yīng)于小的運(yùn)動(dòng)目標(biāo),當(dāng)目標(biāo)運(yùn)動(dòng)速度較大時(shí),即幀間運(yùn)動(dòng)較大,算法誤差比較大,因此在運(yùn)動(dòng)分析中引入基于金字塔的多分辨率約束,通過(guò)二次采樣縮小圖像的尺寸,同時(shí)在二次采樣前對(duì)圖像應(yīng)用低通濾波器以減小強(qiáng)度變化對(duì)運(yùn)動(dòng)估計(jì)的影響,提高算法的魯棒性。該方法通過(guò)提出一種可靠運(yùn)動(dòng)區(qū)域的提取策略,得到一個(gè)密集的速度場(chǎng),并在一系列連續(xù)日冕圖像上進(jìn)行測(cè)試來(lái)檢測(cè)和跟蹤相鄰兩幀圖像之間的運(yùn)動(dòng)信息。

    利用光流法檢測(cè)日冕物質(zhì)拋射運(yùn)動(dòng)目標(biāo)的基本思想在于給圖像中的每個(gè)像素點(diǎn)賦予一個(gè)速度矢量,形成一個(gè)矢量場(chǎng),根據(jù)各個(gè)像素點(diǎn)的速度矢量特征對(duì)圖像進(jìn)行動(dòng)態(tài)分析。如果圖像中沒有運(yùn)動(dòng)目標(biāo),則光流矢量在整個(gè)圖像區(qū)域是連續(xù)變化的;當(dāng)圖像中有運(yùn)動(dòng)物體時(shí),目標(biāo)和背景存在相對(duì)運(yùn)動(dòng),兩者形成的速度矢量必然不一樣,如此便可以計(jì)算出圖像中日冕物質(zhì)拋射的位置。

    2.3.4 基于學(xué)習(xí)的方法

    傳統(tǒng)的檢測(cè)方法主要通過(guò)人為定義特征或利用設(shè)定簡(jiǎn)單的閾值進(jìn)行處理,并不能很準(zhǔn)確地檢測(cè)日冕物質(zhì)拋射現(xiàn)象,而且日冕物質(zhì)拋射通常具有明亮的結(jié)構(gòu)和豐富的紋理特征,僅使用單一特征識(shí)別日冕物質(zhì)拋射區(qū)域也不夠準(zhǔn)確。隨著近年來(lái)機(jī)器學(xué)習(xí)與深度學(xué)習(xí)技術(shù)在計(jì)算機(jī)視覺與圖像處理中的應(yīng)用越來(lái)越廣泛[30-31],機(jī)器學(xué)習(xí)中有很多效果良好的分類器被用來(lái)解決各種分類和檢測(cè)問(wèn)題。同時(shí)日冕物質(zhì)拋射檢測(cè)是一個(gè)二分類問(wèn)題,即識(shí)別圖像中是否存在日冕物質(zhì)拋射,因此基于學(xué)習(xí)的方法可應(yīng)用于日冕物質(zhì)拋射檢測(cè)。

    由于日冕物質(zhì)拋射在形狀和空間尺度上變化很大,因此僅使用單一分類器是不夠的。在機(jī)器學(xué)習(xí)領(lǐng)域,文[32]提出的AdaBoost算法能夠集成同一個(gè)訓(xùn)練集訓(xùn)練的不同弱分類器,構(gòu)成一個(gè)強(qiáng)分類器,該強(qiáng)分類器的性能超過(guò)集成的任一弱分類器?;贏daBoost,文[23]提出了一種日冕物質(zhì)拋射圖像分類算法,與上述大多數(shù)方法在極坐標(biāo)變換下的圖像中檢測(cè)日冕物質(zhì)拋射不同,該方法直接將扇形區(qū)域的差分圖像分割成塊,并分析圖像中最亮塊的灰度。通過(guò)實(shí)驗(yàn)表明:切片大小、灰度閾值和亮點(diǎn)分?jǐn)?shù)閾值都會(huì)影響檢測(cè)結(jié)果,因此該方法通過(guò)獲得學(xué)習(xí)樣本并為其設(shè)計(jì)分類模型以得到這3種影響因子的最佳參數(shù)組來(lái)設(shè)計(jì)弱分類器。最后使用AdaBoost組合這些弱分類器,得到最終的強(qiáng)分類器。

    與文[23]類似,文[25]也認(rèn)為日冕物質(zhì)拋射對(duì)應(yīng)圖像中的最亮塊,因此將日冕物質(zhì)拋射的檢測(cè)建模為差分圖像中最亮塊的分類問(wèn)題,同樣提出了一種基于AdaBoost分類日冕物質(zhì)拋射圖像的算法。不同的是該算法將差分圖像進(jìn)行了極坐標(biāo)變換,提取變換后圖像最亮塊的灰度、紋理和方向梯度直方圖(Histogram of Oriented Gradient,HOG)特征,采用多特征融合的方式構(gòu)造分類器。由于日冕物質(zhì)拋射具有不同的外觀、灰度和紋理特征,因此選擇合適特征和分類器非常重要,而決策樹可以自動(dòng)實(shí)現(xiàn)特征選擇和分類,因而將其設(shè)計(jì)為弱分類器,同時(shí)AdaBoost作為集成分類器可以提高單個(gè)分類器的分類能力。由于集成學(xué)習(xí)的群體智慧和決策樹的特征選擇能力,所提出的基于多特征檢測(cè)算法的集成可以獲得更好的檢測(cè)效果。由于該方法主要考慮最亮塊區(qū)域,而忽略了暗腔區(qū)域及其他類似日冕物質(zhì)拋射的太陽(yáng)結(jié)構(gòu)體,容易漏掉比較微弱的日冕物質(zhì)拋射事件。

    針對(duì)分類器收斂速度慢的問(wèn)題,文[24]還提出了一種基于極限學(xué)習(xí)機(jī)[33]的新檢測(cè)算法,該方法通過(guò)幀間差分、極坐標(biāo)變換、濾波去噪完成預(yù)處理。與上述方法類似,同樣基于塊分割圖像的方法,不同的是之前提出的方法是將最亮的部分分割成塊,而該算法是利用固定大小的塊遍歷預(yù)處理后的整個(gè)圖像。塊尺寸的選擇如下:若塊中包含的所有像素的平均灰度值高于指定閾值,該塊被視為亮塊,通過(guò)計(jì)算每個(gè)塊大小下的分類準(zhǔn)確率(正確檢測(cè)樣本與所有樣本之比)選擇最高準(zhǔn)確率的塊尺寸分割圖像。同時(shí)使用灰度和紋理特征組成特征向量并構(gòu)建基于極限學(xué)習(xí)機(jī)的分類器。該方法的優(yōu)點(diǎn)是可以選擇合適的特征,收斂速度快,同時(shí)通過(guò)使用時(shí)空連續(xù)性的決策規(guī)則,可以去除大多數(shù)具有類似日冕物質(zhì)拋射的太陽(yáng)結(jié)構(gòu)體,提高了分類效果。

    上述幾種基于機(jī)器學(xué)習(xí)的方法都將圖像中的最亮塊視為日冕物質(zhì)拋射,因此都是對(duì)最亮塊進(jìn)行檢測(cè)和分類,同時(shí)選擇的分類器也各有特點(diǎn)。但是這些方法都是通過(guò)人為定義特征的提取,由于日冕物質(zhì)拋射具有多種特征,人工選擇的特征不一定能得到好的檢測(cè)效果。而深度學(xué)習(xí)具有強(qiáng)大的學(xué)習(xí)能力和特征表達(dá)能力,會(huì)自動(dòng)提取特征進(jìn)行分類。文[3]提出了一種基于卷積神經(jīng)網(wǎng)絡(luò)(Convolutional Neural Networks,CNN)的日冕物質(zhì)拋射檢測(cè)方法,卷積神經(jīng)網(wǎng)絡(luò)采用有監(jiān)督的學(xué)習(xí)方式,使用海量具有標(biāo)簽的訓(xùn)練樣本,結(jié)合反向傳播算法更新權(quán)重,訓(xùn)練出一個(gè)具有分類能力的網(wǎng)絡(luò)。該方法自動(dòng)提取適應(yīng)于日冕物質(zhì)拋射檢測(cè)的圖像特征并建立檢測(cè)模型,網(wǎng)絡(luò)結(jié)構(gòu)采用包含卷積層、池化層和全連接層的六層卷積神經(jīng)網(wǎng)絡(luò)。為了網(wǎng)絡(luò)模型能夠在訓(xùn)練階段學(xué)習(xí)到較好的參數(shù)進(jìn)行分類任務(wù),訓(xùn)練采用明顯的日冕物質(zhì)拋射圖像數(shù)據(jù)集。同時(shí)以此參數(shù)為基準(zhǔn)在微調(diào)階段采用微弱的日冕物質(zhì)拋射數(shù)據(jù)集進(jìn)行參數(shù)調(diào)優(yōu),使得網(wǎng)絡(luò)模型能獲得較好的分類效果。相較于以往的檢測(cè)方法,基于卷積神經(jīng)網(wǎng)絡(luò)的分類方法在圖形處理器模式下圖像處理效率遠(yuǎn)勝于其他算法,實(shí)時(shí)性好。

    目前基于學(xué)習(xí)方法檢測(cè)日冕物質(zhì)拋射的研究還不多,文獻(xiàn)也很少,與上述傳統(tǒng)方法相比,該類方法分類效果好,檢測(cè)準(zhǔn)確率高,因此,可以將此方法作為未來(lái)日冕物質(zhì)拋射檢測(cè)的一個(gè)重要研究方向。

    3 總結(jié)與展望

    日冕物質(zhì)拋射是一種劇烈的太陽(yáng)爆發(fā)活動(dòng),噴發(fā)出的物質(zhì)攜帶巨大的能量,從而引起地球空間環(huán)境的強(qiáng)烈擾動(dòng),可能對(duì)通訊、導(dǎo)航、航天器等高技術(shù)系統(tǒng)造成災(zāi)害性影響,因此,日冕物質(zhì)拋射的檢測(cè)對(duì)預(yù)防這些災(zāi)害性空間天氣具有非常重要的作用。

    目前雖然已經(jīng)有很多方法應(yīng)用到日冕物質(zhì)拋射的檢測(cè)領(lǐng)域,使得日冕物質(zhì)拋射的檢測(cè)研究取得了很大進(jìn)展,但這仍然是一個(gè)具有挑戰(zhàn)性的問(wèn)題。(1)日冕物質(zhì)拋射具有不同的尺度、大小和形狀,這使得難以選擇適當(dāng)?shù)奶卣鱽?lái)檢測(cè)所有日冕物質(zhì)拋射;(2)日冕物質(zhì)拋射在噴射過(guò)程中伴隨一些與其結(jié)構(gòu)非常相似的干擾源,因此很難將日冕物質(zhì)拋射與這些結(jié)構(gòu)區(qū)分;(3)對(duì)于一些小的、弱的日冕物質(zhì)拋射,通常會(huì)發(fā)生誤檢與漏檢。

    隨著近年來(lái)機(jī)器學(xué)習(xí)技術(shù)在圖像處理與機(jī)器視覺領(lǐng)域獲得巨大成功,可以利用機(jī)器學(xué)習(xí)技術(shù)替代簡(jiǎn)單的閾值分割技術(shù)提高日冕物質(zhì)拋射檢測(cè)效果。同時(shí)深度神經(jīng)網(wǎng)絡(luò)可以自動(dòng)提取合適的特征,在語(yǔ)音識(shí)別、圖像處理等方面都取得了良好的效果。另外,現(xiàn)階段大數(shù)據(jù)和高性能計(jì)算設(shè)備為深度學(xué)習(xí)在日冕物質(zhì)拋射檢測(cè)中的應(yīng)用提供了條件。后續(xù)日冕物質(zhì)拋射檢測(cè)的主要研究方向可以從以下幾方面開展:

    (1)繼續(xù)利用機(jī)器學(xué)習(xí)的方法選擇合適的特征以及分類器,進(jìn)一步提高日冕物質(zhì)拋射的分類效果;

    (2)由于日冕物質(zhì)拋射的爆發(fā)是一個(gè)動(dòng)態(tài)過(guò)程,可以將基于深度學(xué)習(xí)的運(yùn)動(dòng)目標(biāo)檢測(cè)方法應(yīng)用于該檢測(cè);

    (3)可以認(rèn)為日冕物質(zhì)拋射的檢測(cè)是分離日冕圖像中前景運(yùn)動(dòng)目標(biāo)的過(guò)程,因此可以將基于深度學(xué)習(xí)的前景檢測(cè)與背景減除方法應(yīng)用于該檢測(cè)。

    利用基于深度學(xué)習(xí)方法的一個(gè)難點(diǎn)和重點(diǎn)在于數(shù)據(jù)集的制作,手工標(biāo)記需要大量的人力與時(shí)間,同時(shí)對(duì)一些弱、小的日冕物質(zhì)拋射標(biāo)記時(shí)可能會(huì)有誤差和歧義,還需要不斷地深入研究。

    總而言之,隨著人工智能技術(shù)的發(fā)展和進(jìn)步,面向大數(shù)據(jù)分析和挖掘的各類基于學(xué)習(xí)的技術(shù)在太陽(yáng)物理學(xué)研究和應(yīng)用中,將有更為廣泛的應(yīng)用空間和前景。因此,基于機(jī)器學(xué)習(xí)和深度學(xué)習(xí)等技術(shù)可以很好地應(yīng)用到日冕物質(zhì)拋射的檢測(cè)上,有望得到更快速、更準(zhǔn)確的檢測(cè)結(jié)果。

    猜你喜歡
    日冕分類器灰度
    采用改進(jìn)導(dǎo)重法的拓?fù)浣Y(jié)構(gòu)灰度單元過(guò)濾技術(shù)
    基于Faster R-CNN的日冕物質(zhì)拋射檢測(cè)方法*
    Task 3
    基于灰度拉伸的圖像水位識(shí)別方法研究
    BP-GA光照分類器在車道線識(shí)別中的應(yīng)用
    基于最大加權(quán)投影求解的彩色圖像灰度化對(duì)比度保留算法
    加權(quán)空-譜與最近鄰分類器相結(jié)合的高光譜圖像分類
    結(jié)合模糊(C+P)均值聚類和SP-V-支持向量機(jī)的TSK分類器
    基于灰度線性建模的亞像素圖像抖動(dòng)量計(jì)算
    太陽(yáng)上也下“傾盆大雨”
    知識(shí)窗(2014年9期)2014-12-08 08:29:54
    av福利片在线| 97人妻精品一区二区三区麻豆| 搡老熟女国产l中国老女人| 99久久综合精品五月天人人| 国产成人一区二区三区免费视频网站| 亚洲黑人精品在线| 18禁黄网站禁片免费观看直播| 日韩免费av在线播放| 精品久久久久久久久久免费视频| 国产亚洲欧美98| 国产av又大| 黄片小视频在线播放| 国产蜜桃级精品一区二区三区| 午夜a级毛片| 少妇的丰满在线观看| av片东京热男人的天堂| 叶爱在线成人免费视频播放| av有码第一页| 亚洲欧美日韩无卡精品| 国产亚洲精品综合一区在线观看 | 欧美精品亚洲一区二区| 伦理电影免费视频| 国内久久婷婷六月综合欲色啪| 亚洲国产中文字幕在线视频| 亚洲人与动物交配视频| 少妇裸体淫交视频免费看高清 | 午夜福利高清视频| 亚洲一区二区三区色噜噜| 熟妇人妻久久中文字幕3abv| 99久久综合精品五月天人人| 免费看日本二区| 欧美黄色片欧美黄色片| 亚洲欧美日韩东京热| 午夜久久久久精精品| 日韩大码丰满熟妇| 亚洲精品美女久久av网站| 久久国产乱子伦精品免费另类| 男人舔奶头视频| 日韩欧美在线乱码| 精品国内亚洲2022精品成人| 老司机福利观看| 亚洲av成人av| 亚洲国产欧洲综合997久久,| 午夜福利高清视频| 少妇的丰满在线观看| 免费看日本二区| 国产v大片淫在线免费观看| 少妇人妻一区二区三区视频| 日本免费一区二区三区高清不卡| 精品无人区乱码1区二区| 最近在线观看免费完整版| 宅男免费午夜| 亚洲成人中文字幕在线播放| 天天躁狠狠躁夜夜躁狠狠躁| 一本精品99久久精品77| 桃色一区二区三区在线观看| 久久性视频一级片| 一本久久中文字幕| 俄罗斯特黄特色一大片| 日韩精品青青久久久久久| 波多野结衣高清无吗| 国产精品98久久久久久宅男小说| 亚洲av中文字字幕乱码综合| 亚洲 欧美一区二区三区| 桃色一区二区三区在线观看| 啪啪无遮挡十八禁网站| 老司机午夜福利在线观看视频| 欧美zozozo另类| 男女之事视频高清在线观看| 久久亚洲精品不卡| 久久天堂一区二区三区四区| 最好的美女福利视频网| 亚洲人成网站在线播放欧美日韩| 韩国av一区二区三区四区| 国产亚洲精品综合一区在线观看 | 我的老师免费观看完整版| 看片在线看免费视频| 最新美女视频免费是黄的| 欧美成人一区二区免费高清观看 | 欧美成人一区二区免费高清观看 | 日本 欧美在线| 男人的好看免费观看在线视频 | 精品第一国产精品| 欧美最黄视频在线播放免费| 国产人伦9x9x在线观看| xxx96com| 欧美3d第一页| 亚洲五月天丁香| av有码第一页| 国产一区在线观看成人免费| 99国产极品粉嫩在线观看| 99久久久亚洲精品蜜臀av| 老司机靠b影院| 欧美一区二区精品小视频在线| 国产精品久久久久久久电影 | 91成年电影在线观看| 国产精品一区二区三区四区久久| 老汉色∧v一级毛片| 啦啦啦观看免费观看视频高清| 可以在线观看毛片的网站| 国产探花在线观看一区二区| 少妇的丰满在线观看| 丝袜人妻中文字幕| or卡值多少钱| 五月玫瑰六月丁香| 日韩欧美一区二区三区在线观看| 亚洲精品一区av在线观看| 男人舔女人的私密视频| 男人舔女人的私密视频| 好男人电影高清在线观看| 999精品在线视频| 宅男免费午夜| 女警被强在线播放| 欧美性长视频在线观看| 1024香蕉在线观看| 国产成人一区二区三区免费视频网站| www日本在线高清视频| 搡老岳熟女国产| 亚洲精品中文字幕一二三四区| 成人亚洲精品av一区二区| 老司机在亚洲福利影院| 亚洲av美国av| 久久香蕉激情| 成人18禁高潮啪啪吃奶动态图| 无限看片的www在线观看| 动漫黄色视频在线观看| 午夜福利在线在线| 少妇裸体淫交视频免费看高清 | 日韩有码中文字幕| 国产在线精品亚洲第一网站| 97碰自拍视频| 淫秽高清视频在线观看| 免费在线观看成人毛片| 好看av亚洲va欧美ⅴa在| 免费看a级黄色片| 免费高清视频大片| 99在线视频只有这里精品首页| 伦理电影免费视频| 精品不卡国产一区二区三区| 一区福利在线观看| 少妇裸体淫交视频免费看高清 | 亚洲一卡2卡3卡4卡5卡精品中文| 久久久国产成人精品二区| 久久久久久国产a免费观看| 高清在线国产一区| av视频在线观看入口| 欧美三级亚洲精品| 97人妻精品一区二区三区麻豆| 九色成人免费人妻av| 国产一区在线观看成人免费| 亚洲一区中文字幕在线| 国产免费男女视频| 国产在线观看jvid| 九九热线精品视视频播放| 给我免费播放毛片高清在线观看| 人成视频在线观看免费观看| 丁香欧美五月| 国产成+人综合+亚洲专区| 亚洲 欧美一区二区三区| 国产亚洲av高清不卡| 国产精品av视频在线免费观看| 免费高清视频大片| 青草久久国产| 成人永久免费在线观看视频| 欧美+亚洲+日韩+国产| 亚洲精品美女久久久久99蜜臀| 高清毛片免费观看视频网站| 正在播放国产对白刺激| 中文字幕高清在线视频| 最好的美女福利视频网| or卡值多少钱| 精品久久蜜臀av无| 又黄又爽又免费观看的视频| 色噜噜av男人的天堂激情| 国产精品自产拍在线观看55亚洲| 久久久久久亚洲精品国产蜜桃av| 在线十欧美十亚洲十日本专区| 日本 av在线| 99国产精品一区二区蜜桃av| 久久精品国产99精品国产亚洲性色| 黄色毛片三级朝国网站| 人人妻,人人澡人人爽秒播| 丰满人妻一区二区三区视频av | 特大巨黑吊av在线直播| 亚洲中文字幕一区二区三区有码在线看 | 午夜福利成人在线免费观看| 激情在线观看视频在线高清| 国产精品一及| 桃红色精品国产亚洲av| 午夜精品一区二区三区免费看| 欧美激情久久久久久爽电影| 91在线观看av| 亚洲成a人片在线一区二区| 两性夫妻黄色片| av在线天堂中文字幕| or卡值多少钱| 亚洲欧美一区二区三区黑人| 国产黄色小视频在线观看| 亚洲九九香蕉| 国产高清有码在线观看视频 | 婷婷六月久久综合丁香| 久久婷婷人人爽人人干人人爱| 国产三级中文精品| 亚洲午夜精品一区,二区,三区| 麻豆一二三区av精品| 级片在线观看| 两人在一起打扑克的视频| 国产av又大| 国产精品美女特级片免费视频播放器 | 午夜免费观看网址| 91国产中文字幕| 脱女人内裤的视频| 午夜免费成人在线视频| 18美女黄网站色大片免费观看| 热99re8久久精品国产| 成人欧美大片| 两个人免费观看高清视频| 国模一区二区三区四区视频 | 亚洲精品中文字幕在线视频| 一个人免费在线观看电影 | 久久这里只有精品中国| 亚洲成人精品中文字幕电影| 国产精品久久久久久亚洲av鲁大| 成人三级黄色视频| av中文乱码字幕在线| 国产免费av片在线观看野外av| 日韩精品免费视频一区二区三区| 亚洲熟妇中文字幕五十中出| 免费观看人在逋| 精品免费久久久久久久清纯| 午夜福利在线观看吧| tocl精华| 真人一进一出gif抽搐免费| 亚洲成人久久爱视频| 欧美色欧美亚洲另类二区| 亚洲av电影不卡..在线观看| 亚洲人成网站高清观看| 啦啦啦韩国在线观看视频| 9191精品国产免费久久| 99久久99久久久精品蜜桃| 婷婷亚洲欧美| 悠悠久久av| 天天躁狠狠躁夜夜躁狠狠躁| 国产欧美日韩一区二区三| 19禁男女啪啪无遮挡网站| 精品免费久久久久久久清纯| 久久久久国内视频| 可以在线观看的亚洲视频| 成年免费大片在线观看| 在线播放国产精品三级| 中文字幕久久专区| 国内精品一区二区在线观看| 久久久久免费精品人妻一区二区| 变态另类成人亚洲欧美熟女| 伊人久久大香线蕉亚洲五| 国产在线观看jvid| 十八禁网站免费在线| 九色成人免费人妻av| 性色av乱码一区二区三区2| 日韩大码丰满熟妇| 在线看三级毛片| 亚洲av美国av| 91麻豆精品激情在线观看国产| 搞女人的毛片| 婷婷精品国产亚洲av在线| 岛国视频午夜一区免费看| 欧美成人一区二区免费高清观看 | 亚洲精品国产一区二区精华液| 老熟妇乱子伦视频在线观看| 男女午夜视频在线观看| 一级片免费观看大全| 国产人伦9x9x在线观看| 久久久久亚洲av毛片大全| 一夜夜www| 高清在线国产一区| 神马国产精品三级电影在线观看 | 亚洲午夜精品一区,二区,三区| 两个人免费观看高清视频| a在线观看视频网站| 国产av麻豆久久久久久久| 女人被狂操c到高潮| 国产亚洲精品综合一区在线观看 | 高清毛片免费观看视频网站| 怎么达到女性高潮| 亚洲成人精品中文字幕电影| 人成视频在线观看免费观看| 99国产精品99久久久久| ponron亚洲| 性色av乱码一区二区三区2| 一本精品99久久精品77| 久久久久久人人人人人| 最近最新中文字幕大全免费视频| 国产探花在线观看一区二区| 国产高清激情床上av| 很黄的视频免费| 一二三四在线观看免费中文在| АⅤ资源中文在线天堂| 99热这里只有精品一区 | 国产激情偷乱视频一区二区| 99re在线观看精品视频| 久久这里只有精品19| 婷婷亚洲欧美| 亚洲av第一区精品v没综合| 午夜久久久久精精品| 亚洲av成人不卡在线观看播放网| 欧美zozozo另类| 黄色a级毛片大全视频| 女人高潮潮喷娇喘18禁视频| 男插女下体视频免费在线播放| 99久久99久久久精品蜜桃| 日韩精品免费视频一区二区三区| 久久午夜亚洲精品久久| 亚洲一区中文字幕在线| or卡值多少钱| 精品一区二区三区av网在线观看| 国产91精品成人一区二区三区| 老司机午夜十八禁免费视频| 女同久久另类99精品国产91| 听说在线观看完整版免费高清| 制服人妻中文乱码| 国产精品久久久av美女十八| 制服诱惑二区| 一级黄色大片毛片| a级毛片在线看网站| 午夜日韩欧美国产| 国产精品免费一区二区三区在线| 成人国语在线视频| 身体一侧抽搐| a级毛片在线看网站| 欧美国产日韩亚洲一区| 欧美av亚洲av综合av国产av| 亚洲欧美日韩高清专用| 老司机在亚洲福利影院| 一个人免费在线观看电影 | 国产av一区在线观看免费| 男人舔奶头视频| 国内毛片毛片毛片毛片毛片| 中文字幕久久专区| 波多野结衣高清作品| 亚洲欧美精品综合久久99| 一区二区三区高清视频在线| cao死你这个sao货| 久久 成人 亚洲| 久久婷婷人人爽人人干人人爱| 亚洲欧美激情综合另类| 国产精品久久久久久久电影 | 波多野结衣高清无吗| a在线观看视频网站| 亚洲国产精品合色在线| 狠狠狠狠99中文字幕| 精品国内亚洲2022精品成人| 国产蜜桃级精品一区二区三区| av片东京热男人的天堂| 国产精品亚洲一级av第二区| 香蕉丝袜av| 国内精品久久久久精免费| 色噜噜av男人的天堂激情| 99久久国产精品久久久| 精品国产超薄肉色丝袜足j| 亚洲国产精品成人综合色| 国产亚洲精品综合一区在线观看 | 精品国内亚洲2022精品成人| 亚洲成人国产一区在线观看| 熟女少妇亚洲综合色aaa.| 大型黄色视频在线免费观看| 国产成人影院久久av| 日本熟妇午夜| 久久久精品欧美日韩精品| 波多野结衣巨乳人妻| 床上黄色一级片| 窝窝影院91人妻| 天堂√8在线中文| 淫妇啪啪啪对白视频| 精品乱码久久久久久99久播| 99热这里只有是精品50| 亚洲av成人av| 免费一级毛片在线播放高清视频| 亚洲av日韩精品久久久久久密| 看片在线看免费视频| 变态另类丝袜制服| 精品高清国产在线一区| 欧美性猛交╳xxx乱大交人| 亚洲av片天天在线观看| 99国产综合亚洲精品| 一本综合久久免费| 国产成人精品久久二区二区91| 两性午夜刺激爽爽歪歪视频在线观看 | 九色成人免费人妻av| 成人三级做爰电影| 国产三级在线视频| 久久久国产成人免费| 国产视频一区二区在线看| 黄频高清免费视频| netflix在线观看网站| 亚洲精品国产精品久久久不卡| 欧美久久黑人一区二区| 国产欧美日韩一区二区三| 成人三级做爰电影| 亚洲美女黄片视频| 亚洲人成伊人成综合网2020| 成人国产一区最新在线观看| 成人av在线播放网站| 小说图片视频综合网站| 99国产精品一区二区蜜桃av| 91九色精品人成在线观看| 黄片小视频在线播放| www.精华液| 午夜福利欧美成人| 亚洲成a人片在线一区二区| 国产野战对白在线观看| 久久这里只有精品19| 日韩欧美 国产精品| 一级片免费观看大全| 欧美日韩瑟瑟在线播放| 亚洲一码二码三码区别大吗| 久久国产精品影院| 久久久久久久精品吃奶| 国产午夜精品论理片| 91成年电影在线观看| 99精品在免费线老司机午夜| 岛国视频午夜一区免费看| 国产99久久九九免费精品| 国产精品一及| 法律面前人人平等表现在哪些方面| 午夜影院日韩av| 在线看三级毛片| 午夜视频精品福利| av福利片在线| 全区人妻精品视频| 一级片免费观看大全| 国产精品乱码一区二三区的特点| 欧美日韩亚洲综合一区二区三区_| 国产1区2区3区精品| 亚洲熟妇中文字幕五十中出| 人妻丰满熟妇av一区二区三区| 国产精品av视频在线免费观看| 妹子高潮喷水视频| 亚洲aⅴ乱码一区二区在线播放 | cao死你这个sao货| 一区二区三区激情视频| 狂野欧美白嫩少妇大欣赏| 最好的美女福利视频网| 亚洲免费av在线视频| 亚洲精品在线观看二区| 久久这里只有精品19| 久久久久久九九精品二区国产 | 黑人巨大精品欧美一区二区mp4| 久久久久久久精品吃奶| 午夜激情av网站| 久久久久九九精品影院| 麻豆久久精品国产亚洲av| 99久久国产精品久久久| 中文字幕人成人乱码亚洲影| 午夜久久久久精精品| 亚洲男人天堂网一区| 国内精品久久久久精免费| 19禁男女啪啪无遮挡网站| 黄色毛片三级朝国网站| 欧美日韩乱码在线| www日本在线高清视频| 国产精品久久视频播放| 在线观看www视频免费| 禁无遮挡网站| 大型黄色视频在线免费观看| e午夜精品久久久久久久| 我的老师免费观看完整版| 欧美色欧美亚洲另类二区| 女人高潮潮喷娇喘18禁视频| 老司机午夜福利在线观看视频| 搡老岳熟女国产| 一卡2卡三卡四卡精品乱码亚洲| 在线观看一区二区三区| 校园春色视频在线观看| 美女高潮喷水抽搐中文字幕| 欧美最黄视频在线播放免费| 琪琪午夜伦伦电影理论片6080| 精品久久久久久久人妻蜜臀av| 色综合站精品国产| 成人永久免费在线观看视频| 变态另类丝袜制服| 91在线观看av| 极品教师在线免费播放| 欧美不卡视频在线免费观看 | 久久久久免费精品人妻一区二区| a级毛片a级免费在线| www日本在线高清视频| 母亲3免费完整高清在线观看| 久热爱精品视频在线9| 又大又爽又粗| 国产99白浆流出| 俺也久久电影网| 日韩免费av在线播放| xxxwww97欧美| 丁香欧美五月| 欧美中文综合在线视频| 18禁国产床啪视频网站| 日本一本二区三区精品| 舔av片在线| 欧美绝顶高潮抽搐喷水| 一个人观看的视频www高清免费观看 | av片东京热男人的天堂| 一级a爱片免费观看的视频| 欧美另类亚洲清纯唯美| 日韩欧美一区二区三区在线观看| 久久久久久久久免费视频了| 中文字幕人妻丝袜一区二区| 免费一级毛片在线播放高清视频| 久久久精品欧美日韩精品| 国产精品1区2区在线观看.| 国产精品久久视频播放| 不卡一级毛片| 国产成人av教育| av中文乱码字幕在线| 99国产精品99久久久久| 精品久久久久久久毛片微露脸| 国产成+人综合+亚洲专区| 一夜夜www| 国产黄片美女视频| 亚洲av熟女| 国产一级毛片七仙女欲春2| 精品熟女少妇八av免费久了| 神马国产精品三级电影在线观看 | 日韩欧美精品v在线| 老司机午夜福利在线观看视频| 久久精品夜夜夜夜夜久久蜜豆 | 亚洲欧美精品综合一区二区三区| 久久午夜亚洲精品久久| 欧美+亚洲+日韩+国产| 国产精品亚洲av一区麻豆| 国产在线观看jvid| 美女扒开内裤让男人捅视频| 18禁国产床啪视频网站| 精品国产乱子伦一区二区三区| 国产亚洲精品久久久久5区| 白带黄色成豆腐渣| 99riav亚洲国产免费| 不卡一级毛片| 夜夜夜夜夜久久久久| 香蕉av资源在线| 精品国内亚洲2022精品成人| 欧美性猛交╳xxx乱大交人| 一本大道久久a久久精品| 一个人观看的视频www高清免费观看 | 法律面前人人平等表现在哪些方面| 国产精品,欧美在线| 91九色精品人成在线观看| 99久久无色码亚洲精品果冻| 国产亚洲精品久久久久久毛片| 麻豆国产av国片精品| 麻豆久久精品国产亚洲av| 人妻久久中文字幕网| 色综合亚洲欧美另类图片| 99热这里只有是精品50| 露出奶头的视频| 老司机深夜福利视频在线观看| 免费人成视频x8x8入口观看| 熟女电影av网| 久久久国产成人精品二区| 欧美黄色片欧美黄色片| 又黄又粗又硬又大视频| 久久午夜亚洲精品久久| 一a级毛片在线观看| 精品久久久久久久人妻蜜臀av| 午夜福利在线观看吧| 国产精品影院久久| 日韩欧美国产在线观看| 国产伦人伦偷精品视频| 老汉色∧v一级毛片| 精品午夜福利视频在线观看一区| 变态另类丝袜制服| 亚洲精品一卡2卡三卡4卡5卡| 精品一区二区三区四区五区乱码| 丁香欧美五月| 观看免费一级毛片| 黄频高清免费视频| 日韩中文字幕欧美一区二区| 久99久视频精品免费| 亚洲国产精品sss在线观看| 午夜精品一区二区三区免费看| 免费看美女性在线毛片视频| 精品福利观看| 麻豆成人午夜福利视频| 国产男靠女视频免费网站| 欧美大码av| 淫秽高清视频在线观看| 欧美精品啪啪一区二区三区| 高清毛片免费观看视频网站| www.自偷自拍.com| 特大巨黑吊av在线直播| 精品国产美女av久久久久小说| 欧美高清成人免费视频www| 两个人视频免费观看高清| 日本黄色视频三级网站网址| 老司机午夜十八禁免费视频| 日本精品一区二区三区蜜桃| 亚洲国产日韩欧美精品在线观看 | 国产欧美日韩一区二区三| 免费在线观看日本一区| 俄罗斯特黄特色一大片| 久久香蕉精品热| 女警被强在线播放| 国产免费av片在线观看野外av| 中文字幕久久专区| 午夜亚洲福利在线播放| 亚洲熟女毛片儿| 99精品久久久久人妻精品| 午夜免费成人在线视频| 亚洲国产中文字幕在线视频| 亚洲国产欧洲综合997久久,| 亚洲专区字幕在线| 大型av网站在线播放| 亚洲精品一卡2卡三卡4卡5卡| av福利片在线| 午夜福利欧美成人| 真人做人爱边吃奶动态|