Carroll Liam, Enger Shirin A
Medical Physics Unit, Department of Oncology, Faculty of Medicine, McGill University, Montréal, Quebec, Canada.
Lady Davis Institute for Medical Research, Jewish General Hospital, Montreal, Quebec, Canada.
Med Phys. 2023 Mar;50(3):1647-1659. doi: 10.1002/mp.16055. Epub 2022 Nov 9.
Dynamic positron emission tomography (dPET) is a nuclear medicine imaging technique providing functional images for organs of interest with applications in oncology, cardiology, and drug discovery. This technique requires the acquisition of the time-course arterial plasma activity concentration, called the arterial input function (AIF), which is conventionally acquired via arterial blood sampling.
The aim of this study was to (A) optimize the geometry for a novel and cost efficient non-invasive detector called NID designed to measure the AIF for dPET scans through Monte Carlo simulations and (B) develop a clinical data analysis chain to successfully separate the arterial component of a simulated AIF signal from the venous component.
The NID was optimized by using an in-house Geant4-based software package. The sensitive volume of the NID consists of a band of 10 cm long and 1 mm in diameter scintillating fibers placed over a wrist phantom. The phantom was simulated as a cylinder, 10 cm long and 6.413 cm in diameter comprised of polyethylene with two holes placed through it to simulate the patient's radial artery and vein. This phantom design was chosen to match the wrist phantom used in our previous proof of concept work. Two geometries were simulated with different arrangements of scintillating fibers. The first design used a single layer of 64 fibers. The second used two layers, an inner layer with 29 fibers and an outer layer with 30 fibers. Four positron emitting radioisotopes were simulated: F, C, O, and Ga with 100 million simulated decay events per run. The total and intrinsic efficiencies of both designs were calculated as well as the full width half maximum (FWHM) of the signal. In addition, contribution by the annihilation photons versus positrons to the signal was investigated. The results obtained from the two simulated detector models were compared. A clinical data analysis chain using an expectation maximization maximum likelihood algorithm was tested. This analysis chain will be used to separate arterial counts from the total signal.
The second NID design with two layers of scintillating fibers had a higher efficiency for all simulations with a maximum increase of 17% total efficiency for C simulation. All simulations had a significant annihilation photon contribution. The signal for F and C was almost entirely due to photons. The clinical data analysis chain was within 1% of the true value for 434 out of 440 trials. Further experimental studies to validate these simulations will be required.
The design of the NID was optimized and its efficiency increased through Monte Carlo simulations. A clinical data analysis chain was successfully developed to separate the arterial component of an AIF signal from the venous component. The simulations show that the NID can be used to accurately measure the AIF non-invasively for dPET scans.
动态正电子发射断层扫描(dPET)是一种核医学成像技术,可为感兴趣的器官提供功能图像,应用于肿瘤学、心脏病学和药物研发领域。该技术需要获取时程动脉血浆活度浓度,即动脉输入函数(AIF),传统上是通过动脉采血来获取。
本研究的目的是:(A)通过蒙特卡罗模拟优化一种新型且经济高效的无创探测器(称为NID)的几何结构,该探测器旨在测量dPET扫描的AIF;(B)开发一个临床数据分析链,以成功地将模拟AIF信号的动脉成分与静脉成分分离。
使用基于Geant4的内部软件包对NID进行优化。NID的灵敏体积由放置在手腕模型上的10厘米长、1毫米直径的闪烁光纤带组成。该模型被模拟为一个10厘米长、6.413厘米直径的圆柱体,由聚乙烯制成,有两个孔贯穿其中,以模拟患者的桡动脉和静脉。选择这种模型设计是为了匹配我们之前概念验证工作中使用的手腕模型。模拟了两种闪烁光纤排列不同的几何结构。第一种设计使用单层64根光纤。第二种使用两层,内层29根光纤,外层30根光纤。模拟了四种发射正电子的放射性同位素:F、C、O和Ga,每次运行模拟1亿次衰变事件。计算了两种设计的总效率和本征效率以及信号的半高宽(FWHM)。此外,研究了湮灭光子与正电子对信号的贡献。比较了从两个模拟探测器模型获得的结果。测试了使用期望最大化最大似然算法的临床数据分析链。该分析链将用于从总信号中分离动脉计数。
对于所有模拟,具有两层闪烁光纤的第二种NID设计效率更高,在C模拟中总效率最高提高了17%。所有模拟都有显著的湮灭光子贡献。F和C的信号几乎完全归因于光子。在440次试验中的434次试验中,临床数据分析链与真实值的偏差在1%以内。需要进一步的实验研究来验证这些模拟。
通过蒙特卡罗模拟对NID的设计进行了优化,并提高了其效率。成功开发了一个临床数据分析链,以将AIF信号的动脉成分与静脉成分分离。模拟表明,NID可用于非侵入性地准确测量dPET扫描的AIF。