TWI816171B - 用於基於輻射的成像的擴散場成像準直儀 - Google Patents
用於基於輻射的成像的擴散場成像準直儀 Download PDFInfo
- Publication number
- TWI816171B TWI816171B TW110130688A TW110130688A TWI816171B TW I816171 B TWI816171 B TW I816171B TW 110130688 A TW110130688 A TW 110130688A TW 110130688 A TW110130688 A TW 110130688A TW I816171 B TWI816171 B TW I816171B
- Authority
- TW
- Taiwan
- Prior art keywords
- collimator
- aperture
- detector
- radiation
- apertures
- Prior art date
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 88
- 230000005855 radiation Effects 0.000 title claims abstract description 77
- 125000006850 spacer group Chemical group 0.000 description 30
- 238000000034 method Methods 0.000 description 21
- 239000011148 porous material Substances 0.000 description 18
- 238000013461 design Methods 0.000 description 14
- 239000011295 pitch Substances 0.000 description 13
- 229940121896 radiopharmaceutical Drugs 0.000 description 11
- 239000012217 radiopharmaceutical Substances 0.000 description 11
- 230000002799 radiopharmaceutical effect Effects 0.000 description 11
- 230000035945 sensitivity Effects 0.000 description 11
- 238000000926 separation method Methods 0.000 description 10
- 238000009826 distribution Methods 0.000 description 9
- 238000005286 illumination Methods 0.000 description 9
- 238000009206 nuclear medicine Methods 0.000 description 6
- 238000013459 approach Methods 0.000 description 5
- 239000000463 material Substances 0.000 description 5
- 239000011159 matrix material Substances 0.000 description 5
- QWUZMTJBRUASOW-UHFFFAOYSA-N cadmium tellanylidenezinc Chemical compound [Zn].[Cd].[Te] QWUZMTJBRUASOW-UHFFFAOYSA-N 0.000 description 4
- 230000008878 coupling Effects 0.000 description 4
- 238000010168 coupling process Methods 0.000 description 4
- 238000005859 coupling reaction Methods 0.000 description 4
- 229910001385 heavy metal Inorganic materials 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000001172 regenerating effect Effects 0.000 description 4
- 229910052721 tungsten Inorganic materials 0.000 description 4
- 229910052720 vanadium Inorganic materials 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 230000000903 blocking effect Effects 0.000 description 3
- 230000007423 decrease Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000033001 locomotion Effects 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000013507 mapping Methods 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- MARUHZGHZWCEQU-UHFFFAOYSA-N 5-phenyl-2h-tetrazole Chemical compound C1=CC=CC=C1C1=NNN=N1 MARUHZGHZWCEQU-UHFFFAOYSA-N 0.000 description 2
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 2
- 239000011358 absorbing material Substances 0.000 description 2
- XQPRBTXUXXVTKB-UHFFFAOYSA-M caesium iodide Chemical compound [I-].[Cs+] XQPRBTXUXXVTKB-UHFFFAOYSA-M 0.000 description 2
- 230000015556 catabolic process Effects 0.000 description 2
- 238000006731 degradation reaction Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 239000007769 metal material Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 229910052710 silicon Inorganic materials 0.000 description 2
- 239000010703 silicon Substances 0.000 description 2
- 238000002603 single-photon emission computed tomography Methods 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- WFKWXMTUELFFGS-UHFFFAOYSA-N tungsten Chemical compound [W] WFKWXMTUELFFGS-UHFFFAOYSA-N 0.000 description 2
- 239000010937 tungsten Substances 0.000 description 2
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 206010036618 Premenstrual syndrome Diseases 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 239000000956 alloy Substances 0.000 description 1
- 229910045601 alloy Inorganic materials 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000005251 gamma ray Effects 0.000 description 1
- 229910052732 germanium Inorganic materials 0.000 description 1
- GNPVGFCGXDBREM-UHFFFAOYSA-N germanium atom Chemical compound [Ge] GNPVGFCGXDBREM-UHFFFAOYSA-N 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000009828 non-uniform distribution Methods 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000035790 physiological processes and functions Effects 0.000 description 1
- 238000002600 positron emission tomography Methods 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 230000002285 radioactive effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- FVAUCKIRQBBSSJ-UHFFFAOYSA-M sodium iodide Chemical compound [Na+].[I-] FVAUCKIRQBBSSJ-UHFFFAOYSA-M 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/06—Diaphragms
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/037—Emission tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4266—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a plurality of detector units
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4291—Arrangements for detecting radiation specially adapted for radiation diagnosis the detector being combined with a grid or grating
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Engineering & Computer Science (AREA)
- Radiology & Medical Imaging (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Measurement Of Radiation (AREA)
- Nuclear Medicine (AREA)
- Radiation-Therapy Devices (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
基於輻射的成像系統包括準直儀,該準直儀被配置為過濾從目標物體發射的輻射,該準直儀包括多個不均勻分佈在該準直儀上的孔徑。該多個孔徑的最大接收角不大於15°。基於輻射的成像系統還包括用於檢測已經通過準直儀的輻射的檢測器。準直儀與檢測器間隔開,以使得面向準直儀的檢測器的頂表面上的點同時被該多個孔徑中的兩個或更多個照明。
Description
在基於輻射的成像中,例如分子醫學成像(有時稱為核子醫學成像),可以產生表示放射性藥物分佈的影像以用於醫學診斷。在成像之前,放射性藥劑被注入到目標物體中,例如患者。放射性藥劑發射放射性光子,可以穿透身體被光子檢測器檢測到。基於來自接收到的光子的資訊,光子檢測器隨後可以確定患者體內放射性藥劑的分佈。分佈代表了患者的生理功能,因此分佈的影像為心臟病學、腫瘤學、神經病學等的多種疾病和病症的診斷提供了有價值的臨床資訊。
為了產生影像,準直儀和光子檢測器協同工作。準直儀是引導光子路徑的裝置(即引導光子走一定的路徑)。在基於輻射的成像中,光子可能來自受試者內部的未知位置,這與X射線或CT中的光子從已知來源(或多
個來源)位置發射不同。如果沒有準直儀,來自各個方向的光子可能會被光子檢測器記錄,影像重建可能會變得困難。因此,準直儀被用來引導可能的光子路徑,從而可以重建影像,類似於鏡頭在攝影相機中的作用。儘管現有的基於輻射的成像系統通常足以滿足其預期目的,但它們並非在所有方面都完全令人滿意。例如,現有的準直儀和檢測器通常必須在成像解析度和訊號靈敏度之間進行權衡部署,但不能同時兼顧兩者。因此,需要改進基於輻射的成像系統。
根據各種實施例,本發明提供一種基於輻射的成像系統。該基於輻射的成像系統包括準直儀,該準直儀被配置為過濾從目標物體發出的輻射,該準直儀包括多個不均勻分佈在該準直儀上的孔徑,其中該多個孔徑的最大接收角不大於15°;以及檢測器,用於檢測已經通過準直儀的輻射,其中準直儀與檢測器間隔開,以使得面向準直儀的檢測器的頂表面上的點同時被該多個孔徑中的兩個或更多個孔徑照明。在一些實施例中,檢測器的頂表面上的點同時被多個孔徑的一定百分比照明,其中該百分比約小於25%。在一些實施例中,多個孔徑的數量超過一千。在一些實施例中,多個孔徑至少在孔徑尺寸、孔徑形狀、接收角、長度和孔徑間距中之一方面變化。在一些實施例中,準直儀與檢測器之間的距離為準直儀的厚度約0.5至
約10倍(例如1.5至10倍)。在一些實施例中,多個孔徑形成包括重複的基本圖案的圖案。在一些實施例中,重複的基本圖案是編碼孔徑圖案。在一些實施例中,其中該些孔徑中的一孔徑的照明區域實質上等於重複的基本圖案的區域。在一些實施例中,準直儀包含第一部分和第二部分,其中第一部分包含均勻分佈在第一部分上的通孔,第二部分包含不均勻分佈在第二部分上並且對應於多個孔徑的通孔,以使得第一部分的通孔的一部分被第二部分擋住。在一些實施例中,第一部分的厚度為第二部分的厚度約2至約10倍。在一些實施例中,準直儀是第一準直儀並且檢測器是第一檢測器,該系統還包括第二準直儀和耦合到第二準直儀的第二檢測器,其中第二準直儀附接到第二檢測器。
根據各種實施例,本發明提供一種基於輻射的成像系統。基於輻射的成像系統包括第一準直儀,其被配置為過濾從目標物體發射的輻射,第一準直儀包括形成第一孔徑圖案的第一多個孔徑;與第一準直儀相關聯的第一檢測器,用於檢測已經通過第一準直儀的輻射;第二準直儀被配置為過濾從目標物體發射的輻射,第二準直儀包括形成第二孔徑圖案的第二多個孔徑;以及與第二準直儀相關聯的第二檢測器,用於檢測已經通過第二準直儀的輻射,其中目標物體位於第一準直儀和第二準直儀之間,並且其中第一孔徑圖案與第二孔徑圖案不同。在一些實施例中,第一準直儀與第一檢測器接觸,並且其中第二準直儀
與第二檢測器間隔開。在一些實施例中,第二準直儀與第二檢測器之間的距離為第二準直儀的厚度約0.5至約7倍(例如1.5至7倍)。在一些實施例中,第一多個孔徑均勻地分佈在第一準直儀上,並且其中第二多個孔徑不均勻地分佈在第二準直儀上。在一些實施例中,第二孔徑圖案包括重複的編碼圖案。在一些實施例中,重複的編碼圖案是均勻冗餘陣列(uniformly redundant array,URA)圖案、修改的均勻冗餘陣列(modified uniformly redundant array,MURA)圖案及完美二進制陣列(perfect binary array,PBA)圖案其中之一。在一些實施例中,第一和第二檢測器指向的方向具有偏移角。
根據各種實施例,本發明提供一種用於過濾物體發射的輻射的準直儀。準直儀包括第一部分,第一部分包括均勻分佈的第一多個通孔;以及包括非均勻分佈的第二多個通孔的第二部分,其中第一部分和第二部分以如下方式間隔開並對齊,亦即,物體發射的光子穿過第一多個通孔中的一些通孔並被第二部分擋住。在一些實施例中,第二多個通孔形成編碼圖案。在一些實施例中,第二部分係可操作以相對於第一部分的面向第二部分的表面滑動,從而可以改變第一和第二多個通孔中對齊的通孔。
根據各種實施例,本發明提供了一種成像系統中未對準校正的方法。該方法包括提供準直儀和檢測器,檢測器具有數位化像素網格;用泛光光源照明準直儀;記錄數位化像素網格的每個像素中的訊號強度,其
中,訊號強度是由準直儀的照明引起的;從記錄的訊號強度中導出描述符;產生一系列被引入至數位化像素網格的偏移量;從最佳化描述符的一系列的偏移量中找到一個選定的偏移量;並且基於從所選定的偏移量導出的實際最佳偏移量來重新產生數位化像素網格。在一些實施例中,所選定的偏移量包括分別在兩個正交方向上的一對偏移值。在一些實施例中,描述符是該記錄的訊號強度的對比度值、峰值和谷值其中之一。在一些實施例中,對比度值是由對應於在數位化像素網格處記錄的訊號強度的訊號強度線的峰值和谷值的比值確定。在一些實施例中,選定偏移量的發現包括掃描偏移量系列並挑選最大化或最小化描述符的選定偏移量。在一些實施例中,重新產生數位化像素網格包括將選定的偏移量調整固定偏移量以產生實際最佳偏移量並用實際最佳偏移量重新產生數位化像素網格。
100:成像系統
102:整合高架
104:轉子
108,204:檢測器
110,202:準直儀
112:患者檢查台
114:機台支座系統
116:耦合和屏蔽組件
120:控制台
150:患者
152:圓
154A,154B:法線方向
203:SFI準直儀
206:目標物體
208,208':孔徑
210:準直儀頂表面
212:隔片
214,214a:感測像素
214b:相鄰感測像素
230:入射光子
240:單元網格
242:泛光源
244:屏蔽
262:平行孔區段
262-a和262-b:子平行孔區段
264:薄板
264-a/264-b:圖案化孔徑區段
A1,A2,A1',A2',C1,C2,D1,D2:線
B1,B2:虛線
D:間隔
H:高度/厚度
H1,H2,H1-a,H1-b,H2-a,H2-b:厚度
H3,Z:距離
P:孔徑間距
P':感測像素間距
R0,R1:射線
Rn,Rc:成像解析度
S:孔徑間隔
Vp和Vv:訊號強度
W:寬度
W':較小的直徑
Δx,Δy:未對準
245:訊號強度線
當結合附圖閱讀本案時,從以下詳細描述可以最好地理解本發明。需要強調的是,根據行業標準慣例,各種特徵並非按比例繪製,僅用於說明目的。事實上,為了討論的清楚起見,可以任意增加或減少各種特徵的尺寸。
[圖1A、1B和1C]是根據本發明的各個方面的示例性基於輻射的成像系統的示意圖。
[圖2A、2B、2C、2D、2E和2F]提供了根據
本發明的一些實施例的準直儀的一部分的透視圖、截面圖和俯視圖。
[圖3]是根據本發明的一些實施例的基於輻射的成像系統的一部分的截面圖。
[圖4]是根據本發明的一些其他實施例的基於輻射的成像系統的一部分的截面圖。
[圖5A和5B]示例了根據本發明的各個方面的具有非均勻分佈的孔徑的準直儀的一部分的俯視圖。
[圖6A、6B和6C]示例了根據本發明的一些實施例的具有重複的基本圖案的準直儀的俯視圖。
[圖7]示例了根據本發明的各個方面的在快速影像重建演算法中將原始影像分解為一系列子影像的步驟。
[圖8A、8B、8C、9A、9B、9C、10A和10B]示例了根據本發明的一些實施例的輕量級準直儀的一些實施例的截面圖。
[圖11A、11B、12A、12B和13]示例了具有對準調整的基於輻射的成像系統的一些實施例的俯視圖和截面圖。
以下揭露內容提供了許多不同的實施例或範例,用於實現本發明的不同特徵。下面描述的組件和配置的具體範例以簡化本發明。當然,這些僅僅是範例而不是
限制性的。對所描述的裝置、系統、方法以及本發明的原理的任何進一步的應用的任何改變和進一步修改都被完全考慮,如本發明所涉及的領域的普通技術人士通常會想到的。例如,關於一個實施例所描述的特徵、組件及/或步驟可以與關於本發明的其他實施例所描述的特徵、組件及/或步驟組合以形成根據本發明的裝置、系統或方法的又一個實施例,儘管這樣的組合沒有明確示出。此外,本發明可以在各種範例中重複參考數字及/或字母。這種重複是為了簡單和清楚起見,其本身並不規定所討論的各種實施例及/或配置之間的關係。
此外,在隨後的本發明中的一個特徵上的、連接至及/或耦合至另一個特徵的特徵可以包括該些特徵直接接觸的實施例,並且還可包括附加特徵可插入該些特徵的實施例,這樣該些特徵間可能不會直接接觸。此外,空間相對術語,例如“較低”、“較高”、“水平”、“垂直”、“上面”、“在...的上方”、“下方”、“在...的下方”、“向上”、“向下”、“頂部”、“底部”等,以及其衍生詞(例如,“水平地”、“向下地”、“向上地”等)用於簡化本發明的一個特徵與另一個特徵間的關係。空間相對術語旨在涵蓋包括該些特徵的裝置的不同定向。更進一步地,當用“大約”、“近似”等術語來描述數字或數字之範圍時,該術語旨在涵蓋在包括所描述的數字在內的合理範圍內的數字,例如在+/-10%內之本領域技術人士所理解的所描述的數量或其他值。例如,術語“約5公分”包含從4.5公分到5.5公分的尺寸
範圍。
本發明總體上涉及用於在基於輻射的成像中使用的準直儀,並且更具體地涉及用於在核子醫學(分子)成像系統中使用的具有非均勻分佈的孔徑的準直儀。術語“非均勻分佈的孔徑”是指準直儀的孔徑(通孔)在孔徑輪廓中的至少一個方面變化,包括但不限於孔徑大小、孔徑形狀(包括橫截面和縱向形狀)、孔徑長度、孔徑間距和孔徑方向。
在基於輻射的成像中,例如核子醫學(分子)成像系統,準直儀和檢測器協同工作以產生表示受試者內放射性藥物分佈的影像。許多核子醫學成像系統,例如單光子射出電腦斷層造影(single photon emission computed tomography,SPECT)和正子斷層造影術(positron emission tomography,PET)成像系統,使用一個或多個檢測器來獲取成像資料,例如伽馬射線或光子成像資料。在獲取影像之前,通常將放射性藥物口服或注射到諸如患者之類的物體中。放射性藥物經歷核衰變,透過直接或間接地互毀,以特定速率和特徵能量放射伽馬光子。一個或多個檢測器單元放置在物體周圍以記錄或監測放射。在很多情況下,為了製造和資料處理的方便,檢測器被組織成平面形狀,因此以二維矩陣格式獲取資料,這通常被稱為投影。基於包括位置、能量和此類檢測事件之計數在內的記錄資訊,可以重建放射性藥物分佈的影像以研究物體的某些部分(例如,患者的身體部分)的功能。
然而,現有的準直儀和檢測器的設計存在各種問題。例如,平行孔準直儀通常與檢測器緊密耦合(或附接)以減少檢測器中感測像素之間的串擾,使得只有一個可能的光子路徑到達任何感測像素。通常,準直儀孔徑越大,成像系統的成像解析度越有利,但會降低其訊號靈敏度。因此,必須在成像解析度和訊號靈敏度之間進行權衡來確定準直儀的孔徑長度。
本發明提供了準直儀設計的實施例,其中具有非均勻分佈的孔徑的準直儀被部署在遠離檢測器的距離處,使得檢測器的一些區域同時被穿過一個以上的孔徑的光子照明。換言之,準直儀的至少一些孔徑在檢測器上的各自的照明區域彼此重疊。透過將準直儀與檢測器隔開並利用更大的孔徑尺寸(或更寬的接收角),準直儀的有效孔徑長度增加,從而在不犧牲訊號靈敏度的情況下提高成像解析度。此外,在本發明的一些實施例中,準直儀的孔徑具有重複的圖案。在實施例的進一步中,圖案可以是均勻冗餘陣列(URA)圖案、修改的均勻冗餘陣列(MURA)圖案、完美二進制陣列(PBA)圖案、隨機圖案、偽隨機圖案和其他合適的圖案。重複的圖案允許基於輻射的成像系統採用簡化的影像重建演算法以降低計算複雜度並提高系統性能。在各種實施例中,準直儀還可以透過具有限定孔徑而不完全填充封閉孔的薄板而採用輕質設計。因此,新的準直儀設計從各個方面提高了基於輻射的成像系統的性能。
圖1A~1C示例了結合根據一些實施例的本發明的特徵的範例基於輻射的成像系統100(特別是核子醫學成像系統)。圖1A是系統100在X-Z平面中的一部分的橫截面側視圖。圖1B和1C是兩個實施例中系統100在Y-Z平面中的軸向圖。成像系統100可用於醫學檢查或治療諸如患者的目標物體。成像系統100包括整合高架102,該整合高架102還包括旋轉機構(例如,圍繞高架中心孔定向的轉子),該旋轉機構被配置為支撐一個或多個檢測器108和圍繞一個軸向(例如,如圖所示的X方向)旋轉一個或多個檢測器108(圖示出了處於相反位置的兩個檢測器108)。在如圖1B所示例的一個實施例中,兩個檢測器108可操作以在對立的位置處沿圓152旋轉,兩個檢測器108的各自頂表面的法線方向154A和154B均指向圓152的中心並朝向相反方向。在如圖1C所示的替代實施例中,兩個檢測器108類似地可操作以在對立的位置處沿圓152旋轉,兩個檢測器108的相應頂表面的法線方向154A和154B也指向圓152的中心。不同之處在於圖1B的法線方向154A和154B並不完全相反,而是形成一個小角度。根據一些實施例,取決於系統性能需要,角度可以在從大約0.5°到大約20°的範圍內變化。有關圖1B和1C中的設置的詳細資訊,稍後將會討論。每個檢測器108與準直儀110協同工作。準直儀110是引導光子路徑的裝置。在分子成像中,光子可能來自受試者內部的未知位置,這與X射線或CT中的光子是從已知的來源(或多個來源)位置發射不同。在沒有準直儀110的
情況下,來自所有方向的光子可能被檢測器108記錄,並且影像重建可能變得困難。因此,部署準直儀110來引導光子路徑,從而可以重建影像。成像系統100還包括耦合到機台支座系統114的患者檢查台112,該機台支座系統114可以直接耦合到地板或可以透過底座耦合到高架102。患者檢查台112被配置為能夠相對於機台支座系統114滑動,這有利於患者150進入和離開到與軸向軸線實質上對齊的檢查位置。控制台120例如以本領域已知的任何方式提供成像系統100的操作和控制。例如,操作員或技術員可以使用控制台120來控制機械運動,例如旋轉轉子104、移動、旋轉或傾斜檢測器108和準直儀110,以及滑動患者檢查台112。成像系統100還包括用於獲取資料和重建核子醫學影像的電腦組件(未示出),例如資料儲存單元、影像處理器、影像儲存單元和顯示器。在一些實施例中,一個或多個電腦組件可以部分地或全部地位於遠端位置(例如,使用雲端計算)。在一些實施例中,這些組件中的一個或多個可以存在於本地或遠端。
在一些實施例中,檢測器108是半導體檢測器,例如基於碲化鎘(cadmium telluride,CdTe)、碲鋅鎘(cadmium zinc telluride,CZT)或高純度鍺(high purity germanium,HPGe)的檢測器。在一些實施例中,檢測器108是閃爍體(例如基於碘化鈉(sodium iodide,NaI)或碘化銫(caesium iodide,CsI))檢測器。在一些其他實施例中,檢測器108也可以是與緊湊型光電倍增管(photo multiplier tubes,
PMT)、矽光電倍增管(silicon photomultiplier tubes,SiPMT)或雪崩光電二極體耦合的閃爍體。口服或注射到患者150中的一種或多種放射性藥物經歷核衰變,並且可以透過直接或間接地互毀以特定速率和特徵能量放射輻射(例如,伽馬光子)。檢測器108被放置在患者150附近以記錄或監測放射。基於記錄的資訊,例如位置、能量和此類檢測事件的計數,可以重建放射性藥物分佈的影像以研究患者150的某些身體部位的狀態或功能。
準直儀110包括限定一個或多個孔徑(也稱為開口或通孔)的多個壁(也稱為隔片)。在各種實施例中,隔片由諸如鉛或鎢之類的重金屬製成。取決於光子的能量的隔片的厚度,係足夠大以阻止大部分輻射,從而使光子主要地通過板上的小孔徑。厚度需要更大才能成像更高能量的伽馬射線。準直儀110被放置在檢測器108和成像物體(例如患者150)之間。準直儀的孔徑決定輻射可以穿過的方向和角跨度(接收角)以到達檢測器上的特定位置。根據孔徑的數量和幾何佈局,準直儀110可以是單孔準直儀、多孔準直儀、或編碼孔徑準直儀、或其他合適類型的準直儀。
成像系統100可以包括成像高架的其他必要部件,例如將部件耦合在一起(例如,將檢測器108和準直儀110連接在一起)的連接器、導致部件移動的馬達、光子屏蔽組件、包含其他部件的外殼組件等。例如,耦合和屏蔽組件116可以連接檢測器108和準直儀110,使得兩者一
起移動(例如,旋轉),並且防止輻射(光子)透過除準直儀110之外的路徑到達檢測器108。在其他實施例中,檢測器108和準直儀110可以相對於彼此單獨移動。
圖2A示出了範例多孔(或多孔徑)準直儀202的透視圖,並且圖2B示出了準直儀202的截面圖(沿圖2A中的B-B剖切)。準直儀202具有多個孔208,每個孔具有高度H和直徑(或寬度)W。高度H也是準直儀202的厚度。孔208也稱為孔徑、開口或通孔。孔徑208的接收角θ由從孔徑208上部開口上的邊緣點傳播到孔徑208下部開口上的對立邊緣點的兩條射線R0和R1之間的角度定義。只有入射光子(在接收角θ內傳播的光子或輻射)可以穿過孔徑208(不考慮可能穿透隔片212的一小部分射線)。如上所述,隔片212由吸收輻射的重金屬或合金製成,例如鉛或鎢。隔片212吸收大部分不從感興趣的方向散發出(或沿其行進)的輻射。用於較高能量的輻射的準直儀比用於較低能量的輻射的準直儀具有更厚的隔片。隔片通常設計為使不想要的光子穿透隔片不超過5%,在某些情況下不超過1%。應當注意,被準直儀阻擋或吸收的輻射或光子不需要100%的阻擋率,因為少數百分比的光子(例如,5%或更少)仍然可以穿透輻射吸收材料的厚度。換言之,阻擋(或其他類似術語)意味著絕大多數光子(例如,95%或更多,或99%或更多)被輻射吸收材料吸收。
準直儀202透過限制入射光子接收角θ來提供檢測到的光子的位置資訊。通常,接收角θ越小,準直儀
202提供的成像解析度越高。準直儀202可具有大量並排且彼此平行放置的實質上相同的、長而窄的孔徑208。長而窄的孔徑具有較小的接收角θ和相應的高成像解析度。準直儀202可以具有大於15mm的厚度H,例如20mm至70mm,具有大約1mm至5mm的孔徑直徑(或寬度)W,從而提供大寬高比(H/W)和小接收角θ。平行孔準直儀202的接收角θ被限制在一個小角度,例如不大於15°。如果接收角θ太大(例如大於15°),則特定平行孔準直儀提供的成像解析度被認為是不合需要地低。除非另有說明,本發明中的範例準直儀的接收角θ不大於15°。除了圖2B中的孔在“x-z”截面為矩形,接收角的概念可以擴展到其他形式的孔,例如在聚合式或散焦式準直儀中的孔,其中孔可以是梯形的,也可以是傾斜的,代表可以通過孔的射線的最大角度差。
圖2A所示例的準直儀202是具有均勻分佈的孔徑的平行孔準直儀。在所示例的實施例中,孔徑208實質上相同,除了由製造精度限制引起的變化之外,每個孔徑具有相同的高度H和相同的直徑或寬度W。因此,每個孔徑具有相同的較小的接收角θ(例如不大於15°,例如在特定範例中約5°)。準直儀202上的孔徑間距P也實質相同。換言之,準直儀202上的孔徑間隔S(隔片厚度)是均勻的。在此,孔徑形狀、孔徑尺寸、孔徑長度和孔徑間距均一致的孔徑被稱為“均勻分佈的孔徑”。
圖2C和2D示例了具有均勻分佈的孔徑的平
行孔準直儀(例如準直儀202)的兩種佈置的俯視圖,而其他的佈置也是可能的並且並非脫離本發明的精神和範圍。在圖2C中,孔徑208被佈置成正交的行和列以形成正交的二維網格佈置。每個孔徑208可被視為與網格網路(由圖2C中的虛線表示)中的方形單元網格240對齊。兩個相鄰的孔徑208具有相同的間距P,其中P=W+S,W是孔徑208的寬度(或直徑),S是兩個相鄰孔徑208之間的最小間隔(或隔片厚度)。在圖2D中,孔徑208以交錯的方式佈置,並呈現一連串的行彼此相鄰並且兩個相鄰的行彼此偏移一定距離(或偏移量)。在一個實施例中,偏移量可以是W,或者在替代實施例中可以是其他合適的值。每個孔徑208可被視為與網格網路(在圖2D中由虛線表示)中的六邊形形狀的單元網格240對齊。兩個相鄰孔徑具有相同的間距P,其中P=W+S,W是孔徑208的寬度(或直徑),S是兩個相鄰孔徑208之間的最小間隔(或隔片厚度)。在一些替代實施例中(未示出),孔徑208可以與網格網路的蜂窩形狀的單元網格240對齊。圖2E和2F分別示例了圖2C和2D的變體。單元網格240中的一個或多個可設計成沒有孔徑穿過,這相當於具有藉由隔片材料填滿(或封閉)的孔徑(由圖2E和2F中的虛線圓圈表示)。這種孔徑佈置被認為具有變化的孔徑間距或具有封閉的孔,因此不屬於具有均勻分佈的孔徑的平行孔準直儀的類別。換句話說,圖2E和2F中的孔徑是非均勻分佈的。
圖3是如圖2B所示例的準直儀202以及與其協
同工作的檢測器204以產生表示目標物體206(例如,患者)內的放射性藥物分佈的影像的示意性截面圖。準直儀202位於目標物體206和檢測器204之間並且被配置為透過阻擋某些光子並通過其他光子來過濾輻射。在所示例的實施例中,為了簡單起見,省略了連接準直儀202和檢測器204的耦合和屏蔽組件。在本實施例中,準直儀202為平行孔準直儀,其中準直儀202中的孔徑或孔具有矩形截面並且彼此平行佈置。在替代實施例中,準直儀202中的孔徑或孔可具有發散橫截面、聚合橫截面或一些其他橫截面形狀。
檢測器204包括感測像素214的陣列,例如矩形陣列、方形陣列或其他合適的像素214陣列。在操作中,每個感測像素214單獨記錄或監測入射在其上的輻射量並且產生與輻射量相關的訊號(例如,電壓或電流)。該些感測像素214可以是具有實質相同的尺寸和相同的形狀(例如,圓形、矩形或正方形)。在各種實施例中,感測像素214的尺寸範圍可以從大約1×1mm2到大約5×5mm2。在各種實施例中,陣列的感測像素間距P'的範圍可以從小於1mm到大約6mm。在一個範例中,基於CZT或矽倍增器(silicon multiplier,SiPM)的檢測器204可以被製造為4cmx4cm的尺寸,其包括具有2.5mmx2.5mm的單元像素尺寸的16x16感測像素的陣列。在一些實施例中,檢測器204包括適當的電子電路(例如,ASIC)以收集和處理從感測像素214產生的訊號。在一些其他實施例中,感測像素214在實體上彼此沒有區別並且僅僅是連續檢測器表面的數位化
的結果,就像大多數基於PMT檢測器系統一樣的情形。例如,檢測器表面可以被數位化為與準直儀202中的網格網路相同的網格網路,其中檢測器表面中的每個單元網格對應於準直儀202中的孔徑。
在圖3中,為了簡單起見,假設感測像素間距P'實質上等於孔徑間距P,使得每個感測像素214直接位於對應的孔徑208的下方。此外,準直儀202緊密耦合(或附接)到檢測器204並在準直儀202和檢測器204之間具有最小間隔,使得光子僅被允許沿著孔徑208行進以到達緊鄰其下方的對應感測像素214,儘管最小的或可忽略的串擾是可能的(例如,光子可能穿過隔片或光子可能會穿過一個孔徑,但會直接擊中另一個不同的孔徑的下方的感測像素)。在本發明中,術語“緊密耦合到”和“附接到”均指準直儀和檢測器之間的間隔小於準直儀厚度的一半,包括準直儀202和檢測器204是實體接觸的實例。
對於平行孔準直儀,當成像來源遠離準直儀時,其成像解析度迅速降低。對於孔徑尺寸W和長度(或有效長度)為He的準直儀,對於位於準直儀頂表面210上方垂直距離Z處的目標物體206,準直儀的成像解析度Rc由方程式(1)給定如下:Rc=(He+Z)/He*W=(1+Z/He)*W. (1)有效孔徑長度He是由He=H-2/μ給定,其中H是孔徑長度,μ是準直儀材料的線性衰減係數。對於鉛,在150keV時,μ=22.43cm-1,並且2/μ1mm。孔徑長度H通常大於20mm,
因此HHe。為簡單起見,H和He在本文中可互換使用,除非它們具有有意義的區別。作為以簡單的方式查看方程式(1),成像解析度Rc等於從孔徑底部透過孔徑可以看到的距離Z處的區域,如圖3所示。
從方程式(1)可以看出,成像解析度Rc隨距離Z線性增加。並且在相同距離Z下,如果He越大,Rc越小(亦即具有更好的成像解析度),這意味著較長的孔徑(亦即,較大的H)的成像解析度會提高。然而,較長的孔徑會以平方反比關係對訊號靈敏度產生不利影響。令G表示準直儀效率,被定義為通過準直儀的輻射射線(例如伽馬射線)量與光源發出的輻射射線量之比,則G與He -2成正比(亦即,G He -2)。因此,必須選擇孔徑長度H以平衡成像解析度和訊號靈敏度。
圖4示例了稱為“擴展場成像(spread field imaging,SFI)準直儀”的替代準直儀設計,其提供增加的等效孔徑長度,因此提供更高的成像解析度,而不犧牲訊號靈敏度。與圖3所示的平行孔準直儀202做比較,圖4所示的SFI準直儀203至少有兩個顯著特徵。首先,SFI準直儀203和檢測器204之間的間隔D與圖3中的配置相比顯著地增加。其次,SFI準直儀203的孔徑208並不完全相同並且至少在一個方面不同。這些將在下面做進一步討論。
增加SFI準直儀203和檢測器204之間的間隔D以允許穿過一個孔徑208的光子入射到檢測器204的延伸超出孔徑208的區域上。換言之,來自目標物體206的輻射穿
過孔徑208並傳播和照明檢測器204上的一個區域,該區域實質上大於孔徑208的尺寸並且延伸到其他相鄰孔徑的正下方的區域。相比之下,當準直儀202緊密耦合到檢測器204時,如圖3所示,來自目標物體206的穿過一個孔徑208的輻射僅照明檢測器204上實質上等於孔徑208的尺寸的區域。如圖4所示例,在定義接收角θ的線A1和A2之間入射的入射光子將穿過孔徑208並到達檢測器204並且“照明”檢測器204。由定義孔徑208的接收角θ的線所界定的檢測器204的表面區域被表示為相應孔徑208的照明區域。例如,中心孔徑208的照明區域由與檢測器204相交的線A1和A2之間的檢測器204的表面區域所限定。照明區域的大小和形狀取決於相應孔徑的輪廓以及準直儀202和檢測器204之間的間隔D。在一些實施例中,檢測器204上的孔徑208的照明區域至少是相對應的單元網格的尺寸的兩倍,該單元網格包括相對應的孔徑208和圍繞該相對應的孔徑208的隔片的一半(注意,單元網格的邊界位於介於單元網格與相鄰單元網格之間的隔片的中間)。此外,一個孔徑208的照明區域可以與其相鄰孔徑的照明區域重疊。例如,在各種實施例中,照明區域的超過10%、超過20%或甚至50%(僅考慮穿過孔徑208的輻射並且排除穿過圍繞孔徑208的隔片的輻射)的部分可以與其他的照明區域重疊。圖4示例了介於由孔徑208下方的線A1和A2所限定的照明區域與由相鄰孔徑下方的線A1'和A2'所限定的另一個照明區域之間的重疊區域。換言之,檢測器204的重疊區域中的點(或感測像
素214)可以接收來自不止一個孔徑208的輻射。此外,在一些實施例中,檢測器204中的任何感測像素214被不大於孔徑208的數量的百分比同時照明,其中該百分比小於約25%。在一些實施例中,該百分比小於10%、小於5%或甚至小於2%。在一些其他的實施例中,任何孔徑的照明區域不超過所耦合的檢測器204的總照明區域(藉由所有孔)的百分比,其中該百分比小於10%、小於5%,或者甚至小於2%。相比之下,編碼孔徑或多針孔準直儀中的孔通常比本發明實施例的總照明區域照明更大百分比的總照明區域。例如,那些準直儀中的孔的照明區域可能大於相應準直儀中所有孔的總照明區域的25%甚至大於40%,這將大量串擾引入成像系統。透過將任何孔的照明區域限制為不超過所耦合的檢測器的總照明區域(藉由所有孔)的10%、5%或甚至2%(取決於系統性能需求),來自不同孔徑的串擾被有效減小,並且因此可以實現良好的成像解析度和訊號靈敏度。
沿外圍區域提供屏蔽244(通常具有重金屬材料)以覆蓋(或密封)準直儀203和檢測器204之間的空間以防止輻射通過該空間到達檢測器204。值得注意的是,雖然圖4示出了屏蔽244和SFI準直儀203之間以及屏蔽244和檢測器204之間的間隙,其僅用於說明的目的。在各種實施例中,屏蔽244被部署以提供封閉(或密封)的區域,使得只有通過準直儀203的輻射才能到達檢測器204。
在本實施例中,SFI準直儀203的孔徑208並
不完全相同。它們的孔徑形狀、孔徑大小甚至孔徑長度可能不同。如果孔徑大小不同,則接收角θ可以變化,但在本實施例中最大接收角θ保持不大於15°。如上所述,不大於15°的接收角θ提供了理想的成像解析度。此外,對於給定的準直儀厚度H,較小的接收角θ轉化為較小的孔徑尺寸,從而導致較小的準直儀維度。不大於15°的接收角θ在緊密的準直儀設計和製造困難度(例如,機械公差)之間提供了良好的折衷。此外,大於15°的接收角θ會增加來自相鄰孔徑的串擾,並導致解析度降低以及成像後處理的複雜性增加。在一些應用中,接收角θ可以小於約10°,甚至小於約5°,以便基於系統性能需求實現更高的成像解析度和更少的串擾。在一些實施例中,SFI準直儀203的一些孔徑208可以是關閉的,這意味著其實質上是被阻擋的並且不允許光子通過。因此,孔徑間距也可以跨越準直儀203而變化。舉例來說,如圖5B所示,208'附近的間距和隔片與其他隔片相對地不同。這種在孔徑輪廓(例如,孔徑形狀、孔徑尺寸、孔徑長度及孔徑間距)的至少一個參數中具有非均勻性的孔徑被稱為“非均勻分佈的孔徑”。圖5A示例了具有不同孔徑形狀(例如,圓形和正方形)、不同孔徑尺寸和不同孔徑間距(由於由虛線圓圈表示的被阻擋的孔徑)的SFI準直儀203的實施例。在圖5A中,每個孔徑208與相應的單元網格240對齊(例如,中心到中心對齊)。圖5B示例了SFI準直儀203的又一實施例,其中至少有一個孔徑208'偏離其相應的單元網格240。具有這種孔徑佈置的SFI
準直儀203也被視為具有非均勻分佈的孔徑的平行孔準直儀。注意,圖5A和5B僅示出了SFI準直儀203的一部分。在一個實施例中,準直儀203上的孔徑208的總數超過一千(1,000)個,例如大約兩千(2000)、五千(5000)、一萬(10,000)、三萬(30,000),甚至十萬(100,000)。大量的孔徑有利於更高的靈敏度、良好的視域(field of view,FOV)覆蓋範圍以及更好地利用檢測器區域,因為每個孔徑僅照明檢測器204的表面區域的一小個百分比並且存在著顯著的串擾或在孔徑208之間多工。
在一些實施例中,SFI準直儀203中的孔徑208的一部分或全部可以以重複的圖案佈置,例如尺寸為3×3、4×4、3×5、5×5等等的重複的圖案。重複的圖案被稱為基本圖案。在一些實施例中,基本圖案具有奇數行和/列,使得基本圖案中存在一個中心孔,這在處理資料時是有利的。圖6A示例了SFI準直儀203的實施例的頂視圖(或頂視圖的一部分)。SFI準直儀203的至少一部分包括5×5單元網格的基本圖案的重複部分,具有在圖6B中所示的一個基本圖案,其中白色像素代表開孔徑,灰色像素代表閉孔徑。開孔徑被隔片包圍。圖6A中的虛線表示基本圖案之間的邊界。在所示例的實施例中,孔徑與正方形的單元網格對齊並且可以擴展到整個單元網格的大小範圍。基本圖案可以是編碼孔徑圖案,例如URA、MURA、PBA、隨機或偽隨機圖案。在一些實施例中,重複的圖案的數量大於2×2,例如是3×3、3×5、5×5等。在圖6A中,所示例的基本圖
案是MURA5圖案並且在形成尺寸為15×15的準直儀(或準直儀的一部分)時在X方向上重複3次並且在Y方向上重複3次。在各種實施例中,一個基本圖案可以分別在X方向上和Y方向上重複任意次數,以滿足系統需求。圖6C示例了SFI準直儀203(或其一部分)的另一個實施例的俯視圖,其中在形成尺寸為125×125的準直儀(或準直儀的一部分)時MURA5的基本圖案在X方向上重複25次並且在Y方向上重複25次。在所示例的實施例中,孔徑的形狀和尺寸是相同的。在一些實施例中,孔徑可以在形狀、尺寸或兩者上不同。
返回參考圖4,可以選擇SFI準直儀203和檢測器204之間的間隔D,使得檢測器204上的每個孔徑208的照明區域等於或不大於重複的基本圖案但大於孔徑的間隔(間距)至少兩倍。檢測器204上的孔徑208的照明區域的大小與SFI準直儀203的基本圖案的大小實質相等,如圖4中虛線B1和B2所示例的。如下文將進一步詳細討論的,使照明區域的尺寸實質上等於(或為基本圖案的整數倍)基本圖案提供了一種快速影像重建演算法,該演算法顯著提高了系統效率和準確度。結果,當SFI準直儀203和檢測器204之間的間隔D大於孔徑長度(準直儀厚度)H的一半時,它被優化。在各種實施例中,D可以是小於H的10倍,例如介於H的1.5至10倍、介於H的1至7倍之間、介於H的1.5至7倍之間,或一些其他合適的範圍。如果D太小,例如小於H的1.5倍,則對於解析度提高的影響不顯著。另一方
面,如果D大於H的10倍,來自相鄰孔徑的串擾會很高,並導致影像重建過程中的複雜性增加和解析度下降。在一些範例中,D是H的約1倍、1.5倍、2倍、2.5倍、3倍或3.5倍。
仍參考圖4,對於檢測器204的感測像素214,可以通過其正上方的孔徑看到的區域,也稱為成像解析度Rn,由方程式(2)給定如下:Rn=(H+D+Z)/(H+D)*W. (2)請注意,方程式(1)和方程式(2)之間的差異是方程式(1)中的He被方程式(2)中的(H+D)代替,根本上將孔徑長度延長了D,因此可以顯著提高成像解析度。即使為了更高的訊號靈敏度而增加孔徑的尺寸W,Rn也可以小於Rc。
光子可能會穿過一個孔徑,然後直接撞擊另一個孔徑下方的檢測器區域,這一事實意味著存在串擾或多工。相比之下,當配備平行孔準直儀(例如圖3所示的準直儀202)時,對於檢測器表面上的任何感應像素,只有一個孔徑光子可以穿過並到達它(不包括隔片穿透,這在平行孔準直儀設計中是要被抑制的效應)。這種串擾效應通常會降低解析度,因為在檢測器位置檢測到光子的可能路徑有多個,而對於位於某個點的來源,輻射射線可能穿過多個孔徑到達檢測器。如果所有孔徑與平行孔準直儀中的孔徑本質上地相同,就會出現這種情況。
另一方面,關於圖4所示的SFI準直儀203,孔徑不完全相同。更具體地,檢測器表面上的感測像素的
可見區域內的孔徑並不完全相同(此處的可見區域是指光子可以穿過以到達檢測器表面上的感測像素的所有孔徑的集合)。所以,彼此靠近的光源會在檢測器上投射出不同的陰影,因為它們下方的孔徑圖案不同。這些差異有助於重建演算法以類似於Rn的高解析度恢復原始光源的分佈,即物體影像。檢測器表面通常被數位化為網格網路,並且每個網格代表一個感測像素的邊界。在一些實施例中,由相鄰感測像素所見的孔徑圖案是不同的。
類似於SPECT成像,可以透過旋轉耦合的準直儀(例如準直儀203)和檢測器(例如檢測器204)的組來環繞目標物體從多個角度獲取影像(也稱為投影)。在這種獲取模式下,耦合的準直儀和檢測器一起移動,它們之間沒有相對運動。在某個角度α處的前向投影p可以寫成
其中i和j是檢測器像素和物體影像三維像素的索引,f(j)是三維像素物體j的值,Kα(i,j)是當相機處於角度α時,在像素i處檢測到從三維像素j發出的光子的概率。為了重建原始的物體影像f,使用的一種方法是最大相似度與期望值最大化演算法(Maximum Likelihood Expectation Maximization,MLEM)演算法,其中可以迭代地找到f
或者
其中f(k)(j)是第k次迭代的目標影像,並且p(i,α)是檢測器像素i和角度α處的測量的投影。該演算法分為三個步驟--前向投影、計算比值和反向投影。該演算法的一種變體稱為OSEM,其中投影被分成N個子集,計算從使用一個投影子集的幾次迭代開始,然後移動到下一個子集,直到計算完所有子集。該過程可以以不同的子集劃分並可以以每個子集的不同迭代次數重複。
如方程式(5)所示的演算法,考慮到大量的三維像素及像素並乘以角度的數量,方程式(5)非常耗時並且需要大量的儲存空間來儲存K矩陣。一個角度的投影可以數位化為64x64到512x512像素,通常需要從60或120個角度採集。並且物體影像通常具有128x128x128或256x256x256個三維像素。因此矩陣Kα(i,j)非常地大。
當孔徑由重複的基本圖案組成並且所有重複的圖案中的各個通孔具有相同的形狀和維度時(即使基本圖案中的孔可能不同),矩陣K和演算法可以被簡化。通常,對於大小為X×Y(X和Y可以不同)的重複的基本圖案,僅有些許圖案等於X×Y的乘積的倍數,X×Y係點光源可能投射到檢測器上的大小。如圖6A或6B所示是重複的圖案為5×5(X=Y=5)MURA 5圖案的範例。令m為基本圖案中孔(網格單元,包括開孔和閉孔)的索引,m=1,…25。在該範例中,考慮位於距離準直儀203距離Z處的平面(見圖4),該平面上任意一點的一個點光源將圖案投射到檢測器204上,其中投射的圖案的尺寸遠小於檢測器的大小,並且點
光源可以投射到檢測器上的圖案只有X×Y=25個圖案,這取決於它正下方(邊界區域除外)的孔徑指數,其被指代為psfm,Z(x,y),其中x=-M,…M以及y=-N,…,N。這裡檢測器像素的大小與準直儀間距(孔間間隔)相同或與它的一小部分相同。因此,該25個圖案可以透過基本的5x5圖案中的孔徑而被編索引。同時,距離Z處的物體影像,記為fZ(其像素(也稱為三維像素)大小與準直儀間距(孔間間隔)相同並且像素與準直儀間距對齊)可被分為25個與fZ大小相同的子影像,其中每個子影像中相同索引的孔徑頂部的像素保持不變,其餘像素設置為0,如圖7所示例。將m設為圖案索引,則距離Z處的物體影像切片,fZ,可以被指代為由以下所給定的子影像,fZ(j,m),式子
這樣,方程式(5)中Z處物體影像的前向投影的步驟可以重寫為
其中(Z)是基於物體影像fZ (k)的當前估計值而計算的前向投影,*表示卷積。可以採用類似的方法來計算反向投影步驟。眾所周知,卷積在計算上比使用快速傅立葉變換(Fast Fourier Transform,FFT)演算法的矩陣乘法更有效。這個想法可以擴展到相當於準直儀間距的一部分的三維像素大小,其中可以使用psf圖案和相應子影像的X×Y卷積的倍數。例如,如果三維像素的大小是準直儀間距的一半,
則可以使用psf圖案和相應子影像的25個卷積的4(2×2)倍(即總共100個卷積)。這種方法可能會隨著計算複雜度的增加而帶來更高的準確性。
方程式(7)還指出,如果物體影像被組織成平行於由距離Z所指代的準直儀表面的切片,那麼對於不同的角度α,投影圖案psfm,z是相同的。這進一步簡化了演算法。
一種方法是對每個角度獨立地實現方程式(5),從總和中除去α;然後沿距離Z對估計值fZ(j)求和以獲得該角度的二次投影,p'α。由於從物體影像f導出p'α的投影幾何類似於平行孔準直,平行孔準直影像重建演算法可用於從p'α重建f,例如濾波反投影(filtered back-projection,FBP)、代數重建技術(algebraic reconstruction technique,ART)和OSEM方法。
另一種方法是在每次迭代開始時透過內插法將fZ轉換為fα,z,其中fα,z由在該角度平行於準直儀表面的切片(由Z標記)組成。然後可以在應用內插法之前執行投影-反投影步驟(projection-back-projection steps)來將反投影比值,qa,z,其由下式給定
轉換為f中的原始網格,然後計算fZ的更新比值。
回過頭來參考圖4,對於準直儀203,孔徑的高度H(亦即準直儀203的厚度)有時遠大於其開口寬度或直
徑W。例如,孔徑寬度W可以小於2mm,而孔徑高度H可以大於20mm。然而,對於封閉的孔徑,孔徑不需要完全填充,而是填充到有足夠的厚度以阻擋輻射穿過孔。例如,當放射性同位素Tc-99m用於隔片材料時,3mm厚的鉛可能會阻擋140keV伽馬射線的入射輻射的99%以上。因此,準直儀可以以重量輕的形式製造。參考圖8A和8B,準直儀203是具有開孔和閉孔的準直儀並且包括兩個區段。一個區段類似於具有均勻分佈的通孔的平行孔準直儀(也稱為平行孔區段262)。另一個區段實質上是薄板264(具有足以實質上阻擋輻射的厚度),其具有不均勻分佈的通孔,例如具有比平行孔區段262中的通孔少的通孔(例如,僅為其大約一半的通孔數)。在圖8A和8B中示出的準直儀203也稱為多區段準直儀。薄板264中的通孔與平行孔區段262中的通孔對齊以共同限定準直儀203的孔徑圖案,而平行孔區段262中的其他通孔被薄板264的不透明材料阻擋。因此,薄板264中的通孔限定了準直儀203的最終孔徑的位置。薄板264也被稱為圖案化孔徑區段264。在一些實施例中,平行孔區段262中的所有通孔實質相同,而圖案化孔徑區段264中的通孔可以不同,甚至包括封閉的通孔。此外,在一些實施例中,圖案化孔徑區段264可以是可移動的,例如可操作以相對於平行孔區段262的表面水平滑動。因此,透過以相對於平行孔區段262的方式移動圖案化孔徑區段264可以獲得在一個角度的多個投影。
圖案化孔徑區段264具有足夠厚以阻擋目標
輻射的厚度H1。在各種實施例中,平行孔區段262的厚度H2可以是編碼孔徑區段264的厚度H1的大約2倍到大約10倍。例如,25mm的孔徑長度可以構造成具有平行孔區段262(具有20mm厚度)以及圖案化孔徑區段264(具有5mm厚度)。透過這種方式,與完全填充的封閉孔徑相比,準直儀的總重量可以顯著降低。圖案化孔徑區段264可以位於平行孔區段的面向檢測器(圖8A)或面向目標物體(圖8B)的一側,或者甚至在兩個子平行孔區段的中間(例如,圖8C中的子平行孔區段262-a和262-b,其中H2-a+H2-b=H2)。事實上,封閉的孔徑只需要一層薄薄的重金屬材料就可以阻擋輻射。這個薄層甚至可以根據需要分成多個片段。因此,這種新設計可以從製造正規的平行通孔開始,然後可以在開啟的通孔內的任何地方放置薄層以製造封閉的孔徑。
儘管在圖8A~8C中所示例的孔徑具有相同的橫截面和開口寬度,但還是存在許多不同的設計以允許不同的孔徑輪廓。例如,輕量化設計還可以應用於具有具備相同間距但具有不同孔徑尺寸和隔片的孔徑的準直儀。一些這樣的例子在圖9A~9C中示例出。在圖9A中,準直儀203具有一些直徑(或寬度)W的孔徑和具有一些較小的直徑(或寬度)W'的孔徑。將在準直儀203的整個厚度上對較窄的孔徑的側壁使用恆定寬度的隔片的作法,以隔片可僅在側端處具有較厚的部分代替。側端處的厚度H1-a和H1-b足以阻擋目標輻射並保持與一個較窄的孔徑所被設計為應具
有的較小的接收角θ相同。以上參考圖8A~8C的多部分設計的精神也可以應用於圖9A中的範例準直儀。圖9B示例了具有一個平行孔區段262和兩個限定孔徑尺寸的圖案化孔徑區段264-a/264-b的準直儀203。在某些情況下,僅需要264-a/264-b其中一者以調整接收角θ。此外,孔徑可以在橫截面和垂直截面兩者中具有不同的形狀。在一些實施例中,孔徑可以具有不同的定向。圖案化孔徑區段264-a/264-b可以與開啟和閉合的孔徑圖案一起定義孔徑尺寸,例如圖9C所示。共同參考圖8A~9C,一個平行孔區段262可以與不同的圖案化孔徑區段264a及/或264b配對以創造不同的準直儀配置。因此,根據本發明的多區段準直儀,其除了實現輕量化設計之外還提供了低成本的解決方案。值得注意的是,圖8A~C和9A~C描繪了多區段準直儀的相鄰區段彼此實體接觸的實施例。在替代實施例中,只要準直儀的性能得到保證,相鄰區段係彼此緊密耦合(以防止光子透過一個孔進入而從另一個孔離開)而沒有實體接觸。
現在參考圖10A。不同於圖8A~9C中所示例的範例性準直儀,在圖8A~9C中平行孔區段262和圖案化孔徑區段264彼此附接,圖10A所示的範例性準直儀具有與圖案化孔徑區段264間隔開的平行孔區段262。此外,圖案化孔徑區段264可緊密耦合到檢測器204。在所示例的實施例中,檢測器204的感測像素214具有在圖案化孔徑區段264中的對應的孔徑以接收來自平行孔區段262中更上方的
平行孔的照明。圖案化孔徑區段264中的孔徑有意地在尺寸上不同,從而導致平行孔區段262中可照明像素214的平行孔的數量不同。例如,感測像素214a在圖案化孔徑區段264中具有比相鄰感測像素214b更大的對應孔徑。如此一來,感測像素214a可被平行孔區段262中的多個平行孔照明,平行孔區段262介於與平行孔區段262的底表面相交的線C1與C2之間。作為比較,由於在圖案化孔徑區段264中具有較小的相應孔徑,像素214b只能被位於平行孔區段262正上方的單個平行孔照明,該平行孔區段262由與平行孔區段262的底表面相交的線D1和D2界定。平行孔區段262中可以照明同一個感測像素214的平行孔數量取決於圖案化孔徑區段264中對應的孔徑的接收角。平行孔區段262中能夠照明同一個感測像素214的平行孔數量也取決於介於平行孔區段262與圖案化孔徑區段264之間的距離H3。距離H3是基於對系統性能的需求而決定。一個通常的規則是設置一個能夠允許檢測器204的一些像素214被其上的平行孔208中至少兩個或更多個照明到的距離H3。值得注意的是,在如圖10B所示的替代實施例中,可以將平行孔區段262配置為緊密耦合到檢測器204,而圖案化孔徑區段264卻與平行孔區段262間隔開。
本發明中提出的設計可以適用於具有大量長而窄的孔的其他變體,例如聚合式或散焦式孔準直儀,其用前面提出的兩個特徵修改。
返回參考圖1A~1C,基於輻射的成像系統100
可以採用至少兩組不同的耦合的準直儀和檢測器。在這樣的實施例其中之一,一組是具有與耦合的檢測器間隔開的SFI(spread-field imaging)準直儀203(例如,如圖4、6A~C、8A~C、9A~C或10A~B中所示)並且另一組是具有附接到耦合的檢測器的平行孔準直儀202(例如,如圖3所示)。通常在這種配置中,平行孔準直儀202中的孔徑開口(及/或接收角)小於SFI準直儀203中的孔徑開口(及/或接收角)。如方程式(1)中所示,由於較小的孔徑開口,平行孔準直儀通常在近距離(小Z)處具有良好的解析度,但隨著目標物體遠離準直儀,其成像解析度迅速下降。相反地,當目標物體離準直儀較遠時,SFI準直儀203提供優異的成像解析度。在多角度影像採集過程中,兩組耦合的準直儀和檢測器模組繞著目標物體旋轉以獲取影像(投影)。理所當然地,目標物體的某個部分在某些角度靠近一個準直儀,而在相對的角度則遠離一個準直儀,除了靠近旋轉中心(center of rotation,COR)的部分。透過組合這兩個準直儀,對於較接近平行孔準直儀202的位置,使用者可以分配較多的權重來更新由平行孔準直儀202投影所估計的因子,而對於離較平行孔準直儀202較遠的位置,使用者可以分配較多的權重來更新由SFI準直儀203投影所估計的因子。換句話說,該兩組耦合的準直儀和檢測器模組的估計之間的權重是依賴於距離的。這種配置的另一個優點是可以透過單獨使用平行孔投影得到的重建來獲得物體影像的快速估計,並且可以將其用作使用來自兩個準直儀的投影進行迭代重建的初始估計。有幾種方法可以選擇來重建多
角度平行孔投影,例如濾波反投影(filtered back-projection,FBP)、代數重建技術(algebraic reconstruction technique,ART)和有序子集期望及最大化(ordered subset expectation and maximization,OSEM)。並且可以執行該步驟以在較低解析度下重建影像(透過對投影應用低通濾波(low-pass filtering,LPF),及/或為較大像素重建)以減輕由於以下事實而導致的較低計數:與帶有雙攝影機的傳統掃描相比(其中兩個攝影機皆具有平行孔準直儀),此處的平行孔準直儀投影僅佔採集角度的一半。在另一個實施例中,系統採用至少兩組耦合的準直儀和檢測器模組,它們均包括SFI準直儀203(例如,如圖4、6A~C、8A~C、9A~C或10A~B所示),其與耦合的檢測器間隔開。但是兩個準直儀中長窄的通孔具有不同的接收角。例如,其中一個SFI準直儀203的最大接收角至少是另一個SFI準直儀203的兩倍,例如一個為10°,一個為5°。在另一個實施例中,系統採用至少兩組耦合的準直儀和檢測器模組,至少一個檢測器與具有長窄孔(小接收角)的SFI準直儀203耦合(例如,如圖4、6A~C、8A~C、9A~C,或10A~B所示),其中該長窄孔與耦合的檢測器間隔開,並且另一個檢測器與多針孔準直儀耦合,該多針孔準直儀包括編碼孔徑準直儀,其中最大孔徑接收角實質上大於SFI準直儀203的接收角,比如20°以上。
參考圖1B,如果基於輻射的成像系統100在對立位置處採用兩組相同的耦合的準直儀和檢測器,則每組耦合的準直儀和檢測器可能只需要以相等的掃描步長旋
轉半圓。透過組合來自第一半圓和第二半圓的資料,獲取來自全圓的影像。例如,在掃描步長為3°時,第一組可以依次獲取0°、3°、6°、...、177°角度的影像;第二組可以依次獲取180°、183°、186°、...、357°角度的影像。如果基於輻射的成像系統100採用兩組不同的耦合的準直儀和檢測器,由於每組所“看到”的不同,每組需要旋轉一整圈。為了保持相同的操作時間,每組可能以兩倍的掃描步長掃描一整圈。例如,以6°的掃描步長下,第一組可以依次獲取0°、6°、12°、……、180°、186°、……、354°等角度的影像;第二組可以依次獲取180°、186°、192°、...、354°、0°、...、174°等角度的影像。對於系統掃描的每個給定的角度,將獲取兩個影像,一個來自第一組,另一個來自第二組。換句話說,雖然這兩組所“看到”的不同,但每一組重複從另一組獲取影像的相同角度獲取影像。作為比較,參考圖1C,這兩組不同的耦合的準直儀和檢測器並不完全相反,而是偏移了一半的掃描步長。例如,如果掃描步長是6°,則這兩組指向的方向154A和154B可以形成大約3°的角度。在這兩組指向的方向之間加入偏移允許系統100從不同角度獲取每個影像。例如,以6°的掃描步長和3°的偏移量,第一組可以在0°、6°、12°、...、180°、186°、...、354°等角度依次獲取影像;第二組可以在183°、189°、195°、……、357°、3°、……、177°等角度依次獲取影像。
現在參考圖11A和11B。在基於輻射的成像系統中,當準直儀的孔徑與檢測器的感測像素對齊時,性
能得到優化。圖11A示例了疊加在檢測器204上的範例性SFI準直儀203(如圖6B所示的基本圖案)的一部分在X-Y平面中的俯視圖。圖11B示例了準直儀203和檢測器204的一部分沿圖11A中的A-A剖切在X-Z平面中的示意性截面圖。如圖10A和10B所示例的檢測器204包括離散感測像素214的陣列。離散感測像素214可以基於CZT或SiPM結構。通常,感測像素214排列在以行和列形成的網格網路(也稱為像素網格)中,例如在所示例的實施例中行i=1、2、3、...和列j=1、2、3、...。位於X-Y平面中的位置(x,y)處的任何感測像素214可以由網格網路中的索引(i,j)(i和j為整數)表示並且在準直儀203(也稱為孔徑網格)的網格網路中具有對應的孔徑(或著編碼圖案中的閉孔徑或一個隔片)。由於圖11A和11B中的感測像素214是離散的並且在實體上彼此不同,實現檢測器和準直儀之間的對準是直接的,亦即是將感測像素214直接定位在準直儀203的相應孔徑208下方,使得像素網格和孔徑網格重疊。
現在參考圖12A和12B。如上面結合圖3和4所討論的,在一些實施例中,檢測器的感測像素彼此在實體上並非有所區別並且僅僅是連續檢測器表面的數位化的結果,如在大多數基於PMT檢測器系統的情況下。圖12A示例了疊加在這樣的檢測器204上的範例性SFI準直儀203(如圖6B所示的基本圖案)的一部分在X-Y平面中的俯視圖。圖12B示例了準直儀203和檢測器204的一部分沿圖12A中的A-A剖切在X-Z平面中的示意性截面圖。檢測器
表面被數位化為與準直儀203的相應孔徑網格相同的像素網格,亦即,網格點之間的間隙相同的準直儀203。換言之,檢測器表面上的每個像素網格單元在最近的準直儀203中具有對應的孔徑(或著編碼圖案中的閉孔徑或一個隔片)。檢測器表面上的每個像素網格單元也稱為檢測器的一個像素。檢測器的像素排列在由行和列所形成的網格網路中,例如行i=1、2、3、...和列j=1、2、3、...並且可以標記為像素(i,j)(i和j是整數)。理想情況下,在完美對齊的情況下,檢測器的像素網格和準直儀的孔徑網格應該重疊。然而,如圖12A和12B所示,像素網格和孔徑網格之間經常會發生未對準,例如在設備組裝期間,導致在X方向(指代為Δx)及/或在Y方向(指代為Δy)上的未對準。未對準會導致系統退化。例如,如圖12B所示,在入射光子230穿過準直儀203的孔徑這樣的場景下(這原本應該被記錄為在完美對準下發生在像素(3,1)處的一個事件(例如,計數)),卻反而將被記錄為發生在像素(3,2)處的一個事件,因此降低了在記錄的訊號中孔徑圖案的對比度。
以基於光電倍增管(photomultiplier tubes,PMT)的檢測器系統為範例,其中基於其他結構的檢測器系統是類似的,簡單地說,PMT的作用是將入射光子在光電陰極處轉化為光電子。這些電子從一系列帶電的陰極產生大量的二次電子,在陽極處產生可測量的電流脈衝。基於在相鄰的PMT中的這些可測量的電流脈衝,可以將連續X-Y平面中的位置(x,y)確定為光子入射事件發生的位置。
基於連續X-Y平面和(i,j)索引之間的映射,發生在位置(x,y)的事件隨後被數位化為發生在像素(i,j)中的事件。然而,如上所述,由於未對準(其係依賴於連續X-Y平面和(i,j)索引之間的初始映射),可能會導致性能下降。影像對比度提供了一種調整以抵銷未對準的方法。為了使人眼和頭腦能夠感知影像,所獲取影像的像素陣列必須顯示對比度。關於樣本的某些東西必須使記錄在不同像素的訊號強度發生變化。最簡單的說,傳輸對比度可能是由於部分或完全不透明的結構造成的,例如準直儀的孔徑和封閉的孔徑(或隔片)。影像中存在的對比度數量決定了調整未對準的關鍵精度。透過演算法來抵銷未對準而不實體移動準直儀的設置在圖13中示例出。在圖13中,諸如Co-57或Tc-99m矩形泛光源等的泛光源242在準直儀203上方提供均勻的輻射場。泛光源242可以定位成與準直儀203緊密耦合,亦即兩者間的距離要比準直儀厚度的一半還短,包括準直儀203和泛光源242實體接觸的情況。在泛光源242的照明下,從檢測器204獲取的影像是由準直儀203所遮蔽的泛光源242的影像。由於輻射均勻,所獲取的影像的對比度是由於準直儀203的透明和不透明特徵。換言之,與開放或大孔徑下方的像素相比,閉合孔徑(或隔片)或小孔徑下方的像素產生較弱的影像訊號強度。
對比度可以定義為介於峰值和谷值之間的影像單一強度的變動(例如,比值或絕對差)的量度。作為一個範例,圖13中的訊號強度線245示例了對應於峰值和谷
值的訊號強度Vp和Vv。最大化對比度的一種方法是最大化Vp/Vv的比值,或最大化|Vp-Vv|/Vp的差值。當對比度被最大化時,準直器和檢測器之間的對準被認為已實現。
由於檢測器表面被數位化為像素網格,因此無需實體移動準直儀即可透過演算法來達成抵銷未對準(換句話說,最大化對比度)。在一個範例性方法中,從檢測器表面(例如,如果在基於PMT的檢測器系統中,則來自PMT的電流脈衝)出現的每個訊號被記錄為具有座標(x,y)的事件。因此,資料列表(也稱為列表模式)中的每個條目都是一個屬性向量,其包含有關單個檢測到的光子的資訊。例如,事件流可以標記為(xn,yn)(n=1、2、3、...),其中每一對代表被檢測到的第n個事件的X-坐標和Y-坐標。基於連續X-Y平面和離散(i,j)索引之間的映射,像素(i,j)中報告的總計數被指代為C(i,j),其可以被表示為
,其中D是將(x,y)處發生的事件映射到相應像素(i,j)的數位化函數,如果D(xn,yn)等於(i,j),則δ函數給定1,其餘則給定0。此外,可以可選地應用能量窗口來過濾事件,使得僅當在能量窗口內發生事件時才被計數。根據定義,C(i,j)為像素(i,j)處的影像訊號強度,表示像素(i,j)覆蓋的區域內記錄到的事件數量,其由數位化函數D定義。峰值和谷值,Vp和Vv,處的訊號強度可以從C(i,j)(i=1、2、3、...;j=1、2、3、...)的集合中識別出來並且初始對比度(Vp/Vv或|Vp-Vv|/Vp)可以計算出來。
接下來,在數位化函數中引入偏移對(δx,δy)以“偏移”像素網格並重新計算C(i,j)。因此,C(i,j)可以被表示為,
其中(δx,δy)的引入相當於在X方向上將像素網格移動了距離δx,並且在Y方向上移動了距離δy。結果是,在位置(x,y)發生的事件(原本會記錄在像素(i,j)下),可能會被視為等效地發生在位置(x+δx,y+δy)處,並相應地記錄在相鄰像素(i',j')之下。因此,C(i,j)的曲線將會不同,並且計算更新的對比度(Vp/Vv或|Vp-Vv|/Vp)。可以透過最佳化方法(例如最陡下降(steepest decent)、共軛梯度(conjugate gradient)等)找到使峰值/谷值對比度最大化的最佳(δx,δy)。另一種方法是以小步長創造一系列偏移對(δx,δy)。透過掃描該一系列偏移對(δx,δy),可以選擇具有最大對比度的一個偏移對(δx,δy)作為正確的偏移量來反制未對準,例如透過重新產生像素網格。這個偏移對(δx,δy)被稱為選定偏移對或選定偏移。或者,可以應用選定的偏移對(δx,δy)來偏移X-Y平面的原點(x0,y0)以移動整個像素網格(重新產生像素網格的特殊實例)。理想情況下,選定的對(δx,δy)實質上等於未對準(Δx,Δy),但還有取決於用於產生該一系列偏移對(δx,δy)的掃描步長。在各種實施例中,可以透過分別等分像素間距Px(X方向上的間距)和Py(Y方向上的間距)來產生該一系列偏移對(δx,δy)。例如,對於在X-以及Y-方向上皆為2.0mm的像素間距,偏移
對(δx,δy)的集合可以是從0開始以0.2mm的步長並且δx<Px以及δy<Py進行掃描,亦即{(0,0),(0.2mm,0),(0.4mm,0)...(2.0mm,1.8mm),(2.0mm,2.0mm)}的一個集合。用其他方式產生偏移對(δx,δy)的集合也是可能的,例如透過泡式排序演算法(bubble sort algorithm)、桶排序演算法(bucket sort algorithm)、插入分類演算法(insertion sort algorithm)或其他合適的演算法。
代替同時掃描δx和δy,替代地,該方法可以單獨在X方向掃描(透過固定δy)以首先在X方向獲得最優化的偏移δx,然後再單獨在Y方向掃描(透過使用最優化的δx)以隨後在Y方向上獲得最優化的偏移δy。或者,類似地,該方法可以先在Y方向掃描,然後再在X方向掃描。此外,代替從整個像素網格二維地挑選Vp和Vv來計算對比度,該方法可以挑選一行(透過在C(i,j)的集合中固定i,例如圖13中的訊號強度線245)並掃描(δx,δy)以獲得該行中的最大對比度,並且將選定的偏移對(δx,δy)應用於整個像素網格。或者,類似地,該方法可以挑選一列(透過在C(i,j)的集合中固定j並掃描(δx,δy)以獲得該列中的最大對比度,並且將偏移對(δx,δy)應用於整個像素網格。
值得注意的是,不同的孔圖案可能需要不同的最優化方法。在一些圖案中,記錄的訊號強度的描述符(例如,對比度、峰值訊號強度或谷值訊號強度)的最大值或最小值可用於最大化/最小化,例如最大化對比度、最大化的峰值,或最小化的谷值。當使用重複的孔圖案時,
可以將相應的像素(對應於基本圖案中的相同孔)加在一起以降低像素值的隨機性(亦即,量子雜訊)。此外,在一些實施例中,當使用特定孔圖案時,來自最優化例程的選定的偏移對(δx,δy)可能與實際的最佳偏移間具有固定的偏移,並且該固定的偏移由孔圖案決定。例如,在某些孔圖案中,泛光源影像中的最大像素值位於相鄰的2x2開孔的中間。在這種情況下,固定的偏移是孔間距的一半,並且用於將像素與孔對齊的實際最佳偏移量將透過固定的偏移進行校正。換句話說,在找到選定的偏移對(δx,δy)之後,可以可選地執行進一步的步驟,以向(或從)選定的偏移對(δx,δy)添加(或減去)固定偏移,從而獲得實際的最佳偏移量以便用於未對準調整。
儘管不旨在進行限制,但本發明的一個或多個實施例為目標物體(例如患者)的基於輻射的成像提供了許多益處。例如,具有與相關的檢測器間隔開的重複的編碼圖案的準直儀設計可提供卓越的成像解析度,而不會犧牲訊號靈敏度。因此,提高了系統性能。
前述概述了幾個實施例的特徵,以便本領域普通技術人士可以更好地理解本發明的各方面。本領域普通技術人士應當理解,他們可以容易地使用本發明作為設計或修改其他過程和結構的基礎,以實現與本文在此介紹的實施例相同的目的及/或實現相同的優點。本領域普通技術人士也應該認識到,這樣的等效構造並不脫離本發明的精神和範圍,在不脫離本發明的精神和範圍的情況下,
可以對本文的揭露進行各種變化、替換和變更。因此,所附申請專利範圍被廣泛地並且以與本發明內容一致的方式來解釋是合適的。
100:成像系統
102:整合高架
108:檢測器
110:準直儀
112:患者檢查台
114:機台支座系統
116:耦合和屏蔽組件
120:控制台
150:患者
Claims (9)
- 一種基於輻射的成像系統,包括:準直儀,被配置為過濾從目標物體發出的輻射,該準直儀包含多個不均勻分佈在該準直儀上的孔徑,其中該多個孔徑的最大接收角不大於15°,並且該多個孔徑形成包含基本圖案的孔徑圖案,該基本圖案在該孔徑圖案內重複超過四次;以及檢測器,用於檢測已經通過該準直儀的輻射,其中,該準直儀與該檢測器間隔開,以使得面向該準直儀的該檢測器的頂表面上的點同時被該多個孔徑中的兩個或更多個孔徑照明,以及其中,該孔徑中的一孔徑的照明區域實質上等於該基本圖案的區域。
- 根據請求項1所述的基於輻射的成像系統,其中,該檢測器的該頂表面上的該點同時被該多個孔徑的一定百分比照明,其中該百分比約小於25%。
- 根據請求項1所述的基於輻射的成像系統,其中,該多個孔徑的數量超過一千。
- 根據請求項1所述的基於輻射的成像系統,其中,該多個孔徑至少在孔徑尺寸、孔徑形狀、接收角、長度和孔徑間距中之一方面變化。
- 根據請求項1所述的基於輻射的成像系統,其中,該準直儀與該檢測器之間的距離為該準直儀的厚度約1.5至約10倍。
- 根據請求項1所述的基於輻射的成像系統,其中,該基本圖案是編碼孔徑圖案。
- 根據請求項1所述的基於輻射的成像系統,其中,該準直儀包含第一部分和第二部分,其中,該第一部分包含均勻分佈在該第一部分上的通孔,以及該第二部分包含不均勻分佈在該第二部分上並且對應於該多個孔徑的通孔,以使得該第一部分的該通孔的一部分被該第二部分擋住。
- 根據請求項7所述的基於輻射的成像系統,其中,該第一部分的厚度為該第二部分的厚度約2至約10倍。
- 根據請求項1所述的基於輻射的成像系統,其中,該準直儀是第一準直儀並且該檢測器是第一檢測器,還包括:第二準直儀;以及第二檢測器,耦合到該第二準直儀,其中,該第二準直儀附接到該第二檢測器。
Applications Claiming Priority (6)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US202063071540P | 2020-08-28 | 2020-08-28 | |
US63/071,540 | 2020-08-28 | ||
US202163168778P | 2021-03-31 | 2021-03-31 | |
US63/168,778 | 2021-03-31 | ||
US17/399,768 | 2021-08-11 | ||
US17/399,768 US11723609B2 (en) | 2020-08-28 | 2021-08-11 | Spread field imaging collimators for radiation-based imaging and methods of using the same |
Publications (2)
Publication Number | Publication Date |
---|---|
TW202222258A TW202222258A (zh) | 2022-06-16 |
TWI816171B true TWI816171B (zh) | 2023-09-21 |
Family
ID=80353783
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
TW110130688A TWI816171B (zh) | 2020-08-28 | 2021-08-19 | 用於基於輻射的成像的擴散場成像準直儀 |
Country Status (6)
Country | Link |
---|---|
US (2) | US11723609B2 (zh) |
EP (1) | EP4203795A4 (zh) |
JP (1) | JP2023539256A (zh) |
CN (1) | CN115988994A (zh) |
TW (1) | TWI816171B (zh) |
WO (1) | WO2022046427A1 (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2024090412A1 (ja) * | 2022-10-25 | 2024-05-02 | 日本結晶光学株式会社 | X線撮像システムのx線検出器 |
WO2024090411A1 (ja) * | 2022-10-25 | 2024-05-02 | 日本結晶光学株式会社 | X線撮像システムのx線検出器 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150216488A1 (en) * | 2011-04-05 | 2015-08-06 | Triple Ring Technologies, Inc. | Method and apparatus for advanced x-ray imaging |
US20200146641A1 (en) * | 2018-11-09 | 2020-05-14 | Argospect Technologies Inc. | Collimator and Detector Based Medical Imaging Systems |
US20200261034A1 (en) * | 2019-02-18 | 2020-08-20 | Argospect Technologies Inc. | Collimators for medical imaging systems and image reconstruction methods thereof |
Family Cites Families (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4506374A (en) * | 1982-04-08 | 1985-03-19 | Technicare Corporation | Hybrid collimator |
US6271524B1 (en) * | 1998-08-05 | 2001-08-07 | Elgems, Ltd. | Gamma ray collimator |
US6353227B1 (en) * | 1998-12-18 | 2002-03-05 | Izzie Boxen | Dynamic collimators |
US6392235B1 (en) * | 1999-02-22 | 2002-05-21 | The Arizona Board Of Regents On Behalf Of The University Of Arizona | Coded-aperture system for planar imaging of volumetric sources |
US6580939B1 (en) * | 1999-11-04 | 2003-06-17 | Board Of Regents, The University Of Texas System | Method and system for reducing background artifacts from uniformly redundant array collimators in single photon emission computed tomography |
US6737652B2 (en) | 2000-09-29 | 2004-05-18 | Massachusetts Institute Of Technology | Coded aperture imaging |
US6652143B2 (en) * | 2001-04-12 | 2003-11-25 | Siemens Aktiengesellschaft | Method and apparatus for measuring the position, shape, size and intensity distribution of the effective focal spot of an x-ray tube |
US20030001098A1 (en) * | 2001-05-09 | 2003-01-02 | Stoddart Hugh A. | High resolution photon emission computed tomographic imaging tool |
US20040120464A1 (en) | 2002-12-19 | 2004-06-24 | Hoffman David Michael | Cast collimators for CT detectors and methods of making same |
US20060261278A1 (en) * | 2003-04-30 | 2006-11-23 | Roberto Accorsi | Soft x-ray imager with ten micrometer resolution |
US20050084072A1 (en) * | 2003-10-17 | 2005-04-21 | Jmp Industries, Inc., An Ohio Corporation | Collimator fabrication |
US7394053B2 (en) | 2004-09-09 | 2008-07-01 | Beth Israel Deaconess Medical Center, Inc. | Systems and methods for multi-modal imaging having a spatial relationship in three dimensions between first and second image data |
JP3928647B2 (ja) * | 2004-09-24 | 2007-06-13 | 株式会社日立製作所 | 放射線撮像装置およびそれを用いた核医学診断装置 |
US7683333B2 (en) * | 2004-10-15 | 2010-03-23 | Koninklijke Philips Electronics N.V. | Detector for nuclear medicine |
WO2007109227A2 (en) * | 2006-03-17 | 2007-09-27 | The Trustees Of The University Of Pennsylvania | Slit-slat collimation |
US20080087829A1 (en) * | 2006-04-11 | 2008-04-17 | John Hoppin | Single-photon emission computed tomography (SPECT) using helical scanning with multiplexing multi-pinhole apertures |
EP2224850B1 (en) | 2006-05-23 | 2012-06-27 | University Of The Witwatersrand, Johannesburg | Coded aperture masks for radiation-based medical imaging |
US7692156B1 (en) * | 2006-08-23 | 2010-04-06 | Radiation Monitoring Devices, Inc. | Beam-oriented pixellated scintillators for radiation imaging |
US20080230707A1 (en) * | 2007-03-23 | 2008-09-25 | Verista Imaging, Inc. | High resolution near-field imaging method and apparatus |
US7671340B2 (en) * | 2007-07-16 | 2010-03-02 | General Electric Company | Adjustable-focal-length collimators method and system |
WO2011100044A2 (en) * | 2010-02-10 | 2011-08-18 | The Regents Of The University Of California, San Francisco | Spect targeted volume molecular imaging using multiple pinhole apertures |
EP2707750B1 (en) * | 2011-05-10 | 2016-01-06 | Universiteit Gent | Methods and systems for collimating |
US10488535B2 (en) | 2013-03-12 | 2019-11-26 | Rearden, Llc | Apparatus and method for capturing still images and video using diffraction coded imaging techniques |
US9182507B2 (en) * | 2013-12-20 | 2015-11-10 | General Electric Company | Imaging system using high and low energy collimation |
US9349495B2 (en) | 2014-04-22 | 2016-05-24 | General Electric Company | Systems and methods for improved collimation sensitivity |
-
2021
- 2021-08-11 US US17/399,768 patent/US11723609B2/en active Active
- 2021-08-12 WO PCT/US2021/045770 patent/WO2022046427A1/en unknown
- 2021-08-12 EP EP21862369.2A patent/EP4203795A4/en active Pending
- 2021-08-12 JP JP2023513548A patent/JP2023539256A/ja active Pending
- 2021-08-12 CN CN202180053274.5A patent/CN115988994A/zh active Pending
- 2021-08-19 TW TW110130688A patent/TWI816171B/zh active
-
2023
- 2023-08-14 US US18/449,189 patent/US20230389882A1/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150216488A1 (en) * | 2011-04-05 | 2015-08-06 | Triple Ring Technologies, Inc. | Method and apparatus for advanced x-ray imaging |
US20200146641A1 (en) * | 2018-11-09 | 2020-05-14 | Argospect Technologies Inc. | Collimator and Detector Based Medical Imaging Systems |
US20200261034A1 (en) * | 2019-02-18 | 2020-08-20 | Argospect Technologies Inc. | Collimators for medical imaging systems and image reconstruction methods thereof |
Also Published As
Publication number | Publication date |
---|---|
US20230389882A1 (en) | 2023-12-07 |
US20220061782A1 (en) | 2022-03-03 |
US11723609B2 (en) | 2023-08-15 |
EP4203795A1 (en) | 2023-07-05 |
WO2022046427A1 (en) | 2022-03-03 |
JP2023539256A (ja) | 2023-09-13 |
EP4203795A4 (en) | 2024-02-14 |
TW202222258A (zh) | 2022-06-16 |
CN115988994A (zh) | 2023-04-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10502844B2 (en) | Sparse acquisition gamma cameras | |
US5818050A (en) | Collimator-free photon tomography | |
US6762413B2 (en) | Correction for depth-dependent sensitivity in rotating slat-collimated gamma camera | |
US11213262B2 (en) | Collimator and detector based medical imaging systems | |
TWI816171B (zh) | 用於基於輻射的成像的擴散場成像準直儀 | |
US11501474B2 (en) | Collimators for medical imaging systems and image reconstruction methods thereof | |
JP2011503535A (ja) | 間接放射線検出器 | |
WO2011058891A1 (ja) | 放射線撮像装置およびそれを用いた核医学診断装置 | |
US20050285042A1 (en) | Nuclear imaging system using rotating scintillation bar detectors with slat collimation and method for imaging using the same | |
EP4237876A1 (en) | Radiation detection system and method | |
Li | Novel Gamma-Ray Imaging and Detection Methods | |
WO2023244268A1 (en) | Emission tomography with generalized time encoded aperture imaging | |
CN116917773A (zh) | 辐射检测系统和方法 | |
JP2012137349A (ja) | ピクセル型放射線撮像装置,画像作成方法,プラナー画像の奥行き位置推定方法 |