專利名稱：磁共振成像采樣軌跡優(yōu)化方法
磁共振成像采樣軌跡優(yōu)化方法
技術(shù)領(lǐng)域：
本發(fā)明涉及核磁共振成像領(lǐng)域，尤其涉及一種磁共振成像采樣軌跡優(yōu)化方法。背景技
磁共振成像是當(dāng)前臨床醫(yī)學(xué)影像學(xué)的重要檢查手段之一，具有無創(chuàng)，無輻射，高分辨率，靈活的選擇成像參數(shù)和層面等優(yōu)點(diǎn)。然而，磁共振成像的最大缺點(diǎn)在于檢查時(shí)間較長(zhǎng)，縮短檢查時(shí)間不僅能夠節(jié)約成本，同時(shí)是心血管檢查，介入手術(shù)等實(shí)時(shí)動(dòng)態(tài)成像的關(guān)鍵所在，因此，縮短成像時(shí)間一直是磁共振技術(shù)發(fā)展的重要目標(biāo)。在磁共振成像中，成像時(shí)間長(zhǎng)的主要原因在于K空間中相位編碼線的獲取時(shí)間較長(zhǎng)。因此，需要對(duì)采樣軌跡進(jìn)行優(yōu)化，在較少的相位編碼線的情況下得到高分辨率的圖像是解決成像速度問題的關(guān)鍵。近年來出現(xiàn)的壓縮傳感理論很好的解決的了這一問題。傳統(tǒng)壓縮傳感理論對(duì)是基于K空間的隨機(jī)采樣進(jìn)行重建，但是此方法的下采樣率較低，重建效果差，雖然能夠使磁共振成像的時(shí)間變短，但是不能得到所要求質(zhì)量的重建圖像。
發(fā)明內(nèi)容有鑒于此，有必要提供一種重建質(zhì)量好且用時(shí)少的磁共振成像采樣軌跡優(yōu)化方法。一種磁共振成像采樣軌跡優(yōu)化方法，包括如下步驟由全采樣的K空間中選取多條相位編碼線，作為初始的采樣軌跡，剩余相位編碼線作為候選集；對(duì)采樣軌跡進(jìn)行重建，得到重建圖像和噪聲壓制因子；對(duì)重建圖像進(jìn)行反變換，得到新的K空間；將新的K空間與全采樣的K空間進(jìn)行對(duì)比，得出K空間中每條相位編碼線的均方差;對(duì)候選集中低頻段的相位編碼線進(jìn)行采集，采集將低頻段的相位編碼線按照均方差由大到小排列的前幾條，并加入初始的采樣軌跡中，得到低頻段加入后的采樣軌跡；對(duì)加入的采樣軌跡進(jìn)行重建，得到重建圖像及噪聲壓制因子，判斷噪聲壓制因子是否不變，若不變則低頻段采集完成，獲得低頻段采集完后的采樣軌跡，否則返回對(duì)重建圖像進(jìn)行反變換，得到新的K空間的步驟進(jìn)行循環(huán)，直至低頻段采集完成；對(duì)候選集中高頻段的相位編碼線進(jìn)行采集，將高頻段的相位編碼線分為多個(gè)區(qū)，采集每個(gè)區(qū)中均方差最大的相位編碼線加入低頻段采集完后的采樣軌跡中，得到高頻段加入后的采樣軌跡；對(duì)高頻段加入后的采樣軌跡進(jìn)行重建得到重建圖像，判斷重建圖像是否達(dá)到設(shè)定圖像質(zhì)量，若達(dá)到則優(yōu)化結(jié)束，獲得最終的采樣軌跡，否則返回對(duì)重建圖像進(jìn)行反變換，得到新的K空間的步驟。優(yōu)選的，所述由全采樣的K空間中選取多條相位編碼線，作為初始的采樣軌跡，剩余相位編碼線作為候選集的步驟具體為由全采樣的K空間的低頻段采集多條相位編碼線，作為初始的采樣軌跡。優(yōu)選的，所述由全采樣的K空間中選取多條相位編碼線，作為初始的采樣軌跡，剩余相位編碼線作為候選集的步驟具體為相對(duì)全采樣的K空間的中心，對(duì)稱的采集多條相位編碼線。優(yōu)選的，所述對(duì)采樣軌跡進(jìn)行重建，得到重建圖像和噪聲壓制因子的步驟具體為 采用非線性方法對(duì)采樣軌跡進(jìn)行重建，得到重建圖像。優(yōu)選的，所述對(duì)重建圖像進(jìn)行反變換，得到新的K空間的步驟具體為對(duì)重建圖像進(jìn)行傅里葉反變換，得到新的K空間。優(yōu)選的，所述對(duì)候選集中低頻段的相位編碼線進(jìn)行采集，采集將低頻段的相位編碼線按照均方差由大到小排列的前幾條，并加入初始的采樣軌跡中，得到低頻段加入后的采樣軌跡的步驟具體為采集將低頻段的相位編碼線按照均方差由大到小排列的前10至15條，加入到初始的采樣軌跡中。優(yōu)選的，所述對(duì)候選集中高頻段的相位編碼線進(jìn)行采集，將高頻段的相位編碼線分為多個(gè)區(qū)，采集每個(gè)區(qū)中均方差最大的相位編碼線加入低頻段采集完后的采樣軌跡中，得到高頻段加入后的采樣軌跡的步驟具體為將候選集高頻段的相位編碼線分為不大于18個(gè)區(qū)，由每個(gè)區(qū)中采集均方差最大的相位編碼線，加入到采樣軌跡中。一種磁共振成像采樣軌跡優(yōu)化方法，包括如下步驟由全采樣的K空間中選取多條相位編碼線，作為初始的采樣軌跡，剩余相位編碼線作為候選集；對(duì)采樣軌跡進(jìn)行重建，得到重建圖像和噪聲壓制因子；對(duì)重建圖像進(jìn)行反變換，得到新的K空間；將新的K空間與全采樣的K空間進(jìn)行對(duì)比，得出K空間中每條相位編碼線的均方差;對(duì)候選集中低頻段的相位編碼線進(jìn)行采集，采集將低頻段的相位編碼線按照均方差由大到小排列的前幾條，并加入初始的采樣軌跡中，得到低頻段加入后的采樣軌跡；對(duì)加入的采樣軌跡進(jìn)行重建，得到重建圖像及噪聲壓制因子，判斷噪聲壓制因子是否不變，若不變則低頻段采集完成，獲得低頻段采集完后的采樣軌跡，否則返回對(duì)重建圖像進(jìn)行反變換，得到新的K空間的步驟進(jìn)行循環(huán)，直至低頻段采集完成；對(duì)候選集中高頻段的相位編碼線進(jìn)行采集，將高頻段的相位編碼線分為多個(gè)區(qū)，采集每個(gè)區(qū)中均方差最大的相位編碼線加入低頻段采集完后的采樣軌跡中，得到高頻段加入后的采樣軌跡；判斷高頻段加入后的采樣軌跡是否達(dá)到設(shè)定的加速因子，若達(dá)到則優(yōu)化結(jié)束，獲得最終的采樣軌跡，否則對(duì)高頻段加入后的采樣軌跡進(jìn)行重建，得到重建圖像并返回對(duì)重建圖像進(jìn)行反變換，得到新的K空間的步驟。
上述磁共振成像采樣軌跡優(yōu)化方法，初始選取幾條相位編碼線作為初始的采樣軌跡，利用初始的采樣軌跡進(jìn)行重建得到重建圖像和噪聲壓制因子，然后反變換得到新的K空間，與全采樣的K空間進(jìn)行比較獲取均方差，然后根據(jù)均方差分別對(duì)低頻段的相位編碼線和高頻段的相位編碼線進(jìn)行采集，且低頻段和高頻段采用不同的采集策略，對(duì)采樣軌跡進(jìn)行優(yōu)化，優(yōu)化效率高，保證了較高的下采樣因子和重建圖像質(zhì)量，又大大降低了重建用時(shí)，提高了重建的效率。
圖I是一個(gè)實(shí)施例中磁共振成像采樣軌跡優(yōu)化方法的流程圖；圖2是傳統(tǒng)利用隨機(jī)選取采用軌跡重建的重建圖像；圖3是一個(gè)實(shí)施例中利用初始的采樣軌跡進(jìn)行重建的重建圖像； 圖4是一個(gè)實(shí)施例中利用優(yōu)化后采樣軌跡重建的重建圖像。
具體實(shí)施方式
下面結(jié)合附圖，對(duì)本發(fā)明的具體實(shí)施方式
進(jìn)行詳細(xì)描述。圖I是一個(gè)實(shí)施例中壓縮感知磁共振成像采樣軌跡優(yōu)化方法的流程圖。該方法包括如下步驟SlO :由全采樣的K空間中采集多條相位編碼線，作為初始的采樣軌跡，剩余相位編碼線作為候選集。全采樣的K空間中具有較多的相位編碼線，初始基于壓縮感知磁共振成像原理，采集多條相位編碼線作為初始的采樣軌跡，剩余的相位編碼線作為候選集。由于K空間具有中心對(duì)稱的性質(zhì)，且低頻段包含了圖像的全局特征，大部分能量集中在低頻段，高頻段僅包含了圖像的細(xì)節(jié)信息，該實(shí)施例中，由全采樣的K空間中采集多條相位編碼線，作為初始的采樣軌跡的步驟具體為由全采樣的K空間的低頻段采集多條相位編碼線，作為初始的采樣軌跡。在其他實(shí)施方式中，相對(duì)全采樣的K空間的中心，對(duì)稱的采集多條相位編碼線。S20 :對(duì)采樣軌跡進(jìn)行重建，得到重建圖像以及噪聲壓制因子。采集初始的采樣軌跡后，需要進(jìn)行重建，得到重建圖像以及噪聲壓制因子。噪聲壓住參數(shù)用來判斷K空間低頻段相位編碼線是否已經(jīng)采集完成。若低頻段采集完成則噪聲壓制因子不再變化。該實(shí)施例中，采用非線性方法對(duì)采樣軌跡進(jìn)行重建，得到重建圖像。S30 :對(duì)重建圖像進(jìn)行反變換，得到新的K空間。根據(jù)重建圖像的重建過程，在得到重建圖像后，對(duì)重建圖像進(jìn)行反變換可以得到新的K空間。該實(shí)施例中，對(duì)重建圖像進(jìn)行傅里葉反變換，得到新的K空間。S40 :將新的K空間與全采樣的K空間進(jìn)行對(duì)比，得出候選集中每條相位編碼線的均方差。在得到新的K空間后，與全采樣的K空間進(jìn)行對(duì)比，能夠得出K空間的均方差，即得出K空間中每條相位編碼線的均方差，也即得出了候選集中每條相位編碼線對(duì)應(yīng)的均方差。S50:對(duì)候選集中低頻段的相位編碼線進(jìn)行采集，采集將低頻段的相位編碼線按照均方差由大到小排列的前幾條，并加入初始的采樣軌跡中，得到低頻段加入后的采樣軌跡。
由于K空間低頻段包含了圖像的全局特征，在得到每條相位編碼線的均方差后，首先對(duì)候選集中低頻段的相位編碼線進(jìn)行采集，加入到初始的采樣軌跡中，增加低頻段的相位編碼線到初始的采樣軌跡中，更新初始的采樣軌跡，得到低頻段加入后的采樣軌跡。此時(shí)，初始的采樣軌跡已變成低頻段加入后的采樣軌跡，對(duì)采樣軌跡進(jìn)行了優(yōu)化。該實(shí)施例中，對(duì)于候選集中低頻段的相位編碼線，采集將低頻段的相位編碼線按照均方差由大到小排列的前10至15條，加入到初始的采樣軌 跡中。S60:對(duì)低頻段加入后的采樣軌跡進(jìn)行重建，得到重建圖像及噪聲壓制因子，判斷噪聲壓制因子是否不變，若不變則低頻段采集完成，獲得低頻段采集完后的采樣軌跡，進(jìn)行步驟S70，否則返回S30進(jìn)行循環(huán)，直至低頻段采集完成。得到低頻段加入后的采樣軌跡后，進(jìn)行重建獲取重建圖像及噪聲壓制因子，此時(shí)需要知道低頻段相位編碼線是否采集完成，則要判斷噪聲壓制因子是否保持不變。如果低頻段的相位編碼線采集完成，采樣軌跡中再加入低頻段的相位編碼線進(jìn)行重建后，噪聲壓制因子保持不變，不會(huì)再變化，則低頻段采集就完成，不再對(duì)低頻段進(jìn)行采集。否則的話，需要返回步驟S30，繼續(xù)循環(huán)，對(duì)低頻段的相位編碼線進(jìn)行采集，不斷往上一循環(huán)的低頻段加入后的采樣軌跡中增加相位編碼線，直至采樣軌跡重建后，噪聲壓制因子相比之前不再變化為止。S70:對(duì)候選集中高頻段的相位編碼線進(jìn)行采集，將高頻段的相位編碼線分為多個(gè)區(qū)，采集每個(gè)區(qū)中均方差最大的相位編碼線加入低頻段采集完后的采樣軌跡中，獲得高頻段加入后的采樣軌跡。在對(duì)候選集中低頻段的相位編碼線采集完成后，由于低頻段的包含的是圖像的全局特征，而圖像的細(xì)節(jié)信息包含在高頻段，為更好的重建圖像，還需對(duì)高頻段的相位編碼線進(jìn)行采集，加入到采樣軌跡中，進(jìn)一步優(yōu)化采樣軌跡。該實(shí)施例中，將候選集高頻段的相位編碼線分為不大于18個(gè)區(qū)，由每個(gè)區(qū)中采集均方差最大的相位編碼線，加入到低頻段采集完后的采樣軌跡中。S80 :對(duì)高頻段加入后的采樣軌跡進(jìn)行重建得到重建圖像；判斷重建圖像是否達(dá)到設(shè)定圖像質(zhì)量，若達(dá)到則優(yōu)化結(jié)束，進(jìn)行步驟S90 :獲取最終的采樣軌跡，否則返回步驟S30。對(duì)高頻段加入后的采樣軌跡進(jìn)行重建，得到重建圖像判斷重建圖像是否達(dá)到設(shè)定的圖像質(zhì)量，若達(dá)到則說明采樣軌跡滿足要求，否則返回步驟S30繼續(xù)采集高頻段的相位編碼線，加入到采樣軌跡中，直至采樣軌跡重建后滿足設(shè)定的圖像質(zhì)量。在其他實(shí)施例中，該步驟還可為判斷高頻段加入后的采樣軌跡是否達(dá)到設(shè)定的加速因子，若達(dá)到則進(jìn)行步驟S90，否則對(duì)高頻段加入后的采樣軌跡進(jìn)行重建，得到重建圖像并返回步驟S30。此外，加速因子用來決定整個(gè)采樣軌跡需要采集的相位編碼線的條數(shù)。在初始的時(shí)候進(jìn)行設(shè)定。在低頻段采集完成后，高頻段還需按照設(shè)定的加速因子，采集達(dá)到設(shè)定加速因子條數(shù)的相位編碼線，若一次不能采集夠，則需要回到步驟S30，重復(fù)進(jìn)行采集，直至達(dá)到設(shè)定的加速因子為止。由于K空間具有中心對(duì)稱性質(zhì)，該方法可以先在半K空間中對(duì)低頻段以及高頻段的相位編碼線進(jìn)行采集，對(duì)采樣軌跡進(jìn)行優(yōu)化，然后對(duì)稱的獲取另一半K空間的相位編碼線，對(duì)采樣軌跡進(jìn)行優(yōu)化。圖2是傳統(tǒng)利用隨機(jī)選取采用軌跡重建的重建圖像。圖3是一個(gè)實(shí)施例中利用初始的采樣軌跡進(jìn)行重建的重建圖像；圖4是一個(gè)實(shí)施例中利用優(yōu)化后采樣軌跡重建的重建圖像。結(jié)合圖2至圖4，該方法對(duì)256條相位編碼線的采樣軌跡進(jìn)行優(yōu)化，以第129條相位編碼線為中心，對(duì)稱的選取20條相位編碼線，共21條相位編碼線作為初始的采樣軌跡，剩余的239條相位編碼線為候選集。在低頻段，每次采集10條相位編碼先加入初始的采樣軌跡，循環(huán)兩次至噪聲壓制因子不變。在高頻段，將剩余的相位編碼線劃分為18個(gè)子模塊，選擇每個(gè)模塊中均方差最大的相位編碼線加入采樣軌跡，從而使得采樣軌跡中低頻段的相位編碼線增多，整個(gè)過程僅耗時(shí)56秒，下采樣因子為4. 34，下采樣因子提高，重建圖像質(zhì)量高
且大大節(jié)省了時(shí)間。在相同的下采樣因子下，與基于K空間分布的隨機(jī)采樣相比，圖像的均方差有所下降，用時(shí)更少，圖像質(zhì)量也大大提高。該方法，初始選取幾條相位編碼線作為初始的采樣軌跡，利用初始的采樣軌跡進(jìn)行重建得到重建圖像和噪聲壓制因子，然后反變換得到新的K空間，與全采樣的K空間進(jìn)行比較獲取均方差，然后根據(jù)均方差分別對(duì)低頻段的相位編碼線和高頻段的相位編碼線進(jìn)行采集，且低頻段和高頻段采用不同的采集策略，對(duì)采樣軌跡進(jìn)行優(yōu)化，優(yōu)化效果高，保證了較高的下采樣因子和重建圖像質(zhì)量，又大大降低了重建用時(shí)，提高了重建的效率。以上所述實(shí)施例僅表達(dá)了本發(fā)明的幾種實(shí)施方式，其描述較為具體和詳細(xì)，但并不能因此而理解為對(duì)本發(fā)明專利范圍的限制。應(yīng)當(dāng)指出的是，對(duì)于本領(lǐng)域的普通技術(shù)人員來說，在不脫離本發(fā)明構(gòu)思的前提下，還可以做出若干變形和改進(jìn)，這些都屬于本發(fā)明的保護(hù)范圍。因此，本發(fā)明專利的保護(hù)范圍應(yīng)以所附權(quán)利要求為準(zhǔn)。
權(quán)利要求
1.一種磁共振成像采樣軌跡優(yōu)化方法，包括如下步驟 由全采樣的K空間中選取多條相位編碼線，作為初始的采樣軌跡，剩余相位編碼線作為候選集； 對(duì)采樣軌跡進(jìn)行重建，得到重建圖像和噪聲壓制因子； 對(duì)重建圖像進(jìn)行反變換，得到新的K空間； 將新的K空間與全采樣的K空間進(jìn)行對(duì)比，得出K空間中每條相位編碼線的均方差； 對(duì)候選集中低頻段的相位編碼線進(jìn)行采集，采集將低頻段的相位編碼線按照均方差由大到小排列的前幾條，并加入初始的采樣軌跡中，得到低頻段加入后的采樣軌跡； 對(duì)加入的采樣軌跡進(jìn)行重建，得到重建圖像及噪聲壓制因子，判斷噪聲壓制因子是否不變，若不變則低頻段采集完成，獲得低頻段采集完后的采樣軌跡，否則返回對(duì)重建圖像進(jìn)行反變換，得到新的K空間的步驟進(jìn)行循環(huán)，直至低頻段采集完成； 對(duì)候選集中高頻段的相位編碼線進(jìn)行采集，將高頻段的相位編碼線分為多個(gè)區(qū)，采集每個(gè)區(qū)中均方差最大的相位編碼線加入低頻段采集完后的采樣軌跡中，得到高頻段加入后的采樣軌跡； 對(duì)高頻段加入后的采樣軌跡進(jìn)行重建得到重建圖像，判斷重建圖像是否達(dá)到設(shè)定圖像質(zhì)量，若達(dá)到則優(yōu)化結(jié)束，獲得最終的采樣軌跡，否則返回對(duì)重建圖像進(jìn)行反變換，得到新的K空間的步驟。
2.根據(jù)權(quán)利要求I所述的磁共振成像采樣軌跡優(yōu)化方法，其特征在于，所述由全采樣的K空間中選取多條相位編碼線，作為初始的采樣軌跡，剩余相位編碼線作為候選集的步驟具體為 由全采樣的K空間的低頻段采集多條相位編碼線，作為初始的采樣軌跡。
3.根據(jù)權(quán)利要求I所述的磁共振成像采樣軌跡優(yōu)化方法，其特征在于，所述由全采樣的K空間中選取多條相位編碼線，作為初始的采樣軌跡，剩余相位編碼線作為候選集的步驟具體為 相對(duì)全采樣的K空間的中心，對(duì)稱的采集多條相位編碼線。
4.根據(jù)權(quán)利要求I所述的磁共振成像采樣軌跡優(yōu)化方法，其特征在于，所述對(duì)采樣軌跡進(jìn)行重建，得到重建圖像和噪聲壓制因子的步驟具體為 采用非線性方法對(duì)采樣軌跡進(jìn)行重建，得到重建圖像。
5.根據(jù)權(quán)利要求I所述的磁共振成像采樣軌跡優(yōu)化方法，其特征在于，所述對(duì)重建圖像進(jìn)行反變換，得到新的K空間的步驟具體為 對(duì)重建圖像進(jìn)行傅里葉反變換，得到新的K空間。
6.根據(jù)權(quán)利要求I所述的磁共振成像采樣軌跡優(yōu)化方法，其特征在于，所述對(duì)候選集中低頻段的相位編碼線進(jìn)行采集，采集將低頻段的相位編碼線按照均方差由大到小排列的前幾條，并加入初始的采樣軌跡中，得到低頻段加入后的采樣軌跡的步驟具體為 采集將低頻段的相位編碼線按照均方差由大到小排列的前10至15條，加入到初始的采樣軌跡中。
7.根據(jù)權(quán)利要求I所述的磁共振成像采樣軌跡優(yōu)化方法，其特征在于，所述對(duì)候選集中高頻段的相位編碼線進(jìn)行采集，將高頻段的相位編碼線分為多個(gè)區(qū)，采集每個(gè)區(qū)中均方差最大的相位編碼線加入低頻段采集完后的采樣軌跡中，得到高頻段加入后的采樣軌跡的步驟具體為將候選集高頻段的相位編碼線分為不大于18個(gè)區(qū)，由每個(gè)區(qū)中采集均方差最大的相位編碼線，加入到采樣軌跡中。
8.—種磁共振成像采樣軌跡優(yōu)化方法，包括如下步驟 由全采樣的K空間中選取多條相位編碼線，作為初始的采樣軌跡，剩余相位編碼線作為候選集； 對(duì)采樣軌跡進(jìn)行重建，得到重建圖像和噪聲壓制因子； 對(duì)重建圖像進(jìn)行反變換，得到新的K空間； 將新的K空間與全采樣的K空間進(jìn)行對(duì)比，得出K空間中每條相位編碼線的均方差； 對(duì)候選集中低頻段的相位編碼線進(jìn)行采集，采集將低頻段的相位編碼線按照均方差由大到小排列的前幾條，并加入初始的采樣軌跡中，得到低頻段加入后的采樣軌跡； 對(duì)加入的采樣軌跡進(jìn)行重建，得到重建圖像及噪聲壓制因子，判斷噪聲壓制因子是否不變，若不變則低頻段采集完成，獲得低頻段采集完后的采樣軌跡，否則返回對(duì)重建圖像進(jìn)行反變換，得到新的K空間的步驟進(jìn)行循環(huán)，直至低頻段采集完成； 對(duì)候選集中高頻段的相位編碼線進(jìn)行采集，將高頻段的相位編碼線分為多個(gè)區(qū)，采集每個(gè)區(qū)中均方差最大的相位編碼線加入低頻段采集完后的采樣軌跡中，得到高頻段加入后的采樣軌跡； 判斷高頻段加入后的采樣軌跡是否達(dá)到設(shè)定的加速因子，若達(dá)到則優(yōu)化結(jié)束，獲得最終的采樣軌跡，否則對(duì)高頻段加入后的采樣軌跡進(jìn)行重建，得到重建圖像并返回對(duì)重建圖像進(jìn)行反變換，得到新的K空間的步驟。
全文摘要
本發(fā)明涉及一種磁共振成像采樣軌跡優(yōu)化方法，該方法初始選取幾條相位編碼線作為初始的采樣軌跡，利用初始的采樣軌跡進(jìn)行重建得到重建圖像和噪聲壓制因子，然后反變換得到新的K空間，與全采樣的K空間進(jìn)行比較獲取均方差，然后根據(jù)均方差分別對(duì)低頻段的相位編碼線和高頻段的相位編碼線進(jìn)行采集，且低頻段和高頻段采用不同的采集策略，對(duì)采樣軌跡進(jìn)行優(yōu)化，優(yōu)化效率高，保證了較高的下采樣因子和重建圖像質(zhì)量，又大大降低了重建用時(shí)，提高了重建的效率。
文檔編號(hào)G01R33/56GK102914753SQ20121039028
公開日2013年2月6日 申請(qǐng)日期2012年10月15日 優(yōu)先權(quán)日2011年12月8日
發(fā)明者梁棟, 劉端端, 劉新, 鄭海榮 申請(qǐng)人:中國(guó)科學(xué)院深圳先進(jìn)技術(shù)研究院