ICT并行重建算法设计与实现

分析了三维ICT近似重建FDK快速算法以及算法串行执行与并行执行的复杂度,利用MPI并行编程环境,对原有ICT重建流程进行并行化,在PC集群上实现了并行ICT重建,进行了系统集成。给出了系统应用实现的解决方案与软件集成流程图。实验结果证明串行重建和并行重建结果是一致的,并行重建可以得到比较理想地重建时间结果和比较理想的加速比与效率。     

Monte Carlo simulation of glandular dose in a dedicated breast CT system

A dedicated breast CT system (DBCT) is a new method for breast cancer detection proposed in recent years. In this paper, the glandular dose in the DBCT is simulated using the Monte Carlo method. The phantom shape is half ellipsoid, and a series of phantoms with different sizes, shapes and compositions were constructed. In order to optimize the spectra, monoenergy X-ray beams of 5-80 keV were used in simulation. The dose distribution of a breast phantom was studied: a higher energy beam generated more uniform distribution, and the outer parts got more dose than the inner parts. For polyenergtic spectra, four spectra of Al filters with different thicknesses were simulated, and the polyenergtic glandular dose was calculated as a spectral weighted combination of the monoenergetic dose.     

Beam hardening correction for a cone-beam CT system and its effect on spatial resolution

In this paper, we present a beam hardening correction (BHC) method in three-dimension space for a cone-beam computed tomography (CBCT) system in a mono-material case and investigate its effect on the spatial resolution. Due to the polychromatic character of the X-ray spectrum used, cupping and streak artifacts called beam hardening artifacts arise in the reconstructed CT images, causing reduced image quality. In addition, enhanced edges are introduced in the reconstructed CT images because of the beam hardening effect. The spatial resolution of the CBCT system is calculated from the edge response function (ERF) on different planes in space. Thus, in the CT images with beam hardening artifacts, enhanced ERFs will be extracted to calculate the modulation transfer function (MTF), obtaining a better spatial resolution that deviates from the real value. Reasonable spatial resolution can be obtained after reducing the artifacts. The 10% MTF value and the full width at half maximum (FWHM) of the point spread function with and without BHC are presented.     

一种新的二维CT扫描方式及其重建算法

介绍了一种二维CT扫描方式及其重建算法。它是一种利用近轴X光线,针对大型回转试件的二维CT检测方法。被扫描试件依次绕与探测器平行的多个旋转中心各旋转一周。在旋转过程中,射线源在周向均布的多个角位置发射射线,同时探测器接收射线并记录投影数据。同传统的二代CT相比,对一个具体例子的计算表明,它可以节约3/4的扫描时间。计算机仿真表明了这种扫描方式和重建算法的可行性。     

RT二维CT扫描的精确重建

基于多个旋转中心的旋转/平移扫描方法(RT扫描)可利用近轴X射线对大型回转试件进行二维CT检测。同传统的二代CT平移/旋转扫描(TR)相比,RT二维CT扫描的扫描时间是TR二维CT扫描时间的1/3。但由RT二维CT扫描插值全局重建的图像,在物体内部,其特征边缘有明显的伪影。采用一种精确的分区域重建方法可消除这种伪影。计算机仿真验证了这种重建算法的可行性。     

Measurement of spatial resolution of the micro-CT system

The performance test is an important and necessary work for the micro-CT （computed tomography） system. The focal spot size of the micro focus X-ray tube is measured. The method of measuring the spatial resolution of micro-CT is introduced. A line-pair resolution of 28.2 lp/mm at the 10% modulation transfer function （MTF） level can be achieved with 14.7 μm spot size, 12.3 μm voxel size and a 25 mm field of view. In addition, a tungsten wire with the diameter of 5 μm can be detected by the system.     

Computed laminography and reconstruction algorithm

Computed laminography (CL) is an alternative to computed tomography if large objects are to be inspected with high resolution. This is especially true for planar objects. In this paper, we set up a new scanning geometry for CL, and study the algebraic reconstruction technique (ART) for CL imaging. We compare the results of ART with variant weighted functions by computer simulation with a digital phantom. It proves that ART algorithm is a good choice for the CL system.     

Measurement of the spatial resolution and the relative density resolution in an industrial cone-beam micro computed tomography system

The spatial resolution and the relative density resolution are the two most critical indicators in CT system. The method recommended in the ASTM E1695-95 and GJB 5311-2004 is only suitable to the fan-beam CT system. In this paper, for industrial cone-beam micro CT system, we will adopt the edge response function (ERF) created by the step edges of a steel ball to measure the system 3D PSF and MTF. To describe the contrast discrimination function more accurately, we will first propose to extend the two-dimensional measurement region to the three-dimensional space. Our experimental spatial resolution is (55.56±0.56) lp/mm and the relative density resolution is 1% within 300 μm×300 μm×300 μm according to the 3σ rule.     

A generalized method of converting CT image to PET linear attenuation coefficient distribution in PET/CT imaging

The accuracy of attenuation correction in positron emission tomography scanners depends mainly on deriving the reliable 511-keV linear attenuation coefficient distribution in the scanned objects. In the PET/CT system, the linear attenu- ation distribution is usually obtained from the intensities of the CT image. However, the intensities of the CT image relate to the attenuation of photons in an energy range of 40 keV-140 keV. Before implementing PET attenuation correction, the intensities of CT images must be transformed into the PET 511-keV linear attenuation coefficients. However, the CT scan parameters can affect the effective energy of CT X-ray photons and thus affect the intensities of the CT image. Therefore, for PET/CT attenuation correction, it is crucial to determine the conversion curve with a given set of CT scan parameters and convert the CT image into a PET linear attenuation coefficient distribution. A generalized method is proposed for con- verting a CT image into a PET linear attenuation coefficient distribution. Instead of some parameter-dependent phantom calibration experiments, the conversion curve is calculated directly by employing the consistency conditions to yield the most consistent attenuation map with the measured PET data. The method is evaluated with phantom experiments and small animal experiments. In phantom studies, the estimated conversion curve fits the true attenuation coefficients accurately, and accurate PET attenuation maps are obtained by the estimated conversion curves and provide nearly the same correction results as the true attenuation map. In small animal studies, a more complicated attenuation distribution of the mouse is obtained successfully to remove the attenuation artifact and improve the PET image contrast efficiently.