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.
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
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)
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.
3-D Printed Slit Nozzles for Fourier Transform Microwave Spectroscopy
NASA Astrophysics Data System (ADS)
Dewberry, Chris; Mackenzie, Becca; Green, Susan; Leopold, Ken
2015-06-01
3-D printing is a new technology whose applications are only beginning to be explored. In this report, we describe the application of 3-D printing to the facile design and construction of supersonic nozzles. The efficacy of a variety of designs is assessed by examining rotational spectra OCS and Ar-OCS using a Fourier transform microwave spectrometer with tandem cavity and chirped-pulse capabilities. This work focuses primarily on the use of slit nozzles but other designs have been tested as well. New nozzles can be created for 0.50 or less each, and the ease and low cost should facilitate the optimization of nozzle performance (e.g., jet temperature or cluster size distribution) for the needs of any particular experiment.
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
FastScript3D - A Companion to Java 3D
NASA Technical Reports Server (NTRS)
Koenig, Patti
2005-01-01
FastScript3D is a computer program, written in the Java 3D(TM) programming language, that establishes an alternative language that helps users who lack expertise in Java 3D to use Java 3D for constructing three-dimensional (3D)-appearing graphics. The FastScript3D language provides a set of simple, intuitive, one-line text-string commands for creating, controlling, and animating 3D models. The first word in a string is the name of a command; the rest of the string contains the data arguments for the command. The commands can also be used as an aid to learning Java 3D. Developers can extend the language by adding custom text-string commands. The commands can define new 3D objects or load representations of 3D objects from files in formats compatible with such other software systems as X3D. The text strings can be easily integrated into other languages. FastScript3D facilitates communication between scripting languages [which enable programming of hyper-text markup language (HTML) documents to interact with users] and Java 3D. The FastScript3D language can be extended and customized on both the scripting side and the Java 3D side.
3D palmprint data fast acquisition and recognition
NASA Astrophysics Data System (ADS)
Wang, Xiaoxu; Huang, Shujun; Gao, Nan; Zhang, Zonghua
2014-11-01
This paper presents a fast 3D (Three-Dimension) palmprint capturing system and develops an efficient 3D palmprint feature extraction and recognition method. In order to fast acquire accurate 3D shape and texture of palmprint, a DLP projector triggers a CCD camera to realize synchronization. By generating and projecting green fringe pattern images onto the measured palm surface, 3D palmprint data are calculated from the fringe pattern images. The periodic feature vector can be derived from the calculated 3D palmprint data, so undistorted 3D biometrics is obtained. Using the obtained 3D palmprint data, feature matching test have been carried out by Gabor filter, competition rules and the mean curvature. Experimental results on capturing 3D palmprint show that the proposed acquisition method can fast get 3D shape information of palmprint. Some initial experiments on recognition show the proposed method is efficient by using 3D palmprint data.
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.
3D fast wavelet network model-assisted 3D face recognition
NASA Astrophysics Data System (ADS)
Said, Salwa; Jemai, Olfa; Zaied, Mourad; Ben Amar, Chokri
2015-12-01
In last years, the emergence of 3D shape in face recognition is due to its robustness to pose and illumination changes. These attractive benefits are not all the challenges to achieve satisfactory recognition rate. Other challenges such as facial expressions and computing time of matching algorithms remain to be explored. In this context, we propose our 3D face recognition approach using 3D wavelet networks. Our approach contains two stages: learning stage and recognition stage. For the training we propose a novel algorithm based on 3D fast wavelet transform. From 3D coordinates of the face (x,y,z), we proceed to voxelization to get a 3D volume which will be decomposed by 3D fast wavelet transform and modeled after that with a wavelet network, then their associated weights are considered as vector features to represent each training face . For the recognition stage, an unknown identity face is projected on all the training WN to obtain a new vector features after every projection. A similarity score is computed between the old and the obtained vector features. To show the efficiency of our approach, experimental results were performed on all the FRGC v.2 benchmark.
Fast 3D shape measurements using laser speckle projection
NASA Astrophysics Data System (ADS)
Schaffer, Martin; Grosse, Marcus; Harendt, Bastian; Kowarschik, Richard
2011-05-01
3D measurement setups based on structured light projection are widely used for many industrial applications. Due to intense research in the past the accuracy is comparably high in connection with relatively low cost of the equipment. But facing higher acquisition rates in industries especially for chain assembling lines there are still hurdles to take when accelerating 3D measurements and at the same time retaining accuracies. We developed a projection technique that uses laser speckles to enable fast 3D measurements with statistically structured light patterns. In combination with a temporal correlation technique dense and accurate 3D reconstructions at nearly video rate can be achieved.
Increased Speed: 3D Silicon Sensors. Fast Current Amplifiers
Parker, Sherwood; Kok, Angela; Kenney, Christopher; Jarron, Pierre; Hasi, Jasmine; Despeisse, Matthieu; Da Via, Cinzia; Anelli, Giovanni; /CERN
2012-05-07
The authors describe techniques to make fast, sub-nanosecond time resolution solid-state detector systems using sensors with 3D electrodes, current amplifiers, constant-fraction comparators or fast wave-form recorders, and some of the next steps to reach still faster results.
Fast and precise 3D fluorophore localization by gradient fitting
NASA Astrophysics Data System (ADS)
Ma, Hongqiang; Xu, Jianquan; Jin, Jingyi; Gao, Ying; Lan, Li; Liu, Yang
2016-02-01
Astigmatism imaging is widely used to encode the 3D position of fluorophore in single-particle tracking and super-resolution localization microscopy. Here, we present a fast and precise localization algorithm based on gradient fitting to decode the 3D subpixel position of the fluorophore. This algorithm determines the center of the emitter by finding the position with the best-fit gradient direction distribution to the measured point spread function (PSF), and can retrieve the 3D subpixel position of the emitter in a single iteration. Through numerical simulation and experiments with mammalian cells, we demonstrate that our algorithm yields comparable localization precision to the traditional iterative Gaussian function fitting (GF) based method, while exhibits over two orders-of-magnitude faster execution speed. Our algorithm is a promising online reconstruction method for 3D super-resolution microscopy.
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 Mode Decision for 3D-HEVC Depth Intracoding
Li, Nana; Wu, Qinggang
2014-01-01
The emerging international standard of high efficiency video coding based 3D video coding (3D-HEVC) is a successor to multiview video coding (MVC). In 3D-HEVC depth intracoding, depth modeling mode (DMM) and high efficiency video coding (HEVC) intraprediction mode are both employed to select the best coding mode for each coding unit (CU). This technique achieves the highest possible coding efficiency, but it results in extremely large encoding time which obstructs the 3D-HEVC from practical application. In this paper, a fast mode decision algorithm based on the correlation between texture video and depth map is proposed to reduce 3D-HEVC depth intracoding computational complexity. Since the texture video and its associated depth map represent the same scene, there is a high correlation among the prediction mode from texture video and depth map. Therefore, we can skip some specific depth intraprediction modes rarely used in related texture CU. Experimental results show that the proposed algorithm can significantly reduce computational complexity of 3D-HEVC depth intracoding while maintaining coding efficiency. PMID:24963512
a Fast Method for Measuring the Similarity Between 3d Model and 3d Point Cloud
NASA Astrophysics Data System (ADS)
Zhang, Zongliang; Li, Jonathan; Li, Xin; Lin, Yangbin; Zhang, Shanxin; Wang, Cheng
2016-06-01
This paper proposes a fast method for measuring the partial Similarity between 3D Model and 3D point Cloud (SimMC). It is crucial to measure SimMC for many point cloud-related applications such as 3D object retrieval and inverse procedural modelling. In our proposed method, the surface area of model and the Distance from Model to point Cloud (DistMC) are exploited as measurements to calculate SimMC. Here, DistMC is defined as the weighted distance of the distances between points sampled from model and point cloud. Similarly, Distance from point Cloud to Model (DistCM) is defined as the average distance of the distances between points in point cloud and model. In order to reduce huge computational burdens brought by calculation of DistCM in some traditional methods, we define SimMC as the ratio of weighted surface area of model to DistMC. Compared to those traditional SimMC measuring methods that are only able to measure global similarity, our method is capable of measuring partial similarity by employing distance-weighted strategy. Moreover, our method is able to be faster than other partial similarity assessment methods. We demonstrate the superiority of our method both on synthetic data and laser scanning data.
Fast 3D stereo flood simulations in urban areas
NASA Astrophysics Data System (ADS)
Hoes, O.; de Haan, G.; Stelling, G.; van Leeuwen, E.; van Dam, A.; Pleumeekers, O.; Schuurmans, W.
2012-04-01
Flood propagation models are essential to study floods as it is problematic to collect data during actual floods. Moreover, models are needed to explore the consequences of additional scenarios above the actual flood itself. The results of these model studies are generally graphs with water levels over time for certain locations or maps with the flood extent in an area for different return periods. A main difficulty with these depictions of flood information is that they do not reflect the seriousness of flood impacts well in terms of life-like human experience. Typically, one needs a (near) flood before measures are implemented. Apparently, a graph or map is not the proper material to convince politicians and policy makers, even if they live in the threatened area. The recent introduction of commercially available 3D stereo projectors and high resolution elevation data make it possible to build life-like visualizations of simulations. In our research we explored using 3D stereo, the recently collected elevation data of the Netherlands (20 laser points per m2!) in combination with aerial photographs, and a new fast 2D flood propagation calculation scheme. This scheme (under construction) is able to simulate floods using such high amounts of data points. The model simulates flood propagation on an irregular grid; at locations with large elevation differences (e.g. in urban areas) and fast flowing water, smaller cells are used compared to flat surfaces where the water is not or hardly flowing. The result of our combination is a very detailed flood simulation model that can be used to simulate floods within a fraction of the current calculation time. The opportunities of models and their results increase enormously with fast calculations and visualizations combined. For example, the model allows on the spot exploration of measures during a flood, with the 3D visualization ensuring that flood impacts become clear for decision makers. We will show the preliminary
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.
Holographic display of real existing objects from their 3D Fourier spectrum
NASA Astrophysics Data System (ADS)
Yatagai, Toyohiko; Sando, Yusuke
2005-02-01
A method for synthesizing computer-generated holograms of real-existing objects is described. A series of projection images are recorded both vertically and horizontally with an incoherent light source and a color CCD camera. According to the principle of computer tomography(CT), the 3-D Fourier spectrum is calculated from several projection images of objects and the Fresnel computer-generated hologram(CGH) is synthesized using a part of the 3-D Fourier spectrum. This method has following advantages. At first, no-blur reconstructed images in any direction are obtained owing to two-dimensionally scanning in recording. Secondarily, since not interference fringes but simple projection images of objects are recorded, a coherent light source is not necessary for recording. The use of a color CCD in recording enables us to record and reconstruct colorful objects. Finally, we demonstrate color reconstruction of objects both numerically and optically.
Undersampling k-space using fast progressive 3D trajectories.
Spiniak, Juan; Guesalaga, Andres; Mir, Roberto; Guarini, Marcelo; Irarrazaval, Pablo
2005-10-01
In 3D MRI, sampling k-space with traditional trajectories can be excessively time-consuming. Fast imaging trajectories are used in an attempt to efficiently cover the k-space and reduce the scan time without significantly affecting the image quality. In many applications, further reductions in scan time can be achieved via undersampling of the k-space; however, no clearly optimal method exists. In most 3D trajectories the k-space is divided into regions that are sampled with shots that share a common geometry (e.g., spirals). A different approach is to design trajectories that gradually but uniformly cover the k-space. In the current work, successive shots progressively add sampled regions to the 3D frequency space. By cutting the sequence short, a natural undersampled method is obtained. This can be particularly efficient because in these types of trajectories the contribution of new information by later shots is less significant. In this work the performance of progressive trajectories for different degrees of undersampling is assessed with trajectories based on missile guidance (MG) ideas. The results show that the approach can be efficient in terms of reducing the scan time, and performs better than the stack of spirals (SOS) technique, particularly under nonideal conditions. PMID:16142719
Trillo, Cristina; Doval, Angel F; Mendoza-Santoyo, Fernando; Pérez-López, Carlos; de la Torre-Ibarra, Manuel; Deán, J Luis
2009-09-28
The combination of a high-speed TV holography system and a 3D Fourier-transform data processing is proposed for the analysis of multimode vibrations in plates. The out-of-plane displacement of the object under generic vibrational excitation is resolved in time by the fast acquisition rate of a high-speed camera, and recorded in a sequence of interferograms with spatial carrier. A full-field temporal history of the multimode vibration is thus obtained. The optical phase of the interferograms is extracted and subtracted from the phase of a reference state to yield a sequence of optical phase-change maps. Each map represents the change undergone by the object between any given state and the reference state. The sequence of maps is a 3D array of data (two spatial dimensions plus time) that is processed with a 3D Fourier-transform algorithm. The individual vibration modes are separated in the 3D frequency space due to their different vibration frequencies and, to a lesser extent, to the different spatial frequencies of the mode shapes. The contribution of each individual mode (or indeed the superposition of several modes) to the dynamic behaviour of the object can then be separated by means of a bandpass filter (or filters). The final output is a sequence of complex-valued maps that contain the full-field temporal history of the selected mode (or modes) in terms of its mechanical amplitude and phase. The proof-of-principle of the technique is demonstrated with a rectangular, fully clamped, thin metal plate vibrating simultaneously in several of its natural resonant frequencies under white-noise excitation. PMID:19907591
3-D lookup: Fast protein structure database searches
Holm. L.; Sander, C.
1995-12-31
There are far fewer classes of three-dimensional protein folds than sequence families but the problem of detecting three-dimensional similarities is NP-complete. We present a novel heuristic for identifying 3-D similarities between a query structure and the database of known protein structures. Many methods for structure alignment use a bottom-up approach, identifying first local matches and then solving a combinatorial problem in building up larger clusters of matching substructures. Here the top-down approach is to start with the global comparison and select a rough superimposition using a fast 3-D lookup of secondary structure motifs. The superimposition is then extended to an alignment of C{sup {alpha}} atoms by an iterative dynamic programming step. An all-against-all comparison of 385-representative proteins (150,000 pair comparisons) took 1 day of computer time on a single R8000 processor. In other words, one query structure is scanned against the database in a matter of minutes. The method is rated at 90% reliability at capturing statistically significant similarities. It is useful as a rapid preprocessor to a comprehensive protein structure database search system.
Imaging open-path Fourier transform infrared spectrometer for 3D cloud profiling
NASA Astrophysics Data System (ADS)
Rentz Dupuis, Julia; Mansur, David J.; Vaillancourt, Robert; Carlson, David; Evans, Thomas; Schundler, Elizabeth; Todd, Lori; Mottus, Kathleen
2009-05-01
OPTRA is developing an imaging open-path Fourier transform infrared (I-OP-FTIR) spectrometer for 3D profiling of chemical and biological agent simulant plumes released into test ranges and chambers. An array of I-OP-FTIR instruments positioned around the perimeter of the test site, in concert with advanced spectroscopic algorithms, enables real time tomographic reconstruction of the plume. The approach is intended as a referee measurement for test ranges and chambers. This Small Business Technology Transfer (STTR) effort combines the instrumentation and spectroscopic capabilities of OPTRA, Inc. with the computed tomographic expertise of the University of North Carolina, Chapel Hill.
3D homogeneity study in PMMA layers using a Fourier domain OCT system
NASA Astrophysics Data System (ADS)
Briones-R., Manuel de J.; Torre-Ibarra, Manuel H. De La; Tavera, Cesar G.; Luna H., Juan M.; Mendoza-Santoyo, Fernando
2016-11-01
Micro-metallic particles embedded in polymers are now widely used in several industrial applications in order to modify the mechanical properties of the bulk. A uniform distribution of these particles inside the polymers is highly desired for instance, when a biological backscattering is simulated or a bio-framework is designed. A 3D Fourier domain optical coherence tomography system to detect the polymer's internal homogeneity is proposed. This optical system has a 2D camera sensor array that records a fringe pattern used to reconstruct with a single shot the tomographic image of the sample. The system gathers the full 3D tomographic and optical phase information during a controlled deformation by means of a motion linear stage. This stage avoids the use of expensive tilting stages, which in addition are commonly controlled by piezo drivers. As proof of principle, a series of different deformations were proposed to detect the uniform or non-uniform internal deposition of copper micro particles. The results are presented as images coming from the 3D tomographic micro reconstruction of the samples, and the 3D optical phase information that identifies the in-homogeneity regions within the Poly methyl methacrylate (PMMA) volume.
Research of Fast 3D Imaging Based on Multiple Mode
NASA Astrophysics Data System (ADS)
Chen, Shibing; Yan, Huimin; Ni, Xuxiang; Zhang, Xiuda; Wang, Yu
2016-02-01
Three-dimensional (3D) imaging has received increasingly extensive attention and has been widely used currently. Lots of efforts have been put on three-dimensional imaging method and system study, in order to meet fast and high accurate requirement. In this article, we realize a fast and high quality stereo matching algorithm on field programmable gate array (FPGA) using the combination of time-of-flight (TOF) camera and binocular camera. Images captured from the two cameras own a same spatial resolution, letting us use the depth maps taken by the TOF camera to figure initial disparity. Under the constraint of the depth map as the stereo pairs when comes to stereo matching, expected disparity of each pixel is limited within a narrow search range. In the meanwhile, using field programmable gate array (FPGA, altera cyclone IV series) concurrent computing we can configure multi core image matching system, thus doing stereo matching on embedded system. The simulation results demonstrate that it can speed up the process of stereo matching and increase matching reliability and stability, realize embedded calculation, expand application range.
Reconstruction for 3D PET Based on Total Variation Constrained Direct Fourier Method
Yu, Haiqing; Chen, Zhi; Zhang, Heye; Loong Wong, Kelvin Kian; Chen, Yunmei; Liu, Huafeng
2015-01-01
This paper presents a total variation (TV) regularized reconstruction algorithm for 3D positron emission tomography (PET). The proposed method first employs the Fourier rebinning algorithm (FORE), rebinning the 3D data into a stack of ordinary 2D data sets as sinogram data. Then, the resulted 2D sinogram are ready to be reconstructed by conventional 2D reconstruction algorithms. Given the locally piece-wise constant nature of PET images, we introduce the total variation (TV) based reconstruction schemes. More specifically, we formulate the 2D PET reconstruction problem as an optimization problem, whose objective function consists of TV norm of the reconstructed image and the data fidelity term measuring the consistency between the reconstructed image and sinogram. To solve the resulting minimization problem, we apply an efficient methods called the Bregman operator splitting algorithm with variable step size (BOSVS). Experiments based on Monte Carlo simulated data and real data are conducted as validations. The experiment results show that the proposed method produces higher accuracy than conventional direct Fourier (DF) (bias in BOSVS is 70% of ones in DF, variance of BOSVS is 80% of ones in DF). PMID:26398232
Anisotropic magnification distortion of the 3D galaxy correlation. II. Fourier and redshift space
Hui Lam; Gaztanaga, Enrique; LoVerde, Marilena
2008-03-15
In paper I of this series we discuss how magnification bias distorts the 3D correlation function by enhancing the observed correlation in the line-of-sight (LOS) orientation, especially on large scales. This lensing anisotropy is distinctive, making it possible to separately measure the galaxy-galaxy, galaxy-magnification and magnification-magnification correlations. Here we extend the discussion to the power spectrum and also to redshift space. In real space, pairs oriented close to the LOS direction are not protected against nonlinearity even if the pair separation is large; this is because nonlinear fluctuations can enter through gravitational lensing at a small transverse separation (or i.e. impact parameter). The situation in Fourier space is different: by focusing on a small wave number k, as is usually done, linearity is guaranteed because both the LOS and transverse wave numbers must be small. This is why magnification distortion of the galaxy correlation appears less severe in Fourier space. Nonetheless, the effect is non-negligible, especially for the transverse Fourier modes, and should be taken into account in interpreting precision measurements of the galaxy power spectrum, for instance those that focus on the baryon oscillations. The lensing induced anisotropy of the power spectrum has a shape that is distinct from the more well-known redshift space anisotropies due to peculiar motions and the Alcock-Paczynski effect. The lensing anisotropy is highly localized in Fourier space while redshift space distortions are more spread out. This means that one could separate the magnification bias component in real observations, implying that potentially it is possible to perform a gravitational lensing measurement without measuring galaxy shapes.
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.
Imaging open-path Fourier transform infrared spectrometer for 3D cloud profiling
NASA Astrophysics Data System (ADS)
Rentz Dupuis, Julia; Mansur, David J.; Vaillancourt, Robert; Carlson, David; Evans, Thomas; Schundler, Elizabeth; Todd, Lori; Mottus, Kathleen
2010-04-01
OPTRA has developed an imaging open-path Fourier transform infrared (I-OP-FTIR) spectrometer for 3D profiling of chemical and biological agent simulant plumes released into test ranges and chambers. An array of I-OP-FTIR instruments positioned around the perimeter of the test site, in concert with advanced spectroscopic algorithms, enables real time tomographic reconstruction of the plume. The approach is intended as a referee measurement for test ranges and chambers. This Small Business Technology Transfer (STTR) effort combines the instrumentation and spectroscopic capabilities of OPTRA, Inc. with the computed tomographic expertise of the University of North Carolina, Chapel Hill. In this paper, we summarize the design and build and detail system characterization and test of a prototype I-OP-FTIR instrument. System characterization includes radiometric performance and spectral resolution. Results from a series of tomographic reconstructions of sulfur hexafluoride plumes in a laboratory setting are also presented.
Imaging open-path Fourier transform infrared spectrometer for 3D cloud profiling
NASA Astrophysics Data System (ADS)
Rentz Dupuis, Julia; Mansur, David J.; Engel, James R.; Vaillancourt, Robert; Todd, Lori; Mottus, Kathleen
2008-04-01
OPTRA and University of North Carolina are developing an imaging open-path Fourier transform infrared (I-OP-FTIR) spectrometer for 3D profiling of chemical and biological agent simulant plumes released into test ranges and chambers. An array of I-OP-FTIR instruments positioned around the perimeter of the test site, in concert with advanced spectroscopic algorithms, enables real time tomographic reconstruction of the plume. The approach will be considered as a candidate referee measurement for test ranges and chambers. This Small Business Technology Transfer (STTR) effort combines the instrumentation and spectroscopic capabilities of OPTRA, Inc. with the computed tomographic expertise of the University of North Carolina, Chapel Hill. In this paper, we summarize progress to date and overall system performance projections based on the instrument, spectroscopy, and tomographic reconstruction accuracy. We then present a preliminary optical design of the I-OP-FTIR.
Fast phase-added stereogram algorithm for generation of photorealistic 3D content.
Kang, Hoonjong; Stoykova, Elena; Yoshikawa, Hiroshi
2016-01-20
A new phase-added stereogram algorithm for accelerated computation of holograms from a point cloud model is proposed. The algorithm relies on the hologram segmentation, sampling of directional information, and usage of the fast Fourier transform with a finer grid in the spatial frequency domain than is provided by the segment size. The algorithm gives improved quality of reconstruction due to new phase compensation introduced in the segment fringe patterns. The result is finer beam steering leading to high peak intensity and a large peak signal-to-noise ratio in reconstruction. The feasibility of the algorithm is checked by the generation of 3D contents for a color wavefront printer. PMID:26835945
3D RISM theory with fast reciprocal-space electrostatics
Heil, Jochen; Kast, Stefan M.
2015-03-21
The calculation of electrostatic solute-solvent interactions in 3D RISM (“three-dimensional reference interaction site model”) integral equation theory is recast in a form that allows for a computational treatment analogous to the “particle-mesh Ewald” formalism as used for molecular simulations. In addition, relations that connect 3D RISM correlation functions and interaction potentials with thermodynamic quantities such as the chemical potential and average solute-solvent interaction energy are reformulated in a way that calculations of expensive real-space electrostatic terms on the 3D grid are completely avoided. These methodical enhancements allow for both, a significant speedup particularly for large solute systems and a smoother convergence of predicted thermodynamic quantities with respect to box size, as illustrated for several benchmark systems.
3D RISM theory with fast reciprocal-space electrostatics.
Heil, Jochen; Kast, Stefan M
2015-03-21
The calculation of electrostatic solute-solvent interactions in 3D RISM ("three-dimensional reference interaction site model") integral equation theory is recast in a form that allows for a computational treatment analogous to the "particle-mesh Ewald" formalism as used for molecular simulations. In addition, relations that connect 3D RISM correlation functions and interaction potentials with thermodynamic quantities such as the chemical potential and average solute-solvent interaction energy are reformulated in a way that calculations of expensive real-space electrostatic terms on the 3D grid are completely avoided. These methodical enhancements allow for both, a significant speedup particularly for large solute systems and a smoother convergence of predicted thermodynamic quantities with respect to box size, as illustrated for several benchmark systems. PMID:25796231
3D imaging of dental hard tissues with Fourier domain optical coherence tomography
NASA Astrophysics Data System (ADS)
Chen, Yueli L.; Yang, Yi; Ma, Jing; Yan, Jun; Shou, Yuanxin; Wang, Tianheng; Ramesh, Aruna; Zhao, Jing; Zhu, Quing
2011-03-01
A fiber optical coherence tomography (OCT) probe is used for three dimensional dental imaging. The probe has a lightweight miniaturized design with a size of a pen to facilitate clinic in vivo diagnostics. The probe is interfaced with a swept-source / Fourier domain optical coherence tomography at 20K axial scanning rate. The tooth samples were scanned from occlusal, buccal, lingual, mesial, and distal orientations. Three dimensional imaging covers tooth surface area up to 10 mm x 10 mm with a depth about 5 mm, where a majority of caries affection occurs. OCT image provides better resolution and contrast compared to gold standard dental radiography (X-ray). In particular, the technology is well suited for occlusal caries detection. This is complementary to X-ray as occlusal caries affection is difficult to be detected due to the X-ray projectile scan geometry. The 3D topology of occlusal surface as well as the dentin-enamel junction (DEJ) surface inside the tooth can be visualized. The lesion area appears with much stronger back scattering signal intensity.
Fast vision-based catheter 3D reconstruction.
Moradi Dalvand, Mohsen; Nahavandi, Saeid; Howe, Robert D
2016-07-21
Continuum robots offer better maneuverability and inherent compliance and are well-suited for surgical applications as catheters, where gentle interaction with the environment is desired. However, sensing their shape and tip position is a challenge as traditional sensors can not be employed in the way they are in rigid robotic manipulators. In this paper, a high speed vision-based shape sensing algorithm for real-time 3D reconstruction of continuum robots based on the views of two arbitrary positioned cameras is presented. The algorithm is based on the closed-form analytical solution of the reconstruction of quadratic curves in 3D space from two arbitrary perspective projections. High-speed image processing algorithms are developed for the segmentation and feature extraction from the images. The proposed algorithms are experimentally validated for accuracy by measuring the tip position, length and bending and orientation angles for known circular and elliptical catheter shaped tubes. Sensitivity analysis is also carried out to evaluate the robustness of the algorithm. Experimental results demonstrate good accuracy (maximum errors of ±0.6 mm and ±0.5 deg), performance (200 Hz), and robustness (maximum absolute error of 1.74 mm, 3.64 deg for the added noises) of the proposed high speed algorithms. PMID:27352011
Fast 3D fluid registration of brain magnetic resonance images
NASA Astrophysics Data System (ADS)
Leporé, Natasha; Chou, Yi-Yu; Lopez, Oscar L.; Aizenstein, Howard J.; Becker, James T.; Toga, Arthur W.; Thompson, Paul M.
2008-03-01
Fluid registration is widely used in medical imaging to track anatomical changes, to correct image distortions, and to integrate multi-modality data. Fluid mappings guarantee that the template image deforms smoothly into the target, without tearing or folding, even when large deformations are required for accurate matching. Here we implemented an intensity-based fluid registration algorithm, accelerated by using a filter designed by Bro-Nielsen and Gramkow. We validated the algorithm on 2D and 3D geometric phantoms using the mean square difference between the final registered image and target as a measure of the accuracy of the registration. In tests on phantom images with different levels of overlap, varying amounts of Gaussian noise, and different intensity gradients, the fluid method outperformed a more commonly used elastic registration method, both in terms of accuracy and in avoiding topological errors during deformation. We also studied the effect of varying the viscosity coefficients in the viscous fluid equation, to optimize registration accuracy. Finally, we applied the fluid registration algorithm to a dataset of 2D binary corpus callosum images and 3D volumetric brain MRIs from 14 healthy individuals to assess its accuracy and robustness.
Fast vision-based catheter 3D reconstruction
NASA Astrophysics Data System (ADS)
Moradi Dalvand, Mohsen; Nahavandi, Saeid; Howe, Robert D.
2016-07-01
Continuum robots offer better maneuverability and inherent compliance and are well-suited for surgical applications as catheters, where gentle interaction with the environment is desired. However, sensing their shape and tip position is a challenge as traditional sensors can not be employed in the way they are in rigid robotic manipulators. In this paper, a high speed vision-based shape sensing algorithm for real-time 3D reconstruction of continuum robots based on the views of two arbitrary positioned cameras is presented. The algorithm is based on the closed-form analytical solution of the reconstruction of quadratic curves in 3D space from two arbitrary perspective projections. High-speed image processing algorithms are developed for the segmentation and feature extraction from the images. The proposed algorithms are experimentally validated for accuracy by measuring the tip position, length and bending and orientation angles for known circular and elliptical catheter shaped tubes. Sensitivity analysis is also carried out to evaluate the robustness of the algorithm. Experimental results demonstrate good accuracy (maximum errors of ±0.6 mm and ±0.5 deg), performance (200 Hz), and robustness (maximum absolute error of 1.74 mm, 3.64 deg for the added noises) of the proposed high speed algorithms.
XuvTools: free, fast and reliable stitching of large 3D datasets.
Emmenlauer, M; Ronneberger, O; Ponti, A; Schwarb, P; Griffa, A; Filippi, A; Nitschke, R; Driever, W; Burkhardt, H
2009-01-01
Current biomedical research increasingly requires imaging large and thick 3D structures at high resolution. Prominent examples are the tracking of fine filaments over long distances in brain slices, or the localization of gene expression or cell migration in whole animals like Caenorhabditis elegans or zebrafish. To obtain both high resolution and a large field of view (FOV), a combination of multiple recordings ('tiles') is one of the options. Although hardware solutions exist for fast and reproducible acquisition of multiple 3D tiles, generic software solutions are missing to assemble ('stitch') these tiles quickly and accurately. In this paper, we present a framework that achieves fully automated recombination of tiles recorded at arbitrary positions in 3D space, as long as some small overlap between tiles is provided. A fully automated 3D correlation between all tiles is achieved such that no manual interaction or prior knowledge about tile positions is needed. We use (1) phase-only correlation in a multi-scale approach to estimate the coarse positions, (2) normalized cross-correlation of small patches extracted at salient points to obtain the precise matches, (3) find the globally optimal placement for all tiles by a singular value decomposition and (4) accomplish a nearly seamless stitching by a bleaching correction at the tile borders. If the dataset contains multiple channels, all channels are used to obtain the best matches between tiles. For speedup we employ a heuristic method to prune unneeded correlations, and compute all correlations via the fast Fourier transform (FFT), thereby achieving very good runtime performance. We demonstrate the successful application of the proposed framework to a wide range of different datasets from whole zebrafish embryos and C. elegans, mouse and rat brain slices and fine plant hairs (trichome). Further, we compare our stitching results to those of other commercially and freely available software solutions. The
Fast 3D Surface Extraction 2 pages (including abstract)
Sewell, Christopher Meyer; Patchett, John M.; Ahrens, James P.
2012-06-05
Ocean scientists searching for isosurfaces and/or thresholds of interest in high resolution 3D datasets required a tedious and time-consuming interactive exploration experience. PISTON research and development activities are enabling ocean scientists to rapidly and interactively explore isosurfaces and thresholds in their large data sets using a simple slider with real time calculation and visualization of these features. Ocean Scientists can now visualize more features in less time, helping them gain a better understanding of the high resolution data sets they work with on a daily basis. Isosurface timings (512{sup 3} grid): VTK 7.7 s, Parallel VTK (48-core) 1.3 s, PISTON OpenMP (48-core) 0.2 s, PISTON CUDA (Quadro 6000) 0.1 s.
Scanning fringe projection for fast 3D inspection
NASA Astrophysics Data System (ADS)
Honegger, Marc; Kahl, Michael; Trunz, Sandra; Rinner, Stefan; Ettemeyer, Andreas; Lambelet, Patrick
2013-04-01
In an earlier paper we have described a concept for high speed 3D inspection using fringe projection techniques. We use a special CMOS camera with 300 x 300 px which can calculate the phase on board. The focus of the first step of development had been a fringe projector, which was able to modulate the projected fringes with up to 250 kHz. In the second step the image acquisition part of the system was developed. In case of 3D measurement with a matrix camera, the camera resp. the measuring object has to be moved stepwise in the lateral direction to cover multiple acquisition areas of the measurement object. Between each image the lateral movement has to correspond to the field of view of the camera. At the intended very high image acquisition rates the high acceleration of the system between each image will lead to inacceptable mechanical forces. In order to obtain a continuous scanning procedure and at the same time to use the performance of a matrix camera, a special lens system was developed. A measurement field 120 mm long and 3 mm wide is imaged onto the camera. The width of the measuring field is imaged onto the 300 rows of the camera, giving a lateral resolution of 10 μm. In the longitudinal direction the 120 mm object length is divided into 12'000 lines to reach the same resolution of 10 μm. Due to technical reasons that will be described later only 240 of the 300 pixel rows on the chip were used. Consequently, with each camera shot 240 separated lines are imaged onto the chip. Between each of these 240 lines there is a large empty space, which is not imaged by the camera. In principle, the camera is operating as 240 single line cameras. Therefore, if the camera is shifted in an inclined direction to the camera orientation over the object, the empty spaces can be recorded as well. In an optimum alignment, the complete measuring object can be scanned in a continuous movement, covering the total length of 120 mm. In this paper we will describe this image
Bouallègue, Fayçal Ben; Crouzet, Jean-François; Comtat, Claude; Fourcade, Marjolaine; Mohammadi, Bijan; Mariano-Goulart, Denis
2007-07-01
This paper presents an extended 3-D exact rebinning formula in the Fourier space that leads to an iterative reprojection algorithm (iterative FOREPROJ), which enables the estimation of unmeasured oblique projection data on the basis of the whole set of measured data. In first approximation, this analytical formula also leads to an extended Fourier rebinning equation that is the basis for an approximate reprojection algorithm (extended FORE). These algorithms were evaluated on numerically simulated 3-D positron emission tomography (PET) data for the solution of the truncation problem, i.e., the estimation of the missing portions in the oblique projection data, before the application of algorithms that require complete projection data such as some rebinning methods (FOREX) or 3-D reconstruction algorithms (3DRP or direct Fourier methods). By taking advantage of all the 3-D data statistics, the iterative FOREPROJ reprojection provides a reliable alternative to the classical FOREPROJ method, which only exploits the low-statistics nonoblique data. It significantly improves the quality of the external reconstructed slices without loss of spatial resolution. As for the approximate extended FORE algorithm, it clearly exhibits limitations due to axial interpolations, but will require clinical studies with more realistic measured data in order to decide on its pertinence. PMID:17649913
Temporal-spatial modeling of fast-moving and deforming 3D objects
NASA Astrophysics Data System (ADS)
Wu, Xiaoliang; Wei, Youzhi
1998-09-01
This paper gives a brief description of the method and techniques developed for the modeling and reconstruction of fast moving and deforming 3D objects. A new approach using close-range digital terrestrial photogrammetry in conjunction with high speed photography and videography is proposed. A sequential image matching method (SIM) has been developed to automatically process pairs of images taken continuously of any fast moving and deforming 3D objects. Using the SIM technique a temporal-spatial model (TSM) of any fast moving and deforming 3D objects can be developed. The TSM would include a series of reconstructed surface models of the fast moving and deforming 3D object in the form of 3D images. The TSM allows the 3D objects to be visualized and analyzed in sequence. The SIM method, specifically the left-right matching and forward-back matching techniques are presented in the paper. An example is given which deals with the monitoring of a typical blast rock bench in a major open pit mine in Australia. With the SIM approach and the TSM model it is possible to automatically and efficiently reconstruct the 3D images of the blasting process. This reconstruction would otherwise be impossible to achieve using a labor intensive manual processing approach based on 2D images taken from conventional high speed cameras. The case study demonstrates the potential of the SIM approach and the TSM for the automatic identification, tracking and reconstruction of any fast moving and deforming 3D targets.
Fast 3D subsurface imaging with stepped-frequency GPR
NASA Astrophysics Data System (ADS)
Masarik, Matthew P.; Burns, Joseph; Thelen, Brian T.; Sutter, Lena
2015-05-01
This paper investigates an algorithm for forming 3D images of the subsurface using stepped-frequency GPR data. The algorithm is specifically designed for a handheld GPR and therefore accounts for the irregular sampling pattern in the data and the spatially-variant air-ground interface by estimating an effective "ground-plane" and then registering the data to the plane. The algorithm efficiently solves the 4th-order polynomial for the Snell reflection points using a fully vectorized iterative scheme. The forward operator is implemented efficiently using an accelerated nonuniform FFT (Greengard and Lee, 2004); the adjoint operator is implemented efficiently using an interpolation step coupled with an upsampled FFT. The imaging is done as a linearized version of the full inverse problem, which is regularized using a sparsity constraint to reduce sidelobes and therefore improve image localization. Applying an appropriate sparsity constraint, the algorithm is able to eliminate most the surrounding clutter and sidelobes, while still rendering valuable image properties such as shape and size. The algorithm is applied to simulated data, controlled experimental data (made available by Dr. Waymond Scott, Georgia Institute of Technology), and government-provided data with irregular sampling and air-ground interface.
Fast algorithm of 3D median filter for medical image despeckling
NASA Astrophysics Data System (ADS)
Xiong, Chengyi; Hou, Jianhua; Gao, Zhirong; He, Xiang; Chen, Shaoping
2007-12-01
Three-dimensional (3-D) median filtering is very useful to eliminate speckle noise from a medical imaging source, such as functional magnetic resonance imaging (fMRI) and ultrasonic imaging. 3-D median filtering is characterized by its higher computation complexity. N 3(N 3-1)/2 comparison operations would be required for 3-D median filtering with N×N×N window if the conventional bubble-sorting algorithm is adopted. In this paper, an efficient fast algorithm for 3-D median filtering was presented, which considerably reduced the computation complexity for extracting the median of a 3-D data array. Compared to the state-of-the-art, the proposed method could reduce the computation complexity of 3-D median filtering by 33%. It results in efficiently reducing the system delay of the 3-D median filter by software implementation, and the system cost and power consumption by hardware implementation.
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.
NASA Astrophysics Data System (ADS)
Fernandez, Sergio; Gdeisat, Munther A.; Salvi, Joaquim; Burton, David
2011-06-01
Fringe pattern analysis in coded structured light constitutes an active field of research. Techniques based on first projecting a sinusoidal pattern and then recovering the phase deviation permit the computation of the phase map and its corresponding depth map, leading to a dense acquisition of the measuring object. Among these techniques, the ones based on time-frequency analysis permit to extract the depth map from a single image, thus having potential applications measuring moving objects. The main techniques are Fourier Transform (FT), Windowed Fourier Transform (WFT) and Wavelet Transform (WT). This paper first analyzes the pros and cons of these three techniques, then a new algorithm for the automatic selection of the window size in WFT is proposed. This algorithm is compared to the traditional WT using adapted mother wavelet signals both with simulated and real objects, showing the performance results for quantitative and qualitative evaluations of the new method.
A fast 3D reconstruction system with a low-cost camera accessory
Zhang, Yiwei; Gibson, Graham M.; Hay, Rebecca; Bowman, Richard W.; Padgett, Miles J.; Edgar, Matthew P.
2015-01-01
Photometric stereo is a three dimensional (3D) imaging technique that uses multiple 2D images, obtained from a fixed camera perspective, with different illumination directions. Compared to other 3D imaging methods such as geometry modeling and 3D-scanning, it comes with a number of advantages, such as having a simple and efficient reconstruction routine. In this work, we describe a low-cost accessory to a commercial digital single-lens reflex (DSLR) camera system allowing fast reconstruction of 3D objects using photometric stereo. The accessory consists of four white LED lights fixed to the lens of a commercial DSLR camera and a USB programmable controller board to sequentially control the illumination. 3D images are derived for different objects with varying geometric complexity and results are presented, showing a typical height error of <3 mm for a 50 mm sized object. PMID:26057407
A fast 3D reconstruction system with a low-cost camera accessory
NASA Astrophysics Data System (ADS)
Zhang, Yiwei; Gibson, Graham M.; Hay, Rebecca; Bowman, Richard W.; Padgett, Miles J.; Edgar, Matthew P.
2015-06-01
Photometric stereo is a three dimensional (3D) imaging technique that uses multiple 2D images, obtained from a fixed camera perspective, with different illumination directions. Compared to other 3D imaging methods such as geometry modeling and 3D-scanning, it comes with a number of advantages, such as having a simple and efficient reconstruction routine. In this work, we describe a low-cost accessory to a commercial digital single-lens reflex (DSLR) camera system allowing fast reconstruction of 3D objects using photometric stereo. The accessory consists of four white LED lights fixed to the lens of a commercial DSLR camera and a USB programmable controller board to sequentially control the illumination. 3D images are derived for different objects with varying geometric complexity and results are presented, showing a typical height error of <3 mm for a 50 mm sized object.
In Vivo 3D Meibography of the Human Eyelid Using Real Time Imaging Fourier-Domain OCT
Hwang, Ho Sik; Shin, Jun Geun; Lee, Byeong Ha; Eom, Tae Joong; Joo, Choun-Ki
2013-01-01
Recently, we reported obtaining tomograms of meibomian glands from healthy volunteers using commercial anterior segment optical coherence tomography (AS-OCT), which is widely employed in clinics for examination of the anterior segment. However, we could not create 3D images of the meibomian glands, because the commercial OCT does not have a 3D reconstruction function. In this study we report the creation of 3D images of the meibomian glands by reconstructing the tomograms of these glands using high speed Fourier-Domain OCT (FD-OCT) developed in our laboratory. This research was jointly undertaken at the Department of Ophthalmology, Seoul St. Mary's Hospital (Seoul, Korea) and the Advanced Photonics Research Institute of Gwangju Institute of Science and Technology (Gwangju, Korea) with two healthy volunteers and seven patients with meibomian gland dysfunction. A real time imaging FD-OCT system based on a high-speed wavelength swept laser was developed that had a spectral bandwidth of 100 nm at the 1310 nm center wavelength. The axial resolution was 5 µm and the lateral resolution was 13 µm in air. Using this device, the meibomian glands of nine subjects were examined. A series of tomograms from the upper eyelid measuring 5 mm (from left to right, B-scan) × 2 mm (from upper part to lower part, C-scan) were collected. Three-D images of the meibomian glands were then reconstructed using 3D “data visualization, analysis, and modeling software”. Established infrared meibography was also performed for comparison. The 3D images of healthy subjects clearly showed the meibomian glands, which looked similar to bunches of grapes. These results were consistent with previous infrared meibography results. The meibomian glands were parallel to each other, and the saccular acini were clearly visible. Here we report the successful production of 3D images of human meibomian glands by reconstructing tomograms of these glands with high speed FD-OCT. PMID:23805297
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.
Real-time 3D Fourier-domain optical coherence tomography guided microvascular anastomosis
NASA Astrophysics Data System (ADS)
Huang, Yong; Ibrahim, Zuhaib; Lee, W. P. A.; Brandacher, Gerald; Kang, Jin U.
2013-03-01
Vascular and microvascular anastomosis is considered to be the foundation of plastic and reconstructive surgery, hand surgery, transplant surgery, vascular surgery and cardiac surgery. In the last two decades innovative techniques, such as vascular coupling devices, thermo-reversible poloxamers and suture-less cuff have been introduced. Intra-operative surgical guidance using a surgical imaging modality that provides in-depth view and 3D imaging can improve outcome following both conventional and innovative anastomosis techniques. Optical coherence tomography (OCT) is a noninvasive high-resolution (micron level), high-speed, 3D imaging modality that has been adopted widely in biomedical and clinical applications. In this work we performed a proof-of-concept evaluation study of OCT as an assisted intraoperative and post-operative imaging modality for microvascular anastomosis of rodent femoral vessels. The OCT imaging modality provided lateral resolution of 12 μm and 3.0 μm axial resolution in air and 0.27 volume/s imaging speed, which could provide the surgeon with clearly visualized vessel lumen wall and suture needle position relative to the vessel during intraoperative imaging. Graphics processing unit (GPU) accelerated phase-resolved Doppler OCT (PRDOCT) imaging of the surgical site was performed as a post-operative evaluation of the anastomosed vessels and to visualize the blood flow and thrombus formation. This information could help surgeons improve surgical precision in this highly challenging anastomosis of rodent vessels with diameter less than 0.5 mm. Our imaging modality could not only detect accidental suture through the back wall of lumen but also promptly diagnose and predict thrombosis immediately after reperfusion. Hence, real-time OCT can assist in decision-making process intra-operatively and avoid post-operative complications.
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.
NASA Astrophysics Data System (ADS)
Lanusse, F.; Rassat, A.; Starck, J.-L.
2015-06-01
Context. Upcoming spectroscopic galaxy surveys are extremely promising to help in addressing the major challenges of cosmology, in particular in understanding the nature of the dark universe. The strength of these surveys, naturally described in spherical geometry, comes from their unprecedented depth and width, but an optimal extraction of their three-dimensional information is of utmost importance to best constrain the properties of the dark universe. Aims: Although there is theoretical motivation and novel tools to explore these surveys using the 3D spherical Fourier-Bessel (SFB) power spectrum of galaxy number counts Cℓ(k,k'), most survey optimisations and forecasts are based on the tomographic spherical harmonics power spectrum C(ij)_ℓ. The goal of this paper is to perform a new investigation of the information that can be extracted from these two analyses in the context of planned stage IV wide-field galaxy surveys. Methods: We compared tomographic and 3D SFB techniques by comparing the forecast cosmological parameter constraints obtained from a Fisher analysis. The comparison was made possible by careful and coherent treatment of non-linear scales in the two analyses, which makes this study the first to compare 3D SFB and tomographic constraints on an equal footing. Nuisance parameters related to a scale- and redshift-dependent galaxy bias were also included in the computation of the 3D SFB and tomographic power spectra for the first time. Results: Tomographic and 3D SFB methods can recover similar constraints in the absence of systematics. This requires choosing an optimal number of redshift bins for the tomographic analysis, which we computed to be N = 26 for zmed ≃ 0.4, N = 30 for zmed ≃ 1.0, and N = 42 for zmed ≃ 1.7. When marginalising over nuisance parameters related to the galaxy bias, the forecast 3D SFB constraints are less affected by this source of systematics than the tomographic constraints. In addition, the rate of increase of the
Frank, Lawrence R.; Jung, Youngkyoo; Inati, Souheil; Tyszka, J. Michael; Wong, Eric C.
2009-01-01
We present an acquisition and reconstruction method designed to acquire high resolution 3D fast spin echo diffusion tensor images while mitigating the major sources of artifacts in DTI - field distortions, eddy currents and motion. The resulting images, being 3D, are of high SNR, and being fast spin echoes, exhibit greatly reduced field distortions. This sequence utilizes variable density spiral acquisition gradients, which allow for the implementation of a self-navigation scheme by which both eddy current and motion artifacts are removed. The result is that high resolution 3D DTI images are produced without the need for eddy current compensating gradients or B0 field correction. In addition, a novel method for fast and accurate reconstruction of the non-Cartesian data is employed. Results are demonstrated in the brains of normal human volunteers. PMID:19778618
A fast rebinning algorithm for 3D positron emission tomography using John's equation
NASA Astrophysics Data System (ADS)
Defrise, Michel; Liu, Xuan
1999-08-01
Volume imaging in positron emission tomography (PET) requires the inversion of the three-dimensional (3D) x-ray transform. The usual solution to this problem is based on 3D filtered-backprojection (FBP), but is slow. Alternative methods have been proposed which factor the 3D data into independent 2D data sets corresponding to the 2D Radon transforms of a stack of parallel slices. Each slice is then reconstructed using 2D FBP. These so-called rebinning methods are numerically efficient but are approximate. In this paper a new exact rebinning method is derived by exploiting the fact that the 3D x-ray transform of a function is the solution to the second-order partial differential equation first studied by John. The method is proposed for two sampling schemes, one corresponding to a pair of infinite plane detectors and another one corresponding to a cylindrical multi-ring PET scanner. The new FORE-J algorithm has been implemented for this latter geometry and was compared with the approximate Fourier rebinning algorithm FORE and with another exact rebinning algorithm, FOREX. Results with simulated data demonstrate a significant improvement in accuracy compared to FORE, while the reconstruction time is doubled. Compared to FOREX, the FORE-J algorithm is slightly less accurate but more than three times faster.
A fast method to measure the 3D surface of the human heart
NASA Astrophysics Data System (ADS)
Cao, Yiping; Su, Xianyu; Xiang, Liqun; Chen, Wenjing; Zhang, Qican
2003-12-01
Three-dimensional (3-D) automatic measurement of an object is widely used in many fields. In Biology and Medicine society, it can be applicable for surgery, orthopedics, viscera disease analysis and diagnosis etc. Here a new fast method to measure the 3D surface of human heart is proposed which can provide doctors a lot of information, such as the size of heart profile, the sizes of the left or right heart ventricle, and the curvature center and radius of heart ventricle, to fully analyze and diagnose pathobiology of human heart. The new fast method is optically and noncontacted and based upon the Phase Measurement Profilometry (PMP), which has higher measuring precision. A human heart specimen experiment has verified our method.
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.
Fast Probabilistic Fusion of 3d Point Clouds via Occupancy Grids for Scene Classification
NASA Astrophysics Data System (ADS)
Kuhn, Andreas; Huang, Hai; Drauschke, Martin; Mayer, Helmut
2016-06-01
High resolution consumer cameras on Unmanned Aerial Vehicles (UAVs) allow for cheap acquisition of highly detailed images, e.g., of urban regions. Via image registration by means of Structure from Motion (SfM) and Multi View Stereo (MVS) the automatic generation of huge amounts of 3D points with a relative accuracy in the centimeter range is possible. Applications such as semantic classification have a need for accurate 3D point clouds, but do not benefit from an extremely high resolution/density. In this paper, we, therefore, propose a fast fusion of high resolution 3D point clouds based on occupancy grids. The result is used for semantic classification. In contrast to state-of-the-art classification methods, we accept a certain percentage of outliers, arguing that they can be considered in the classification process when a per point belief is determined in the fusion process. To this end, we employ an octree-based fusion which allows for the derivation of outlier probabilities. The probabilities give a belief for every 3D point, which is essential for the semantic classification to consider measurement noise. For an example point cloud with half a billion 3D points (cf. Figure 1), we show that our method can reduce runtime as well as improve classification accuracy and offers high scalability for large datasets.
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.
3D WKB solution for fast magnetoacoustic wave behaviour around an X-line
NASA Astrophysics Data System (ADS)
McLaughlin, J. A.; Botha, G. J. J.; Régnier, S.; Spoors, D. L.
2016-06-01
Context. We study the propagation of a fast magnetoacoustic wave in a 3D magnetic field created from two magnetic dipoles. The magnetic topology contains an X-line. Aims: We aim to contribute to the overall understanding of MHD wave propagation within inhomogeneous media, specifically around X-lines. Methods: We investigate the linearised, 3D MHD equations under the assumptions of ideal and cold plasma. We utilise the WKB approximation and Charpit's method during our investigation. Results: It is found that the behaviour of the fast magnetoacoustic wave is entirely dictated by the local, inhomogeneous, equilibrium Alfvén speed profile. All parts of the wave experience refraction during propagation, where the magnitude of the refraction effect depends on the location of an individual wave element within the inhomogeneous magnetic field. The X-line, along which the Alfvén speed is identically zero, acts as a focus for the refraction effect. There are two main types of wave behaviour: part of the wave is either trapped by the X-line or escapes the system, and there exists a critical starting region around the X-line that divides these two types of behaviour. For the set-up investigated, it is found that 15.5% of the fast wave energy is trapped by the X-line. Conclusions: We conclude that linear, β = 0 fast magnetoacoustic waves can accumulate along X-lines and thus these will be specific locations of fast wave energy deposition and thus preferential heating. The work here highlights the importance of understanding the magnetic topology of a system. We also demonstrate how the 3D WKB technique described in this paper can be applied to other magnetic configurations.
Fast extraction of minimal paths in 3D images and applications to virtual endoscopy.
Deschamps, T; Cohen, L D
2001-12-01
The aim of this article is to build trajectories for virtual endoscopy inside 3D medical images, using the most automatic way. Usually the construction of this trajectory is left to the clinician who must define some points on the path manually using three orthogonal views. But for a complex structure such as the colon, those views give little information on the shape of the object of interest. The path construction in 3D images becomes a very tedious task and precise a priori knowledge of the structure is needed to determine a suitable trajectory. We propose a more automatic path tracking method to overcome those drawbacks: we are able to build a path, given only one or two end points and the 3D image as inputs. This work is based on previous work by Cohen and Kimmel [Int. J. Comp. Vis. 24 (1) (1997) 57] for extracting paths in 2D images using Fast Marching algorithm. Our original contribution is twofold. On the first hand, we present a general technical contribution which extends minimal paths to 3D images and gives new improvements of the approach that are relevant in 2D as well as in 3D to extract linear structures in images. It includes techniques to make the path extraction scheme faster and easier, by reducing the user interaction. We also develop a new method to extract a centered path in tubular structures. Synthetic and real medical images are used to illustrate each contribution. On the other hand, we show that our method can be efficiently applied to the problem of finding a centered path in tubular anatomical structures with minimum interactivity, and that this path can be used for virtual endoscopy. Results are shown in various anatomical regions (colon, brain vessels, arteries) with different 3D imaging protocols (CT, MR). PMID:11731307
Fast 3D Spatial EPR Imaging Using Spiral Magnetic Field Gradient
Deng, Yuanmu; Petryakov, Sergy; He, Guanglong; Kesselring, Eric; Kuppusamy, Periannan; Zweier, Jay L.
2007-01-01
Electron paramagnetic resonance imaging (EPRI) provides direct detection and mapping of free radicals. The continuous wave (CW) EPRI technique, in particular, has been widely used in a variety of applications in the fields of biology and medicine due to its high sensitivity and applicability to a wide range of free radicals and paramagnetic species. However, the technique requires long image acquisition periods, and this limits its use for many in vivo applications where relatively rapid changes occur in the magnitude and distribution of spins. Therefore, there has been a great need to develop fast EPRI techniques. We report the development of a fast 3D CW EPRI technique using spiral magnetic field gradient. By spiraling the magnetic field gradient and stepping the main magnetic field, this approach acquires a 3D image in one sweep of the main magnetic field, enabling significant reduction of the imaging time. A direct one-stage 3D image reconstruction algorithm, modified for reconstruction of the EPR images from the projections acquired with the spiral magnetic field gradient, was used. We demonstrated using a home-built L-band EPR system that the spiral magnetic field gradient technique enabled a 4 to 7-fold accelerated acquisition of projections. This technique has great potential for in vivo studies of free radicals and their metabolism. PMID:17267252
Comparison of Parallel MRI Reconstruction Methods for Accelerated 3D Fast Spin-Echo Imaging
Xiao, Zhikui; Hoge, W. Scott; Mulkern, R.V.; Zhao, Lei; Hu, Guangshu; Kyriakos, Walid E.
2014-01-01
Parallel MRI (pMRI) achieves imaging acceleration by partially substituting gradient-encoding steps with spatial information contained in the component coils of the acquisition array. Variable-density subsampling in pMRI was previously shown to yield improved two-dimensional (2D) imaging in comparison to uniform subsampling, but has yet to be used routinely in clinical practice. In an effort to reduce acquisition time for 3D fast spin-echo (3D-FSE) sequences, this work explores a specific nonuniform sampling scheme for 3D imaging, subsampling along two phase-encoding (PE) directions on a rectilinear grid. We use two reconstruction methods—2D-GRAPPA-Operator and 2D-SPACE RIP—and present a comparison between them. We show that high-quality images can be reconstructed using both techniques. To evaluate the proposed sampling method and reconstruction schemes, results via simulation, phantom study, and in vivo 3D human data are shown. We find that fewer artifacts can be seen in the 2D-SPACE RIP reconstructions than in 2D-GRAPPA-Operator reconstructions, with comparable reconstruction times. PMID:18727083
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).
Algorithms for Accurate and Fast Plotting of Contour Surfaces in 3D Using Hexahedral Elements
NASA Astrophysics Data System (ADS)
Singh, Chandan; Saini, Jaswinder Singh
2016-07-01
In the present study, Fast and accurate algorithms for the generation of contour surfaces in 3D are described using hexahedral elements which are popular in finite element analysis. The contour surfaces are described in the form of groups of boundaries of contour segments and their interior points are derived using the contour equation. The locations of contour boundaries and the interior points on contour surfaces are as accurate as the interpolation results obtained by hexahedral elements and thus there are no discrepancies between the analysis and visualization results.
Algorithms for Accurate and Fast Plotting of Contour Surfaces in 3D Using Hexahedral Elements
NASA Astrophysics Data System (ADS)
Singh, Chandan; Saini, Jaswinder Singh
2016-05-01
In the present study, Fast and accurate algorithms for the generation of contour surfaces in 3D are described using hexahedral elements which are popular in finite element analysis. The contour surfaces are described in the form of groups of boundaries of contour segments and their interior points are derived using the contour equation. The locations of contour boundaries and the interior points on contour surfaces are as accurate as the interpolation results obtained by hexahedral elements and thus there are no discrepancies between the analysis and visualization results.
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.
Thermal analysis of 3D composites by a new fast multipole hybrid boundary node method
NASA Astrophysics Data System (ADS)
Miao, Yu; Wang, Qiao; Zhu, Hongping; Li, Yinping
2014-01-01
This paper applies the hybrid boundary node method (Hybrid BNM) for the thermal analysis of 3D composites. A new formulation is derived for the inclusion-based composites. In the new formulation, the unknowns of the interfaces are assembled only once in the final system equation, which can reduce nearly one half of degrees of freedom (DOFs) compared with the conventional multi-domain solver when there are lots of inclusions. A new version of the fast multipole method (FMM) is also coupled with the new formulation and the technique is applied to thermal analysis of composites with many inclusions. In the new fast multipole hybrid boundary node method (FM-HBNM), a diagonal form for translation operators is used and the method presented can be applied to the computation of more than 1,000,000 DOFs on a personal computer. Numerical examples are presented to analyze the thermal behavior of composites with many inclusions.
3D Fast Automatic Segmentation of Kidney Based on Modified AAM and Random Forest.
Jin, Chao; Shi, Fei; Xiang, Dehui; Jiang, Xueqing; Zhang, Bin; Wang, Ximing; Zhu, Weifang; Gao, Enting; Chen, Xinjian
2016-06-01
In this paper, a fully automatic method is proposed to segment the kidney into multiple components: renal cortex, renal column, renal medulla and renal pelvis, in clinical 3D CT abdominal images. The proposed fast automatic segmentation method of kidney consists of two main parts: localization of renal cortex and segmentation of kidney components. In the localization of renal cortex phase, a method which fully combines 3D Generalized Hough Transform (GHT) and 3D Active Appearance Models (AAM) is applied to localize the renal cortex. In the segmentation of kidney components phase, a modified Random Forests (RF) method is proposed to segment the kidney into four components based on the result from localization phase. During the implementation, a multithreading technology is applied to speed up the segmentation process. The proposed method was evaluated on a clinical abdomen CT data set, including 37 contrast-enhanced volume data using leave-one-out strategy. The overall true-positive volume fraction and false-positive volume fraction were 93.15%, 0.37% for renal cortex segmentation; 83.09%, 0.97% for renal column segmentation; 81.92%, 0.55% for renal medulla segmentation; and 80.28%, 0.30% for renal pelvis segmentation, respectively. The average computational time of segmenting kidney into four components took 20 seconds. PMID:26742124
Fast 3D modeling in complex environments using a single Kinect sensor
NASA Astrophysics Data System (ADS)
Yue, Haosong; Chen, Weihai; Wu, Xingming; Liu, Jingmeng
2014-02-01
Three-dimensional (3D) modeling technology has been widely used in inverse engineering, urban planning, robot navigation, and many other applications. How to build a dense model of the environment with limited processing resources is still a challenging topic. A fast 3D modeling algorithm that only uses a single Kinect sensor is proposed in this paper. For every color image captured by Kinect, corner feature extraction is carried out first. Then a spiral search strategy is utilized to select the region of interest (ROI) that contains enough feature corners. Next, the iterative closest point (ICP) method is applied to the points in the ROI to align consecutive data frames. Finally, the analysis of which areas can be walked through by human beings is presented. Comparative experiments with the well-known KinectFusion algorithm have been done and the results demonstrate that the accuracy of the proposed algorithm is the same as KinectFusion but the computing speed is nearly twice of KinectFusion. 3D modeling of two scenes of a public garden and traversable areas analysis in these regions further verified the feasibility of our algorithm.
Fast and robust 3D ultrasound registration--block and game theoretic matching.
Banerjee, Jyotirmoy; Klink, Camiel; Peters, Edward D; Niessen, Wiro J; Moelker, Adriaan; van Walsum, Theo
2015-02-01
Real-time 3D US has potential for image guidance in minimally invasive liver interventions. However, motion caused by patient breathing makes it hard to visualize a localized area, and to maintain alignment with pre-operative information. In this work we develop a fast affine registration framework to compensate in real-time for liver motion/displacement due to breathing. The affine registration of two consecutive ultrasound volumes in time is performed using block-matching. For a set of evenly distributed points in one volume and their correspondences in the other volume, we propose a robust outlier rejection method to reject false matches. The inliers are then used to determine the affine transformation. The approach is evaluated on 13 4D ultrasound sequences acquired from 8 subjects. For 91 pairs of 3D ultrasound volumes selected from these sequences, a mean registration error of 1.8mm is achieved. A graphics processing unit (GPU) implementation runs the 3D US registration at 8 Hz. PMID:25484018
Capabilities of a Global 3D MHD Model for Monitoring Extremely Fast CMEs
NASA Astrophysics Data System (ADS)
Wu, C. C.; Plunkett, S. P.; Liou, K.; Socker, D. G.; Wu, S. T.; Wang, Y. M.
2015-12-01
Since the start of the space era, spacecraft have recorded many extremely fast coronal mass ejections (CMEs) which have resulted in severe geomagnetic storms. Accurate and timely forecasting of the space weather effects of these events is important for protecting expensive space assets and astronauts and avoiding communications interruptions. Here, we will introduce a newly developed global, three-dimensional (3D) magnetohydrodynamic (MHD) model (G3DMHD). The model takes the solar magnetic field maps at 2.5 solar radii (Rs) and intepolates the solar wind plasma and field out to 18 Rs using the algorithm of Wang and Sheeley (1990, JGR). The output is used as the inner boundary condition for a 3D MHD model. The G3DMHD model is capable of simulating (i) extremely fast CME events with propagation speeds faster than 2500 km/s; and (ii) multiple CME events in sequence or simultaneously. We will demonstrate the simulation results (and comparison with in-situ observation) for the fastest CME in record on 23 July 2012, the shortest transit time in March 1976, and the well-known historic Carrington 1859 event.
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
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.
Park, Sung-Hee; Ryu, Keun Ho; Gilbert, David
2005-08-01
Similarity search for protein 3D structures become complex and computationally expensive due to the fact that the size of protein structure databases continues to grow tremendously. Recently, fast structural similarity search systems have been required to put them into practical use in protein structure classification whilst existing comparison systems do not provide comparison results on time. Our approach uses multi-step processing that composes of a preprocessing step to represent geometry of protein structures with spatial objects, a filter step to generate a small candidate set using approximate topological string matching, and a refinement step to compute a structural alignment. This paper describes the preprocessing and filtering for fast similarity search using the discovery of topological patterns of secondary structure elements based on spatial relations. Our system is fully implemented by using Oracle 8i spatial. We have previously shown that our approach has the advantage of speed of performance compared with other approach such as DALI. This work shows that the discovery of topological relations of secondary structure elements in protein structures by using spatial relations of spatial databases is practical for fast structural similarity search for proteins. PMID:16187404
Miyazaki, Jun; Iida, Tadatsune; Tanaka, Shinji; Hayashi-Takagi, Akiko; Kasai, Haruo; Okabe, Shigeo; Kobayashi, Takayoshi
2016-05-01
A fast, high-sensitivity photothermal microscope was developed by implementing a spatially segmented balanced detection scheme into a laser scanning microscope. We confirmed a 4.9 times improvement in signal-to-noise ratio in the spatially segmented balanced detection compared with that of conventional detection. The system demonstrated simultaneous bi-modal photothermal and confocal fluorescence imaging of transgenic mouse brain tissue with a pixel dwell time of 20 μs. The fluorescence image visualized neurons expressing yellow fluorescence proteins, while the photothermal signal detected endogenous chromophores in the mouse brain, allowing 3D visualization of the distribution of various features such as blood cells and fine structures probably due to lipids. This imaging modality was constructed using compact and cost-effective laser diodes, and will thus be widely useful in the life and medical sciences. PMID:27231615
Fast-response switchable lens for 3D and wearable displays.
Lee, Yun-Han; Peng, Fenglin; Wu, Shin-Tson
2016-01-25
We report a switchable lens in which a twisted nematic (TN) liquid crystal cell is utilized to control the input polarization. Different polarization state leads to different path length in the proposed optical system, which in turn results in different focal length. This type of switchable lens has advantages in fast response time, low operation voltage, and inherently lower chromatic aberration. Using a pixelated TN panel, we can create depth information to the selected pixels and thus add depth information to a 2D image. By cascading three such device structures together, we can generate 8 different focuses for 3D displays, wearable virtual/augmented reality, and other head mounted display devices. PMID:26832545
Miyazaki, Jun; Iida, Tadatsune; Tanaka, Shinji; Hayashi-Takagi, Akiko; Kasai, Haruo; Okabe, Shigeo; Kobayashi, Takayoshi
2016-01-01
A fast, high-sensitivity photothermal microscope was developed by implementing a spatially segmented balanced detection scheme into a laser scanning microscope. We confirmed a 4.9 times improvement in signal-to-noise ratio in the spatially segmented balanced detection compared with that of conventional detection. The system demonstrated simultaneous bi-modal photothermal and confocal fluorescence imaging of transgenic mouse brain tissue with a pixel dwell time of 20 μs. The fluorescence image visualized neurons expressing yellow fluorescence proteins, while the photothermal signal detected endogenous chromophores in the mouse brain, allowing 3D visualization of the distribution of various features such as blood cells and fine structures probably due to lipids. This imaging modality was constructed using compact and cost-effective laser diodes, and will thus be widely useful in the life and medical sciences. PMID:27231615
Fall, Mandiaye; Boutami, Salim; Glière, Alain; Stout, Brian; Hazart, Jerome
2013-06-01
A combination of the multilevel fast multipole method (MLFMM) and boundary element method (BEM) can solve large scale photonics problems of arbitrary geometry. Here, MLFMM-BEM algorithm based on a scalar and vector potential formulation, instead of the more conventional electric and magnetic field formulations, is described. The method can deal with multiple lossy or lossless dielectric objects of arbitrary geometry, be they nested, in contact, or dispersed. Several examples are used to demonstrate that this method is able to efficiently handle 3D photonic scatterers involving large numbers of unknowns. Absorption, scattering, and extinction efficiencies of gold nanoparticle spheres, calculated by the MLFMM, are compared with Mie's theory. MLFMM calculations of the bistatic radar cross section (RCS) of a gold sphere near the plasmon resonance and of a silica coated gold sphere are also compared with Mie theory predictions. Finally, the bistatic RCS of a nanoparticle gold-silver heterodimer calculated with MLFMM is compared with unmodified BEM calculations. PMID:24323115
ALE3D Simulation of Heating and Violence in a Fast Cookoff Experiment with LX-10
McClelland, M A; Maienschein, J L; Howard, W M; Nichols, A L; deHaven, M R; Strand, O T
2006-06-26
We performed a computational and experimental analysis of fast cookoff of LX-10 (94.7% HMX, 5.3% Viton A) confined in a 2 kbar steel tube with reinforced end caps. A Scaled-Thermal-Explosion-eXperiment (STEX) was completed in which three radiant heaters were used to heat the vessel until ignition, resulting in a moderately violent explosion after 20.4 minutes. Thermocouple measurements showed tube temperatures as high as 340 C at ignition and LX-10 surface temperatures as high as 279 C, which is near the melting point of HMX. Three micro-power radar systems were used to measure mean fragment velocities of 840 m/s. Photonics Doppler Velocimeters (PDVs) showed a rapid acceleration of fragments over 80 {micro}s. A one-dimensional ALE3D cookoff model at the vessel midplane was used to simulate the heating, thermal expansion, LX-10 decomposition composition, and closing of the gap between the HE (High Explosive) and vessel wall. Although the ALE3D simulation terminated before ignition, the model provided a good representation of heat transfer through the case and across the dynamic gap to the explosive.
ALE3D Simulation and Measurement of Violence in a Fast Cookoff Experiment with LX-10
McClelland, M A; Maienschein, J L; Howard, W M; deHaven, M R
2006-11-22
We performed a computational and experimental analysis of fast cookoff of LX-10 (94.7% HMX, 5.3% Viton A) confined in a 2 kbar steel tube with reinforced end caps. A Scaled-Thermal-Explosion-eXperiment (STEX) was completed in which three radiant heaters were used to heat the vessel until ignition, resulting in a moderately violent explosion after 20.4 minutes. Thermocouple measurements showed tube temperatures as high as 340 C at ignition and LX-10 surface temperatures as high as 279 C, which is near the melting point of HMX. Three micro-power radar systems were used to measure mean fragment velocities of 840 m/s. Photonics Doppler Velocimeters (PDVs) showed a rapid acceleration of fragments over 80 {micro}s. A one-dimensional ALE3D cookoff model at the vessel midplane was used to simulate the heating, thermal expansion, LX-10 decomposition composition, and closing of the gap between the HE (High Explosive) and vessel wall. Although the ALE3D simulation terminated before ignition, the model provided a good representation of heat transfer through the case and across the dynamic gap to the explosive.
Fast Wave Trains Associated with Solar Eruptions: Insights from 3D Thermodynamic MHD Simulations
NASA Astrophysics Data System (ADS)
Downs, C.; Liu, W.; Torok, T.; Linker, J.; Mikic, Z.; Ofman, L.
2015-12-01
EUV imaging observations during the SDO/AIA era have provided new insights into a variety of wave phenomena occurring in the low solar corona. One example is the observation of quasi-periodic, fast-propagating wave trains that are associated with solar eruptions, including flares and CMEs. While there has been considerable progress in understanding such waves from both an observational and theoretical perspective, it remains a challenge to pin down their physical origin. In this work, we detail our results from a case-study 3D thermodynamic MHD simulation of a coronal mass ejection where quasi-periodic wave trains are generated during the simulated eruption. We find a direct correlation between the onset of non-steady reconnection in the flare current sheet and the generation of quasi-periodic wave train signatures when patchy, collimated downflows interact with the flare arcade. Via forward modeling of SDO/AIA observables, we explore how the appearance of the wave trains is affected by line-of-sight integration and the multi-thermal nature of the coronal medium. We also examine how the wave trains themselves are channeled by natural waveguides formed in 3D by the non-uniform background magnetic field. While the physical association of the reconnection dynamics to the generation of quasi-periodic wave trains appears to be a compelling result, unanswered questions posed from recent observations as well as future prospects will be discussed.
NASA Astrophysics Data System (ADS)
Wang, Zhongren; Li, Bo; Zhou, Yuebin
2014-11-01
Fast 3D reconstruction of tool wear from 2D images has great importance to 3D measuring and objective evaluating tool wear condition, determining accurate tool change and insuring machined part's quality. Extracting 3D information of tool wear zone based on monocular multi-color structured light can realize fast recovery of surface topography of tool wear, which overcomes the problems of traditional methods such as solution diversity and slow convergence when using SFS method and stereo match when using 3D reconstruction from multiple images. In this paper, a kind of new multi-color structured light illuminator was put forward. An information mapping model was established among illuminator's structure parameters, surface morphology and color images. The mathematical model to reconstruct 3D morphology based on monocular multi-color structured light was presented. Experimental results show that this method is effective and efficient to reconstruct the surface morphology of tool wear zone.
ALE3D Simulation and Measurement of Violence in a Fast Cookoff Experiment for LX-10
McClelland, M A; Maienschein, J L; Howard, W M; deHaven, M R
2006-05-23
Fast cookoff is of interest in the areas of fire hazard reduction and the development of directed energy systems for defense. During a fast cookoff (thermal explosion), high heat fluxes cause rapid temperature increases and ignition in thin boundary layers. We are developing ALE3D models to describe the thermal, chemical, and mechanical behavior during the heating, ignition, and explosive phases. The candidate models and numerical strategies are being evaluated using benchmark cookoff experiments. Fast cookoff measurements were made in a Scaled-Thermal-Explosion-eXperiment (STEX) for LX-10 (94.7% HMX, 5.3% Viton A) confined in a 4130 steel tube with reinforced end caps. Gaps were present at the side and top of the explosive charge to allow for thermal expansion. The explosive was heated until explosion using radiant heaters. Temperatures were measured using thermocouples positioned on the tube wall and in the explosive. During the explosion, the tube expansion and fragment velocities were measured with strain gauges, Photonic-Doppler-Velocimeters (PDVs), and micropower radar units. A fragment size distribution was constructed from fragments captured in Lexan panels. ALE3D models for chemical, thermal, and mechanical behavior were developed for the heating and explosive processes. A multi-step chemical kinetics model is employed for the HMX while a one-step model is used for the Viton. A pressure-dependent deflagration model is employed during the expansion. A Steinberg-Guinan model represents the mechanical behavior of the solid constituents while polynomial and gamma-law expressions are used for the equation of state of the solid and gas species, respectively. Parameters for the kinetics model were specified using measurements of the One-Dimensional-Time-to-Explosion (ODTX), while measurements for burn rate were employed to determine parameters in the burn front model. The simulations include radiative and conductive transport across the dynamic gaps between the
Fast 3D Net Expeditions: Tools for Effective Scientific Collaboration on the World Wide Web
NASA Technical Reports Server (NTRS)
Watson, Val; Chancellor, Marisa K. (Technical Monitor)
1996-01-01
Two new technologies, the FASTexpedition and Remote FAST, have been developed that provide remote, 3D (three dimensional), high resolution, dynamic, interactive viewing of scientific data. The FASTexpedition permits one to access scientific data from the World Wide Web, take guided expeditions through the data, and continue with self controlled expeditions through the data. Remote FAST permits collaborators at remote sites to simultaneously view an analysis of scientific data being controlled by one of the collaborators. Control can be transferred between sites. These technologies are now being used for remote collaboration in joint university, industry, and NASA projects. Also, NASA Ames Research Center has initiated a project to make scientific data and guided expeditions through the data available as FASTexpeditions on the World Wide Web for educational purposes. Previously, remote visualization of dynamic data was done using video format (transmitting pixel information) such as video conferencing or MPEG (Motion Picture Expert Group) movies on the Internet. The concept for this new technology is to send the raw data (e.g., grids, vectors, and scalars) along with viewing scripts over the Internet and have the pixels generated by a visualization tool running on the viewers local workstation. The visualization tool that is currently used is FAST (Flow Analysis Software Toolkit). The advantages of this new technology over using video format are: (1) The visual is much higher in resolution (1280x1024 pixels with 24 bits of color) than typical video format transmitted over the network. (2) The form of the visualization can be controlled interactively (because the viewer is interactively controlling the visualization tool running on his workstation). (3) A rich variety of guided expeditions through the data can be included easily. (4) A capability is provided for other sites to see a visual analysis of one site as the analysis is interactively performed. Control of
3-D Ground Displacement Monitoring of very fast-moving Landslides in Emergency Scenario
NASA Astrophysics Data System (ADS)
Casu, Francesco; Manconi, Andrea; Bonano, Manuela; De Luca, Claudio; Elefante, Stefano
2014-05-01
On December 3rd, 2013, a large and fast-moving landslide phenomena, which occurred in South-West of Montescaglioso town (southern Italy) after some days of intense raining, caused ground displacements on the order of several meters. The mass wasting involved an important freeway connection disrupting more than 500 meters of the route and some isolated buildings. In this work we present a case study of application of SAR remote sensing techniques for retrieving ground displacement field in a landslide emergency scenario. To this aim, thanks to the availability of ascending and descending COSMO-SkyMed (CSK) satellite acquisitions, we first applied the DInSAR technique (Massonnet et al., 1993) to both datasets, for generating differential interferograms across the investigated event. In particular, two data pairs (one ascending and one descending) involving pre- and post-event epochs and approximately spanning the same time interval were identified. Unfortunately, the DInSAR analysis produced unsatisfactory results, because of the excessive phase noise within the area of interest, mainly related to the fast-moving deformation pattern (several meters) and also to the presence of vegetation. To overcome the above mentioned limitations, the amplitude-based Pixel Offset (PO) technique (Fialko and Simons, 2001) was applied to the previous identified CSK data pairs. In this case, the PO technique allowed us to retrieve the projection of the surface displacements across and along the satellite's track (range and azimuth, respectively) for both the ascending and descending orbits. Then, by properly combining these 2-D maps of the measured surface movements, we also retrieved the 3-D ground deformation pattern, i.e. the North, East and Vertical displacement components. The ground displacements have a main SSE component, with values exceeding 10 meters. Moreover, large subsidence values were identified in those areas experiencing the largest damages, as well as a clear uplift
Fast 3D inversion of airborne gravity-gradiometry data using Lanczos bidiagonalization method
NASA Astrophysics Data System (ADS)
Meng, Zhaohai; Li, Fengting; Zhang, Dailei; Xu, Xuechun; Huang, Danian
2016-09-01
We developed a new fast inversion method for to process and interpret airborne gravity gradiometry data, which was based on Lanczos bidiagonalization algorithm. Here, we describe the application of this new 3D gravity gradiometry inversion method to recover a subsurface density distribution model from the airborne measured gravity gradiometry anomalies. For this purpose, the survey area is divided into a large number of rectangular cells with each cell possessing a constant unknown density. It is well known that the solution of large linear gravity gradiometry is an ill-posed problem since using the smoothest inversion method is considerably time consuming. We demonstrate that the Lanczos bidiagonalization method can be an appropriate algorithm to solve a Tikhonov solver time cost function for resolving the large equations within a short time. Lanczos bidiagonalization is designed to make the very large gravity gradiometry forward modeling matrices to become low-rank, which will considerably reduce the running time of the inversion method. We also use a weighted generalized cross validation method to choose the appropriate Tikhonov parameter to improve inversion results. The inversion incorporates a model norm that allows us to attain the smoothing and depth of the solution; in addition, the model norm counteracts the natural decay of the kernels, which concentrate at shallow depths. The method is applied on noise-contaminated synthetic gravity gradiometry data to demonstrate its suitability for large 3D gravity gradiometry data inversion. The airborne gravity gradiometry data from the Vinton Salt Dome, USE, were considered as a case study. The validity of the new method on real data is discussed with reference to the Vinton Dome inversion result. The intermediate density values in the constructed model coincide well with previous results and geological information. This demonstrates the validity of the gravity gradiometry inversion method.
Numerical 3D models support two distinct hydrothermal circulation systems at fast spreading ridges
NASA Astrophysics Data System (ADS)
Hasenclever, Jörg; Theissen-Krah, Sonja; Rüpke, Lars
2013-04-01
We present 3D numerical calculations of hydrothermal fluid flow at fast spreading ridges. The setup of the 3D models is based our previous 2D studies, in which we have coupled numerical models for crustal accretion and hydrothermal fluid flow. One result of these calculations is a crustal permeability field that leads to a thermal structure in the crust that matches seismic tomography data of the East Pacific Rise (EPR). The 1000°C isotherm obtained from the 2D results is now used as the lower boundary of the 3D model domain, while the upper boundary is a smoothed bathymetry of the EPR. The same permeability field as in the 2D models is used, with the highest permeability at the ridge axis and a decrease with both depth and distance to the ridge. Permeability is also reduced linearly between 600 and 1000°C. Using a newly developed parallel finite element code written in Matlab that solves for thermal evolution, fluid pressure and Darcy flow, we simulate the flow patterns of hydrothermal circulation in a segment of 5000m along-axis, 10000m across-axis and up to 5000m depth. We observe two distinct hydrothermal circulation systems: An on-axis system forming a series of vents with a spacing ranging from 100 to 500m that is recharged by nearby (100-200m) downflows on both sides of the ridge axis. Simultaneously a second system with much broader extensions both laterally and vertically exists off-axis. It is recharged by fluids intruding between 1500m to 5000m off-axis and sampling both upper and lower crust. These fluids are channeled in the deepest and hottest regions with high permeability and migrate up-slope following the 600°C isotherm until reaching the edge of the melt lens. Depending on the width of the melt lens these off-axis fluids either merge with the on-axis hydrothermal system or form separate vents. We observe separate off-axis vent fields if the magma lens half-width exceeds 1000m and confluence of both systems for half-widths smaller than 500m. For
Fast Semantic Segmentation of 3d Point Clouds with Strongly Varying Density
NASA Astrophysics Data System (ADS)
Hackel, Timo; Wegner, Jan D.; Schindler, Konrad
2016-06-01
We describe an effective and efficient method for point-wise semantic classification of 3D point clouds. The method can handle unstructured and inhomogeneous point clouds such as those derived from static terrestrial LiDAR or photogammetric reconstruction; and it is computationally efficient, making it possible to process point clouds with many millions of points in a matter of minutes. The key issue, both to cope with strong variations in point density and to bring down computation time, turns out to be careful handling of neighborhood relations. By choosing appropriate definitions of a point's (multi-scale) neighborhood, we obtain a feature set that is both expressive and fast to compute. We evaluate our classification method both on benchmark data from a mobile mapping platform and on a variety of large, terrestrial laser scans with greatly varying point density. The proposed feature set outperforms the state of the art with respect to per-point classification accuracy, while at the same time being much faster to compute.
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
Ultra-fast hybrid CPU-GPU multiple scatter simulation for 3-D PET.
Kim, Kyung Sang; Son, Young Don; Cho, Zang Hee; Ra, Jong Beom; Ye, Jong Chul
2014-01-01
Scatter correction is very important in 3-D PET reconstruction due to a large scatter contribution in measurements. Currently, one of the most popular methods is the so-called single scatter simulation (SSS), which considers single Compton scattering contributions from many randomly distributed scatter points. The SSS enables a fast calculation of scattering with a relatively high accuracy; however, the accuracy of SSS is dependent on the accuracy of tail fitting to find a correct scaling factor, which is often difficult in low photon count measurements. To overcome this drawback as well as to improve accuracy of scatter estimation by incorporating multiple scattering contribution, we propose a multiple scatter simulation (MSS) based on a simplified Monte Carlo (MC) simulation that considers photon migration and interactions due to photoelectric absorption and Compton scattering. Unlike the SSS, the MSS calculates a scaling factor by comparing simulated prompt data with the measured data in the whole volume, which enables a more robust estimation of a scaling factor. Even though the proposed MSS is based on MC, a significant acceleration of the computational time is possible by using a virtual detector array with a larger pitch by exploiting that the scatter distribution varies slowly in spatial domain. Furthermore, our MSS implementation is nicely fit to a parallel implementation using graphic processor unit (GPU). In particular, we exploit a hybrid CPU-GPU technique using the open multiprocessing and the compute unified device architecture, which results in 128.3 times faster than using a single CPU. Overall, the computational time of MSS is 9.4 s for a high-resolution research tomograph (HRRT) system. The performance of the proposed MSS is validated through actual experiments using an HRRT. PMID:24403412
Loh, Ne-Te Duane
2011-08-01
These 2000 single-shot diffraction patterns include were either background-scattering only or hits (background-scattering plus diffraction signal from sub-micron ellipsoidal particles at random, undetermined orientations). Candidate hits were identified by eye, and the remainder were presumed as background. 54 usable, background-subtracted hits in this set (procedure in referenced article) were used to reconstruct the 3D diffraction intensities of the average ellipsoidal particle.
Fast Susceptibility-Weighted Imaging (SWI) with 3D Short-Axis Propeller (SAP)-EPI
Holdsworth, Samantha J.; Yeom, Kristen W.; Moseley, Michael E.; Skare, S.
2014-01-01
Purpose Susceptibility-Weighted Imaging (SWI) in neuroimaging can be challenging due to long scan times of 3D Gradient Recalled Echo (GRE), while faster techniques such as 3D interleaved EPI (iEPI) are prone to motion artifacts. Here we outline and implement a 3D Short-Axis Propeller Echo-Planar Imaging (SAP-EPI) trajectory as a faster, motion-correctable approach for SWI. Methods Experiments were conducted on a 3T MRI system. 3D SAP-EPI, 3D iEPI, and 3D GRE SWI scans were acquired on two volunteers. Controlled motion experiments were conducted to test the motion-correction capability of 3D SAP-EPI. 3D SAP-EPI SWI data were acquired on two pediatric patients as a potential alternative to 2D GRE used clinically. Results 3D GRE images had a better target resolution (0.47 × 0.94 × 2mm, scan time = 5min), iEPI and SAP-EPI images (resolution = 0.94 × 0.94 × 2mm) were acquired in a faster scan time (1:52min) with twice the brain coverage. SAP-EPI showed motion-correction capability and some immunity to undersampling from rejected data. Conclusion While 3D SAP-EPI suffers from some geometric distortion, its short scan time and motion-correction capability suggest that SAP-EPI may be a useful alternative to GRE and iEPI for use in SWI, particularly in uncooperative patients. PMID:24956237
Smalyuk, V.A.; Sadot, O.; Delettrez, J.A.; Meyerhofer, D.D.; Regan, S.P.; Sangster, T.C.
2005-12-05
Nonlinear growth of 3-D broadband nonuniformities was measured near saturation levels using x-ray radiography in planar foils accelerated by laser light. The initial target modulations were seeded by laser nonuniformities and later amplified during acceleration by Rayleigh-Taylor instability. The nonlinear saturation velocities are measured for the first time and are found to be in excellent agreement with Haan predictions. The measured growth of long-wavelength modes is consistent with enhanced, nonlinear, long-wavelength generation in ablatively driven targets.
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.
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
Fast and Precise 3D Fluorophore Localization based on Gradient Fitting
NASA Astrophysics Data System (ADS)
Ma, Hongqiang; Xu, Jianquan; Jin, Jingyi; Gao, Ying; Lan, Li; Liu, Yang
2015-09-01
Astigmatism imaging approach has been widely used to encode the fluorophore’s 3D position in single-particle tracking and super-resolution localization microscopy. Here, we present a new high-speed localization algorithm based on gradient fitting to precisely decode the 3D subpixel position of the fluorophore. This algebraic algorithm determines the center of the fluorescent emitter by finding the position with the best-fit gradient direction distribution to the measured point spread function (PSF), and can retrieve the 3D subpixel position of the fluorophore in a single iteration. Through numerical simulation and experiments with mammalian cells, we demonstrate that our algorithm yields comparable localization precision to the traditional iterative Gaussian function fitting (GF) based method, while exhibits over two orders-of-magnitude faster execution speed. Our algorithm is a promising high-speed analyzing method for 3D particle tracking and super-resolution localization microscopy.
Fast and Precise 3D Fluorophore Localization based on Gradient Fitting
Ma, Hongqiang; Xu, Jianquan; Jin, Jingyi; Gao, Ying; Lan, Li; Liu, Yang
2015-01-01
Astigmatism imaging approach has been widely used to encode the fluorophore’s 3D position in single-particle tracking and super-resolution localization microscopy. Here, we present a new high-speed localization algorithm based on gradient fitting to precisely decode the 3D subpixel position of the fluorophore. This algebraic algorithm determines the center of the fluorescent emitter by finding the position with the best-fit gradient direction distribution to the measured point spread function (PSF), and can retrieve the 3D subpixel position of the fluorophore in a single iteration. Through numerical simulation and experiments with mammalian cells, we demonstrate that our algorithm yields comparable localization precision to the traditional iterative Gaussian function fitting (GF) based method, while exhibits over two orders-of-magnitude faster execution speed. Our algorithm is a promising high-speed analyzing method for 3D particle tracking and super-resolution localization microscopy. PMID:26390959
Song, Liang; Maslov, Konstantin; Bitton, Rachel; Shung, K. Kirk; Wang, Lihong V.
2009-01-01
We present an in vivo dark-field reflection-mode photoacoustic microscopy system that performs cross-sectional (B-scan) imaging at 50 Hz with realtime beamforming and 3D imaging consisting of 166 B-scan frames at 1 Hz with post-beamforming. To our knowledge, this speed is currently the fastest in photoacoustic imaging. A custom-designed light delivery system is integrated with a 30-MHz ultrasound linear array to realize dark-field reflection-mode imaging. Linear mechanical scanning of the array produces 3D images. The system has axial, lateral, and elevational resolutions of 25, 70, and 200 μm, respectively, and can image 3 mm deep in scattering biological tissues. Volumetric images of subcutaneous vasculature in rats are demonstrated in vivo. Fast 3D photoacoustic microscopy is anticipated to facilitate applications of photoacoustic imaging in biomedical studies that involve dynamics and clinical procedures that demand immediate diagnosis. PMID:19021408
Boulanger, Jérôme; Gueudry, Charles; Münch, Daniel; Cinquin, Bertrand; Paul-Gilloteaux, Perrine; Bardin, Sabine; Guérin, Christophe; Senger, Fabrice; Blanchoin, Laurent; Salamero, Jean
2014-01-01
Total internal reflection fluorescence microscopy (TIRFM) is the method of choice to visualize a variety of cellular processes in particular events localized near the plasma membrane of live adherent cells. This imaging technique not relying on particular fluorescent probes provides a high sectioning capability. It is, however, restricted to a single plane. We present here a method based on a versatile design enabling fast multiwavelength azimuthal averaging and incidence angles scanning to computationally reconstruct 3D images sequences. We achieve unprecedented 50-nm axial resolution over a range of 800 nm above the coverslip. We apply this imaging modality to obtain structural and dynamical information about 3D actin architectures. We also temporally decipher distinct Rab11a-dependent exocytosis events in 3D at a rate of seven stacks per second. PMID:25404337
Courau, T.; Moustafa, S.; Plagne, L.; Poncot, A.
2013-07-01
As part of its activity, EDF R and D is developing a new nuclear core simulation code named COCAGNE. This code relies on DIABOLO, a Simplified PN (SPN) method to compute the neutron flux inside the core for eigenvalue calculations. In order to assess the accuracy of SPN calculations, we have developed DOMINO, a new 3D Cartesian SN solver. The parallel implementation of DOMINO is very efficient and allows to complete an eigenvalue calculation involving around 300 x 10{sup 9} degrees of freedom within a few hours on a single shared-memory supercomputing node. This computation corresponds to a 26-group S{sub 8} 3D PWR core model used to assess the SPN accuracy. At the pin level, the maximal error for the SP{sub 5} DIABOLO fission production rate is lower than 0.2% compared to the S{sub 8} DOMINO reference for this 3D PWR core model. (authors)
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.
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.
The fast and accurate 3D-face scanning technology based on laser triangle sensors
NASA Astrophysics Data System (ADS)
Wang, Jinjiang; Chang, Tianyu; Ge, Baozhen; Tian, Qingguo; Chen, Yang; Kong, Bin
2013-08-01
A laser triangle scanning method and the structure of 3D-face measurement system were introduced. In presented system, a liner laser source was selected as an optical indicated signal in order to scanning a line one times. The CCD image sensor was used to capture image of the laser line modulated by human face. The system parameters were obtained by system calibrated calculated. The lens parameters of image part of were calibrated with machine visual image method and the triangle structure parameters were calibrated with fine wire paralleled arranged. The CCD image part and line laser indicator were set with a linear motor carry which can achieve the line laser scanning form top of the head to neck. For the nose is ledge part and the eyes are sunk part, one CCD image sensor can not obtain the completed image of laser line. In this system, two CCD image sensors were set symmetric at two sides of the laser indicator. In fact, this structure includes two laser triangle measure units. Another novel design is there laser indicators were arranged in order to reduce the scanning time for it is difficult for human to keep static for longer time. The 3D data were calculated after scanning. And further data processing include 3D coordinate refine, mesh calculate and surface show. Experiments show that this system has simply structure, high scanning speed and accurate. The scanning range covers the whole head of adult, the typical resolution is 0.5mm.
Fast error simulation of optical 3D measurements at translucent objects
NASA Astrophysics Data System (ADS)
Lutzke, P.; Kühmstedt, P.; Notni, G.
2012-09-01
The scan results of optical 3D measurements at translucent objects deviate from the real objects surface. This error is caused by the fact that light is scattered in the objects volume and is not exclusively reflected at its surface. A few approaches were made to separate the surface reflected light from the volume scattered. For smooth objects the surface reflected light is dominantly concentrated in specular direction and could only be observed from a point in this direction. Thus the separation either leads to measurement results only creating data for near specular directions or provides data from not well separated areas. To ensure the flexibility and precision of optical 3D measurement systems for translucent materials it is necessary to enhance the understanding of the error forming process. For this purpose a technique for simulating the 3D measurement at translucent objects is presented. A simple error model is shortly outlined and extended to an efficient simulation environment based upon ordinary raytracing methods. In comparison the results of a Monte-Carlo simulation are presented. Only a few material and object parameters are needed for the raytracing simulation approach. The attempt of in-system collection of these material and object specific parameters is illustrated. The main concept of developing an error-compensation method based on the simulation environment and the collected parameters is described. The complete procedure is using both, the surface reflected and the volume scattered light for further processing.
WaveQ3D: Fast and accurate acoustic transmission loss (TL) eigenrays, in littoral environments
NASA Astrophysics Data System (ADS)
Reilly, Sean M.
This study defines a new 3D Gaussian ray bundling acoustic transmission loss model in geodetic coordinates: latitude, longitude, and altitude. This approach is designed to lower the computation burden of computing accurate environmental effects in sonar training application by eliminating the need to transform the ocean environment into a collection of Nx2D Cartesian radials. This approach also improves model accuracy by incorporating real world 3D effects, like horizontal refraction, into the model. This study starts with derivations for a 3D variant of Gaussian ray bundles in this coordinate system. To verify the accuracy of this approach, acoustic propagation predictions of transmission loss, time of arrival, and propagation direction are compared to analytic solutions and other models. To validate the model's ability to predict real world phenomena, predictions of transmission loss and propagation direction are compared to at-sea measurements, in an environment where strong horizontal refraction effect have been observed. This model has been integrated into U.S. Navy active sonar training system applications, where testing has demonstrated its ability to improve transmission loss calculation speed without sacrificing accuracy.
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.
Fast imaging of laboratory core floods using 3D compressed sensing RARE MRI.
Ramskill, N P; Bush, I; Sederman, A J; Mantle, M D; Benning, M; Anger, B C; Appel, M; Gladden, L F
2016-09-01
Three-dimensional (3D) imaging of the fluid distributions within the rock is essential to enable the unambiguous interpretation of core flooding data. Magnetic resonance imaging (MRI) has been widely used to image fluid saturation in rock cores; however, conventional acquisition strategies are typically too slow to capture the dynamic nature of the displacement processes that are of interest. Using Compressed Sensing (CS), it is possible to reconstruct a near-perfect image from significantly fewer measurements than was previously thought necessary, and this can result in a significant reduction in the image acquisition times. In the present study, a method using the Rapid Acquisition with Relaxation Enhancement (RARE) pulse sequence with CS to provide 3D images of the fluid saturation in rock core samples during laboratory core floods is demonstrated. An objective method using image quality metrics for the determination of the most suitable regularisation functional to be used in the CS reconstructions is reported. It is shown that for the present application, Total Variation outperforms the Haar and Daubechies3 wavelet families in terms of the agreement of their respective CS reconstructions with a fully-sampled reference image. Using the CS-RARE approach, 3D images of the fluid saturation in the rock core have been acquired in 16min. The CS-RARE technique has been applied to image the residual water saturation in the rock during a water-water displacement core flood. With a flow rate corresponding to an interstitial velocity of vi=1.89±0.03ftday(-1), 0.1 pore volumes were injected over the course of each image acquisition, a four-fold reduction when compared to a fully-sampled RARE acquisition. Finally, the 3D CS-RARE technique has been used to image the drainage of dodecane into the water-saturated rock in which the dynamics of the coalescence of discrete clusters of the non-wetting phase are clearly observed. The enhancement in the temporal resolution that has
Fast imaging of laboratory core floods using 3D compressed sensing RARE MRI
NASA Astrophysics Data System (ADS)
Ramskill, N. P.; Bush, I.; Sederman, A. J.; Mantle, M. D.; Benning, M.; Anger, B. C.; Appel, M.; Gladden, L. F.
2016-09-01
Three-dimensional (3D) imaging of the fluid distributions within the rock is essential to enable the unambiguous interpretation of core flooding data. Magnetic resonance imaging (MRI) has been widely used to image fluid saturation in rock cores; however, conventional acquisition strategies are typically too slow to capture the dynamic nature of the displacement processes that are of interest. Using Compressed Sensing (CS), it is possible to reconstruct a near-perfect image from significantly fewer measurements than was previously thought necessary, and this can result in a significant reduction in the image acquisition times. In the present study, a method using the Rapid Acquisition with Relaxation Enhancement (RARE) pulse sequence with CS to provide 3D images of the fluid saturation in rock core samples during laboratory core floods is demonstrated. An objective method using image quality metrics for the determination of the most suitable regularisation functional to be used in the CS reconstructions is reported. It is shown that for the present application, Total Variation outperforms the Haar and Daubechies3 wavelet families in terms of the agreement of their respective CS reconstructions with a fully-sampled reference image. Using the CS-RARE approach, 3D images of the fluid saturation in the rock core have been acquired in 16 min. The CS-RARE technique has been applied to image the residual water saturation in the rock during a water-water displacement core flood. With a flow rate corresponding to an interstitial velocity of vi = 1.89 ± 0.03 ft day-1, 0.1 pore volumes were injected over the course of each image acquisition, a four-fold reduction when compared to a fully-sampled RARE acquisition. Finally, the 3D CS-RARE technique has been used to image the drainage of dodecane into the water-saturated rock in which the dynamics of the coalescence of discrete clusters of the non-wetting phase are clearly observed. The enhancement in the temporal resolution
Fast quantitative susceptibility mapping using 3D EPI and total generalized variation.
Langkammer, Christian; Bredies, Kristian; Poser, Benedikt A; Barth, Markus; Reishofer, Gernot; Fan, Audrey Peiwen; Bilgic, Berkin; Fazekas, Franz; Mainero, Caterina; Ropele, Stefan
2015-05-01
Quantitative susceptibility mapping (QSM) allows new insights into tissue composition and organization by assessing its magnetic property. Previous QSM studies have already demonstrated that magnetic susceptibility is highly sensitive to myelin density and fiber orientation as well as to para- and diamagnetic trace elements. Image resolution in QSM with current approaches is limited by the long acquisition time of 3D scans and the need for high signal to noise ratio (SNR) to solve the dipole inversion problem. We here propose a new total-generalized-variation (TGV) based method for QSM reconstruction, which incorporates individual steps of phase unwrapping, background field removal and dipole inversion in a single iteration, thus yielding a robust solution to the reconstruction problem. This approach has beneficial characteristics for low SNR data, allowing for phase data to be rapidly acquired with a 3D echo planar imaging (EPI) sequence. The proposed method was evaluated with a numerical phantom and in vivo at 3 and 7 T. Compared to total variation (TV), TGV-QSM enforced higher order smoothness which yielded solutions closer to the ground truth and prevented stair-casing artifacts. The acquisition time for images with 1mm isotropic resolution and whole brain coverage was 10s on a clinical 3 Tesla scanner. In conclusion, 3D EPI acquisition combined with single-step TGV reconstruction yields reliable QSM images of the entire brain with 1mm isotropic resolution in seconds. The short acquisition time combined with the robust reconstruction may enable new QSM applications in less compliant populations, clinical susceptibility tensor imaging, and functional resting state examinations. PMID:25731991
NASA Astrophysics Data System (ADS)
Chakraborty, Bidisha; Heyde, Brecht; Alessandrini, Martino; D'hooge, Jan
2016-04-01
Image registration techniques using free-form deformation models have shown promising results for 3D myocardial strain estimation from ultrasound. However, the use of this technique has mostly been limited to research institutes due to the high computational demand, which is primarily due to the computational load of the regularization term ensuring spatially smooth cardiac strain estimates. Indeed, this term typically requires evaluating derivatives of the transformation field numerically in each voxel of the image during every iteration of the optimization process. In this paper, we replace this time-consuming step with a closed-form solution directly associated with the transformation field resulting in a speed up factor of ~10-60,000, for a typical 3D B-mode image of 2503 and 5003 voxels, depending upon the size and the parametrization of the transformation field. The performance of the numeric and the analytic solutions was contrasted by computing tracking and strain accuracy on two realistic synthetic 3D cardiac ultrasound sequences, mimicking two ischemic motion patterns. Mean and standard deviation of the displacement errors over the cardiac cycle for the numeric and analytic solutions were 0.68+/-0.40 mm and 0.75+/-0.43 mm respectively. Correlations for the radial, longitudinal and circumferential strain components at end-systole were 0.89, 0.83 and 0.95 versus 0.90, 0.88 and 0.92 for the numeric and analytic regularization respectively. The analytic solution matched the performance of the numeric solution as no statistically significant differences (p>0.05) were found when expressed in terms of bias or limits-of-agreement.
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)
Cheng, Cheng; Bodin, Thomas; Allen, Richard M.
2016-02-01
We present a novel 3D pre-stack Kirchhoff depth migration (PKDM) method for teleseismic receiver functions. The proposed algorithm considers the effects of diffraction, scattering, and travel time alteration caused by 3D volumetric heterogeneities. It is therefore particularly useful for imaging complex 3D structures such as dipping discontinuities, which is hard to accomplish with traditional methods. The scheme is based on the acoustic wave migration principle, where at each time step of the receiver function, the energy is migrated back to the ensemble of potential conversion points in the image, given a smooth 3D reference model. Travel times for P and S waves are computed with an efficient Eikonal solver, the Fast Marching Method. We also consider elastic scattering patterns, where the amplitude of converted S waves depends on the angle between the incident P wave, and the scattered S wave. Synthetic experiments demonstrate the validity of the method for a variety of dipping angle discontinuities. Comparison with the widely used Common Conversion Point (CCP) stacking method reveals that our migration shows considerable improvement. For example, the effect of multiple reflections that usually produce apparent discontinuities is avoided. The proposed approach is practical, computationally efficient, and is therefore a potentially powerful alternative to standard CCP methods for imaging large-scale continental structure under dense networks.
Fast and efficient particle reconstruction on a 3D grid using sparsity
NASA Astrophysics Data System (ADS)
Cornic, P.; Champagnat, F.; Cheminet, A.; Leclaire, B.; Le Besnerais, G.
2015-03-01
We propose an approach for efficient localization and intensity reconstruction of particles on a 3D grid based on sparsity principles. The computational complexity of the method is limited by using the particle volume reconstruction paradigm (Champagnat et al. in Meas Sci Technol 25, 2014) and a reduction in the problem dimension. Tests on synthetic and experimental data show that the proposed method leads to more efficient detections and to reconstructions of higher quality than classical tomoPIV approaches on a large range of seeding densities, up to ppp ≈ 0.12.
Atlas Toolkit: Fast registration of 3D morphological datasets in the absence of landmarks
Grocott, Timothy; Thomas, Paul; Münsterberg, Andrea E.
2016-01-01
Image registration is a gateway technology for Developmental Systems Biology, enabling computational analysis of related datasets within a shared coordinate system. Many registration tools rely on landmarks to ensure that datasets are correctly aligned; yet suitable landmarks are not present in many datasets. Atlas Toolkit is a Fiji/ImageJ plugin collection offering elastic group-wise registration of 3D morphological datasets, guided by segmentation of the interesting morphology. We demonstrate the method by combinatorial mapping of cell signalling events in the developing eyes of chick embryos, and use the integrated datasets to predictively enumerate Gene Regulatory Network states. PMID:26864723
3D Multislab, Multishot Acquisition for Fast, Whole-Brain MR Elastography with High SNR Efficiency
Johnson, Curtis L; Holtrop, Joseph L; McGarry, Matthew DJ; Weaver, John B; Paulsen, Keith D; Georgiadis, John G; Sutton, Bradley P
2014-01-01
Purpose To develop an acquisition scheme for generating magnetic resonance elastography (MRE) displacement data with whole-brain coverage, high spatial resolution, and adequate signal-to-noise ratio (SNR) in a short scan time. Theory and Methods A 3D multislab, multishot acquisition for whole-brain MRE with 2.0 mm isotropic spatial resolution is proposed. The multislab approach allowed for the use of short repetition time to achieve very high SNR efficiency. High SNR efficiency allowed for a reduced acquisition time of only six minutes while the minimum SNR needed for inversion was maintained. Results The mechanical property maps estimated from whole-brain displacement data with nonlinear inversion (NLI) demonstrated excellent agreement with neuroanatomical features, including the cerebellum and brainstem. A comparison with an equivalent 2D acquisition illustrated the improvement in SNR efficiency of the 3D multislab acquisition. The flexibility afforded by the high SNR efficiency allowed for higher resolution with a 1.6 mm isotropic voxel size, which generated higher estimates of brainstem stiffness compared with the 2.0 mm isotropic acquisition. Conclusions The acquisition presented allows for the capture of whole-brain MRE displacement data in a short scan time, and may be used to generate local mechanical property estimates of neuroanatomical features throughout the brain. PMID:24347237
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
Contrast-Based 3D/2D Registration of the Left Atrium: Fast versus Consistent
Kowalewski, Christopher; Kurzidim, Klaus; Strobel, Norbert; Hornegger, Joachim
2016-01-01
For augmented fluoroscopy during cardiac ablation, a preoperatively acquired 3D model of a patient's left atrium (LA) can be registered to X-ray images recorded during a contrast agent (CA) injection. An automatic registration method that works also for small amounts of CA is desired. We propose two similarity measures: The first focuses on edges of the patient anatomy. The second computes a contrast agent distribution estimate (CADE) inside the 3D model and rates its consistency with the CA as seen in biplane fluoroscopic images. Moreover, temporal filtering on the obtained registration results of a sequence is applied using a Markov chain framework. Evaluation was performed on 11 well-contrasted clinical angiographic sequences and 10 additional sequences with less CA. For well-contrasted sequences, the error for all 73 frames was 7.9 ± 6.3 mm and it dropped to 4.6 ± 4.0 mm when registering to an automatically selected, well enhanced frame in each sequence. Temporal filtering reduced the error for all frames from 7.9 ± 6.3 mm to 5.7 ± 4.6 mm. The error was typically higher if less CA was used. A combination of both similarity measures outperforms a previously proposed similarity measure. The mean accuracy for well contrasted sequences is in the range of other proposed manual registration methods. PMID:27051412
Radial subsampling for fast cost function computation in intensity-based 3D image registration
NASA Astrophysics Data System (ADS)
Boettger, Thomas; Wolf, Ivo; Meinzer, Hans-Peter; Celi, Juan Carlos
2007-03-01
Image registration is always a trade-off between accuracy and speed. Looking towards clinical scenarios the time for bringing two or more images into registration should be around a few seconds only. We present a new scheme for subsampling 3D-image data to allow for efficient computation of cost functions in intensity-based image registration. Starting from an arbitrary center point voxels are sampled along scan lines which do radially extend from the center point. We analyzed the characteristics of different cost functions computed on the sub-sampled data and compared them to known cost functions with respect to local optima. Results show the cost functions are smooth and give high peaks at the expected optima. Furthermore we investigated capture range of cost functions computed under the new subsampling scheme. Capture range was remarkably better for the new scheme compared to metrics using all voxels or different subsampling schemes and high registration accuracy was achieved as well. The most important result is the improvement in terms of speed making this scheme very interesting for clinical scenarios. We conclude using the new subsampling scheme intensity-based 3D image registration can be performed much faster than using other approaches while maintaining high accuracy. A variety of different extensions of the new approach is conceivable, e.g. non-regular distribution of the scan lines or not to let the scan lines start from a center point only, but from the surface of an organ model for example.
Contrast-Based 3D/2D Registration of the Left Atrium: Fast versus Consistent.
Hoffmann, Matthias; Kowalewski, Christopher; Maier, Andreas; Kurzidim, Klaus; Strobel, Norbert; Hornegger, Joachim
2016-01-01
For augmented fluoroscopy during cardiac ablation, a preoperatively acquired 3D model of a patient's left atrium (LA) can be registered to X-ray images recorded during a contrast agent (CA) injection. An automatic registration method that works also for small amounts of CA is desired. We propose two similarity measures: The first focuses on edges of the patient anatomy. The second computes a contrast agent distribution estimate (CADE) inside the 3D model and rates its consistency with the CA as seen in biplane fluoroscopic images. Moreover, temporal filtering on the obtained registration results of a sequence is applied using a Markov chain framework. Evaluation was performed on 11 well-contrasted clinical angiographic sequences and 10 additional sequences with less CA. For well-contrasted sequences, the error for all 73 frames was 7.9 ± 6.3 mm and it dropped to 4.6 ± 4.0 mm when registering to an automatically selected, well enhanced frame in each sequence. Temporal filtering reduced the error for all frames from 7.9 ± 6.3 mm to 5.7 ± 4.6 mm. The error was typically higher if less CA was used. A combination of both similarity measures outperforms a previously proposed similarity measure. The mean accuracy for well contrasted sequences is in the range of other proposed manual registration methods. PMID:27051412
Ruh, Dominic; Tränkle, Benjamin; Rohrbach, Alexander
2011-10-24
Multi-dimensional, correlated particle tracking is a key technology to reveal dynamic processes in living and synthetic soft matter systems. In this paper we present a new method for tracking micron-sized beads in parallel and in all three dimensions - faster and more precise than existing techniques. Using an acousto-optic deflector and two quadrant-photo-diodes, we can track numerous optically trapped beads at up to tens of kHz with a precision of a few nanometers by back-focal plane interferometry. By time-multiplexing the laser focus, we can calibrate individually all traps and all tracking signals in a few seconds and in 3D. We show 3D histograms and calibration constants for nine beads in a quadratic arrangement, although trapping and tracking is easily possible for more beads also in arbitrary 2D arrangements. As an application, we investigate the hydrodynamic coupling and diffusion anomalies of spheres trapped in a 3 × 3 arrangement. PMID:22109012
NASA Astrophysics Data System (ADS)
Huang, Qun-xing; Yan, Fei Wang Jian-hua; Chi, Yong
2013-04-01
A new approach to inverse radiation analysis is presented for non-intrusive 3-D flame temperature reconstruction using flame emission images from four CCD camera detectors installed on the furnace wall. The scattering from participating medium in the flame was considered by combining the discrete radiative transfer method with the discrete ordinate method. A modified minimum residual algorithm was employed to calculate the least squares solution of the ill-conditioned inverse problem. A numerical test problem simulating real temperature measurements in an industrial furnace was used to assess the performance of the proposed method. These assessments indicate that this method is capable of reconstructing 3-D temperature distributions fast and accurately, even with noisy flame emission data. Such a capability has potential in real-time temperature measurement for combustion optimization and pollution emission control.
Jia, Jia; Liu, Juan; Jin, Guofan; Wang, Yongtian
2014-09-20
Occlusion culling is an important process that produces correct depth cues for observers in holographic displays, whereas current methods suffer from occlusion errors or high computational loads. We propose a fast and effective method for occlusion culling based on multiple light-point sampling planes and an inverse orthographic projection technique. Multiple light-point sampling planes are employed to remove the hidden surfaces for each direction of the view of the three-dimensional (3D) scene by forward orthographic projection, and the inverse orthographic projection technique is used to determine the effective sampling points of the 3D scene. A numerical simulation and an optical experiment are performed. The results show that this approach can realize accurate occlusion effects, smooth motion parallax, and continuous depth using low angular sampling without any extra computation costs. PMID:25322109
Very fast road database verification using textured 3D city models obtained from airborne imagery
NASA Astrophysics Data System (ADS)
Bulatov, Dimitri; Ziems, Marcel; Rottensteiner, Franz; Pohl, Melanie
2014-10-01
Road databases are known to be an important part of any geodata infrastructure, e.g. as the basis for urban planning or emergency services. Updating road databases for crisis events must be performed quickly and with the highest possible degree of automation. We present a semi-automatic algorithm for road verification using textured 3D city models, starting from aerial or even UAV-images. This algorithm contains two processes, which exchange input and output, but basically run independently from each other. These processes are textured urban terrain reconstruction and road verification. The first process contains a dense photogrammetric reconstruction of 3D geometry of the scene using depth maps. The second process is our core procedure, since it contains various methods for road verification. Each method represents a unique road model and a specific strategy, and thus is able to deal with a specific type of roads. Each method is designed to provide two probability distributions, where the first describes the state of a road object (correct, incorrect), and the second describes the state of its underlying road model (applicable, not applicable). Based on the Dempster-Shafer Theory, both distributions are mapped to a single distribution that refers to three states: correct, incorrect, and unknown. With respect to the interaction of both processes, the normalized elevation map and the digital orthophoto generated during 3D reconstruction are the necessary input - together with initial road database entries - for the road verification process. If the entries of the database are too obsolete or not available at all, sensor data evaluation enables classification of the road pixels of the elevation map followed by road map extraction by means of vectorization and filtering of the geometrically and topologically inconsistent objects. Depending on the time issue and availability of a geo-database for buildings, the urban terrain reconstruction procedure has semantic models
Mapping of the spontaneous deletion in the Ap3d1 gene of mocha mice: fast and reliable genotyping
Drasbek, Kim Ryun; Holm, Mai Marie; Delenclos, Marion; Jensen, Kimmo
2008-01-01
Background The mocha mouse carries a spontaneous deletion in the Ap3d1 gene, encoding the delta 1 subunit of the adaptor related protein complex 3, (Ap3d1), and subsequently lack the expression of functional AP-3. This leads to a deficiency in vesicle transport and storage, which affects neurotransmitter vesicle turnover and release in the central nervous system. Since the genomic sequence of the Ap3d1 gene of mocha mouse is not known, precise mapping of the deletion as well as reliable genotyping protocols are lacking. Findings We sequenced the Ap3d1 gene (HGNC GeneID: 8943) around the deletion site in the mocha mouse and revealed a 10639 bp deletion covering exon 2 to 6. Subsequently, new PCR primers were designed yielding a reliable genotyping protocol of both newborn and adult tissue. To examine the genotypes further, hippocampal neurons were cultured from mocha and control mice. Patch-clamp recordings showed that mocha neurons had a higher input resistance, and that autaptic EPSC in mocha cultures depressed faster and stronger as compared with control cultures. Conclusion Our study reports the sequence of the deleted part of the Ap3d1 gene in mocha mice, as well as a reliable PCR-based genotyping protocol. We cultured hippocampal neurons from control and mocha mice, and found a difference in input resistance of the neurons, and in the synaptic short-term plasticity of glutamatergic autapses showing a larger synaptic depression than controls. The described procedures may be useful for the future utilization of the mocha mouse as a model of defective vesicle biogenesis. Importantly, as genotyping by eye color is complicated in newborn mice, the designed protocol is so fast and reliable that newborn mice could rapidly be genotyped and hippocampal neurons dissociated and cultured, which is normally best done at P0-P2. PMID:19032734
Fast and Robust Sixth Order Multigrid Computation for 3D Convection Diffusion Equation.
Wang, Yin; Zhang, Jun
2010-10-15
We present a sixth order explicit compact finite difference scheme to solve the three dimensional (3D) convection diffusion equation. We first use multiscale multigrid method to solve the linear systems arising from a 19-point fourth order discretization scheme to compute the fourth order solutions on both the coarse grid and the fine grid. Then an operator based interpolation scheme combined with an extrapolation technique is used to approximate the sixth order accurate solution on the fine grid. Since the multigrid method using a standard point relaxation smoother may fail to achieve the optimal grid independent convergence rate for solving convection diffusion equation with a high Reynolds number, we implement the plane relaxation smoother in the multigrid solver to achieve better grid independency. Supporting numerical results are presented to demonstrate the efficiency and accuracy of the sixth order compact scheme (SOC), compared with the previously published fourth order compact scheme (FOC). PMID:21151737
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.
NASA Technical Reports Server (NTRS)
Hung, Stephen H. Y.
1989-01-01
A fast 3-D object recognition algorithm that can be used as a quick-look subsystem to the vision system for the Special-Purpose Dexterous Manipulator (SPDM) is described. Global features that can be easily computed from range data are used to characterize the images of a viewer-centered model of an object. This algorithm will speed up the processing by eliminating the low level processing whenever possible. It may identify the object, reject a set of bad data in the early stage, or create a better environment for a more powerful algorithm to carry the work further.
A fast technique applied to the analysis of Resistive Wall Modes with 3D conducting structures
Rubinacci, Guglielmo Liu, Yueqiang
2009-03-20
This paper illustrates the development of a 'fast' technique for the analysis of Resistive Wall Modes (RWMs) in fusion devices with three-dimensional conducting structures, by means of the recently developed CarMa code. Thanks to its peculiar features, the computational cost scales almost linearly with the number of discrete unknowns. Some large scale problems are solved in configurations of interest for the International Thermonuclear Experimental Reactor (ITER)
A fast nested dissection solver for Cartesian 3D elliptic problems using hierarchical matrices
NASA Astrophysics Data System (ADS)
Schmitz, Phillip G.; Ying, Lexing
2014-02-01
We present a fast algorithm for solutions to linear systems arising from three dimensional elliptic problems on a regular Cartesian mesh. We follow the approach of Schmitz and Ying (2012) on combining the nested dissection matrix factorization method with hierarchical matrices in two dimensions and extend it to the three dimensional case. A theoretical linear time complexity is derived and a more practical variant with slightly worse scaling is demonstrated.
Poulin, Eric; Racine, Emmanuel; Beaulieu, Luc; Binnekamp, Dirk
2015-03-15
Purpose: In high dose rate brachytherapy (HDR-B), current catheter reconstruction protocols are relatively slow and error prone. The purpose of this technical note is to evaluate the accuracy and the robustness of an electromagnetic (EM) tracking system for automated and real-time catheter reconstruction. Methods: For this preclinical study, a total of ten catheters were inserted in gelatin phantoms with different trajectories. Catheters were reconstructed using a 18G biopsy needle, used as an EM stylet and equipped with a miniaturized sensor, and the second generation Aurora{sup ®} Planar Field Generator from Northern Digital Inc. The Aurora EM system provides position and orientation value with precisions of 0.7 mm and 0.2°, respectively. Phantoms were also scanned using a μCT (GE Healthcare) and Philips Big Bore clinical computed tomography (CT) system with a spatial resolution of 89 μm and 2 mm, respectively. Reconstructions using the EM stylet were compared to μCT and CT. To assess the robustness of the EM reconstruction, five catheters were reconstructed twice and compared. Results: Reconstruction time for one catheter was 10 s, leading to a total reconstruction time inferior to 3 min for a typical 17-catheter implant. When compared to the μCT, the mean EM tip identification error was 0.69 ± 0.29 mm while the CT error was 1.08 ± 0.67 mm. The mean 3D distance error was found to be 0.66 ± 0.33 mm and 1.08 ± 0.72 mm for the EM and CT, respectively. EM 3D catheter trajectories were found to be more accurate. A maximum difference of less than 0.6 mm was found between successive EM reconstructions. Conclusions: The EM reconstruction was found to be more accurate and precise than the conventional methods used for catheter reconstruction in HDR-B. This approach can be applied to any type of catheters and applicators.
Development of a 3D numerical methodology for fast prediction of gun blast induced loading
NASA Astrophysics Data System (ADS)
Costa, E.; Lagasco, F.
2014-05-01
In this paper, the development of a methodology based on semi-empirical models from the literature to carry out 3D prediction of pressure loading on surfaces adjacent to a weapon system during firing is presented. This loading is consequent to the impact of the blast wave generated by the projectile exiting the muzzle bore. When exceeding a pressure threshold level, loading is potentially capable to induce unwanted damage to nearby hard structures as well as frangible panels or electronic equipment. The implemented model shows the ability to quickly predict the distribution of the blast wave parameters over three-dimensional complex geometry surfaces when the weapon design and emplacement data as well as propellant and projectile characteristics are available. Considering these capabilities, the use of the proposed methodology is envisaged as desirable in the preliminary design phase of the combat system to predict adverse effects and then enable to identify the most appropriate countermeasures. By providing a preliminary but sensitive estimate of the operative environmental loading, this numerical means represents a good alternative to more powerful, but time consuming advanced computational fluid dynamics tools, which use can, thus, be limited to the final phase of the design.
Fast simulation of solid tumors thermal ablation treatments with a 3D reaction diffusion model.
Bertaccini, Daniele; Calvetti, Daniela
2007-08-01
An efficient computational method for near real-time simulation of thermal ablation of tumors via radio frequencies is proposed. Model simulations of the temperature field in a 3D portion of tissue containing the tumoral mass for different patterns of source heating can be used to design the ablation procedure. The availability of a very efficient computational scheme makes it possible to update the predicted outcome of the procedure in real time. In the algorithms proposed here a discretization in space of the governing equations is followed by an adaptive time integration based on implicit multistep formulas. A modification of the ode15s MATLAB function which uses Krylov space iterative methods for the solution of the linear systems arising at each integration step makes it possible to perform the simulations on standard desktop for much finer grids than using the built-in ode15s. The proposed algorithm can be applied to a wide class of nonlinear parabolic differential equations. PMID:17173888
FAST SIMULATION OF SOLID TUMORS THERMAL ABLATION TREATMENTS WITH A 3D REACTION DIFFUSION MODEL *
BERTACCINI, DANIELE; CALVETTI, DANIELA
2007-01-01
An efficient computational method for near real-time simulation of thermal ablation of tumors via radio frequencies is proposed. Model simulations of the temperature field in a 3D portion of tissue containing the tumoral mass for different patterns of source heating can be used to design the ablation procedure. The availability of a very efficient computational scheme makes it possible update the predicted outcome of the procedure in real time. In the algorithms proposed here a discretization in space of the governing equations is followed by an adaptive time integration based on implicit multistep formulas. A modification of the ode15s MATLAB function which uses Krylov space iterative methods for the solution of for the linear systems arising at each integration step makes it possible to perform the simulations on standard desktop for much finer grids than using the built-in ode15s. The proposed algorithm can be applied to a wide class of nonlinear parabolic differential equations. PMID:17173888
Sakhalkar, H S; Oldham, M
2008-01-01
This study introduces a charge coupled device (CCD) area detector based optical-computed tomography (optical-CT) scanner for comprehensive verification of radiation dose distributions recorded in nonscattering radiochromic dosimeters. Defining characteristics include: (i) a very fast scanning time of approximately 5 min to acquire a complete three-dimensional (3D) dataset, (ii) improved image formation through the use of custom telecentric optics, which ensures accurate projection images and minimizes artifacts from scattered and stray-light sources, and (iii) high resolution (potentially 50 microm) isotropic 3D dose readout. The performance of the CCD scanner for 3D dose readout was evaluated by comparison with independent 3D readout from the single laser beam OCTOPUS-scanner for the same PRESAGE dosimeters. The OCTOPUS scanner was considered the "gold standard" technique in light of prior studies demonstrating its accuracy. Additional comparisons were made against calculated dose distributions from the ECLIPSE treatment-planning system. Dose readout for the following treatments were investigated: (i) a single rectangular beam irradiation to investigate small field and very steep dose gradient dosimetry away from edge effects, (ii) a 2-field open beam parallel-opposed irradiation to investigate dosimetry along steep dose gradients, and (iii) a 7-field intensity modulated radiation therapy (IMRT) irradiation to investigate dosimetry for complex treatment delivery involving modulation of fluence and for dosimetry along moderate dose gradients. Dose profiles, dose-difference plots, and gamma maps were employed to evaluate quantitative estimates of agreement between independently measured and calculated dose distributions. Results indicated that dose readout from the CCD scanner was in agreement with independent gold-standard readout from the OCTOPUS-scanner as well as the calculated ECLIPSE dose distribution for all treatments, except in regions within a few
NCC-RANSAC: A Fast Plane Extraction Method for 3-D Range Data Segmentation
Qian, Xiangfei; Ye, Cang
2015-01-01
This paper presents a new plane extraction (PE) method based on the random sample consensus (RANSAC) approach. The generic RANSAC-based PE algorithm may over-extract a plane, and it may fail in case of a multistep scene where the RANSAC procedure results in multiple inlier patches that form a slant plane straddling the steps. The CC-RANSAC PE algorithm successfully overcomes the latter limitation if the inlier patches are separate. However, it fails if the inlier patches are connected. A typical scenario is a stairway with a stair wall where the RANSAC plane-fitting procedure results in inliers patches in the tread, riser, and stair wall planes. They connect together and form a plane. The proposed method, called normal-coherence CC-RANSAC (NCC-RANSAC), performs a normal coherence check to all data points of the inlier patches and removes the data points whose normal directions are contradictory to that of the fitted plane. This process results in separate inlier patches, each of which is treated as a candidate plane. A recursive plane clustering process is then executed to grow each of the candidate planes until all planes are extracted in their entireties. The RANSAC plane-fitting and the recursive plane clustering processes are repeated until no more planes are found. A probabilistic model is introduced to predict the success probability of the NCC-RANSAC algorithm and validated with real data of a 3-D time-of-flight camera–SwissRanger SR4000. Experimental results demonstrate that the proposed method extracts more accurate planes with less computational time than the existing RANSAC-based methods. PMID:24771605
NASA Astrophysics Data System (ADS)
Zhang, Lisha
We present fast and robust numerical algorithms for 3-D scattering from perfectly electrical conducting (PEC) and dielectric random rough surfaces in microwave remote sensing. The Coifman wavelets or Coiflets are employed to implement Galerkin's procedure in the method of moments (MoM). Due to the high-precision one-point quadrature, the Coiflets yield fast evaluations of the most off-diagonal entries, reducing the matrix fill effort from O(N2) to O( N). The orthogonality and Riesz basis of the Coiflets generate well conditioned impedance matrix, with rapid convergence for the conjugate gradient solver. The resulting impedance matrix is further sparsified by the matrix-formed standard fast wavelet transform (SFWT). By properly selecting multiresolution levels of the total transformation matrix, the solution precision can be enhanced while matrix sparsity and memory consumption have not been noticeably sacrificed. The unified fast scattering algorithm for dielectric random rough surfaces can asymptotically reduce to the PEC case when the loss tangent grows extremely large. Numerical results demonstrate that the reduced PEC model does not suffer from ill-posed problems. Compared with previous publications and laboratory measurements, good agreement is observed.
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.
Fast particles identification in programmable form at level-0 trigger by means of the 3D-Flow system
Crosetto, Dario B.
1998-10-30
The 3D-Flow Processor system is a new, technology-independent concept in very fast, real-time system architectures. Based on either an FPGA or an ASIC implementation, it can address, in a fully programmable manner, applications where commercially available processors would fail because of throughput requirements. Possible applications include filtering-algorithms (pattern recognition) from the input of multiple sensors, as well as moving any input validated by these filtering-algorithms to a single output channel. Both operations can easily be implemented on a 3D-Flow system to achieve a real-time processing system with a very short lag time. This system can be built either with off-the-shelf FPGAs or, for higher data rates, with CMOS chips containing 4 to 16 processors each. The basic building block of the system, a 3D-Flow processor, has been successfully designed in VHDL code written in ''Generic HDL'' (mostly made of reusable blocks that are synthesizable in different technologies, or FPGAs), to produce a netlist for a four-processor ASIC featuring 0.35 micron CBA (Ceil Base Array) technology at 3.3 Volts, 884 mW power dissipation at 60 MHz and 63.75 mm sq. die size. The same VHDL code has been targeted to three FPGA manufacturers (Altera EPF10K250A, ORCA-Lucent Technologies 0R3T165 and Xilinx XCV1000). A complete set of software tools, the 3D-Flow System Manager, equally applicable to ASIC or FPGA implementations, has been produced to provide full system simulation, application development, real-time monitoring, and run-time fault recovery. Today's technology can accommodate 16 processors per chip in a medium size die, at a cost per processor of less than $5 based on the current silicon die/size technology cost.
Fast spatial beam shaping by acousto-optic diffraction for 3D non-linear microscopy.
Akemann, Walther; Léger, Jean-François; Ventalon, Cathie; Mathieu, Benjamin; Dieudonné, Stéphane; Bourdieu, Laurent
2015-11-01
Acousto-optic deflection (AOD) devices offer unprecedented fast control of the entire spatial structure of light beams, most notably their phase. AOD light modulation of ultra-short laser pulses, however, is not straightforward to implement because of intrinsic chromatic dispersion and non-stationarity of acousto-optic diffraction. While schemes exist to compensate chromatic dispersion, non-stationarity remains an obstacle. In this work we demonstrate an efficient AOD light modulator for stable phase modulation using time-locked generation of frequency-modulated acoustic waves at the full repetition rate of a high power laser pulse amplifier of 80 kHz. We establish the non-local relationship between the optical phase and the generating acoustic frequency function and verify the system for temporal stability, phase accuracy and generation of non-linear two-dimensional phase functions. PMID:26561090
LBM-EP: Lattice-Boltzmann method for fast cardiac electrophysiology simulation from 3D images.
Rapaka, S; Mansi, T; Georgescu, B; Pop, M; Wright, G A; Kamen, A; Comaniciu, Dorin
2012-01-01
Current treatments of heart rhythm troubles require careful planning and guidance for optimal outcomes. Computational models of cardiac electrophysiology are being proposed for therapy planning but current approaches are either too simplified or too computationally intensive for patient-specific simulations in clinical practice. This paper presents a novel approach, LBM-EP, to solve any type of mono-domain cardiac electrophysiology models at near real-time that is especially tailored for patient-specific simulations. The domain is discretized on a Cartesian grid with a level-set representation of patient's heart geometry, previously estimated from images automatically. The cell model is calculated node-wise, while the transmembrane potential is diffused using Lattice-Boltzmann method within the domain defined by the level-set. Experiments on synthetic cases, on a data set from CESC'10 and on one patient with myocardium scar showed that LBM-EP provides results comparable to an FEM implementation, while being 10 - 45 times faster. Fast, accurate, scalable and requiring no specific meshing, LBM-EP paves the way to efficient and detailed models of cardiac electrophysiology for therapy planning. PMID:23286029
Field calibration of binocular stereo vision based on fast reconstruction of 3D control field
NASA Astrophysics Data System (ADS)
Zhang, Haijun; Liu, Changjie; Fu, Luhua; Guo, Yin
2015-08-01
Construction of high-speed railway in China has entered a period of rapid growth. To accurately and quickly obtain the dynamic envelope curve of high-speed vehicle is an important guarantee for safe driving. The measuring system is based on binocular stereo vision. Considering the difficulties in field calibration such as environmental changes and time limits, carried out a field calibration method based on fast reconstruction of three-dimensional control field. With the rapid assembly of pre-calibrated three-dimensional control field, whose coordinate accuracy is guaranteed by manufacture accuracy and calibrated by V-STARS, two cameras take a quick shot of it at the same time. The field calibration parameters are then solved by the method combining linear solution with nonlinear optimization. Experimental results showed that the measurement accuracy can reach up to +/- 0.5mm, and more importantly, in the premise of guaranteeing accuracy, the speed of the calibration and the portability of the devices have been improved considerably.
Comparison of fast 3D simulation and actinic inspection for EUV masks with buries defects
Clifford, C. H.; Wiraatmadja, S.; Chan, T. T.; Neureuther, A. R.; Goldberg, K. A.; Mochi, I.; Liang, T.
2009-02-23
Aerial images for isolated defects and the interactions of defects with features are compared between the Actinic Inspection Tool (AIT) at Lawrence Berkeley National Laboratory (LBNL) and the fast EUV simulation program RADICAL. Comparisons between AIT images from August 2007 and RADICAL simulations are used to extract aberrations. At this time astigmatism was the dominant aberration with a value of 0.55 waves RMS. Significant improvements in the imaging performance of the AIT were made between August 2007 and December 2008. A good match will be shown between the most recent AIT images and RADICAL simulations without aberrations. These comparisons will demonstrate that a large defect, in this case 7nm tall on the surface, is still printable even if it is centered under the absorber line. These comparisons also suggest that the minimum defect size is between 1.5nm and 0.8nm surface height because a 1.5nm defect was printable but a 0.8nm was not. Finally, the image of a buried defect near an absorber line through focus will demonstrate an inversion in the effect of the defect from a protrusion of the dark line into the space to a protrusion of the space into the line.
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.
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.
Sanz-Requena, Roberto; Moratal, David; García-Sánchez, Diego Ramón; Bodí, Vicente; Rieta, José Joaquín; Sanchis, Juan Manuel
2007-03-01
Intravascular ultrasound (IVUS) imaging is used along with X-ray coronary angiography to detect vessel pathologies. Manual analysis of IVUS images is slow and time-consuming and it is not feasible for clinical purposes. A semi-automated method is proposed to generate 3D reconstructions from IVUS video sequences, so that a fast diagnose can be easily done, quantifying plaque length and severity as well as plaque volume of the vessels under study. The methodology described in this work has four steps: a pre-processing of IVUS images, a segmentation of media-adventitia contour, a detection of intima and plaque and a 3D reconstruction of the vessel. Preprocessing is intended to remove noise from the images without blurring the edges. Segmentation of media-adventitia contour is achieved using active contours (snakes). In particular, we use the gradient vector flow (GVF) as external force for the snakes. The detection of lumen border is obtained taking into account gray-level information of the inner part of the previously detected contours. A knowledge-based approach is used to determine which level of gray corresponds statistically to the different regions of interest: intima, plaque and lumen. The catheter region is automatically discarded. An estimate of plaque type is also given. Finally, 3D reconstruction of all detected regions is made. The suitability of this methodology has been verified for the analysis and visualization of plaque length, stenosis severity, automatic detection of the most problematic regions, calculus of plaque volumes and a preliminary estimation of plaque type obtaining for automatic measures of lumen and vessel area an average error smaller than 1mm(2) (equivalent aproximately to 10% of the average measure), for calculus of plaque and lumen volume errors smaller than 0.5mm(3) (equivalent approximately to 20% of the average measure) and for plaque type estimates a mismatch of less than 8% in the analysed frames. PMID:17215103
Fast, large field-of-view, telecentric optical-CT scanning system for 3D radiochromic dosimetry
Thomas, A; Oldham, M
2010-01-01
We describe initial experiences with an in-house, fast, large field-of-view optical-CT telecentric scanner (the Duke Large field of view Optical-CT Scanner (DLOS)). The DLOS system is designed to enable telecentric optical-CT imaging of dosimeters up to 24 cm in diameter with a spatial resolution of 1 mm3, in approximately 10 minutes. These capabilities render the DLOS system a unique device at present. The system is a scaled up version of early prototypes in our lab. This scaling introduces several challenges, including the accurate measurement of a greatly increased range of light attenuation within the dosimeter, and the need to reduce even minor reflections and scattered light within the imaging chain. We present several corrections and techniques that enable accurate, low noise, 3D dosimetery with the DLOS system. PMID:21218169
Fast, Accurate and Precise Mid-Sagittal Plane Location in 3D MR Images of the Brain
NASA Astrophysics Data System (ADS)
Bergo, Felipe P. G.; Falcão, Alexandre X.; Yasuda, Clarissa L.; Ruppert, Guilherme C. S.
Extraction of the mid-sagittal plane (MSP) is a key step for brain image registration and asymmetry analysis. We present a fast MSP extraction method for 3D MR images, based on automatic segmentation of the brain and on heuristic maximization of the cerebro-spinal fluid within the MSP. The method is robust to severe anatomical asymmetries between the hemispheres, caused by surgical procedures and lesions. The method is also accurate with respect to MSP delineations done by a specialist. The method was evaluated on 64 MR images (36 pathological, 20 healthy, 8 synthetic), and it found a precise and accurate approximation of the MSP in all of them with a mean time of 60.0 seconds per image, mean angular variation within a same image (precision) of 1.26o and mean angular difference from specialist delineations (accuracy) of 1.64o.
NASA Astrophysics Data System (ADS)
Lee, Kyungmoo; Zhang, Li; Abramoff, Michael D.; Sonka, Milan
2015-03-01
Image segmentation is important for quantitative analysis of medical image data. Recently, our research group has introduced a 3-D graph search method which can simultaneously segment optimal interacting surfaces with respect to the cost function in volumetric images. Although it provides excellent segmentation accuracy, it is computationally demanding (both CPU and memory) to simultaneously segment multiple surfaces from large volumetric images. Therefore, we propose a new, fast, and memory-efficient graph search method for intraretinal layer segmentation of 3-D macular optical coherence tomograpy (OCT) scans. The key idea is to reduce the size of a graph by combining the nodes with high costs based on the multiscale approach. The new approach requires significantly less memory and achieves significantly faster processing speeds (p < 0.01) with only small segmentation differences compared to the original graph search method. This paper discusses sub-optimality of this approach and assesses trade-off relationships between decreasing processing speed and increasing segmentation differences from that of the original method as a function of employed scale of the underlying graph construction.
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
The ATLAS3D project - XIX. The hot gas content of early-type galaxies: fast versus slow rotators
NASA Astrophysics Data System (ADS)
Sarzi, Marc; Alatalo, Katherine; Blitz, Leo; Bois, Maxime; Bournaud, Frédéric; Bureau, Martin; Cappellari, Michele; Crocker, Alison; Davies, Roger L.; Davis, Timothy A.; de Zeeuw, P. T.; Duc, Pierre-Alain; Emsellem, Eric; Khochfar, Sadegh; Krajnović, Davor; Kuntschner, Harald; Lablanche, Pierre-Yves; McDermid, Richard M.; Morganti, Raffaella; Naab, Thorsten; Oosterloo, Tom; Scott, Nicholas; Serra, Paolo; Young, Lisa M.; Weijmans, Anne-Marie
2013-07-01
For early-type galaxies, the ability to sustain a corona of hot, X-ray-emitting gas could have played a key role in quenching their star formation history. A halo of hot gas may act as an effective shield against the acquisition of cold gas and can quickly absorb stellar mass loss material. Yet, since the discovery by the Einstein Observatory of such X-ray haloes around early-type galaxies, the precise amount of hot gas around these galaxies still remains a matter of debate. By combining homogeneously derived photometric and spectroscopic measurements for the early-type galaxies observed as part of the ATLAS3D integral field survey with measurements of their X-ray luminosity based on X-ray data of both low and high spatial resolution (for 47 and 19 objects, respectively) we conclude that the hot gas content of early-type galaxies can depend on their dynamical structure. Specifically, whereas slow rotators generally have X-ray haloes with luminosity LX, gas and temperature T values that are well in line with what is expected if the hot gas emission is sustained by the thermalization of the kinetic energy carried by the stellar mass loss material, fast rotators tend to display LX, gas values that fall consistently below the prediction of this model, with similar T values that do not scale with the stellar kinetic energy (traced by the stellar velocity dispersion) as observed in the case of slow rotators. Such a discrepancy between the hot gas content of slow and fast rotators would appear to reduce, or even disappear, for large values of the dynamical mass (above ˜3 × 1011 M⊙), with younger fast rotators displaying also somewhat larger LX, gas values possibly owing to the additional energy input from recent supernovae explosions. Considering that fast rotators are likely to be intrinsically flatter than slow rotators, and that the few LX, gas-deficient slow rotators also happen to be relatively flat, the observed LX, gas deficiency in these objects would support
NASA Astrophysics Data System (ADS)
Cooper, Wilfred A.; Pfefferle, David; Graves, Jonathan P.
2014-10-01
Designed to accurately solve the motion of energetic particles in the presence of 3D magnetic fields, the VENUS-LEVIS code leans on a non-canonical general coordinate Lagrangian formulation of the equations of motion. It switches between full-orbit particle following and guiding-centre tracing by verifying the perpendicular variation of magnetic vector field, not only including gradients and curvature terms but also the shearing of field-lines. The criteria is particularly relevant for the study of fast ion redistribution in the kinked core of hybrid plasmas, where the compression of flux-surfaces against the axisymmetric outer mantle creates strongly varying magnetic field-lines and large parallel currents. Slowing-down simulations of NBI fast ions show that co-passing particles helically align in the opposite side of the plasma deformation whereas counter-passing particles are barely affected by the kinked structure. Results are compared with experimental neutron camera traces and FIDA measurements during long-lived modes (LLM).
NASA Astrophysics Data System (ADS)
Giles, B. L.; Pollock, C. J.; Avanov, L. A.; Barrie, A. C.; Burch, J. L.; Chandler, M. O.; Clark, G. B.; Coffey, V. N.; Dickson, C.; Dorelli, J.; Ergun, R. E.; Fuselier, S. A.; Gershman, D. J.; Gliese, U.; Holland, M. P.; Jacques, A. D.; Kreisler, S.; Lavraud, B.; MacDonald, E.; Mauk, B.; Moore, T. E.; Mukai, T.; Nakamura, R.; Paterson, W. R.; Rager, A. C.; Saito, Y.; Salo, C.; Sauvaud, J. A.; Torbert, R. B.; Vinas, A. F.; Yokota, S.
2015-12-01
The primary focus of the Magnetospheric Multiscale (MMS) mission, launched in March 2015, is magnetic reconnection and associated processes. Understanding hinges critically on the kinetic physics that allows reconnection to take place. The Fast Plasma Investigation (FPI) provides electron and ion distribution functions at 4.5s cadence and, for select periods of time, at cadences of 30ms for electrons and 150ms for ions. These select time periods are chosen after in situ acquisition based on inspection of the low resolution data. Thus the FPI provides, independent of spacecraft spin rate, the time resolution needed to resolve the small, fast-moving reconnection diffusion regions. The first mission phase focuses on the dayside magnetopause and this presentation is intended to demonstrate the capabilities of FPI to resolve the important spatial scales relevant to the reconnection process. Magnetopause and other boundary crossings will be examined and the phase-space trajectories identified at the tetrahedral satellite locations through analysis of the 3D distribution functions.
NASA Astrophysics Data System (ADS)
Czarske, J. W.; Kuschmierz, R.; Günther, P.
2013-06-01
Precise measurements of distance, eccentricity and 3D-shape of fast moving objects such as turning parts of lathes, gear shafts, magnetic bearings, camshafts, crankshafts and rotors of vacuum pumps are on the one hand important tasks. On the other hand they are big challenges, since contactless precise measurement techniques are required. Optical techniques are well suitable for distance measurements of non-moving surfaces. However, measurements of laterally fast moving surfaces are still challenging. For such tasks the laser Doppler distance sensor technique was invented by the TU Dresden some years ago. This technique has been realized by two mutually tilted interference fringe systems, where the distance is coded in the phase difference between the generated interference signals. However, due to the speckle effect different random envelopes and phase jumps of the interference signals occur. They disturb the phase difference estimation between the interference signals. In this paper, we will report on a scientific breakthrough on the measurement uncertainty budget which has been achieved recently. Via matching of the illumination and receiving optics the measurement uncertainty of the displacement and distance can be reduced by about one magnitude. For displacement measurements of a recurring rough surface a standard deviation of 110 nm were attained at lateral velocities of 5 m / s. Due to the additionally measured lateral velocity and the rotational speed, the two-dimensional shape of rotating objects is calculated. The three-dimensional shape can be conducted by employment of a line camera. Since the measurement uncertainty of the displacement, vibration, distance, eccentricity, and shape is nearly independent of the lateral surface velocity, this technique is predestined for fast-rotating objects. Especially it can be advantageously used for the quality control of workpieces inside of a lathe towards the reduction of process tolerances, installation times and
Mitsouras, Dimitris; Mulkern, Robert V.; Rybicki, Frank J.
2006-01-01
Fast Spin Echo (FSE) trains elicited by non-selective “hard” refocusing radio frequency (RF) pulses have been proposed as a means to enable application of FSE methods for high resolution 3D magnetic resonance imaging (MRI). Hard-pulse FSE (HPFSE) trains offer short (3–4 ms) echo spacings, but are unfortunately limited to imaging the entire sample within the coil sensitivity thus requiring lengthy imaging times, consequently limiting clinical application. In this work we formulate and analyze two general purpose combinations of 3D HPFSE with Inner Volume (IV) MR imaging to circumvent this limitation. The first method employs a 2D selective RF excitation followed by the HPFSE train, and focuses on required properties of the spatial excitation profile with respect to limiting RF pulse duration in the 5–6 ms range. The second method employs two orthogonally selective 1D RF excitations (a 90x°– 180y° pair) to generate an echo from magnetization within the volume defined by their intersection. Subsequent echoes are formed via the HPFSE train, placing the focus of the method on (a) avoiding spurious echoes that may arise from transverse magnetization located outside the slab intersection when it is unavoidably affected by the non-selective refocusing pulses, and (b) avoiding signal losses due to the necessarily different spacing (in time) of the RF pulse applications. The performance of each method is experimentally measured using Carr-Purcell-Meiboom-Gill (CPMG) multi-echo imaging, enabling examination of the magnetization evolution throughout the echo train. The methods as implemented achieve 95% to 97% outer volume signal suppression, and higher suppression appears to be well within reach, by further refinement of the selective RF excitations. Example images of the human brain and spine are presented with each technique. We conclude that the SNR effciency of volume imaging in conjunction with the short echo spacing afforded by hard pulse trains enable high
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.
Optimized design of a LED-array-based TOF range imaging sensor for fast 3-D shape measurement
NASA Astrophysics Data System (ADS)
Wang, Huanqin; Wang, Ying; Xu, Jun; He, Deyong; Zhao, Tianpeng; Ming, Hai; Kong, Deyi
2011-06-01
A LED-array-based range imaging sensor using Time-of-Flight (TOF) distance measurement was developed to capture the depth information of three-dimensional (3-D) object. By time-division electronic scanning of the LED heterodyne phase-shift TOF range finders in array, the range images were fast obtained without any mechanical moving parts. The design of LED-array-based range imaging sensor was adequately described and a range imaging theoretical model based on photoelectric signal processing was built, which showed there was mutual restriction relationship among the measurement time of a depth pixel, the bandwidth of receiver and the sensor's signal-to-noise ratio (SNR). In order to improve the key parameters of sensor such as range resolution and measurement speed simultaneously, some optimized designs needed to be done for the proposed range imaging sensor, including choosing proper parameters for the filters in receiver, adopting special structure feedback automatic gain control (AGC) circuit with short response time, etc. The final experiment results showed the sensor after optimization could acquire the range images at a rate of 10 frames per second with a range resolution as high as +/-2mm in the range of 50-1200mm. The essential advantages of the proposed range imaging sensor were construction with simple structure, high range resolution, short measurement time and low cost, which was sufficient for many robotic and industrial automation applications.
Medical applications of fast 3D cameras in real-time image-guided radiotherapy (IGRT) of cancer
NASA Astrophysics Data System (ADS)
Li, Shidong; Li, Tuotuo; Geng, Jason
2013-03-01
Dynamic volumetric medical imaging (4DMI) has reduced motion artifacts, increased early diagnosis of small mobile tumors, and improved target definition for treatment planning. High speed cameras for video, X-ray, or other forms of sequential imaging allow a live tracking of external or internal movement useful for real-time image-guided radiation therapy (IGRT). However, none of 4DMI can track real-time organ motion and no camera has correlated with 4DMI to show volumetric changes. With a brief review of various IGRT techniques, we propose a fast 3D camera for live-video stereovision, an automatic surface-motion identifier to classify body or respiratory motion, a mechanical model for synchronizing the external surface movement with the internal target displacement by combination use of the real-time stereovision and pre-treatment 4DMI, and dynamic multi-leaf collimation for adaptive aiming the moving target. Our preliminary results demonstrate that the technique is feasible and efficient in IGRT of mobile targets. A clinical trial has been initiated for validation of its spatial and temporal accuracies and dosimetric impact for intensity-modulated RT (IMRT), volumetric-modulated arc therapy (VMAT), and stereotactic body radiotherapy (SBRT) of any mobile tumors. The technique can be extended for surface-guided stereotactic needle insertion in biopsy of small lung nodules.
Comparison of 3-D Modeling With Experimental Results on Fast Wave Antenna Loading in DIII-D
NASA Astrophysics Data System (ADS)
Pinsker, R. I.; Ryan, P. M.; Goulding, R. H.; Hanson, G. R.; Milanesio, D.; Maggiora, R.; Hosea, J. C.; Nagy, A.; Porkolab, M.; Zeng, L.
2011-10-01
In DIII-D and other tokamaks, with a fixed system voltage limit, the parameter that limits the ICRF power that can be coupled to H-mode plasmas is the antenna loading resistance RL. For a fixed antenna geometry and excitation (phasing), RL is determined by the electron density profile in the antenna near-field region. Quantitative understanding of the coupling physics is obtained by comparing the resistive (RL) and reactive components of the antenna loading, without and with plasma, to predictions of 3-D models of the antenna and the edge plasma (Microwave Studio and TOPICA). When measured density profiles from reflectometers are used, good agreement between predicted and measured values of RL is obtained without any adjustable parameters in the model. The improved understanding is applied to enhancement of RL in advanced scenarios in DIII-D to increase the coupled fast wave power. Supported in part by US DOE under DE-FC02-04ER54698, DE-AC05-00OR22725, DE-AC02-09CH11466, DE-FG02-08ER54984.
2016-01-01
Objective. Our objective is to develop a computerized scheme for liver tumor segmentation in MR images. Materials and Methods. Our proposed scheme consists of four main stages. Firstly, the region of interest (ROI) image which contains the liver tumor region in the T1-weighted MR image series was extracted by using seed points. The noise in this ROI image was reduced and the boundaries were enhanced. A 3D fast marching algorithm was applied to generate the initial labeled regions which are considered as teacher regions. A single hidden layer feedforward neural network (SLFN), which was trained by a noniterative algorithm, was employed to classify the unlabeled voxels. Finally, the postprocessing stage was applied to extract and refine the liver tumor boundaries. The liver tumors determined by our scheme were compared with those manually traced by a radiologist, used as the “ground truth.” Results. The study was evaluated on two datasets of 25 tumors from 16 patients. The proposed scheme obtained the mean volumetric overlap error of 27.43% and the mean percentage volume error of 15.73%. The mean of the average surface distance, the root mean square surface distance, and the maximal surface distance were 0.58 mm, 1.20 mm, and 6.29 mm, respectively. PMID:27597960
Le, Trong-Ngoc; Bao, Pham The; Huynh, Hieu Trung
2016-01-01
Objective. Our objective is to develop a computerized scheme for liver tumor segmentation in MR images. Materials and Methods. Our proposed scheme consists of four main stages. Firstly, the region of interest (ROI) image which contains the liver tumor region in the T1-weighted MR image series was extracted by using seed points. The noise in this ROI image was reduced and the boundaries were enhanced. A 3D fast marching algorithm was applied to generate the initial labeled regions which are considered as teacher regions. A single hidden layer feedforward neural network (SLFN), which was trained by a noniterative algorithm, was employed to classify the unlabeled voxels. Finally, the postprocessing stage was applied to extract and refine the liver tumor boundaries. The liver tumors determined by our scheme were compared with those manually traced by a radiologist, used as the "ground truth." Results. The study was evaluated on two datasets of 25 tumors from 16 patients. The proposed scheme obtained the mean volumetric overlap error of 27.43% and the mean percentage volume error of 15.73%. The mean of the average surface distance, the root mean square surface distance, and the maximal surface distance were 0.58 mm, 1.20 mm, and 6.29 mm, respectively. PMID:27597960
Nugroho, Hendro; Widiyantoro, Sri; Nugraha, Andri Dian
2013-09-09
Determination of earthquake hypocenter in Indonesia conducted by the Meteorological, Climatological, and Geophysical Agency (MCGA) has still used a 1-D seismic velocity model. In this research, we have applied a Fast Grid Search (FGM) method and a 3-D velocity model resulting from tomographic imaging to relocate earthquakes in the Sumatran region. The data were taken from the MCGA data catalog from 2009 to 2011 comprising of subduction zone and on land fault earthquakes with magnitude greater than 4 Mw. Our preliminary results show some significant changes in the depths of the relocated earthquakes which are in general deeper than the depths of hypocenters from the MCGA data catalog. The residual times resulting from the relocation process are smaller than those prior to the relocation. Encouraged by these results, we will continue to conduct hypocenter relocation for all events from the MCGA data catalog periodically in order to produce a new data catalog with good quality. We hope that the new data catalog will be useful for further studies.
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.
MO-G-BRF-07: Anomalously Fast Diffusion of Carbon Nanotubes Carriers in 3D Tissue Model
Wang, Y; Bahng, J; Kotov, N
2014-06-15
Purpose: We aim to investigate and understand diffusion process of carbon nanotubes (CNTs) and other nanoscale particles in tissue and organs. Methods: In this research, we utilized a 3D model tissue of hepatocellular carcinoma (HCC)cultured in inverted colloidal crystal (ICC) scaffolds to compare the diffusivity of CNTs with small molecules such as Rhodamine and FITC in vitro, and further investigated the transportation of CNTs with and without targeting ligand, TGFβ1. The real-time permeation profiles of CNTs in HCC tissue model with high temporal and spatial resolution was demonstrated by using standard confocal microscopy. Quantitative analysis of the diffusion process in 3D was carried out using luminescence intensity in a series of Z-stack images obtained for different time points of the diffusion process after initial addition of CNTs or small molecules to the cell culture and the image data was analyzed by software ImageJ and Mathematica. Results: CNTs display diffusion rate in model tissues substantially faster than small molecules of the similar charge such as FITC, and the diffusion rate of CNTs are significantly enhanced with targeting ligand, TGFβ1. Conclusion: In terms of the advantages of in-vitro model, we were able to have access to measuring the rate of CNT penetration at designed conditions with variable parameters. And the findings by using this model, changed our understanding about advantages of CNTs as nanoscale drug carriers and provides design principles for making new drug carriers for both treatment and diagnostics. Additionally the fast diffusion opens the discussion of the best possible drug carriers to reach deep parts of cancerous tissues, which is often a prerequisite for successful cancer treatment. This work was supported by the Center for Photonic and Multiscale Nanomaterials funded by National Science Foundation Materials Research Science and Engineering Center program DMR 1120923. The work was also partially supported by NSF
NASA Astrophysics Data System (ADS)
Canavesi, Cristina; Cogliati, Andrea; Hayes, Adam; Santhanam, Anand P.; Tankam, Patrice; Rolland, Jannick P.
2015-10-01
Fast, robust, nondestructive 3D imaging is needed for characterization of microscopic structures in industrial and clinical applications. A custom micro-electromechanical system (MEMS)-based 2D scanner system was developed to achieve 55 kHz A-scan acquisition in a Gabor-domain optical coherence microscopy (GD-OCM) instrument with a novel multilevel GPU architecture for high-speed imaging. GD-OCM yields high-definition volumetric imaging with dynamic depth of focusing through a bio-inspired liquid lens-based microscope design, which has no moving parts and is suitable for use in a manufacturing setting or in a medical environment. A dual-axis MEMS mirror was chosen to replace two single-axis galvanometer mirrors; as a result, the astigmatism caused by the mismatch between the optical pupil and the scanning location was eliminated and a 12x reduction in volume of the scanning system was achieved. Imaging at an invariant resolution of 2 μm was demonstrated throughout a volume of 1 × 1 × 0.6 mm3, acquired in less than 2 minutes. The MEMS-based scanner resulted in improved image quality, increased robustness and lighter weight of the system - all factors that are critical for on-field deployment. A custom integrated feedback system consisting of a laser diode and a position-sensing detector was developed to investigate the impact of the resonant frequency of the MEMS and the driving signal of the scanner on the movement of the mirror. Results on the metrology of manufactured materials and characterization of tissue samples with GD-OCM are presented.
None, None
2015-09-28
Coulomb interaction between charged particles inside a bunch is one of the most importance collective effects in beam dynamics, becoming even more significant as the energy of the particle beam is lowered to accommodate analytical and low-Z material imaging purposes such as in the time resolved Ultrafast Electron Microscope (UEM) development currently underway at Michigan State University. In addition, space charge effects are the key limiting factor in the development of ultrafast atomic resolution electron imaging and diffraction technologies and are also correlated with an irreversible growth in rms beam emittance due to fluctuating components of the nonlinear electron dynamics. In the short pulse regime used in the UEM, space charge effects also lead to virtual cathode formation in which the negative charge of the electrons emitted at earlier times, combined with the attractive surface field, hinders further emission of particles and causes a degradation of the pulse properties. Space charge and virtual cathode effects and their remediation are core issues for the development of the next generation of high-brightness UEMs. Since the analytical models are only applicable for special cases, numerical simulations, in addition to experiments, are usually necessary to accurately understand the space charge effect. In this paper we will introduce a grid-free differential algebra based multiple level fast multipole algorithm, which calculates the 3D space charge field for n charged particles in arbitrary distribution with an efficiency of O(n), and the implementation of the algorithm to a simulation code for space charge dominated photoemission processes.
None, None
2015-09-28
Coulomb interaction between charged particles inside a bunch is one of the most importance collective effects in beam dynamics, becoming even more significant as the energy of the particle beam is lowered to accommodate analytical and low-Z material imaging purposes such as in the time resolved Ultrafast Electron Microscope (UEM) development currently underway at Michigan State University. In addition, space charge effects are the key limiting factor in the development of ultrafast atomic resolution electron imaging and diffraction technologies and are also correlated with an irreversible growth in rms beam emittance due to fluctuating components of the nonlinear electron dynamics.more » In the short pulse regime used in the UEM, space charge effects also lead to virtual cathode formation in which the negative charge of the electrons emitted at earlier times, combined with the attractive surface field, hinders further emission of particles and causes a degradation of the pulse properties. Space charge and virtual cathode effects and their remediation are core issues for the development of the next generation of high-brightness UEMs. Since the analytical models are only applicable for special cases, numerical simulations, in addition to experiments, are usually necessary to accurately understand the space charge effect. In this paper we will introduce a grid-free differential algebra based multiple level fast multipole algorithm, which calculates the 3D space charge field for n charged particles in arbitrary distribution with an efficiency of O(n), and the implementation of the algorithm to a simulation code for space charge dominated photoemission processes.« less
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.
Ai, H; Pan, T; Hwang, K
2014-06-15
Purpose: To determine the feasibility of identifying cortical bone on MR images with a short-TE 3D fast-GRE sequence for attenuation correction of PET data in PET/MR. Methods: A water-fat-bone phantom was constructed with two pieces of beef shank. MR scans were performed on a 3T MR scanner (GE Discovery™ MR750). A 3D GRE sequence was first employed to measure the level of residual signal in cortical bone (TE{sub 1}/TE{sub 2}/TE{sub 3}=2.2/4.4/6.6ms, TR=20ms, flip angle=25°). For cortical bone segmentation, a 3D fast-GRE sequence (TE/TR=0.7/1.9ms, acquisition voxel size=2.5×2.5×3mm{sup 3}) was implemented along with a 3D Dixon sequence (TE{sub 1}/TE{sub 2}/TR=1.2/2.3/4.0ms, acquisition voxel size=1.25×1.25×3mm{sup 3}) for water/fat imaging. Flip angle (10°), acquisition bandwidth (250kHz), FOV (480×480×144mm{sup 3}) and reconstructed voxel size (0.94×0.94×1.5mm{sup 3}) were kept the same for both sequences. Soft tissue and fat tissue were first segmented on the reconstructed water/fat image. A tissue mask was created by combining the segmented water/fat masks, which was then applied on the fast-GRE image (MRFGRE). A second mask was created to remove the Gibbs artifacts present in regions in close vicinity to the phantom. MRFGRE data was smoothed with a 3D anisotropic diffusion filter for noise reduction, after which cortical bone and air was separated using a threshold determined from the histogram. Results: There is signal in the cortical bone region in the 3D GRE images, indicating the possibility of separating cortical bone and air based on signal intensity from short-TE MR image. The acquisition time for the 3D fast-GRE sequence was 17s, which can be reduced to less than 10s with parallel imaging. The attenuation image created from water-fat-bone segmentation is visually similar compared to reference CT. Conclusion: Cortical bone and air can be separated based on intensity in MR image with a short-TE 3D fast-GRE sequence. Further research is required
The polyGeVero® software for fast and easy computation of 3D radiotherapy dosimetry data
NASA Astrophysics Data System (ADS)
Kozicki, Marek; Maras, Piotr
2015-01-01
The polyGeVero® software package was elaborated for calculations of 3D dosimetry data such as the polymer gel dosimetry. It comprises four workspaces designed for: i) calculating calibrations, ii) storing calibrations in a database, iii) calculating dose distribution 3D cubes, iv) comparing two datasets e.g. a measured one with a 3D dosimetry with a calculated one with the aid of a treatment planning system. To accomplish calculations the software was equipped with a number of tools such as the brachytherapy isotopes database, brachytherapy dose versus distance calculation based on the line approximation approach, automatic spatial alignment of two 3D dose cubes for comparison purposes, 3D gamma index, 3D gamma angle, 3D dose difference, Pearson's coefficient, histograms calculations, isodoses superimposition for two datasets, and profiles calculations in any desired direction. This communication is to briefly present the main functions of the software and report on the speed of calculations performed by polyGeVero®.
3D multiple optical tweezers based on time-shared scanning with a fast focus tunable lens
NASA Astrophysics Data System (ADS)
Tanaka, Yoshio
2013-02-01
Three-dimensional controlled manipulation of individual micro-objects requires multiple optical tweezers that can be independently controlled in a 3D working space with high spatiotemporal resolution. Here, the author presents 3D multiple optical tweezers based on a time-shared scanning technique with an electrically focus tunable lens for axial steering and a two-axis steering mirror for lateral steering. Four typical examples of 3D controlled manipulation, including the rotation of a single bead on its axis, are demonstrated in real time. The optical system design and the control method are also described.
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
NASA Astrophysics Data System (ADS)
Jo, Hong Li; Song, Yo Han; Park, Jinho; Jo, Eun-Jung; Goh, Yeongchang; Shin, Kyujin; Kim, Min-Gon; Lee, Kang Taek
2015-11-01
We report on the development of a three-dimensional (3D) live-cell imaging technique with high spatiotemporal resolution using lanthanide-doped upconverting nanoparticles (UCNPs). It employs the sectioning capability of confocal microscopy except that the two-dimensional (2D) section images are acquired by wide-field epi-fluorescence microscopy. Although epi-fluorescence images are contaminated with the out-of-focus background in general, the near-infrared (NIR) excitation used for the excitation of UCNPs does not generate any autofluorescence, which helps to lower the background. Moreover, the image blurring due to defocusing was naturally eliminated in the image reconstruction process. The 3D images were used to investigate the cellular dynamics such as nuclear uptake and single-particle tracking that require 3D description.We report on the development of a three-dimensional (3D) live-cell imaging technique with high spatiotemporal resolution using lanthanide-doped upconverting nanoparticles (UCNPs). It employs the sectioning capability of confocal microscopy except that the two-dimensional (2D) section images are acquired by wide-field epi-fluorescence microscopy. Although epi-fluorescence images are contaminated with the out-of-focus background in general, the near-infrared (NIR) excitation used for the excitation of UCNPs does not generate any autofluorescence, which helps to lower the background. Moreover, the image blurring due to defocusing was naturally eliminated in the image reconstruction process. The 3D images were used to investigate the cellular dynamics such as nuclear uptake and single-particle tracking that require 3D description. Electronic supplementary information (ESI) available: Methods for synthetic procedure, characterization, live-cell imaging with UCNPs, and two video files for 3D imaging. See DOI: 10.1039/c5nr05875a
Rupprecht, Peter; Prendergast, Andrew; Wyart, Claire; Friedrich, Rainer W
2016-01-01
There is a high demand for 3D multiphoton imaging in neuroscience and other fields but scanning in axial direction presents technical challenges. We developed a focusing technique based on a remote movable mirror that is conjugate to the specimen plane and translated by a voice coil motor. We constructed cost-effective z-scanning modules from off-the-shelf components that can be mounted onto standard multiphoton laser scanning microscopes to extend scan patterns from 2D to 3D. Systems were designed for large objectives and provide high resolution, high speed and a large z-scan range (>300 μm). We used these systems for 3D multiphoton calcium imaging in the adult zebrafish brain and measured odor-evoked activity patterns across >1500 neurons with single-neuron resolution and high signal-to-noise ratio. PMID:27231612
Rupprecht, Peter; Prendergast, Andrew; Wyart, Claire; Friedrich, Rainer W
2016-05-01
There is a high demand for 3D multiphoton imaging in neuroscience and other fields but scanning in axial direction presents technical challenges. We developed a focusing technique based on a remote movable mirror that is conjugate to the specimen plane and translated by a voice coil motor. We constructed cost-effective z-scanning modules from off-the-shelf components that can be mounted onto standard multiphoton laser scanning microscopes to extend scan patterns from 2D to 3D. Systems were designed for large objectives and provide high resolution, high speed and a large z-scan range (>300 μm). We used these systems for 3D multiphoton calcium imaging in the adult zebrafish brain and measured odor-evoked activity patterns across >1500 neurons with single-neuron resolution and high signal-to-noise ratio. PMID:27231612
David L. Alumbaugh
2006-07-03
This report addresses the effects of electrical anisotropy on the 3D inversion of single-well induction logging data when anisotropy is not considered. Of concern are possible artifacts that may lead to an incorrect interpretation of the formation about the borehole. Comparison is made of 3D isotropic inversion on a suite of model data, with and without anisotropy, consisting of an infinite layer and layer terminated at the borehole. In both cases, the layer dip (or well deviation) is varied. Inversion of the anisotropic data result in an overestimate of the layer conductivity, and the lateral extent of the layer about the borehole.
NASA Astrophysics Data System (ADS)
Zhang, M. Q.
1989-09-01
A new Monte Carlo algorithm for 3D Kawasaki spin-exchange simulations and its implementation on a CDC CYBER 205 is presented. This approach is applicable to lattices with sizes between 4×4×4 and 256×L2×L3 ((L2+2)(L3+4)/4⩽65535) and periodic boundary conditions. It is adjustable to various kinetic models in which the total magnetization is conserved. Maximum speed on 10 million steps per second can be reached for 3-D Ising model with Metropolis rate.
NASA Astrophysics Data System (ADS)
Stover, E.; Berger, G.; Wendel, M.; Petter, J.
2015-10-01
A method for non-contact 3D form testing of aspheric surfaces including determination of decenter and wedge errors and lens thickness is presented. The principle is based on the absolute measurement capability of multi-wavelength interferometry (MWLI). The approach produces high density 3D shape information and geometric parameters at high accuracy in short measurement times. The system allows inspection of aspheres without restrictions in terms of spherical departures, of segmented and discontinuous optics. The optics can be polished or ground and made of opaque or transparent materials.
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
VIRO 3D: fast three-dimensional full-body scanning for humans and other living objects
NASA Astrophysics Data System (ADS)
Stein, Norbert; Minge, Bernhard
1998-03-01
The development of a family of partial and whole body scanners provides a complete technology for fully three-dimensional and contact-free scans on human bodies or other living objects within seconds. This paper gives insight into the design and the functional principles of the whole body scanner VIRO 3D operating on the basis of the laser split-beam method. The arrangement of up to 24 camera/laser combinations, thus dividing the area into different camera fields and an all- around sensor configuration travelling in vertical direction allow the complete 360-degree-scan of an object within 6 - 20 seconds. Due to a special calibration process the different sensors are matched and the measured data are combined. Up to 10 million 3D measuring points with a resolution of approximately 1 mm are processed in all coordinate axes to generate a 3D model. By means of high-performance processors in combination with real-time image processing chips the image data from almost any number of sensors can be recorded and evaluated synchronously in video real-time. VIRO 3D scanning systems have already been successfully implemented in various applications and will open up new perspectives in different other fields, ranging from industry, orthopaedic medicine, plastic surgery to art and photography.
3D in the Fast Lane: Render as You Go with the Latest OpenGL Boards.
ERIC Educational Resources Information Center
Sauer, Jeff; Murphy, Sam
1997-01-01
NT OpenGL hardware allows modelers and animators to work at relatively inexpensive NT workstations in their own offices or homes previous to shared space and workstation time in expensive studios. Rates seven OpenGL boards and two QuickDraw 3D accelerator boards for Mac users on overall value, wireframe and texture rendering, 2D acceleration, and…
NASA Astrophysics Data System (ADS)
Bates, Jason; Schmitt, Andrew; Klapisch, Marcel; Karasik, Max; Obenschain, Steve
2013-10-01
Modifications to the FAST3D code have been made to enhance its ability to simulate the dynamics of plastic ICF targets with high-Z overcoats. This class of problems is challenging computationally due in part to plasma conditions that are not in a state of local thermodynamic equilibrium and to the presence of mixed computational cells containing more than one material. Recently, new opacity tables for gold, palladium and plastic have been generated with an improved version of the STA code. These improved tables provide smoother, higher-fidelity opacity data over a wider range of temperature and density states than before, and contribute to a more accurate treatment of radiative transfer processes in FAST3D simulations. Furthermore, a new, more efficient subroutine known as ``MMEOS'' has been installed in the FAST3D code for determining pressure and temperature equilibrium conditions within cells containing multiple materials. We will discuss these topics, and present new simulation results for high-Z planar-target experiments performed recently on the NIKE Laser Facility. Work supported by DOE/NNSA.
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
3D pyCloudy modelling of bipolar planetary nebulae: Evidence for fast fading of the lobes
NASA Astrophysics Data System (ADS)
Gesicki, K.; Zijlstra, A. A.; Morisset, C.
2016-01-01
Aims: The origin and evolution of the shapes of bipolar planetary nebulae are poorly understood. We postulate that their history can be traced through their internal velocity fields in a procedure similar to the one well established for spherical objects. Such an analysis requires 3D photoionization and kinematical modelling that is computationally very time consuming. We apply an axially symmetric pseudo-3D photoionization model, pyCloudy, to derive the structures of six bipolar nebulae and two suggested post-bipolars in an attempt to constrain the bipolar planetary nebulae evolution. Methods: HST images and VLT/UVES spectroscopy are used for the modelling. The targets are located in the direction of the Galactic bulge. A 3D model structure is used as input to the photoionization code in order to fit the HST images. Line profiles of different ions constrain the velocity field. The model and associated velocity fields allow us to derive masses, velocities, and ages. Results: The 3D models find much lower ionized masses than required in 1D models: ionized masses are reduced by factors of 2-7. The selected bi-lobed planetary nebulae show a narrow range of ages: the averaged radii and velocities result in values between 1300 and 2000 yr. The lobes are fitted well with velocities that increase linearly with radius. These Hubble-type flows have been found before, and suggest that the lobes form at a defined point in time. The lobes appear to be slightly younger, by ~500 yr, than the main (host) nebulae; they seem to form at an early phase of PN evolution and fade after 1-2 kyr. We find that 30-35% of bulge PNe pass through a bipolar phase. Based on observations collected at the European Organisation for Astronomical Research in the Southern Hemisphere, Chile (proposal 075.D-0104) and HST (program 9356).
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.
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.
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.
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
Bao Yidong; Hu Sibo; Lang Zhikui; Hu Ping
2005-08-05
A fast simulation scheme for 3D curved binder flanging and blank shape prediction of sheet metal based on one-step inverse finite element method is proposed, in which the total plasticity theory and proportional loading assumption are used. The scheme can be actually used to simulate 3D flanging with complex curve binder shape, and suitable for simulating any type of flanging model by numerically determining the flanging height and flanging lines. Compared with other methods such as analytic algorithm and blank sheet-cut return method, the prominent advantage of the present scheme is that it can directly predict the location of the 3D flanging lines when simulating the flanging process. Therefore, the prediction time of flanging lines will be obviously decreased. Two typical 3D curve binder flanging including stretch and shrink characters are simulated in the same time by using the present scheme and incremental FE non-inverse algorithm based on incremental plasticity theory, which show the validity and high efficiency of the present scheme.
NASA Astrophysics Data System (ADS)
Huang, Guo-Jiao; Bai, Chao-Ying; Greenhalgh, Stewart
2013-09-01
The traditional grid/cell-based wavefront expansion algorithms, such as the shortest path algorithm, can only find the first arrivals or multiply reflected (or mode converted) waves transmitted from subsurface interfaces, but cannot calculate the other later reflections/conversions having a minimax time path. In order to overcome the above limitations, we introduce the concept of a stationary minimax time path of Fermat's Principle into the multistage irregular shortest path method. Here we extend it from Cartesian coordinates for a flat earth model to global ray tracing of multiple phases in a 3-D complex spherical earth model. The ray tracing results for 49 different kinds of crustal, mantle and core phases show that the maximum absolute traveltime error is less than 0.12 s and the average absolute traveltime error is within 0.09 s when compared with the AK135 theoretical traveltime tables for a 1-D reference model. Numerical tests in terms of computational accuracy and CPU time consumption indicate that the new scheme is an accurate, efficient and a practical way to perform 3-D multiphase arrival tracking in regional or global traveltime tomography.
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 Astrophysics Data System (ADS)
Khochfar, Sadegh; Emsellem, Eric; Serra, Paolo; Bois, Maxime; Alatalo, Katherine; Bacon, R.; Blitz, Leo; Bournaud, Frédéric; Bureau, M.; Cappellari, Michele; Davies, Roger L.; Davis, Timothy A.; de Zeeuw, P. T.; Duc, Pierre-Alain; Krajnović, Davor; Kuntschner, Harald; Lablanche, Pierre-Yves; McDermid, Richard M.; Morganti, Raffaella; Naab, Thorsten; Oosterloo, Tom; Sarzi, Marc; Scott, Nicholas; Weijmans, Anne-Marie; Young, Lisa M.
2011-10-01
We propose a simple model for the origin of fast and slow rotator early-type galaxies (ETG) within the hierarchical Λcold dark matter (ΛCDM) scenario, that is based on the assumption that the mass fraction of stellar discs in ETGs is a proxy for the specific angular momentum expressed via λR. Within our model we reproduce the fraction of fast and slow rotators as a function of magnitude in the ATLAS3D survey, assuming that fast-rotating ETGs have at least 10 per cent of their total stellar mass in a disc component. In agreement with ATLAS3D observations we find that slow rotators are predominantly galaxies with M* > 1010.5 M⊙ contributing ˜20 per cent to the overall ETG population. We show in detail that the growth histories of fast and slow rotators are different, supporting the classification of ETGs into these two categories. Slow rotators accrete between ˜50 and 90 per cent of their stellar mass from satellites and their most massive progenitors have on average up to three major mergers during their evolution. Fast rotators in contrast accrete less than 50 per cent and have on average less than one major merger in their past. We find that the underlying physical reason for the different growth histories is the slowing down and ultimately complete shut-down of gas cooling in massive galaxies. Once cooling and associated star formation in disc stop, galaxies grow via infall from satellites. Frequent minor mergers thereby destroy existing stellar discs via violent relaxation and also tend to lower the specific angular momentum of the main stellar body, lowering λR into the slow rotator regime. On average, the last gas-rich major merger interaction in slow rotators happens at z > 1.5, followed by a series of minor mergers. These results support the idea that kinematically decoupled cores (KDC) form during gas-rich major mergers at high z followed by minor mergers, which build-up the outer layers of the remnant, and make remnants that are initially too flat
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 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.
NASA Astrophysics Data System (ADS)
Aghdasi, Nava; Li, Yangming; Berens, Angelique; Moe, Kris S.; Bly, Randall A.; Hannaford, Blake
2015-03-01
Minimally invasive neuroendoscopic surgery provides an alternative to open craniotomy for many skull base lesions. These techniques provides a great benefit to the patient through shorter ICU stays, decreased post-operative pain and quicker return to baseline function. However, density of critical neurovascular structures at the skull base makes planning for these procedures highly complex. Furthermore, additional surgical portals are often used to improve visualization and instrument access, which adds to the complexity of pre-operative planning. Surgical approach planning is currently limited and typically involves review of 2D axial, coronal, and sagittal CT and MRI images. In addition, skull base surgeons manually change the visualization effect to review all possible approaches to the target lesion and achieve an optimal surgical plan. This cumbersome process relies heavily on surgeon experience and it does not allow for 3D visualization. In this paper, we describe a rapid pre-operative planning system for skull base surgery using the following two novel concepts: importance-based highlight and mobile portal. With this innovation, critical areas in the 3D CT model are highlighted based on segmentation results. Mobile portals allow surgeons to review multiple potential entry portals in real-time with improved visualization of critical structures located inside the pathway. To achieve this we used the following methods: (1) novel bone-only atlases were manually generated, (2) orbits and the center of the skull serve as features to quickly pre-align the patient's scan with the atlas, (3) deformable registration technique was used for fine alignment, (4) surgical importance was assigned to each voxel according to a surgical dictionary, and (5) pre-defined transfer function was applied to the processed data to highlight important structures. The proposed idea was fully implemented as independent planning software and additional
SU-E-J-251: Fast MR-Based DRR Generation Using Highly Undersampled 3D Radial Trajectories
Pereira, G; Traughber, B; Traughber, M; Hu, L; Su, K; Muzic, R
2014-06-01
Purpose: The construction of a digitally reconstructed radiograph (DRR) from a magnetic resonance image (MRI) is possible if the cortical bone signal can be acquired and separated from air and soft tissue. This may be accomplished by subtracting a long echo-time, in-phase, gradient echo (GRE) image volume from an ultra-short echo time free induction decay (FID) image to produce a bone-enhanced (BE) image that reveals cortical bone. One limitation of this approach is the length of time required for data acquisition, which can limit the quality of the DRRs due to patient and organ motion. This study aimed to significantly reduce the acquisition time without compromising DRR quality. Methods: Brain data were acquired from two volunteers using a 3T MR scanner (Ingenia, Philips Healthcare). The FID and GRE images were acquired in a single acquisition using a 3D radial readout sequence with the following parameters: TE1=0.142ms (ultra-short), TE2=2.197ms (nearly in-phase), 2*2*2mm3 isotropic voxels, 250*250*250mm3 FOV. To reduce the acquisition time, k-space was sampled at 75, 50 and 25% of a full 3D sphere . The TE2 image was subtracted from the TE1 image to generate the BE images. The BE images were used to generate DRRs using the Pinnacle treatment planning system (Philips-version 9.2). The quality of the DRRs was evaluated qualitatively by 5 board certified medical physicists for clinical usefulness. Results: The acquisition time for 75, 50 and 25% sampling schemes were 219s, 146s, and 73s, respectively, the latter of which was a four-fold reduction in scan time compared to a 300s fully-sampled acquisition. All DRRs obtained were of acceptable quality and were shown to have sufficient information for clinical 2D image matching. Conclusion: Undersampling k-space while maintaining the same range of frequency information results in significantly reduced scan time and clinically acceptable DRR image quality. Drs. B Traughber and R Muzic have research support from Philips
Jiang, Chao; Jia, Shuhai; Dong, Jun; Bao, Qingchen; Yang, Jia; Lian, Qin; Li, Dichen
2015-09-21
We propose a novel multi-frequency color-marked fringe projection profilometry approach to measure the 3D shape of objects with depth discontinuities. A digital micromirror device projector is used to project a color map consisting of a series of different-frequency color-marked fringe patterns onto the target object. We use a chromaticity curve to calculate the color change caused by the height of the object. The related algorithm to measure the height is also described in this paper. To improve the measurement accuracy, a chromaticity curve correction method is presented. This correction method greatly reduces the influence of color fluctuations and measurement error on the chromaticity curve and the calculation of the object height. The simulation and experimental results validate the utility of our method. Our method avoids the conventional phase shifting and unwrapping process, as well as the independent calculation of the object height required by existing techniques. Thus, it can be used to measure complex and dynamic objects with depth discontinuities. These advantages are particularly promising for industrial applications. PMID:26406621
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.''
NASA Astrophysics Data System (ADS)
Chixiao, Chen; Jixuan, Xiang; Huabin, Chen; Jun, Xu; Fan, Ye; Ning, Li; Junyan, Ren
2015-05-01
Asynchronous successive approximation register (SAR) analog-to-digital converters (ADC) feature high energy efficiency but medium performance. From the point of view of speed, the key bottleneck is the unit capacitor size. In this paper, a small size three-dimensional (3-D) metal—oxide—metal (MOM) capacitor is proposed. The unit capacitor has a capacitance of 1-fF. It shapes as an umbrella, which is designed for fast settling consideration. A comparison among the proposed capacitor with other 3-D MOM capacitors is also given in the paper. To demonstrate the effectiveness of the MOM capacitor, a 6-b capacitive DAC is implemented in TSMC 1P9M 65 nm LP CMOS technology. The DAC consumes a power dissipation of 0.16 mW at the rate of 100 MS/s, excluding a source-follower based output buffer. Static measurement result shows that INL is less than ±1 LSB and DNL is less than ±0.5 LSB. In addition, a 100 MS/s 9-bit SAR ADC with the proposed 3-D capacitor is simulated.
NASA Astrophysics Data System (ADS)
Chai, Kil-Byoung; Zhai, Xiang; Bellan, Paul
2015-11-01
The Caltech astrophysical jet experiment provides a highly resolved demonstration of the interaction between single fluid and 2-fluid scales and possibly kinetic scales as well. The jet evolves through the following sequence: (i) a current-carrying MHD-driven plasma jet self-forms, (ii) the jet undergoes a kink instability, (iii) the kink provides the environment for a secondary, Rayleigh-Taylor (RT) instability, (iv) the RT instability erodes the current channel radius to a scale smaller than ion skin depth to cause fast magnetic reconnection, (v) the reconnection emits broadband obliquely-propagating, right-hand circularly polarized whistler waves, and (vi) the reconnection energizes electrons and ions. The observation of the whistler waves confirms that the reconnection is in the Hall MHD regime (i.e., 2-fluid). An energetic extreme ultraviolet burst is observed at the location of reconnection indicating strong, localized electron heating. Spectroscopic measurement shows simultaneous fast ion heating. The analysis shows that electrons are plausibly heated by Ohmic dissipation, and that ions are plausibly heated stochastically, i.e., the guiding center approximation fails, a kinetic effect.
Lukas, Carsten; Hahn, Horst K; Bellenberg, Barbara; Rexilius, Jan; Schmid, Gebhard; Schimrigk, Sebastian K; Przuntek, Horst; Köster, Odo; Peitgen, Heinz-Otto
2004-11-01
Fast, reliable and easy-to-use methods to quantify brain atrophy are of increasing importance in clinical studies on neuro-degenerative diseases. Here, ILAB 4, a new volumetry software that uses a fast semi-automated 3D segmentation of thin-slice T1-weighted 3D MR images based on a modified watershed transform and an automatic histogram analysis was evaluated. It provides the cerebral volumes: whole brain, white matter, gray matter and intracranial cavity. Inter- and intra-rater reliability and scan-rescan reproducibility were excellent in measuring whole brain volumes (coefficients of variation below 0.5%) of volunteers and patients. However, gray and white matter volumes were more susceptible to image quality. High accuracy of the absolute volume results (+/-5 ml) were shown by phantom and preparation measurements. Analysis times were 6 min for processing of 128 slices. The proposed technique is reliable and highly suitable for quantitative studies of brain atrophy, e.g., in multiple sclerosis. PMID:15536555
NASA Astrophysics Data System (ADS)
Bates, Jason; Schmitt, Andrew; Zalesak, Steve
2015-11-01
The ablative Rayleigh-Taylor (RT) instability is a key factor in the performance of directly-drive inertial-confinement-fusion (ICF) targets. Although this subject has been studied for quite some time, the accurate simulation of the ablative RT instability has proven to be a challenging task for many radiation hydrodynamics codes, particularly when it comes to capturing the ablatively-stabilized region of the linear dispersion spectrum and modeling ab initio perturbations. In this poster, we present results from recent two-dimensional numerical simulations of the ablative RT instability that were performed using the Eulerian code FastRad3D at the U.S. Naval Research Laboratory. We consider both planar and spherical geometries, low and moderate-Z target materials, different laser wavelengths and where possible, compare our findings with experiment data, linearized theory and/or results from other radiation hydrodynamics codes. Overall, we find that FastRad3D is capable of simulating the ablative RT instability quite accurately, although some uncertainties/discrepancies persist. We discuss these issues, as well as some of the numerical challenges associated with modeling this class of problems. Work supported by U.S. DOE/NNSA.
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.
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.
A fast rigid-registration method of inferior limb X-ray image and 3D CT images for TKA surgery
NASA Astrophysics Data System (ADS)
Ito, Fumihito; O. D. A, Prima; Uwano, Ikuko; Ito, Kenzo
2010-03-01
In this paper, we propose a fast rigid-registration method of inferior limb X-ray films (two-dimensional Computed Radiography (CR) images) and three-dimensional Computed Tomography (CT) images for Total Knee Arthroplasty (TKA) surgery planning. The position of the each bone, such as femur and tibia (shin bone), in X-ray film and 3D CT images is slightly different, and we must pay attention how to use the two different images, since X-ray film image is captured in the standing position, and 3D CT is captured in decubitus (face up) position, respectively. Though the conventional registration mainly uses cross-correlation function between two images,and utilizes optimization techniques, it takes enormous calculation time and it is difficult to use it in interactive operations. In order to solve these problems, we calculate the center line (bone axis) of femur and tibia (shin bone) automatically, and we use them as initial positions for the registration. We evaluate our registration method by using three patient's image data, and we compare our proposed method and a conventional registration, which uses down-hill simplex algorithm. The down-hill simplex method is an optimization algorithm that requires only function evaluations, and doesn't need the calculation of derivatives. Our registration method is more effective than the downhill simplex method in computational time and the stable convergence. We have developed the implant simulation system on a personal computer, in order to support the surgeon in a preoperative planning of TKA. Our registration method is implemented in the simulation system, and user can manipulate 2D/3D translucent templates of implant components on X-ray film and 3D CT images.
NASA Astrophysics Data System (ADS)
Liang, Xian-hua; Sun, Wei-dong
2011-06-01
Inventory checking is one of the most significant parts for grain reserves, and plays a very important role on the macro-control of food and food security. Simple, fast and accurate method to obtain internal structure information and further to estimate the volume of the grain storage is needed. Here in our developed system, a special designed multi-site laser scanning system is used to acquire the range data clouds of the internal structure of the grain storage. However, due to the seriously uneven distribution of the range data, this data should firstly be preprocessed by an adaptive re-sampling method to reduce the data redundancy as well as noise. Then the range data is segmented and useful features, such as plane and cylinder information, are extracted. With these features a coarse registration between all of these single-site range data is done, and then an Iterative Closest Point (ICP) algorithm is carried out to achieve fine registration. Taking advantage of the structure of the grain storage being well defined and the types of them are limited, a fast automatic registration method based on the priori model is proposed to register the multi-sites range data more efficiently. Then after the integration of the multi-sites range data, the grain surface is finally reconstructed by a delaunay based algorithm and the grain volume is estimated by a numerical integration method. This proposed new method has been applied to two common types of grain storage, and experimental results shown this method is more effective and accurate, and it can also avoids the cumulative effect of errors when registering the overlapped area pair-wisely.
NASA Astrophysics Data System (ADS)
Zhang, He
2013-01-01
The space charge effect is one of the most important collective effects in beam dynamic studies. In many cases, numerical simulations are inevitable in order to get a clear understanding of this effect. The particle-particle interaction algorithms and the article-in-cell algorithms are widely used in space charge effect simulations. But they both have difficulties in dealing with highly correlated beams with abnormal distributions or complicated geometries. We developed a new algorithm to calculate the three dimensional self-field between charged particles by combining the differential algebra (DA) techniques with the fast multi-pole method (FMM). The FMM hierarchically decomposes the whole charged domain into many small regions. For each region it uses multipole expansions to represent the potential/field contributions from the particles far away from the region and then converts the multipole expansions into a local expansion inside the region. The potential/field due to the far away particles is calculated from the expansions and the potential/field due to the nearby particles is calculated from the Coulomb force law. The DA techniques are used in the calculation, translation and converting of the expansions. The new algorithm scales linearly with the total number of particles and it is suitable for any arbitrary charge distribution. Using the DA techniques, we can calculate both the potential/field and its high order derivatives, which will be useful for the purpose of including the space charge effect into transfer maps in the future. We first present the single level FMM, which decomposes the whole domain into boxes of the same size. It works best for charge distributions that are not overly non-uniform. Then we present the multilevel fast multipole algorithm (MLFMA), which decomposes the whole domain into different sized boxes according to the charge density. Finer boxes are generated where the higher charge density exists; thus the algorithm works for any
NASA Astrophysics Data System (ADS)
Spuhler, Christoph; Harders, Matthias; Székely, Gábor
2006-03-01
We present a fast and robust approach for automatic centerline extraction of tubular structures. The underlying idea is to cut traditional snakes into a set of shorter, independent segments - so-called snakelets. Following the same variational principles, each snakelet acts locally and extracts a subpart of the overall structure. After a parallel optimization step, outliers are detected and the remaining segments then form an implicit centerline. No manual initialization of the snakelets is necessary, which represents one advantage of the method. Moreover, computational complexity does not directly depend on dataset size, but on the number of snake segments necessary to cover the structure of interest, resulting in short computation times. Lastly, the approach is robust even for very complex datasets such as the small intestine. Our approach was tested on several medical datasets (CT datasets of colon, small bowel, and blood vessels) and yielded smooth, connected centerlines with few or no branches. The computation time needed is less than a minute using standard computing hardware.
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
The ATLAS3D project - XXVI. H I discs in real and simulated fast and slow rotators
NASA Astrophysics Data System (ADS)
Serra, Paolo; Oser, Ludwig; Krajnović, Davor; Naab, Thorsten; Oosterloo, Tom; Morganti, Raffaella; Cappellari, Michele; Emsellem, Eric; Young, Lisa M.; Blitz, Leo; Davis, Timothy A.; Duc, Pierre-Alain; Hirschmann, Michaela; Weijmans, Anne-Marie; Alatalo, Katherine; Bayet, Estelle; Bois, Maxime; Bournaud, Frédéric; Bureau, Martin; Crocker, Alison F.; Davies, Roger L.; de Zeeuw, P. T.; Khochfar, Sadegh; Kuntschner, Harald; Lablanche, Pierre-Yves; McDermid, Richard M.; Sarzi, Marc; Scott, Nicholas
2014-11-01
One quarter of all nearby early-type galaxies (ETGs) outside Virgo host a disc/ring of H I with size from a few to tens of kpc and mass up to ˜109 M⊙. Here we investigate whether this H I is related to the presence of a stellar disc within the host making use of the classification of ETGs in fast and slow rotators (FR/SR). We find a large diversity of H I masses and morphologies within both families. Surprisingly, SRs are detected as often, host as much H I and have a similar rate of H I discs/rings as FRs. Accretion of H I is therefore not always linked to the growth of an inner stellar disc. The weak relation between H I and stellar disc is confirmed by their frequent kinematical misalignment in FRs, including cases of polar and counterrotating gas. In SRs the H I is usually polar. This complex picture highlights a diversity of ETG formation histories which may be lost in the relative simplicity of their inner structure and emerges when studying their outer regions. We find that Λ CDM hydrodynamical simulations have difficulties reproducing the H I properties of ETGs. The gas discs formed in simulations are either too massive or too small depending on the star formation feedback implementation. Kinematical misalignments match the observations only qualitatively. The main point of conflict is that nearly all simulated FRs and a large fraction of all simulated SRs host corotating H I. This establishes the H I properties of ETGs as a novel challenge to simulations.
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.
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.
Merino-Caviedes, Susana; Cordero-Grande, Lucilio; Revilla-Orodea, Ana; Sevilla-Ruiz, Teresa; Pérez, M Teresa; Martín-Fernández, Marcos; Alberola-López, Carlos
2014-01-01
We propose a fully 3-D methodology for the computation of myocardial nonviable tissue transmurality in contrast enhanced magnetic resonance images. The outcome is a continuous map defined within the myocardium where not only current state-of-the-art measures of transmurality can be calculated, but also information on the location of nonviable tissue is preserved. The computation is done by means of a partial differential equation framework we have called multi-stencil streamline fast marching. Using it, the myocardial and scarred tissue thickness is simultaneously computed. Experimental results show that the proposed 3-D method allows for the computation of transmurality in myocardial regions where current 2-D methods are not able to as conceived, and it also provides more robust and accurate results in situations where the assumptions on which current 2-D methods are based-i.e., there is a visible endocardial contour and its corresponding epicardial points lie on the same slice-, are not met. PMID:24235299
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.