The Radon transform and its adjoint, the back-projection operator, can both be expressed as convolutions in log-polar coordinates. Hence, fast algorithms for the application of these operators can be constructed by us...
详细信息
The Radon transform and its adjoint, the back-projection operator, can both be expressed as convolutions in log-polar coordinates. Hence, fast algorithms for the application of these operators can be constructed by using the FFT, if data is resampled at log-polar coordinates. Radon data is typically measured on an equally spaced grid in polar coordinates, and reconstructions are represented (as images) in Cartesian coordinates. Therefore, in addition to FFT, several steps of interpolation have to be conducted in order to apply the Radon transform and the back-projection operator by means of convolutions. However, in comparison to the interpolation conducted in Fourier-based gridding methods, the interpolation performed in the Radon and image domains will typically deal with functions that are substantially less oscillatory. Reasonable reconstruction results can thus be expected using interpolation schemes of moderate order. The approach also provides better control over the artifacts that can appear due to measurement errors. Both the interpolation and the FFT operations can be efficiently implemented on graphical processor units (GPUs). For the interpolation, it is possible to make use of the fact that linear interpolation is hard-wired on GPUs, meaning that it has the same computational cost as direct memory access. Cubic order interpolation schemes can be constructed by combining linear interpolation steps, and this provides important computation speedup. We provide details about how the Radon transform and the back-projection can be implemented efficiently as convolution operators on GPUs. For large data sizes, these algorithms are several times faster than those of other software packages based on GPU implementations of the Radon transform and the back-projection operator. Moreover, the gain in computational speed is substantially higher when comparing against other CPU-based algorithms.
To explore the redundancy in large ESPRIT arrays, we propose fast algorithms for direction-of-arrival (DOA) finding without computing (partial) eigendecompositions. For batch processing, the algorithms have computatio...
详细信息
To explore the redundancy in large ESPRIT arrays, we propose fast algorithms for direction-of-arrival (DOA) finding without computing (partial) eigendecompositions. For batch processing, the algorithms have computational complexity O(M(2)d) for the covariance matrix case and O(MNd) for the data matrix case, where M is the number of sensors, d the number of signals and N the number of snapshots. For a real-time implementation of the proposed algorithms, the computational complexity is O(Md)+O(d(3)) per data vector update. We also present numerical simulation results to illustrate the efficiency of the proposed algorithms.
Efficient algorithms are developed for area morphology. As opposed to traditional morphological operations that alter grayscale images via a concatenation of order statistic filters, the area morphological operators m...
详细信息
Efficient algorithms are developed for area morphology. As opposed to traditional morphological operations that alter grayscale images via a concatenation of order statistic filters, the area morphological operators manipulate connected components within the image level sets. Essentially, the area morphology filters are capable of removing objects based on the object area solely. These operators can then be effectively used in important multiscale and scale space tasks such as object-based coding and hierarchical image searches. Unfortunately, the traditional implementation of these filters based on level set theory precludes real-time implementation. This paper reviews previous fast algorithms and introduces a pyramidal approach. The full pyramidal algorithm is over 1000 times faster than the standard algorithm for typical image sizes. The paper provides supporting simulation results in terms of computational complexity and solution quality. (C) 2001 Academic Press.
An algorithm is introduced for the rapid evaluation at appropriately chosen nodes on the two-dimensional sphere S-2 in R-3 of functions specified by their spherical harmonic expansions ( known as the inverse spherical...
详细信息
An algorithm is introduced for the rapid evaluation at appropriately chosen nodes on the two-dimensional sphere S-2 in R-3 of functions specified by their spherical harmonic expansions ( known as the inverse spherical harmonic transform), and for the evaluation of the coefficients in spherical harmonic expansions of functions specified by their values at appropriately chosen points on S-2 ( known as the forward spherical harmonic transform). The procedure is numerically stable and requires an amount of CPU time proportional to N-2( logN) log( 1/epsilon), where N-2 is the number of nodes in the discretization of S-2, and e is the precision of computations. The performance of the algorithm is illustrated via several numerical examples.
fast algorithms for generalized predictive control (GPC) are derived by adopting an approach whereby dynamic programming and a polynomial formulation are jointly exploited. They consist of a set of coupled linear poly...
详细信息
fast algorithms for generalized predictive control (GPC) are derived by adopting an approach whereby dynamic programming and a polynomial formulation are jointly exploited. They consist of a set of coupled linear polynomial recursions by which the dynamic output feedback GPC law is recursively computed with only O(Nn) computations for an n-th order plant and N-steps prediction horizon.
We consider the problem of estimating the uncertainty in large-scale linear statistical inverse problems with high-dimensional parameter spaces within the framework of Bayesian inference. When the noise and prior prob...
详细信息
We consider the problem of estimating the uncertainty in large-scale linear statistical inverse problems with high-dimensional parameter spaces within the framework of Bayesian inference. When the noise and prior probability densities are Gaussian, the solution to the inverse problem is also Gaussian and is thus characterized by the mean and covariance matrix of the posterior probability density. Unfortunately, explicitly computing the posterior covariance matrix requires as many forward solutions as there are parameters and is thus prohibitive when the forward problem is expensive and the parameter dimension is large. However, for many ill-posed inverse problems, the Hessian matrix of the data misfit term has a spectrum that collapses rapidly to zero. We present a fast method for computation of an approximation to the posterior covariance that exploits the low-rank structure of the preconditioned (by the prior covariance) Hessian of the data misfit. Analysis of an infinite-dimensional model convection-diffusion problem, and numerical experiments on large-scale three-dimensional convection-diffusion inverse problems with up to 1.5 million parameters, demonstrate that the number of forward PDE solves required for an accurate low-rank approximation is independent of the problem dimension. This permits scalable estimation of the uncertainty in large-scale ill-posed linear inverse problems at a small multiple (independent of the problem dimension) of the cost of solving the forward problem.
This paper presents fast algorithms for type-II, type-III, and type-IV generalized discrete Hartley transform. In particular, new odd-factor algorithms are derived to support transforms whose sequence length contains ...
详细信息
This paper presents fast algorithms for type-II, type-III, and type-IV generalized discrete Hartley transform. In particular, new odd-factor algorithms are derived to support transforms whose sequence length contains multiple odd factors. By jointly using the odd-factor and radix-2 algorithms, fast computation for arbitrarily composite sequence length can be achieved. Compared to other reported algorithms, the proposed ones have a regular computational structure, achieve a substantial reduction of computational complexity, and support a wider range of choices on the sequence length.
In this paper, algorithms for fast implementations of Montgomery's modular multiplication algorithm are proposed. These algorithms use nonredundant multibit recoding. Two techniques, one based on multiplication fo...
详细信息
In this paper, algorithms for fast implementations of Montgomery's modular multiplication algorithm are proposed. These algorithms use nonredundant multibit recoding. Two techniques, one based on multiplication followed by estimation of scaled residue and the other involving integrated multiplication with scaled residue estimation, have been considered in detail. Techniques for simplifying the computation to estimate the multiple of the modulus to be added have been described. The area and computation time requirements of the proposed techniques are estimated for a general radix in order to highlight the time-space trade-offs.
We reformulate the original component-by-component algorithm for rank-1 lattices in a matrix-vector notation so as to highlight its structural properties. For function spaces similar to a weighted Korobov space, we de...
详细信息
We reformulate the original component-by-component algorithm for rank-1 lattices in a matrix-vector notation so as to highlight its structural properties. For function spaces similar to a weighted Korobov space, we derive a technique which has construction cost O(sn log(n)), in contrast with the original algorithm which has construction cost O(sn(2)). Herein s is the number of dimensions and n the number of points ( taken prime). In contrast to other approaches to speed up construction, our fast algorithm computes exactly the same quantity as the original algorithm. The presented algorithm can also be used to construct randomly shifted lattice rules in weighted Sobolev spaces.
We survey algorithms for computing isogenies between elliptic curves defined over a field of characteristic either 0 or a large prime. We introduce a new algorithm that computes an isogeny of degree l (l different fro...
详细信息
We survey algorithms for computing isogenies between elliptic curves defined over a field of characteristic either 0 or a large prime. We introduce a new algorithm that computes an isogeny of degree l (l different from the characteristic) in time quasi-linear with respect to l. This is based in particular on fast algorithms for power series expansion of the Weierstrass P-function and related functions.
暂无评论