Iterative restoration of SPECT projection images

Photon attenuation and the limited nonstationary spatial resolution of the detector can reduce both qualitative and quantitative image quality in single photon emission computed tomography (SPECT). In this paper, a reconstruction approach is described which can compensate for both of these degradations. The approach involves processing the projection data with Bellini's method for attenuation compensation followed by an iterative deconvolution technique which uses the frequency distance principle (FDP) to model the distance-dependent camera blur. Modeling of the camera blur with the FDP allows an efficient implementation using fast Fourier transform (FFT) methods. After processing of the projection data, reconstruction is performed using filtered backprojection. Simulation studies using two different brain phantoms show that this approach gives reconstructions with a favorable bias versus noise tradeoff, provides no visually undesirable noise artifacts, and requires a low computational load     

大型集装箱检测系统管理站的任务需求分析与结构规划

由加速器周期性产生的Ｘ射线经前准直器成形为扇形脉冲束流，透过集装箱后再经过后准直器进入阵列探测器，再由前置放大器放大成形为电压信号送入采样保持电路，经多路选择器送去变换成为一组数字量，再经过一系列的预处理就可形成一列象素。由于集装箱是载于拖动装置上的，随着拖动装置的运动，这一列列的象素就可组成一幅完整的集装箱透视图，这样的一个完整过程由中心管理站统一控制和同步。本文介绍实时扫描数据获取与处理，并给     

An Iterative Reconstruction from Truncated Projection Data

Various methods have been proposed for tomographic reconstruction from truncated projection data. In this paper, a reconstructive method is discussed which consists of iterations of filtered back-projection, reprojection and some nonlinear processings. First, the method is so constructed that it converges to a fixed point. Then, to examine its effectiveness, comparisons are made by computer experiments with two existing reconstructive methods for truncated projection data, that is, the method of extrapolation based on the smooth assumption followed by filtered back-projection, and modified additive ART.     

Limited Angular View Planar-Integral Projection Reconstruction in NMR Imaging

Limited angular view planar-integral projection reconstruction for fast multislice imaging in NMR is proposed and its efficacy is studied. Alghough it is known that the planar-integral projection reconstruction (PPR) method has advantages in signal to noise ratio improvement in NMR imaging, necessary data acquisition time is too long to be of practical value. In this paper, a study is undertaken on the variations of the PPR method which would lead to a shorter data acquisition time than the orginal PPR method. Also a new method which provides high quality multislice images through the use of limited angular view data acquisition mode is introduced. Computer simulation results and experimentally obtained images are presented for the demonstration of the potential of the method for clinical applications.     

Dynamic fragility concepts for equipment design and qualification

Fragility concepts are explored for use in the design and qualification of nuclear plant equipment and for relating the ultimate capability of equipment to that of the overall plant. In the most general sense, the fragility level of a device may depend on several different types of environmental stress or challenge factors (i.e., heat, nuclear radiation, vibration, etc.) that influence its operation. However, emphasis is concentrated on the dynamic and particularly the seismic fragility levels of equipment. A general definition of dynamic fragility and various methods for its measurement are described. The state of published data on nuclear equipment fragility is discussed, and limitations on its use are noted. From there, the concept of a standardized seismic fragility data base and its potential uses are considered. Various gaps in the methodology are identified, and recommendations for further research are suggested.     

基于SART算法的CL硬化伪影校正方法研究

X射线薄板层析成像(CL)系统是一种针对板状构件的新型断层成像系统,该系统中存在射束硬化现象。本工作研究了一种基于SART算法的硬化伪影校正算法（SART-BHC算法）,通过对该算法的研究及对CL系统矩阵的计算,使该方法能应用于CL中。通过模拟与实验相结合的方式,利用蒙特卡罗方法获得CL投影的模拟数据,在实际CL中获得实验数据,然后利用SART-BHC与FPB算法对两种数据分别重建。结果表明,SART-BHC算法可应用于CL中,不需要任何先验知识,能很好地校正硬化伪影,并对层间混叠有抑制作用。     

基于阵列式吸收发光CT法的辐射场实时测量系统的研究

分析了已有的各种核辐射剂量场测量方法的基础上, 利用闪烁体吸收发光特点并结合计算机实时成像处理技术,提出了对剂量场分布进行实时成像测量的新方法——阵列式吸收发光ＣＴ 法。研制出闪烁光纤阵列构成的探测器及其伺服控制系统, 用高灵敏度的电荷耦合器件（ＣＣＤ）拾取探测器产生的微弱闪烁光信号,并采用定点采集的方法对视频信号进行数据的快速采集。在图像重建方法上, 提出了迭代滤波反投影重建方法和利用非完全投影进行数据修复,对获取的投影数据进行变换和处理,并通过选择滤波函数及其适当的参数, 可获得最佳滤波效果,以重建剂量场的二维场分布,从而建立了剂量场的实时、高精度的成像测量系统。     

Thermal creep modeling of HT9 steel for fast reactor applications

The ferritic martensitic steel HT9 is a primary candidate material for the fuel cladding of liquid-metal-cooled fast reactors (LMFRs) owing to its excellent stability under irradiation. Thermal creep of fuel cladding is a potential life-limiting factor in the long-life fuel design of LMFRs. Using the measured data available in the literature, such as creep strain data, stress rupture data, and steady-state creep rate data, a generalized creep correlation was developed. The new correlation is based on the Garofalo equation and the modified Monkman-Grant equation, and it shows better agreement with the experimental data than existing correlations that use either the theta projection method or the minimum commitment method, making it more appropriate for use in long-life applications.     

离散Radon变换卷积分配性及在从投影直接提取CT图像边缘的应用

研究了离散Radon变换(DRT)在二维卷积运算方面的分配性。通过对这一性质离散条件下的计算方法的研究,提出了在CT重建过程之前直接对投影数据进行处理的方法,并对图像边缘提取这一具体工作进行了应用研究．给出了计算机模拟计算的结果。通过对传统图像边缘提取与从投影进行图像边缘提取这两种方式进行比较,分析了二者的运算时间的影响参数,论证了该方法在运算时间上的优越性,以及在CT图像处理中的应用前景。     

CT image correction for beam hardening using simulated projectiondata

A general beam-hardening correction technique is presented. Knowledge of the CT (computerized tomography) scanner X-ray spectrum is necessary. This postreconstruction method does not require the original projection data. Simulated projections through an uncorrected reconstructed image are used to correct for beam hardening. Errors in the mean linear attenuation coefficient are decreased from 30% to 5% with virtual elimination of the visual streaking artifact. The theoretical image improvement equals that of projection linearization postreconstruction methods using the original projection data. The correction is limited to cases where the material causing the beam hardening is contained within the reconstruction space     

基于硬件实现的锥束CT图像重建系统的存储机制设计

为了实现三维锥束CT图像重建加速系统的小型化,建立了基于FPGA的三维图像重建系统。并对该系统中所采用的FDK重建算法所需的数据存储量和数据传输量以及由SDRAM、SRAM和FPGA内部缓存组成的存储系统的数据吞吐率进行了研究。首先,根据FDK算法的滤波与反投影两个步骤介绍了三维锥束CT图像重建系统的数据规模。接着,介绍了一种以SDRAM为外部主存,以SRAM为外部缓存和以FPGA内部SRAM资源作为内部高速缓存的存储机制。然后,介绍了该存储机制的实现方法以及测试方法。最后对该三维图像重建系统的数据吞吐能力进行了测试,并将之与FDK算法所需的数据传输量进行了对比分析。试验结果表明:该存储机制的数据连续存取速度为151.9MB/s,数据随机存取速度为100MB/s,基本满足小规模的三维图像重建的数据存储与传输带宽的要求。     

Characterizing the MTF in 3D for a Quantized SPECT Camera Having Arbitrary Trajectories

The emergence of application-specific 3D tomographic small animal and dedicated breast imaging systems has stimulated the development of simple methods to quantify the spatial resolution or Modulation Transfer Function (MTF) of the system in three dimensions. Locally determined MTFs, obtained from line source measurements at specific locations, can characterize spatial variations in the system resolution and can help correct for such variations. In this study, a method is described to measure the MTF in 3D for a compact SPECT system that uses a 16 × 20 cm(2) CZT-based compact gamma camera and 3D positioning gantry capable of moving in different trajectories. Image data are acquired for a novel phantom consisting of three radioactivity-filled capillary tubes, positioned nearly orthogonally to each other. These images provide simultaneous measurements of the local MTF along three dimensions of the reconstructed imaged volume. The usefulness of this approach is shown by characterizing the MTF at different locations in the reconstructed imaged 3D volume using various (1) energy windows; (2) iterative reconstruction parameters including number of iterations, voxel size, and number of projection views; (3) simple and complex 3D orbital trajectories including simple vertical axis of rotation, simple tilt, complex circle-plus-arc, and complex sinusoids projected onto a hemisphere; and (4) object shapes in the camera's field of view. Results indicate that the method using the novel phantom can provide information on spatial resolution effects caused by system design, sampling, energy windows, reconstruction parameters, novel 3D orbital trajectories, and object shapes. Based on these measurements that are useful for dedicated tomographic breast imaging, it was shown that there were small variations in the MTF in 3D for various energy windows and reconstruction parameters. However, complex trajectories that uniformly sample the breast volume of interest were quantitatively shown to have slightly better spatial resolution performance than more simple orbits.     

Derivative-Hilbert-Backprojection based image reconstruction from truncated projections in helical cone-beam CT

In helical cone-beam computed tomography(CT), Feldkamp-Davis-Kress(FDK) based image reconstruction algorithms are by far the most popular. However, artifacts are commonly met in the presence of lateral projection truncation. The reason is that the ramp filter is global. To restrain the truncation artifacts, an approximate reconstruction formula is proposed based on the Derivative-Hilbert-Backprojection(DHB) framework. In the method, the first order derivative filter is followed by the Hilbert transform. Since the filtered projection values are almost zero by the first order derivative filter, the following Hilbert transform has little influence on the projection values, even though the projections are laterally truncated. The proposed method has two main advantages. First, it has comparable computational efficiency and image quality as well as the conventional helical FDK algorithm for non-truncated projections. The second advantage is that images can be reconstructed with acceptable quality and much lower computational cost in comparison to the Laplace operator based algorithm in cases with truncated projections. To point out the advantages of our method, simulations on the computer and real data experiments on our laboratory industrial cone-beam CT are conducted. The simulated and experimental results demonstrate that the method is feasible for image reconstruction in the case of projection truncation.     

基于平板探测器投影数据预处理的三维CT重构算法研究

大面积非晶硅平板探测器(flat panel detector-FPD)在工业透射射线成像领域得到越来越多的应用。基于FPD的圆轨道FDK型三维计算机层析成像技术(Computed Tomography,3D-CT)检测速度快,但成像质量不及2D-CT。为发挥3D-CT速度优势,抑制噪声同时提高图像重建的质量。以由平板探测器得到的二维投影数据为研究对象,提出了一种对投影数据分别进行高频和低频两次滤波的投影预处理方法,然后将分别重建的图像叠加最终得到高质量的重建图像。由算法的仿真实验结果表明,重建图像质量得到明显改善,系统噪声得到抑制。     

Study of R.f., Gradient Pulse and Magnet Instability Effect in NMR Tomography

Newly emerged NMR imaging requires careful studies on r. f. pulse shapes and sequencing for the selection of the region, gradient pulsing for the 2- or 3-D spatial coding, and suitable signal handling technique for the compensation of the inherent instability of the system, especially fluctuation of the static magnetic field. Above subjects are discussed in detail and a new method which would be useful for the line integral projection reconstruction is proposed. The method could equally be applied to other 3-D NMR imaging techniques such as KWE (Kumar-Welti-Ernst) direct Fourier reconstruction or planar integral projection type reconstruction.     

扇形束投影域重建—再投影迭代法

本文所提出的投影域重建——再投影迭代算法可以恢复不完整的扇形束投影数据,从而实现从有限的投影数据重建图象。这种算法把滤波——反投影图象重建过程和再投影的过程纳入一个公式,仅在投影域上估算未知的投影数据,因此避免了重建过程和反投影过程中的插值运算,改善了收敛性,提高了计算精度,缩短了运算时间。模拟实验表明:即使有70％的投影数据丢失,经过两次迭代就能以令人满意的精度恢复丢失的投影数据。     

Reconstruction from Limited Angular Projection Data Using Constrained Optimization

In this paper, a new type of pixel decomposition is introduced to suppress modeling noise and then an algebraic method is developed for reconstruction from limited angular projection data. The basic idea with this method is to directly control the tradeoff between resolution and noise, taking into account the uncertainty of the discrepancy criterion generally used in an algebraic method for reconstruction. The results of a simulation using a simple mathematical phantom are also presented, demonstrating the effectiveness of this method.     

Strain distribution around a crack tip in high temperature environment using picture processing

A new optical strain measuring method developed by the authors is described. Here, we employ a simple method which uses point recognition picture processing and the least square approximation method.The surface of the specimen for which we are going to obtain the strain distribution is marked with points using a punch or a paint marker beforehand. At each loading stage, the digital data of the black and white picture of the specimen's marked surface are taken with the television camera under the control of the computer to which the data are sent. From these data, the locations of the marked points at each loading stage are obtained by the picture processing technique. Using these locations decided, we obtain the least square approximation of the displacement distribution, and then the strain distribution.To show the effectiveness of the present technique, the method is applied to the strain measurement near the crack tip in the CT specimen of Type 304 stainless steel. Because of the high ductility of this material, the crack tip blunting is notable so that the crack tip opening displacement observed by this method is as large as 2.5 mm at the specimen surface before initiation.Next, the strain distributions of a center-notched panel made of Type 304 stainless steel under high temperature creep condition are obtained.It is found from the above two experiments that the present strain measuring technique is powerful for the measurement of high strain, say more than 1%, whereas it is unsuitable to apply to the case of small strain at the present developing stage.     

用薄膜~(252)Cf源测定ST451快中子探测器的相对效率

本文叙述了用飞行时间测量~(262)Cf瞬发裂变中子能谱来刻度ST 451探测器的快中子相对效率响应的方法。在飞行时间谱仪中,一个流气式的微型电离室作为裂变碎块的探测器。文中给出了对实验数据的处理和修正的过程。得到了阈值分别为0.420、0.625、0.825、1.168、1.565和1.882MeV中子能量从几百keV到10MeV范围内的相对探测效率。实验结果与用Monte Carlo方法计算的效率作了比较。     

基于半扇束扫描的CT图像迭代重建算法

为解决探测器尺寸受限导致的CT图像重建区域面积受限以及X射线辐射剂量较大问题,提出一种半扇束扫描模式下的CT图像迭代重建算法。该算法将探测器水平偏置,稀疏角度下采集物体投影数据,利用ART-TV迭代算法对投影数据进行重建,达到预设迭代次数后,输出重建图像。实验结果表明:半扫描可扩大重建面积,降低探测器成本;利用迭代重建算法,降低剂量的同时重建优质CT图像。