多層容積ct灌注成像源影像的自動化減影處理方法
【專利摘要】本發明公開了一種多層容積CT灌注成像源影像的自動化減影處理方法，該方法首先根據CT灌注影像序列計算CT成像設備的球管每轉一周所掃描的層數n，并確定掃描時間中每一掃描時刻容積成像數據的n幅圖像，獲取每一時刻所有圖像的容積CT值集合，然后確定CT減影的被減時刻，由被減時刻后的容積CT值集合減去被減時刻的容積CT值集合得到減影容積CT值集合，最后根據減影容積CT值集合與其對應文件元信息集合，獲得減影DICOM文件。通過該方法實現容積CTP源影像的自動化減影處理與成像，提高了減影的效率與可靠性，為形成一站式的減影CTP軟件奠定了基礎。
【專利說明】多層容積CT灌注成像源影像的自動化減影處理方法
【技術領域】
[0001]本發明涉及容積減影成像的計算機自動化處理【技術領域】，具體涉及一種多層容積CT灌注成像源影像的自動化減影處理方法。
【背景技術】
[0002]CT灌注成像(CT Perfusion，CTP)是在靜脈注射對比劑的同時，對感興趣區層面或容積層塊進行連續動態CT掃描，以獲得所選層面或容積層塊內每一個像素的時間一密度曲線(Time Density Curve, TDC),并根據此曲線通過不同的數學模型轉換和計算機偽彩處理得到局部組織的血流量(Blood Flow, BF)、血容量(Blood Volume, BV)、平均通過時間(Mean Transit Time, MTT)、達峰時間(Time to Peak, TTP)和表面通透性(permeabilitysurface area, PS)等血流動力學參數和灌注偽彩圖。
[0003]CTP可以從形態學、功能代謝諸方面提供疾病的診斷信息，但若檢查部位組織在CT圖像上表現出高密度影(如骨組織、血腫、鈣化、成骨性病變等)，就限制了 CTP技術的應用。這是因為高密度影本身，往往掩蓋了對比劑血流灌注前、后局部組織密度及微循環血流動力學的變化，使得CTP難以甚至無法提供診病信息。近年有文獻報道，使用減影CTP(Subtraction CT Perfusion, SCTP)技術從軟件的圖像后處理角度能夠有效地克服因組織本身高密度影對CTP診斷的 影響[1，2]。但SCTP技術沒有整合，后處理時間過長，自動化程度低，限制了其臨床推廣應用，需要研發一站式全自動計算機分析系統。SCTP —站式全自動計算機分析系統需要解決兩個關鍵問題，一是CTP源影像容積數據的全自動計算機減影，生成SCTP源影像；二是找到合適數學模型，根據層面或容積層塊內圖像減運算后的每一個像素的TDC，計算機自動轉換和生成局部組織BF、BV、MTT、TTP和PS等血流動力學參數值和灌注圖像。
[0004]圖像的減運算(即減影技術)在醫學上已有廣泛的應用，典型技術莫過于數字減影血管造影(Digital Subtraction Angiography, DSA) [3],近些年來 CT、MRI 領域也有涉獵。Hanna等[4]使用減影技術對肌肉與骨骼腫瘤MR灌注影像序列的對比劑前、后的的圖像進行了減影處理，發現血管惡性組織、脂肪和出血在減影后的圖像上對比度得到了增強。Terai等[5]使用動態對比劑增強減影磁共振方法診斷睪丸扭轉。李征宇等[6]使用MR動態增強時間減影灌注成像技術對周圍型肺癌肺動脈血供定性與定量評估研究時在工作站上將強化前后的圖像行減影處理。Huang和Chamberli [7]使用計算機化的自動方法對CT灌注影像序列的對比劑增強前后的圖像進行減處理，在減影后的圖像上發現了對比劑注入前和注入后的CT圖像上肉眼觀察不到的病灶。但這些CT、MRI應用與已報道的SCTP —樣，大多基于影像工作站上的減影軟件[1，2，4-6]對單層圖像(未涉及容積層塊圖像)進行半自動的減影處理，操作繁瑣，費時且容易出錯。
【發明內容】
[0005]針對現有技術中存在的缺陷，本發明的目的在于提供一種多層容積CT灌注成像源影像的自動化減影處理方法，通過該方法實現容積CTP源影像的自動化減影處理與成像。
[0006]為實現上述目的，本發明采用的技術方案如下:
[0007]—種多層容積CT灌注成像源影像的自動化減影處理方法，包括以下步驟:
[0008](1)進行容積層CT灌注成像掃描，獲取容積成像數據，根據掃描時間和容積成像數據計算CT成像設備的球管每轉一周所掃描的層數η;所述的容積成像數據包括掃描得到的CT灌注影像序列中的DICOM圖像；
[0009](2)確定掃描時間中每一掃描時刻容積成像數據的η幅圖像，獲取每一掃描時刻η幅圖像的容積CT值集合；
[0010](3)確定CT減影的被減時刻，根據被減時刻的容積CT值集合得到其它掃描時刻的減影容積CT值集合；
[0011](4)根據減影容積CT值集合與其對應的文件元信息集合，獲取容積成像數據的減影DICOM文件。
[0012]進一步，如上所述的一種多層容積CT灌注成像源影像的自動化減影處理方法，步驟(I)中，計算CT成像設備的球管每轉一周所掃描的層數η的具體方式為:
[0013](1-1)計算CT灌注影像序列中所有圖像的圖像平面與原點的最小距離，并存儲在一個一維數組中；
[0014](1-2)刪除所述一維數組中的重復的元素，數組中剩余元素的個數為CT成像設備的球管每轉一周所掃描的層數η。
[0015]進一步，如上所述的一種多層容積CT灌注成像源影像的自動化減影處理方法，步驟(2)中，獲取每一時刻所有圖像的容積CT值集合的具體方式為:
[0016](2-1)根據所獲取的CT灌注影像序列中DICOM圖像的圖像序列號InstanceNunber確定某一掃描時刻timei容積成像數據的η幅圖像；
[0017](2-2)根據η幅圖的圖像序列號Instance Nunber的順序依次讀入對應圖像的像素數值，并轉換成對應圖像的CT值，得到η個一維CT值數組；
[0018](2-3)將所得到的η個一維CT值數組按照其對應圖像的圖像序列號InstanceNunber串接得到timei時刻的容積CT集合。
[0019]進一步，如上所述的一種多層容積CT灌注成像源影像的自動化減影處理方法，步驟(3)中，所述的其它掃描時刻的減影容積CT值集合包括被減時刻之前的掃描時刻的減影容積CT值集合和被減時刻之后的掃描時刻的減影容積CT值集合。
[0020]進一步，如上所述的一種多層容積CT灌注成像源影像的自動化減影處理方法，步驟(3)中，獲取被減時刻之后的掃描時刻的減影容積CT值集合的具體方式為:
[0021](3-1)確定被減時刻subtractedi,獲取subtractedi時刻的容積CT集合HUsetO ；
[0022](3-2)獲取被減時刻subtractedi之后的某一掃描時刻t的容積容積CT集合HUset ；
[0023](3-3)計算某一掃描時刻t的減影容積CT值集合deltaHUset,計算公式為:deltaHUset=HUset_HUsetO。
[0024]進一步，如上所述的一種多層容積CT灌注成像源影像的自動化減影處理方法，步驟(3)中，被減時刻subtractedi之前的所有掃描時刻的減影容積CT值集合等于subtractedi+1時刻的減影容積CT集合。
[0025]再進一步，如上所述的一種多層容積CT灌注成像源影像的自動化減影處理方法，步驟(3)中，如果減影容積CT值集合中某一元素的數值小于0，將該元素的數值置為O。
[0026]更進一步，如上所述的一種多層容積CT灌注成像源影像的自動化減影處理方法，步驟(4)中，獲取容積成像數據的減影DICOM文件的具體方式為:
[0027]假設需要獲取t時刻的容積成像數據中第j幅圖像的DICOM文件，首先獲取t時刻的減影容積CT值集合deltaHUset,對集合deltaHUset中的第j個元素deltaHUset [j]進行線性變換得到第j幅圖像的像素值，根據轉換得到的第j幅圖像的像素值與t時刻容積成像數據中第j幅圖像的DICOM文件元信息得到subtractedi時刻容積成像數據第j幅圖像的DICOM文件；其中，I≤j≤η。
[0028]本發明的有益效果在于:本發明所述的方法實現了容積灌注源影像的減影CTP源影像的全自動生成，為形成一站式的減影CTP軟件奠定了基礎。該方法全自動減影，速度快，不會因為人為的因素誤操作，而生成錯誤的圖像序列，提高了減影的效率與可靠性。
【專利附圖】
【附圖說明】
[0029]圖1為本發明一種多層容積CT灌注成像源影像的自動化減影處理方法的流程圖； [0030]圖2為【具體實施方式】中計算某一時刻的減影容積CT值集合的流程圖。
【具體實施方式】
[0031]下面結合說明書附圖與【具體實施方式】對本發明做進一步的詳細說明。
[0032]圖1示出了本發明一種多層容積CT灌注成像源影像的自動化減影處理方法的流程圖，該方法主要包括以下步驟:
[0033]步驟Sll:根據CT灌注影像序列計算CT成像設備的球管每轉一周所掃描的層數；
[0034]進行容積層CT灌注成像掃描，獲取容積成像數據,根據掃描時間timenumber和掃描得到的容積成像數據計算CT成像設備的球管每轉一周所掃描的層數η;所述的容積成像數據包括掃描得到的CT灌注影像序列中的DICOM圖像。
[0035]CT灌注影像序列中的圖像的存儲格式為DICOM文件格式，DICOM文件一般由文件兀信息(File Meta Information)和DICOM數據集合(DICOM Data Set)組成。數據兀素(Data Element)是DICOM格式文件最基本的組成單元。它由4個部分組成:標簽(Tag)、數據描述 VR (Value Representation)、數據長度 VL (Value Lengh)以及數據域(Value)。使用數據元素標簽Tag為(0020,0032)和(0020,0037)的數值就可以計算出圖像平面與原點間的最小距離，用以表示圖像位置。在CT灌注影像序列中，兩幅圖像的圖像平面與原點間的最小距離相同，就認為這兩幅圖像處于同一個層面。計算CT灌注影像序列中所有圖像的圖像平面與原點間的最小距離，并存儲在I個一維數組中。剔除這個數組中重復的元素，此時剩下元素的個數η就是CT成像設備的球管每轉一周獲得的層數。而這個影像序列的數據掃描時間timenumber是影像序列的總共擁有的圖像數除以這個層數η。
[0036]步驟S12:根據每一掃描時刻容積成像數據的圖像獲取每一時刻所有圖像的容積CT值集合；[0037]確定掃描時間中每一掃描時刻容積成像數據的η幅圖像，獲取每一掃描時刻所得到的η幅圖像的容積CT值集合。若多幅DICOM圖像同屬于一個醫學圖像序列，那么它們一般有一個唯一標識的Series Instance UID(0020, 000E)(圖像序列的唯一標識符)和StudyInstance UID (0020, 000D)(研究的唯一標識符)，可以根據每幅DICOM圖像的圖像序列號Instance Number (0020, 0013)(標識圖象的號碼)確定某一掃描時刻容積成像數據是由哪幾幅圖像組成的。例如，進行8層CT灌注成像掃描時，即CT成像設備的球管每轉一周同時掃描8層，假設數據采集(掃描時間)30s產生240幅圖像，則每一掃描時刻的容積成像數據由8幅圖像組成，而掃描得到的每一層的圖像也可以確定，設獲取到的240幅DICOM的文件名即圖像序列號為頂1，IM2,……，頂239，頂240，因為一共是8層，每一層30幅圖像(文件)，所以每一層包含的圖像(文件)情況如下:
[0038]第I 層:頂1，頂9，......, IM233 ；
[0039]第2 層:IM2, IM10,......，頂234 ；
[0040]第3 層:IM3, IM11,......，頂235 ；
[0041]第
4 層:IM4, IM12,......，頂236 ；
[0042]第5 層:IM5, IM13,......，頂237 ；
[0043]第6 層:IM6, IM14,......，頂238 ；
[0044]第7 層:IM7, IM15,......，頂239 ；
[0045]第8 層:IM8, IM16,......，頂240 ；
[0046]共計240幅圖像(文件)。
[0047]由步驟Sll可知，CT成像設備的球管每轉一周所掃描的層數n，因此每一掃描時刻的容積成像數據都有η幅圖像，如第I秒獲得的容積圖像包括第I幅圖像到第η幅圖像，第I秒獲得的容積圖像包括第η+1幅圖像到第2η幅圖像。第i秒獲得的容積圖像包括第i*n+l幅圖像到第(i+l)*n幅圖像，這樣，第i秒獲得的容積的CT值集合是依次讀取第i*n+l幅圖像到第(i+1) *n幅圖像的CT值,形成η個I維數組,并按圖像對應的Instance Number從小到大的順序串接成一 個一維數組。計算容積CT值集合的具體方式如下:
[0048]例如，第timei時刻容積成像數據由η幅圖像組成，其中timei≥1且timei≤timenumber。這η幅圖像的Instance Number形成一個等差數列。這個等差數列的起始值是(time1-l)Xn+l，而公差是1，一共有η項。其中η是CT成像設備的球管每轉一周獲得的層數。確定了某掃描時刻容積成像數據是由哪幾幅圖像組成后，就可以根據這些圖像Instance Number從小到大的順序讀入對應圖像的像素數值,再變換成CT值。每讀入一幅圖像就得到一個一維CT值數組，這樣就可以得到η個一維CT值數組。把這η個一維CT值數組按照對應圖像的Instance Number從小到大的順序串接成一個一維數組。這個一維數組存儲的就是第timei時刻容積CT值集合。
[0049]步驟S13:計算CT減影容積CT值集合；
[0050]確定CT減影的被減時刻，根據被減時刻的容積CT值集合得到其它掃描時刻的減影容積CT值集合。假設被減時刻記為subtractedi,被減時刻subtractedi之后的某一掃描時刻記為t，計算被減時刻之外的其它掃描時刻的減影容積CT值集合包括被減時刻之前的掃描時刻的減影容積CT值集合和被減時刻之后的掃描時刻的減影容積CT值集合。本實施方式中計算被減時刻之后某一掃描時刻t的減影容積CT值集合的具體方式如圖2所示:[0051]首先確定被減時刻subtractedi,然后使用步驟S12中計算容積CT值集合的方法獲得被減時刻的CT值集合HUsetO，再使用步驟S12中計算出的第t時刻的容積CT值集合HUset，然后用HUset減去HUsetO得到t時刻的減影容積CT值集合deltaHUset，即被減時刻之后的t時刻的減影容積CT值集合deltaHUset = HUset_HUsetO。按此方法,可以計算出大于subtractedi的所有掃描時刻(被減時刻之后)的減影容積CT值集合。而小于等于被減時刻subtractedi的所有掃描時刻的減影容積CT值集合等于第subtractedi+Ι時刻的減影容積CT值集合。通過上述計算即可得到所有掃描時刻的減影容積CT值集合。由于deltaHUset、HUset 和 HUsetO 都是一維數組，那么:
[0052]deltaHUset[i] = HUset[i]-HUsetO[i]
[0053]其中，deltaHUset[i]、HUset [i]和 HUsetO[i]分別表不 deltaHUset、HUset 和HUsetO下標為i的元素的數值。為了便于計算，deltaHUset中小于O的數值都置為O。
[0054]步驟S14:根據減影容積CT值集合與其對應的文件元信息集合，獲取容積成像數據的減影DICOM文件。
[0055]根據減影容積CT值集合與其對應的文件元信息集合，獲取容積成像數據的減影DICOM文件，具體方式如下:
[0056]假設某時刻t的減影容積CT值集合是deltaHUset，那么，該時刻對應的原始容積數據第j幅圖像貢獻的減影CT值數據是一個一維數組deltaHUset [ j],其中1≤ j ≤ η。而deltaHUset [ j]等于 deltaHUset 中從下標(j-1) Xheight X width 到 j Xheight X width 之間的數值，其中，height和width分別是影像序列中任意一幅圖像的高和寬，參見文獻《基于DICOM標準的醫學圖像的顯示方法研究》。再根據這個第j幅圖像的DICOM文件元信息以及對deltaHUset [j]進行線性變換，就可以把deltaHUset [j]存儲為一幅DICOM文件。這樣就可以得到所有時刻容積成像數據對應的減影DICOM文件。
[0057]其中，對是deltaHUset [j]進行線性變換為現有技術，具體變化方法可參見文獻《基于DICOM標準的醫學圖像的顯示方法研究》。DICOM文件格式一般包括兩部分:文件元信息(File Meta Information)和DICOM數據集合(DICOM Data Set)。文件兀信息包含DICOM數據集合的識別信息，可分為3部分:文件前言(File Preamble)、DICOM前綴(DICOMPref ix)和文件兀要素(File Meta Elements), DICOM文件兀信息可以通過Gdcm、DCMTK、CTN、ITK,VTK開源免費下載的DICOM開發工具包獲取，也可以使用Matlab和MITK獲取。將DICOM文件元信息與線性變換后的deltaHUset [j]結合得到DICOM文件也是現有技術，可以使用Gdcm、DCMTK、CTN、ITK、VTK等開源免費下載的DICOM開發工具包可以得到，也可以使用Matlab和MITK得到。例如，使用Matlab的dicomwrite函數即可辦到。把DICOM文件元信息和對deltaHUset [j]線性變換的結果作為dicomwrite函數的輸入參數,就可以得到DICOM文件。
[0058]本發明所述的減影處理方法與現有的半自動手工技術相比，提供了一種在PC機上對容積CT灌注(CTP)源影像全自動生成減影CTP源影像的方法，全自動減影，速度快，減影過程按照指定的程序進行，不會因為人為的因素誤操作，而生成錯誤的圖像序列。因此，減影的可靠性得到了保證。促進SCTP的臨床應用以及理論研究，填補國產醫學影像軟件領域的一些空白點。
[0059]下面結合具體的實施例對本發明所述的方法進行進一步說明。[0060]實施例
[0061]本實施例中CT掃描CT機進行2層CT灌注成像掃描，即球管每轉一周同時掃描2層，數據采集4s產生8幅圖像。每幅圖像的存儲格式為DICOM文件格式。8幅圖像分別記為頂1、IM2-頂8，則第I秒得到的容積成像數據是I Ml和頂2，第2秒得到的容積成像數據是頂3和頂4，第3秒得到的容積成像數據是頂5和頂6，第4秒得到的容積成像數據是IM7和頂8。其中，頂1、頂3、頂5、頂7是第一層隨時間變化得到的數據信息，頂2、頂4、頂6、IM8是第二層隨時間變化得到的數據信息。
[0062]如果令第I秒是被減時刻，那么第I秒掃描得到的兩幅圖像為頂1和頂2，根據2幅圖的Instance Nunber的順序依次讀入IMl和IM2對應的像素數值,并轉換成相應的CT值，得到2個一維CT值數組；讀取頂1得到以為CT數組記為CTHU1，讀取頂2得到以為CT數組記為CTHU2，將CTHUl和CTHU2串接成更大的一維數組VolumeCTHUl得到第I秒的容積CT集合。假設CTHUl有nl個元素，CTHU2有n2個元素，那么VolumeCTHUl有nl+n2個元素。VolumeCTHUl表示第I秒容積成像數據的CT值集合。同理，可以得到第2秒至第4秒的CT值集合 VolumeCTHU2、VolumeCTHU3、VolumeCTHU4。在線性代數中 VolumeCTHUl、VolumeCTHU2、VolumeCTHU3 和 VolumeCTHU4 可以看做向量，用向量 VolumeCTHU2 減去 VolumeCTHUl 得到deltaVolumeCTHU2，用向量 VolumeCTHU3 減去 VolumeCTHUl 得到 deltaVolumeCTHU3，用向量VolumeCTHU4 減去 VolumeCTHUl 得到 deltaVolumeCTHU4。對于 deltaVolumeCTHU2 前 nl 個元素可以組成向量deItaCTHU3，而后n2個元素可以組成向量deItaCTHU4。把deItaCTHU3和deltaCTHU4進行線性變換得到deltaX3和deltaX4。在得到VolumeCTHU2的同時，通過DICOM開發工具包可以得到IM3的文件元信息metadata3和IM4的文件元信息metadata4。使用DICOM開發工具包(例如Matlab的dicomwrite函數)把IM3的文件元信息metadata3和deltaX3作為輸入參數可以得到DICOM文件。使用DICOM開發工具包(例如Matlab的dicomwrite函數)把IM4的文 件兀信息metadata4和deltaX4作為輸入參數可以得到DICOM文件。同理，可以得到2層8幅CT灌注影像序列對應的減影CT灌注影像序列文件。
[0063]顯然，本領域的技術人員可以對本發明進行各種改動和變型而不脫離本發明的精神和范圍。這樣，倘若本發明的這些修改和變型屬于本發明權利要求及其同等技術的范圍之內，則本發明也意圖包含這些改動和變型在內。
【權利要求】
1.一種多層容積CT灌注成像源影像的自動化減影處理方法，包括以下步驟: (1)進行容積層CT灌注成像掃描，獲取容積成像數據，根據掃描時間和容積成像數據計算CT成像設備的球管每轉一周所掃描的層數η;所述的容積成像數據包括掃描得到的CT灌注影像序列中的DICOM圖像； (2)確定掃描時間中每一掃描時刻容積成像數據的η幅圖像，獲取每一掃描時刻η幅圖像的容積CT值集合； (3)確定CT減影的被減時刻，根據被減時刻的容積CT值集合得到其它掃描時刻的減影容積CT值集合； (4)根據減影容積CT值集合與其對應的文件元信息集合，獲取容積成像數據的減影DICOM文件。
2.如權利要求1所述的一種多層容積CT灌注成像源影像的自動化減影處理方法，其特征在于:步驟(1)中，計算CT成像設備的球管每轉一周所掃描的層數η的具體方式為: (1-1)計算CT灌注影像序列中所有圖像的圖像平面與原點的最小距離，并存儲在一個一維數組中； (1-2)刪除所述一維數組中的重復的元素，數組中剩余元素的個數為CT成像設備的球管每轉一周所掃描的層數η。
3.如權利要求2所述的一種多層容積CT灌注成像源影像的自動化減影處理方法，其特征在于:步驟(2)中，獲取每一時刻所有圖像的容積CT值集合的具體方式為: (2-1)根據所獲取的CT灌注影像序列中DICOM圖像的圖像序列號Instance Nunber確定某一掃描時刻timei容積成像數據的η幅圖像； (2-2)根據η幅圖的圖像序列號Instance Nunber的順序依次讀入對應圖像的像素數值，并轉換成對應圖像的CT值，得到η個一維CT值數組； (2-3)將所得到的η個一維CT值數組按照其對應圖像的Instance Nunber串接得到timei時刻的容積CT集合。
4.如權利要求3所述的一種多層容積CT灌注成像源影像的自動化減影處理方法，其特征在于:步驟(3)中，所述的其它掃描時刻的減影容積CT值集合包括被減時刻之前的掃描時刻的減影容積CT值集合和被減時刻之后的掃描時刻的減影容積CT值集合。
5.如權利要求4所述的一種多層容積CT灌注成像源影像的自動化減影處理方法，其特征在于:步驟(3)中，獲取被減時刻之后的掃描時刻的減影容積CT值集合的具體方式為: (3-1)確定被減時刻subtractedi,獲取subtractedi時刻的容積CT集合HUsetO ； (3-2)獲取被減時刻subtractedi之后的某一掃描時刻t的容積容積CT集合HUset ； (3-3)計算某一掃描時刻t的減影容積CT值集合deltaHUset,計算公式為:deltaHUset=HUset_HUsetO。
6.如權利要求5所述的一種多層容積CT灌注成像源影像的自動化減影處理方法，其特征在于:步驟(3)中，被減時刻subtractedi之前的所有掃描時刻的減影容積CT值集合等于subtractedi+Ι時刻的減影容積CT集合。
7.如權利要求6所述的一種多層容積CT灌注成像源影像的自動化減影處理方法，其特征在于:步驟(3)中，如果減影容積CT值集合中某一元素的數值小于0，將該元素的數值置為0。
8.如權利要求7所述的一種多層容積CT灌注成像源影像的自動化減影處理方法，其特征在于:步驟(4)中，獲取容積成像數據的減影DICOM文件的具體方式為: 假設需要獲取t時刻的容積成像數據中第j幅圖像的DICOM文件，首先獲取t時刻的減影容積CT值集合deltaHUset,對集合deltaHUset中的第j個元素deltaHUset [j]進行線性變換得到第j幅圖像的像素值，根據轉換得到的第j幅圖像的像素值與t時刻容積成像數據中第j幅圖像的DICOM文件元信息得到subtractedi時刻容積成像數據第j幅圖像的DICOM文件；其中，1≤j≤η。
【文檔編號】A61B6/03GK103908276SQ201210592536
【公開日】2014年7月9日 申請日期:2012年12月29日 優先權日:2012年12月29日
【發明者】楊朝輝, 楊秀軍, 汪雪林, 佟慶彬, 雷大偉 申請人:北京握奇數據系統有限公司