專(zhuān)利名稱(chēng)：基于螺旋馬鞍線的錐形束ct精確濾波反投影重建方法
技術(shù)領(lǐng)域：
本發(fā)明涉及的是一種生物醫(yī)學(xué)成像技術(shù)領(lǐng)域的方法，具體地說(shuō)，是一種基于 螺旋馬鞍線的錐形束CT精確濾波反投影重建方法。
背景技術(shù)：
錐形束CT的X射線球管在一點(diǎn)發(fā)出錐形束X射線，穿過(guò)被測(cè)對(duì)象后，衰減 的X射線被對(duì)側(cè)的檢測(cè)器檢測(cè)，并將X光轉(zhuǎn)化為電信號(hào)，送至計(jì)算機(jī)的模數(shù)轉(zhuǎn)換 器和數(shù)據(jù)采集器，轉(zhuǎn)化為數(shù)字信號(hào)并存儲(chǔ)在計(jì)算機(jī)中。這只是完成一次的曝光和 采集過(guò)程，為了得到重建被測(cè)對(duì)象所需的投影數(shù)據(jù)，X射線源需沿一定的軌跡掃 描，并在每個(gè)位置完成一次曝光和數(shù)據(jù)采集過(guò)程。目前X射線源常用掃描方式有 兩種， 一種以被測(cè)物體為參考系觀察，X射線源的軌跡為單螺旋，螺旋錐形束CT 能夠較快地采集投影數(shù)據(jù)，實(shí)現(xiàn)長(zhǎng)物體的體積重建；另外一種以被測(cè)物體為參考 系觀察，X射線源的軌跡為馬鞍線，基于馬鞍線的錐形束CT具有周期性的特點(diǎn)， 能夠?qū)崿F(xiàn)動(dòng)態(tài)物體如心臟的快速成像和連續(xù)觀察。
經(jīng)對(duì)現(xiàn)有技術(shù)的文獻(xiàn)檢索發(fā)現(xiàn)，A. Katsevich在Advances in A卯lied Mathematics [應(yīng)用數(shù)學(xué)進(jìn)展](2004年第32期，第681頁(yè)至697頁(yè))上發(fā)表的 "Improved exact FBP algorithm for spiral CT"[改進(jìn)的螺旋CT精確濾波 反投影算法],提出了一種基于螺旋線的錐形束CT精確濾波反投影重建方法。具 體的實(shí)現(xiàn)方法是讓X射線源和檢測(cè)器同時(shí)沿著螺旋線運(yùn)動(dòng)，對(duì)檢測(cè)器上獲得的投 影數(shù)據(jù)先進(jìn)行濾波，然后沿著PI線進(jìn)行反投影重建。所述PI線是指螺旋線上兩 點(diǎn)的連線，表示這兩點(diǎn)位置的角度參數(shù)之差小于360° 。這個(gè)方法的優(yōu)點(diǎn)是成像 速度快，成像精度高，垂直軸Z覆蓋范圍大且成像精度與螺距無(wú)關(guān)；這個(gè)方法的 缺點(diǎn)是反投影重建必須沿著PI線進(jìn)行，重建的圖像還需要進(jìn)行坐標(biāo)變換才能被 普通人所識(shí)別，而且這個(gè)方法只適合于螺旋錐形束CT，不能用于沿著馬鞍線掃 描的錐形束CT。

發(fā)明內(nèi)容
本發(fā)明的目的在于克服現(xiàn)有技術(shù)中的不足，提出一種基于螺旋馬鞍線的錐形 束CT精確濾波反投影重建方法，使其具有成像速度快，成像精度高，重建的圖 像不需要經(jīng)過(guò)坐標(biāo)變換就可以直接顯示，而且既可適用于長(zhǎng)物體的成像，也適用 于動(dòng)態(tài)物體的快速成像和連續(xù)觀察。
本發(fā)明是通過(guò)以下技術(shù)方法實(shí)現(xiàn)的，本發(fā)明采用單個(gè)X射線源，首先X射線 源以螺旋馬鞍線路徑進(jìn)行掃描，然后選擇欲重建的平面，計(jì)算該平面與螺旋馬鞍 線掃描路徑的交點(diǎn)，選擇滿足一定條件的交點(diǎn)，確定濾波的方向，對(duì)掃描獲得的 投影數(shù)據(jù)進(jìn)行濾波，最后對(duì)濾波過(guò)的數(shù)據(jù)進(jìn)行反投影重建，獲得重建的圖像。
本發(fā)明方法包括以下步驟
第一步，X射線源在一點(diǎn)發(fā)出錐形束X射線，穿過(guò)被測(cè)對(duì)象后，衰減的x
射線被對(duì)側(cè)的檢測(cè)器檢測(cè)。
第二步，x射線源與檢測(cè)器中心的連線垂直于垂直軸z。
第三步，x射線源與檢測(cè)器同時(shí)繞垂直軸z轉(zhuǎn)動(dòng)，相對(duì)于被測(cè)物體作螺旋馬
鞍線運(yùn)動(dòng)。檢測(cè)器及相對(duì)應(yīng)的投影數(shù)據(jù)采集系統(tǒng)采集投影數(shù)據(jù)。所選螺旋馬鞍線
的參數(shù)要保證任一經(jīng)過(guò)被測(cè)物體并且垂直于垂直軸z的平面與螺旋馬鞍線的交
點(diǎn)不小于4個(gè)。
所述的螺旋馬鞍線，是指螺旋線和馬鞍線在垂直軸Z方向上的線性疊加， 即滿足=cosO)，sin(",^0)COS(2" +力的函數(shù)曲線，其中/7是
螺旋馬鞍線上任一點(diǎn)的三維坐標(biāo)，s是角度參數(shù)，i ("是螺旋馬鞍線上角度參數(shù) 為s的那一點(diǎn)在X-Y平面(垂直于垂直軸Z的平面)上的投影到原點(diǎn)間的距離， A"是馬鞍線的振幅，//("是螺旋線的螺距，s的范圍超過(guò)2;r，及("要取得足
夠大使得螺旋馬鞍線所在的圓柱面包圍被測(cè)物體，A"、 //("為非負(fù)的實(shí)數(shù)。
所述的馬鞍線的振幅，是指馬鞍線在垂直軸Z方向偏離X-Y平面的距離。 所述的螺旋線的螺距，是指螺旋線在垂直軸Z方向上相鄰兩個(gè)螺圈間的距離。
所述的投影數(shù)據(jù)采集系統(tǒng)是指把檢測(cè)器上的數(shù)據(jù)記錄下來(lái)并轉(zhuǎn)化成數(shù)字信 號(hào)傳給計(jì)算機(jī)處理的系統(tǒng)。
第四步，選擇一經(jīng)過(guò)被測(cè)物體并且垂直于垂直軸Z的平面，計(jì)算該平面與螺
旋馬鞍線的交點(diǎn)，記為Pl， P2， P3,…，PN， N〉=4。
第五步，在這N個(gè)交點(diǎn)中，選擇4個(gè)滿足如下條件的交點(diǎn)交點(diǎn)之間連線所 圍成的區(qū)域必須覆蓋被測(cè)物體的斷層層面，另一個(gè)是交點(diǎn)之間角度參數(shù)的差必須 小于360° 。按照這4個(gè)交點(diǎn)角度參數(shù)的大小，分別記為PA， PB, PC， PD。
第六步，這4個(gè)交點(diǎn)將螺旋馬鞍線分為5段，分別是角度參數(shù)小于PA的 部分，角度參數(shù)在PA與PB之間的部分，角度參數(shù)在PB與PC之間的部分，角度 參數(shù)在PC與PD之間的部分和角度參數(shù)大于PD的部分。其中位于PA與PB之間 的投影數(shù)據(jù)的濾波方向是以PA為起點(diǎn)，PB為終點(diǎn)的向量；位于PB與PC之間的 投影數(shù)據(jù)的濾波方向是以PB為起點(diǎn)，PC為終點(diǎn)的向量；位于PC與TO之間的投 影數(shù)據(jù)的濾波方向是以PC為起點(diǎn)，PD為終點(diǎn)的向量；位于PA與PD之間的投影 數(shù)據(jù)的濾波方向是以PD為起點(diǎn)，PA為終點(diǎn)的向量。
第七步，在螺旋馬鞍線上，從點(diǎn)PA開(kāi)始的每一個(gè)點(diǎn)P，存在一系列的平面 經(jīng)過(guò)點(diǎn)P并且平行于點(diǎn)P所對(duì)應(yīng)的濾波方向，做這些平面與檢測(cè)器平面的交線， 沿著這些交線對(duì)檢測(cè)器上獲得的投影數(shù)據(jù)進(jìn)行濾波希爾伯特變換運(yùn)算。
第八步，對(duì)濾波過(guò)的數(shù)據(jù)進(jìn)行反投影重建，得到直角坐標(biāo)系下的重建的二維 斷層圖像。
所述的對(duì)濾波過(guò)的數(shù)據(jù)進(jìn)行反投影重建，是指選擇角度參數(shù)最小的交點(diǎn)作
為積分下限，角度參數(shù)最大的交點(diǎn)作為積分上限，對(duì)濾波后的數(shù)據(jù)進(jìn)行積分的過(guò) 程。
第九步，改變欲重建的平面的位置，重復(fù)第五步到第六步，得到一系列的二 維斷層圖像。
第十步，用這一系列直角坐標(biāo)系下的二維重建圖像，生成直角坐標(biāo)系下的被 測(cè)對(duì)象的三維圖像，在指定的范圍內(nèi)顯示出重建后的高分辨率的三維圖像或二維 斷層圖像。
本發(fā)明原理是對(duì)投影數(shù)據(jù)進(jìn)行濾波時(shí)，濾波的方向在同一個(gè)平面內(nèi)，而且 各個(gè)濾波方向首尾相接，矢量和是0，因此可以直接重建該平面；而基于螺旋線 的精確濾波反投影重建方法對(duì)投影數(shù)據(jù)進(jìn)行濾波時(shí)，濾波的方向由PI線決定， 重建時(shí)也必須沿著PI線，因此得到的是在PI坐標(biāo)系下的重建圖像，還需要進(jìn)行 坐標(biāo)轉(zhuǎn)換才能得到直角坐標(biāo)系下的重建圖像。
本發(fā)明相對(duì)于基于螺旋線的錐形束CT精確濾波反投影重建方法的優(yōu)勢(shì)在 于(1)重建的圖像不需要經(jīng)過(guò)坐標(biāo)變換就可以直接顯示，時(shí)間耗費(fèi)減少30%以 上；(2)掃描方式靈活，通過(guò)改變馬鞍線的振幅和螺旋線的螺距，使這種重建方 法既可適用于長(zhǎng)物體的成像，也適用于動(dòng)態(tài)物體的快速成像和連續(xù)觀察。


圖1為本發(fā)明實(shí)施例中投影數(shù)據(jù)濾波方向的示意圖
具體實(shí)施例方式
下面結(jié)合附圖對(duì)本發(fā)明的實(shí)施例作詳細(xì)說(shuō)明本實(shí)施例在以本發(fā)明技術(shù)方案
為前提下進(jìn)行實(shí)施，給出了詳細(xì)的實(shí)施方式和具體的操作過(guò)程，但本發(fā)明的保護(hù) 范圍不限于下述的實(shí)施例。
實(shí)施例按照以下步驟實(shí)施
1. 把一個(gè)生物樣本置于可沿垂直軸Z方向作勻速直線運(yùn)動(dòng)床上，X射線源
在一點(diǎn)發(fā)出錐形束x射線，經(jīng)準(zhǔn)直器，穿過(guò)活鼠后，衰減的x射線被對(duì)側(cè)的平
板檢測(cè)器檢測(cè)，檢測(cè)器為長(zhǎng)方形的平面，由512x128個(gè)檢測(cè)單元組成。
2. X射線源距離垂直軸Z75cm，檢測(cè)器距離垂直軸Z60cm， X射線源與檢測(cè) 器的中心垂直于垂直軸Z。
3. X射線源和檢測(cè)器繞垂直軸Z沿著馬鞍線做勻速運(yùn)動(dòng)，每隔0.3。采集一 次投影數(shù)據(jù)。生物樣本隨著床沿垂直軸Z做直線運(yùn)動(dòng)。初始時(shí)馬鞍線的振幅為 10cm，螺旋線的螺距為10cm。在記錄投影數(shù)據(jù)時(shí)，還應(yīng)記錄掃描裝置的空間位 置以及其他在之后圖像重建中必要的相關(guān)信息，包括X射線發(fā)射源的垂直軸Z 位置；獲取數(shù)據(jù)時(shí)每個(gè)檢測(cè)器組中的每個(gè)檢測(cè)器單元的位置等。
4. 選擇經(jīng)過(guò)生物樣本并與垂直軸Z垂直的平面Z=0，計(jì)算它與螺旋馬鞍線 的交點(diǎn)，記為P1，P2，P3，…，PN， N>=4。
5. 在這N個(gè)交點(diǎn)中，選擇4個(gè)點(diǎn)，按照這4個(gè)交點(diǎn)角度參數(shù)的大小，分別 記為PA,PB，PC，PD。這4個(gè)點(diǎn)要滿足以下條件4個(gè)點(diǎn)圍成的四邊形區(qū)域必須覆 蓋生物樣本位于Z=0的斷層層面，PD與PA角度參數(shù)的差必須小于360° 。
6. 這4個(gè)交點(diǎn)將螺旋馬鞍線分為5段，分別是角度參數(shù)小于PA的部分，角 度參數(shù)在PA與PB之間的部分，角度參數(shù)在PB與PC之間的部分，角度參數(shù)在PC 與PD之間的部分和角度參數(shù)大于PD的部分。其中位于PA與PB之間的投影數(shù)據(jù)
的濾波方向是以PA為起點(diǎn)，PB為終點(diǎn)的向量；位于PB與PC之間的投影數(shù)據(jù)的 濾波方向是以PB為起點(diǎn)，PC為終點(diǎn)的向量；位于PC與PD之間的投影數(shù)據(jù)的濾 波方向是以PC為起點(diǎn)，PD為終點(diǎn)的向量；位于PA與PD之間的投影數(shù)據(jù)的濾波 方向是以PD為起點(diǎn)，PA為終點(diǎn)的向量。圖1中實(shí)線就是X射線源與檢測(cè)器掃描 時(shí)的螺旋馬鞍線軌跡，箭頭所指的方向即每一部分的濾波方向，
7. 在螺旋馬鞍線上，從點(diǎn)PA開(kāi)始，角度參數(shù)每增加0.3度，就選擇一個(gè)點(diǎn) P，同時(shí)在垂直軸Z上選定一個(gè)點(diǎn)Q，點(diǎn)P與點(diǎn)Q在垂直軸Z方向上的坐標(biāo)相同。 判斷點(diǎn)Q與欲重建的平面Z二O的關(guān)系(1)若點(diǎn)Q在平面Z^之上，則將點(diǎn)Q 往下移動(dòng)，每次移動(dòng)lmm，作經(jīng)過(guò)點(diǎn)P、點(diǎn)Q并且平行于點(diǎn)P所對(duì)應(yīng)的濾波方向 的平面W，平面W與平板檢測(cè)器有交線L，沿著交線L對(duì)檢測(cè)器上獲得的投影數(shù) 據(jù)進(jìn)行希爾伯特變換， 一直移動(dòng)點(diǎn)Q，直到平面W與平板檢測(cè)器沒(méi)有交線為止；
(2)若點(diǎn)Q在平面Z:O之下，則將點(diǎn)Q往上移動(dòng)，每次移動(dòng)lmm，作經(jīng)過(guò)點(diǎn)P、 點(diǎn)Q并且平行于點(diǎn)P所對(duì)應(yīng)的濾波方向的平面W，平面W與平板檢測(cè)器有交線L， 沿著交線L對(duì)檢測(cè)器上獲得的投影數(shù)據(jù)進(jìn)行濾波希爾伯特變換， 一直移動(dòng)點(diǎn)Q， 直到平面W與平板檢測(cè)器沒(méi)有交線為止。按照如上的方法分別對(duì)從點(diǎn)PA到點(diǎn)PB、 點(diǎn)PB到點(diǎn)PC、點(diǎn)PC到點(diǎn)PD、點(diǎn)PA到點(diǎn)PD間的投影數(shù)據(jù)進(jìn)行濾波。
8. 對(duì)濾波后的投影數(shù)據(jù)進(jìn)行反投影，重建出直角坐標(biāo)系下生物樣本斷層 512X512點(diǎn)的圖像。具體地說(shuō)，就是先計(jì)算4個(gè)積分，第一個(gè)積分以點(diǎn)PA為積 分下限、點(diǎn)PB為積分上限，第二個(gè)積分以點(diǎn)PB為積分下限、點(diǎn)PC為積分上限， 第三個(gè)積分以點(diǎn)PC為積分下限、點(diǎn)PD為積分上限，第四個(gè)積分以點(diǎn)PD為積分 下限、點(diǎn)PA為積分上限，這4個(gè)積分都沿著螺旋馬鞍線對(duì)濾波后的投影數(shù)據(jù)進(jìn) 行積分運(yùn)算，然后將這4個(gè)積分線性疊加起來(lái)即獲得重建的體數(shù)據(jù)。
9. 重建的圖像通過(guò)可視化軟件顯示在高分辨率的醫(yī)用顯示器上。 這樣的實(shí)施方式，避免了A. Katsevich精確重建方法中需要將體數(shù)據(jù)重建
完全，然后將體數(shù)據(jù)從PI坐標(biāo)系變換到笛卡兒坐標(biāo)系后，重建圖像才能顯示的 缺陷，不需要對(duì)重建圖像的進(jìn)行坐標(biāo)變換即可直接顯示，大大加快了成像速度， 實(shí)施效果是時(shí)間耗費(fèi)降低了 45%。
權(quán)利要求
1、一種基于螺旋馬鞍線的錐形束CT精確濾波反投影重建方法，其特征在于，包括以下步驟第一步，X射線源在一點(diǎn)發(fā)出錐形束X射線，穿過(guò)被測(cè)對(duì)象后，衰減的X射線被對(duì)側(cè)的檢測(cè)器檢測(cè)；第二步，X射線源與檢測(cè)器中心的連線垂直于垂直軸Z；第三步，X射線源與檢測(cè)器同時(shí)繞垂直軸Z轉(zhuǎn)動(dòng)，相對(duì)于被測(cè)物體作螺旋馬鞍線運(yùn)動(dòng)，檢測(cè)器及相對(duì)應(yīng)的投影數(shù)據(jù)采集系統(tǒng)采集投影數(shù)據(jù)，所選螺旋馬鞍線的參數(shù)要保證任一經(jīng)過(guò)被測(cè)物體并且垂直于垂直軸Z的平面與螺旋馬鞍線的交點(diǎn)大于或等于4個(gè)；第四步，選擇一經(jīng)過(guò)被測(cè)物體并且垂直于垂直軸Z的平面，計(jì)算該平面與螺旋馬鞍線的交點(diǎn)，記為P1，P2，P3，…，PN，N＞＝4；第五步，在這N個(gè)交點(diǎn)中，選擇4個(gè)滿足如下條件的交點(diǎn)交點(diǎn)之間連線所圍成的區(qū)域必須覆蓋被測(cè)物體的斷層層面，另一個(gè)是交點(diǎn)之間角度參數(shù)的差必須小于360°，按照這4個(gè)交點(diǎn)角度參數(shù)的大小，分別記為PA，PB，PC，PD；第六步，這4個(gè)交點(diǎn)將螺旋馬鞍線分為5段，分別是角度參數(shù)小于PA的部分，角度參數(shù)在PA與PB之間的部分，角度參數(shù)在PB與PC之間的部分，角度參數(shù)在PC與PD之間的部分和角度參數(shù)大于PD的部分；第七步，在螺旋馬鞍線上，從點(diǎn)角度參數(shù)PA開(kāi)始的每一個(gè)點(diǎn)P，存在一系列的平面經(jīng)過(guò)點(diǎn)P并且平行于點(diǎn)P所對(duì)應(yīng)的濾波方向，做這些平面與檢測(cè)器平面的交線，沿著這些交線對(duì)檢測(cè)器上獲得的投影數(shù)據(jù)進(jìn)行濾波希爾伯特變換運(yùn)算；第八步，對(duì)濾波過(guò)的數(shù)據(jù)進(jìn)行反投影重建，得到直角坐標(biāo)系下的重建的二維斷層圖像；第九步，改變欲重建的平面的位置，重復(fù)第五步到第六步，得到一系列的二維斷層圖像；第十步，用這一系列直角坐標(biāo)系下的二維重建圖像，生成直角坐標(biāo)系下的被測(cè)對(duì)象的三維圖像，在指定的范圍內(nèi)顯示出重建后的三維圖像或二維斷層圖像。
2、 根據(jù)權(quán)利要求1所述的基于螺旋馬鞍線的錐形束CT精確濾波反投影重建 方法，其特征是，所述的螺旋馬鞍線，是指螺旋線和馬鞍線在垂直軸Z方向上的線性疊加，即滿足= cos(",及O)sin(",」(5)cos(2s) + 的函數(shù)曲線，其中p是螺旋馬鞍線上任一點(diǎn)的三維坐標(biāo)，s是角度參數(shù)，/ 00是螺旋馬鞍線上角度參數(shù)為的那一點(diǎn)在垂直于垂直軸Z的X-Y平面上的投影到原點(diǎn)間的距 離，J(力是馬鞍線的振幅，//("是螺旋線的螺距，s的范圍超過(guò)2;r， i (力要取得足夠大使得螺旋馬鞍線所在的圓柱面包圍被測(cè)物體，A"、 //(力為非負(fù)的實(shí)數(shù)。
3、 根據(jù)權(quán)利要求1所述的基于螺旋馬鞍線的錐形束CT精確濾波反投影重建 方法，其特征是，所述的馬鞍線的振幅，是指馬鞍線在垂直軸Z方向偏離垂直 于垂直軸Z的X-Y平面的距離。
4、 根據(jù)權(quán)利要求1所述的基于螺旋馬鞍線的錐形束CT精確濾波反投影重建 方法，其特征是，所述的螺旋線的螺距，是指螺旋線在垂直軸Z方向上相鄰兩 個(gè)螺圈間的距離。
5、 根據(jù)權(quán)利要求1所述的基于螺旋馬鞍線的錐形束CT精確濾波反投影重建 方法，其特征是，所述的投影數(shù)據(jù)采集系統(tǒng)是指把檢測(cè)器上的數(shù)據(jù)記錄下來(lái)并轉(zhuǎn) 化成數(shù)字信號(hào)傳給計(jì)算機(jī)處理的系統(tǒng)。
6、 根據(jù)權(quán)利要求1所述的基于螺旋馬鞍線的錐形束CT精確濾波反投影重建 方法，其特征是，所述4個(gè)交點(diǎn)將螺旋馬鞍線分為5段，其中位于PA與PB 之間的投影數(shù)據(jù)的濾波方向是以PA為起點(diǎn)，PB為終點(diǎn)的向量；位于PB與PC之 間的投影數(shù)據(jù)的濾波方向是以PB為起點(diǎn)，PC為終點(diǎn)的向量；位于PC與PD之間 的投影數(shù)據(jù)的濾波方向是以PC為起點(diǎn)，PD為終點(diǎn)的向量；位于PA與PD之間的 投影數(shù)據(jù)的濾波方向是以PD為起點(diǎn)，PA為終點(diǎn)的向量。
7、 根據(jù)權(quán)利要求1所述的基于螺旋馬鞍線的錐形束CT精確濾波反投影重建 方法，其特征是，所述的對(duì)濾波過(guò)的數(shù)據(jù)進(jìn)行反投影重建，是指選擇角度參數(shù) 最小的交點(diǎn)作為積分下限，角度參數(shù)最大的交點(diǎn)作為積分上限，對(duì)濾波后的數(shù)據(jù) 進(jìn)行積分的過(guò)程。
全文摘要
本發(fā)明涉及的是一種生物醫(yī)學(xué)成像技術(shù)領(lǐng)域的基于螺旋馬鞍線的錐形束CT精確濾波反投影重建方法。本發(fā)明采用單個(gè)X射線源，首先X射線源以螺旋馬鞍線路徑進(jìn)行掃描，然后選擇欲重建的平面，計(jì)算該平面與螺旋馬鞍線掃描路徑的交點(diǎn)，選擇交點(diǎn)，確定濾波的方向，對(duì)掃描獲得的投影數(shù)據(jù)進(jìn)行濾波，最后對(duì)濾波過(guò)的數(shù)據(jù)進(jìn)行反投影重建，獲得重建的圖像。本發(fā)明具有成像速度快，成像精度高，重建的圖像不需要經(jīng)過(guò)坐標(biāo)變換就可以直接顯示，而且既可適用于長(zhǎng)物體的成像，也適用于動(dòng)態(tài)物體的快速成像和連續(xù)觀察。
文檔編號(hào)A61B6/03GK101194844SQ20071017327
公開(kāi)日2008年6月11日 申請(qǐng)日期2007年12月27日 優(yōu)先權(quán)日2007年12月27日
發(fā)明者楊 呂, 俊 趙 申請(qǐng)人:上海交通大學(xué)