We study the statistical convergence properties of two Cast iterative reconstruction algorithms, the rescaled block-iterative (RBI) and ordered subset (OS) EM algorithms, in the context of cardiac SPECT with 3D detect...
详细信息
We study the statistical convergence properties of two Cast iterative reconstruction algorithms, the rescaled block-iterative (RBI) and ordered subset (OS) EM algorithms, in the context of cardiac SPECT with 3D detector response modeling. The Monte Carlo method was used to generate nearly noise-free projection data modeling the effects of attenuation, detector response, and scatter from the MCAT phantom. One thousand noise realizations were generated with an average count level approximating a typical Tl-201 cardiac study. Each noise realization was reconstructed using the RBI and OS algorithms for cases with and without detector response modeling. For each iteration up to twenty, we generated mean and variance images, as well as covariance images for six specific locations. Both OS and RBI converged in the mean to results that were close to the noise-free ML-EM result using the same projection model. When detector response was not modeled in the reconstruction, RBI exhibited considerably lower noise variance than OS for the same resolution. When 3D detector response was modeled, the RBI-EM provided a small improvement in the tradeoff between noise level and resolution recovery, primarily in the axial direction, while OS required about half the number of iterations of RBI to reach the same resolution. We conclude that OS is faster than RBI, but may be sensitive to errors in the projection model. Both OS-EM and RBI-EM are effective alternatives to the ML-EM algorithm, but noise level and speed of convergence depend on the projection model used.
In this paper we will present numerical studies of the performance of a 3D Compton camera being developed at the University of Michigan. We present a physical model of the camera which exploits symmetries by using an ...
详细信息
In this paper we will present numerical studies of the performance of a 3D Compton camera being developed at the University of Michigan. We present a physical model of the camera which exploits symmetries by using an adapted spatial sampling pattern in the object domain. This model increases the sparsity of the transition matrix to reduce the very high storage and computation requirements. This model also allows the decomposition of the transition matrix into several small blocks that are easy to store. We also present the transition matrix of our system analytically computed with a real time algorithm. We geometrically optimize our hemispherical sampling and develop a 3D volumetric interpolation. Finally, we present a 3D image reconstruction method which uses the Gauss-Seidel algorithm to minimize a penalized least square objective.
We conducted an extensive computational experiment, lasting multiple CPU-years, to optimally select parameters for two important classes of algorithms for finding sparse solutions of underdetermined systems of linear ...
详细信息
We conducted an extensive computational experiment, lasting multiple CPU-years, to optimally select parameters for two important classes of algorithms for finding sparse solutions of underdetermined systems of linear equations. We make the optimally tuned implementations available at ***;they run "out of the box" with no user tuning: it is not necessary to select thresholds or know the likely degree of sparsity. Our class of algorithms includes iterative hard and soft thresholding with or without relaxation, as well as CoSaMP, subspace pursuit and some natural extensions. As a result, our optimally tuned algorithms dominate such proposals. Our notion of optimality is defined in terms of phase transitions, i.e., we maximize the number of nonzeros at which the algorithm can successfully operate. We show that the phase transition is a well-defined quantity with our suite of random underdetermined linear systems. Our tuning gives the highest transition possible within each class of algorithms. We verify by extensive computation the robustness of our recommendations to the amplitude distribution of the nonzero coefficients as well as the matrix ensemble defining the underdetermined system. Our findings include the following. 1) For all algorithms, the worst amplitude distribution for nonzeros is generally the constant-amplitude random-sign distribution, where all nonzeros are the same amplitude. 2) Various random matrix ensembles give the same phase transitions;random partial isometries may give different transitions and require different tuning. 3) Optimally tuned subspace pursuit dominates optimally tuned CoSaMP, particularly so when the system is almost square.
The authors introduce two algorithms in quadratic and sign gradient-type versions for spikey deconvolution based on a double (detection-estimation) formulation of the problem. Both algorithms are based on the applicat...
详细信息
The authors introduce two algorithms in quadratic and sign gradient-type versions for spikey deconvolution based on a double (detection-estimation) formulation of the problem. Both algorithms are based on the application of a selection operator to detect the nonzero values of a solution provided by a standard approach. In the first algorithm, a gradient-type iterative reduction of the error is used. The second algorithm minimizes the signal to be obtained. Their performances using l/sub 1/ and l/sub 2/ norms in the minimization are compared in terms of resolution, robustness, speed of convergence, and computational load.< >
Presented high-level synthesis describes scheduling for wide class of DSP algorithms. Several FPGA vendors or even ASIC designs are targeted via Handel-C compiled by Celoxica DK3.1 compiler. Using the authors' app...
详细信息
Presented high-level synthesis describes scheduling for wide class of DSP algorithms. Several FPGA vendors or even ASIC designs are targeted via Handel-C compiled by Celoxica DK3.1 compiler. Using the authors' approach, the designer can easily change type of used pipelined arithmetic modules and then check new performance. The optimal time schedule is found by cyclic scheduling using integer linear programming while minimizing the schedule period in the terms of clock cycles. Experimental results in HW implementation, performed on logarithmic arithmetic and floating-point arithmetic, confirm significant influence of the period on the resulting performance of DSP algorithms.
The Poisson model based image reconstruction algorithms in PET require the estimation of the mean value of random coincidences. The random data are often acquired with a delayed coincidence technique, and expected ran...
详细信息
The Poisson model based image reconstruction algorithms in PET require the estimation of the mean value of random coincidences. The random data are often acquired with a delayed coincidence technique, and expected randoms are estimated through variance reduction (VR) of measured delayed coincidences. In the past we have developed iterative VR algorithms that estimate singles rates from random compressed sinograms. The final estimation of singles rates can be noisy in short scans. In the present work, we apply commonly used and specific regularization methods. A sequential coordinate ascent algorithm is used to maximize the penalized Poisson Likelihood objective function. Measured data from a Siemens TruePoint clinical scanner are used to validate the algorithm performance.
In this paper we present an application of some iterative algorithms to solve a deconvolution problem. The aim is to solve an integral equation between an impulse response of a sensor and a desired and not accessible ...
详细信息
In this paper we present an application of some iterative algorithms to solve a deconvolution problem. The aim is to solve an integral equation between an impulse response of a sensor and a desired and not accessible input signal. Experimentally, we want to remove the distortion effects of the measuring system on an input signal. The applications concern the measurement of the distance between a flat eddy-current proximity sensor and a metallic surface in front of it. Numerical results are presented as an illustration. Experiments have been carried out to ascertain the performance of the conjugate gradient algorithm outlined in this paper. We also compare its performance to that of the iterative sequential algorithm of Placko et al. (C) 1997 Elsevier Science S.A.
algorithms and architectures for performing modular multiplication operations, which is central to crypto-system and authentication schemes, are important in today's needs of secure communications. This paper pres...
详细信息
algorithms and architectures for performing modular multiplication operations, which is central to crypto-system and authentication schemes, are important in today's needs of secure communications. This paper presents two new iterative algorithms for modular multiplication. The implementation of these algorithms yields to scalable architectures that can be used for any modulus without altering the design. In addition, the Radix-2 algorithm shows almost similar features when compared with similar architectures available in the literature. Furthermore, the radix-4 algorithm can be used to develop higher radix algorithms since it only requires the use of powers of two of the modulus.
We review the notion of factor graphs and give an overview of recent developments in the field. Factor graphs provide a framework in which iterative algorithms and iteratively decodable codes are easily described. The...
详细信息
We review the notion of factor graphs and give an overview of recent developments in the field. Factor graphs provide a framework in which iterative algorithms and iteratively decodable codes are easily described. They are obtained as a slight modification of the graphs, which were introduced by Wiberg, Loeliger and Koetter (1995). The notion of factor graphs is related to other graphical models, such as Bayesian networks and Markov random fields, and allows an easy integration of the decoding task and related problems in a communication system, such as channel estimation and equalization. Recently a number of research topics involving factor graphs have developed. These include the construction and representation of codes on graphs and an analysis of iterative decoding algorithm stemming from the theory of dynamical systems.
A class of very efficient iterative algorithms is presented for the analysis of planar scatterers at the very fine detail regime. The mechanisms causing slow convergence in these cases have been identified and their e...
详细信息
A class of very efficient iterative algorithms is presented for the analysis of planar scatterers at the very fine detail regime. The mechanisms causing slow convergence in these cases have been identified and their effect corrected in a newly introduced algorithm. Further convergence enhancement is achieved by one or two additional minimization steps. The resultant new algorithms converge very fast for problems with fine details sampled faster than 0.001/spl lambda/, where other methods fail.
暂无评论