[發明專利]電鏡三維圖像重構的方法及系統

有效

• 申請號: 201210326514.7
• 申請日: 2012-09-05
• 公開（公告）號: CN102881042A
• 公開（公告）日: 2013-01-16
• 發明（設計）人: 盧曉偉;張清
• 申請（專利權）人: 浪潮（北京）電子信息產業有限公司
• 主分類號: G06T17/00
• 分類號: G06T17/00;G06T1/20
• 代理公司: 北京安信方達知識產權代理有限公司 11262
• 代理人: 栗若木;曲鵬

技術領域

本發明涉及圖像重構技術，尤其涉及一種電鏡三維圖像重構的方法及系統。

背景技術

英特爾(Intel)公司推出的眾核(Many?Integrated?Core，MIC)處理器，跟通用的多核至強處理器相比，MIC眾核架構具有更小的內核和硬件線程，眾核處理器計算資源密度更高，片上通信開銷顯著降低，更多的晶體管和能量，能夠勝任更為復雜的并行應用。Intel?MIC產品基于X86架構，基于重核的眾核處理器，包含50個以上的核心，以及512bit的向量位寬，雙精性能超過1TFlops。

MIC擁有極其靈活的編程方式，MIC卡可以作為一個處理器存在，也可以被看作是一個獨立的節點。基本的MIC編程模型是將MIC看作一個處理器，中央處理單元(CPU)根據程序的指令，將一部分代碼運行在MIC端。此時存在兩類設備，即CPU端和MIC眾核處理器端。

電子斷層三維重構技術(Electron?Tomography，ET)是從一個物體的二維投影圖經過重構獲取物體內部結構的技術，通過獲取多個不同角度的二維投影圖進行反向重構出所研究對象的三維結構。圖像重構無論是在醫學診斷，還是在生物樣品的成像，以及在地表層析成像等諸多領域中都有著廣泛的應用。

目前，采用迭代重構圖像方法實現圖像重構，該迭代法可以表述成：Ax＝p，其中A為非奇異矩陣，x為未知數向量，p為已知向量，迭代法求解目標在于確定x的數值。在迭代重構圖像方法中，A為加權因子，p為投影得到的圖像像素值，x則為需要重構的圖像的像素值。

圖1說明了采用迭代重構圖像方法實現圖像重構的過程，該方法通過賦0值或隨機值或背投影方法獲取用于迭代過程中的初值X(0)，然后經過多次的重投影和背投影的迭代過程，來求得最后滿足精度要求的重構圖像的像素值。

聯合迭代重構(simultaneous?Iterative?Reconstruction?Technique，SIRT)算法在所有圖像重構算法中重構效果較好，但其運算量巨大，因此需要提高運算效率，來滿足快速圖像重構的需求。

發明內容

本發明提供了一種電鏡三維圖像重構的方法及系統，以克服現有圖像重構速度較慢的缺陷。

本發明提供了一種電鏡三維圖像重構的方法，該方法包括：

中央處理器(CPU)將要進行迭代重構圖像中的物體分成切片，構建所述CPU與眾核(MIC)處理器協同計算的框架；

所述CPU調用所述框架中的所有設備多線程并行地獲得各自當前切片的重構圖像值和測量投影圖像值，根據所述重構圖像值和所述測量投影圖像值依次通過重投影算法和背投影算法計算出當前切片的最終重構圖像值，直至計算出所有切片的最終重構圖像值。

優選地，所述CPU將要進行迭代重構圖像中的物體分成切片，包括：

所述CPU根據所述物體的大小將所述物體分成切片，切片的數目為該物體寬度的像素值。

優選地，所述CPU構建所述CPU與眾核(MIC)處理器協同計算的框架，包括：

所述CPU將連接在同一單節點服務器的所述CPU與多個MIC處理器作為協調計算的框架，該框架中所有設備的數目為CPU的數目和MIC處理器的數目之和。

優選地，所述方法還包括：所述所有設備分別在獲得各自當前切片的所述重構圖像值和所述測量投影圖像值之后，執行上鎖操作。

優選地，所述CPU調用所述框架中的所有設備多線程并行地獲得各自當前切片的重構圖像值和測量投影圖像值，根據所述重構圖像值和所述測量投影圖像值依次通過重投影算法和背投影算法計算出當前切片的最終重構圖像值，包括：

所述MIC處理器或所述CPU根據獲得的所述當前切片的初始重構圖像值X(0)或上次迭代得到的重構圖像值X(k-1)進行重投影得到計算投影圖像值；根據所述X(0)或X(k-1)得到所述測量投影圖像值，根據所述計算投影圖像值和所述測量投影圖像值進行背投影得到重構圖像值X(k)，通過多次的重投影和背投影操作得到滿足精度要求的所述當前切片的重構圖像值，其中，所述初始重構圖像值X(0)通過賦初值或隨機值或背投影方法獲得，k為大于1的整數。

本發明還提供了一種電鏡三維圖像重構的系統，該系統包括：