Department of Radiation Oncology (Medical Physics), Nova Scotia Cancer Centre, Halifax, NS, B3H 4R2, Canada.
Department of Physics and Atmospheric Science, Dalhousie University, Halifax, NS, B3H 4J5, Canada.
Med Phys. 2019 Feb;46(2):528-543. doi: 10.1002/mp.13354. Epub 2019 Jan 16.
The purpose of this study was to develop a novel patient-specific pixel-based weighting factor dual-energy (PP-DE) algorithm to effectively suppress bone throughout the image and overcome the limitation of the conventional DE algorithm with constant weighting factor which is restricted to regions with uniform patient thickness. Additionally, to derive theoretical expressions to describe the dependence of the weighting factors on several imaging parameters and validate them with measurement.
A step phantom was constructed consisting of slabs of solid water and bone materials. Thicknesses of bone ranged [0-6] cm in one direction and solid water [5-30] cm in the other direction. Projection images at 60 and 140 kVp were acquired using a clinical imaging system. Optimal weighting factors were found by iteratively varying it in the range [0-1.4], where bone and soft-tissue contrast-to-noise ratio (CNR) reached zero. Bone and soft-tissue digitally reconstructed thicknesses were created using computed tomography (CT) images of a Rando phantom and ray tracing techniques. A weighting factor image (ω) was calculated using digitally reconstructed thicknesses (DRTs) and precalculated weighting factors from the step phantom. This ω image was then used to generate a PP-DE image. The PP-DE image was compared to the conventional DE image which uses a constant weighting factor throughout the image. The effect of the misaligned ω image on PP-DE images was investigated by acquiring LE and HE images at various shifts of Rando phantom. A rigid registration was used based on mutual information algorithm in Matlab. The signal-to-noise ratios (SNR) were calculated in the step phantom for the PP-DE image and compared to that of conventional DE technique. Analytical expressions for theoretical weighting factors were derived which included various effects such as beam hardening, scatter, and detector response. The analytical expressions were simulated in Spektr3.0 for different bone and solid water thicknesses as per the step phantom. A tray of steel pins was constructed and used with the step phantom to remove the scattered radiation. The simulated theoretical weighting factors were validated by comparing to those from the step phantom measurement.
Optimal weighting factor values for the step phantom varied from 0.633 to 1.372 depending on region thickness. Thicker regions required larger weighting factors for bone cancellation. The PP-DE image of the Rando phantom favorably cancelled both ribs and spine, whereas in the conventional DE image, only one could be cancelled at a time. The misaligned ω image was less effective in removing all bones indicating the importance of alignment as part of the PP-DE algorithm implementation. The SNRs for the PP-DE image was larger than those of the conventional DE images for regions which required smaller weighting factors for bone suppression. Comparisons of measured and simulated weighting factors demonstrated a 3% agreement for all bone overlapped regions except for the thickest region with 30 cm of solid water overlapped with 6 cm bone where the signal was lost due to excess attenuation.
A novel PP-DE algorithm was developed which can create higher quality DE images with enhanced bone cancellation and improved noise characteristics compared to conventional DE technique. In addition, theoretical weighting factor expressions were derived and validated against measurement.
本研究旨在开发一种新的基于像素的加权因子双能(PP-DE)算法,以有效抑制整个图像中的骨骼,并克服传统 DE 算法中恒定加权因子的局限性,该算法仅限于患者厚度均匀的区域。此外,推导出描述加权因子依赖于几个成像参数的理论表达式,并通过测量进行验证。
使用临床成像系统,构建了一个由实心水和骨材料板组成的阶梯式幻影。骨的厚度在一个方向上从 0 到 6cm,实心水的厚度在另一个方向上从 5 到 30cm。在 60 和 140kVp 下获取投影图像。通过在 0 到 1.4 的范围内迭代地改变它来找到最佳的加权因子,在该范围内,骨和软组织的对比度噪声比(CNR)达到零。使用 Rando 幻影的 CT 图像和射线追踪技术创建骨和软组织的数字重建厚度。使用数字重建厚度(DRTs)和从台阶幻影中预先计算的加权因子计算加权因子图像(ω)。然后,使用该 ω 图像生成 PP-DE 图像。将 PP-DE 图像与传统 DE 图像进行比较,后者在整个图像中使用恒定的加权因子。通过在 Rando 幻影的各种移位处获取 LE 和 HE 图像来研究错位 ω 图像对 PP-DE 图像的影响。在 Matlab 中,使用基于互信息算法的刚性配准。在台阶幻影中计算了 PP-DE 图像的信噪比(SNR),并与传统 DE 技术进行了比较。推导了包括束硬化、散射和探测器响应等各种影响的理论加权因子的解析表达式。在 Spektr3.0 中,根据台阶幻影中的不同骨和实心水厚度对理论加权因子进行了模拟。构建了一个钢针托盘并与台阶幻影一起使用,以去除散射辐射。通过将模拟的理论加权因子与台阶幻影的测量值进行比较,验证了其有效性。
台阶幻影的最佳加权因子值因区域厚度而异,范围从 0.633 到 1.372。较厚的区域需要更大的加权因子来消除骨。Rando 幻影的 PP-DE 图像能够很好地消除肋骨和脊柱,而在传统的 DE 图像中,一次只能消除一个。错位的 ω 图像在去除所有骨骼方面效果较差,表明作为 PP-DE 算法实施一部分的对准的重要性。对于需要较小加权因子来抑制骨骼的区域,PP-DE 图像的 SNR 大于传统 DE 图像的 SNR。测量和模拟加权因子的比较表明,除了重叠有 30cm 实心水和 6cm 骨的最厚区域外,所有骨重叠区域的一致性为 3%,在该区域,由于过度衰减,信号丢失。
开发了一种新的 PP-DE 算法,与传统的 DE 技术相比,它可以创建更高质量的 DE 图像,具有更好的骨骼消除和改善的噪声特性。此外,推导出了理论加权因子表达式,并通过测量进行了验证。