Tensor representation of color images and fast 2D quaternion discrete Fourier transform
NASA Astrophysics Data System (ADS)
Grigoryan, Artyom M.; Agaian, Sos S.
2015-03-01
In this paper, a general, efficient, split algorithm to compute the two-dimensional quaternion discrete Fourier transform (2-D QDFT), by using the special partitioning in the frequency domain, is introduced. The partition determines an effective transformation, or color image representation in the form of 1-D quaternion signals which allow for splitting the N × M-point 2-D QDFT into a set of 1-D QDFTs. Comparative estimates revealing the efficiency of the proposed algorithms with respect to the known ones are given. In particular, a proposed method of calculating the 2r × 2r -point 2-D QDFT uses 18N2 less multiplications than the well-known column-row method and method of calculation based on the symplectic decomposition. The proposed algorithm is simple to apply and design, which makes it very practical in color image processing in the frequency domain.
Real-time 2D floating-point fast Fourier transforms for seeker simulation
NASA Astrophysics Data System (ADS)
Chamberlain, Richard; Lord, Eric; Shand, David J.
2002-07-01
The floating point Fast Fourier Transform (FFT) is one of the most useful basic functions available to the image and signal processing engineer allowing many complex and detailed special functions to be implemented more simply in the frequency domain. In the Hardware-in-the-Loop field an image transformed using FFT would allow the designer to think about accurate frequency based simulation of seeker lens effects, motion blur, detector transfer functions and much more. Unfortunately, the transform requires many hundreds of thousands or millions of floating point operations on a single modest sized image making it impractical for realtime Hardware-in-the-Loop systems. .until now. This paper outlines the development, by Nallatech, of an FPGA based IEEE floating point core. It traces the subsequent use of this core to develop a full 256 X 256 FFT and filter process implemented on COTS hardware at frame rates up to 150Hz. This transform can be demonstrated to model optical transfer functions at a far greater accuracy than the current spatial models. Other applications and extensions of this technique will be discussed such as filtering for image tracking algorithms and in the simulation of radar processing in the frequency domain.
2D Regimes of Non-Fourier Convection
NASA Astrophysics Data System (ADS)
Papanicolaou, N. C.
2010-11-01
In this work, we investigate the 2D flow in a rectangular cavity subject to both vertical and horizontal temperature gradients. The linearized model is studied and the effect of thermal relaxation, as described by the Maxwell-Cattaneo law of heat conduction is examined. To this end, a spectral numerical model is created based on a Galerkin expansion. The basis is the Cartesian product of systems of beam functions and trigonometric functions. The natural modes of the system are derived for both the Fourier and non-Fourier models. The results are compared to earlier works for the plain Fourier law. Our computations show that for the same set of parameters, the Maxwell-Cattaneo law yields modes which are quantitatively different from the Fourier. It is found that the real parts of the eigenvalues increase with the Straughan number Sg, which quantifies the non-Fourier effects. This confirms the destabilizing effect of the MC-law on the convective flow.
A novel sliding window algorithm for 2D discrete Fourier transform
NASA Astrophysics Data System (ADS)
Dong, Zhifang; Wu, Jiasong; Gui, Jiyong
2015-12-01
Discrete Fourier transform (DFT) is one of the most wildly used tools for signal processing. In this paper, a novel sliding window algorithm is presented for fast computing 2D DFT when sliding window shifts more than one-point. The propose algorithm computing the DFT of the current window using that of the previous window. For fast computation, we take advantage of the recursive process of 2D SDFT and butterfly-based algorithm. So it can be directly applied to 2D signal processing. The theoretical analysis shows that the computational complexity is equal to 2D SDFT when one sample comes into current window. As well, the number of additions and multiplications of our proposed algorithm are less than those of 2D vector radix FFT when sliding window shifts mutiple-point.
FFTLog: Fast Fourier or Hankel transform
NASA Astrophysics Data System (ADS)
Hamilton, Andrew J. S.
2015-12-01
FFTLog is a set of Fortran subroutines that compute the fast Fourier or Hankel (= Fourier-Bessel) transform of a periodic sequence of logarithmically spaced points. FFTLog can be regarded as a natural analogue to the standard Fast Fourier Transform (FFT), in the sense that, just as the normal FFT gives the exact (to machine precision) Fourier transform of a linearly spaced periodic sequence, so also FFTLog gives the exact Fourier or Hankel transform, of arbitrary order m, of a logarithmically spaced periodic sequence.
Fast Fourier Transform algorithm design and tradeoffs
NASA Technical Reports Server (NTRS)
Kamin, Ray A., III; Adams, George B., III
1988-01-01
The Fast Fourier Transform (FFT) is a mainstay of certain numerical techniques for solving fluid dynamics problems. The Connection Machine CM-2 is the target for an investigation into the design of multidimensional Single Instruction Stream/Multiple Data (SIMD) parallel FFT algorithms for high performance. Critical algorithm design issues are discussed, necessary machine performance measurements are identified and made, and the performance of the developed FFT programs are measured. Fast Fourier Transform programs are compared to the currently best Cray-2 FFT program.
Fast Fourier Transforms of Piecewise Constant Functions
NASA Astrophysics Data System (ADS)
Sorets, Eugene
1995-02-01
We present an algorithm for the evaluation of the Fourier transform of piecewise constant functions of two variables. The algorithm overcomes the accuracy problems associated with computing the Fourier transform of discontinuous functions; in fact, its time complexity is O (N2 logN + NP log2 (1/ε) + V log3 (1/ε)), where ε is the accuracy, N is the size of the problem, P is the perimeter of the set of discontinuities, and V is its number of vertices. The algorithm is based on the Lagrange interpolation formula and the Green's theorem, which are used to preprocess the data before applying the fast Fourier transform. It readily generalizes to higher dimensions and to piecewise smooth functions.
Fast Computation of Wideband Beam Pattern for Designing Large-Scale 2-D Arrays.
Chi, Cheng; Li, Zhaohui
2016-06-01
For real-time and high-resolution 3-D ultrasound imaging, the design of sparse distribution and weights of elements of a large-scale wideband 2-D array is needed to reduce hardware cost and achieve better directivity. However, due to the high time consumption of computing the wideband beam pattern, the design methods that need massive iterations have rarely been applied to design large-scale wideband 2-D arrays by directly computing the wideband beam pattern. In this paper, a fast method is proposed to realize the computation of a wideband beam pattern of arbitrary 2-D arrays in the far field in order to design large-scale wideband 2-D arrays. The proposed fast method exploits two important techniques: 1) nonuniform fast Fourier transform (FFT) and 2) short inverse FFT. Compared with the commonly used ultrasound simulator Field II, two orders of magnitude improvement in computation speed is achieved with comparable accuracy. The proposed fast method enables massive iterations of direct wideband beam pattern computation of arbitrary large-scale 2-D arrays. A design example in this paper demonstrates that the proposed fast method can help achieve better performance in designing large-scale wideband 2-D arrays. PMID:27046870
Building a symbolic computer algebra toolbox to compute 2D Fourier transforms in polar coordinates.
Dovlo, Edem; Baddour, Natalie
2015-01-01
The development of a symbolic computer algebra toolbox for the computation of two dimensional (2D) Fourier transforms in polar coordinates is presented. Multidimensional Fourier transforms are widely used in image processing, tomographic reconstructions and in fact any application that requires a multidimensional convolution. By examining a function in the frequency domain, additional information and insights may be obtained. The advantages of our method include: •The implementation of the 2D Fourier transform in polar coordinates within the toolbox via the combination of two significantly simpler transforms.•The modular approach along with the idea of lookup tables implemented help avoid the issue of indeterminate results which may occur when attempting to directly evaluate the transform.•The concept also helps prevent unnecessary computation of already known transforms thereby saving memory and processing time. PMID:26150988
2D Fourier series representation of gravitational functionals in spherical coordinates
NASA Astrophysics Data System (ADS)
Ghobadi-Far, Khosro; Sharifi, Mohammad Ali; Sneeuw, Nico
2016-05-01
2D Fourier series representation of a scalar field like gravitational potential is conventionally derived by making use of the Fourier series of the Legendre functions in the spherical harmonic representation. This representation has been employed so far only in the case of a scalar field or the functionals that are related to it through a radial derivative. This paper provides a unified scheme to represent any gravitational functional in terms of spherical coordinates using a 2D Fourier series representation. The 2D Fourier series representation for each individual point is derived by transforming the spherical harmonics from the geocentric Earth-fixed frame to a rotated frame so that its equator coincides with the local meridian plane of that point. In the obtained formulation, each functional is linked to the potential in the spectral domain using a spectral transfer. We provide the spectral transfers of the first-, second- and third-order gradients of the gravitational potential in the local north-oriented reference frame and also those of some functionals of frequent use in the physical geodesy. The obtained representation is verified numerically. Moreover, spherical harmonic analysis of anisotropic functionals and contribution analysis of the third-order gradient tensor are provided as two numerical examples to show the power of the formulation. In conclusion, the 2D Fourier series representation on the sphere is generalized to functionals of the potential. In addition, the set of the spectral transfers can be considered as a pocket guide that provides the spectral characteristics of the functionals. Therefore, it extends the so-called Meissl scheme.
VLSI Implementation Of The Fast Fourier Transform
NASA Astrophysics Data System (ADS)
Chau, Paul M.; Ku, Walter H.
1986-03-01
A VLSI implementation of a Fast Fourier Transform (FFT) processor consisting of a mesh interconnection of complex floating-point butterfly units is presented. The Cooley-Tukey radix-2 Decimation-In-Frequency (DIF) formulation of the FFT was chosen since it offered the best overall compromise between the need for fast and efficient algorithmic computation and the need for a structure amenable to VLSI layout. Thus the VLSI implementation is modular, regular, expandable to various problem sizes and has a simple systolic flow of data and control. To evaluate the FFT architecture, VLSI area-time complexity concepts are used, but are now adapted to a complex floating-point number system rather than the usual integer ring representation. We show by our construction that the Thompson area-time optimum bound for the VLSI computation of an N-point FFT, area-time2oc = ORNlogN)1+a] can be attained by an alternative number representation, and hence the theoretical bound is a tight bound regardless of number system representation.
Fast modular data acquisition system for GEM-2D detector
NASA Astrophysics Data System (ADS)
Kasprowicz, G.; Byszuk, Adrian; Wojeński, A.; Zienkiewicz, P.; Czarski, T.; Chernyshova, M.; Poźniak, K.; Rzadkiewicz, J.; Zabolotny, W.; Juszczyk, B.
2014-11-01
A novel approach to two dimensional Gas Electron Multiplier (GEM) detector readout is presented. Unlike commonly used methods, based on discriminators and analogue FIFOs, the method developed uses simulta- neously sampling high speed ADCs with fast hybrid integrator and advanced FPGA-based processing logic to estimate the energy of every single photon. Such a method is applied to every GEM strip / pixel signal. It is especially useful in case of crystal-based spectrometers for soft X-rays, 2D imaging for plasma tomography and all these applications where energy resolution of every single photon is required. For the purpose of the detector readout, a novel, highly modular and extendable conception of the measurement platform was developed. It is evolution of already deployed measurement system for JET Spectrometer.
Fourier based methodology for simulating 2D-random shapes in heterogeneous materials
NASA Astrophysics Data System (ADS)
Mattrand, C.; Béakou, A.; Charlet, K.
2015-08-01
Gaining insights into the effects of microstructural details on materials behavior may be achieved by incorporating their attributes into numerical modeling. This requires us to make considerable efforts to feature heterogeneity morphology distributions and their spatial arrangement. This paper focuses on modeling the scatter observed in materials heterogeneity geometry. The proposed strategy is based on the development of a 1D-shape signature function representing the 2D-section of a given shape, on Fourier basis functions. The Fourier coefficients are then considered as random variables. This methodology has been applied to flax fibers which are gradually introduced into composite materials as a potential alternative to synthetic reinforcements. In this contribution, the influence of some underlying assumptions regarding the choice of one 1D-shape signature function, its discretization scheme and truncation level, and the best way of modeling the associated random variables is also investigated. Some configurations coming from the combination of these tuning parameters are found to be sufficiently relevant to render efficiently the morphometric factors of the observed fibers statistically speaking.
NASA Technical Reports Server (NTRS)
Alfano, Robert R. (Inventor); Cai, Wei (Inventor)
2007-01-01
A reconstruction technique for reducing computation burden in the 3D image processes, wherein the reconstruction procedure comprises an inverse and a forward model. The inverse model uses a hybrid dual Fourier algorithm that combines a 2D Fourier inversion with a 1D matrix inversion to thereby provide high-speed inverse computations. The inverse algorithm uses a hybrid transfer to provide fast Fourier inversion for data of multiple sources and multiple detectors. The forward model is based on an analytical cumulant solution of a radiative transfer equation. The accurate analytical form of the solution to the radiative transfer equation provides an efficient formalism for fast computation of the forward model.
Development of ultra-fast 2D ion Doppler tomography using image intensified CMOS fast camera
NASA Astrophysics Data System (ADS)
Tanabe, Hiroshi; Kuwahata, Akihiro; Yamanaka, Haruki; Inomoto, Michiaki; Ono, Yasushi; TS-group Team
2015-11-01
The world fastest novel time-resolved 2D ion Doppler tomography diagnostics has been developed using fast camera with high-speed gated image intensifier (frame rate: 200kfps. phosphor decay time: ~ 1 μ s). Time evolution of line-integrated spectra are diffracted from a f=1m, F/8.3 and g=2400L/mm Czerny-Turner polychromator, whose output is intensified and recorded to a high-speed camera with spectral resolution of ~0.005nm/pixel. The system can accommodate up to 36 (9 ×4) spatial points recorded at 5 μs time resolution, tomographic reconstruction is applied for the line-integrated spectra, time-resolved (5 μs/frame) local 2D ion temperature measurement has been achieved without any assumption of shot repeatability. Ion heating during intermittent reconnection event which tends to happen during high guide field merging tokamak was measured around diffusion region in UTST. The measured 2D profile shows ion heating inside the acceleration channel of reconnection outflow jet, stagnation point and downstream region where reconnected field forms thick closed flux surface as in MAST. Achieved maximum ion temperature increases as a function of Brec2 and shows good fit with MAST experiment, demonstrating promising CS-less startup scenario for spherical tokamak. This work is supported by JSPS KAKENHI Grant Number 15H05750 and 15K20921.
Fast Fourier transform analysis of rotor-bearing systems
NASA Technical Reports Server (NTRS)
Choy, K. C.; Gunter, E. J.; Allaire, P. E.
1978-01-01
Nonlinear transient analysis of rotor-bearing systems is becoming increasingly important in the analysis of modern-day rotating machinery to model such phenomena as oil film whirl. This paper develops an analysis technique incorporating modal analysis and fast Fourier transform techniques to analyze rotors with residual shaft bow and realistic nonlinear bearings. The technique is demonstrated on single-mass and three-mass rotor examples. Comparisons of the theoretical results with experimental data give excellent agreement.
Fast Fourier transformation results from gamma-ray burst profiles
NASA Technical Reports Server (NTRS)
Kouveliotou, Chryssa; Norris, Jay P.; Fishman, Gerald J.; Meegan, Charles A.; Wilson, Robert B.; Paciesas, W. S.
1992-01-01
Several gamma-ray bursts in the BATSE data have sufficiently long durations and complex temporal structures with pulses that appear to be spaced quasi-periodically. In order to test and quantify these periods we have applied fast Fourier transformations (FFT) to all these events. We have also performed cross spectral analyses of the FFT of the two extreme (high-low) energy bands in each case to determine the lead/lag of the pulses in different energies.
Preconditioning 2D Integer Data for Fast Convex Hull Computations.
Cadenas, José Oswaldo; Megson, Graham M; Luengo Hendriks, Cris L
2016-01-01
In order to accelerate computing the convex hull on a set of n points, a heuristic procedure is often applied to reduce the number of points to a set of s points, s ≤ n, which also contains the same hull. We present an algorithm to precondition 2D data with integer coordinates bounded by a box of size p × q before building a 2D convex hull, with three distinct advantages. First, we prove that under the condition min(p, q) ≤ n the algorithm executes in time within O(n); second, no explicit sorting of data is required; and third, the reduced set of s points forms a simple polygonal chain and thus can be directly pipelined into an O(n) time convex hull algorithm. This paper empirically evaluates and quantifies the speed up gained by preconditioning a set of points by a method based on the proposed algorithm before using common convex hull algorithms to build the final hull. A speedup factor of at least four is consistently found from experiments on various datasets when the condition min(p, q) ≤ n holds; the smaller the ratio min(p, q)/n is in the dataset, the greater the speedup factor achieved. PMID:26938221
Preconditioning 2D Integer Data for Fast Convex Hull Computations
2016-01-01
In order to accelerate computing the convex hull on a set of n points, a heuristic procedure is often applied to reduce the number of points to a set of s points, s ≤ n, which also contains the same hull. We present an algorithm to precondition 2D data with integer coordinates bounded by a box of size p × q before building a 2D convex hull, with three distinct advantages. First, we prove that under the condition min(p, q) ≤ n the algorithm executes in time within O(n); second, no explicit sorting of data is required; and third, the reduced set of s points forms a simple polygonal chain and thus can be directly pipelined into an O(n) time convex hull algorithm. This paper empirically evaluates and quantifies the speed up gained by preconditioning a set of points by a method based on the proposed algorithm before using common convex hull algorithms to build the final hull. A speedup factor of at least four is consistently found from experiments on various datasets when the condition min(p, q) ≤ n holds; the smaller the ratio min(p, q)/n is in the dataset, the greater the speedup factor achieved. PMID:26938221
Fast triangulated vortex methods for the 2D Eulen equations
NASA Astrophysics Data System (ADS)
Russo, Giovanni; Strain, John A.
1994-04-01
Vortex methods for inviscid incompressible two-dimensional fluid flow are usually based on blob approximations. This paper presents a vortex method in which the vorticity is approximated by a piecewise polynomial interpolant on a Delaunay triangulation of the vortices. An efficient reconstruction of the Delaunay triangulation at each step makes the method accurate for long times. The vertices of the triangulation move with the fluid velocity, which is reconstructed from the vorticity via a simplified fast multipole method for the Biot-Savart law with a continuous source distribution. The initial distribution of vortices is constructed from the initial vorticity field by an adaptive approximation method which produces good accuracy even for discontinuous initial data. Numerical results show that the method is highly accurate over long time intervals. Experiments with single and multiple circular and elliptical rotating patches of both piecewise constant and smooth vorticity indicate that the method produces much smaller errors than blob methods with the same number of degrees of freedom, at little additional cost. Generalizations to domains with boundaries, viscous flow, and three space dimensions are discussed.
Fast triangulated vortex methods for the 2D Euler equations
Russo, G. ); Strain, J.A. )
1994-04-01
Vortex methods for inviscid incompressible two-dimensional fluid flow are usually based on blob approximations. This paper presents a vortex method in which the vorticity is approximated by a piecewise polynomial interpolant on a Delaunay triangulation of the vortices. An efficient reconstruction of the Delaunay triangulation at each step makes the method accurate for long times. The vertices of the triangulation move with the fluid velocity, which is reconstructed from the vorticity via a simplified fast multipole method for the Biot-Savart law with a continuous source distribution. The initial distribution of vortices is constructed from the initial vorticity field by an adaptive approximation method which produces good accuracy even for discontinuous initial data. Numerical results show that the method is highly accurate over long time intervals. Experiments with single and multiple circular and elliptical rotating patches of both piecewise constant and smooth vorticity indicate that the method produces much smaller errors than blob methods with the same number of degrees of freedom, at little additional cost. Generalizations to domains with boundaries, viscous flow, and three space dimensions are discussed. 52 refs., 28 figs., 2 tabs.
Computing the Fast Fourier Transform on a vector computer
NASA Technical Reports Server (NTRS)
Korn, D. G.; Lambiotte, J. J., Jr.
1979-01-01
Two algorithms are presented for performing a Fast Fourier Transform on a vector computer and are compared on the Control Data Corporation STAR-100. The relative merits of the two algorithms are shown to depend upon whether only a few or many independent transforms are desired. A theorem is proved which shows that a set of independent transforms can be computed by performing a partial transformation on a single vector. The results of this theorem also apply to nonvector machines and have reduced the average time per transform by a factor of two on the CDC 6600 computer.
The Empirical Mode Decomposition algorithm via Fast Fourier Transform
NASA Astrophysics Data System (ADS)
Myakinin, Oleg O.; Zakharov, Valery P.; Bratchenko, Ivan A.; Kornilin, Dmitry V.; Artemyev, Dmitry N.; Khramov, Alexander G.
2014-09-01
In this paper we consider a problem of implementing a fast algorithm for the Empirical Mode Decomposition (EMD). EMD is one of the newest methods for decomposition of non-linear and non-stationary signals. A basis of EMD is formed "on-the-fly", i.e. it depends from a distribution of the signal and not given a priori in contrast on cases Fourier Transform (FT) or Wavelet Transform (WT). The EMD requires interpolating of local extrema sets of signal to find upper and lower envelopes. The data interpolation on an irregular lattice is a very low-performance procedure. A classical description of EMD by Huang suggests doing this through splines, i.e. through solving of a system of equations. Existence of a fast algorithm is the main advantage of the FT. A simple description of an algorithm in terms of Fast Fourier Transform (FFT) is a standard practice to reduce operation's count. We offer a fast implementation of EMD (FEMD) through FFT and some other cost-efficient algorithms. Basic two-stage interpolation algorithm for EMD is composed of a Upscale procedure through FFT and Downscale procedure through a selection procedure for signal's points. First we consider the local maxima (or minima) set without reference to the axis OX, i.e. on a regular lattice. The Upscale through the FFT change the signal's length to the Least Common Multiple (LCM) value of all distances between neighboring extremes on the axis OX. If the LCM value is too large then it is necessary to limit local set of extrema. In this case it is an analog of the spline interpolation. A demo for FEMD in noise reduction task for OCT has been shown.
NASA Astrophysics Data System (ADS)
Aizenberg, Evgeni; Bigio, Irving J.; Rodriguez-Diaz, Eladio
2012-03-01
The Fourier descriptors paradigm is a well-established approach for affine-invariant characterization of shape contours. In the work presented here, we extend this method to images, and obtain a 2D Fourier representation that is invariant to image rotation. The proposed technique retains phase uniqueness, and therefore structural image information is not lost. Rotation-invariant phase coefficients were used to train a single multi-valued neuron (MVN) to recognize satellite and human face images rotated by a wide range of angles. Experiments yielded 100% and 96.43% classification rate for each data set, respectively. Recognition performance was additionally evaluated under effects of lossy JPEG compression and additive Gaussian noise. Preliminary results show that the derived rotation-invariant features combined with the MVN provide a promising scheme for efficient recognition of rotated images.
NASA Astrophysics Data System (ADS)
Nash, J. Gregory
2002-07-01
A specialized CAD tool is described that will take a user's high level code description of a non-uniform affinely indexed algorithm and automatically generate abstract latency-optimal systolic arrays. Emphasis has been placed on ease of use and the ability to either force conformation to specific design criteria or perform unconstrained explorations. How such design goals are achieved is illustrated in the context of LU decomposition and the matrix Lyapunov equation. The tool is then used to generate new 1-D and 2-D hardware efficient systolic arrays for the discreet Fourier transform that take advantage of the use of the radix-4 matrix decomposition.
On the abstracted dataflow complexity of Fast Fourier Transforms
Boehm, A.P.W.; Hiromoto, R.E.; Kelly, K.A.; Ashley, J.M.
1992-05-01
In this paper we develop and analyze the simulated performance of codes for the Fast Fourier Transform written in If and targeted for execution on Motorola`s dataflow machine Monsoon. The FFT application is of interest because of its computational parallelism, its requirement for global communications, and its array element data dependences. We use the parallel profiling simulator Id World to study the dataflow performance of various implementations. Our approach is comparative. We study two approaches, a recursive and an iterative one, and in each version we examine the effect of a variety of implementations. We contend that only through such comparative evaluations can significant insight be gained in understanding the computational and structural details of functional algorithms.
Solar radiometry at millimeter wavelengths. [Fast Fourier Transformation solutions
NASA Technical Reports Server (NTRS)
Henze, W.
1974-01-01
In the area of resolution enhancement, the use of Fast Fourier Transform programs was investigated for possible application to millimeter wavelength maps of the sun. A difficulty arises with the La Posta maps in that they are limited to 35 arc-minutes square while the smeared out solar image is larger than that. A list of possible cometary emission lines near 13 millimeters is presented. Although preparation of the list was inspired by the appearance of Comet Kohoutek, the results are applicable to any future comet. The brightness temperature of the sun at 8.6 millimeters was measured using the moon as a calibration source. The result does not confirm a deep absorption feature as apparently observed by earlier workers.
On the abstracted dataflow complexity of Fast Fourier Transforms
Boehm, A.P.W. . Dept. of Computer Science); Hiromoto, R.E.; Kelly, K.A. ); Ashley, J.M. . Dept. of Computer Science)
1992-01-01
In this paper we develop and analyze the simulated performance of codes for the Fast Fourier Transform written in If and targeted for execution on Motorola's dataflow machine Monsoon. The FFT application is of interest because of its computational parallelism, its requirement for global communications, and its array element data dependences. We use the parallel profiling simulator Id World to study the dataflow performance of various implementations. Our approach is comparative. We study two approaches, a recursive and an iterative one, and in each version we examine the effect of a variety of implementations. We contend that only through such comparative evaluations can significant insight be gained in understanding the computational and structural details of functional algorithms.
Parameterized Spectral Bathymetric Roughness Using the Nonequispaced Fast Fourier Transform
NASA Astrophysics Data System (ADS)
Fabre, David Hanks
The ocean and acoustic modeling community has specifically asked for roughness from bathymetry. An effort has been undertaken to provide what can be thought of as the high frequency content of bathymetry. By contrast, the low frequency content of bathymetry is the set of contours. The two-dimensional amplitude spectrum calculated with the nonequispaced fast Fourier transform (Kunis, 2006) is exploited as the statistic to provide several parameters of roughness following the method of Fox (1996). When an area is uniformly rough, it is termed isotropically rough. When an area exhibits lineation effects (like in a trough or a ridge line in the bathymetry), the term anisotropically rough is used. A predominant spatial azimuth of lineation summarizes anisotropic roughness. The power law model fit produces a roll-off parameter that also provides insight into the roughness of the area. These four parameters give rise to several derived parameters. Algorithmic accomplishments include reviving Fox's method (1985, 1996) and improving the method with the possibly geophysically more appropriate nonequispaced fast Fourier transform. A new composite parameter, simply the overall integral length of the nonlinear parameterizing function, is used to make within-dataset comparisons. A synthetic dataset and six multibeam datasets covering practically all depth regimes have been analyzed with the tools that have been developed. Data specific contributions include possibly discovering an aspect ratio isotropic cutoff level (less than 1.2), showing a range of spectral fall-off values between about -0.5 for a sandybottomed Gulf of Mexico area, to about -1.8 for a coral reef area just outside of the Saipan harbor. We also rank the targeted type of dataset, the best resolution gridded datasets, from smoothest to roughest using a factor based on the kernel dimensions, a percentage from the windowing operation, all multiplied by the overall integration length.
Bi-planar 2D-to-3D registration in Fourier domain for stereoscopic x-ray motion tracking
NASA Astrophysics Data System (ADS)
Zosso, Dominique; Le Callennec, Benoît; Bach Cuadra, Meritxell; Aminian, Kamiar; Jolles, Brigitte M.; Thiran, Jean-Philippe
2008-03-01
In this paper we present a new method to track bone movements in stereoscopic X-ray image series of the knee joint. The method is based on two different X-ray image sets: a rotational series of acquisitions of the still subject knee that allows the tomographic reconstruction of the three-dimensional volume (model), and a stereoscopic image series of orthogonal projections as the subject performs movements. Tracking the movements of bones throughout the stereoscopic image series means to determine, for each frame, the best pose of every moving element (bone) previously identified in the 3D reconstructed model. The quality of a pose is reflected in the similarity between its theoretical projections and the actual radiographs. We use direct Fourier reconstruction to approximate the three-dimensional volume of the knee joint. Then, to avoid the expensive computation of digitally rendered radiographs (DRR) for pose recovery, we develop a corollary to the 3-dimensional central-slice theorem and reformulate the tracking problem in the Fourier domain. Under the hypothesis of parallel X-ray beams, the heavy 2D-to-3D registration of projections in the signal domain is replaced by efficient slice-to-volume registration in the Fourier domain. Focusing on rotational movements, the translation-relevant phase information can be discarded and we only consider scalar Fourier amplitudes. The core of our motion tracking algorithm can be implemented as a classical frame-wise slice-to-volume registration task. Results on both synthetic and real images confirm the validity of our approach.
Fast large scale structure perturbation theory using one-dimensional fast Fourier transforms
NASA Astrophysics Data System (ADS)
Schmittfull, Marcel; Vlah, Zvonimir; McDonald, Patrick
2016-05-01
The usual fluid equations describing the large-scale evolution of mass density in the universe can be written as local in the density, velocity divergence, and velocity potential fields. As a result, the perturbative expansion in small density fluctuations, usually written in terms of convolutions in Fourier space, can be written as a series of products of these fields evaluated at the same location in configuration space. Based on this, we establish a new method to numerically evaluate the 1-loop power spectrum (i.e., Fourier transform of the 2-point correlation function) with one-dimensional fast Fourier transforms. This is exact and a few orders of magnitude faster than previously used numerical approaches. Numerical results of the new method are in excellent agreement with the standard quadrature integration method. This fast model evaluation can in principle be extended to higher loop order where existing codes become painfully slow. Our approach follows by writing higher order corrections to the 2-point correlation function as, e.g., the correlation between two second-order fields or the correlation between a linear and a third-order field. These are then decomposed into products of correlations of linear fields and derivatives of linear fields. The method can also be viewed as evaluating three-dimensional Fourier space convolutions using products in configuration space, which may also be useful in other contexts where similar integrals appear.
Fast Fourier Transform Co-processor (FFTC), towards embedded GFLOPs
NASA Astrophysics Data System (ADS)
Kuehl, Christopher; Liebstueckel, Uwe; Tejerina, Isaac; Uemminghaus, Michael; Witte, Felix; Kolb, Michael; Suess, Martin; Weigand, Roland; Kopp, Nicholas
2012-10-01
Many signal processing applications and algorithms perform their operations on the data in the transform domain to gain efficiency. The Fourier Transform Co-Processor has been developed with the aim to offload General Purpose Processors from performing these transformations and therefore to boast the overall performance of a processing module. The IP of the commercial PowerFFT processor has been selected and adapted to meet the constraints of the space environment. In frame of the ESA activity "Fast Fourier Transform DSP Co-processor (FFTC)" (ESTEC/Contract No. 15314/07/NL/LvH/ma) the objectives were the following: • Production of prototypes of a space qualified version of the commercial PowerFFT chip called FFTC based on the PowerFFT IP. • The development of a stand-alone FFTC Accelerator Board (FTAB) based on the FFTC including the Controller FPGA and SpaceWire Interfaces to verify the FFTC function and performance. The FFTC chip performs its calculations with floating point precision. Stand alone it is capable computing FFTs of up to 1K complex samples in length in only 10μsec. This corresponds to an equivalent processing performance of 4.7 GFlops. In this mode the maximum sustained data throughput reaches 6.4Gbit/s. When connected to up to 4 EDAC protected SDRAM memory banks the FFTC can perform long FFTs with up to 1M complex samples in length or multidimensional FFT-based processing tasks. A Controller FPGA on the FTAB takes care of the SDRAM addressing. The instructions commanded via the Controller FPGA are used to set up the data flow and generate the memory addresses. The paper will give an overview on the project, including the results of the validation of the FFTC ASIC prototypes.
Xu, Min; Beck, Martin; Alber, Frank
2013-01-01
Cryo-electron tomography allows the visualization of macromolecular complexes in their cellular environments in close-to-live conditions. The nominal resolution of subtomograms can be significantly increased when individual subtomograms of the same kind are aligned and averaged. A vital step for such a procedure are algorithms that speedup subtomogram alignment and improve accuracy for reference-free subtomogram classification, which will facilitate automation of tomography analysis and overall high throughput in the data processing. In this paper, we propose a fast rotational alignment method that uses the Fourier equivalent form of a popular constrained correlation measure that considers missing wedge corrections and density variances in the subtomograms. The fast rotational search is based on 3D volumetric matching, which significantly improves the rotational alignment accuracy in particular for highly distorted subtomograms with low SNR and tilt angle ranges in comparison to a fast rotational alignment based on matching of projected 2D spherical images. We further integrate our fast rotational alignment method in a reference free iterative subtomogram classification scheme, and propose a local feature enhancement strategy in the classification process. We can demonstrate that the automatic method can be used to successfully classify a large number of experimental subtomograms without the need of a reference structure. PMID:22420977
Gao, Wenjing; Huyen, Nguyen Thi Thanh; Loi, Ho Sy; Kemao, Qian
2009-12-01
In optical interferometers, fringe projection systems, and synthetic aperture radars, fringe patterns are common outcomes and usually degraded by unavoidable noises. The presence of noises makes the phase extraction and phase unwrapping challenging. Windowed Fourier transform (WFT) based algorithms have been proven to be effective for fringe pattern analysis to various applications. However, the WFT-based algorithms are computationally expensive, prohibiting them from real-time applications. In this paper, we propose a fast parallel WFT-based library using graphics processing units and computer unified device architecture. Real-time WFT-based algorithms are achieved with 4 frames per second in processing 256x256 fringe patterns. Up to 132x speedup is obtained for WFT-based algorithms using NVIDIA GTX295 graphics card than sequential C in quad-core 2.5GHz Intel(R)Xeon(R) CPU E5420. PMID:20052242
Advances in fast 2D camera data handling and analysis on NSTX
Davis, W. M.; Patel, R. I.; Boeglin, W. U.; Roquemore, A. L.; Maqueda, R. J.; Zweben, S. J.
2010-07-01
The use of fast 2D cameras on NSTX continues to grow. There are 6 cameras with the capability of taking up to 1–2 gigabytes (GBs) of data apiece during each plasma shot on the National Spherical Torus Experiment (NSTX). Efficient storage and retrieval of this data remains a challenge. Performance comparisons are presented for reading data stored in MDSplus, using both compressed data and segmented records, and direct access I/O with different read sizes. Encouragingly, fast 2D camera data provides considerable insight into plasma complexities, such as small-scale turbulence and particle transport. The last part of this paper is an example of more recent uses: dual cameras looking at the same region of the plasma from different angles, which can provide trajectories of incandescent particles in 3D. A laboratory simulation of the 3D trajectories is presented, as well as corresponding data from NSTX plasma where glowing dust particles can be followed.
Cryo-EM Image Alignment Based on Nonuniform Fast Fourier Transform
Yang, Zhengfan; Penczek, Pawel A.
2008-01-01
In single particle analysis, two-dimensional (2-D) alignment is a fundamental step intended to put into register various particle projections of biological macromolecules collected at the electron microscope. The efficiency and quality of three-dimensional (3-D) structure reconstruction largely depends on the computational speed and alignment accuracy of this crucial step. In order to improve the performance of alignment, we introduce a new method that takes advantage of the highly accurate interpolation scheme based on the gridding method, a version of the nonuniform Fast Fourier Transform, and utilizes a multi-dimensional optimization algorithm for the refinement of the orientation parameters. Using simulated data, we demonstrate that by using less than half of the sample points and taking twice the runtime, our new 2-D alignment method achieves dramatically better alignment accuracy than that based on quadratic interpolation. We also apply our method to image to volume registration, the key step in the single particle EM structure refinement protocol. We find that in this case the accuracy of the method not only surpasses the accuracy of the commonly used real-space implementation, but results are achieved in much shorter time, making gridding-based alignment a perfect candidate for efficient structure determination in single particle analysis. PMID:18499351
NASA Astrophysics Data System (ADS)
Kwiecien, Pavel; Richter, Ivan; Čtyroký, Jiří
2012-02-01
Frequency-domain Fourier modal methods have recently evolved into efficient tools for rigorous numerical modeling of a wide class of photonic and plasmonic structures and devices. In this contribution we describe the application of our 2D and 3D in-house tools, namely aperiodic rigorous coupled wave analysis (aRCWA) and bi-directional mode expansion propagation method using harmonic expansion (BEXX), on a recently described novel type of subwavelength grating (SWG) waveguides. They are created by means of periodically interlacing silicon segments with a superstrate material with a lower refractive index. It has been shown recently, both theoretically and experimentally, that for a suitable choice of SWG parameters such as grating period and duty cycle, the structure can support low-loss guided (Bloch) mode. Its effective index, mode profile and dispersion characteristics can thus be tailored to specific needs without the necessity of changing material composition. In our methods, either complex coordinate transformation or uniaxial anisotropic perfectly matched layers have been applied as efficient absorption boundary conditions. In order to reduce the number of expansion terms needed to reach required accuracy, the adaptive spatial resolution technique has been implemented. Structural symmetries of the devices can be fully utilized to this aim, too. Propagation constants of Bloch modes are also compared with those obtained with a full-vector film mode matching (FiMM) mode solver using the very simple effective medium theory (EMT).
NASA Astrophysics Data System (ADS)
Kwiecien, Pavel; Richter, Ivan; Čtyroký, Jiří
2011-09-01
Frequency-domain Fourier modal methods have recently evolved into efficient tools for rigorous numerical modeling of a wide class of photonic and plasmonic structures and devices. In this contribution we describe the application of our 2D and 3D in-house tools, namely aperiodic rigorous coupled wave analysis (aRCWA) and bi-directional mode expansion propagation method using harmonic expansion (BEXX), on a recently described novel type of subwavelength grating (SWG) waveguides. They are created by means of periodically interlacing silicon segments with a superstrate material with a lower refractive index. It has been shown recently, both theoretically and experimentally, that for a suitable choice of SWG parameters such as grating period and duty cycle, the structure can support low-loss guided (Bloch) mode. Its effective index, mode profile and dispersion characteristics can thus be tailored to specific needs without the necessity of changing material composition. In our methods, either complex coordinate transformation or uniaxial anisotropic perfectly matched layers have been applied as efficient absorption boundary conditions. In order to reduce the number of expansion terms needed to reach required accuracy, the adaptive spatial resolution technique has been implemented. Structural symmetries of the devices can be fully utilized to this aim, too. Propagation constants of Bloch modes are also compared with those obtained with a full-vector film mode matching (FiMM) mode solver using the very simple effective medium theory (EMT).
Optimal Padding for the Two-Dimensional Fast Fourier Transform
NASA Technical Reports Server (NTRS)
Dean, Bruce H.; Aronstein, David L.; Smith, Jeffrey S.
2011-01-01
One-dimensional Fast Fourier Transform (FFT) operations work fastest on grids whose size is divisible by a power of two. Because of this, padding grids (that are not already sized to a power of two) so that their size is the next highest power of two can speed up operations. While this works well for one-dimensional grids, it does not work well for two-dimensional grids. For a two-dimensional grid, there are certain pad sizes that work better than others. Therefore, the need exists to generalize a strategy for determining optimal pad sizes. There are three steps in the FFT algorithm. The first is to perform a one-dimensional transform on each row in the grid. The second step is to transpose the resulting matrix. The third step is to perform a one-dimensional transform on each row in the resulting grid. Steps one and three both benefit from padding the row to the next highest power of two, but the second step needs a novel approach. An algorithm was developed that struck a balance between optimizing the grid pad size with prime factors that are small (which are optimal for one-dimensional operations), and with prime factors that are large (which are optimal for two-dimensional operations). This algorithm optimizes based on average run times, and is not fine-tuned for any specific application. It increases the amount of times that processor-requested data is found in the set-associative processor cache. Cache retrievals are 4-10 times faster than conventional memory retrievals. The tested implementation of the algorithm resulted in faster execution times on all platforms tested, but with varying sized grids. This is because various computer architectures process commands differently. The test grid was 512 512. Using a 540 540 grid on a Pentium V processor, the code ran 30 percent faster. On a PowerPC, a 256x256 grid worked best. A Core2Duo computer preferred either a 1040x1040 (15 percent faster) or a 1008x1008 (30 percent faster) grid. There are many industries that
Parallel fast Fourier transforms for non power of two data
Semeraro, B.D.
1994-09-01
This report deals with parallel algorithms for computing discrete Fourier transforms of real sequences of length N not equal to a power of two. The method described is an extension of existing power of two transforms to sequences with N a product of small primes. In particular, this implementation requires N = 2{sup p}3{sup q}5{sup r}. The communication required is the same as for a transform of length N = 2{sup p}. The algorithm presented is intended for use in the solution of partial differential equations, or in any situation in which a large number of forward and backward transforms must be performed and in which the Fourier Coefficients need not be ordered. This implementation is a one dimensional FFT but the techniques are applicable to multidimensional transforms as well. The algorithm has been implemented on a 128 node Intel Ipsc/860.
Fast acquisition of high-resolution 2D NMR spectroscopy in inhomogeneous magnetic fields
NASA Astrophysics Data System (ADS)
Lin, Liangjie; Wei, Zhiliang; Zeng, Qing; Yang, Jian; Lin, Yanqin; Chen, Zhong
2016-05-01
High-resolution nuclear magnetic resonance (NMR) spectroscopy plays an important role in chemical and biological analyses. In this study, we combine the J-coupling coherence transfer module with the echo-train acquisition technique for fast acquisition of high-resolution 2D NMR spectra in magnetic fields with unknown spatial variations. The proposed method shows satisfactory performance on a 5 mM ethyl 3-bromopropionate sample, under a 5-kHz (10 ppm at 11.7 T) B0 inhomogeneous field, as well as under varying degrees of pulse-flip-angle deviations. Moreover, a simulative ex situ NMR measurement is also conducted to show the effectiveness of the proposed pulse sequence.
A feasibility study using radiochromic films for fast neutron 2D passive dosimetry
Brady, Samuel L; Gunasingha, Rathnayaka; Yoshizumi, Terry T; Howell, Calvin R; Crowell, Alexander S; Fallin, Brent; Tonchev, Anton P; Dewhirst, Mark W
2013-01-01
The objective of this paper is threefold: (1) to establish sensitivity of XRQA and EBT radiochromic films to fast neutron exposure; (2) to develop a film response to radiation dose calibration curve and (3) to investigate a two-dimensional (2D) film dosimetry technique for use in establishing an experimental setup for a radiobiological irradiation of mice and to assess the dose to the mice in this setup. The films were exposed to a 10 MeV neutron beam via the 2H(d,n)3He reaction. The XRQA film response was a factor of 1.39 greater than EBT film response to the 10 MeV neutron beam when exposed to a neutron dose of 165 cGy. A film response-to-soft tissue dose calibration function was established over a range of 0–10 Gy and had a goodness of fit of 0.9926 with the calibration data. The 2D film dosimetry technique estimated the neutron dose to the mice by measuring the dose using a mouse phantom and by placing a piece of film on the exterior of the experimental mouse setup. The film results were benchmarked using Monte Carlo and aluminum (Al) foil activation measurements. The radiochromic film, Monte Carlo and Al foil dose measurements were strongly correlated, and the film within the mouse phantom agreed to better than 7% of the externally mounted films. These results demonstrated the potential application of radiochromic films for passive 2D neutron dosimetry. PMID:20693612
Lerche, Ernesto; Van Eester, Dirk
2011-12-23
Fourier analysis in the poloidal direction is a standard ingredient in present-day 2D wave equation solvers describing radio frequency waves in hot tokamak plasmas. Although a powerful and elegant technique, Fourier analysis has the disadvantage that a large number of modes is needed to describe the field pattern on a magnetic surface if a short wavelength mode exists on any - even very small - subpart of the particle trajectory. The present paper examines the potential of a method that does not suffer from this drawback: a finite element technique relying on simple linear or cubic area base functions that are defined on irregular elementary surfaces of triangular shape. The wave equation is solved in its weak Galerkin variational form and for realistic 2D tokamak geometry, accounting for the toroidal curvature but assuming the toroidal angle is ignorable, allowing to study the wave pattern for each of the independent toroidal modes excited by the antenna individually.The locally uniform full hot plasma dielectric tensor to all orders in finite Larmor radius was adopted. As the main intended application is the study of fast wave behavior (heating and current drive) at arbitrary harmonics, the wave vector complex amplitude appearing in the dielectric tensor is determined through a local dispersion root evaluation. High frequency fast wave propagation and damping is provided as an illustration in view of possible application of this type of current drive in future high density reactor-like tokamaks.
A new hybrid algorithm for computing a fast discrete Fourier transform
NASA Technical Reports Server (NTRS)
Reed, I. S.; Truong, T. K.
1979-01-01
In this paper for certain long transform lengths, Winograd's algorithm for computing the discrete Fourier transform (DFT) is extended considerably. This is accomplished by performing the cyclic convolution, required by Winograd's method, with the Mersenne prime number-theoretic transform developed originally by Rader. This new algorithm requires fewer multiplications than either the standard fast Fourier transform (FFT) or Winograd's more conventional algorithm. However, more additions are required.
Fast Confocal Raman Imaging Using a 2-D Multifocal Array for Parallel Hyperspectral Detection.
Kong, Lingbo; Navas-Moreno, Maria; Chan, James W
2016-01-19
We present the development of a novel confocal hyperspectral Raman microscope capable of imaging at speeds up to 100 times faster than conventional point-scan Raman microscopy under high noise conditions. The microscope utilizes scanning galvomirrors to generate a two-dimensional (2-D) multifocal array at the sample plane, generating Raman signals simultaneously at each focus of the array pattern. The signals are combined into a single beam and delivered through a confocal pinhole before being focused through the slit of a spectrometer. To separate the signals from each row of the array, a synchronized scan mirror placed in front of the spectrometer slit positions the Raman signals onto different pixel rows of the detector. We devised an approach to deconvolve the superimposed signals and retrieve the individual spectra at each focal position within a given row. The galvomirrors were programmed to scan different focal arrays following Hadamard encoding patterns. A key feature of the Hadamard detection is the reconstruction of individual spectra with improved signal-to-noise ratio. Using polystyrene beads as test samples, we demonstrated not only that our system images faster than a conventional point-scan method but that it is especially advantageous under noisy conditions, such as when the CCD detector operates at fast read-out rates and high temperatures. This is the first demonstration of multifocal confocal Raman imaging in which parallel spectral detection is implemented along both axes of the CCD detector chip. We envision this novel 2-D multifocal spectral detection technique can be used to develop faster imaging spontaneous Raman microscopes with lower cost detectors. PMID:26654100
Unger, Miriam; Pfeifer, Frank; Siesler, Heinz W
2016-07-01
The main objective of this communication is to compare the performance of a miniaturized handheld near-infrared (NIR) spectrometer with a benchtop Fourier transform near-infrared (FT-NIR) spectrometer. Generally, NIR spectroscopy is an extremely powerful analytical tool to study hydrogen-bonding changes of amide functionalities in solid and liquid materials and therefore variable temperature NIR measurements of polyamide II (PAII) have been selected as a case study. The information content of the measurement data has been further enhanced by exploiting the potential of two-dimensional correlation spectroscopy (2D-COS) and the perturbation correlation moving window two-dimensional (PCMW2D) evaluation technique. The data provide valuable insights not only into the changes of the hydrogen-bonding structure and the recrystallization of the hydrocarbon segments of the investigated PAII but also in their sequential order. Furthermore, it has been demonstrated that the 2D-COS and PCMW2D results derived from the spectra measured with the miniaturized NIR instrument are equivalent to the information extracted from the data obtained with the high-performance FT-NIR instrument. PMID:27287846
An algorithm for computing the 2D structure of fast rotating stars
NASA Astrophysics Data System (ADS)
Rieutord, Michel; Espinosa Lara, Francisco; Putigny, Bertrand
2016-08-01
Stars may be understood as self-gravitating masses of a compressible fluid whose radiative cooling is compensated by nuclear reactions or gravitational contraction. The understanding of their time evolution requires the use of detailed models that account for a complex microphysics including that of opacities, equation of state and nuclear reactions. The present stellar models are essentially one-dimensional, namely spherically symmetric. However, the interpretation of recent data like the surface abundances of elements or the distribution of internal rotation have reached the limits of validity of one-dimensional models because of their very simplified representation of large-scale fluid flows. In this article, we describe the ESTER code, which is the first code able to compute in a consistent way a two-dimensional model of a fast rotating star including its large-scale flows. Compared to classical 1D stellar evolution codes, many numerical innovations have been introduced to deal with this complex problem. First, the spectral discretization based on spherical harmonics and Chebyshev polynomials is used to represent the 2D axisymmetric fields. A nonlinear mapping maps the spheroidal star and allows a smooth spectral representation of the fields. The properties of Picard and Newton iterations for solving the nonlinear partial differential equations of the problem are discussed. It turns out that the Picard scheme is efficient on the computation of the simple polytropic stars, but Newton algorithm is unsurpassed when stellar models include complex microphysics. Finally, we discuss the numerical efficiency of our solver of Newton iterations. This linear solver combines the iterative Conjugate Gradient Squared algorithm together with an LU-factorization serving as a preconditioner of the Jacobian matrix.
Fast 3D shape measurement using Fourier transform profilometry without phase unwrapping
NASA Astrophysics Data System (ADS)
Song, Kechen; Hu, Shaopeng; Wen, Xin; Yan, Yunhui
2016-09-01
This paper presents a novel, simple, yet fast 3D shape measurement method using Fourier transform profilometry. Different from the conventional Fourier transform profilometry, this proposed method introduces the binocular stereo vision and employs two image pairs (i.e., original image pairs and fringe image pairs) to restructure 3D shape. In this proposed method, instead of phase unwrapping algorithm, a coarse disparity map is adopted as a constraint condition to realize phase matching using wrapped phase. Since the local phase matching and sub-pixel disparity refinement are proposed to obtain high measuring accuracy, high-quality phase is not required. The validity of the proposed method is verified by experiments.
Gridding and fast Fourier transformation on non-uniformly sparse sampled multidimensional NMR data.
Jiang, Bin; Jiang, Xianwang; Xiao, Nan; Zhang, Xu; Jiang, Ling; Mao, Xi-an; Liu, Maili
2010-05-01
For multidimensional NMR method, indirect dimensional non-uniform sparse sampling can dramatically shorten acquisition time of the experiments. However, the non-uniformly sampled NMR data cannot be processed directly using fast Fourier transform (FFT). We show that the non-uniformly sampled NMR data can be reconstructed to Cartesian grid with the gridding method that has been wide applied in MRI, and sequentially be processed using FFT. The proposed gridding-FFT (GFFT) method increases the processing speed sharply compared with the previously proposed non-uniform Fourier Transform, and may speed up application of the non-uniform sparse sampling approaches. PMID:20236843
Gridding and fast Fourier transformation on non-uniformly sparse sampled multidimensional NMR data
NASA Astrophysics Data System (ADS)
Jiang, Bin; Jiang, Xianwang; Xiao, Nan; Zhang, Xu; Jiang, Ling; Mao, Xi-an; Liu, Maili
2010-05-01
For multidimensional NMR method, indirect dimensional non-uniform sparse sampling can dramatically shorten acquisition time of the experiments. However, the non-uniformly sampled NMR data cannot be processed directly using fast Fourier transform (FFT). We show that the non-uniformly sampled NMR data can be reconstructed to Cartesian grid with the gridding method that has been wide applied in MRI, and sequentially be processed using FFT. The proposed gridding-FFT (GFFT) method increases the processing speed sharply compared with the previously proposed non-uniform Fourier Transform, and may speed up application of the non-uniform sparse sampling approaches.
NASA Astrophysics Data System (ADS)
Blacic, Tanya M.; Jun, Hyunggu; Rosado, Hayley; Shin, Changsoo
2016-02-01
In seismic oceanography, processed images highlight small temperature changes, but inversion is needed to obtain absolute temperatures. Local search-based full waveform inversion has a lower computational cost than global search but requires accurate starting models. Unfortunately, most marine seismic data have little associated hydrographic data and the band-limited nature of seismic data makes extracting the long wavelength sound speed trend directly from seismic data inherently challenging. Laplace and Laplace-Fourier domain inversion (LDI) can use rudimentary starting models without prior information about the medium. Data are transformed to the Laplace domain, and a smooth sound speed model is extracted by examining the zero and low frequency components of the damped wavefield. We applied LDI to five synthetic data sets based on oceanographic features and recovered smoothed versions of our synthetic models, showing the viability of LDI for creating starting models suitable for more detailed inversions.
NASA Astrophysics Data System (ADS)
Tan, T. L.; Gabona, M. G.; Woo, J. Q.; Ng, L. L.; Wong, Andy; McNaughton, Don
2016-03-01
The Fourier transform infrared (FTIR) rovibrational spectrum of the b-type ν9 band of 13C2D4 was recorded at a unapodized resolution of 0.0063 cm-1 in the 2230-2450 cm-1 region. By assigning and fitting a total of 1171 rotationally resolved infrared transitions of the ν9 band and using the Watson's A-reduced Hamiltonian in the Ir representation, rovibrational constants for the upper state (v9 = 1) up to five quartic centrifugal distortion terms were derived for the first time. The root-mean-square (rms) deviation of the infrared fit was 0.00043 cm-1. The ground state constants of 13C2D4 were determined with higher accuracy than previously by fitting 1485 ground state combination-differences (GSCDs) from the present and previous infrared measurements, with rms deviation of 0.00034 cm-1. The band center of ν9 band of 13C2D4 was determined to be at 2324.3593 cm-1. The equilibrium state rovibrational constants up to all 5 quartic terms were derived from theoretical harmonic calculations at three levels of theory: B3LYP/cc-pVTZ, MP2/cc-pVTZ, and CCSD/cc-pVTZ using the principal axis coordinate system. These constants agreed reasonably well with the ground state constants of 13C2D4 derived from the present experimental GSCD fit. Furthermore, all 3 rotational constants of the upper state (v9 = 1) and of the ground state of 13C2D4 were obtained from anharmonic calculations using B3LYP and MP2 levels with the cc-pVTZ basis set. The calculated rotational constants were found to agree with those derived experimentally within 0.40%.
Program for the analysis of time series. [by means of fast Fourier transform algorithm
NASA Technical Reports Server (NTRS)
Brown, T. J.; Brown, C. G.; Hardin, J. C.
1974-01-01
A digital computer program for the Fourier analysis of discrete time data is described. The program was designed to handle multiple channels of digitized data on general purpose computer systems. It is written, primarily, in a version of FORTRAN 2 currently in use on CDC 6000 series computers. Some small portions are written in CDC COMPASS, an assembler level code. However, functional descriptions of these portions are provided so that the program may be adapted for use on any facility possessing a FORTRAN compiler and random-access capability. Properly formatted digital data are windowed and analyzed by means of a fast Fourier transform algorithm to generate the following functions: (1) auto and/or cross power spectra, (2) autocorrelations and/or cross correlations, (3) Fourier coefficients, (4) coherence functions, (5) transfer functions, and (6) histograms.
NASA Astrophysics Data System (ADS)
Vilardy, Juan M.; Giacometto, F.; Torres, C. O.; Mattos, L.
2011-01-01
The two-dimensional Fast Fourier Transform (FFT 2D) is an essential tool in the two-dimensional discrete signals analysis and processing, which allows developing a large number of applications. This article shows the description and synthesis in VHDL code of the FFT 2D with fixed point binary representation using the programming tool Simulink HDL Coder of Matlab; showing a quick and easy way to handle overflow, underflow and the creation registers, adders and multipliers of complex data in VHDL and as well as the generation of test bench for verification of the codes generated in the ModelSim tool. The main objective of development of the hardware architecture of the FFT 2D focuses on the subsequent completion of the following operations applied to images: frequency filtering, convolution and correlation. The description and synthesis of the hardware architecture uses the XC3S1200E family Spartan 3E FPGA from Xilinx Manufacturer.
Simulation of micromechanical behavior of polycrystals: finite elements vs. fast Fourier transforms
Lebensohn, Ricardo A; Prakash, Arun
2009-01-01
In this work, we compare finite element and fast Fourier transform approaches for the prediction of micromechanical behavior of polycrystals. Both approaches are full-field approaches and use the same visco-plastic single crystal constitutive law. We investigate the texture and the heterogeneity of the inter- and intragranular, stress and strain fields obtained from the two models. Additionally, we also look into their computational performance. Two cases - rolling of aluminium and wire drawing of tungsten - are used to evaluate the predictions of the two mode1s. Results from both the models are similar, when large grain distortions do not occur in the polycrystal. The finite element simulations were found to be highly computationally intensive, in comparison to the fast Fourier transform simulations.
Using single buffers and data reorganization to implement a multi-megasample fast Fourier transform
NASA Technical Reports Server (NTRS)
Brown, R. D.
1992-01-01
Data ordering in large fast Fourier transforms (FFT's) is both conceptually and implementationally difficult. Discribed here is a method of visualizing data orderings as vectors of address bits, which enables the engineer to use more efficient data orderings and reduce double-buffer memory designs. Also detailed are the difficulties and algorithmic solutions involved in FFT lengths up to 4 megasamples (Msamples) and sample rates up to 80 MHz.
Isotropic Spin Trap EPR Spectra Simulation by Fast Fourier Transform (FFT)
NASA Astrophysics Data System (ADS)
Laachir, S.; Moussetad, M.; Adhiri, R.; Fahli, A.
2005-03-01
The detection and investigation of free radicals forming in living systems became possible due to the introduction of the method of spin traps. In this work, the electron spin resonance (ESR) spectra of DMPO/HO(.) and MGD-Fe-NO adducts are reproduced by simulation, based on the Fast Fourier Transform (FFT). The calculated spectral parameters as the hyperfine coupling constants, agree reasonably with the experimental data and the results are discussed.
NASA Technical Reports Server (NTRS)
Powers, E. J.; Kim, Y. C.; Hong, J. Y.; Roth, J. R.; Krawczonek, W. M.
1978-01-01
A diagnostic, based on fast Fourier-transform spectral analysis techniques, that provides experimental insight into the relationship between the experimentally observable spectral characteristics of the fluctuations and the fluctuation-induced plasma transport is described. The model upon which the diagnostic technique is based and its experimental implementation is discussed. Some characteristic results obtained during the course of an experimental study of fluctuation-induced transport in the electric field dominated NASA Lewis bumpy torus plasma are presented.
Fast 2D flood modelling using GPU technology - recent applications and new developments
NASA Astrophysics Data System (ADS)
Crossley, Amanda; Lamb, Rob; Waller, Simon; Dunning, Paul
2010-05-01
In recent years there has been considerable interest amongst scientists and engineers in exploiting the potential of commodity graphics hardware for desktop parallel computing. The Graphics Processing Units (GPUs) that are used in PC graphics cards have now evolved into powerful parallel co-processors that can be used to accelerate the numerical codes used for floodplain inundation modelling. We report in this paper on experience over the past two years in developing and applying two dimensional (2D) flood inundation models using GPUs to achieve significant practical performance benefits. Starting with a solution scheme for the 2D diffusion wave approximation to the 2D Shallow Water Equations (SWEs), we have demonstrated the capability to reduce model run times in ‘real-world' applications using GPU hardware and programming techniques. We then present results from a GPU-based 2D finite volume SWE solver. A series of numerical test cases demonstrate that the model produces outputs that are accurate and consistent with reference results published elsewhere. In comparisons conducted for a real world test case, the GPU-based SWE model was over 100 times faster than the CPU version. We conclude with some discussion of practical experience in using the GPU technology for flood mapping applications, and for research projects investigating use of Monte Carlo simulation methods for the analysis of uncertainty in 2D flood modelling.
NASA Astrophysics Data System (ADS)
Li, Dong Jun; Huang, Zhegang; Hwang, Tae Hoon; Narayan, Rekha; Choi, Jang Wook; Kim, Sang Ouk
2016-03-01
Realizing practical high performance materials and devices using the properties of 2D materials is of key research interest in the materials science field. In particular, building well-defined heterostructures using more than two different 2D components in a rational way is highly desirable. In this paper, a 2D heterostructure consisting of atomic thin titania nanosheets densely grown on reduced graphene oxide surface is successfully prepared through incorporating polymer functionalized graphene oxide into the novel TiO2 nanosheets synthesis scheme. As a result of the synergistic combination of a highly accessible surface area and abundant interface, which can modulate the physicochemical properties, the resultant heterostructure can be used in high efficiency visible light photocatalysis as well as fast energy storage with a long lifecycle. [Figure not available: see fulltext.
Analysis of the High-Resolution Fourier Spectrum of the ν6 Band of the cis-C2h2d2 Molecule
NASA Astrophysics Data System (ADS)
Konov, I. A.; Chertavskikh, Yu. V.; Fomchenko, A. L.; Aslapovskaya, Yu. S.; Zhdanovich, S. A.; Sydow, C.
2016-03-01
The spectrum of the ν6 band of the cis-ethylene-d2 molecule (cis-C2H2D2) is recorded with a Bruker IFS 125 HR Fourier spectrometer in the range 580-1210 cm-1 with resolution of 0.0021 cm-1. An analysis of the experimental spectrum allows more than 1500 transitions belonging to this band to be assigned that by more than 2.5 times greater than it has been known in the literature so far. The obtained experimental data are then used to determine the model parameters of the molecule (the effective Hamiltonian in the A-reduction and I'- representation). Strong resonance interaction with the band ν4 forbidden in absorption by the symmetry of a molecule is taken into account. 10 parameters of the Hamiltonian obtained by solving inverse spectroscopic problem reproduce 427 initial experimental energies (more than 1500 transitions) with accuracy close to the experimental uncertainty.
Suppression law of quantum states in a 3D photonic fast Fourier transform chip.
Crespi, Andrea; Osellame, Roberto; Ramponi, Roberta; Bentivegna, Marco; Flamini, Fulvio; Spagnolo, Nicolò; Viggianiello, Niko; Innocenti, Luca; Mataloni, Paolo; Sciarrino, Fabio
2016-01-01
The identification of phenomena able to pinpoint quantum interference is attracting large interest. Indeed, a generalization of the Hong-Ou-Mandel effect valid for any number of photons and optical modes would represent an important leap ahead both from a fundamental perspective and for practical applications, such as certification of photonic quantum devices, whose computational speedup is expected to depend critically on multi-particle interference. Quantum distinctive features have been predicted for many particles injected into multimode interferometers implementing the Fourier transform over the optical modes. Here we develop a scalable approach for the implementation of the fast Fourier transform algorithm using three-dimensional photonic integrated interferometers, fabricated via femtosecond laser writing technique. We observe the suppression law for a large number of output states with four- and eight-mode optical circuits: the experimental results demonstrate genuine quantum interference between the injected photons, thus offering a powerful tool for diagnostic of photonic platforms. PMID:26843135
Suppression law of quantum states in a 3D photonic fast Fourier transform chip
Crespi, Andrea; Osellame, Roberto; Ramponi, Roberta; Bentivegna, Marco; Flamini, Fulvio; Spagnolo, Nicolò; Viggianiello, Niko; Innocenti, Luca; Mataloni, Paolo; Sciarrino, Fabio
2016-01-01
The identification of phenomena able to pinpoint quantum interference is attracting large interest. Indeed, a generalization of the Hong–Ou–Mandel effect valid for any number of photons and optical modes would represent an important leap ahead both from a fundamental perspective and for practical applications, such as certification of photonic quantum devices, whose computational speedup is expected to depend critically on multi-particle interference. Quantum distinctive features have been predicted for many particles injected into multimode interferometers implementing the Fourier transform over the optical modes. Here we develop a scalable approach for the implementation of the fast Fourier transform algorithm using three-dimensional photonic integrated interferometers, fabricated via femtosecond laser writing technique. We observe the suppression law for a large number of output states with four- and eight-mode optical circuits: the experimental results demonstrate genuine quantum interference between the injected photons, thus offering a powerful tool for diagnostic of photonic platforms. PMID:26843135
Dynamic UltraFast 2D EXchange SpectroscopY (UF-EXSY) of hyperpolarized substrates
NASA Astrophysics Data System (ADS)
Leon Swisher, Christine; Koelsch, Bertram; Sukumar, Subramianam; Sriram, Renuka; Santos, Romelyn Delos; Wang, Zhen Jane; Kurhanewicz, John; Vigneron, Daniel; Larson, Peder
2015-08-01
In this work, we present a new ultrafast method for acquiring dynamic 2D EXchange SpectroscopY (EXSY) within a single acquisition. This technique reconstructs two-dimensional EXSY spectra from one-dimensional spectra based on the phase accrual during echo times. The Ultrafast-EXSY acquisition overcomes long acquisition times typically needed to acquire 2D NMR data by utilizing sparsity and phase dependence to dramatically undersample in the indirect time dimension. This allows for the acquisition of the 2D spectrum within a single shot. We have validated this method in simulations and hyperpolarized enzyme assay experiments separating the dehydration of pyruvate and lactate-to-pyruvate conversion. In a renal cell carcinoma cell (RCC) line, bidirectional exchange was observed. This new technique revealed decreased conversion of lactate-to-pyruvate with high expression of monocarboxylate transporter 4 (MCT4), known to correlate with aggressive cancer phenotypes. We also showed feasibility of this technique in vivo in a RCC model where bidirectional exchange was observed for pyruvate-lactate, pyruvate-alanine, and pyruvate-hydrate and were resolved in time. Broadly, the technique is well suited to investigate the dynamics of multiple exchange pathways and applicable to hyperpolarized substrates where chemical exchange has shown great promise across a range of disciplines.
Fast and Memory-Efficient Topological Denoising of 2D and 3D Scalar Fields.
Günther, David; Jacobson, Alec; Reininghaus, Jan; Seidel, Hans-Peter; Sorkine-Hornung, Olga; Weinkauf, Tino
2014-12-01
Data acquisition, numerical inaccuracies, and sampling often introduce noise in measurements and simulations. Removing this noise is often necessary for efficient analysis and visualization of this data, yet many denoising techniques change the minima and maxima of a scalar field. For example, the extrema can appear or disappear, spatially move, and change their value. This can lead to wrong interpretations of the data, e.g., when the maximum temperature over an area is falsely reported being a few degrees cooler because the denoising method is unaware of these features. Recently, a topological denoising technique based on a global energy optimization was proposed, which allows the topology-controlled denoising of 2D scalar fields. While this method preserves the minima and maxima, it is constrained by the size of the data. We extend this work to large 2D data and medium-sized 3D data by introducing a novel domain decomposition approach. It allows processing small patches of the domain independently while still avoiding the introduction of new critical points. Furthermore, we propose an iterative refinement of the solution, which decreases the optimization energy compared to the previous approach and therefore gives smoother results that are closer to the input. We illustrate our technique on synthetic and real-world 2D and 3D data sets that highlight potential applications. PMID:26356972
Dynamic UltraFast 2D EXchange SpectroscopY (UF-EXSY) of hyperpolarized substrates
Swisher, Christine Leon; Koelsch, Bertram; Sukumar, Subramianam; Sriram, Renuka; Santos, Romelyn Delos; Wang, Zhen Jane; Kurhanewicz, John; Vigneron, Daniel; Larson, Peder
2015-01-01
In this work, we present a new ultrafast method for acquiring dynamic 2D EXchange SpectroscopY (EXSY) within a single acquisition. This technique reconstructs two-dimensional EXSY spectra from one-dimensional spectra based on the phase accrual during echo times. The Ultrafast-EXSY acquisition overcomes long acquisition times typically needed to acquire 2D NMR data by utilizing sparsity and phase dependence to dramatically undersample in the indirect time dimension. This allows for the acquisition of the 2D spectrum within a single shot. We have validated this method in simulations and hyperpolarized enzyme assay experiments separating the dehydration of pyruvate and lactate-to-pyruvate conversion. In a renal cell carcinoma cell (RCC) line, bidirectional exchange was observed. This new technique revealed decreased conversion of lactate-to-pyruvate with high expression of monocarboxylate transporter 4 (MCT4), known to correlate with aggressive cancer phenotypes. We also showed feasibility of this technique in vivo in a RCC model where bidirectional exchange was observed for pyruvate–lactate, pyruvate–alanine, and pyruvate–hydrate and were resolved in time. Broadly, the technique is well suited to investigate the dynamics of multiple exchange pathways and applicable to hyperpolarized substrates where chemical exchange has shown great promise across a range of disciplines. PMID:26117655
Progressive attenuation fields: Fast 2D-3D image registration without precomputation
Rohlfing, Torsten; Russakoff, Daniel B.; Denzler, Joachim; Mori, Kensaku; Maurer, Calvin R. Jr.
2005-09-15
Computation of digitally reconstructed radiograph (DRR) images is the rate-limiting step in most current intensity-based algorithms for the registration of three-dimensional (3D) images to two-dimensional (2D) projection images. This paper introduces and evaluates the progressive attenuation field (PAF), which is a new method to speed up DRR computation. A PAF is closely related to an attenuation field (AF). A major difference is that a PAF is constructed on the fly as the registration proceeds; it does not require any precomputation time, nor does it make any prior assumptions of the patient pose or limit the permissible range of patient motion. A PAF effectively acts as a cache memory for projection values once they are computed, rather than as a lookup table for precomputed projections like standard AFs. We use a cylindrical attenuation field parametrization, which is better suited for many medical applications of 2D-3D registration than the usual two-plane parametrization. The computed attenuation values are stored in a hash table for time-efficient storage and access. Using clinical gold-standard spine image data sets from five patients, we demonstrate consistent speedups of intensity-based 2D-3D image registration using PAF DRRs by a factor of 10 over conventional ray casting DRRs with no decrease of registration accuracy or robustness.
A fast non-Fourier method for Landau-fluid operators
Dimits, A. M. Joseph, I.; Umansky, M. V.
2014-05-15
An efficient and versatile non-Fourier method for the computation of Landau-fluid (LF) closure operators [Hammett and Perkins, Phys. Rev. Lett. 64, 3019 (1990)] is presented, based on an approximation by a sum of modified-Helmholtz-equation solves (SMHS) in configuration space. This method can yield fast-Fourier-like scaling of the computational time requirements and also provides a very compact data representation of these operators, even for plasmas with large spatial nonuniformity. As a result, the method can give significant savings compared with direct application of “delocalization kernels” [e.g., Schurtz et al., Phys. Plasmas 7, 4238 (2000)], both in terms of computational cost and memory requirements. The method is of interest for the implementation of Landau-fluid models in situations where the spatial nonuniformity, particular geometry, or boundary conditions render a Fourier implementation difficult or impossible. Systematic procedures have been developed to optimize the resulting operators for accuracy and computational cost. The four-moment Landau-fluid model of Hammett and Perkins has been implemented in the BOUT++ code using the SMHS method for LF closure. Excellent agreement has been obtained for the one-dimensional plasma density response function between driven initial-value calculations using this BOUT++ implementation and matrix eigenvalue calculations using both Fourier and SMHS non-Fourier implementations of the LF closures. The SMHS method also forms the basis for the implementation, which has been carried out in the BOUT++ code, of the parallel and toroidal drift-resonance LF closures. The method is a key enabling tool for the extension of gyro-Landau-fluid models [e.g., Beer and Hammett, Phys. Plasmas 3, 4046 (1996)] to codes that treat regions with strong profile variation, such as the tokamak edge and scrapeoff-layer.
Radar cross-section reduction based on an iterative fast Fourier transform optimized metasurface
NASA Astrophysics Data System (ADS)
Song, Yi-Chuan; Ding, Jun; Guo, Chen-Jiang; Ren, Yu-Hui; Zhang, Jia-Kai
2016-07-01
A novel polarization insensitive metasurface with over 25 dB monostatic radar cross-section (RCS) reduction is introduced. The proposed metasurface is comprised of carefully arranged unit cells with spatially varied dimension, which enables approximate uniform diffusion of incoming electromagnetic (EM) energy and reduces the threat from bistatic radar system. An iterative fast Fourier transform (FFT) method for conventional antenna array pattern synthesis is innovatively applied to find the best unit cell geometry parameter arrangement. Finally, a metasurface sample is fabricated and tested to validate RCS reduction behavior predicted by full wave simulation software Ansys HFSSTM and marvelous agreement is observed.
Fast Acceleration of 2D Wave Propagation Simulations Using Modern Computational Accelerators
Wang, Wei; Xu, Lifan; Cavazos, John; Huang, Howie H.; Kay, Matthew
2014-01-01
Recent developments in modern computational accelerators like Graphics Processing Units (GPUs) and coprocessors provide great opportunities for making scientific applications run faster than ever before. However, efficient parallelization of scientific code using new programming tools like CUDA requires a high level of expertise that is not available to many scientists. This, plus the fact that parallelized code is usually not portable to different architectures, creates major challenges for exploiting the full capabilities of modern computational accelerators. In this work, we sought to overcome these challenges by studying how to achieve both automated parallelization using OpenACC and enhanced portability using OpenCL. We applied our parallelization schemes using GPUs as well as Intel Many Integrated Core (MIC) coprocessor to reduce the run time of wave propagation simulations. We used a well-established 2D cardiac action potential model as a specific case-study. To the best of our knowledge, we are the first to study auto-parallelization of 2D cardiac wave propagation simulations using OpenACC. Our results identify several approaches that provide substantial speedups. The OpenACC-generated GPU code achieved more than speedup above the sequential implementation and required the addition of only a few OpenACC pragmas to the code. An OpenCL implementation provided speedups on GPUs of at least faster than the sequential implementation and faster than a parallelized OpenMP implementation. An implementation of OpenMP on Intel MIC coprocessor provided speedups of with only a few code changes to the sequential implementation. We highlight that OpenACC provides an automatic, efficient, and portable approach to achieve parallelization of 2D cardiac wave simulations on GPUs. Our approach of using OpenACC, OpenCL, and OpenMP to parallelize this particular model on modern computational accelerators should be applicable to other computational models of wave propagation in
A fast partial Fourier transform (FPFT) for data compression and filtering.
Smith, Mark William
2010-07-01
A discrete Fourier transform (DFT) or the closely related discrete cosine transform (DCT) is often employed as part of a data compression scheme. This paper presents a fast partial Fourier transform (FPFT) algorithm that is useful for calculating a subset of M Fourier transform coefficients for a data set comprised of N points (M < N). This algorithm reduces to the standard DFT when M = 1 and it reduces to the radix-2, decimation-in-time FFT when M = N and N is a power of 2. The DFT requires on the order of MN complex floating point multiplications to calculate M coefficients for N data points, a complete FFT requires on the order of (N/2)log{sub 2}N multiplications independent of M, and the new FPFT algorithm requires on the order of (N/2)log{sub 2}M + N multiplications. The FPFT algorithm introduced in this paper could be readily adapted to parallel processing. In addition to data compression, the FPFT algorithm described in this paper might be useful for very narrow band filter operations that pass only a small number of non-zero frequency coefficients such that M << N.
Fast 2D DOA Estimation Algorithm by an Array Manifold Matching Method with Parallel Linear Arrays.
Yang, Lisheng; Liu, Sheng; Li, Dong; Jiang, Qingping; Cao, Hailin
2016-01-01
In this paper, the problem of two-dimensional (2D) direction-of-arrival (DOA) estimation with parallel linear arrays is addressed. Two array manifold matching (AMM) approaches, in this work, are developed for the incoherent and coherent signals, respectively. The proposed AMM methods estimate the azimuth angle only with the assumption that the elevation angles are known or estimated. The proposed methods are time efficient since they do not require eigenvalue decomposition (EVD) or peak searching. In addition, the complexity analysis shows the proposed AMM approaches have lower computational complexity than many current state-of-the-art algorithms. The estimated azimuth angles produced by the AMM approaches are automatically paired with the elevation angles. More importantly, for estimating the azimuth angles of coherent signals, the aperture loss issue is avoided since a decorrelation procedure is not required for the proposed AMM method. Numerical studies demonstrate the effectiveness of the proposed approaches. PMID:26907301
A fast and accurate method to predict 2D and 3D aerodynamic boundary layer flows
NASA Astrophysics Data System (ADS)
Bijleveld, H. A.; Veldman, A. E. P.
2014-12-01
A quasi-simultaneous interaction method is applied to predict 2D and 3D aerodynamic flows. This method is suitable for offshore wind turbine design software as it is a very accurate and computationally reasonably cheap method. This study shows the results for a NACA 0012 airfoil. The two applied solvers converge to the experimental values when the grid is refined. We also show that in separation the eigenvalues remain positive thus avoiding the Goldstein singularity at separation. In 3D we show a flow over a dent in which separation occurs. A rotating flat plat is used to show the applicability of the method for rotating flows. The shown capabilities of the method indicate that the quasi-simultaneous interaction method is suitable for design methods for offshore wind turbine blades.
Fast 2D FWI on a multi and many-cores workstation.
NASA Astrophysics Data System (ADS)
Thierry, Philippe; Donno, Daniela; Noble, Mark
2014-05-01
Following the introduction of x86 co-processors (Xeon Phi) and the performance increase of standard 2-socket workstations using the latest 12 cores E5-v2 x86-64 CPU, we present here a MPI + OpenMP implementation of an acoustic 2D FWI (full waveform inversion) code which simultaneously runs on the CPUs and on the co-processors installed in a workstation. The main advantage of running a 2D FWI on a workstation is to be able to quickly evaluate new features such as more complicated wave equations, new cost functions, finite-difference stencils or boundary conditions. Since the co-processor is made of 61 in-order x86 cores, each of them having up to 4 threads, this many-core can be seen as a shared memory SMP (symmetric multiprocessing) machine with its own IP address. Depending on the vendor, a single workstation can handle several co-processors making the workstation as a personal cluster under the desk. The original Fortran 90 CPU version of the 2D FWI code is just recompiled to get a Xeon Phi x86 binary. This multi and many-core configuration uses standard compilers and associated MPI as well as math libraries under Linux; therefore, the cost of code development remains constant, while improving computation time. We choose to implement the code with the so-called symmetric mode to fully use the capacity of the workstation, but we also evaluate the scalability of the code in native mode (i.e running only on the co-processor) thanks to the Linux ssh and NFS capabilities. Usual care of optimization and SIMD vectorization is used to ensure optimal performances, and to analyze the application performances and bottlenecks on both platforms. The 2D FWI implementation uses finite-difference time-domain forward modeling and a quasi-Newton (with L-BFGS algorithm) optimization scheme for the model parameters update. Parallelization is achieved through standard MPI shot gathers distribution and OpenMP for domain decomposition within the co-processor. Taking advantage of the 16
Fast 2D DOA Estimation Algorithm by an Array Manifold Matching Method with Parallel Linear Arrays
Yang, Lisheng; Liu, Sheng; Li, Dong; Jiang, Qingping; Cao, Hailin
2016-01-01
In this paper, the problem of two-dimensional (2D) direction-of-arrival (DOA) estimation with parallel linear arrays is addressed. Two array manifold matching (AMM) approaches, in this work, are developed for the incoherent and coherent signals, respectively. The proposed AMM methods estimate the azimuth angle only with the assumption that the elevation angles are known or estimated. The proposed methods are time efficient since they do not require eigenvalue decomposition (EVD) or peak searching. In addition, the complexity analysis shows the proposed AMM approaches have lower computational complexity than many current state-of-the-art algorithms. The estimated azimuth angles produced by the AMM approaches are automatically paired with the elevation angles. More importantly, for estimating the azimuth angles of coherent signals, the aperture loss issue is avoided since a decorrelation procedure is not required for the proposed AMM method. Numerical studies demonstrate the effectiveness of the proposed approaches. PMID:26907301
Chaotically Spiking Canards in an Excitable System with 2D Inertial Fast Manifolds
NASA Astrophysics Data System (ADS)
Marino, Francesco; Marin, Francesco; Balle, Salvador; Piro, Oreste
2007-02-01
We introduce a new class of excitable systems with two-dimensional fast dynamics that includes inertia. A novel transition from excitability to relaxation oscillations is discovered where the usual Hopf bifurcation is followed by a cascade of period doubled and chaotic small excitable attractors and, as they grow, by a new type of canard explosion where a small chaotic background erratically but deterministically triggers excitable spikes. This scenario is also found in a model for a nonlinear Fabry-Perot cavity with one pendular mirror.
Katoh, Marcus Opitz, Armin; Minko, Peter; Massmann, Alexander; Berlich, Joachim; Buecker, Arno
2011-06-15
Purpose: To investigate the value of rotational digital subtraction angiography (rDSA) for evaluation of peripheral and visceral artery stenoses compared to conventional digital subtraction angiography (cDSA). Methods: A phantom study was performed comparing the radiation dose of cDSA with two projections and rDSA by means of the 2D Dynavision technique (Siemens Medical Solutions, Forchheim, Germany). Subsequently, 33 consecutive patients (18 women, 15 men; mean {+-} SD age 67 {+-} 15 years) were examined by both techniques. In total, 63 vessel segments were analyzed by two observers with respect to stenoses, image contrast, and vessel sharpness. Results: Radiation dose was significantly lower with rDSA. cDSA and rDSA revealed 21 and 24 flow-relevant stenotic lesions and vessel occlusions (70-100%), respectively. The same stenosis grade was assessed in 45 segments. By means of rDSA, 10 lesions were judged to have a higher and 8 lesions a lower stenosis grade compared to cDSA. rDSA yielded additive information regarding the vessel anatomy and pathology in 29 segments. However, a tendency toward better image quality and sharper vessel visualization was seen with cDSA. Conclusion: rDSA allows for multiprojection assessment of peripheral and visceral arteries and provides additional clinically relevant information after a single bolus of contrast medium. At the same time, radiation dose can be significantly reduced compared to cDSA.
Fast DRR generation for 2D to 3D registration on GPUs
Tornai, Gabor Janos; Cserey, Gyoergy
2012-08-15
Purpose: The generation of digitally reconstructed radiographs (DRRs) is the most time consuming step on the CPU in intensity based two-dimensional x-ray to three-dimensional (CT or 3D rotational x-ray) medical image registration, which has application in several image guided interventions. This work presents optimized DRR rendering on graphical processor units (GPUs) and compares performance achievable on four commercially available devices. Methods: A ray-cast based DRR rendering was implemented for a 512 Multiplication-Sign 512 Multiplication-Sign 72 CT volume. The block size parameter was optimized for four different GPUs for a region of interest (ROI) of 400 Multiplication-Sign 225 pixels with different sampling ratios (1.1%-9.1% and 100%). Performance was statistically evaluated and compared for the four GPUs. The method and the block size dependence were validated on the latest GPU for several parameter settings with a public gold standard dataset (512 Multiplication-Sign 512 Multiplication-Sign 825 CT) for registration purposes. Results: Depending on the GPU, the full ROI is rendered in 2.7-5.2 ms. If sampling ratio of 1.1%-9.1% is applied, execution time is in the range of 0.3-7.3 ms. On all GPUs, the mean of the execution time increased linearly with respect to the number of pixels if sampling was used. Conclusions: The presented results outperform other results from the literature. This indicates that automatic 2D to 3D registration, which typically requires a couple of hundred DRR renderings to converge, can be performed quasi on-line, in less than a second or depending on the application and hardware in less than a couple of seconds. Accordingly, a whole new field of applications is opened for image guided interventions, where the registration is continuously performed to match the real-time x-ray.
Hill, N.C.; Limbach, P.A.; Shomo, R.E. II; Marshall, A.G. ); Appelhans, A.D.; Delmore, J.E. )
1991-11-01
The coupling of an autoneutralizing SF{sup {minus}}{sub 6} fast ion-beam gun to a Fourier transform ion cyclotron resonance (FT/ICR) mass spectrometer is described. The fast neutral beam provides for secondary-ion-type FT/ICR mass analysis (e.g., production of abundant pseudomolecular (M+H){sup +} ions) of involatile samples without the need for external ion injection, since ions are formed at the entrance to the ICR ion trap. The design, construction, and testing of the hybrid instrument are described. The feasibility of the experiment (for both broadband and high-resolution FT/ICR positive-ion mass spectra) is demonstrated with {ital tetra}-butylammonium bromide and a Tylenol{sup ( )} sample. The ability to analyze high molecular weight polymers with high mass resolution is demonstrated for Teflon{sup ( )}. All of the advantages of the fast neutral beam ion source previously demonstrated with quadrupole mass analysis are preserved, and the additional advantages of FT/ICR mass analysis (e.g., high mass resolving power, ion trapping) are retained.
NASA Astrophysics Data System (ADS)
Hill, Nicholas C.; Limbach, Patrick A.; Shomo, Ronald E., II; Marshall, Alan G.; Appelhans, Anthony D.; Delmore, James E.
1991-11-01
The coupling of an autoneutralizing SF-6 fast ion-beam gun to a Fourier transform ion cyclotron resonance (FT/ICR) mass spectrometer is described. The fast neutral beam provides for secondary-ion-type FT/ICR mass analysis [e.g., production of abundant pseudomolecular (M+H)+ ions] of involatile samples without the need for external ion injection, since ions are formed at the entrance to the ICR ion trap. The design, construction, and testing of the hybrid instrument are described. The feasibility of the experiment (for both broadband and high-resolution FT/ICR positive-ion mass spectra) is demonstrated with tetra-butylammonium bromide and a Tylenol■ sample. The ability to analyze high molecular weight polymers with high mass resolution is demonstrated for Teflon■. All of the advantages of the fast neutral beam ion source previously demonstrated with quadrupole mass analysis are preserved, and the additional advantages of FT/ICR mass analysis (e.g., high mass resolving power, ion trapping) are retained.
Liu, Derek Sloboda, Ron S.
2014-05-15
Purpose: Boyer and Mok proposed a fast calculation method employing the Fourier transform (FT), for which calculation time is independent of the number of seeds but seed placement is restricted to calculation grid points. Here an interpolation method is described enabling unrestricted seed placement while preserving the computational efficiency of the original method. Methods: The Iodine-125 seed dose kernel was sampled and selected values were modified to optimize interpolation accuracy for clinically relevant doses. For each seed, the kernel was shifted to the nearest grid point via convolution with a unit impulse, implemented in the Fourier domain. The remaining fractional shift was performed using a piecewise third-order Lagrange filter. Results: Implementation of the interpolation method greatly improved FT-based dose calculation accuracy. The dose distribution was accurate to within 2% beyond 3 mm from each seed. Isodose contours were indistinguishable from explicit TG-43 calculation. Dose-volume metric errors were negligible. Computation time for the FT interpolation method was essentially the same as Boyer's method. Conclusions: A FT interpolation method for permanent prostate brachytherapy TG-43 dose calculation was developed which expands upon Boyer's original method and enables unrestricted seed placement. The proposed method substantially improves the clinically relevant dose accuracy with negligible additional computation cost, preserving the efficiency of the original method.
Peak velocity determination using fast Fourier velocity encoding with minimal spatial encoding.
Galea, Daniela; Lauzon, M Louis; Drangova, Maria
2002-08-01
For quantitative peak velocity determination, a technique was developed that uses Fourier velocity encoding (FVE) for the fast acquisition of images of velocity with no spatial encoding other than slice selection. The technique produces images of velocity versus temporal frequency. In applications where the quantity of interest is the peak velocity and in-plane spatial localization is not required, high SNR images are produced with reduced sensitivity to errors due to slice thickness and motion. The technique was validated using steady and pulsatile flow in a straight tube, and compared to both phase contrast measurements and numerical models using steady flow in a 50% and a 75% cosinusoidal stenosis phantom. Results show that for slices as large as 2 cm and/or undergoing periodic motion, FVE can accurately measure the peak velocity in cases where a distribution of velocities exist. PMID:12201419
NASA Astrophysics Data System (ADS)
Raynolds, James E.; Mullin, Lenore R.
2005-07-01
The techniques of Conformal Computing are introduced with an application to the Fast Fourier Transform. Conformal Computing is a design methodology, based on a rigorous mathematical foundation, which provides a systematic approach to the most efficient organization of all levels of the software and hardware design hierarchy from high-level software constructs all the way down to the design of the integrated circuits. We show that using these general design principles, without any specialized optimization, leads to portable, scalable, code that is competitive with other well-tuned machine specific routines. Further improvements are straightforward within our formalism by taking into account specific hardware details (e.g., cache loops) in a portable parametric way. We also argue that the present theory constitutes a uniform way of reasoning about physics and the data structures that define physics on computers.
Nanowire humidity optical sensor system based on fast Fourier transform technique
NASA Astrophysics Data System (ADS)
Rota-Rodrigo, S.; Pérez-Herrera, R.; Lopez-Aldaba, A.; López Bautista, M. C.; Esteban, O.; López-Amo, M.
2015-09-01
In this paper, a new sensor system for relative humidity measurements based on its interaction with the evanescent field of a nanowire is presented. The interrogation of the sensing head is carried out by monitoring the fast Fourier transform phase variations of one of the nanowire interference frequencies. This method is independent of the signal amplitude and also avoids the necessity of tracking the wavelength evolution in the spectrum, which can be a handicap when there are multiple interference frequency components with different sensitivities. The sensor is operated within a wide humidity range (20%-70% relative humidity) with a maximum sensitivity achieved of 0.14rad/% relative humidity. Finally, due to the system uses an optical interrogator as unique active element, the system presents a cost-effective feature.
A general purpose subroutine for fast fourier transform on a distributed memory parallel machine
NASA Technical Reports Server (NTRS)
Dubey, A.; Zubair, M.; Grosch, C. E.
1992-01-01
One issue which is central in developing a general purpose Fast Fourier Transform (FFT) subroutine on a distributed memory parallel machine is the data distribution. It is possible that different users would like to use the FFT routine with different data distributions. Thus, there is a need to design FFT schemes on distributed memory parallel machines which can support a variety of data distributions. An FFT implementation on a distributed memory parallel machine which works for a number of data distributions commonly encountered in scientific applications is presented. The problem of rearranging the data after computing the FFT is also addressed. The performance of the implementation on a distributed memory parallel machine Intel iPSC/860 is evaluated.
NASA Astrophysics Data System (ADS)
Liu, B.; Raabe, D.; Roters, F.; Eisenlohr, P.; Lebensohn, R. A.
2010-12-01
We compare two full-field formulations, i.e. a crystal plasticity fast Fourier transform-based (CPFFT) model and the crystal plasticity finite element model (CPFEM) in terms of the deformation textures predicted by both approaches. Plane-strain compression of a 1024-grain ensemble is simulated with CPFFT and CPFEM to assess the models in terms of their predictions of texture evolution for engineering applications. Different combinations of final textures and strain distributions are obtained with the CPFFT and CPFEM models for this 1024-grain polycrystal. To further understand these different predictions, the correlation between grain rotations and strain gradients is investigated through the simulation of plane-strain compression of bicrystals. Finally, a study of the influence of the initial crystal orientation and the crystallographic neighborhood on grain rotations and grain subdivisions is carried out by means of plane-strain compression simulations of a 64-grain cluster.
Fast Fourier transform analysis of sounds made while swallowing various foods.
Taniwaki, Mitsuru; Kohyama, Kaoru
2012-10-01
The cervical auscultation method was applied to investigate sounds generated while swallowing various foods with unique physical properties, including liquid (water), semiliquid (yogurt), and solid (konjac jelly). To study the differences among swallowing sounds for various foods, fast Fourier transform (FFT) analysis was applied to signals that were attributed to the flow of a food bolus, which is a swallowable soft mass of chewed food. An FFT program was developed that enabled the calculation of a spectrum for a specified region of time domain swallowing sound signals. The intensity of spectra in the frequency range between 400 and 1000 Hz significantly differed: liquid > semiliquid > solid. The FFT spectrum in this range was suggested to represent the frequency characteristics of the swallowing sounds of various foods. PMID:23039442
Murali, S; Vladimir, Kulish V
2007-10-01
Fast Fourier transform (FFT) and fractal dimension (FD) are computed to discover the evoked potentials (EPs) in the base human electroencephalograms (EEGs) corresponding to six types of odor inhalation, one at a time. During the first 20 s of the EEG recording, the patient breathes normally and the rest of the recording time 21-64 s, inhales an odor, when both the eye open and closed modes. FFT spectral analysis cannot discriminate the potential due to olfactory stimuli from the base EEG. However, fractal spectra analysis predicts the evoked potential due to inhalation of English rose, jasmine, lemon, orange, rose garden, and sandalwood concentrates. Similar results are obtained for three of the patients examined and on computing for five EEG recordings of two patients during inhalation of six types of odor. PMID:17729151
Directional 2D functions as models for fast layout pattern transfer verification
NASA Astrophysics Data System (ADS)
Torres, J. Andres; Hofmann, Mark; Otto, Oberdan
2009-03-01
As advanced manufacturing processes become more stable, the need to adapt new designs to fully utilize the available manufacturing technology becomes a key technologic differentiator. However, many times such gains can only be realized and evaluated during full chip analysis. It has been demonstrated that the most accurate layout verification methods require application of the actual OPC recipes along with most of the mask data preparation that defines the pattern transfer characteristics of the process. Still, this method in many instances is not sufficiently fast to be used in a layout creation environment which undergoes constant updates. By doing an analysis of typical mask data processing, it is possible to determine that the most CPUintensive computations are the OPC and contour simulation steps needed to perform layout printability checks. Several researchers have tried to reduce the time it takes to compute the OPC mask by introducing matrix convolutions of the layout with empirically calibrated two-dimensional functions. However, most of these approaches do not provide a sufficient speed-up since they only replace the OPC computation and still require a full contour computation. Another alternative is to try to find effective ways of pattern matching those topologies that will exhibit transfer difficulties4, but such methods lack the ability to be predictive beyond their calibration data. In this paper we present a methodology that includes common resolution enhancement techniques, such as retargeting and sub-resolution assist feature insertion, and which replaces the OPC computation and subsequent contour calculation with an edge bias function based on an empirically-calibrated, directional, two-dimensional function. Because the edge bias function does not provide adequate control over the corner locations, a spline-based smoothing process is applied. The outcome is a piecewise-linear curve similar to those obtained by full lithographic simulations. Our
Sasaoka, Sara; Saito, Koichi; Higashi, Kenjirou; Limwikrant, Waree; Moribe, Kunikazu; Suzuki, Shinichi; Yamamoto, Keiji
2015-12-01
A novel type of spectrum, the one-dimensional power spectrum (1D-PS), was designed for the discrimination of adhesive packing tapes, i.e., kraft tapes. The 1D-PS offered complementary information to that provided by the improved two-dimensional PS (2D-PS), which was calculated using our previously established image processes combined with a two-dimensional fast Fourier transform (2D-FFT) to obtain information about the spatial periodicity within kraft tapes. The 1D-PS was calculated using a three-step image process: (i) the 2D-FFT was applied to 50 randomly selected areas in a transmitted light image; (ii) the obtained 2D-PSs were accumulated without applying a logarithmic transform; (iii) the wavenumber and the maximum intensity were plotted on the x-axis and y-axis, respectively. Through an intra-roll comparison, the 1D-PSs collected from single rolls showed similar profiles. In an inter-roll comparison, the 1D-PSs from 50 commercially available brand-name products were classified into 29 groups. The 1D-PSs contained other useful information than that provided by the improved 2D-PSs: they presented more peaks and absolute intensity with a wider range. The 1D-PSs enabled us to compare the spectra quickly and easily, owing to their unchanging profiles regardless of the orientation of the scanned images. A combined use of the 1D-PSs with the improved 2D-PSs-both spectrum types being convenient, rapid, non-destructive, and applicable to dirty and/or damaged samples-could further improve the identification of kraft tapes. PMID:26461031
Robust 2D/3D registration for fast-flexion motion of the knee joint using hybrid optimization.
Ohnishi, Takashi; Suzuki, Masahiko; Kobayashi, Tatsuya; Naomoto, Shinji; Sukegawa, Tomoyuki; Nawata, Atsushi; Haneishi, Hideaki
2013-01-01
Previously, we proposed a 2D/3D registration method that uses Powell's algorithm to obtain 3D motion of a knee joint by 3D computed-tomography and bi-plane fluoroscopic images. The 2D/3D registration is performed consecutively and automatically for each frame of the fluoroscopic images. This method starts from the optimum parameters of the previous frame for each frame except for the first one, and it searches for the next set of optimum parameters using Powell's algorithm. However, if the flexion motion of the knee joint is fast, it is likely that Powell's algorithm will provide a mismatch because the initial parameters are far from the correct ones. In this study, we applied a hybrid optimization algorithm (HPS) combining Powell's algorithm with the Nelder-Mead simplex (NM-simplex) algorithm to overcome this problem. The performance of the HPS was compared with the separate performances of Powell's algorithm and the NM-simplex algorithm, the Quasi-Newton algorithm and hybrid optimization algorithm with the Quasi-Newton and NM-simplex algorithms with five patient data sets in terms of the root-mean-square error (RMSE), target registration error (TRE), success rate, and processing time. The RMSE, TRE, and the success rate of the HPS were better than those of the other optimization algorithms, and the processing time was similar to that of Powell's algorithm alone. PMID:23138929
NASA Astrophysics Data System (ADS)
Sotiropoulou, C. L.; Gkaitatzis, S.; Annovi, A.; Beretta, M.; Kordas, K.; Nikolaidis, S.; Petridou, C.; Volpi, G.
2014-10-01
The parallel 2D pixel clustering FPGA implementation used for the input system of the ATLAS Fast TracKer (FTK) processor is presented. The input system for the FTK processor will receive data from the Pixel and micro-strip detectors from inner ATLAS read out drivers (RODs) at full rate, for total of 760Gbs, as sent by the RODs after level-1 triggers. Clustering serves two purposes, the first is to reduce the high rate of the received data before further processing, the second is to determine the cluster centroid to obtain the best spatial measurement. For the pixel detectors the clustering is implemented by using a 2D-clustering algorithm that takes advantage of a moving window technique to minimize the logic required for cluster identification. The cluster detection window size can be adjusted for optimizing the cluster identification process. Additionally, the implementation can be parallelized by instantiating multiple cores to identify different clusters independently thus exploiting more FPGA resources. This flexibility makes the implementation suitable for a variety of demanding image processing applications. The implementation is robust against bit errors in the input data stream and drops all data that cannot be identified. In the unlikely event of missing control words, the implementation will ensure stable data processing by inserting the missing control words in the data stream. The 2D pixel clustering implementation is developed and tested in both single flow and parallel versions. The first parallel version with 16 parallel cluster identification engines is presented. The input data from the RODs are received through S-Links and the processing units that follow the clustering implementation also require a single data stream, therefore data parallelizing (demultiplexing) and serializing (multiplexing) modules are introduced in order to accommodate the parallelized version and restore the data stream afterwards. The results of the first hardware tests of
NASA Astrophysics Data System (ADS)
McLaughlin, J. A.; De Moortel, I.; Hood, A. W.; Brady, C. S.
2009-01-01
Context: This paper extends the models of Craig & McClymont (1991, ApJ, 371, L41) and McLaughlin & Hood (2004, A&A, 420, 1129) to include finite β and nonlinear effects. Aims: We investigate the nature of nonlinear fast magnetoacoustic waves about a 2D magnetic X-point. Methods: We solve the compressible and resistive MHD equations using a Lagrangian remap, shock capturing code (Arber et al. 2001, J. Comp. Phys., 171, 151) and consider an initial condition in {v}×{B} \\cdot {hat{z}} (a natural variable of the system). Results: We observe the formation of both fast and slow oblique magnetic shocks. The nonlinear wave deforms the X-point into a “cusp-like” point which in turn collapses to a current sheet. The system then evolves through a series of horizontal and vertical current sheets, with associated changes in connectivity, i.e. the system exhibits oscillatory reconnection. Our final state is non-potential (but in force balance) due to asymmetric heating from the shocks. Larger amplitudes in our initial condition correspond to larger values of the final current density left in the system. Conclusions: The inclusion of nonlinear terms introduces several new features to the system that were absent from the linear regime. A movie is available in electronic form at http://www.aanda.org
Protein-protein docking by fast generalized Fourier transforms on 5D rotational manifolds.
Padhorny, Dzmitry; Kazennov, Andrey; Zerbe, Brandon S; Porter, Kathryn A; Xia, Bing; Mottarella, Scott E; Kholodov, Yaroslav; Ritchie, David W; Vajda, Sandor; Kozakov, Dima
2016-07-26
Energy evaluation using fast Fourier transforms (FFTs) enables sampling billions of putative complex structures and hence revolutionized rigid protein-protein docking. However, in current methods, efficient acceleration is achieved only in either the translational or the rotational subspace. Developing an efficient and accurate docking method that expands FFT-based sampling to five rotational coordinates is an extensively studied but still unsolved problem. The algorithm presented here retains the accuracy of earlier methods but yields at least 10-fold speedup. The improvement is due to two innovations. First, the search space is treated as the product manifold [Formula: see text], where [Formula: see text] is the rotation group representing the space of the rotating ligand, and [Formula: see text] is the space spanned by the two Euler angles that define the orientation of the vector from the center of the fixed receptor toward the center of the ligand. This representation enables the use of efficient FFT methods developed for [Formula: see text] Second, we select the centers of highly populated clusters of docked structures, rather than the lowest energy conformations, as predictions of the complex, and hence there is no need for very high accuracy in energy evaluation. Therefore, it is sufficient to use a limited number of spherical basis functions in the Fourier space, which increases the efficiency of sampling while retaining the accuracy of docking results. A major advantage of the method is that, in contrast to classical approaches, increasing the number of correlation function terms is computationally inexpensive, which enables using complex energy functions for scoring. PMID:27412858
AN OPTIMIZED 64X64 POINT TWO-DIMENSIONAL FAST FOURIER TRANSFORM
NASA Technical Reports Server (NTRS)
Miko, J.
1994-01-01
Scientists at Goddard have developed an efficient and powerful program-- An Optimized 64x64 Point Two-Dimensional Fast Fourier Transform-- which combines the performance of real and complex valued one-dimensional Fast Fourier Transforms (FFT's) to execute a two-dimensional FFT and its power spectrum coefficients. These coefficients can be used in many applications, including spectrum analysis, convolution, digital filtering, image processing, and data compression. The program's efficiency results from its technique of expanding all arithmetic operations within one 64-point FFT; its high processing rate results from its operation on a high-speed digital signal processor. For non-real-time analysis, the program requires as input an ASCII data file of 64x64 (4096) real valued data points. As output, this analysis produces an ASCII data file of 64x64 power spectrum coefficients. To generate these coefficients, the program employs a row-column decomposition technique. First, it performs a radix-4 one-dimensional FFT on each row of input, producing complex valued results. Then, it performs a one-dimensional FFT on each column of these results to produce complex valued two-dimensional FFT results. Finally, the program sums the squares of the real and imaginary values to generate the power spectrum coefficients. The program requires a Banshee accelerator board with 128K bytes of memory from Atlanta Signal Processors (404/892-7265) installed on an IBM PC/AT compatible computer (DOS ver. 3.0 or higher) with at least one 16-bit expansion slot. For real-time operation, an ASPI daughter board is also needed. The real-time configuration reads 16-bit integer input data directly into the accelerator board, operating on 64x64 point frames of data. The program's memory management also allows accumulation of the coefficient results. The real-time processing rate to calculate and accumulate the 64x64 power spectrum output coefficients is less than 17.0 mSec. Documentation is included
Diffuse correlation spectroscopy with a fast Fourier transform-based software autocorrelator
NASA Astrophysics Data System (ADS)
Dong, Jing; Bi, Renzhe; Ho, Jun Hui; Thong, Patricia S. P.; Soo, Khee-Chee; Lee, Kijoon
2012-09-01
Diffuse correlation spectroscopy (DCS) is an emerging noninvasive technique that probes the deep tissue blood flow, by using the time-averaged intensity autocorrelation function of the fluctuating diffuse reflectance signal. We present a fast Fourier transform (FFT)-based software autocorrelator that utilizes the graphical programming language LabVIEW (National Instruments) to complete data acquisition, recording, and processing tasks. The validation and evaluation experiments were conducted on an in-house flow phantom, human forearm, and photodynamic therapy (PDT) on mouse tumors under the acquisition rate of ˜400 kHz. The software autocorrelator in general has certain advantages, such as flexibility in raw photon count data preprocessing and low cost. In addition to that, our FFT-based software autocorrelator offers smoother starting and ending plateaus when compared to a hardware correlator, which could directly benefit the fitting results without too much sacrifice in speed. We show that the blood flow index (BFI) obtained by using a software autocorrelator exhibits better linear behavior in a phantom control experiment when compared to a hardware one. The results indicate that an FFT-based software autocorrelator can be an alternative solution to the conventional hardware ones in DCS systems with considerable benefits.
Big Data in Reciprocal Space: Sliding Fast Fourier Transforms for Determining Periodicity
Vasudevan, Rama K.; Belianinov, Alex; Gianfrancesco, Anthony G.; Baddorf, Arthur P.; Tselev, Alexander; Kalinin, Sergei V.; Jesse, Stephen
2015-03-03
Significant advances in atomically resolved imaging of crystals and surfaces have occurred in the last decade allowing unprecedented insight into local crystal structures and periodicity. Yet, the analysis of the long-range periodicity from the local imaging data, critical to correlation of functional properties and chemistry to the local crystallography, remains a challenge. Here, we introduce a Sliding Fast Fourier Transform (FFT) filter to analyze atomically resolved images of in-situ grown La5/8Ca3/8MnO3 films. We demonstrate the ability of sliding FFT algorithm to differentiate two sub-lattices, resulting from a mixed-terminated surface. Principal Component Analysis (PCA) and Independent Component Analysis (ICA) of themore » Sliding FFT dataset reveal the distinct changes in crystallography, step edges and boundaries between the multiple sub-lattices. The method is universal for images with any periodicity, and is especially amenable to atomically resolved probe and electron-microscopy data for rapid identification of the sub-lattices present.« less
NASA Technical Reports Server (NTRS)
Logan, T. L.; Huning, J. R.; Glackin, D. L.
1983-01-01
The use of two dimensional Fast Fourier Transforms (FFTs) subjected to pattern recognition technology for the identification and classification of low altitude stratus cloud structure from Geostationary Operational Environmental Satellite (GOES) imagery was examined. The development of a scene independent pattern recognition methodology, unconstrained by conventional cloud morphological classifications was emphasized. A technique for extracting cloud shape, direction, and size attributes from GOES visual imagery was developed. These attributes were combined with two statistical attributes (cloud mean brightness, cloud standard deviation), and interrogated using unsupervised clustering amd maximum likelihood classification techniques. Results indicate that: (1) the key cloud discrimination attributes are mean brightness, direction, shape, and minimum size; (2) cloud structure can be differentiated at given pixel scales; (3) cloud type may be identifiable at coarser scales; (4) there are positive indications of scene independence which would permit development of a cloud signature bank; (5) edge enhancement of GOES imagery does not appreciably improve cloud classification over the use of raw data; and (6) the GOES imagery must be apodized before generation of FFTs.
Diffuse correlation spectroscopy with a fast Fourier transform-based software autocorrelator.
Dong, Jing; Bi, Renzhe; Ho, Jun Hui; Thong, Patricia S P; Soo, Khee-Chee; Lee, Kijoon
2012-09-01
Diffuse correlation spectroscopy (DCS) is an emerging noninvasive technique that probes the deep tissue blood flow, by using the time-averaged intensity autocorrelation function of the fluctuating diffuse reflectance signal. We present a fast Fourier transform (FFT)-based software autocorrelator that utilizes the graphical programming language LabVIEW (National Instruments) to complete data acquisition, recording, and processing tasks. The validation and evaluation experiments were conducted on an in-house flow phantom, human forearm, and photodynamic therapy (PDT) on mouse tumors under the acquisition rate of ∼400 kHz. The software autocorrelator in general has certain advantages, such as flexibility in raw photon count data preprocessing and low cost. In addition to that, our FFT-based software autocorrelator offers smoother starting and ending plateaus when compared to a hardware correlator, which could directly benefit the fitting results without too much sacrifice in speed. We show that the blood flow index (BFI) obtained by using a software autocorrelator exhibits better linear behavior in a phantom control experiment when compared to a hardware one. The results indicate that an FFT-based software autocorrelator can be an alternative solution to the conventional hardware ones in DCS systems with considerable benefits. PMID:23085922
Big Data in Reciprocal Space: Sliding Fast Fourier Transforms for Determining Periodicity
Vasudevan, Rama K.; Belianinov, Alex; Gianfrancesco, Anthony G.; Baddorf, Arthur P.; Tselev, Alexander; Kalinin, Sergei V.; Jesse, Stephen
2015-03-03
Significant advances in atomically resolved imaging of crystals and surfaces have occurred in the last decade allowing unprecedented insight into local crystal structures and periodicity. Yet, the analysis of the long-range periodicity from the local imaging data, critical to correlation of functional properties and chemistry to the local crystallography, remains a challenge. Here, we introduce a Sliding Fast Fourier Transform (FFT) filter to analyze atomically resolved images of in-situ grown La5/8Ca3/8MnO3 films. We demonstrate the ability of sliding FFT algorithm to differentiate two sub-lattices, resulting from a mixed-terminated surface. Principal Component Analysis (PCA) and Independent Component Analysis (ICA) of the Sliding FFT dataset reveal the distinct changes in crystallography, step edges and boundaries between the multiple sub-lattices. The method is universal for images with any periodicity, and is especially amenable to atomically resolved probe and electron-microscopy data for rapid identification of the sub-lattices present.
Fast Fourier transform to measure pressure coefficient of muons in the GRAPES-3 experiment
NASA Astrophysics Data System (ADS)
Mohanty, P. K.; Ahmad, S.; Antia, H. M.; Arunbabu, K. P.; Chandra, A.; Dugad, S. R.; Gupta, S. K.; Hariharan, B.; Hayashi, Y.; Jagadeesan, P.; Jain, A.; Kawakami, S.; Kojima, H.; Morris, S. D.; Nayak, P. K.; Oshima, A.; Rao, B. S.; Reddy, L. V.; Shibata, S.
2016-06-01
The GRAPES-3 large area (560 m2) tracking muon telescope is operating at Ooty in India since 2001. It records 4 × 109 muons of energy ≥ 1 GeV every day. These high statistics data have enabled extremely sensitive measurements of solar phenomena, including the solar anisotropies, Forbush decreases, coronal mass ejections etc. to be made. However, prior to such studies, the variation in observed muon rate caused by changes in atmospheric pressure needs to be corrected. Traditionally, the pressure coefficient (β) for the muon rate was derived from the observed data. But the influence of various solar effects makes the measurement of β somewhat difficult. In the present work, a different approach to circumvent this difficulty was used to measure β, almost independent of the solar activity. This approach exploits a small amplitude (∼1 hPa) periodic (12 h) variation of atmospheric pressure at Ooty that introduces a synchronous variation in the muon rate. By using the fast Fourier transform technique the spectral power distributions at 12 h from the atmospheric pressure, and muon rate were used to measure β. The value of pressure coefficient was found to be β =(- 0.128 ± 0.005) % hPa-1.
A Fourier-series-based kernel-independent fast multipole method
Zhang Bo; Huang Jingfang; Pitsianis, Nikos P.; Sun Xiaobai
2011-07-01
We present in this paper a new kernel-independent fast multipole method (FMM), named as FKI-FMM, for pairwise particle interactions with translation-invariant kernel functions. FKI-FMM creates, using numerical techniques, sufficiently accurate and compressive representations of a given kernel function over multi-scale interaction regions in the form of a truncated Fourier series. It provides also economic operators for the multipole-to-multipole, multipole-to-local, and local-to-local translations that are typical and essential in the FMM algorithms. The multipole-to-local translation operator, in particular, is readily diagonal and does not dominate in arithmetic operations. FKI-FMM provides an alternative and competitive option, among other kernel-independent FMM algorithms, for an efficient application of the FMM, especially for applications where the kernel function consists of multi-physics and multi-scale components as those arising in recent studies of biological systems. We present the complexity analysis and demonstrate with experimental results the FKI-FMM performance in accuracy and efficiency.
Iterative Image Reconstruction for PROPELLER-MRI using the NonUniform Fast Fourier Transform
Tamhane, Ashish A.; Anastasio, Mark A.; Gui, Minzhi; Arfanakis, Konstantinos
2013-01-01
Purpose To investigate an iterative image reconstruction algorithm using the non-uniform fast Fourier transform (NUFFT) for PROPELLER (Periodically Rotated Overlapping parallEL Lines with Enhanced Reconstruction) MRI. Materials and Methods Numerical simulations, as well as experiments on a phantom and a healthy human subject were used to evaluate the performance of the iterative image reconstruction algorithm for PROPELLER, and compare it to that of conventional gridding. The trade-off between spatial resolution, signal to noise ratio, and image artifacts, was investigated for different values of the regularization parameter. The performance of the iterative image reconstruction algorithm in the presence of motion was also evaluated. Results It was demonstrated that, for a certain range of values of the regularization parameter, iterative reconstruction produced images with significantly increased SNR, reduced artifacts, for similar spatial resolution, compared to gridding. Furthermore, the ability to reduce the effects of motion in PROPELLER-MRI was maintained when using the iterative reconstruction approach. Conclusion An iterative image reconstruction technique based on the NUFFT was investigated for PROPELLER MRI. For a certain range of values of the regularization parameter the new reconstruction technique may provide PROPELLER images with improved image quality compared to conventional gridding. PMID:20578028
NASA Technical Reports Server (NTRS)
Rider, D.; Blavier, J-F.; Cunningham, T.; Hancock, B.; Key, R.; Pannell, Z.; Sander, S.; Seshadri, S.; Sun, C.; Wrigley, C.
2011-01-01
Focal plane arrays (FPAs) with high frame rates and many pixels benefit several upcoming Earth science missions including GEO-CAPE, GACM, and ACE by enabling broader spatial coverage and higher spectral resolution. FPAs for the PanFTS, a high spatial resolution Fourier transform spectrometer and a candidate instrument for the GEO-CAPE mission are the focus of the developments reported here, but this FPA technology has the potential to enable a variety of future measurements and instruments. The ESTO ACT Program funded the developed of a fast readout integrated circuit (ROIC) based on an innovative in-pixel analog-to-digital converter (ADC). The 128 X 128 pixel ROIC features 60 ?m pixels, a 14-bit ADC in each pixel and operates at a continuous frame rate of 14 kHz consuming only 1.1 W of power. The ROIC outputs digitized data completely eliminating the bulky, power consuming signal chains needed by conventional FPAs. The 128 X 128 pixel ROIC has been fabricated in CMOS and tested at the Jet Propulsion Laboratory. The current version is designed to be hybridized with PIN photodiode arrays via indium bump bonding for light detection in the visible and ultraviolet spectral regions. However, the ROIC design incorporates a small photodiode in each cell to permit detailed characterization of the ROICperformance without the need for hybridization. We will describe the essential features of the ROIC design and present results of ROIC performance measurements.
Support vector machine based classification of fast Fourier transform spectroscopy of proteins
NASA Astrophysics Data System (ADS)
Lazarevic, Aleksandar; Pokrajac, Dragoljub; Marcano, Aristides; Melikechi, Noureddine
2009-02-01
Fast Fourier transform spectroscopy has proved to be a powerful method for study of the secondary structure of proteins since peak positions and their relative amplitude are affected by the number of hydrogen bridges that sustain this secondary structure. However, to our best knowledge, the method has not been used yet for identification of proteins within a complex matrix like a blood sample. The principal reason is the apparent similarity of protein infrared spectra with actual differences usually masked by the solvent contribution and other interactions. In this paper, we propose a novel machine learning based method that uses protein spectra for classification and identification of such proteins within a given sample. The proposed method uses principal component analysis (PCA) to identify most important linear combinations of original spectral components and then employs support vector machine (SVM) classification model applied on such identified combinations to categorize proteins into one of given groups. Our experiments have been performed on the set of four different proteins, namely: Bovine Serum Albumin, Leptin, Insulin-like Growth Factor 2 and Osteopontin. Our proposed method of applying principal component analysis along with support vector machines exhibits excellent classification accuracy when identifying proteins using their infrared spectra.
Big data in reciprocal space: Sliding fast Fourier transforms for determining periodicity
Vasudevan, Rama K. Belianinov, Alex; Baddorf, Arthur P.; Tselev, Alexander; Jesse, S.; Gianfrancesco, Anthony G.
2015-03-02
Significant advances in atomically resolved imaging of crystals and surfaces have occurred in the last decade allowing unprecedented insight into local crystal structures and periodicity. Yet, the analysis of the long-range periodicity from the local imaging data, critical to correlation of functional properties and chemistry to the local crystallography, remains a challenge. Here, we introduce a Sliding Fast Fourier Transform (FFT) filter to analyze atomically resolved images of in-situ grown La{sub 5/8}Ca{sub 3/8}MnO{sub 3} (LCMO) films. We demonstrate the ability of sliding FFT algorithm to differentiate two sub-lattices, resulting from a mixed-terminated surface. Principal Component Analysis and Independent Component Analysis of the Sliding FFT dataset reveal the distinct changes in crystallography, step edges, and boundaries between the multiple sub-lattices. The implications for the LCMO system are discussed. The method is universal for images with any periodicity, and is especially amenable to atomically resolved probe and electron-microscopy data for rapid identification of the sub-lattices present.
Lebensohn, Ricardo A; Lee, Sukbin; Rollett, Anthony D
2009-01-01
A viscoplastic approach using the Fast Fourier Transform (FFT) method for obtaining local mechanical response is utilized to study microstructure-property relationships in composite materials. Specifically, three-dimensional, two-phase digital materials containing isotropically coarsened particles surrounded by a matrix phase, generated through a Kinetic Monte Carlo Potts model for Ostwald ripening, are used as instantiations in order to calculate the stress and strain rate fields under uniaxial tension. The effects of the morphology of the matrix phase, the volume fraction and the contiguity of particles, and the polycrystallinity of matrix phase, on the stress and strain rate fields under uniaxial tension are examined. It is found that the first moments of the stress and strain rate fields have a different dependence on the particle volume fraction and the particle contiguity from their second moments. The average stresses and average strain rates of both phases and of the overall composite have rather simple relationships with the particle volume fraction whereas their standard deviations vary strongly, especially when the particle volume fraction is high, and the contiguity of particles has a noticeable effect on the mechanical response. It is also found that the shape of stress distribution in the BCC hard particle phase evolves as the volume fraction of particles in the composite varies, such that it agrees with the stress field in the BCC polycrystal as the volume of particles approaches unity. Finally, it is observed that the stress and strain rate fields in the microstructures with a polycrystalline matrix are less sensitive to changes in volume fraction and contiguity of particles.
Fast Fourier transformation analysis of kindling-induced afterdischarge in the rabbit hippocampus.
Tsuchiya, Komei; Kogure, Shinichi
2011-06-01
Kindling is a widely used animal model of intractable temporal lobe epilepsy. In the present study, we performed fast Fourier transformation (FFT) analysis of kindling-induced afterdischarge (AD) in the rabbit hippocampus. Ten adult rabbits were used. Kindling stimulation to the right hippocampus was delivered as a train of biphasic pulses (1 ms duration each) of 50 Hz for 1s, with suprathreshold intensity for AD. Motor responses were classified into five stages according to the conventional criteria. Of 10 animals, five developed stage 5 convulsions with a mean of 21 stimulations (kindled (K) group), while the remaining five animals did not (incomplete kindling (IK) group). We standardized each ratio of power spectral density of lower frequency band component (LFB: 0-9 Hz) and the higher frequency band (HFB: 12-30 Hz) in the initial stage as 1.0. The IK group exhibited small decrements (0.99 and 0.94 times) in LFB and HFB components at the final stage. In contrast, the K group exhibited a significantly (p<0.05) large decrement (0.49 times) in the LFB component and a very large increment (4.45 times) of HFB component at the final stage. Correlation analyses were performed between alteration of power spectral density ratio of the HFB component and AD duration, interictal discharge frequency, and behavioral stage during kindling progression. Fairly strong positive correlations were found in all cases in the K group. FFT analysis of kindling-induced AD demonstrated an important role of the HFB component: enhancement of the HFB component is associated with kindled stage, while decrement of it is associated with incomplete kindling stage. These findings suggest that FFT analysis of stimulus-induced and spontaneous seizure discharges is useful for examination of the progression of epileptic disorders. PMID:21498048
NASA Technical Reports Server (NTRS)
Tooms, S.; Attenborough, K.
1990-01-01
Using a Fast Fourier integration method and a global matrix method for solution of the boundary condition equations at all interfaces simultaneously, a useful tool for predicting acoustic propagation in a stratified fluid over a stratified porous-elastic solid was developed. The model for the solid is a modified Biot-Stoll model incorporating four parameters describing the pore structure corresponding to the Rayleigh-Attenborough rigid-porous structure model. The method is also compared to another Fast Fourier code (CERL-FFP) which models the ground as an impedance surface under a horizontally stratified air. Agreement with the CERL FFP is good. The effects on sound propagation of a combination of ground elasticity, complex ground structure, and atmospheric conditions are demonstrated by theoretical results over a snow layer, and experimental results over a model ground surface.
Fast Fourier and Wavelet Transforms for Wavefront Reconstruction in Adaptive Optics
Dowla, F U; Brase, J M; Olivier, S S
2000-07-28
Wavefront reconstruction techniques using the least-squares estimators are computationally quite expensive. We compare wavelet and Fourier transforms techniques in addressing the computation issues of wavefront reconstruction in adaptive optics. It is shown that because the Fourier approach is not simply a numerical approximation technique unlike the wavelet method, the Fourier approach might have advantages in terms of numerical accuracy. However, strictly from a numerical computations viewpoint, the wavelet approximation method might have advantage in terms of speed. To optimize the wavelet method, a statistical study might be necessary to use the best basis functions or ''approximation tree.''
On the Estimation of T-Wave Alternans Using the Spectral Fast Fourier Transform Method
Armoundas, Antonis A; Mela, Theofanie; Merchant, Faisal M
2012-01-01
BACKGROUND T-wave alternans (TWA), has been associated with increased vulnerability to ventricular tachyarrhythmias and sudden cardiac death (SCD). However, both random (white) noise and (patho)physiologic processes (i.e. premature ventricular contractions [PVCs], heart and respiration rates) may hamper TWA estimation and therefore, lessen its clinical utility for risk stratification. OBJECTIVE To investigate the effect of random noise and certain (patho)physiologic processes on the estimation of TWA using the Fast Fourier Transform (FFT) method and to develop methods to overcome these potential sources of error. METHODS We used a combination of human electrocardiogram data and computer simulations to assess the effects of a PVC, random and colored noise on the accuracy of TWA estimation. RESULTS We quantitatively demonstrate that replacing a “bad” beat with an odd/even median beat is a more accurate approach than replacing it with the overall average or the overall median beat. We also show that phase resetting may have a significant effect on alternans estimation and that estimation of alternans using frequencies greater than 0.4922 cycles/beat in a 128-point FFT provides the most accurate approach for estimating the alternans when phase resetting is likely to occur. Additionally, our data demonstrate that the number of indeterminate TWA tests due to high levels of noise can be reduced when the alternans voltage exceeds a new higher threshold. Also, the amplitude of random noise has a significant effect on alternans estimation and should be considered to adjust the alternans voltage threshold for noise levels greater than 1.8 μV. Finally, we quantitatively demonstrate that colored noise may lead to a false positive or a false negative result. We propose methods to estimate the effect of these (patho)physiologic processes on the alternans estimation in order to determine whether a TWA test is likely to be a true positive or a true negative. CONCLUSION This
A Fourier analysis for a fast simulation algorithm. [for switching converters
NASA Technical Reports Server (NTRS)
King, Roger J.
1988-01-01
This paper presents a derivation of compact expressions for the Fourier series analysis of the steady-state solution of a typical switching converter. The modeling procedure for the simulation and the steady-state solution is described, and some desirable traits for its matrix exponential subroutine are discussed. The Fourier analysis algorithm was tested on a phase-controlled parallel-loaded resonant converter, providing an experimental confirmation.
Using nonequispaced fast Fourier transformation to process optical coherence tomography signals
NASA Astrophysics Data System (ADS)
Hillmann, Dierck; Hüttmann, Gereon; Koch, Peter
2009-07-01
In OCT imaging the spectra that are used for Fourier transformation are in general not acquired linearly in k-space. Therefore one needs to apply an algorithm to re-sample the data and finally do the Fourier Transformation to gain depth information. We compare three algorithms (Non-Equispaced DFT, interpolated FFT and Non-Equispaced FFT) for this purpose in terms of speed and accuracy. The optimal algorithm depends on the OCT device (speed, SNR) and the object.
HIFI-C: a robust and fast method for determining NMR couplings from adaptive 3D to 2D projections.
Cornilescu, Gabriel; Bahrami, Arash; Tonelli, Marco; Markley, John L; Eghbalnia, Hamid R
2007-08-01
We describe a novel method for the robust, rapid, and reliable determination of J couplings in multi-dimensional NMR coupling data, including small couplings from larger proteins. The method, "High-resolution Iterative Frequency Identification of Couplings" (HIFI-C) is an extension of the adaptive and intelligent data collection approach introduced earlier in HIFI-NMR. HIFI-C collects one or more optimally tilted two-dimensional (2D) planes of a 3D experiment, identifies peaks, and determines couplings with high resolution and precision. The HIFI-C approach, demonstrated here for the 3D quantitative J method, offers vital features that advance the goal of rapid and robust collection of NMR coupling data. (1) Tilted plane residual dipolar couplings (RDC) data are collected adaptively in order to offer an intelligent trade off between data collection time and accuracy. (2) Data from independent planes can provide a statistical measure of reliability for each measured coupling. (3) Fast data collection enables measurements in cases where sample stability is a limiting factor (for example in the presence of an orienting medium required for residual dipolar coupling measurements). (4) For samples that are stable, or in experiments involving relatively stronger couplings, robust data collection enables more reliable determinations of couplings in shorter time, particularly for larger biomolecules. As a proof of principle, we have applied the HIFI-C approach to the 3D quantitative J experiment to determine N-C' RDC values for three proteins ranging from 56 to 159 residues (including a homodimer with 111 residues in each subunit). A number of factors influence the robustness and speed of data collection. These factors include the size of the protein, the experimental set up, and the coupling being measured, among others. To exhibit a lower bound on robustness and the potential for time saving, the measurement of dipolar couplings for the N-C' vector represents a realistic
NASA Astrophysics Data System (ADS)
Song, Qiyuan; Isobe, Keisuke; Hirosawa, Kenichi; Midorikawa, Katsumi; Kannari, Fumihiko
2015-03-01
Simultaneous spatial and temporal focusing (SSTF) multiphoton microscopy offers us widefield imaging with sectioning ability. As extending the idea to 2D SSTF, people can utilize a 2D spectral disperser. In this study, we use a 2D spectral disperser via a virtually-imaged phased-array (VIPA) and a diffraction grating to fulfill the back aperture of objective lens with a spectrum matrix. This offers us an axial resolution enhanced by a factor of ~1.7 compared with conventional SSTF microscopy. Furthermore, the small free spectral range (FSR) of VIPA will reduce the temporal self-imaging effect around out-of-focus region and thus will reduce the out-of-focus multiphoton excited fluorescence (MPEF) signal of 2D SSTF microscopy. We experimentally show that inside a sample with dense MPEF, the contrast of the sectioning image is increased in our 2D SSTF microscope compared with SSTF microscope. In our microscope, we use a 1 kHz chirped amplification laser, a piezo stage and a sCMOS camera integrated with 2D SSTF to realize high speed volume imaging at a speed of 50 volumes per second as well as improved sectioning ability. Volume imaging of Brownian motions of fluorescent beads as small as 1μm has been demonstrated. Not only the lateral motion but also the axial motion could be traced.
Xu, Daguang; Huang, Yong; Kang, Jin U
2014-06-16
We implemented the graphics processing unit (GPU) accelerated compressive sensing (CS) non-uniform in k-space spectral domain optical coherence tomography (SD OCT). Kaiser-Bessel (KB) function and Gaussian function are used independently as the convolution kernel in the gridding-based non-uniform fast Fourier transform (NUFFT) algorithm with different oversampling ratios and kernel widths. Our implementation is compared with the GPU-accelerated modified non-uniform discrete Fourier transform (MNUDFT) matrix-based CS SD OCT and the GPU-accelerated fast Fourier transform (FFT)-based CS SD OCT. It was found that our implementation has comparable performance to the GPU-accelerated MNUDFT-based CS SD OCT in terms of image quality while providing more than 5 times speed enhancement. When compared to the GPU-accelerated FFT based-CS SD OCT, it shows smaller background noise and less side lobes while eliminating the need for the cumbersome k-space grid filling and the k-linear calibration procedure. Finally, we demonstrated that by using a conventional desktop computer architecture having three GPUs, real-time B-mode imaging can be obtained in excess of 30 fps for the GPU-accelerated NUFFT based CS SD OCT with frame size 2048(axial) × 1,000(lateral). PMID:24977582
Xu, Daguang; Huang, Yong; Kang, Jin U.
2014-01-01
We implemented the graphics processing unit (GPU) accelerated compressive sensing (CS) non-uniform in k-space spectral domain optical coherence tomography (SD OCT). Kaiser-Bessel (KB) function and Gaussian function are used independently as the convolution kernel in the gridding-based non-uniform fast Fourier transform (NUFFT) algorithm with different oversampling ratios and kernel widths. Our implementation is compared with the GPU-accelerated modified non-uniform discrete Fourier transform (MNUDFT) matrix-based CS SD OCT and the GPU-accelerated fast Fourier transform (FFT)-based CS SD OCT. It was found that our implementation has comparable performance to the GPU-accelerated MNUDFT-based CS SD OCT in terms of image quality while providing more than 5 times speed enhancement. When compared to the GPU-accelerated FFT based-CS SD OCT, it shows smaller background noise and less side lobes while eliminating the need for the cumbersome k-space grid filling and the k-linear calibration procedure. Finally, we demonstrated that by using a conventional desktop computer architecture having three GPUs, real-time B-mode imaging can be obtained in excess of 30 fps for the GPU-accelerated NUFFT based CS SD OCT with frame size 2048(axial)×1000(lateral). PMID:24977582
New generation of Fourier optics instruments for fast multispectral BRDF characterization
NASA Astrophysics Data System (ADS)
Boher, Pierre; Leroux, Thierry; Collomb-Patton, Véronique; Bignon, Thibault
2015-03-01
A new generation of Fourier optics multispectral instruments that allow rapid full diffused or collimated beam spectral BRDF measurements is presented. Light detection is made simultaneously at all angular locations including the illumination direction. Backscattering effect in the fields of cosmetics and parasitic reflection of mobile displays are reported as examples.
Yang, Nuo; Hu, Shiqian; Ma, Dengke; Lu, Tingyu; Li, Baowen
2015-01-01
In this Paper, we investigate numerically and analytically the thermal conductivity of nanoscale graphene disks (NGDs), and discussed the possibility to realize functionally graded material (FGM) with only one material, NGDs. Different from previous studies on divergence/non-diffusive of thermal conductivity in nano-structures with different size, we found a novel non-homogeneous (graded) thermal conductivity along the radius direction in a single nano-disk structure. We found that, instead of a constant value, the NGD has a graded thermal conductivity along the radius direction. That is, Fourier's law of heat conduction is not valid in two dimensional graphene disk structures Moreover, we show the dependent of NGDs' thermal conductivity on radius and temperature. Our study might inspire experimentalists to develop NGD based versatile FGMs, improve understanding of the heat removal of hot spots on chips, and enhance thermoelectric energy conversion efficiency by two dimensional disk with a graded thermal conductivity. PMID:26443206
NASA Astrophysics Data System (ADS)
Maccone, Claudio
2007-04-01
The present article describes that the range of any radiotelescope (and radar in general) may be increased by virtue of software, if one replaces the fast Fourier transform by the Karhunen Loève transform. The range increases with the inverse of the fourth root of the signal-to-noise ratio when this ratio decreases. Thus, the range on any radiotelescope (and radar) may be increased without changing the hardware at all, but by changing the software only. This improvement in the range of the radiotelescope is currently implemented at the 32-m antenna located at Medicina, near Bologna, in Italy, for both SETI and general radioastronomy.
Sasaki, Keiji; Inayoshi, Naruhiko; Tashiro, Kohji
2008-12-01
To observe the time-dependent two-dimensional (2D) images of spatial distribution of chemically/physically modified lubricant molecules on the metal surface during friction motion, a new in situ technique has been developed by combining the 2D fast-imaging Fourier-transform infrared-attenuated total reflection spectrometer with the temperature-controlled friction equipment containing lubricant agent. Using this new instrument, the time-dependent changes in lubricant molecules, for example, cis-trans isomerization, stress-induced molecular deformation, etc., can be detected successfully. The characteristic features of this instrument have been demonstrated in a detailed and concrete manner by demonstrating the experimental data measured for oleic acid and tricresyl phosphate. PMID:19123564
Short time Fourier analysis of the electromyogram - Fast movements and constant contraction
NASA Technical Reports Server (NTRS)
Hannaford, Blake; Lehman, Steven
1986-01-01
Short-time Fourier analysis was applied to surface electromyograms (EMG) recorded during rapid movements, and during isometric contractions at constant forces. A portion of the data to be transformed by multiplying the signal by a Hamming window was selected, and then the discrete Fourier transform was computed. Shifting the window along the data record, a new spectrum was computed each 10 ms. The transformed data were displayed in spectograms or 'voiceprints'. This short-time technique made it possible to see time-dependencies in the EMG that are normally averaged in the Fourier analysis of these signals. Spectra of EMGs during isometric contractions at constant force vary in the short (10-20 ms) term. Short-time spectra from EMGs recorded during rapid movements were much less variable. The windowing technique picked out the typical 'three-burst pattern' in EMG's from both wrist and head movements. Spectra during the bursts were more consistent than those during isometric contractions. Furthermore, there was a consistent shift in spectral statistics in the course of the three bursts. Both the center frequency and the variance of the spectral energy distribution grew from the first burst to the second burst in the same muscle. The analogy between EMGs and speech signals is extended to argue for future applicability of short-time spectral analysis of EMG.
NASA Astrophysics Data System (ADS)
Windhari, Ayuty; Handayani, Gunawan
2015-04-01
The 3D inversion gravity anomaly to estimate topographical density using a matlab source code from gridded data provided by Parker Oldenburg algorithm based on fast Fourier transform was computed. We extend and improved the source code of 3DINVERT.M invented by Gomez Ortiz and Agarwal (2005) using the relationship between Fourier transform of the gravity anomaly and the sum of the Fourier transform from the topography density. We gave density contrast between the two media to apply the inversion. FFT routine was implemented to construct amplitude spectrum to the given mean depth. The results were presented as new graphics of inverted topography density, the gravity anomaly due to the inverted topography and the difference between the input gravity data and the computed ones. It terminates when the RMS error is lower than pre-assigned value used as convergence criterion or until maximum of iterations is reached. As an example, we used the matlab program on gravity data of Banten region, Indonesia.
NASA Astrophysics Data System (ADS)
Chang, Chenliang; Xia, Jun; Lei, Wei
2015-03-01
We proposed a new method to calculate the color computer generated hologram of three-dimensional object in holographic display. The three-dimensional object is composed of several tilted planes which are tilted from the hologram. The diffraction from each tilted plane to the hologram plane is calculated based on the coordinate rotation in Fourier spectrum domains. We used the nonuniform fast Fourier transformation (NUFFT) to calculate the nonuniform sampled Fourier spectrum on the tilted plane after coordinate rotation. By using the NUFFT, the diffraction calculation from tilted plane to the hologram plane with variable sampling rates can be achieved, which overcomes the sampling restriction of FFT in the conventional angular spectrum based method. The holograms of red, green and blue component of the polygon-based object are calculated separately by using our NUFFT based method. Then the color hologram is synthesized by placing the red, green and blue component hologram in sequence. The chromatic aberration caused by the wavelength difference can be solved effectively by restricting the sampling rate of the object in the calculation of each wavelength component. The computer simulation shows the feasibility of our method in calculating the color hologram of polygon-based object. The 3D object can be displayed in color with adjustable size and no chromatic aberration in holographic display system, which can be considered as an important application in the colorful holographic three-dimensional display.
NASA Astrophysics Data System (ADS)
Thoen, D. J.; Bongers, W. A.; Westerhof, E.; Oosterbeek, J. W.; de Baar, M. R.; van den Berg, M. A.; van Beveren, V.; Bürger, A.; Goede, A. P. H.; Graswinckel, M. F.; Hennen, B. A.; Schüller, F. C.
2009-10-01
A fast Fourier transform (FFT) based wide range millimeter wave diagnostics for spectral characterization of scattered millimeter waves in plasmas has been successfully brought into operation. The scattered millimeter waves are heterodyne downconverted and directly digitized using a fast analog-digital converter and a compact peripheral component interconnect computer. Frequency spectra are obtained by FFT in the time domain of the intermediate frequency signal. The scattered millimeter waves are generated during high power electron cyclotron resonance heating experiments on the TEXTOR tokamak and demonstrate the performance of the diagnostics and, in particular, the usability of direct digitizing and Fourier transformation of millimeter wave signals. The diagnostics is able to acquire 4 GHz wide spectra of signals in the range of 136-140 GHz. The rate of spectra is tunable and has been tested between 200 000 spectra/s with a frequency resolution of 100 MHz and 120 spectra/s with a frequency resolution of 25 kHz. The respective dynamic ranges are 52 and 88 dB. Major benefits of the new diagnostics are a tunable time and frequency resolution due to postdetection, near-real time processing of the acquired data. This diagnostics has a wider application in astrophysics, earth observation, plasma physics, and molecular spectroscopy for the detection and analysis of millimeter wave radiation, providing high-resolution spectra at high temporal resolution and large dynamic range.
Thoen, D J; Bongers, W A; Westerhof, E; Oosterbeek, J W; de Baar, M R; van den Berg, M A; van Beveren, V; Bürger, A; Goede, A P H; Graswinckel, M F; Hennen, B A; Schüller, F C
2009-10-01
A fast Fourier transform (FFT) based wide range millimeter wave diagnostics for spectral characterization of scattered millimeter waves in plasmas has been successfully brought into operation. The scattered millimeter waves are heterodyne downconverted and directly digitized using a fast analog-digital converter and a compact peripheral component interconnect computer. Frequency spectra are obtained by FFT in the time domain of the intermediate frequency signal. The scattered millimeter waves are generated during high power electron cyclotron resonance heating experiments on the TEXTOR tokamak and demonstrate the performance of the diagnostics and, in particular, the usability of direct digitizing and Fourier transformation of millimeter wave signals. The diagnostics is able to acquire 4 GHz wide spectra of signals in the range of 136-140 GHz. The rate of spectra is tunable and has been tested between 200,000 spectra/s with a frequency resolution of 100 MHz and 120 spectra/s with a frequency resolution of 25 kHz. The respective dynamic ranges are 52 and 88 dB. Major benefits of the new diagnostics are a tunable time and frequency resolution due to postdetection, near-real time processing of the acquired data. This diagnostics has a wider application in astrophysics, earth observation, plasma physics, and molecular spectroscopy for the detection and analysis of millimeter wave radiation, providing high-resolution spectra at high temporal resolution and large dynamic range. PMID:19895061
Addendum to 'A new hybrid algorithm for computing a fast discrete Fourier transform'
NASA Technical Reports Server (NTRS)
Reed, I. S.; Truong, T. K.; Benjauthrit, B.
1981-01-01
The reported investigation represents a continuation of a study conducted by Reed and Truong (1979), who proposed a hybrid algorithm for computing the discrete Fourier transform (DFT). The proposed technique employs a Winograd-type algorithm in conjunction with the Mersenne prime-number theoretic transform to perform a DFT. The implementation of the technique involves a considerable number of additions. The new investigation shows an approach which can reduce the number of additions significantly. It is proposed to use Winograd's algorithm for computing the Mersenne prime-number theoretic transform in the transform portion of the hybrid algorithm.
Fast multi-core based multimodal registration of 2D cross-sections and 3D datasets
2010-01-01
Background Solving bioinformatics tasks often requires extensive computational power. Recent trends in processor architecture combine multiple cores into a single chip to improve overall performance. The Cell Broadband Engine (CBE), a heterogeneous multi-core processor, provides power-efficient and cost-effective high-performance computing. One application area is image analysis and visualisation, in particular registration of 2D cross-sections into 3D image datasets. Such techniques can be used to put different image modalities into spatial correspondence, for example, 2D images of histological cuts into morphological 3D frameworks. Results We evaluate the CBE-driven PlayStation 3 as a high performance, cost-effective computing platform by adapting a multimodal alignment procedure to several characteristic hardware properties. The optimisations are based on partitioning, vectorisation, branch reducing and loop unrolling techniques with special attention to 32-bit multiplies and limited local storage on the computing units. We show how a typical image analysis and visualisation problem, the multimodal registration of 2D cross-sections and 3D datasets, benefits from the multi-core based implementation of the alignment algorithm. We discuss several CBE-based optimisation methods and compare our results to standard solutions. More information and the source code are available from http://cbe.ipk-gatersleben.de. Conclusions The results demonstrate that the CBE processor in a PlayStation 3 accelerates computational intensive multimodal registration, which is of great importance in biological/medical image processing. The PlayStation 3 as a low cost CBE-based platform offers an efficient option to conventional hardware to solve computational problems in image processing and bioinformatics. PMID:20064262
Lin, Lin; Yang, Chao; Lu, Jiangfeng; Ying, Lexing; E, Weinan
2009-09-25
We present an efficient parallel algorithm and its implementation for computing the diagonal of $H^-1$ where $H$ is a 2D Kohn-Sham Hamiltonian discretized on a rectangular domain using a standard second order finite difference scheme. This type of calculation can be used to obtain an accurate approximation to the diagonal of a Fermi-Dirac function of $H$ through a recently developed pole-expansion technique \\cite{LinLuYingE2009}. The diagonal elements are needed in electronic structure calculations for quantum mechanical systems \\citeHohenbergKohn1964, KohnSham 1965,DreizlerGross1990. We show how elimination tree is used to organize the parallel computation and how synchronization overhead is reduced by passing data level by level along this tree using the technique of local buffers and relative indices. We analyze the performance of our implementation by examining its load balance and communication overhead. We show that our implementation exhibits an excellent weak scaling on a large-scale high performance distributed parallel machine. When compared with standard approach for evaluating the diagonal a Fermi-Dirac function of a Kohn-Sham Hamiltonian associated a 2D electron quantum dot, the new pole-expansion technique that uses our algorithm to compute the diagonal of $(H-z_i I)^-1$ for a small number of poles $z_i$ is much faster, especially when the quantum dot contains many electrons.
Mori, Shinichiro; Kumagai, Motoki; Miki, Kentaro; Fukuhara, Riki; Haneishi, Hideaki
2015-09-01
To improve treatment workflow, we developed a graphic processing unit (GPU)-based patient positional verification software application and integrated it into carbon-ion scanning beam treatment. Here, we evaluated the basic performance of the software. The algorithm provides 2D/3D registration matching using CT and orthogonal X-ray flat panel detector (FPD) images. The participants were 53 patients with tumors of the head and neck, prostate or lung receiving carbon-ion beam treatment. 2D/3D-ITchi-Gime (ITG) calculation accuracy was evaluated in terms of computation time and registration accuracy. Registration calculation was determined using the similarity measurement metrics gradient difference (GD), normalized mutual information (NMI), zero-mean normalized cross-correlation (ZNCC), and their combination. Registration accuracy was dependent on the particular metric used. Representative examples were determined to have target registration error (TRE) = 0.45 ± 0.23 mm and angular error (AE) = 0.35 ± 0.18° with ZNCC + GD for a head and neck tumor; TRE = 0.12 ± 0.07 mm and AE = 0.16 ± 0.07° with ZNCC for a pelvic tumor; and TRE = 1.19 ± 0.78 mm and AE = 0.83 ± 0.61° with ZNCC for lung tumor. Calculation time was less than 7.26 s.The new registration software has been successfully installed and implemented in our treatment process. We expect that it will improve both treatment workflow and treatment accuracy. PMID:26081313
Mori, Shinichiro; Kumagai, Motoki; Miki, Kentaro; Fukuhara, Riki; Haneishi, Hideaki
2015-01-01
To improve treatment workflow, we developed a graphic processing unit (GPU)-based patient positional verification software application and integrated it into carbon-ion scanning beam treatment. Here, we evaluated the basic performance of the software. The algorithm provides 2D/3D registration matching using CT and orthogonal X-ray flat panel detector (FPD) images. The participants were 53 patients with tumors of the head and neck, prostate or lung receiving carbon-ion beam treatment. 2D/3D-ITchi-Gime (ITG) calculation accuracy was evaluated in terms of computation time and registration accuracy. Registration calculation was determined using the similarity measurement metrics gradient difference (GD), normalized mutual information (NMI), zero-mean normalized cross-correlation (ZNCC), and their combination. Registration accuracy was dependent on the particular metric used. Representative examples were determined to have target registration error (TRE) = 0.45 ± 0.23 mm and angular error (AE) = 0.35 ± 0.18° with ZNCC + GD for a head and neck tumor; TRE = 0.12 ± 0.07 mm and AE = 0.16 ± 0.07° with ZNCC for a pelvic tumor; and TRE = 1.19 ± 0.78 mm and AE = 0.83 ± 0.61° with ZNCC for lung tumor. Calculation time was less than 7.26 s.The new registration software has been successfully installed and implemented in our treatment process. We expect that it will improve both treatment workflow and treatment accuracy. PMID:26081313
NASA Astrophysics Data System (ADS)
Chae, Kyu-Hyun
2002-04-01
Fourier series solutions to the deflection and magnification by a family of three-dimensional cusped two-power-law ellipsoidal mass distributions are presented. The cusped two-power-law ellipsoidal mass distributions are characterized by inner and outer power-law radial indices and a break (or transition) radius. The model family includes mass models mimicking Jaffe, Hernquist, and η models and dark matter halo profiles from numerical simulations. The Fourier series solutions for the cusped two-power-law mass distributions are relatively simple and allow a very fast calculation, even for a chosen small fractional calculational error (e.g., 10-5). These results will be particularly useful for studying lensed systems that provide a number of accurate lensing constraints and for systematic analyses of large numbers of lenses. Subroutines employing these results for the two-power-law model and the results by Chae, Khersonsky, & Turnshek for the generalized single-power-law mass model are made publicly available.
Mackie, David M.; Jahnke, Justin P.; Benyamin, Marcus S.; Sumner, James J.
2016-01-01
The standard methodologies for quantitative analysis (QA) of mixtures using Fourier transform infrared (FTIR) instruments have evolved until they are now more complicated than necessary for many users’ purposes. We present a simpler methodology, suitable for widespread adoption of FTIR QA as a standard laboratory technique across disciplines by occasional users.•Algorithm is straightforward and intuitive, yet it is also fast, accurate, and robust.•Relies on component spectra, minimization of errors, and local adaptive mesh refinement.•Tested successfully on real mixtures of up to nine components. We show that our methodology is robust to challenging experimental conditions such as similar substances, component percentages differing by three orders of magnitude, and imperfect (noisy) spectra. As examples, we analyze biological, chemical, and physical aspects of bio-hybrid fuel cells. PMID:26977411
Talhaoui, Hicham; Menacer, Arezki; Kessal, Abdelhalim; Kechida, Ridha
2014-09-01
This paper presents new techniques to evaluate faults in case of broken rotor bars of induction motors. Procedures are applied with closed-loop control. Electrical and mechanical variables are treated using fast Fourier transform (FFT), and discrete wavelet transform (DWT) at start-up and steady state. The wavelet transform has proven to be an excellent mathematical tool for the detection of the faults particularly broken rotor bars type. As a performance, DWT can provide a local representation of the non-stationary current signals for the healthy machine and with fault. For sensorless control, a Luenberger observer is applied; the estimation rotor speed is analyzed; the effect of the faults in the speed pulsation is compensated; a quadratic current appears and used for fault detection. PMID:25004798
SnO2-MOF-Fabry-Perot humidity optical sensor system based on fast Fourier transform technique
NASA Astrophysics Data System (ADS)
Lopez-Aldaba, A.; Lopez-Torres, D.; Ascorbe, J.; Rota-Rodrigo, S.; Elosua, C.; Lopez-Amo, M.; Arregui, F. J.; Corres, J. M.; Auguste, J.-L.; Jamier, R.; Roy, P.
2016-05-01
In this paper, a new sensor system for relative humidity measurements based on a SnO2 sputtering deposition on a microstructured optical fiber (MOF) low-finesse Fabry-Perot (FP) sensing head is presented and characterized. The interrogation of the sensing head is carried out by monitoring the Fast Fourier Transform phase variations of the FP interference frequency. This method is low-sensitive to signal amplitude variations and also avoids the necessity of tracking the evolution of peaks and valleys in the spectrum. The sensor is operated within a wide humidity range (20%-90% relative humidity) with a maximum sensitivity achieved of 0.14rad/%. The measurement method uses a commercial optical interrogator as the only active element, this compact solution allows real time analysis of the data.
Bettadapura, Radhakrishna; Rasheed, Muhibur; Vollrath, Antje; Bajaj, Chandrajit
2015-10-01
There continue to be increasing occurrences of both atomistic structure models in the PDB (possibly reconstructed from X-ray diffraction or NMR data), and 3D reconstructed cryo-electron microscopy (3D EM) maps (albeit at coarser resolution) of the same or homologous molecule or molecular assembly, deposited in the EMDB. To obtain the best possible structural model of the molecule at the best achievable resolution, and without any missing gaps, one typically aligns (match and fits) the atomistic structure model with the 3D EM map. We discuss a new algorithm and generalized framework, named PF(2) fit (Polar Fast Fourier Fitting) for the best possible structural alignment of atomistic structures with 3D EM. While PF(2) fit enables only a rigid, six dimensional (6D) alignment method, it augments prior work on 6D X-ray structure and 3D EM alignment in multiple ways: Scoring. PF(2) fit includes a new scoring scheme that, in addition to rewarding overlaps between the volumes occupied by the atomistic structure and 3D EM map, rewards overlaps between the volumes complementary to them. We quantitatively demonstrate how this new complementary scoring scheme improves upon existing approaches. PF(2) fit also includes two scoring functions, the non-uniform exterior penalty and the skeleton-secondary structure score, and implements the scattering potential score as an alternative to traditional Gaussian blurring. Search. PF(2) fit utilizes a fast polar Fourier search scheme, whose main advantage is the ability to search over uniformly and adaptively sampled subsets of the space of rigid-body motions. PF(2) fit also implements a new reranking search and scoring methodology that considerably improves alignment metrics in results obtained from the initial search. PMID:26469938
Bettadapura, Radhakrishna; Rasheed, Muhibur; Vollrath, Antje; Bajaj, Chandrajit
2015-01-01
There continue to be increasing occurrences of both atomistic structure models in the PDB (possibly reconstructed from X-ray diffraction or NMR data), and 3D reconstructed cryo-electron microscopy (3D EM) maps (albeit at coarser resolution) of the same or homologous molecule or molecular assembly, deposited in the EMDB. To obtain the best possible structural model of the molecule at the best achievable resolution, and without any missing gaps, one typically aligns (match and fits) the atomistic structure model with the 3D EM map. We discuss a new algorithm and generalized framework, named PF2 fit (Polar Fast Fourier Fitting) for the best possible structural alignment of atomistic structures with 3D EM. While PF2 fit enables only a rigid, six dimensional (6D) alignment method, it augments prior work on 6D X-ray structure and 3D EM alignment in multiple ways: Scoring. PF2 fit includes a new scoring scheme that, in addition to rewarding overlaps between the volumes occupied by the atomistic structure and 3D EM map, rewards overlaps between the volumes complementary to them. We quantitatively demonstrate how this new complementary scoring scheme improves upon existing approaches. PF2 fit also includes two scoring functions, the non-uniform exterior penalty and the skeleton-secondary structure score, and implements the scattering potential score as an alternative to traditional Gaussian blurring. Search. PF2 fit utilizes a fast polar Fourier search scheme, whose main advantage is the ability to search over uniformly and adaptively sampled subsets of the space of rigid-body motions. PF2 fit also implements a new reranking search and scoring methodology that considerably improves alignment metrics in results obtained from the initial search. PMID:26469938
NASA Astrophysics Data System (ADS)
Otake, Y.; Schafer, S.; Stayman, J. W.; Zbijewski, W.; Kleinszig, G.; Graumann, R.; Khanna, A. J.; Siewerdsen, J. H.
2012-02-01
Localization of target vertebrae is an essential step in minimally invasive spine surgery, with conventional methods relying on "level counting" - i.e., manual counting of vertebrae under fluoroscopy starting from readily identifiable anatomy (e.g., the sacrum). The approach requires an undesirable level of radiation, time, and is prone to counting errors due to the similar appearance of vertebrae in projection images; wrong-level surgery occurs in 1 of every ~3000 cases. This paper proposes a method to automatically localize target vertebrae in x-ray projections using 3D-2D registration between preoperative CT (in which vertebrae are preoperatively labeled) and intraoperative fluoroscopy. The registration uses an intensity-based approach with a gradient-based similarity metric and the CMA-ES algorithm for optimization. Digitally reconstructed radiographs (DRRs) and a robust similarity metric are computed on GPU to accelerate the process. Evaluation in clinical CT data included 5,000 PA and LAT projections randomly perturbed to simulate human variability in setup of mobile intraoperative C-arm. The method demonstrated 100% success for PA view (projection error: 0.42mm) and 99.8% success for LAT view (projection error: 0.37mm). Initial implementation on GPU provided automatic target localization within about 3 sec, with further improvement underway via multi-GPU. The ability to automatically label vertebrae in fluoroscopy promises to streamline surgical workflow, improve patient safety, and reduce wrong-site surgeries, especially in large patients for whom manual methods are time consuming and error prone.
NASA Astrophysics Data System (ADS)
Rawlinson, N.; Sambridge, M.
2003-12-01
The accurate prediction of seismic traveltimes in layered media is required in many areas of seismology. In addition to simple refractions and reflections, complex phases comprising numerous transmission and reflection branches may exist; for instance, the so-called ``multiples" frequently identified in marine reflection seismology. We present a grid-based method for the accurate determination of multi-phase traveltimes in layered media of significant complexity. A finite difference eikonal solver known as the Fast Marching Method (FMM) is used to track wavefronts within a layer. FMM is a fast and unconditionally stable upwind scheme that is well suited to complex models, and can be used sequentially to track the multiple refraction and/or reflection branches of virtually any required phase. Although FMM was initially introduced as a first-order scheme, higher order operators can be used. A mixed-order scheme that preferentially uses second-order operators, but reverts to first-order operators when the required upwind traveltimes are unavailable, is one possibility. Despite improved accuracy, this scheme still suffers from first-order convergence due to high wavefront curvature and first-order accuracy in the vicinity of the source. To overcome this problem, we implement local grid refinement about the source. In order to retain stability, the edge of the refined grid conforms to the shape of the wavefront, so that information only flows out of the refined grid, and never back into it. Application of our new scheme to complex velocity media shows that grid refinement typically improves accuracy by an order of magnitude, with only a small increase in computation time ( ˜5%). Significantly, first-order convergence is replaced by near second-order convergence, even in media with velocity contrasts as large as 8:1. In one example, with a velocity grid defined by 257,121 nodes, reflection traveltimes from a strongly undulating interface were calculated with an error of
Abrishami, V; Bilbao-Castro, J R; Vargas, J; Marabini, R; Carazo, J M; Sorzano, C O S
2015-10-01
We describe a fast and accurate method for the reconstruction of macromolecular complexes from a set of projections. Direct Fourier inversion (in which the Fourier Slice Theorem plays a central role) is a solution for dealing with this inverse problem. Unfortunately, the set of projections provides a non-equidistantly sampled version of the macromolecule Fourier transform in the single particle field (and, therefore, a direct Fourier inversion) may not be an optimal solution. In this paper, we introduce a gridding-based direct Fourier method for the three-dimensional reconstruction approach that uses a weighting technique to compute a uniform sampled Fourier transform. Moreover, the contrast transfer function of the microscope, which is a limiting factor in pursuing a high resolution reconstruction, is corrected by the algorithm. Parallelization of this algorithm, both on threads and on multiple CPU's, makes the process of three-dimensional reconstruction even faster. The experimental results show that our proposed gridding-based direct Fourier reconstruction is slightly more accurate than similar existing methods and presents a lower computational complexity both in terms of time and memory, thereby allowing its use on larger volumes. The algorithm is fully implemented in the open-source Xmipp package and is downloadable from http://xmipp.cnb.csic.es. PMID:26094203
NASA Astrophysics Data System (ADS)
Weese, Juergen; Goecke, Roland; Penney, Graeme P.; Desmedt, Paul; Buzug, Thorsten M.; Schumann, Heidrun
1999-05-01
2D/3D registration makes it possible to use pre-operative CT scans for navigation purposes during X-ray fluoroscopy guided interventions. We present a fast voxel-based method for this registration task, which uses a recently introduced similarity measure (pattern intensity). This measure is especially suitable for 2D/3D registration, because it is robust with respect to structures such as a stent visible in the X-ray fluoroscopy image but not in the CT scan. The method uses only a part of the CT scan for the generation of digitally reconstructed radiographs (DRRs) to accelerate their computation. Nevertheless, computation time is crucial for intra-operative application and a further speed-up is required, because numerous DRRs must be computed. For that reason, the suitability of different volume rendering methods for 2D/3D registration has been investigated. A method based on the shear-warp factorization of the viewing transformation turned out to be especially suitable and builds the basis of the registration algorithm. The algorithm has been applied to images of a spine phantom and to clinical images. For comparison, registration results have been calculated using ray-casting. The shear-warp factorization based rendering method accelerates registration by a factor of up to seven compared to ray-casting without degrading registration accuracy. Using a vertebra as feature for registration, computation time is in the range of 3-4s (Sun UltraSparc, 300 MHz) which is acceptable for intra-operative application.
Kasinathan, N.; Rajakumar, A.; Vaidyanathan, G.; Chetal, S.C.
1995-09-01
Post shutdown decay heat removal is an important safety requirement in any nuclear system. In order to improve the reliability of this function, Liquid metal (sodium) cooled fast breeder reactors (LMFBR) are equipped with redundant hot pool dipped immersion coolers connected to natural draught air cooled heat exchangers through intermediate sodium circuits. During decay heat removal, flow through the core, immersion cooler primary side and in the intermediate sodium circuits are also through natural convection. In order to establish the viability and validate computer codes used in making predictions, a 1:20 scale experimental model called RAMONA with water as coolant has been built and experimental simulation of decay heat removal situation has been performed at KfK Karlsruhe. Results of two such experiments have been compiled and published as benchmarks. This paper brings out the results of the numerical simulation of one of the benchmark case through a 1D/2D coupled code system, DHDYN-1D/THYC-2D and the salient features of the comparisons. Brief description of the formulations of the codes are also included.
NASA Astrophysics Data System (ADS)
Seidel, A.; Wagner, S.; Dreizler, A.; Ebert, V.
2015-05-01
We have developed a fast, spatially scanning direct tunable diode laser absorption spectrometer (dTDLAS) that combines four polygon-mirror based scanning units with low-cost retro-reflective foils. With this instrument, tomographic measurements of absolute 2-D water vapor concentration profiles are possible without any calibration using a reference gas. A spatial area of 0.8 m x 0.8 m was covered, which allows for application in soil physics, where greenhouse gas emission from certain soil structures shall be monitored. The whole concentration field was measured with up to 2.5 Hz. In this paper, we present the setup and spectroscopic performance of the instrument regarding the influence of the polygon rotation speed and mode on the absorption signal. Homogeneous H2O distributions were measured and compared to a single channel, bi-static reference TDLAS spectrometer for validation of the instrument. Good accuracy and precision with errors of less than 6% of the absolute concentration and length and bandwidth normalized detection limits of up to 1.1 ppmv . m (Hz)-0.5 were achieved. The spectrometer is a robust and easy to set up instrument for tomographic reconstructions of 2-D-concentration fields that can be considered as a good basis for future field measurements in environmental research.
NASA Astrophysics Data System (ADS)
Seidel, A.; Wagner, S.; Dreizler, A.; Ebert, V.
2014-12-01
We have developed a fast, spatially direct scanning tunable diode laser absorption spectrometer (dTDLAS) that combines four polygon-mirror based scanning units with low-cost retro-reflective foils. With this instrument, tomographic measurements of absolute 2-D water vapour concentration profiles are possible without any calibration using a reference gas. A spatial area of 0.8 m × 0.8 m was covered, which allows for application in soil physics, where greenhouse gas emission from certain soil structures shall be monitored. The whole concentration field was measured with up to 2.5 Hz. In this paper, we present the setup and spectroscopic performance of the instrument regarding the influence of the polygon rotation speed and mode on the absorption signal. Homogeneous H2O distributions were measured and compared to a single channel, bi-static reference TDLAS spectrometer for validation of the instrument. Good accuracy and precision with errors of less than 6% of the absolute concentration and length and bandwidth normalized detection limits of up to 1.1 ppmv · m · √Hz-1 were achieved. The spectrometer is a robust and easy to set up instrument for tomographic reconstructions of 2-D-concentration fields that can be considered a good basis for future field measurements in environmental research.
Konrad-Martin, D; Neely, S T; Keefe, D H; Dorn, P A; Gorga, M P
2001-06-01
Primary and secondary sources combine to produce the 2f1-f2 distortion product otoacoustic emission (DPOAE) measured in the ear canals of humans. DPOAEs were obtained in nine normal-hearing subjects using a fixed-f2 paradigm in which f1 was varied. The f2 was 2 or 4 kHz, and absolute and relative primary levels were varied. Data were obtained with and without a third tone (f3) placed 15.6 Hz below 2f1-f2. The level of f3 was varied in order to suppress the stimulus frequency otoacoustic emission (SFOAE) coming from the 2f1-f2 place. These data were converted from the complex frequency domain into an equivalent time representation using an inverse fast Fourier transform (IFFT). IFFTs of unsuppressed DPOAE data were characterized by two or more peaks. Relative amplitudes of these peaks depended on overall primary level and on primary-level differences. The suppressor eliminated later peaks, but early peaks remained relatively unaltered. Results are interpreted to mean that the DPOAE measured in humans includes components from the f2 place (intermodulation distortion) and DP place (in the form of a SFOAE). These findings build on previous work by providing evidence that multiple peaks in the IFFT are due to a secondary source at the DP place. PMID:11425129
NASA Astrophysics Data System (ADS)
Hillerkuss, D.; Schmogrow, R.; Schellinger, T.; Jordan, M.; Winter, M.; Huber, G.; Vallaitis, T.; Bonk, R.; Kleinow, P.; Frey, F.; Roeger, M.; Koenig, S.; Ludwig, A.; Marculescu, A.; Li, J.; Hoh, M.; Dreschmann, M.; Meyer, J.; Ben Ezra, S.; Narkiss, N.; Nebendahl, B.; Parmigiani, F.; Petropoulos, P.; Resan, B.; Oehler, A.; Weingarten, K.; Ellermeyer, T.; Lutz, J.; Moeller, M.; Huebner, M.; Becker, J.; Koos, C.; Freude, W.; Leuthold, J.
2011-06-01
Optical transmission systems with terabit per second (Tbit s-1) single-channel line rates no longer seem to be too far-fetched. New services such as cloud computing, three-dimensional high-definition television and virtual-reality applications require unprecedented optical channel bandwidths. These high-capacity optical channels, however, are fed from lower-bitrate signals. The question then is whether the lower-bitrate tributary information can viably, energy-efficiently and effortlessly be encoded to and extracted from terabit per second data streams. We demonstrate an optical fast Fourier transform scheme that provides the necessary computing power to encode lower-bitrate tributaries into 10.8 and 26.0 Tbit s-1 line-rate orthogonal frequency-division multiplexing (OFDM) data streams and to decode them from fibre-transmitted OFDM data streams. Experiments show the feasibility and ease of handling terabit per second data with low energy consumption. To the best of our knowledge, this is the largest line rate ever encoded onto a single light source.
Bhanot, Gyan V.; Chen, Dong; Gara, Alan G.; Giampapa, Mark E.; Heidelberger, Philip; Steinmacher-Burow, Burkhard D.; Vranas, Pavlos M.
2008-01-01
The present in invention is directed to a method, system and program storage device for efficiently implementing a multidimensional Fast Fourier Transform (FFT) of a multidimensional array comprising a plurality of elements initially distributed in a multi-node computer system comprising a plurality of nodes in communication over a network, comprising: distributing the plurality of elements of the array in a first dimension across the plurality of nodes of the computer system over the network to facilitate a first one-dimensional FFT; performing the first one-dimensional FFT on the elements of the array distributed at each node in the first dimension; re-distributing the one-dimensional FFT-transformed elements at each node in a second dimension via "all-to-all" distribution in random order across other nodes of the computer system over the network; and performing a second one-dimensional FFT on elements of the array re-distributed at each node in the second dimension, wherein the random order facilitates efficient utilization of the network thereby efficiently implementing the multidimensional FFT. The "all-to-all" re-distribution of array elements is further efficiently implemented in applications other than the multidimensional FFT on the distributed-memory parallel supercomputer.
Bhanot, Gyan V.; Chen, Dong; Gara, Alan G.; Giampapa, Mark E.; Heidelberger, Philip; Steinmacher-Burow, Burkhard D.; Vranas, Pavlos M.
2012-01-10
The present in invention is directed to a method, system and program storage device for efficiently implementing a multidimensional Fast Fourier Transform (FFT) of a multidimensional array comprising a plurality of elements initially distributed in a multi-node computer system comprising a plurality of nodes in communication over a network, comprising: distributing the plurality of elements of the array in a first dimension across the plurality of nodes of the computer system over the network to facilitate a first one-dimensional FFT; performing the first one-dimensional FFT on the elements of the array distributed at each node in the first dimension; re-distributing the one-dimensional FFT-transformed elements at each node in a second dimension via "all-to-all" distribution in random order across other nodes of the computer system over the network; and performing a second one-dimensional FFT on elements of the array re-distributed at each node in the second dimension, wherein the random order facilitates efficient utilization of the network thereby efficiently implementing the multidimensional FFT. The "all-to-all" re-distribution of array elements is further efficiently implemented in applications other than the multidimensional FFT on the distributed-memory parallel supercomputer.
NASA Astrophysics Data System (ADS)
Polkowski, Marcin
2016-04-01
Seismic wave travel time calculation is the most common numerical operation in seismology. The most efficient is travel time calculation in 1D velocity model - for given source, receiver depths and angular distance time is calculated within fraction of a second. Unfortunately, in most cases 1D is not enough to encounter differentiating local and regional structures. Whenever possible travel time through 3D velocity model has to be calculated. It can be achieved using ray calculation or time propagation in space. While single ray path calculation is quick it is complicated to find the ray path that connects source with the receiver. Time propagation in space using Fast Marching Method seems more efficient in most cases, especially when there are multiple receivers. In this presentation a Python module pySeismicFMM is presented - simple and very efficient tool for calculating travel time from sources to receivers. Calculation requires regular 2D or 3D velocity grid either in Cartesian or geographic coordinates. On desktop class computer calculation speed is 200k grid cells per second. Calculation has to be performed once for every source location and provides travel time to all receivers. pySeismicFMM is free and open source. Development of this tool is a part of authors PhD thesis. National Science Centre Poland provided financial support for this work via NCN grant DEC-2011/02/A/ST10/00284.
NASA Astrophysics Data System (ADS)
Frankowski, Marek; Chęciński, Jakub; Czapkiewicz, Maciej
2015-04-01
We present a tool for calculations of Fourier transform spatial distribution taken from magnetization dynamics simulated in Object Oriented Micromagnetic Framework (OOMMF). In OOMMF, as well as in other popular micromagnetic software, output data is organized as magnetization vectors from each simulation cell written down to separate file for each simulation step. Therefore, we use parallel computations to reorganize data in files containing time evolution for each cell. Fast Fourier transform is obtained for selected time period by parallel computations using Matlab. The output is a spatial distribution of the magnitude for the selected frequency in the sample cross-section. It allows for analysis of spin waves localization and therefore helps to understand their origin in investigated sample.
NASA Technical Reports Server (NTRS)
Cecil, R. W.; White, R. A.; Szczur, M. R.
1972-01-01
The IDAMS Processor is a package of task routines and support software that performs convolution filtering, image expansion, fast Fourier transformation, and other operations on a digital image tape. A unique task control card for that program, together with any necessary parameter cards, selects each processing technique to be applied to the input image. A variable number of tasks can be selected for execution by including the proper task and parameter cards in the input deck. An executive maintains control of the run; it initiates execution of each task in turn and handles any necessary error processing.
Komorowski, Dariusz; Pietraszek, Stanislaw
2016-01-01
This paper presents the analysis of multi-channel electrogastrographic (EGG) signals using the continuous wavelet transform based on the fast Fourier transform (CWTFT). The EGG analysis was based on the determination of the several signal parameters such as dominant frequency (DF), dominant power (DP) and index of normogastria (NI). The use of continuous wavelet transform (CWT) allows for better visible localization of the frequency components in the analyzed signals, than commonly used short-time Fourier transform (STFT). Such an analysis is possible by means of a variable width window, which corresponds to the scale time of observation (analysis). Wavelet analysis allows using long time windows when we need more precise low-frequency information, and shorter when we need high frequency information. Since the classic CWT transform requires considerable computing power and time, especially while applying it to the analysis of long signals, the authors used the CWT analysis based on the fast Fourier transform (FFT). The CWT was obtained using properties of the circular convolution to improve the speed of calculation. This method allows to obtain results for relatively long records of EGG in a fairly short time, much faster than using the classical methods based on running spectrum analysis (RSA). In this study authors indicate the possibility of a parametric analysis of EGG signals using continuous wavelet transform which is the completely new solution. The results obtained with the described method are shown in the example of an analysis of four-channel EGG recordings, performed for a non-caloric meal. PMID:26573647
NASA Astrophysics Data System (ADS)
Mao, Xi-an; Li, Ning; Mao, Jiezhen; Li, Qiurong; Xiao, Nan; Jiang, Bin; Jiang, Ling; Wang, Xu-xia; Liu, Maili
2012-01-01
Detection and quantification of total choline-containing metabolites (CCMs) in tissues by magnetic resonance spectroscopy (MRS) has received considerable attention as a biomarker of cancer. Tissue CCMs are mainly choline (Cho), phosphocholine (PCho), and glycerophosphocholine (GPCho). Because the methyl 1H resonances of tissue CCMs exhibit small chemical shift differences and overlap significantly in 1D 1H MRS, quantification of individual components is precluded. Development of a MRS method capably of resolving individual components of tissue CCMs would be a significant advance. Herein, a modification of the 2D 1H- 14N HSQC technique is targeted on the two methylene 1H in the CH 2O group ( 3J1H14N = 2.7 Hz) and applied to ex vivo mouse and human liver samples at physiological temperature (37 °C). Specifically, the 1H- 14N HSQC technique is modified into a 2D 1H- 14N three-bond correlation (HN3BC) experiment, which selectively detects the 1H of CH 2O coupled to 14N in CCMs. Separate signals from Cho, PCho, and GPCho components are resolved with high detection sensitivity. A 2D HN3BC spectrum can be recorded from mouse liver in only 1.5 min and from human carcinoma liver tissue in less than 3 min with effective sample volume of 0.2 ml at 14.1 T.
NASA Technical Reports Server (NTRS)
Boriakoff, Valentin
1994-01-01
The goal of this project was the feasibility study of a particular architecture of a digital signal processing machine operating in real time which could do in a pipeline fashion the computation of the fast Fourier transform (FFT) of a time-domain sampled complex digital data stream. The particular architecture makes use of simple identical processors (called inner product processors) in a linear organization called a systolic array. Through computer simulation the new architecture to compute the FFT with systolic arrays was proved to be viable, and computed the FFT correctly and with the predicted particulars of operation. Integrated circuits to compute the operations expected of the vital node of the systolic architecture were proven feasible, and even with a 2 micron VLSI technology can execute the required operations in the required time. Actual construction of the integrated circuits was successful in one variant (fixed point) and unsuccessful in the other (floating point).
NASA Astrophysics Data System (ADS)
Paliwal, Deepak; Choudhur, Achintya; Govandhan, T.
2014-06-01
Fault diagnosis of rolling element bearings requires efficient signal processing techniques. For this purpose, the performances of envelope detection with fast Fourier transform (FFT) and continuous wavelet transform (CWT) of vibration signals produced from a bearing with defects on inner race and rolling element, have been examined at low signal to noise ratio. Both simulated and experimental signals from identical bearings have been considered for the purpose of analysis. The bearings have been modeled as spring-mass-dashpot systems and the simulated signals have been obtained considering transfer functions for the bearing systems subjected to impulsive loads due to the defects. Frequency B spline wavelets have been applied for CWT and a discussion on wavelet selection has been presented for better effectiveness. Results show that use of CWT with the proposed wavelets overcomes the short coming of FFT while processing a noisy vibration signals for defect detection of bearings.
Texture Analysis of the Epidermis Based on Fast Fourier Transformation in Sjögren-Larsson Syndrome
Auada, Mariam P; Adam, Randall L.; Leite, Neucimar J.; Puzzi, Maria B.; Cintra, Maria L.; Rizzo, William B.; Metze, Konradin
2011-01-01
OBJECTIVE To investigate whether image analysis of routine hematoxylin-eosin (H-E) skin sections usingfast Fourier transformation (FFT) could detect structural alterations in patients with Sjögren-Larsson syndrome (SLS) diagnosed by molecular biology. STUDY DESIGN Skin punch biopsies of 9 patients with SLS and 17 healthy volunteers were obtained. Digital images of routine histologic sections were taken, and their gray scale luminance was analyzed by FFT. The inertia values were determined for different ranges of the spatial frequencies in the vertical and horizontal direction. To get an estimation of anisotropy, we calculated the resultant vector of the designated frequency ranges. RESULTS In the prickle cell layer, SLS patients showed more intense amplitudes in spatial structures with periods between 1.2 and 3.6 µm in the vertical direction, which correlated in part with accentuated nuclei and nucleoli and perinucleolar halos in the H-E sections. In a linear discriminant analysis, the variables derived from the FFT images correctly discriminated 84.6% of the patients. Texture features derived from the gray level co-occurrence matrix were not able to separate the groups. CONCLUSION Exploratory texture analysis by FFT was able to detect discrete alterations in the prickle cell layer in routine light microscopy slides of SLS patients. The structural changes identified by FFT may be related to abnormal cellular components associated with aberrant lipid metabolism. PMID:16927642
Utgenannt, André; Maspero, Ross; Fortini, Andrea; Turner, Rebecca; Florescu, Marian; Jeynes, Christopher; Kanaras, Antonios G; Muskens, Otto L; Sear, Richard P; Keddie, Joseph L
2016-02-23
When fabricating photonic crystals from suspensions in volatile liquids using the horizontal deposition method, the conventional approach is to evaporate slowly to increase the time for particles to settle in an ordered, periodic close-packed structure. Here, we show that the greatest ordering of 10 nm aqueous gold nanoparticles (AuNPs) in a template of larger spherical polymer particles (mean diameter of 338 nm) is achieved with very fast water evaporation rates obtained with near-infrared radiative heating. Fabrication of arrays over areas of a few cm(2) takes only 7 min. The assembly process requires that the evaporation rate is fast relative to the particles' Brownian diffusion. Then a two-dimensional colloidal crystal forms at the falling surface, which acts as a sieve through which the AuNPs pass, according to our Langevin dynamics computer simulations. With sufficiently fast evaporation rates, we create a hybrid structure consisting of a two-dimensional AuNP nanoarray (or "nanogrid") on top of a three-dimensional polymer opal. The process is simple, fast, and one-step. The interplay between the optical response of the plasmonic Au nanoarray and the microstructuring of the photonic opal results in unusual optical spectra with two extinction peaks, which are analyzed via finite-difference time-domain method simulations. Comparison between experimental and modeling results reveals a strong interplay of plasmonic modes and collective photonic effects, including the formation of a high-order stopband and slow-light-enhanced plasmonic absorption. The structures, and hence their optical signatures, are tuned by adjusting the evaporation rate via the infrared power density. PMID:26767891
Westphalen, Antonio C.; Noworolski, Susan M.; Harisinghani, Mukesh; Jhaveri, Kartik S.; Raman, Steve S.; Rosenkrantz, Andrew B.; Wang, Zhen J.; Zagoria, Ronald J.; Kurhanewicz, John
2016-01-01
OBJECTIVE The goal of this study was to compare the perceived quality of 3-T axial T2-weighted high-resolution 2D and high-resolution 3D fast spin-echo (FSE) endorectal MR images of the prostate. MATERIALS AND METHODS Six radiologists independently reviewed paired 3-T axial T2-weighted high-resolution 2D and 3D FSE endorectal MR images of the prostates of 85 men in two sessions. In the first session (n = 85), each reader selected his or her preferred images; in the second session (n = 28), they determined their confidence in tumor identification and compared the depiction of the prostatic anatomy, tumor conspicuity, and subjective intrinsic image quality of images. A meta-analysis using a random-effects model, logistic regression, and the paired Wilcoxon rank-sum test were used for statistical analyses. RESULTS Three readers preferred the 2D acquisition (67–89%), and the other three preferred the 3D images (70–80%). The option for one of the techniques was not associated with any of the predictor variables. The 2D FSE images were significantly sharper than 3D FSE (p < 0.001) and significantly more likely to exhibit other (nonmotion) artifacts (p = 0.002). No other statistically significant differences were found. CONCLUSION Our results suggest that there are strong individual preferences for the 2D or 3D FSE MR images, but there was a wide variability among radiologists. There were differences in image quality (image sharpness and presence of artifacts not related to motion) but not in the sequences’ ability to delineate the glandular anatomy and depict a cancerous tumor. PMID:26491891
PixFEL: developing a fine pitch, fast 2D X-ray imager for the next generation X-FELs
NASA Astrophysics Data System (ADS)
Ratti, L.; Comotti, D.; Fabris, L.; Grassi, M.; Lodola, L.; Malcovati, P.; Manghisoni, M.; Re, V.; Traversi, G.; Vacchi, C.; Bettarini, S.; Casarosa, G.; Forti, F.; Morsani, F.; Paladino, A.; Paoloni, E.; Rizzo, G.; Benkechkache, M. A.; Dalla Betta, G.-F.; Mendicino, R.; Pancheri, L.; Verzellesi, G.; Xu, H.
2015-10-01
The PixFEL project is conceived as the first stage of a long term research program aiming at the development of advanced X-ray imaging instrumentation for applications at the free electron laser (FEL) facilities. The project aims at substantially advancing the state-of-the-art in the field of 2D X-ray imaging by exploring cutting-edge solutions for sensor development, for integration processes and for readout channel architectures. The main focus is on the development of the fundamental microelectronic building blocks for detector readout and on the technologies for the assembly of a multilayer module with minimum dead area. This work serves the purpose of introducing the main features of the project, together with the simulation results leading to the first prototyping run.
Nikolic, Goran; Zlatkovic, Sasa; Cakic, Milorad; Cakic, Suzana; Lacnjevac, Caslav; Rajic, Zoran
2010-01-01
The use of fast FT-IR spectroscopy as a sensitive method to estimate a change of the crosslinking kinetics of epoxy resin with polyamine adducts is described in this study. A new epoxy formulation based on the use of polyamine adducts as the hardeners was analyzed. Crosslinking reactions of the different stoichiometric mixtures of the unmodified GY250 epoxy resin with the aliphatic EH606 and the cycloaliphatic EH637 polyamine adducts were studied using mid FT-IR spectroscopic techniques. As the crosslinking proceeded, the primary amine groups in polyamine adduct are converted to secondary and the tertiary amines. The decrease in the IR band intensity of epoxy groups at about 915 cm−1, as well as at about 3,056 cm−1, was observed due to process. Mid IR spectral analysis was used to calculate the content of the epoxy groups as a function of crosslinking time and the crosslinking degree of resin. The amount of all the epoxy species was estimated from IR spectra to changes during the crosslinking kinetics of epichlorhydrin. PMID:22315562
Kuo-Petravic, G.
1988-12-01
Electron cyclotron emission measurements taken at 20 locations in the horizontal midplane during a sawtooth crash have been analysed based on the assumption of fast rigid rotation of the plasma. Due to this fast rotation (approx.100..mu..sec), which remains fairly constant throughout the sawtooth crash, we have been able to make time-to-space reconstructions of half the poloidal plane using points which are separated in time by not more than 40..mu..sec. The existence of a temperature flattening in the precursor phase, which we interpret as an m = 1 temperature island, is clearly demonstrated, and its location and width agree well with local emissivity measurements from soft x-ray tomography viewing the same poloidal plane. The rotating temperature island in the precursor phase, the outward movement of the region of high T/sub c/ during the crash phase, and the shape of T/sub e/ during the crash phase, and the shape of T/sub e/ distribution after the crash during the successor phase have all been documented in a time sequence of color contours. 4 refs., 10 figs.
Adams, Robert; Zboray, Robert; Cortesi, Marco; Prasser, Horst-Michael
2014-04-01
A conceptual design optimization of a fast neutron tomography system was performed. The system is based on a compact deuterium-deuterium fast neutron generator and an arc-shaped array of individual neutron detectors. The array functions as a position sensitive one-dimensional detector allowing tomographic reconstruction of a two-dimensional cross section of an object up to 10 cm across. Each individual detector is to be optically isolated and consists of a plastic scintillator and a Silicon Photomultiplier for measuring light produced by recoil protons. A deterministic geometry-based model and a series of Monte Carlo simulations were used to optimize the design geometry parameters affecting the reconstructed image resolution. From this, it is expected that with an array of 100 detectors a reconstructed image resolution of ~1.5mm can be obtained. Other simulations were performed in order to optimize the scintillator depth (length along the neutron path) such that the best ratio of direct to scattered neutron counts is achieved. This resulted in a depth of 6-8 cm and an expected detection efficiency of 33-37%. Based on current operational capabilities of a prototype neutron generator being developed at the Paul Scherrer Institute, planned implementation of this detector array design should allow reconstructed tomograms to be obtained with exposure times on the order of a few hours. PMID:24495568
NASA Astrophysics Data System (ADS)
Torabi, M.; Mousavi G., S. M.; Younesian, D.
2015-03-01
Registration of the point clouds is a conventional challenge in computer vision related applications. As an application, matching of train wheel profiles extracted from two viewpoints is studied in this paper. The registration problem is formulated into an optimization problem. An error minimization function for registration of the two partially overlapping point clouds is presented. The error function is defined as the sum of the squared distance between the source points and their corresponding pairs which should be minimized. The corresponding pairs are obtained thorough Iterative Closest Point (ICP) variants. Here, a point-to-plane ICP variant is employed. Principal Component Analysis (PCA) is used to obtain tangent planes. Thus it is shown that minimization of the proposed objective function diminishes point-to-plane ICP variant. We utilized this algorithm to register point clouds of two partially overlapping profiles of wheel train extracted from two viewpoints in 2D. Also, a number of synthetic point clouds and a number of real point clouds in 3D are studied to evaluate the reliability and rate of convergence in our method compared with other registration methods.
Reyhan, Meral; Kim, Hyun J.; Brown, Matthew S.; Ennis, Daniel B.
2013-01-01
Purpose To assess the intra- and inter-scan reproducibility of LV twist using FAST. Assessing the reproducibility of the measurement of new magnetic resonance imaging (MRI) biomarkers is an important part of validation. Fourier Analysis of STimulated Echoes (FAST) is a new MRI tissue tagging method that has recently been shown to compare favorably to conventional estimates of left ventricular (LV) twist from cardiac tagged images, but with significantly reduced user interaction time. Materials and Methods Healthy volunteers (N=10) were scanned twice using FAST over one week. On Day-1 two measurements of LV twist were collected for intra-scan comparisons. Measurements for LV twist were again collected on Day-8 for inter-scan assessment. LV short-axis tagged images were acquired on a 3T scanner in order to ensure detectability of tags during early and mid-diastole. Peak LV twist is reported as mean±SD. Reproducibility was assessed using the concordance correlation coefficient (CCC) and the repeatability coefficient (RC) (95%-CI range). Results Mean peak twist measurements were 13.4±4.3° (Day-1, Scan-1), 13.6±3.7° (Day-1, Scan-2), and 13.0±2.7° (Day-8). Bland-Altman analysis resulted in intra- and inter-scan bias and 95%-CI of −0.6° [−1.0°, 1.6°] and 1.4° [−1.0°, 3.0°], respectively. The Bland-Altman RC for peak LV twist was 2.6° and 4.0° for intra- and inter-scan respectively. The CCC was 0.9 and 0.6 for peak LV twist for intra- and inter-scan respectively. Conclusion FAST is a semi-automated method that provides a quick and quantitative assessment of LV systolic and diastolic twist that demonstrates high intra-scan and moderate inter-scan reproducibility in preliminary studies. PMID:23633244
Beaunier, L.; Frydman, J.; Gabrielli, C.; Huet, F.; Keddam, M.
1996-12-31
A comparison of a spectral analysis using the fast Fourier transform (FFT) and the maximum entropy method (MEM) was carried out in the case in which both methods can be performed, that is, when several time acquisitions can be recorded. A summary of the principles of the MEM is given. Then the main properties of this method are investigated, that is, influence of the MEM order on the spectrum accuracy, validity of the low-frequency plateau usually given by this technique, overlapping of spectra measured for different frequency bandwidths, and influence of a slow evolution of the amplitude of the signal fluctuations. The susceptibility to pitting corrosion of type 304 stainless steel and type 304 modified by molybdenum (Mo) by means of ion implantation was studied. The power spectral densities (PSD) measured with the FFT and MEM techniques are in reasonable agreement, except for low electrochemical current noises (ECN) buried in the parasitic noise generated by the power supply. In that case, the FFT technique is more appropriate than the MEM, which gave qualitative results only. The type 304 stainless steel showed a large metastable pitting leading to only a few macroscopic pits, whereas the type 304 Mo-implanted specimen showed a very low metastable pitting leading to many hemispheric pits covered by the Mo-implanted layer, under which localized corrosion occurred.
NASA Astrophysics Data System (ADS)
Hu, Hongda; Shu, Hong; Hu, Zhiyong; Xu, Jianhui
2016-04-01
Kriging interpolation provides the best linear unbiased estimation for unobserved locations, but its heavy computation limits the manageable problem size in practice. To address this issue, an efficient interpolation procedure incorporating the fast Fourier transform (FFT) was developed. Extending this efficient approach, we propose an FFT-based parallel algorithm to accelerate regression Kriging interpolation on an NVIDIA® compute unified device architecture (CUDA)-enabled graphic processing unit (GPU). A high-performance cuFFT library in the CUDA toolkit was introduced to execute computation-intensive FFTs on the GPU, and three time-consuming processes were redesigned as kernel functions and executed on the CUDA cores. A MODIS land surface temperature 8-day image tile at a resolution of 1 km was resampled to create experimental datasets at eight different output resolutions. These datasets were used as the interpolation grids with different sizes in a comparative experiment. Experimental results show that speedup of the FFT-based regression Kriging interpolation accelerated by GPU can exceed 1000 when processing datasets with large grid sizes, as compared to the traditional Kriging interpolation running on the CPU. These results demonstrate that the combination of FFT methods and GPU-based parallel computing techniques greatly improves the computational performance without loss of precision.
NASA Technical Reports Server (NTRS)
Boriakoff, Valentin; Chen, Wei
1990-01-01
The NASA-Cornell Univ.-Worcester Polytechnic Institute Fast Fourier Transform (FFT) chip based on the architecture of the systolic FFT computation as presented by Boriakoff is implemented into an operating device design. The kernel of the system, a systolic inner product floating point processor, was designed to be assembled into a systolic network that would take incoming data streams in pipeline fashion and provide an FFT output at the same rate, word by word. It was thoroughly simulated for proper operation, and it has passed a comprehensive set of tests showing no operational errors. The black box specifications of the chip, which conform to the initial requirements of the design as specified by NASA, are given. The five subcells are described and their high level function description, logic diagrams, and simulation results are presented. Some modification of the Read Only Memory (ROM) design were made, since some errors were found in it. Because a four stage pipeline structure was used, simulating such a structure is more difficult than an ordinary structure. Simulation methods are discussed. Chip signal protocols and chip pinout are explained.
NASA Astrophysics Data System (ADS)
Quinten, M.; Houta, F.; Fries, T.
2015-06-01
Thin film thickness determination with a reflectometer suffers from two problems. One problem is the leakage in the Fast Fourier Transform caused by the fact that the two variables wavenumber 1/λ and optical thickness nṡd are not really independent, since the refractive index n of the film material itself depends upon the wavenumber. This causes uncertainties in the thickness determination in the order of up to 5% for highly refractive materials like semiconductors. We present a simple but effective improvement of this contribution of the leakage that reduces the uncertainty to less than 2% for highly refractive materials. Another problem that mainly affects thin films below about 2 μm arises if one uses measuring heads collimators or even microscope headers to obtain high lateral resolutions in the thickness determination. The use of a header introduces angles of incidence different from the default angle α = 0° in reflectometry. Then, the measured reflectance becomes polarization-dependent and the angle must be explicitly considered in the evaluation algorithm. For a microscope header however, all angles between 0° and the angle of aperture must be considered. We will present a solution that allows to reduce the work for each header on taking into account the polarization of the reflected light and a corresponding effective angle αeff.
NASA Astrophysics Data System (ADS)
Itoh, Taketsugu; Yamauchi, Noriyoshi
2007-05-01
Surface morphology of pentacene thin films and their substrates with under-layers is characterized by using atomic force microscopy (AFM). The power values of power spectral density (PSD) for the AFM digital data were determined by the fast Fourier transform (FFT) algorithms instead of the root-mean-square (rms) and peak-to-valley value. The PSD plots of pentacene films on glass substrate are successfully approximated by the k-correlation model. The pentacene film growth is interpreted the intermediation of the bulk and surface diffusion by parameter C of k-correlation model. The PSD plots of pentacene film on Au under-layer is approximated by using the linear continuum model (LCM) instead of the combination model of the k-correlation model and Gaussian function. The PSD plots of SiO 2 layer on Au under-layer as a gate insulator on a gate electrode of organic thin film transistors (OTFTs) have three power values of PSD. It is interpreted that the specific three PSD power values are caused by the planarization of the smooth SiO 2 layer to rough Au under-layer.
Jafari, Safiye; Faridbod, Farnoush; Norouzi, Parviz; Dezfuli, Amin Shiralizadeh; Ajloo, Davood; Mohammadipanah, Fatemeh; Ganjali, Mohammad Reza
2015-10-01
A new strategy was introduced for ssDNA immobilization on a modified glassy carbon electrode. The electrode surface was modified using polyaniline and chemically reduced graphene oxide decorated cerium oxide nanoparticles (CeO2NPs-RGO). A single-stranded DNA (ssDNA) probe was immobilized on the modified electrode surface. Fast Fourier transform square wave voltammetry (FFT-SWV) was applied as detection technique and [Ru(bpy)3](2+/3+) redox signal was used as electrochemical marker. The hybridization of ssDNA with its complementary target caused a dramatic decrease in [Ru(bpy)3](2+/3+) FFT-SW signal. The proposed electrochemical biosensor was able to detect Aeromonas hydrophila DNA oligonucleotide sequence encoding aerolysin protein. Under optimal conditions, the biosensor showed excellent selectivity toward complementary sequence in comparison with noncomplementary and two-base mismatch sequences. The dynamic linear range of this electrochemical DNA biosensor for detecting 20-mer oligonucleotide sequence of A. hydrophila was from 1 × 10(-15) to 1 × 10(-8) mol L(-1). The proposed biosensor was successfully applied for the detection of DNA extracted from A. hydrophila in fish pond water up to 0.01 μg mL(-1) with RSD of 5%. Besides, molecular docking was applied to consider the [Ru(bpy)3](2+/3+) interaction with ssDNA before and after hybridization. PMID:26454462
Inertial solvation in femtosecond 2D spectra
NASA Astrophysics Data System (ADS)
Hybl, John; Albrecht Ferro, Allison; Farrow, Darcie; Jonas, David
2001-03-01
We have used 2D Fourier transform spectroscopy to investigate polar solvation. 2D spectroscopy can reveal molecular lineshapes beneath ensemble averaged spectra and freeze molecular motions to give an undistorted picture of the microscopic dynamics of polar solvation. The transition from "inhomogeneous" to "homogeneous" 2D spectra is governed by both vibrational relaxation and solvent motion. Therefore, the time dependence of the 2D spectrum directly reflects the total response of the solvent-solute system. IR144, a cyanine dye with a dipole moment change upon electronic excitation, was used to probe inertial solvation in methanol and propylene carbonate. Since the static Stokes' shift of IR144 in each of these solvents is similar, differences in the 2D spectra result from solvation dynamics. Initial results indicate that the larger propylene carbonate responds more slowly than methanol, but appear to be inconsistent with rotational estimates of the inertial response. To disentangle intra-molecular vibrations from solvent motion, the 2D spectra of IR144 will be compared to the time-dependent 2D spectra of the structurally related nonpolar cyanine dye HDITCP.
Energy Science and Technology Software Center (ESTSC)
2005-07-01
Aniso2d is a two-dimensional seismic forward modeling code. The earth is parameterized by an X-Z plane in which the seismic properties Can have monoclinic with x-z plane symmetry. The program uses a user define time-domain wavelet to produce synthetic seismograms anrwhere within the two-dimensional media.
NASA Astrophysics Data System (ADS)
Cai, Ailong; Wang, Linyuan; Yan, Bin; Zhang, Hanming; Li, Lei; Xi, Xiaoqi; Li, Jianxin
2015-03-01
In this study, we consider a novel form of computed tomography (CT), that is, linear scan CT (LCT), which applies a straight line trajectory. Furthermore, an iterative algorithm is proposed for pseudo-polar Fourier reconstruction through total variation minimization (PPF-TVM). Considering that the sampled Fourier data are distributed in pseudo-polar coordinates, the reconstruction model minimizes the TV of the image subject to the constraint that the estimated 2D Fourier data for the image are consistent with the 1D Fourier transform of the projection data. PPF-TVM employs the alternating direction method (ADM) to develop a robust and efficient iteration scheme, which ensures stable convergence provided that appropriate parameter values are given. In the ADM scheme, PPF-TVM applies the pseudo-polar fast Fourier transform and its adjoint to iterate back and forth between the image and frequency domains. Thus, there is no interpolation in the Fourier domain, which makes the algorithm both fast and accurate. PPF-TVM is particularly useful for limited angle reconstruction in LCT and it appears to be robust against artifacts. The PPF-TVM algorithm was tested with the FORBILD head phantom and real data in comparisons with state-of-the-art algorithms. Simulation studies and real data verification suggest that PPF-TVM can reconstruct higher accuracy images with lower time consumption.
Fast computer simulation of reconstructed image from rainbow hologram based on GPU
NASA Astrophysics Data System (ADS)
Shuming, Jiao; Yoshikawa, Hiroshi
2015-10-01
A fast computer simulation solution for rainbow hologram reconstruction based on GPU is proposed. In the commonly used segment Fourier transform method for rainbow hologram reconstruction, the computation of 2D Fourier transform on each hologram segment is very time consuming. GPU-based parallel computing can be applied to improve the computing speed. Compared with CPU computing, simulation results indicate that our proposed GPU computing can effectively reduce the computation time by as much as eight times.
Greg Flach, Frank Smith
2011-12-31
Mesh2d is a Fortran90 program designed to generate two-dimensional structured grids of the form [x(i),y(i,j)] where [x,y] are grid coordinates identified by indices (i,j). The x(i) coordinates alone can be used to specify a one-dimensional grid. Because the x-coordinates vary only with the i index, a two-dimensional grid is composed in part of straight vertical lines. However, the nominally horizontal y(i,j0) coordinates along index i are permitted to undulate or otherwise vary. Mesh2d also assigns an integer material type to each grid cell, mtyp(i,j), in a user-specified manner. The complete grid is specified through three separate input files defining the x(i), y(i,j), and mtyp(i,j) variations.
Energy Science and Technology Software Center (ESTSC)
2011-12-31
Mesh2d is a Fortran90 program designed to generate two-dimensional structured grids of the form [x(i),y(i,j)] where [x,y] are grid coordinates identified by indices (i,j). The x(i) coordinates alone can be used to specify a one-dimensional grid. Because the x-coordinates vary only with the i index, a two-dimensional grid is composed in part of straight vertical lines. However, the nominally horizontal y(i,j0) coordinates along index i are permitted to undulate or otherwise vary. Mesh2d also assignsmore » an integer material type to each grid cell, mtyp(i,j), in a user-specified manner. The complete grid is specified through three separate input files defining the x(i), y(i,j), and mtyp(i,j) variations.« less
NASA Astrophysics Data System (ADS)
Lotsch, Bettina V.
2015-07-01
Graphene's legacy has become an integral part of today's condensed matter science and has equipped a whole generation of scientists with an armory of concepts and techniques that open up new perspectives for the postgraphene area. In particular, the judicious combination of 2D building blocks into vertical heterostructures has recently been identified as a promising route to rationally engineer complex multilayer systems and artificial solids with intriguing properties. The present review highlights recent developments in the rapidly emerging field of 2D nanoarchitectonics from a materials chemistry perspective, with a focus on the types of heterostructures available, their assembly strategies, and their emerging properties. This overview is intended to bridge the gap between two major—yet largely disjunct—developments in 2D heterostructures, which are firmly rooted in solid-state chemistry or physics. Although the underlying types of heterostructures differ with respect to their dimensions, layer alignment, and interfacial quality, there is common ground, and future synergies between the various assembly strategies are to be expected.
Efficient single pixel imaging in Fourier space
NASA Astrophysics Data System (ADS)
Bian, Liheng; Suo, Jinli; Hu, Xuemei; Chen, Feng; Dai, Qionghai
2016-08-01
Single pixel imaging (SPI) is a novel technique capturing 2D images using a bucket detector with a high signal-to-noise ratio, wide spectrum range and low cost. Conventional SPI projects random illumination patterns to randomly and uniformly sample the entire scene’s information. Determined by Nyquist sampling theory, SPI needs either numerous projections or high computation cost to reconstruct the target scene, especially for high-resolution cases. To address this issue, we propose an efficient single pixel imaging technique (eSPI), which instead projects sinusoidal patterns for importance sampling of the target scene’s spatial spectrum in Fourier space. Specifically, utilizing the centrosymmetric conjugation and sparsity priors of natural images’ spatial spectra, eSPI sequentially projects two \\tfrac{π }{2}-phase-shifted sinusoidal patterns to obtain each Fourier coefficient in the most informative spatial frequency bands. eSPI can reduce requisite patterns by two orders of magnitude compared to conventional SPI, which helps a lot for fast and high-resolution SPI.
NASA Astrophysics Data System (ADS)
Li, Jinghe; Song, Linping; Liu, Qing Huo
2016-02-01
A simultaneous multiple frequency contrast source inversion (CSI) method is applied to reconstructing hydrocarbon reservoir targets in a complex multilayered medium in two dimensions. It simulates the effects of a salt dome sedimentary formation in the context of reservoir monitoring. In this method, the stabilized biconjugate-gradient fast Fourier transform (BCGS-FFT) algorithm is applied as a fast solver for the 2D volume integral equation for the forward computation. The inversion technique with CSI combines the efficient FFT algorithm to speed up the matrix-vector multiplication and the stable convergence of the simultaneous multiple frequency CSI in the iteration process. As a result, this method is capable of making quantitative conductivity image reconstruction effectively for large-scale electromagnetic oil exploration problems, including the vertical electromagnetic profiling (VEP) survey investigated here. A number of numerical examples have been demonstrated to validate the effectiveness and capacity of the simultaneous multiple frequency CSI method for a limited array view in VEP.
Dohet-Eraly, Jérôme; Yourassowsky, Catherine; Dubois, Frank
2016-09-01
The knowledge of the complex amplitude of optical fields, that is, both quantitative phase and intensity, enables numeric reconstruction along the optical axis. Nonetheless, a criterion is required for autofocusing. This Letter presents a robust and rapid refocusing criterion suitable for color interferometric digital holographic microscopy, and, more generally, for applications where complex amplitude is known for at least two different wavelengths. This criterion uses the phase in the Fourier domain, which is compared among wavelengths. It is applicable whatever the nature of the observed object: opaque, refractive, or both mixed. The method is validated with simulated and experimental holograms. PMID:27607975
NASA Astrophysics Data System (ADS)
Thomas, Daniel D.; Zhang, Yuhong
2014-06-01
Along-track interferometry (ATI) is used to detect ground moving targets against a stationary background in synthetic aperture radar (SAR) imagery. In this paper, we present a novel approach to multi-channel ATI wherein clutter cancellation is applied to each pixel of the multiple SAR images, followed by a Fourier transform to estimate range rate (Doppler). Range rate estimates allow us to compensate for the cross-range offset of the target, thus geo-locating the targets. We then present a number of benefits to this approach.
Fourier Lucas-Kanade algorithm.
Lucey, Simon; Navarathna, Rajitha; Ashraf, Ahmed Bilal; Sridharan, Sridha
2013-06-01
In this paper, we propose a framework for both gradient descent image and object alignment in the Fourier domain. Our method centers upon the classical Lucas & Kanade (LK) algorithm where we represent the source and template/model in the complex 2D Fourier domain rather than in the spatial 2D domain. We refer to our approach as the Fourier LK (FLK) algorithm. The FLK formulation is advantageous when one preprocesses the source image and template/model with a bank of filters (e.g., oriented edges, Gabor, etc.) as 1) it can handle substantial illumination variations, 2) the inefficient preprocessing filter bank step can be subsumed within the FLK algorithm as a sparse diagonal weighting matrix, 3) unlike traditional LK, the computational cost is invariant to the number of filters and as a result is far more efficient, and 4) this approach can be extended to the Inverse Compositional (IC) form of the LK algorithm where nearly all steps (including Fourier transform and filter bank preprocessing) can be precomputed, leading to an extremely efficient and robust approach to gradient descent image matching. Further, these computational savings translate to nonrigid object alignment tasks that are considered extensions of the LK algorithm, such as those found in Active Appearance Models (AAMs). PMID:23599053
A unified Fourier theory for time-of-flight PET data
NASA Astrophysics Data System (ADS)
Li, Yusheng; Matej, Samuel; Metzler, Scott D.
2016-01-01
Fully 3D time-of-flight (TOF) PET scanners offer the potential of previously unachievable image quality in clinical PET imaging. TOF measurements add another degree of redundancy for cylindrical PET scanners and make photon-limited TOF-PET imaging more robust than non-TOF PET imaging. The data space for 3D TOF-PET data is five-dimensional with two degrees of redundancy. Previously, consistency equations were used to characterize the redundancy of TOF-PET data. In this paper, we first derive two Fourier consistency equations and Fourier-John equation for 3D TOF PET based on the generalized projection-slice theorem; the three partial differential equations (PDEs) are the dual of the sinogram consistency equations and John’s equation. We then solve the three PDEs using the method of characteristics. The two degrees of entangled redundancy of the TOF-PET data can be explicitly elicited and exploited by the solutions of the PDEs along the characteristic curves, which gives a complete understanding of the rich structure of the 3D x-ray transform with TOF measurement. Fourier rebinning equations and other mapping equations among different types of PET data are special cases of the general solutions. We also obtain new Fourier rebinning and consistency equations (FORCEs) from other special cases of the general solutions, and thus we obtain a complete scheme to convert among different types of PET data: 3D TOF, 3D non-TOF, 2D TOF and 2D non-TOF data. The new FORCEs can be used as new Fourier-based rebinning algorithms for TOF-PET data reduction, inverse rebinnings for designing fast projectors, or consistency conditions for estimating missing data. Further, we give a geometric interpretation of the general solutions—the two families of characteristic curves can be obtained by respectively changing the azimuthal and co-polar angles of the biorthogonal coordinates in Fourier space. We conclude the unified Fourier theory by showing that the Fourier consistency equations
A unified Fourier theory for time-of-flight PET data.
Li, Yusheng; Matej, Samuel; Metzler, Scott D
2016-01-21
Fully 3D time-of-flight (TOF) PET scanners offer the potential of previously unachievable image quality in clinical PET imaging. TOF measurements add another degree of redundancy for cylindrical PET scanners and make photon-limited TOF-PET imaging more robust than non-TOF PET imaging. The data space for 3D TOF-PET data is five-dimensional with two degrees of redundancy. Previously, consistency equations were used to characterize the redundancy of TOF-PET data. In this paper, we first derive two Fourier consistency equations and Fourier-John equation for 3D TOF PET based on the generalized projection-slice theorem; the three partial differential equations (PDEs) are the dual of the sinogram consistency equations and John's equation. We then solve the three PDEs using the method of characteristics. The two degrees of entangled redundancy of the TOF-PET data can be explicitly elicited and exploited by the solutions of the PDEs along the characteristic curves, which gives a complete understanding of the rich structure of the 3D x-ray transform with TOF measurement. Fourier rebinning equations and other mapping equations among different types of PET data are special cases of the general solutions. We also obtain new Fourier rebinning and consistency equations (FORCEs) from other special cases of the general solutions, and thus we obtain a complete scheme to convert among different types of PET data: 3D TOF, 3D non-TOF, 2D TOF and 2D non-TOF data. The new FORCEs can be used as new Fourier-based rebinning algorithms for TOF-PET data reduction, inverse rebinnings for designing fast projectors, or consistency conditions for estimating missing data. Further, we give a geometric interpretation of the general solutions--the two families of characteristic curves can be obtained by respectively changing the azimuthal and co-polar angles of the biorthogonal coordinates in Fourier space. We conclude the unified Fourier theory by showing that the Fourier consistency equations are
Fourier rebinning and consistency equations for time-of-flight PET planograms
NASA Astrophysics Data System (ADS)
Li, Yusheng; Defrise, Michel; Matej, Samuel; Metzler, Scott D.
2016-09-01
Due to the unique geometry, dual-panel PET scanners have many advantages in dedicated breast imaging and on-board imaging applications since the compact scanners can be combined with other imaging and treatment modalities. The major challenges of dual-panel PET imaging are the limited-angle problem and data truncation, which can cause artifacts due to incomplete data sampling. The time-of-flight (TOF) information can be a promising solution to reduce these artifacts. The TOF planogram is the native data format for dual-panel TOF PET scanners, and the non-TOF planogram is the 3D extension of linogram. The TOF planograms is five-dimensional while the objects are three-dimensional, and there are two degrees of redundancy. In this paper, we derive consistency equations and Fourier-based rebinning algorithms to provide a complete understanding of the rich structure of the fully 3D TOF planograms. We first derive two consistency equations and John’s equation for 3D TOF planograms. By taking the Fourier transforms, we obtain two Fourier consistency equations (FCEs) and the Fourier–John equation (FJE), which are the duals of the consistency equations and John’s equation, respectively. We then solve the FCEs and FJE using the method of characteristics. The two degrees of entangled redundancy of the 3D TOF data can be explicitly elicited and exploited by the solutions along the characteristic curves. As the special cases of the general solutions, we obtain Fourier rebinning and consistency equations (FORCEs), and thus we obtain a complete scheme to convert among different types of PET planograms: 3D TOF, 3D non-TOF, 2D TOF and 2D non-TOF planograms. The FORCEs can be used as Fourier-based rebinning algorithms for TOF-PET data reduction, inverse rebinnings for designing fast projectors, or consistency conditions for estimating missing data. As a byproduct, we show the two consistency equations are necessary and sufficient for 3D TOF planograms. Finally, we give
Optical modulators with 2D layered materials
NASA Astrophysics Data System (ADS)
Sun, Zhipei; Martinez, Amos; Wang, Feng
2016-04-01
Light modulation is an essential operation in photonics and optoelectronics. With existing and emerging technologies increasingly demanding compact, efficient, fast and broadband optical modulators, high-performance light modulation solutions are becoming indispensable. The recent realization that 2D layered materials could modulate light with superior performance has prompted intense research and significant advances, paving the way for realistic applications. In this Review, we cover the state of the art of optical modulators based on 2D materials, including graphene, transition metal dichalcogenides and black phosphorus. We discuss recent advances employing hybrid structures, such as 2D heterostructures, plasmonic structures, and silicon and fibre integrated structures. We also take a look at the future perspectives and discuss the potential of yet relatively unexplored mechanisms, such as magneto-optic and acousto-optic modulation.
Microholography by Numerical Wavefront Reconstruction in the Fourier Transform Geometry.
NASA Astrophysics Data System (ADS)
Haddad, Waleed Sami
New techniques for pulsed holographic microimaging with x-rays and visible light based on the Fourier transform geometry are presented. A charge-coupled device (CCD) is used to record the holographic information. Images are reconstructed numerically by fast Fourier transform (FFT) and by a novel "stigmatic" aberration-free technique. It is shown by numerical simulation that reconstruction by FFT inherently suffers from coma. An innovative design for a Fourier transform x -ray holographic microscope incorporating a metallic microsphere as the reference wave producing element is presented. The design is optimized for use with an x-ray laser source operating in the range lambda~eq43 A; nickel is found to be the best material for the reference sphere in this range. A visible light analog of the x-ray microscope was built and tested yielding an image of Ascaris larva. Two unique visible light Fourier transform holographic microscopes are also presented. The first incorporates a glycerol microdrop reference lens. This microscope was operated at N.A. ~ 0.25 with an argon ion laser at lambda = 514.5 nm. The transverse point spread function of the system was measured to be 1.40 mum, near the theoretical diffraction limit of 1.29 mu m. The second is a design based on a Mach-Zehnder interferometer. This Fourier Mach Zehnder holographic microscope (FMZHM) offers great flexibility and many advantages. The FMZHM was operated at N.A. _sp{~ }{>} 0.25 with lambda = 514.5 nm and yielded many 2-D and 3-D images of biological samples. A study of potential biological applications of visible and x-ray holographic microscopy is also presented.
Mainali, Dipak; Seelenbinder, John
2016-05-01
Quick and presumptive identification of seized drug samples without destroying evidence is necessary for law enforcement officials to control the trafficking and abuse of drugs. This work reports an automated screening method to detect the presence of cocaine in seized samples using portable Fourier transform infrared (FT-IR) spectrometers. The method is based on the identification of well-defined characteristic vibrational frequencies related to the functional group of the cocaine molecule and is fully automated through the use of an expert system. Traditionally, analysts look for key functional group bands in the infrared spectra and characterization of the molecules present is dependent on user interpretation. This implies the need for user expertise, especially in samples that likely are mixtures. As such, this approach is biased and also not suitable for non-experts. The method proposed in this work uses the well-established "center of gravity" peak picking mathematical algorithm and combines it with the conditional reporting feature in MicroLab software to provide an automated method that can be successfully employed by users with varied experience levels. The method reports the confidence level of cocaine present only when a certain number of cocaine related peaks are identified by the automated method. Unlike library search and chemometric methods that are dependent on the library database or the training set samples used to build the calibration model, the proposed method is relatively independent of adulterants and diluents present in the seized mixture. This automated method in combination with a portable FT-IR spectrometer provides law enforcement officials, criminal investigators, or forensic experts a quick field-based prescreening capability for the presence of cocaine in seized drug samples. PMID:27006022
Fast Steerable Principal Component Analysis
Zhao, Zhizhen; Shkolnisky, Yoel; Singer, Amit
2016-01-01
Cryo-electron microscopy nowadays often requires the analysis of hundreds of thousands of 2-D images as large as a few hundred pixels in each direction. Here, we introduce an algorithm that efficiently and accurately performs principal component analysis (PCA) for a large set of 2-D images, and, for each image, the set of its uniform rotations in the plane and their reflections. For a dataset consisting of n images of size L × L pixels, the computational complexity of our algorithm is O(nL3 + L4), while existing algorithms take O(nL4). The new algorithm computes the expansion coefficients of the images in a Fourier–Bessel basis efficiently using the nonuniform fast Fourier transform. We compare the accuracy and efficiency of the new algorithm with traditional PCA and existing algorithms for steerable PCA. PMID:27570801
2D signature for detection and identification of drugs
NASA Astrophysics Data System (ADS)
Trofimov, Vyacheslav A.; Varentsova, Svetlana A.; Shen, Jingling; Zhang, Cunlin; Zhou, Qingli; Shi, Yulei
2011-06-01
The method of spectral dynamics analysis (SDA-method) is used for obtaining the2D THz signature of drugs. This signature is used for the detection and identification of drugs with similar Fourier spectra by transmitted THz signal. We discuss the efficiency of SDA method for the identification problem of pure methamphetamine (MA), methylenedioxyamphetamine (MDA), 3, 4-methylenedioxymethamphetamine (MDMA) and Ketamine.
The fractional Fourier transform and applications
NASA Technical Reports Server (NTRS)
Bailey, David H.; Swarztrauber, Paul N.
1991-01-01
This paper describes the 'fractional Fourier transform', which admits computation by an algorithm that has complexity proportional to the fast Fourier transform algorithm. Whereas the discrete Fourier transform (DFT) is based on integral roots of unity e exp -2(pi)i/n, the fractional Fourier transform is based on fractional roots of unity e exp -2(pi)i(alpha), where alpha is arbitrary. The fractional Fourier transform and the corresponding fast algorithm are useful for such applications as computing DFTs of sequences with prime lengths, computing DFTs of sparse sequences, analyzing sequences with noninteger periodicities, performing high-resolution trigonometric interpolation, detecting lines in noisy images, and detecting signals with linearly drifting frequencies. In many cases, the resulting algorithms are faster by arbitrarily large factors than conventional techniques.
NASA Astrophysics Data System (ADS)
Wang, Jin; Ma, Jianyong; Zhou, Changhe
2014-11-01
A 3×3 high divergent 2D-grating with period of 3.842μm at wavelength of 850nm under normal incidence is designed and fabricated in this paper. This high divergent 2D-grating is designed by the vector theory. The Rigorous Coupled Wave Analysis (RCWA) in association with the simulated annealing (SA) is adopted to calculate and optimize this 2D-grating.The properties of this grating are also investigated by the RCWA. The diffraction angles are more than 10 degrees in the whole wavelength band, which are bigger than the traditional 2D-grating. In addition, the small period of grating increases the difficulties of fabrication. So we fabricate the 2D-gratings by direct laser writing (DLW) instead of traditional manufacturing method. Then the method of ICP etching is used to obtain the high divergent 2D-grating.
Reducing acquisition times in multidimensional NMR with a time-optimized Fourier encoding algorithm.
Zhang, Zhiyong; Smith, Pieter E S; Frydman, Lucio
2014-11-21
Speeding up the acquisition of multidimensional nuclear magnetic resonance (NMR) spectra is an important topic in contemporary NMR, with central roles in high-throughput investigations and analyses of marginally stable samples. A variety of fast NMR techniques have been developed, including methods based on non-uniform sampling and Hadamard encoding, that overcome the long sampling times inherent to schemes based on fast-Fourier-transform (FFT) methods. Here, we explore the potential of an alternative fast acquisition method that leverages a priori knowledge, to tailor polychromatic pulses and customized time delays for an efficient Fourier encoding of the indirect domain of an NMR experiment. By porting the encoding of the indirect-domain to the excitation process, this strategy avoids potential artifacts associated with non-uniform sampling schemes and uses a minimum number of scans equal to the number of resonances present in the indirect dimension. An added convenience is afforded by the fact that a usual 2D FFT can be used to process the generated data. Acquisitions of 2D heteronuclear correlation NMR spectra on quinine and on the anti-inflammatory drug isobutyl propionic phenolic acid illustrate the new method's performance. This method can be readily automated to deal with complex samples such as those occurring in metabolomics, in in-cell as well as in in vivo NMR applications, where speed and temporal stability are often primary concerns. PMID:25416883
Reducing acquisition times in multidimensional NMR with a time-optimized Fourier encoding algorithm
Zhang, Zhiyong; Frydman, Lucio
2014-01-01
Speeding up the acquisition of multidimensional nuclear magnetic resonance (NMR) spectra is an important topic in contemporary NMR, with central roles in high-throughput investigations and analyses of marginally stable samples. A variety of fast NMR techniques have been developed, including methods based on non-uniform sampling and Hadamard encoding, that overcome the long sampling times inherent to schemes based on fast-Fourier-transform (FFT) methods. Here, we explore the potential of an alternative fast acquisition method that leverages a priori knowledge, to tailor polychromatic pulses and customized time delays for an efficient Fourier encoding of the indirect domain of an NMR experiment. By porting the encoding of the indirect-domain to the excitation process, this strategy avoids potential artifacts associated with non-uniform sampling schemes and uses a minimum number of scans equal to the number of resonances present in the indirect dimension. An added convenience is afforded by the fact that a usual 2D FFT can be used to process the generated data. Acquisitions of 2D heteronuclear correlation NMR spectra on quinine and on the anti-inflammatory drug isobutyl propionic phenolic acid illustrate the new method's performance. This method can be readily automated to deal with complex samples such as those occurring in metabolomics, in in-cell as well as in in vivo NMR applications, where speed and temporal stability are often primary concerns. PMID:25416883
Reducing acquisition times in multidimensional NMR with a time-optimized Fourier encoding algorithm
Zhang, Zhiyong; Smith, Pieter E. S.; Frydman, Lucio
2014-11-21
Speeding up the acquisition of multidimensional nuclear magnetic resonance (NMR) spectra is an important topic in contemporary NMR, with central roles in high-throughput investigations and analyses of marginally stable samples. A variety of fast NMR techniques have been developed, including methods based on non-uniform sampling and Hadamard encoding, that overcome the long sampling times inherent to schemes based on fast-Fourier-transform (FFT) methods. Here, we explore the potential of an alternative fast acquisition method that leverages a priori knowledge, to tailor polychromatic pulses and customized time delays for an efficient Fourier encoding of the indirect domain of an NMR experiment. By porting the encoding of the indirect-domain to the excitation process, this strategy avoids potential artifacts associated with non-uniform sampling schemes and uses a minimum number of scans equal to the number of resonances present in the indirect dimension. An added convenience is afforded by the fact that a usual 2D FFT can be used to process the generated data. Acquisitions of 2D heteronuclear correlation NMR spectra on quinine and on the anti-inflammatory drug isobutyl propionic phenolic acid illustrate the new method's performance. This method can be readily automated to deal with complex samples such as those occurring in metabolomics, in in-cell as well as in in vivo NMR applications, where speed and temporal stability are often primary concerns.
Tomosynthesis imaging with 2D scanning trajectories
NASA Astrophysics Data System (ADS)
Khare, Kedar; Claus, Bernhard E. H.; Eberhard, Jeffrey W.
2011-03-01
Tomosynthesis imaging in chest radiography provides volumetric information with the potential for improved diagnostic value when compared to the standard AP or LAT projections. In this paper we explore the image quality benefits of 2D scanning trajectories when coupled with advanced image reconstruction approaches. It is intuitively clear that 2D trajectories provide projection data that is more complete in terms of Radon space filling, when compared with conventional tomosynthesis using a linearly scanned source. Incorporating this additional information for obtaining improved image quality is, however, not a straightforward problem. The typical tomosynthesis reconstruction algorithms are based on direct inversion methods e.g. Filtered Backprojection (FBP) or iterative algorithms that are variants of the Algebraic Reconstruction Technique (ART). The FBP approach is fast and provides high frequency details in the image but at the same time introduces streaking artifacts degrading the image quality. The iterative methods can reduce the image artifacts by using image priors but suffer from a slow convergence rate, thereby producing images lacking high frequency details. In this paper we propose using a fast converging optimal gradient iterative scheme that has advantages of both the FBP and iterative methods in that it produces images with high frequency details while reducing the image artifacts. We show that using favorable 2D scanning trajectories along with the proposed reconstruction method has the advantage of providing improved depth information for structures such as the spine and potentially producing images with more isotropic resolution.
Reconstruction in Fourier space
NASA Astrophysics Data System (ADS)
Burden, A.; Percival, W. J.; Howlett, C.
2015-10-01
We present a fast iterative fast Fourier transform (FFT) based reconstruction algorithm that allows for non-parallel redshift-space distortions (RSDs). We test our algorithm on both N-body dark matter simulations and mock distributions of galaxies designed to replicate galaxy survey conditions. We compare solenoidal and irrotational components of the redshift distortion and show that an approximation of this distortion leads to a better estimate of the real-space potential (and therefore faster convergence) than ignoring the RSD when estimating the displacement field. Our iterative reconstruction scheme converges in two iterations for the mock samples corresponding to Baryon Oscillation Spectroscopic Survey CMASS Data Release 11 when we start with an approximation of the RSD. The scheme takes six iterations when the initial estimate, measured from the redshift-space overdensity, has no RSD correction. Slower convergence would be expected for surveys covering a larger angle on the sky. We show that this FFT based method provides a better estimate of the real-space displacement field than a configuration space method that uses finite difference routines to compute the potential for the same grid resolution. Finally, we show that a lognormal transform of the overdensity, used as a proxy for the linear overdensity, is beneficial in estimating the full displacement field from a dense sample of tracers. However, the lognormal transform of the overdensity does not perform well when estimating the displacements from sparser simulations with a more realistic galaxy density.
Bingi, Jayachandra; Murukeshan, Vadakke Matham
2015-01-01
Laser speckle pattern is a granular structure formed due to random coherent wavelet interference and generally considered as noise in optical systems including photolithography. Contrary to this, in this paper, we use the speckle pattern to generate predictable and controlled Gaussian random structures and quasi-random structures photo-lithographically. The random structures made using this proposed speckle lithography technique are quantified based on speckle statistics, radial distribution function (RDF) and fast Fourier transform (FFT). The control over the speckle size, density and speckle clustering facilitates the successful fabrication of black silicon with different surface structures. The controllability and tunability of randomness makes this technique a robust method for fabricating predictable 2D Gaussian random structures and black silicon structures. These structures can enhance the light trapping significantly in solar cells and hence enable improved energy harvesting. Further, this technique can enable efficient fabrication of disordered photonic structures and random media based devices. PMID:26679513
Bingi, Jayachandra; Murukeshan, Vadakke Matham
2015-01-01
Laser speckle pattern is a granular structure formed due to random coherent wavelet interference and generally considered as noise in optical systems including photolithography. Contrary to this, in this paper, we use the speckle pattern to generate predictable and controlled Gaussian random structures and quasi-random structures photo-lithographically. The random structures made using this proposed speckle lithography technique are quantified based on speckle statistics, radial distribution function (RDF) and fast Fourier transform (FFT). The control over the speckle size, density and speckle clustering facilitates the successful fabrication of black silicon with different surface structures. The controllability and tunability of randomness makes this technique a robust method for fabricating predictable 2D Gaussian random structures and black silicon structures. These structures can enhance the light trapping significantly in solar cells and hence enable improved energy harvesting. Further, this technique can enable efficient fabrication of disordered photonic structures and random media based devices. PMID:26679513
Ultrafast Fourier-transform parallel processor
Greenberg, W.L.
1980-04-01
A new, flexible, parallel-processing architecture is developed for a high-speed, high-precision Fourier transform processor. The processor is intended for use in 2-D signal processing including spatial filtering, matched filtering and image reconstruction from projections.
Aperture scanning Fourier ptychographic microscopy
Ou, Xiaoze; Chung, Jaebum; Horstmeyer, Roarke; Yang, Changhuei
2016-01-01
Fourier ptychographic microscopy (FPM) is implemented through aperture scanning by an LCOS spatial light modulator at the back focal plane of the objective lens. This FPM configuration enables the capturing of the complex scattered field for a 3D sample both in the transmissive mode and the reflective mode. We further show that by combining with the compressive sensing theory, the reconstructed 2D complex scattered field can be used to recover the 3D sample scattering density. This implementation expands the scope of application for FPM and can be beneficial for areas such as tissue imaging and wafer inspection. PMID:27570705
NASA Technical Reports Server (NTRS)
Alexandrov, Mikhail D.; Cairns, Brian; Mishchenko, Michael I.
2012-01-01
We present a novel technique for remote sensing of cloud droplet size distributions. Polarized reflectances in the scattering angle range between 135deg and 165deg exhibit a sharply defined rainbow structure, the shape of which is determined mostly by single scattering properties of cloud particles, and therefore, can be modeled using the Mie theory. Fitting the observed rainbow with such a model (computed for a parameterized family of particle size distributions) has been used for cloud droplet size retrievals. We discovered that the relationship between the rainbow structures and the corresponding particle size distributions is deeper than it had been commonly understood. In fact, the Mie theory-derived polarized reflectance as a function of reduced scattering angle (in the rainbow angular range) and the (monodisperse) particle radius appears to be a proxy to a kernel of an integral transform (similar to the sine Fourier transform on the positive semi-axis). This approach, called the rainbow Fourier transform (RFT), allows us to accurately retrieve the shape of the droplet size distribution by the application of the corresponding inverse transform to the observed polarized rainbow. While the basis functions of the proxy-transform are not exactly orthogonal in the finite angular range, this procedure needs to be complemented by a simple regression technique, which removes the retrieval artifacts. This non-parametric approach does not require any a priori knowledge of the droplet size distribution functional shape and is computationally fast (no look-up tables, no fitting, computations are the same as for the forward modeling).
Energy Science and Technology Software Center (ESTSC)
2004-08-01
AnisWave2D is a 2D finite-difference code for a simulating seismic wave propagation in fully anisotropic materials. The code is implemented to run in parallel over multiple processors and is fully portable. A mesh refinement algorithm has been utilized to allow the grid-spacing to be tailored to the velocity model, avoiding the over-sampling of high-velocity materials that usually occurs in fixed-grid schemes.
WFR-2D: an analytical model for PWAS-generated 2D ultrasonic guided wave propagation
NASA Astrophysics Data System (ADS)
Shen, Yanfeng; Giurgiutiu, Victor
2014-03-01
This paper presents WaveFormRevealer 2-D (WFR-2D), an analytical predictive tool for the simulation of 2-D ultrasonic guided wave propagation and interaction with damage. The design of structural health monitoring (SHM) systems and self-aware smart structures requires the exploration of a wide range of parameters to achieve best detection and quantification of certain types of damage. Such need for parameter exploration on sensor dimension, location, guided wave characteristics (mode type, frequency, wavelength, etc.) can be best satisfied with analytical models which are fast and efficient. The analytical model was constructed based on the exact 2-D Lamb wave solution using Bessel and Hankel functions. Damage effects were inserted in the model by considering the damage as a secondary wave source with complex-valued directivity scattering coefficients containing both amplitude and phase information from wave-damage interaction. The analytical procedure was coded with MATLAB, and a predictive simulation tool called WaveFormRevealer 2-D was developed. The wave-damage interaction coefficients (WDICs) were extracted from harmonic analysis of local finite element model (FEM) with artificial non-reflective boundaries (NRB). The WFR-2D analytical simulation results were compared and verified with full scale multiphysics finite element models and experiments with scanning laser vibrometer. First, Lamb wave propagation in a pristine aluminum plate was simulated with WFR-2D, compared with finite element results, and verified by experiments. Then, an inhomogeneity was machined into the plate to represent damage. Analytical modeling was carried out, and verified by finite element simulation and experiments. This paper finishes with conclusions and suggestions for future work.
Prozy Visco-plasticity, Fast Fourier Transform
Energy Science and Technology Software Center (ESTSC)
2012-10-25
ProxyVPFFT computes the micromechanical reaponse of polycriystalline samples directly from images of their microstructures, usually captured using SEM. [2] ProxyVPFFT is proxy application, or a somplified version of existing software by Ricardo Lebensohn (VPFFT).
Efficient 2D MRI relaxometry using compressed sensing
NASA Astrophysics Data System (ADS)
Bai, Ruiliang; Cloninger, Alexander; Czaja, Wojciech; Basser, Peter J.
2015-06-01
Potential applications of 2D relaxation spectrum NMR and MRI to characterize complex water dynamics (e.g., compartmental exchange) in biology and other disciplines have increased in recent years. However, the large amount of data and long MR acquisition times required for conventional 2D MR relaxometry limits its applicability for in vivo preclinical and clinical MRI. We present a new MR pipeline for 2D relaxometry that incorporates compressed sensing (CS) as a means to vastly reduce the amount of 2D relaxation data needed for material and tissue characterization without compromising data quality. Unlike the conventional CS reconstruction in the Fourier space (k-space), the proposed CS algorithm is directly applied onto the Laplace space (the joint 2D relaxation data) without compressing k-space to reduce the amount of data required for 2D relaxation spectra. This framework is validated using synthetic data, with NMR data acquired in a well-characterized urea/water phantom, and on fixed porcine spinal cord tissue. The quality of the CS-reconstructed spectra was comparable to that of the conventional 2D relaxation spectra, as assessed using global correlation, local contrast between peaks, peak amplitude and relaxation parameters, etc. This result brings this important type of contrast closer to being realized in preclinical, clinical, and other applications.
2D materials and van der Waals heterostructures.
Novoselov, K S; Mishchenko, A; Carvalho, A; Castro Neto, A H
2016-07-29
The physics of two-dimensional (2D) materials and heterostructures based on such crystals has been developing extremely fast. With these new materials, truly 2D physics has begun to appear (for instance, the absence of long-range order, 2D excitons, commensurate-incommensurate transition, etc.). Novel heterostructure devices--such as tunneling transistors, resonant tunneling diodes, and light-emitting diodes--are also starting to emerge. Composed from individual 2D crystals, such devices use the properties of those materials to create functionalities that are not accessible in other heterostructures. Here we review the properties of novel 2D crystals and examine how their properties are used in new heterostructure devices. PMID:27471306
Fast fully 3-D image reconstruction in PET using planograms.
Brasse, D; Kinahan, P E; Clackdoyle, R; Defrise, M; Comtat, C; Townsend, D W
2004-04-01
We present a method of performing fast and accurate three-dimensional (3-D) backprojection using only Fourier transform operations for line-integral data acquired by planar detector arrays in positron emission tomography. This approach is a 3-D extension of the two-dimensional (2-D) linogram technique of Edholm. By using a special choice of parameters to index a line of response (LOR) for a pair of planar detectors, rather than the conventional parameters used to index a LOR for a circular tomograph, all the LORs passing through a point in the field of view (FOV) lie on a 2-D plane in the four-dimensional (4-D) data space. Thus, backprojection of all the LORs passing through a point in the FOV corresponds to integration of a 2-D plane through the 4-D "planogram." The key step is that the integration along a set of parallel 2-D planes through the planogram, that is, backprojection of a plane of points, can be replaced by a 2-D section through the origin of the 4-D Fourier transform of the data. Backprojection can be performed as a sequence of Fourier transform operations, for faster implementation. In addition, we derive the central-section theorem for planogram format data, and also derive a reconstruction filter for both backprojection-filtering and filtered-backprojection reconstruction algorithms. With software-based Fourier transform calculations we provide preliminary comparisons of planogram backprojection to standard 3-D backprojection and demonstrate a reduction in computation time by a factor of approximately 15. PMID:15084067
NASA Astrophysics Data System (ADS)
Mayor, Louise
2016-05-01
Graphene might be the most famous example, but there are other 2D materials and compounds too. Louise Mayor explains how these atomically thin sheets can be layered together to create flexible “van der Waals heterostructures”, which could lead to a range of novel applications.
Fourier Series Operating Package
NASA Technical Reports Server (NTRS)
Charnow, Milton L.
1961-01-01
This report presents a computer program for multiplying, adding, differentiating, integrating, "barring" and scalarly multiplying "literal" Fourier series as such, and for extracting the coefficients of specified terms.
Wideband aperture array using RF channelizers and massively parallel digital 2D IIR filterbank
NASA Astrophysics Data System (ADS)
Sengupta, Arindam; Madanayake, Arjuna; Gómez-García, Roberto; Engeberg, Erik D.
2014-05-01
Wideband receive-mode beamforming applications in wireless location, electronically-scanned antennas for radar, RF sensing, microwave imaging and wireless communications require digital aperture arrays that offer a relatively constant far-field beam over several octaves of bandwidth. Several beamforming schemes including the well-known true time-delay and the phased array beamformers have been realized using either finite impulse response (FIR) or fast Fourier transform (FFT) digital filter-sum based techniques. These beamforming algorithms offer the desired selectivity at the cost of a high computational complexity and frequency-dependant far-field array patterns. A novel approach to receiver beamforming is the use of massively parallel 2-D infinite impulse response (IIR) fan filterbanks for the synthesis of relatively frequency independent RF beams at an order of magnitude lower multiplier complexity compared to FFT or FIR filter based conventional algorithms. The 2-D IIR filterbanks demand fast digital processing that can support several octaves of RF bandwidth, fast analog-to-digital converters (ADCs) for RF-to-bits type direct conversion of wideband antenna element signals. Fast digital implementation platforms that can realize high-precision recursive filter structures necessary for real-time beamforming, at RF radio bandwidths, are also desired. We propose a novel technique that combines a passive RF channelizer, multichannel ADC technology, and single-phase massively parallel 2-D IIR digital fan filterbanks, realized at low complexity using FPGA and/or ASIC technology. There exists native support for a larger bandwidth than the maximum clock frequency of the digital implementation technology. We also strive to achieve More-than-Moore throughput by processing a wideband RF signal having content with N-fold (B = N Fclk/2) bandwidth compared to the maximum clock frequency Fclk Hz of the digital VLSI platform under consideration. Such increase in bandwidth is
Sehgal, Akansha Ashvani; Pelupessy, Philippe; Rolando, Christian; Bodenhausen, Geoffrey
2016-04-01
Two-dimensional (2D) Fourier transform ion cyclotron resonance (FT-ICR) offers an approach to mass spectrometry (MS) that pursuits similar objectives as MS/MS experiments. While the latter must focus on one ion species at a time, 2D FT ICR can examine all possible correlations due to ion fragmentation in a single experiment: correlations between precursors, charged and neutral fragments. We revisited the original 2D FT-ICR experiment that has hitherto fallen short of stimulating significant analytical applications, probably because it is technically demanding. These shortcomings can now be overcome by improved FT-ICR instrumentation and computer hard- and software. We seek to achieve a better understanding of the intricacies of the behavior of ions during a basic two-dimensional ICR sequence comprising three simple monochromatic pulses. Through simulations based on Lorentzian equations, we have mapped the ion trajectories for different pulse durations and phases. PMID:26974979
Energy Science and Technology Software Center (ESTSC)
2001-01-31
This software reduces the data from two-dimensional kSA MOS program, k-Space Associates, Ann Arbor, MI. Initial MOS data is recorded without headers in 38 columns, with one row of data per acquisition per lase beam tracked. The final MOSS 2d data file is reduced, graphed, and saved in a tab-delimited column format with headers that can be plotted in any graphing software.
NASA Astrophysics Data System (ADS)
Sterken, C.
2003-03-01
This paper gives a short account of some key elements in the life of Jean Baptiste Joseph Fourier (1768-1830), specifically his relation to Napoleon Bonaparte. The mathematical approach to Fourier series and the original scepticism by French mathematicians are briefly illustrated.
Fourier Series Optimization Opportunity
ERIC Educational Resources Information Center
Winkel, Brian
2008-01-01
This note discusses the introduction of Fourier series as an immediate application of optimization of a function of more than one variable. Specifically, it is shown how the study of Fourier series can be motivated to enrich a multivariable calculus class. This is done through discovery learning and use of technology wherein students build the…
Hyun, Eugin; Jin, Young-Seok; Lee, Jong-Hun
2016-01-01
For an automotive pedestrian detection radar system, fast-ramp based 2D range-Doppler Frequency Modulated Continuous Wave (FMCW) radar is effective for distinguishing between moving targets and unwanted clutter. However, when a weak moving target such as a pedestrian exists together with strong clutter, the pedestrian may be masked by the side-lobe of the clutter even though they are notably separated in the Doppler dimension. To prevent this problem, one popular solution is the use of a windowing scheme with a weighting function. However, this method leads to a spread spectrum, so the pedestrian with weak signal power and slow Doppler may also be masked by the main-lobe of clutter. With a fast-ramp based FMCW radar, if the target is moving, the complex spectrum of the range- Fast Fourier Transform (FFT) is changed with a constant phase difference over ramps. In contrast, the clutter exhibits constant phase irrespective of the ramps. Based on this fact, in this paper we propose a pedestrian detection for highly cluttered environments using a coherent phase difference method. By detecting the coherent phase difference from the complex spectrum of the range-FFT, we first extract the range profile of the moving pedestrians. Then, through the Doppler FFT, we obtain the 2D range-Doppler map for only the pedestrian. To test the proposed detection scheme, we have developed a real-time data logging system with a 24 GHz FMCW transceiver. In laboratory tests, we verified that the signal processing results from the proposed method were much better than those expected from the conventional 2D FFT-based detection method. PMID:26805835
Hyun, Eugin; Jin, Young-Seok; Lee, Jong-Hun
2016-01-01
For an automotive pedestrian detection radar system, fast-ramp based 2D range-Doppler Frequency Modulated Continuous Wave (FMCW) radar is effective for distinguishing between moving targets and unwanted clutter. However, when a weak moving target such as a pedestrian exists together with strong clutter, the pedestrian may be masked by the side-lobe of the clutter even though they are notably separated in the Doppler dimension. To prevent this problem, one popular solution is the use of a windowing scheme with a weighting function. However, this method leads to a spread spectrum, so the pedestrian with weak signal power and slow Doppler may also be masked by the main-lobe of clutter. With a fast-ramp based FMCW radar, if the target is moving, the complex spectrum of the range- Fast Fourier Transform (FFT) is changed with a constant phase difference over ramps. In contrast, the clutter exhibits constant phase irrespective of the ramps. Based on this fact, in this paper we propose a pedestrian detection for highly cluttered environments using a coherent phase difference method. By detecting the coherent phase difference from the complex spectrum of the range-FFT, we first extract the range profile of the moving pedestrians. Then, through the Doppler FFT, we obtain the 2D range-Doppler map for only the pedestrian. To test the proposed detection scheme, we have developed a real-time data logging system with a 24 GHz FMCW transceiver. In laboratory tests, we verified that the signal processing results from the proposed method were much better than those expected from the conventional 2D FFT-based detection method. PMID:26805835
Fehler, M.C.; Huang, L.-J.
1998-12-10
During the past few years, there has been interest in developing migration and forward modeling approaches that are both fast and reliable particularly in regions that have rapid spatial variations in structure. The authors have been investigating a suite of modeling and migration methods that are implemented in the wavenumber-space domains and operate on data in the frequency domain. The best known example of these methods is the split-step Fourier method (SSF). Two of the methods that the authors have developed are the extended local Born Fourier (ELBF) approach and the extended local Rytov Fourier (ELRF) approach. Both methods are based on solutions of the scalar (constant density) wave equation, are computationally fast and can reliably model effects of both deterministic and random structures. The authors have investigated their reliability for migrating both 2D synthetic data and real 2D field data. The authors have found that the methods give images that are better than those that can be obtained using other methods like the SSF and Kirchhoff migration approaches. More recently, the authors have developed an approach for solving the acoustic (variable density) wave equation and have begun to investigate its applicability for modeling one-way wave propagation. The methods will be introduced and their ability to model seismic wave propagation and migrate seismic data will be investigated. The authors will also investigate their capability to model forward wave propagation through random media and to image zones of small scale heterogeneity such as those associated with zones of high permeability.
Estimating 2-D vector velocities using multidimensional spectrum analysis.
Oddershede, Niels; Løvstakken, Lasse; Torp, Hans; Jensen, Jørgen Arendt
2008-08-01
Wilson (1991) presented an ultrasonic wideband estimator for axial blood flow velocity estimation through the use of the 2-D Fourier transform. It was shown how a single velocity component was concentrated along a line in the 2-D Fourier space, where the slope was given by the axial velocity. Later, it was shown that this approach could also be used for finding the lateral velocity component by also including a lateral sampling. A single velocity component would then be concentrated along a plane in the 3-D Fourier space, tilted according to the 2 velocity components. This paper presents 2 new velocity estimators for finding both the axial and lateral velocity components. The estimators essentially search for the plane in the 3- D Fourier space, where the integrated power spectrum is largest. The first uses the 3-D Fourier transform to find the power spectrum, while the second uses a minimum variance approach. Based on this plane, the axial and lateral velocity components are estimated. Several phantom measurements, for flow-to-depth angles of 60, 75, and 90 degrees, were performed. Multiple parallel lines were beamformed simultaneously, and 2 different receive apodization schemes were tried. The 2 estimators were then applied to the data. The axial velocity component was estimated with an average standard deviation below 2.8% of the peak velocity, while the average standard deviation of the lateral velocity estimates was between 2.0% and 16.4%. The 2 estimators were also tested on in vivo data from a transverse scan of the common carotid artery, showing the potential of the vector velocity estimation method under in vivo conditions. PMID:18986918
Nanoimprint lithography: 2D or not 2D? A review
NASA Astrophysics Data System (ADS)
Schift, Helmut
2015-11-01
Nanoimprint lithography (NIL) is more than a planar high-end technology for the patterning of wafer-like substrates. It is essentially a 3D process, because it replicates various stamp topographies by 3D displacement of material and takes advantage of the bending of stamps while the mold cavities are filled. But at the same time, it keeps all assets of a 2D technique being able to pattern thin masking layers like in photon- and electron-based traditional lithography. This review reports about 20 years of development of replication techniques at Paul Scherrer Institut, with a focus on 3D aspects of molding, which enable NIL to stay 2D, but at the same time enable 3D applications which are "more than Moore." As an example, the manufacturing of a demonstrator for backlighting applications based on thermally activated selective topography equilibration will be presented. This technique allows generating almost arbitrary sloped, convex and concave profiles in the same polymer film with dimensions in micro- and nanometer scale.
Fourier analysis and signal processing by use of the Moebius inversion formula
NASA Technical Reports Server (NTRS)
Reed, Irving S.; Yu, Xiaoli; Shih, Ming-Tang; Tufts, Donald W.; Truong, T. K.
1990-01-01
A novel Fourier technique for digital signal processing is developed. This approach to Fourier analysis is based on the number-theoretic method of the Moebius inversion of series. The Fourier transform method developed is shown also to yield the convolution of two signals. A computer simulation shows that this method for finding Fourier coefficients is quite suitable for digital signal processing. It competes with the classical FFT (fast Fourier transform) approach in terms of accuracy, complexity, and speed.
Imaging Fourier transform spectrometer
Bennett, C.L.
1993-09-13
This invention is comprised of an imaging Fourier transform spectrometer having a Fourier transform infrared spectrometer providing a series of images to a focal plane array camera. The focal plane array camera is clocked to a multiple of zero crossing occurrences as caused by a moving mirror of the Fourier transform infrared spectrometer and as detected by a laser detector such that the frame capture rate of the focal plane array camera corresponds to a multiple of the zero crossing rate of the Fourier transform infrared spectrometer. The images are transmitted to a computer for processing such that representations of the images as viewed in the light of an arbitrary spectral ``fingerprint`` pattern can be displayed on a monitor or otherwise stored and manipulated by the computer.
Fourier Transform Mass Spectrometry
Scigelova, Michaela; Hornshaw, Martin; Giannakopulos, Anastassios; Makarov, Alexander
2011-01-01
This article provides an introduction to Fourier transform-based mass spectrometry. The key performance characteristics of Fourier transform-based mass spectrometry, mass accuracy and resolution, are presented in the view of how they impact the interpretation of measurements in proteomic applications. The theory and principles of operation of two types of mass analyzer, Fourier transform ion cyclotron resonance and Orbitrap, are described. Major benefits as well as limitations of Fourier transform-based mass spectrometry technology are discussed in the context of practical sample analysis, and illustrated with examples included as figures in this text and in the accompanying slide set. Comparisons highlighting the performance differences between the two mass analyzers are made where deemed useful in assisting the user with choosing the most appropriate technology for an application. Recent developments of these high-performing mass spectrometers are mentioned to provide a future outlook. PMID:21742802
Generalized fiber Fourier optics.
Cincotti, Gabriella
2011-06-15
A twofold generalization of the optical schemes that perform the discrete Fourier transform (DFT) is given: new passive planar architectures are presented where the 2 × 2 3 dB couplers are replaced by M × M hybrids, reducing the number of required connections and phase shifters. Furthermore, the planar implementation of the discrete fractional Fourier transform (DFrFT) is also described, with a waveguide grating router (WGR) configuration and a properly modified slab coupler. PMID:21686007
ERIC Educational Resources Information Center
Debnath, Lokenath
2012-01-01
This article deals with a brief biographical sketch of Joseph Fourier, his first celebrated work on analytical theory of heat, his first great discovery of Fourier series and Fourier transforms. Included is a historical development of Fourier series and Fourier transforms with their properties, importance and applications. Special emphasis is made…
O'Brien, Travis A.; Kashinath, Karthik
2015-05-22
This software implements the fast, self-consistent probability density estimation described by O'Brien et al. (2014, doi: ). It uses a non-uniform fast Fourier transform technique to reduce the computational cost of an objective and self-consistent kernel density estimation method.
Fourier plane imaging microscopy
Dominguez, Daniel Peralta, Luis Grave de; Alharbi, Nouf; Alhusain, Mdhaoui; Bernussi, Ayrton A.
2014-09-14
We show how the image of an unresolved photonic crystal can be reconstructed using a single Fourier plane (FP) image obtained with a second camera that was added to a traditional compound microscope. We discuss how Fourier plane imaging microscopy is an application of a remarkable property of the obtained FP images: they contain more information about the photonic crystals than the images recorded by the camera commonly placed at the real plane of the microscope. We argue that the experimental results support the hypothesis that surface waves, contributing to enhanced resolution abilities, were optically excited in the studied photonic crystals.
2D imaging of functional structures in perfused pig heart
NASA Astrophysics Data System (ADS)
Kessler, Manfred D.; Cristea, Paul D.; Hiller, Michael; Trinks, Tobias
2002-06-01
In 2000 by 2D-imaging we were able for the first time to visualize in subcellular space functional structures of myocardium. For these experiments we used hemoglobin-free perfused pig hearts in our lab. Step by step we learned to understand the meaning of subcellular structures. Principally, the experiment revealed that in subcellular space very fast changes of light scattering can occur. Furthermore, coefficients of different parameters were determined on the basis of multicomponent system theory.
Introduction to Fourier Optics
ERIC Educational Resources Information Center
Huggins, Elisha
2007-01-01
Much like a physical prism, which displays the frequency components of a light wave, Fourier analysis can be thought of as a mathematical prism that can tell us what harmonics or frequency components are contained in a recording of a sound wave. We wrote the MacScope II program so that the user could not only see a plot of the harmonic amplitudes…
Fourier Transform Mass Spectrometry.
ERIC Educational Resources Information Center
Gross, Michael L.; Rempel, Don L.
1984-01-01
Discusses the nature of Fourier transform mass spectrometry and its unique combination of high mass resolution, high upper mass limit, and multichannel advantage. Examines its operation, capabilities and limitations, applications (ion storage, ion manipulation, ion chemistry), and future applications and developments. (JN)
Fourier Transform Spectrometer
NASA Technical Reports Server (NTRS)
1998-01-01
Understanding the global atmospheric changes is difficult with today's current technology. However, with high resolution and nearly continuous observations from a satellite, it's possible to transform our understanding of the atmosphere. To enable the next generation of atmospheric science, a new class of orbiting atmospheric sensors is being developed. The foundation of this advanced concept is the Fourier Transform Spectrometer, or FTS.
NASA Astrophysics Data System (ADS)
Sandirasegaram, Nicholas; English, Ryan
2005-05-01
The performance of several combinations of feature extraction and target classification algorithms is analyzed for Synthetic Aperture Radar (SAR) imagery using the standard Moving and Stationary Target Acquisition and Recognition (MSTAR) evaluation method. For feature extraction, 2D Fast Fourier Transform (FFT) is used to extract Fourier coefficients (frequency information) while 2D wavelet decomposition is used to extract wavelet coefficients (time-frequency information), from which subsets of characteristic in-class "invariant" coefficients are developed. Confusion matrices and Receiver Operating Characteristic (ROC) curves are used to evaluate and compare combinations of these characteristic coefficients with several classification methods, including Lp metric distances, a Multi Layer Perceptron (MLP) Neural Network (NN) and AND Corporation's Holographic Neural Technology (HNeT) classifier. The evaluation method examines the trade-off between correct detection rate and false alarm rate for each combination of feature-classifier systems. It also measures correct classification, misclassification and rejection rates for a 90% detection rate. Our analysis demonstrates the importance of feature and classifier selection in accurately classifying new target images.
NKG2D ligands as therapeutic targets
Spear, Paul; Wu, Ming-Ru; Sentman, Marie-Louise; Sentman, Charles L.
2013-01-01
The Natural Killer Group 2D (NKG2D) receptor plays an important role in protecting the host from infections and cancer. By recognizing ligands induced on infected or tumor cells, NKG2D modulates lymphocyte activation and promotes immunity to eliminate ligand-expressing cells. Because these ligands are not widely expressed on healthy adult tissue, NKG2D ligands may present a useful target for immunotherapeutic approaches in cancer. Novel therapies targeting NKG2D ligands for the treatment of cancer have shown preclinical success and are poised to enter into clinical trials. In this review, the NKG2D receptor and its ligands are discussed in the context of cancer, infection, and autoimmunity. In addition, therapies targeting NKG2D ligands in cancer are also reviewed. PMID:23833565
Independent task Fourier filters
NASA Astrophysics Data System (ADS)
Caulfield, H. John
2001-11-01
Since the early 1960s, a major part of optical computing systems has been Fourier pattern recognition, which takes advantage of high speed filter changes to enable powerful nonlinear discrimination in `real time.' Because filter has a task quite independent of the tasks of the other filters, they can be applied and evaluated in parallel or, in a simple approach I describe, in sequence very rapidly. Thus I use the name ITFF (independent task Fourier filter). These filters can also break very complex discrimination tasks into easily handled parts, so the wonderful space invariance properties of Fourier filtering need not be sacrificed to achieve high discrimination and good generalizability even for ultracomplex discrimination problems. The training procedure proceeds sequentially, as the task for a given filter is defined a posteriori by declaring it to be the discrimination of particular members of set A from all members of set B with sufficient margin. That is, we set the threshold to achieve the desired margin and note the A members discriminated by that threshold. Discriminating those A members from all members of B becomes the task of that filter. Those A members are then removed from the set A, so no other filter will be asked to perform that already accomplished task.
Perspectives for spintronics in 2D materials
NASA Astrophysics Data System (ADS)
Han, Wei
2016-03-01
The past decade has been especially creative for spintronics since the (re)discovery of various two dimensional (2D) materials. Due to the unusual physical characteristics, 2D materials have provided new platforms to probe the spin interaction with other degrees of freedom for electrons, as well as to be used for novel spintronics applications. This review briefly presents the most important recent and ongoing research for spintronics in 2D materials.
A discrete Fourier transform for virtual memory machines
NASA Technical Reports Server (NTRS)
Galant, David C.
1992-01-01
An algebraic theory of the Discrete Fourier Transform is developed in great detail. Examination of the details of the theory leads to a computationally efficient fast Fourier transform for the use on computers with virtual memory. Such an algorithm is of great use on modern desktop machines. A FORTRAN coded version of the algorithm is given for the case when the sequence of numbers to be transformed is a power of two.
NASA Astrophysics Data System (ADS)
Wu, Yan-Wen; Sun, Su-Qin; Zhou, Qun; Tao, Jia-Xun; Noda, Isao
2008-06-01
As a traditional Chinese medicine (TCM), 'Red Flower Oil' preparation is widely used as a household remedy in China and Southeast Asia. Usually, the preparation is a mixture of several plant essential oils with different volatile features, such as wintergreen oil, turpentine oil and clove oil. The proportions of these plant essential oils in 'Red Flower Oil' vary from different manufacturers. Thus, it is important to develop a simple and rapid evaluation method for quality assurance of the preparations. Fourier transform infrared (FT-IR) was applied and two-dimensional correlation infrared spectroscopy (2D IR) based on the volatile characteristic of samples was used to enhance the resolution of FT-IR spectra. 2D IR technique could, not only easily provide the composition and their volatile sequences in 'Red flower Oil' preparations, but also rapidly discriminate the subtle differences in products from different manufacturers. Therefore, FT-IR combined with volatility-dependent 2D IR correlation analysis provides a very fast and effective method for the quality control of essential oil mixtures in TCM.
Fourier technique for studying ammonoid sutures
Gildner, R.F.; Ackerly, S.C.
1985-01-01
Suture patterns have long been recognized as being of primary importance in the study of ammonoids. The authors have developed a technique to use Fourier analysis to study these structures by using a simple transformation: x-y data of a digitized suture are transformed to angle of slope versus position along the suture's length. A Fast Fourier Transform applied to the data produces a power spectrum (amplitude versus wave number) providing a precise and accurate measure of suture shape. The authors have applied this technique to the analysis of ontogenetic change in suture morphology. In goniatitic, ceratitic and preadult ammonitic patterns most of the change is exhibited in the amplitudes of the lowest ten wave numbers. Their Fourier coefficients clearly show trends not readily apparent by visual inspection. The more complex ammonitic patterns are reflected in increased amplitudes of higher wave numbers (a broader peak of the power spectrum) and their analysis is necessarily more complex. The Fourier approach presents the opportunity to quantitatively measure and describe the tempo and mode of evolution in the Ammonoidea. Potential applications of the new technique, as well as limitations, are discussed with special attention to investigations of ammonoid ontogeny and phylogeny.
Astronomical Fourier spectrometer.
Connes, P; Michel, G
1975-09-01
A high resolution near ir Fourier spectrometer with the same general design as previously described laboratory instruments has been built for astronomical observations at a coudé focus. Present spectral range is 0.8-3.5 microm with PbS and Ge detectors and maximum path difference 1 m. The servo system can accommodate various recording modes: stepping or continuous scan, path difference modulation, sky chopping. A real time computer is incorporated into the system, which has been set up at the Hale 500-cm telescope on Mount Palomar. Samples of the results are given. PMID:20154966
NASA Technical Reports Server (NTRS)
Feldkhun, Daniel (Inventor); Wagner, Kelvin H. (Inventor)
2013-01-01
Methods and systems are disclosed of sensing an object. A first radiation is spatially modulated to generate a structured second radiation. The object is illuminated with the structured second radiation such that the object produces a third radiation in response. Apart from any spatially dependent delay, a time variation of the third radiation is spatially independent. With a single-element detector, a portion of the third radiation is detected from locations on the object simultaneously. At least one characteristic of a sinusoidal spatial Fourier-transform component of the object is estimated from a time-varying signal from the detected portion of the third radiation.
Single-shot parallel full range complex Fourier-domain optical coherence tomography
NASA Astrophysics Data System (ADS)
Huang, Bingjie; Bu, Peng; Nan, Nan; Wang, Xiangzhao
2011-01-01
We present a method of parallel full range complex Fourier-domain optical coherence tomography (FDOCT) that is capable of acquiring an artifacts-free two-dimensional (2-D) cross-sectional image, i.e. a full range B-scan tomogram, by a single shot of 2-D CCD camera. This method is based on a spatial carrier technique, in which the spatial carrier-frequency is instantaneously introduced into the 2-D spectral interferogram registered in parallel FDOCT by using a grating-generated reference beam. The spatial-carrier-contained 2-D spectral interferogram is processed through Fourier transformation to obtain a complex 2-D spectral interferogram. From the 2-D complex spectral interferomgram, a full range B-scan tomogram is reconstructed. The principle of our method is confirmed by imaging an onion sample.
Annotated Bibliography of EDGE2D Use
J.D. Strachan and G. Corrigan
2005-06-24
This annotated bibliography is intended to help EDGE2D users, and particularly new users, find existing published literature that has used EDGE2D. Our idea is that a person can find existing studies which may relate to his intended use, as well as gain ideas about other possible applications by scanning the attached tables.
Staring 2-D hadamard transform spectral imager
Gentry, Stephen M.; Wehlburg, Christine M.; Wehlburg, Joseph C.; Smith, Mark W.; Smith, Jody L.
2006-02-07
A staring imaging system inputs a 2D spatial image containing multi-frequency spectral information. This image is encoded in one dimension of the image with a cyclic Hadamarid S-matrix. The resulting image is detecting with a spatial 2D detector; and a computer applies a Hadamard transform to recover the encoded image.
Hackel, L.A.; Hermann, M.R.; Dane, C.B.; Tiszauer, D.H.
1995-12-12
A solid state laser is frequency tripled to 0.3 {micro}m. A small portion of the laser is split off and generates a Stokes seed in a low power oscillator. The low power output passes through a mask with the appropriate hole pattern. Meanwhile, the bulk of the laser output is focused into a larger stimulated Brillouin scattering (SBS) amplifier. The low power beam is directed through the same cell in the opposite direction. The majority of the amplification takes place at the focus which is the fourier transform plane of the mask image. The small holes occupy large area at the focus and thus are preferentially amplified. The amplified output is now imaged onto the multichip module where the holes are drilled. Because of the fourier plane amplifier, only about 1/10th the power of a competitive system is needed. This concept allows less expensive masks to be used in the process and requires much less laser power. 1 fig.
Hackel, Lloyd A.; Hermann, Mark R.; Dane, C. Brent; Tiszauer, Detlev H.
1995-01-01
A solid state laser is frequency tripled to 0.3 .mu.m. A small portion of the laser is split off and generates a Stokes seed in a low power oscillator. The low power output passes through a mask with the appropriate hole pattern. Meanwhile, the bulk of the laser output is focused into a larger stimulated Brillouin scattering (SBS) amplifier. The low power beam is directed through the same cell in the opposite direction. The majority of the amplification takes place at the focus which is the fourier transform plane of the mask image. The small holes occupy large area at the focus and thus are preferentially amplified. The amplified output is now imaged onto the multichip module where the holes are drilled. Because of the fourier plane amplifier, only .about.1/10th the power of a competitive system is needed. This concept allows less expensive masks to be used in the process and requires much less laser power.
2-D soft x-ray arrays in the EAST.
Chen, Kaiyun; Xu, Liqing; Hu, Liqun; Duan, Yanmin; Li, Xueqin; Yuan, Yi; Mao, Songtao; Sheng, Xiuli; Zhao, Jinlong
2016-06-01
A high spatial and temporal resolution soft x-ray (SXR) imaging diagnostic has been installed in EAST for the study of magnetohydrodynamics activities and core high-Z impurity transport. Up to 122 lines of sight view the poloidal plasma from three directions (two up-down symmetrical horizontal arrays and one vertical array), which renders the diagnostic able to provide detailed tomographic reconstructions under various conditions. Fourier-Bessel method based on flux coordinates was employed for 2-D SXR tomographic reconstruction. Examples of several events measured by SXR diagnostic in EAST are shown, namely the crash patterns of sawtooth, periodical burst of edge localized modes, and the transport of high-Z intrinsic impurities. PMID:27370451
2-D soft x-ray arrays in the EAST
NASA Astrophysics Data System (ADS)
Chen, Kaiyun; Xu, Liqing; Hu, Liqun; Duan, Yanmin; Li, Xueqin; Yuan, Yi; Mao, Songtao; Sheng, Xiuli; Zhao, Jinlong
2016-06-01
A high spatial and temporal resolution soft x-ray (SXR) imaging diagnostic has been installed in EAST for the study of magnetohydrodynamics activities and core high-Z impurity transport. Up to 122 lines of sight view the poloidal plasma from three directions (two up-down symmetrical horizontal arrays and one vertical array), which renders the diagnostic able to provide detailed tomographic reconstructions under various conditions. Fourier-Bessel method based on flux coordinates was employed for 2-D SXR tomographic reconstruction. Examples of several events measured by SXR diagnostic in EAST are shown, namely the crash patterns of sawtooth, periodical burst of edge localized modes, and the transport of high-Z intrinsic impurities.
A study of geometric phase topology using Fourier transform method
NASA Astrophysics Data System (ADS)
Samlan, C. T.; Naik, Dinesh N.; Viswanathan, Nirmal K.
2016-07-01
Topological aspect of the geometric phase (GP) due to pure polarization projection is studied using the 2D Fourier transform (2D-FT) method. Projection of orthogonal polarization state results in a phase singularity in the 2D parameter space of ellipticity and orientation of polarization ellipse. Projection of its surrounding states results in an accumulation of GP in different amount that form a spiral structure. A half wave plate–quarter wave plate combination is used to generate different polarization states which are projected using a polarizer. The accumulated phase for each orientation of the wave plate is extracted from 2D-FT of the interferogram, obtained by interfering it with a reference beam in a Mach–Zehnder like interferometer.
Light field morphing using 2D features.
Wang, Lifeng; Lin, Stephen; Lee, Seungyong; Guo, Baining; Shum, Heung-Yeung
2005-01-01
We present a 2D feature-based technique for morphing 3D objects represented by light fields. Existing light field morphing methods require the user to specify corresponding 3D feature elements to guide morph computation. Since slight errors in 3D specification can lead to significant morphing artifacts, we propose a scheme based on 2D feature elements that is less sensitive to imprecise marking of features. First, 2D features are specified by the user in a number of key views in the source and target light fields. Then the two light fields are warped view by view as guided by the corresponding 2D features. Finally, the two warped light fields are blended together to yield the desired light field morph. Two key issues in light field morphing are feature specification and warping of light field rays. For feature specification, we introduce a user interface for delineating 2D features in key views of a light field, which are automatically interpolated to other views. For ray warping, we describe a 2D technique that accounts for visibility changes and present a comparison to the ideal morphing of light fields. Light field morphing based on 2D features makes it simple to incorporate previous image morphing techniques such as nonuniform blending, as well as to morph between an image and a light field. PMID:15631126
2D materials for nanophotonic devices
NASA Astrophysics Data System (ADS)
Xu, Renjing; Yang, Jiong; Zhang, Shuang; Pei, Jiajie; Lu, Yuerui
2015-12-01
Two-dimensional (2D) materials have become very important building blocks for electronic, photonic, and phononic devices. The 2D material family has four key members, including the metallic graphene, transition metal dichalcogenide (TMD) layered semiconductors, semiconducting black phosphorous, and the insulating h-BN. Owing to the strong quantum confinements and defect-free surfaces, these atomically thin layers have offered us perfect platforms to investigate the interactions among photons, electrons and phonons. The unique interactions in these 2D materials are very important for both scientific research and application engineering. In this talk, I would like to briefly summarize and highlight the key findings, opportunities and challenges in this field. Next, I will introduce/highlight our recent achievements. We demonstrated atomically thin micro-lens and gratings using 2D MoS2, which is the thinnest optical component around the world. These devices are based on our discovery that the elastic light-matter interactions in highindex 2D materials is very strong. Also, I would like to introduce a new two-dimensional material phosphorene. Phosphorene has strongly anisotropic optical response, which creates 1D excitons in a 2D system. The strong confinement in phosphorene also enables the ultra-high trion (charged exciton) binding energies, which have been successfully measured in our experiments. Finally, I will briefly talk about the potential applications of 2D materials in energy harvesting.