用於截斷補償的方法和設備的製作方法
2024-02-09 14:11:15 1
專利名稱：用於截斷補償的方法和設備的製作方法
技術領域：
本發明一般涉及計算機斷層掃描(CT)圖像重建中的方法和設備，尤其涉及用於截斷補償的方法和設備。
背景技術：
在某些掃描條件下，病人部分身體可能會超出檢測器所能測量的區域之外，這會導致圖像偽影和物體成像不完整。已經公布了一些已知方法，其解決減少偽影，而不是超出視野(FOV)的病人部分身體成像。然而，有必要對超出FOV的病人部分身體成像，例如在2003年6月24日提交的、具有代理人案卷號為129993的、且要求2002年10月4日所提交的、臨時申請序號為60/416,072的申請的利益的、題為「用於截斷補償的方法和設備」的共同待審申請序號##/###,###中所描述的。另外，有必要改進超出視野的病人部分身體的圖像質量。
發明內容
一方面，提供了一種確定投影是否截斷的方法。該方法包括計算掃描物體的每個投影視圖中的所有採樣的和；確定所計算的和的最大值；當索引k的所有採樣的和小於最大值的預先確定的百分比時，求投影視圖索引k的多個採樣m的均值；將均值和門限t進行比較，當均值大於t時，確定投影為截斷，當均值不大於t時，確定投影為未截斷。
另一方面，提供了一種成像設備。這種成像設備包括射線源；檢測器，其對射線響應，且置於接收該源所放射射線的位置上；以及計算機，其操作地耦合到該射線源和該檢測器。配置該計算機以計算掃描物體的每個投影視圖中的所有採樣的和；確定所計算的和的最大值；當投影視圖索引k的所有採樣的和小於該最大值的預先確定的百分比時，求索引k中多個採樣m的均值；將均值和門限t進行比較，當均值大於t時，確定投影為截斷，當均值不大於t時，確定投影為未截斷。
再一個方面，提供了一種計算機可讀介質，包含配置的編碼程序，以指示計算機計算掃描物體的每個投影視圖中的所有採樣的和；確定所計算的和的最大值；當索引k的所有採樣的和小於所述大值的預先確定的百分比時，求投影視圖索引k的多個採樣m的均值；將均值和門限t進行比較，當均值大於t時，確定投影為截斷，當均值不大於t時，確定投影為未截斷；並按(k)=k2-kk2-k1(k1)+k-k1k2-k1(k2)]]>估計總衰減τ(k)，其中k1和k2是與含至少一個確定為截斷投影的截斷區域相鄰的未截斷的視圖的多個視圖位置的均值；按(k)=i=1Np(i,k)]]>計算ξ(k)。進一步配置該程序以指示計算機按λ(k)＝τ(k)-ξ(k)確定衰減差額λ(k)；按(k)=2Rl2(k)-xl(k)pl(k)-Rl2(k)arcsin(xl(k)Rl(k))+2Rr2(k)-xr(k)pr(k)-Rr2(k)arcsin(xr(k)Rr(k))]]>計算要增加的衰減總量η(k)，其中pl(k)和pr(k)分別是左右投影邊界上的數個視圖的平均採樣值幅值，xl(k)、xr(k)、Rl(k)和Rr(k)分別是左右擬合柱面的位置和半徑；通過計算比率(k)=(k)w(k),]]>將η(k)和λ(k)進行比較，其中μw是水的衰減係數；將ε(k)和門限q進行比較；當ε(k)不大於q時，至少使用η(k)和λ(k)中的一個校正圖像。進一步配置該程序以指示計算機當ε(k)大於q時，基於和k1-n視圖以及k2+n視圖有關的數據，計算ηn(k)，其中n是整數；以及使用ηn(k)校正圖像。
圖1是CT成像系統實施例圖；圖2是圖1所示系統的方框示意圖；圖3是僅出現在部分投影角度上的投影截斷圖；圖4是不進行截斷校正時的重建圖像；圖5是根據代理人案卷號為129993描述的方法作截斷校正得到的圖像；圖6是按照投影視圖的函數(function of projection view)繪製的所測量的每個投影的衰減總量圖；圖7是擬合(fitting)過程的圖示；圖8說明多個圖像，且包括圖像，其僅用單個視圖來估計左截斷投影的幅值和斜率pl，sl，且也用單個視圖來估計右截斷投影的幅值和斜率pr，sr；圖9說明如何基於所估計的柱面計算缺少的投影採樣；圖10是用這裡所述方法和設備對圖4和圖5所示的同樣的體模(phantom)掃描進行校正得到的圖像。
具體實施例方式
於此提供的截斷補償方法和設備用於擴展CT系統的視野。參考附圖，對該設備和方法進行詳細的描述，所有附圖中的相同的附圖標記都指相同的元件。附圖的用意在於例證而不用作限制，且附於此的目的是有助於對本發明中示例性方法和設備的實施例進行解釋。
在一些已知的CT成像系統配置中，射線源投射的扇形光束被校準在位於笛卡兒坐標系的X-Y平面內，該平面通常稱為成像平面。射線光束穿過被成像的物體，如病人。經物體衰減之後，光束照射射線檢測器陣列。檢測器陣列接收到的衰減後的射線光束的強度取決於物體對射線光束的衰減程度。陣列中的每個檢測器元件產生單獨的電信號，它是該檢測器處的光束衰減測量值。從所有檢測器獨立地獲得衰減測量值以產生透射輪廓(transmissionprofile)。
在第三代CT系統中，射線源和檢測器陣列隨臺架在成像平面內圍繞被成像的物體旋轉，從而射線光束和物體交叉的角度以常量變化。檢測器陣列在一個臺架角度上獲得的一組射線衰減測量(即，投影數據)被稱為一個視圖。對物體的掃描包含射線源和檢測器旋轉一周期間在不同臺架角度或角度上獲得的一組視圖。
在軸向掃描中，處理投影數據以重建相應於穿過物體的二維層(twodimensional slice)圖像。從一組投影數據重建圖像的方法技術上稱之為過濾背景投影技術。該處理將掃描所獲得的衰減測量值轉換成稱為「CT數」或「Hounsfield單位」的整數，用以控制顯示設備上相應像素的亮度。
為減少總掃描時間，可作螺旋狀掃描。為進行螺旋狀掃描，在獲取規定量的層數據的時候，要移動病人。這種系統從一次扇形光束螺旋狀掃描中產生單個螺旋。由扇形光束映射的螺旋產生投影數據，根據該數據重建每個規定層中的圖像。
如這裡的用法，在所引述的單數形式元素或步驟之前的「a」或「an」應理解成不排除複數形式的所述元素或步驟，除非明確說明。另外，本發明所稱的一個實施例也不要解釋成排除包含所述特性的其它實施例的存在。
如這裡的用法，短語「重建圖像」也不要解釋成排除本發明中只產生表示圖像的數據而不產生可視圖像的實施例。因此，如這裡的用法，術語「圖像」泛指可視圖像和表示可視圖像的數據。然而，許多實施例產生(或配置來產生)至少一種可視圖像。
圖1是CT成像系統(10)的圖。圖2是圖1所示系統(10)的方框示意圖。在這個示例性實施例中，CT成像系統10包括臺架12，表示第三代CT成像系統。臺架12有射線源14，將錐形X射線光束16投射到臺架12另一端的檢測器陣列18的方向上。
檢測器組18由多個檢測器行(圖中未顯示)組成，包括多個檢測器元件20，都用於感知穿過物體(如病人22)的投射出的X射線光束。每個檢測器元件20產生電信號，以表示照射的射線光束的強度，也即光束穿過物體或病人22時的衰減。成像系統10有多層檢測器18，能提供多個表示物體22的圖像。這些多個圖像中的每一幅圖像相應於其中的單獨層。層的厚度或縫隙取決於檢測器行的厚度。
在掃描獲取射線投影數據期間，臺架12和其上安裝的組件圍繞旋轉中心24旋轉。圖2僅顯示單行檢測器元件20(即檢測器行)。然而，多層檢測器陣列18包括由檢測器元件20構成的多個平行的檢測器行，因此相應於多個準平行或平行層的投影數據可以在一次掃描中同時獲得。
臺架12的旋轉和射線源14的操作由CT系統10中的控制機械裝置26管理。控制機械裝置26包括射線控制器28，其為射線源14提供電源和定時信號；以及臺架電動機控制器30，其控制臺架12的旋轉速度和位置。控制機械裝置26中的數據採集系統DAS 32從檢測器元件20採樣模擬數據，並把這些數據轉換成數位訊號，用於後續處理。圖像重建器34從DAS 32接收採樣的數位化的射線數據，並進行高速圖像重建。重建圖像作為計算機36的輸入，將圖像存儲於海量存儲設備38中。
計算機36也通過具有鍵盤的控制臺40接收來自操作員的命令和掃描參數。相關的陰極射線管顯示器42使操作員能觀察來自計算機36的重建圖像和其它數據。計算機36使用操作員提供的命令和參數為DAS 32、射線控制器28和臺架電動機控制器30提供控制信號和信息。此外，計算機36操縱工作檯電動機控制器44，以控制機動化的工作檯46，把病人22放到臺架12中。尤其是，工作檯46將病人22的部分身體移動通過臺架開口48。
在一個實施例中，計算機36包括設備50，如磁碟驅動器或CD-ROM驅動器，以讀取計算機可讀介質52(如磁碟或CD-ROM)中的指令或數據。在另一個實施例中，計算機36執行存儲在固件(未示出)中的指令。通常，對圖2的DAS 32、重建器34和計算機36中至少一個的所含的處理器進行編程，以執行以下描述的程序。當然，這種方法不僅限於在CT系統10中使用，並可用於其它類型或種類的成像系統中。在一個實施例中，對計算機36編程，以執行上述功能，相應地，如這裡的用法，術語「計算機」不限於技術上所指的集成電路意義上的計算機，而泛指計算機、處理器、微控制器、微機、可編程控制器、ASIC和其它可編程電路。儘管這裡描述的方法針對醫療場合，可預期本發明的益處對非醫用成像系統，如工業場合或運輸場合中常使用的系統，如但不局限於飛機場或其它運輸中心的行李掃描CT系統，也是十分關鍵的。
已經公布了一些已知方法，其解決減少偽影，而不是超出視野的病人部分身體成像。然而，有必要對超出視野的病人部分身體成像。這在許多領域內如腫瘤學、旋轉血管造影術、融合成像系統以及商用CT掃描器十分有用。至少有一種已知的多層CT掃描器，其現有硬體將重建視野限制在大約50cm範圍內。儘管這對大多數醫療應用是足夠的，也有必要對視野進行擴展，以對該視野之外的物體成像。這對諸如腫瘤學或CT/PET等應用非常有益。對腫瘤學應用來說，需要一個較大的視野。這主要是因為對射線診療方案而言，病人的四肢常常放在掃描視野之外，以更好的定位瘤的位置。已知的CT重建算法不考慮截斷投影，且產生的圖像有嚴重的偽影。從治療方案講，這些偽影可能影響精確估計衰減路徑。這裡所描述的是一種算法途徑，以將重建視野(FOV)擴大到由檢測器硬體所限制的FOV之外。校正算法可應用於不同的重建算法，包括但不受限於全掃描、半掃描/分割、螺旋狀掃描、以及基於心臟部分(cardiac sector)的算法。另外，系統10的配置採用了這裡描述的算法。
圖3是出現在部分投影角度上的投影截斷圖。在某些掃描條件下，病人22的部分身體超出檢測器18所能測量的區域，這會導致圖像偽影和成像物體不完整。X射線管和檢測器18被牢牢固定在繞病人22旋轉的框架上。在旋轉過程中，位於完全採樣視野60內，會獲得不間斷的測量。穿過物體22的完全採樣視野60之外的任何區域的射線衰減只能在有限的旋轉角範圍內測量到，且該區域被稱為部分採樣視野區域。換句話說，位於完全採樣視野60內的部分被定位在扇形16之內，因此從所有臺架角度都能測量到，且所測得的數據被定義成完全採樣視野數據。然而，一些部分在某些角度上位於扇形16之內，而在另一些角度上卻位於扇形16之外，與這些部分有關的採集數據被定義成部分採樣視野數據。在這個例子中，時針3所指位置上的投影不存在截斷，但時針12所指位置上的投影嚴重截斷。因此，可依賴未截斷的投影(如圖3中大約時針3所指位置上)估計截斷視圖的截斷總量(如圖3中大約時針12所指位置上)。校正處理初期步驟是在處理投影之前用軟體完成扇形光束到平行光束的轉換。在一個實施例中，該初期步驟是第一步。該處理技術眾所周知，且不需要專門的數據收集。當轉換完成之後，將所有檢測器通道上的投影集中起來，以獲得下面所詳細描述的總衰減曲線。
圖4是不進行截斷校正時的重建圖像80。容易觀察到，圖像80包含截斷偽影82。
代理人案卷號為129993描述了一種針對大多數情況的截斷補償算法。然而，在某些情況下，代理人案卷號為129993中描述的補償算法不能獲得最優結果。圖5是根據代理人案卷號為129993中的方法作截斷校正後的圖像。儘管圖像84比圖像80有很大改進，但校正區域86(恢復手臂)內仍存在失真。在該例子中，胸模置於最外端區域距離相同中心29cm的位置上(在58cm FOV邊緣)。儘管上述校正消除了50cm視野以內的偽影，但所恢復手臂86的形狀和強度從某種程度上而言存在失真。
相信失真主要由所估計的缺少衰減不精確引起。為了說明起見，圖6是根據投影視圖的函數繪製的所測得的每個投影的總衰減90的圖。為清楚起見，對該圖進行了放大。在衰減測量值90中，下降區92相應於截斷視圖。從該圖中可清楚看到，對多個未截斷投影(衰減測量值90中的近似平坦區域)而言，為每個投影所測得的投影測量值的和在所有投影角度上不是常數。這是由非理想化校準和扇形到平行轉換處理(parallel rebinning process)引起的。儘管變化較小(小於2％)，它也能表示整個截斷投影採樣中的主要部分。換句話說，和所估計的缺少投影相比，該變化百分比要大得多。因此，把其他的信息用於截斷估計，而不是用於代理人案卷號為129993中所描述的補償，將會帶來更少的、甚至沒有失真。
在多個體模試驗期間，發現僅投影邊界上的幅值和斜率有時不足以估計擬合投影的柱面物體的大小和位置，因為有時邊界採樣的斜率偏離所掃描物體的內部結構。例如，如果一個氣穴碰巧位於截斷投影邊界採樣的切線上，在外插投影採樣中會出現明顯的不連續(在投影角度上)。這會引起陰影和拖尾偽影。
當投影採樣有噪聲時(由於有限的x射線光子統計)，發生另一種複雜情況。在這種情況下，投影邊界的幅值和斜率都變得不可靠。高頻變化(從一個視圖到另一個視圖)可導致嚴重的拖尾以及物體形狀高度失真。臨床應用上，這常發生在掃描體大的病人時，或當使用低流量(low-flux)獲取協議時。這裡描述的方法和設備至少部分解決上述所描述的難題。
更明確地講，投影被聲明為截斷或未截斷。要聲明為截斷投影，必須滿足一個條件，即投影邊界採樣的任一端(在幾個通道上的均值)必須小於預先確定的門限，除先前代理人案卷號為129993描述的條件之外。多個採樣的均值是指多個邊緣採樣的均值。即如果投影包含N個採樣，索引從0到N-1，使用從0到m-1左面m個投影採樣的均值以及從N-m到N-1右面m個採樣的均值。除上述測試之外，仍使用代理人案卷號為129993中的測試。即計算所有視圖的所有N個採樣的和。計算所有投影視圖中的最大值，且當一個特定視圖的和小於預先確定的最大值的百分比時，該視圖被聲明為截斷的。門限值非零以考慮非理想校準。從算式上講，如果下列條件存在，投影被聲明為截斷pl=1mi=1mp(i,k)>t]]>或pr=1mi=1mp(N-i,k)>t---(1)]]>其中N是檢測器通道的數量，k是投影視圖索引，m是採樣數，t是門限值。在示例性的實施例中，t約等於0.42時，從經驗上講，性能最好。在其它實施例中，t約在0.25和0.58之間，t約在0.33和0.5之間，t約在0.375和0.46之間。
用上面描述的確定方法，識別的未截斷視圖在圖6中以雙箭頭顯示。為了估計缺少投影的量，相鄰未截斷投影的總衰減用於估計。尤其，截斷視圖的總衰減通過對位於截斷視圖區域兩端的兩個未截斷投影線性內插獲得。對圖6中所示的例子來說，第一組截斷視圖位於投影視圖號290(投射角106°)和445(投射角163°)之間。從視圖號289和446獲得的總衰減用於估計截斷視圖的總衰減。通過對兩個未截斷視圖抽樣線性內插進行估計，如圖6中的點線所示。從算式上而言，所估計的總衰減τ(k)是(k)=k2-kk2-k1(k1)+k-k1k2-k1(k2)---(2)]]>其中k1和k2是和截斷區域相鄰的未截斷視圖的視圖位置，ξ(k)是通過將視圖k的投影測量值求總和而獲得的總衰減。從算式上而言，它由下面的等式描述(k)=i=1Np(i,k)---(3)]]>為改進算法的穩健性，可以使用幾個未截斷投影的投影視圖的均值。例如，可以使用視圖285到289的平均總衰減來代替只根據視圖289獲得的總衰減，且k1為視圖285到289的平均值。在示例性的實施例中，取五個視圖的平均值，但在其它實施例中，取非五個視圖的平均值(如3-7個視圖或2-8個視圖)。另外，k2可指單個視圖或與k1無關的幾個視圖的平均值。例如，參照圖6中的例子，若k1指視圖285到289的平均值，則在一個實施例中，k2指視圖446到450的平均值，在另一個實施例中，k2指446。另外，k1和k2無關，且k1可指單個視圖或與k2無關的幾個視圖。
用於截斷補償的投影增加量λ(k)是所估計的總衰減和所測量的總衰減之間的差值λ(k)＝τ(k)-ξ(k)(4)下一步，估計可擬合每個截斷投影邊界的水柱面物體的大小和位置以推導待估計的投影。基於m個投影視圖的截斷投影邊界採樣的平均值的幅值和斜率計算擬合柱面的位置和半徑xl(k)=-sl(k)pl(k)4w2,]]>和Rl(k)=pl2(k)4w2+xl2(k)---(5)]]>xr(k)=-sr(k)pr(k)4w2,]]>和Rr(k)=pr2(k)4w2+xr2(k)---(6)]]>其中pl，sl和pr，sr分別是左右數個視圖的投影邊界的採樣平均值幅值和斜率，xl、xr、Rl和Rr分別是左右擬合柱面的位置和半徑，μw是水的衰減係數。擬合過程如圖7所示。
此前，表明在缺少x射線光子的條件下，參數估計過程會變得不可靠，因為所估計的參數從一個視圖到另一個視圖會有極大波動。圖8說明包含圖像102的多個圖像100，其只使用單個視圖估計左截斷投影的幅值和斜率pl，sl，也只使用單個視圖估計右截斷投影的幅值和斜率pr，sr。圖像102中清楚可見由參數波動引起的拖尾偽影。為克服這種可能的不足，在一個實施例中，使用幾個視圖的投影採樣的均值，以估計截斷投影的幅值和斜率pl，slpr，sr。在示例性的實施例中，使用7個視圖的投影採樣估計幅值和斜率。多個圖像100中的平均圖像104使用7個視圖估計幅值和斜率。在另一個實施例中，5到9個視圖用於估計幅值和斜率。之後，用平均值確定擬合柱面的位置和半徑。基於擬合柱面，對於投影k增加的柱面的總衰減η(k)(在最初投影之外的柱面部分)按下式計算(k)=2Rl2(k)-xl(k)pl(k)-Rl2(k)arcsin(xl(k)Rl(k))+2Rr2(k)-xr(k)pr(k)-Rr2(k)arcsin(xr(k)Rr(k))---(7)]]>在一個實施例中，將η(k)與所估計的增量衰減測量值λ(k)(如等式4描述)進行核對。在理想情況下，這兩個量應該相等。但在實際情況下，因為物體的缺少部分極少為一個理想柱面，所以這兩個量之間總是有一些差別。因此，計算這兩個量的比值ε(k)(k)=(k)w(k)---(8)]]>當ε(k)完全不等於1時，說明用擬合的柱面參數有問題。這可能是由被掃描物體的不規則形狀引起的，正如上面所描述的。因此，當ε(k)＞q時，不能使用所估計的參數作投影補償。q是一個預先確定的門限，且在示例性實施例中設為2。在另一個示例性實施例中，q在1.5和2.5之間。在另一個實施例中，q在1.75和2.25之間。在另一個示例性實施例中，q在1.9和2.1之間。
對ε(k)＞q的投影來說，可使用參數擬合較好的兩個投影(即ε(k)≤q)的柱面參數的內插。即如果視圖k1到k2的ε(k)值超過門限，可使用視圖k1-n和k2+n的擬合柱面參數來估計這些視圖的柱面參數，其中n是參數，n＝1表示第一個較好擬合的視圖。在示例性的實施例中，選擇n＝5，採用線性內插足夠好。即xl(k)、xr(k)、Rl(k)和Rr(k)是基於xl(k1-n)、xl(k2+n)、xr(k1-n)、xr(k2+n)、Rl(k1-n)、Rl(k2+n)、Rr(k1-n)和Rr(k2+n)的線性內插，且使用內插按等式7計算ηn(k)。比值ε(k)也從同樣的視圖組作線性內插。當然，也可使用其它內插方法，如高階插值。另外，沿視圖k的額外低通過濾可被應用於所估計的參數以改進參數的連續性。在其它實施例中，n在2和8之間，或n在3和7之間。
最後，基於所估計的柱面計算缺少的投影採樣。圖9圖示這個過程。因為這些柱面的總衰減不能匹配按一致性條件(等式4)估計的整個缺少衰減，要作進一步調整。例如，x軸(檢測器通道刻度i)按比率縮放，使初始估計投影(圖9中點線所示)在x上進一步擴展(如圖9中實粗線所示)或壓縮。這可通過改變檢測器通道i的刻度距離，或對初始計算出的投影進行內插來實現。
用這裡描述的方法和設備，代理人案卷號為129993所描述的方法和設備的性能和穩健性大大提高。圖10是用這裡所述方法和設備對圖4和圖5所示同樣的體模掃描進行校正得到的圖像。位於掃描視野之外的物體106形狀和密度的精確性明顯好轉，這是本發明的技術效果和貢獻。當然，其它技術效果也存在。
儘管上面所描述的系統和方法只使用了總衰減、邊界採樣的幅值和斜率來估計缺少投影的分布，其他信息也能用於估計。例如，對X線斷層攝影術而言，可使用Helgason-Ludwig條件(HL條件)進一步提煉上述技術。另外，設置不同的門限以確保算法在錯誤測量條件下正常工作。例如，可設置圖9中拉伸率的上下限，以防止由於不可靠的測量造成誤差增加。另外，為所計算的斜率sl，sr設置一個合理的範圍。如果被掃描物體的材料特徵和水完全不同，也可使用已知材料的衰減係數來計算等式3和4中的大小和位置。
因為內插數據不可能有和完全採樣視野內的數據一樣的圖像質量，所以對FOV外推的圖像加以標記就十分有用。在一個實施例中，提供了在重建圖像中描繪代表完全採樣視野數據和部分採樣視野數據區域的方法。也可以在CT號中用色彩代碼和變換的CT號來標記。因為標記可能影響圖像數據的可視性，提供了打開和關閉標記的簡易方式。因此，允許系統10的用戶打開和關閉該標記。
儘管本發明參照不同的特定實施例進行描述，本領域的技術人員理解本發明能夠在權利要求的實質和範圍內加以修改。
權利要求
1.一種成像設備(10)，包括射線源(14)；檢測器(18)，其對射線響應，並置於接收從所述源放射的射線的位置；以及計算機(36)，其操作地耦合到所述射線源和所述檢測器，且所述計算機配置為計算掃描物體(22)的每個投影視圖的所有採樣的和；確定所計算的和的最大值；當索引k的所有採樣的和小於所述最大值的預先確定的百分比時，求投影視圖索引k的多個採樣m的均值；將均值和門限t進行比較；當均值大於t時，確定投影為截斷；以及當均值不大於t時，確定投影為未截斷。
2.根據權利要求1所述的設備(10)，其中進一步配置所述計算機(36)以將均值和門限t進行比較，其中t大約在0.25和0.58之間。
3.根據權利要求1所述的設備(10)，其中進一步配置所述計算機(36)按下式估計總衰減(k)=k2-kk2-k1(k1)+k-k1k2-k1(k2)]]>其中k1和k2是與含至少一個確定為截斷的投影的截斷區域相鄰的未截斷視圖的視圖位置，且按(k)=i=1Np(i,k)]]>計算ξ(k)。
4.根據權利要求1所述的設備(10)，其中進一步配置所述計算機(36)以按下式估計總衰減(k)=k2-kk2-k1(k1)+k-k1k2-k1(k2)]]>其中k1和k2是與含至少一個確定為截斷的投影的截斷區域相鄰的未截斷的視圖的多個視圖位置的均值，且按(k)=i=1Np(i,k)]]>計算ξ(k)。
5.根據權利要求4所述的設備(10)，其中進一步配置所述計算機(36)為按λ(k)＝τ(k)-ξ(k)確定衰減差額λ(k)；按(k)=2Rl2(k)-xl(k)pl(k)-Rl2(k)arcsin(xl(k)Rl(k))+2Rr2(k)-xr(k)pr(k)-Rr2(k)arcsin(xr(k)Rr(k))]]>計算要增加的衰減總量η(k)，其中pl(k)和pr(k)分別是左右投影邊界上的數個視圖的平均採樣值幅值，xl(k)、xr(k)、Rl(k)和Rr(k)分別是左右擬合柱面的位置和半徑；通過計算比率(k)=(k)w(k),]]>將η(k)和λ(k)進行比較，其中μw是水的衰減係數；將ε(k)和門限q進行比較；當ε(k)不大於q時，至少使用η(k)和λ(k)中的一個校正圖像；以及當ε(k)大於q時基於和k1-n視圖以及k2+n視圖有關的數據，計算ηn(k)，其中n是整數；以及使用ηn(k)校正圖像。
全文摘要
一種方法包括計算掃描物體的每個投影視圖中的所有採樣的和；確定所計算的和的最大值；當索引k的所有採樣的和小於最大值的預先確定的百分比時，求投影視圖索引k的多個採樣m的均值；將均值和門限t進行比較；當均值大於t時，確定投影為截斷，當均值不大於t時，確定投影為未截斷。
文檔編號A61B6/03GK1575760SQ200410061719
公開日2005年2月9日 申請日期2004年6月30日 優先權日2003年7月3日
發明者謝強, 布賴恩·格雷科維茨, 愛德華·H·趙, 斯科特·M·麥科拉什, 埃米·L·霍斯特 申請人:Ge醫藥系統環球科技公司