Linearly Scaling 3D Fragment Method for Large-Scale Electronic Structure Calculations
Wang, Lin-Wang; Lee, Byounghak; Shan, Hongzhang; Zhao, Zhengji; Meza, Juan; Strohmaier, Erich; Bailey, David H.
2008-07-01
We present a new linearly scaling three-dimensional fragment (LS3DF) method for large scale ab initio electronic structure calculations. LS3DF is based on a divide-and-conquer approach, which incorporates a novel patching scheme that effectively cancels out the artificial boundary effects due to the subdivision of the system. As a consequence, the LS3DF program yields essentially the same results as direct density functional theory (DFT) calculations. The fragments of the LS3DF algorithm can be calculated separately with different groups of processors. This leads to almost perfect parallelization on tens of thousands of processors. After code optimization, we were able to achieve 35.1 Tflop/s, which is 39percent of the theoretical speed on 17,280 Cray XT4 processor cores. Our 13,824-atom ZnTeO alloy calculation runs 400 times faster than a direct DFTcalculation, even presuming that the direct DFT calculation can scale well up to 17,280 processor cores. These results demonstrate the applicability of the LS3DF method to material simulations, the advantage of using linearly scaling algorithms over conventional O(N3) methods, and the potential for petascale computation using the LS3DF method.
The linearly scaling 3D fragment method for large scale electronic structure calculations
Zhao, Zhengji; Meza, Juan; Lee, Byounghak; Shan, Hongzhang; Strohmaier, Erich; Bailey, David; Wang, Lin-Wang
2009-07-28
The Linearly Scaling three-dimensional fragment (LS3DF) method is an O(N) ab initio electronic structure method for large-scale nano material simulations. It is a divide-and-conquer approach with a novel patching scheme that effectively cancels out the artificial boundary effects, which exist in all divide-and-conquer schemes. This method has made ab initio simulations of thousand-atom nanosystems feasible in a couple of hours, while retaining essentially the same accuracy as the direct calculation methods. The LS3DF method won the 2008 ACM Gordon Bell Prize for algorithm innovation. Our code has reached 442 Tflop/s running on 147,456 processors on the Cray XT5 (Jaguar) at OLCF, and has been run on 163,840 processors on the Blue Gene/P (Intrepid) at ALCF, and has been applied to a system containing 36,000 atoms. In this paper, we will present the recent parallel performance results of this code, and will apply the method to asymmetric CdSe/CdS core/shell nanorods, which have potential applications in electronic devices and solar cells.
The Linearly Scaling 3D Fragment Method for Large Scale Electronic Structure Calculations
Zhao, Zhengji; Meza, Juan; Lee, Byounghak; Shan, Hongzhang; Strohmaier, Erich; Bailey, David; Wang, Lin-Wang
2009-06-26
The Linearly Scaling three-dimensional fragment (LS3DF) method is an O(N) ab initio electronic structure method for large-scale nano material simulations. It is a divide-and-conquer approach with a novel patching scheme that effectively cancels out the artificial boundary effects, which exist in all divide-and-conquer schemes. This method has made ab initio simulations of thousand-atom nanosystems feasible in a couple of hours, while retaining essentially the same accuracy as the direct calculation methods. The LS3DF method won the 2008 ACM Gordon Bell Prize for algorithm innovation. Our code has reached 442 Tflop/s running on 147,456 processors on the Cray XT5 (Jaguar) at OLCF, and has been run on 163,840 processors on the Blue Gene/P (Intrepid) at ALCF, and has been applied to a system containing 36,000 atoms. In this paper, we will present the recent parallel performance results of this code, and will apply the method to asymmetric CdSe/CdS core/shell nanorods, which have potential applications in electronic devices and solar cells.
The linearly scaling 3D fragment method for large scale electronic structure calculations
NASA Astrophysics Data System (ADS)
Zhao, Zhengji; Meza, Juan; Lee, Byounghak; Shan, Hongzhang; Strohmaier, Erich; Bailey, David; Wang, Lin-Wang
2009-07-01
The linearly scaling three-dimensional fragment (LS3DF) method is an O(N) ab initio electronic structure method for large-scale nano material simulations. It is a divide-and-conquer approach with a novel patching scheme that effectively cancels out the artificial boundary effects, which exist in all divide-and-conquer schemes. This method has made ab initio simulations of thousand-atom nanosystems feasible in a couple of hours, while retaining essentially the same accuracy as the direct calculation methods. The LS3DF method won the 2008 ACM Gordon Bell Prize for algorithm innovation. Our code has reached 442 Tflop/s running on 147,456 processors on the Cray XT5 (Jaguar) at OLCF, and has been run on 163,840 processors on the Blue Gene/P (Intrepid) at ALCF, and has been applied to a system containing 36,000 atoms. In this paper, we will present the recent parallel performance results of this code, and will apply the method to asymmetric CdSe/CdS core/shell nanorods, which have potential applications in electronic devices and solar cells.
3D puzzle reconstruction for archeological fragments
NASA Astrophysics Data System (ADS)
Jampy, F.; Hostein, A.; Fauvet, E.; Laligant, O.; Truchetet, F.
2015-03-01
The reconstruction of broken artifacts is a common task in archeology domain; it can be supported now by 3D data acquisition device and computer processing. Many works have been dedicated in the past to reconstructing 2D puzzles but very few propose a true 3D approach. We present here a complete solution including a dedicated transportable 3D acquisition set-up and a virtual tool with a graphic interface allowing the archeologists to manipulate the fragments and to, interactively, reconstruct the puzzle. The whole lateral part is acquired by rotating the fragment around an axis chosen within a light sheet thanks to a step-motor synchronized with the camera frame clock. Another camera provides a top view of the fragment under scanning. A scanning accuracy of 100μm is attained. The iterative automatic processing algorithm is based on segmentation into facets of the lateral part of the fragments followed by a 3D matching providing the user with a ranked short list of possible assemblies. The device has been applied to the reconstruction of a set of 1200 fragments from broken tablets supporting a Latin inscription dating from the first century AD.
Automated building of organometallic complexes from 3D fragments.
Foscato, Marco; Venkatraman, Vishwesh; Occhipinti, Giovanni; Alsberg, Bjørn K; Jensen, Vidar R
2014-07-28
A method for the automated construction of three-dimensional (3D) molecular models of organometallic species in design studies is described. Molecular structure fragments derived from crystallographic structures and accurate molecular-level calculations are used as 3D building blocks in the construction of multiple molecular models of analogous compounds. The method allows for precise control of stereochemistry and geometrical features that may otherwise be very challenging, or even impossible, to achieve with commonly available generators of 3D chemical structures. The new method was tested in the construction of three sets of active or metastable organometallic species of catalytic reactions in the homogeneous phase. The performance of the method was compared with those of commonly available methods for automated generation of 3D models, demonstrating higher accuracy of the prepared 3D models in general, and, in particular, a much wider range with respect to the kind of chemical structures that can be built automatically, with capabilities far beyond standard organic and main-group chemistry. PMID:24998944
New method of 3-D object recognition
NASA Astrophysics Data System (ADS)
He, An-Zhi; Li, Qun Z.; Miao, Peng C.
1991-12-01
In this paper, a new method of 3-D object recognition using optical techniques and a computer is presented. We perform 3-D object recognition using moire contour to obtain the object's 3- D coordinates, projecting drawings of the object in three coordinate planes to describe it and using a method of inquiring library of judgement to match objects. The recognition of a simple geometrical entity is simulated by computer and studied experimentally. The recognition of an object which is composed of a few simple geometrical entities is discussed.
Methods for comparing 3D surface attributes
NASA Astrophysics Data System (ADS)
Pang, Alex; Freeman, Adam
1996-03-01
A common task in data analysis is to compare two or more sets of data, statistics, presentations, etc. A predominant method in use is side-by-side visual comparison of images. While straightforward, it burdens the user with the task of discerning the differences between the two images. The user if further taxed when the images are of 3D scenes. This paper presents several methods for analyzing the extent, magnitude, and manner in which surfaces in 3D differ in their attributes. The surface geometry are assumed to be identical and only the surface attributes (color, texture, etc.) are variable. As a case in point, we examine the differences obtained when a 3D scene is rendered progressively using radiosity with different form factor calculation methods. The comparison methods include extensions of simple methods such as mapping difference information to color or transparency, and more recent methods including the use of surface texture, perturbation, and adaptive placements of error glyphs.
NASA Technical Reports Server (NTRS)
Scruggs, T.; Moraguez, M.; Patankar, K.; Fitz-Coy, N.; Liou, J.-C.; Sorge, M.; Huynh, T.
2016-01-01
Debris fragments from the hypervelocity impact testing of DebriSat are being collected and characterized for use in updating existing satellite breakup models. One of the key parameters utilized in these models is the ballistic coefficient of the fragment which is directly related to its area-to-mass ratio. However, since the attitude of fragments varies during their orbital lifetime, it is customary to use the average cross-sectional area in the calculation of the area-to-mass ratio. The average cross-sectional area is defined as the average of the projected surface areas perpendicular to the direction of motion and has been shown to be equal to one-fourth of the total surface area of a convex object. Unfortunately, numerous fragments obtained from the DebriSat experiment show significant concavity (i.e., shadowing) and thus we have explored alternate methods for computing the average cross-sectional area of the fragments. An imaging system based on the volumetric reconstruction of a 3D object from multiple 2D photographs of the object was developed for use in determining the size characteristic (i.e., characteristics length) of the DebriSat fragments. For each fragment, the imaging system generates N number of images from varied azimuth and elevation angles and processes them using a space-carving algorithm to construct a 3D point cloud of the fragment. This paper describes two approaches for calculating the average cross-sectional area of debris fragments based on the 3D imager. Approach A utilizes the constructed 3D object to generate equally distributed cross-sectional area projections and then averages them to determine the average cross-sectional area. Approach B utilizes a weighted average of the area of the 2D photographs to directly compute the average cross-sectional area. A comparison of the accuracy and computational needs of each approach is described as well as preliminary results of an analysis to determine the "optimal" number of images needed for
Recognition methods for 3D textured surfaces
NASA Astrophysics Data System (ADS)
Cula, Oana G.; Dana, Kristin J.
2001-06-01
Texture as a surface representation is the subject of a wide body of computer vision and computer graphics literature. While texture is always associated with a form of repetition in the image, the repeating quantity may vary. The texture may be a color or albedo variation as in a checkerboard, a paisley print or zebra stripes. Very often in real-world scenes, texture is instead due to a surface height variation, e.g. pebbles, gravel, foliage and any rough surface. Such surfaces are referred to here as 3D textured surfaces. Standard texture recognition algorithms are not appropriate for 3D textured surfaces because the appearance of these surfaces changes in a complex manner with viewing direction and illumination direction. Recent methods have been developed for recognition of 3D textured surfaces using a database of surfaces observed under varied imaging parameters. One of these methods is based on 3D textons obtained using K-means clustering of multiscale feature vectors. Another method uses eigen-analysis originally developed for appearance-based object recognition. In this work we develop a hybrid approach that employs both feature grouping and dimensionality reduction. The method is tested using the Columbia-Utrecht texture database and provides excellent recognition rates. The method is compared with existing recognition methods for 3D textured surfaces. A direct comparison is facilitated by empirical recognition rates from the same texture data set. The current method has key advantages over existing methods including requiring less prior information on both the training and novel images.
Analysis of simple 2-D and 3-D metal structures subjected to fragment impact
NASA Technical Reports Server (NTRS)
Witmer, E. A.; Stagliano, T. R.; Spilker, R. L.; Rodal, J. J. A.
1977-01-01
Theoretical methods were developed for predicting the large-deflection elastic-plastic transient structural responses of metal containment or deflector (C/D) structures to cope with rotor burst fragment impact attack. For two-dimensional C/D structures both, finite element and finite difference analysis methods were employed to analyze structural response produced by either prescribed transient loads or fragment impact. For the latter category, two time-wise step-by-step analysis procedures were devised to predict the structural responses resulting from a succession of fragment impacts: the collision force method (CFM) which utilizes an approximate prediction of the force applied to the attacked structure during fragment impact, and the collision imparted velocity method (CIVM) in which the impact-induced velocity increment acquired by a region of the impacted structure near the impact point is computed. The merits and limitations of these approaches are discussed. For the analysis of 3-d responses of C/D structures, only the CIVM approach was investigated.
Puso, M A; Laursen, T A
2002-05-02
Smoothing of contact surfaces can be used to eliminate the chatter typically seen with node on facet contact and give a better representation of the actual contact surface. The latter affect is well demonstrated for problems with interference fits. In this work we present two methods for the smoothing of contact surfaces for 3D finite element contact. In the first method, we employ Gregory patches to smooth the faceted surface in a node on facet implementation. In the second method, we employ a Bezier interpolation of the faceted surface in a mortar method implementation of contact. As is well known, node on facet approaches can exhibit locking due to the failure of the Babuska-Brezzi condition and in some instances fail the patch test. The mortar method implementation is stable and provides optimal convergence in the energy of error. In the this work we demonstrate the superiority of the smoothed versus the non-smoothed node on facet implementations. We also show where the node on facet method fails and some results from the smoothed mortar method implementation.
3D Face Modeling Using the Multi-Deformable Method
Hwang, Jinkyu; Yu, Sunjin; Kim, Joongrock; Lee, Sangyoun
2012-01-01
In this paper, we focus on the problem of the accuracy performance of 3D face modeling techniques using corresponding features in multiple views, which is quite sensitive to feature extraction errors. To solve the problem, we adopt a statistical model-based 3D face modeling approach in a mirror system consisting of two mirrors and a camera. The overall procedure of our 3D facial modeling method has two primary steps: 3D facial shape estimation using a multiple 3D face deformable model and texture mapping using seamless cloning that is a type of gradient-domain blending. To evaluate our method's performance, we generate 3D faces of 30 individuals and then carry out two tests: accuracy test and robustness test. Our method shows not only highly accurate 3D face shape results when compared with the ground truth, but also robustness to feature extraction errors. Moreover, 3D face rendering results intuitively show that our method is more robust to feature extraction errors than other 3D face modeling methods. An additional contribution of our method is that a wide range of face textures can be acquired by the mirror system. By using this texture map, we generate realistic 3D face for individuals at the end of the paper. PMID:23201976
Building virtual 3D bone fragment models to control diaphyseal fracture reduction
NASA Astrophysics Data System (ADS)
Leloup, Thierry; Schuind, Frederic; Lasudry, Nadine; Van Ham, Philippe
1999-05-01
Most fractures of the long bones are displaced and need to be surgically reduced. External fixation is often used but the crucial point of this technique is the control of reduction, which is effected with a brilliance amplifier. This system, giving instantly a x-ray image, has many disadvantages. It implies frequent irradiation to the patient and the surgical team, the visual field is limited, the supplied images are distorted and it only gives 2D information. Consequently, the reduction is occasionally imperfect although intraoperatively it appears acceptable. Using the pains inserted in each fragment as markers and an optical tracker, it is possible to build a virtual 3D model for each principal fragment and to follow its movement during the reduction. This system will supply a 3D image of the fracture in real time and without irradiation. The brilliance amplifier could then be replaced by such a virtual reality system to provide the surgeon with an accurate tool facilitating the reduction of the fracture. The purpose of this work is to show how to build the 3D model for each principal bone fragment.
NASA Astrophysics Data System (ADS)
Ke, Zhipeng; Lu, Tao; Liu, Haichun; Yuan, Haoliang; Ran, Ting; Zhang, Yanmin; Yao, Sihui; Xiong, Xiao; Xu, Jinxing; Xu, Anyang; Chen, Yadong
2014-06-01
Over expression of anaplastic lymphoma kinase (ALK) has been found in many types of cancer, and ALK is a promising therapeutic target for the treatment of cancer. To obtain new potent inhibitors of ALK, we conducted lead optimization using 3D-QSAR modeling and molecular docking investigation of 2,4-diaminopyrimidines and 2,7-disubstituted-pyrrolo[2,1-f][1,2,4]triazine-based compounds. Three favorable 3D-QSAR models (CoMFA with q2, 0.555; r2, 0.939; CoMSIA with q2, 0.625; r2, 0.974; Topomer CoMFA with q2, 0.557; r2 0.756) have been developed to predict the biological activity of novel compounds. Topomer Search was utilized for virtual screening to obtain suitable fragments. The novel compounds generated by molecular fragment replacement (MFR) were evaluated by Topomer CoMFA prediction, Glide (docking) and further evaluated with CoMFA and CoMSIA prediction. 25 novel 2,7-disubstituted-pyrrolo[2,1-f][1,2,4]triazine derivatives as potential ALK inhibitors were finally obtained. In this paper, a combination of CoMFA, CoMSIA and Topomer CoMFA could obtain favorable 3D-QSAR models and suitable fragments for ALK inhibitors optimization. The work flow which comprised 3D-QSAR modeling, Topomer Search, MFR, molecular docking and evaluating criteria could be applied to de novo drug design and the resulted compounds initiate us to further optimize and design new potential ALK inhibitors.
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.
3D scanning modeling method application in ancient city reconstruction
NASA Astrophysics Data System (ADS)
Ren, Pu; Zhou, Mingquan; Du, Guoguang; Shui, Wuyang; Zhou, Pengbo
2015-07-01
With the development of optical engineering technology, the precision of 3D scanning equipment becomes higher, and its role in 3D modeling is getting more distinctive. This paper proposed a 3D scanning modeling method that has been successfully applied in Chinese ancient city reconstruction. On one hand, for the existing architectures, an improved algorithm based on multiple scanning is adopted. Firstly, two pieces of scanning data were rough rigid registered using spherical displacers and vertex clustering method. Secondly, a global weighted ICP (iterative closest points) method is used to achieve a fine rigid registration. On the other hand, for the buildings which have already disappeared, an exemplar-driven algorithm for rapid modeling was proposed. Based on the 3D scanning technology and the historical data, a system approach was proposed for 3D modeling and virtual display of ancient city.
Light Attenuation Method for 3D data acquisition (LAM3D) of bottom particle deposits
NASA Astrophysics Data System (ADS)
Er, Jenn Wei; Law, Adrian W. K.; Adams, E. Eric; Yang, Yang
2015-11-01
We have developed a novel experimental technique, Light Attenuation Method for 3D data acquisition (LAM3D), to acquire three-dimensional spatial characteristics and temporal development of bottom particle deposits. The new technique performs data acquisition with higher spatial and temporal resolution than existing approaches with laser and ultrasonic 3D profilers, and is therefore ideal for laboratory investigations with fast varying changes in the sediment bed, such as the developing deposition profile from sediment clouds commonly formed during dredging or land reclamation projects and the dynamic evolution in movable bed processes in rivers. The principle of the technique is based on the analysis of the light attenuation due to multiple light scattering through the particle deposits layer compared to the clear water column. With appropriate calibration, the particles size and distribution thickness can be quantified by the transmitted light spectrum. In the presentation, we will first show our measurement setup with a light panel for calibrated illumination and a system of DSLR cameras for the light capturing. Subsequently, we shall present the experimental results of fast evolving deposition profile of a barge-disposed sediment cloud upon its bottom impact on the sea bed.
A 3D Level Set Method for Microwave Breast Imaging
Colgan, Timothy J.; Hagness, Susan C.; Van Veen, Barry D.
2015-01-01
Objective Conventional inverse-scattering algorithms for microwave breast imaging result in moderate resolution images with blurred boundaries between tissues. Recent 2D numerical microwave imaging studies demonstrate that the use of a level set method preserves dielectric boundaries, resulting in a more accurate, higher resolution reconstruction of the dielectric properties distribution. Previously proposed level set algorithms are computationally expensive and thus impractical in 3D. In this paper we present a computationally tractable 3D microwave imaging algorithm based on level sets. Methods We reduce the computational cost of the level set method using a Jacobian matrix, rather than an adjoint method, to calculate Frechet derivatives. We demonstrate the feasibility of 3D imaging using simulated array measurements from 3D numerical breast phantoms. We evaluate performance by comparing full 3D reconstructions to those from a conventional microwave imaging technique. We also quantitatively assess the efficacy of our algorithm in evaluating breast density. Results Our reconstructions of 3D numerical breast phantoms improve upon those of a conventional microwave imaging technique. The density estimates from our level set algorithm are more accurate than those of conventional microwave imaging, and the accuracy is greater than that reported for mammographic density estimation. Conclusion Our level set method leads to a feasible level of computational complexity for full 3D imaging, and reconstructs the heterogeneous dielectric properties distribution of the breast more accurately than conventional microwave imaging methods. Significance 3D microwave breast imaging using a level set method is a promising low-cost, non-ionizing alternative to current breast imaging techniques. PMID:26011863
3D face recognition based on a modified ICP method
NASA Astrophysics Data System (ADS)
Zhao, Kankan; Xi, Jiangtao; Yu, Yanguang; Chicharo, Joe F.
2011-11-01
3D face recognition technique has gained much more attention recently, and it is widely used in security system, identification system, and access control system, etc. The core technique in 3D face recognition is to find out the corresponding points in different 3D face images. The classic partial Iterative Closest Point (ICP) method is iteratively align the two point sets based on repetitively calculate the closest points as the corresponding points in each iteration. After several iterations, the corresponding points can be obtained accurately. However, if two 3D face images with different scale are from the same person, the classic partial ICP does not work. In this paper we propose a modified partial Iterative Closest Point (ICP) method in which the scaling effect is considered to achieve 3D face recognition. We design a 3x3 diagonal matrix as the scale matrix in each iteration of the classic partial ICP. The probing face image which is multiplied by the scale matrix will keep the similar scale with the reference face image. Therefore, we can accurately determine the corresponding points even the scales of probing image and reference image are different. 3D face images in our experiments are acquired by a 3D data acquisition system based on Digital Fringe Projection Profilometry (DFPP). A 3D database consists of 30 group images, three images with the same scale, which are from the same person with different views, are included in each group. And in different groups, the scale of the 3 images may be different from other groups. The experiment results show that our proposed method can achieve 3D face recognition, especially in the case that the scales of probing image and referent image are different.
3-D UNSTRUCTURED HEXAHEDRAL-MESH Sn TRANSPORT METHODS
J. MOREL; J. MCGHEE; ET AL
2000-11-01
This is the final report of a three-year, Laboratory-Directed Research and Development (LDRD) project at the Los Alamos National Laboratory (LANL). We have developed a method for solving the neutral-particle transport equation on 3-D unstructured hexahedral meshes using a S{sub n} discretization in angle in conjunction with a discontinuous finite-element discretization in space and a multigroup discretization in energy. Previous methods for solving this equation in 3-D have been limited to rectangular meshes. The unstructured-mesh method that we have developed is far more efficient for solving problems with complex 3-D geometric features than rectangular-mesh methods. In spite of having to make several compromises in our spatial discretization technique and our iterative solution technique, our method has been found to be both accurate and efficient for a broad class of problems.
[An integrated segmentation method for 3D ultrasound carotid artery].
Yang, Xin; Wu, Huihui; Liu, Yang; Xu, Hongwei; Liang, Huageng; Cai, Wenjuan; Fang, Mengjie; Wang, Yujie
2013-07-01
An integrated segmentation method for 3D ultrasound carotid artery was proposed. 3D ultrasound image was sliced into transverse, coronal and sagittal 2D images on the carotid bifurcation point. Then, the three images were processed respectively, and the carotid artery contours and thickness were obtained finally. This paper tries to overcome the disadvantages of current computer aided diagnosis method, such as high computational complexity, easily introduced subjective errors et al. The proposed method could get the carotid artery overall information rapidly, accurately and completely. It could be transplanted into clinical usage for atherosclerosis diagnosis and prevention. PMID:24195385
Improving automated 3D reconstruction methods via vision metrology
NASA Astrophysics Data System (ADS)
Toschi, Isabella; Nocerino, Erica; Hess, Mona; Menna, Fabio; Sargeant, Ben; MacDonald, Lindsay; Remondino, Fabio; Robson, Stuart
2015-05-01
This paper aims to provide a procedure for improving automated 3D reconstruction methods via vision metrology. The 3D reconstruction problem is generally addressed using two different approaches. On the one hand, vision metrology (VM) systems try to accurately derive 3D coordinates of few sparse object points for industrial measurement and inspection applications; on the other, recent dense image matching (DIM) algorithms are designed to produce dense point clouds for surface representations and analyses. This paper strives to demonstrate a step towards narrowing the gap between traditional VM and DIM approaches. Efforts are therefore intended to (i) test the metric performance of the automated photogrammetric 3D reconstruction procedure, (ii) enhance the accuracy of the final results and (iii) obtain statistical indicators of the quality achieved in the orientation step. VM tools are exploited to integrate their main functionalities (centroid measurement, photogrammetric network adjustment, precision assessment, etc.) into the pipeline of 3D dense reconstruction. Finally, geometric analyses and accuracy evaluations are performed on the raw output of the matching (i.e. the point clouds) by adopting a metrological approach. The latter is based on the use of known geometric shapes and quality parameters derived from VDI/VDE guidelines. Tests are carried out by imaging the calibrated Portable Metric Test Object, designed and built at University College London (UCL), UK. It allows assessment of the performance of the image orientation and matching procedures within a typical industrial scenario, characterised by poor texture and known 3D/2D shapes.
Toward Rational Fragment-Based Lead Design without 3D Structures
2012-01-01
Fragment-based lead discovery (FBLD) has become a prime component of the armamentarium of modern drug design programs. FBLD identifies low molecular weight ligands that weakly bind to important biological targets. Three-dimensional structural information about the binding mode is provided by X-ray crystallography or NMR spectroscopy and is subsequently used to improve the lead compounds. Despite tremendous success rates, FBLD relies on the availability of high-resolution structural information, still a bottleneck in drug discovery programs. To overcome these limitations, we recently demonstrated that the meta-structure approach provides an alternative route to rational lead identification in cases where no 3D structure information about the biological target is available. Combined with information-rich NMR data, this strategy provides valuable information for lead development programs. We demonstrate with several examples the feasibility of the combined NMR and meta-structure approach to devise a rational strategy for fragment evolution without resorting to highly resolved protein complex structures. PMID:22889313
MR image denoising method for brain surface 3D modeling
NASA Astrophysics Data System (ADS)
Zhao, De-xin; Liu, Peng-jie; Zhang, De-gan
2014-11-01
Three-dimensional (3D) modeling of medical images is a critical part of surgical simulation. In this paper, we focus on the magnetic resonance (MR) images denoising for brain modeling reconstruction, and exploit a practical solution. We attempt to remove the noise existing in the MR imaging signal and preserve the image characteristics. A wavelet-based adaptive curve shrinkage function is presented in spherical coordinates system. The comparative experiments show that the denoising method can preserve better image details and enhance the coefficients of contours. Using these denoised images, the brain 3D visualization is given through surface triangle mesh model, which demonstrates the effectiveness of the proposed method.
A method to fabricate disconnected silver nanostructures in 3D.
Vora, Kevin; Kang, SeungYeon; Mazur, Eric
2012-01-01
The standard nanofabrication toolkit includes techniques primarily aimed at creating 2D patterns in dielectric media. Creating metal patterns on a submicron scale requires a combination of nanofabrication tools and several material processing steps. For example, steps to create planar metal structures using ultraviolet photolithography and electron-beam lithography can include sample exposure, sample development, metal deposition, and metal liftoff. To create 3D metal structures, the sequence is repeated multiple times. The complexity and difficulty of stacking and aligning multiple layers limits practical implementations of 3D metal structuring using standard nanofabrication tools. Femtosecond-laser direct-writing has emerged as a pre-eminent technique for 3D nanofabrication.(1,2) Femtosecond lasers are frequently used to create 3D patterns in polymers and glasses.(3-7) However, 3D metal direct-writing remains a challenge. Here, we describe a method to fabricate silver nanostructures embedded inside a polymer matrix using a femtosecond laser centered at 800 nm. The method enables the fabrication of patterns not feasible using other techniques, such as 3D arrays of disconnected silver voxels.(8) Disconnected 3D metal patterns are useful for metamaterials where unit cells are not in contact with each other,(9) such as coupled metal dot(10,11)or coupled metal rod(12,13) resonators. Potential applications include negative index metamaterials, invisibility cloaks, and perfect lenses. In femtosecond-laser direct-writing, the laser wavelength is chosen such that photons are not linearly absorbed in the target medium. When the laser pulse duration is compressed to the femtosecond time scale and the radiation is tightly focused inside the target, the extremely high intensity induces nonlinear absorption. Multiple photons are absorbed simultaneously to cause electronic transitions that lead to material modification within the focused region. Using this approach, one can
Novel 3D Compression Methods for Geometry, Connectivity and Texture
NASA Astrophysics Data System (ADS)
Siddeq, M. M.; Rodrigues, M. A.
2016-06-01
A large number of applications in medical visualization, games, engineering design, entertainment, heritage, e-commerce and so on require the transmission of 3D models over the Internet or over local networks. 3D data compression is an important requirement for fast data storage, access and transmission within bandwidth limitations. The Wavefront OBJ (object) file format is commonly used to share models due to its clear simple design. Normally each OBJ file contains a large amount of data (e.g. vertices and triangulated faces, normals, texture coordinates and other parameters) describing the mesh surface. In this paper we introduce a new method to compress geometry, connectivity and texture coordinates by a novel Geometry Minimization Algorithm (GM-Algorithm) in connection with arithmetic coding. First, each vertex ( x, y, z) coordinates are encoded to a single value by the GM-Algorithm. Second, triangle faces are encoded by computing the differences between two adjacent vertex locations, which are compressed by arithmetic coding together with texture coordinates. We demonstrate the method on large data sets achieving compression ratios between 87 and 99 % without reduction in the number of reconstructed vertices and triangle faces. The decompression step is based on a Parallel Fast Matching Search Algorithm (Parallel-FMS) to recover the structure of the 3D mesh. A comparative analysis of compression ratios is provided with a number of commonly used 3D file formats such as VRML, OpenCTM and STL highlighting the performance and effectiveness of the proposed method.
Color dithering methods for LEGO-like 3D printing
NASA Astrophysics Data System (ADS)
Sun, Pei-Li; Sie, Yuping
2015-01-01
Color dithering methods for LEGO-like 3D printing are proposed in this study. The first method is work for opaque color brick building. It is a modification of classic error diffusion. Many color primaries can be chosen. However, RGBYKW is recommended as its image quality is good and the number of color primary is limited. For translucent color bricks, multi-layer color building can enhance the image quality significantly. A LUT-based method is proposed to speed the dithering proceeding and make the color distribution even smoother. Simulation results show the proposed multi-layer dithering method can really improve the image quality of LEGO-like 3D printing.
SAMA: A Method for 3D Morphological Analysis
Cerruti, Florent; Sonnenschein, Carlos; Soto, Ana M.
2016-01-01
Three-dimensional (3D) culture models are critical tools for understanding tissue morphogenesis. A key requirement for their analysis is the ability to reconstruct the tissue into computational models that allow quantitative evaluation of the formed structures. Here, we present Software for Automated Morphological Analysis (SAMA), a method by which epithelial structures grown in 3D cultures can be imaged, reconstructed and analyzed with minimum human intervention. SAMA allows quantitative analysis of key features of epithelial morphogenesis such as ductal elongation, branching and lumen formation that distinguish different hormonal treatments. SAMA is a user-friendly set of customized macros operated via FIJI (http://fiji.sc/Fiji), an open-source image analysis platform in combination with a set of functions in R (http://www.r-project.org/), an open-source program for statistical analysis. SAMA enables a rapid, exhaustive and quantitative 3D analysis of the shape of a population of structures in a 3D image. SAMA is cross-platform, licensed under the GPLv3 and available at http://montevil.theobio.org/content/sama. PMID:27035711
A method for building 3D models of barchan dunes
NASA Astrophysics Data System (ADS)
Nai, Yang; Li-lan, Su; Lin, Wan; Jie, Yang; Shi-yi, Chen; Wei-lu, Hu
2016-01-01
The distributions of barchan dunes are usually represented by digital terrain models (DTMs) overlaid with digital orthophoto maps. Given that most regions with barchan dues have low relief, a 3D map obtained from a DTM may ineffectively show the stereoscopic shape of each dune. The method of building 3D models of barchan dunes using existing modeling software seldom considers the geographical environment. As a result, barchan dune models are often inconsistent with actual DTMs and incompletely express the morphological characteristics of dunes. Manual construction of barchan dune models is also costly and time consuming. Considering these problems, the morphological characteristics of barchan dunes and the mathematical relationships between the morphological parameters of the dunes, such as length, height, and width, are analyzed in this study. The methods of extracting the morphological feature points of barchan dunes, calculating their morphological parameters and building dune outlines and skeleton lines based on the medial axes, are also presented. The dune outlines, skeleton lines, and part of the medial axes of dunes are used to construct a constrained triangulated irregular network. C# and ArcEngine are employed to build 3D models of barchan dunes automatically. Experimental results of a study conducted in Tengger Desert show that the method can be used to approximate the morphological characteristics of barchan dunes and is less time consuming than manual methods.
Breast tumour visualization using 3D quantitative ultrasound methods
NASA Astrophysics Data System (ADS)
Gangeh, Mehrdad J.; Raheem, Abdul; Tadayyon, Hadi; Liu, Simon; Hadizad, Farnoosh; Czarnota, Gregory J.
2016-04-01
Breast cancer is one of the most common cancer types accounting for 29% of all cancer cases. Early detection and treatment has a crucial impact on improving the survival of affected patients. Ultrasound (US) is non-ionizing, portable, inexpensive, and real-time imaging modality for screening and quantifying breast cancer. Due to these attractive attributes, the last decade has witnessed many studies on using quantitative ultrasound (QUS) methods in tissue characterization. However, these studies have mainly been limited to 2-D QUS methods using hand-held US (HHUS) scanners. With the availability of automated breast ultrasound (ABUS) technology, this study is the first to develop 3-D QUS methods for the ABUS visualization of breast tumours. Using an ABUS system, unlike the manual 2-D HHUS device, the whole patient's breast was scanned in an automated manner. The acquired frames were subsequently examined and a region of interest (ROI) was selected in each frame where tumour was identified. Standard 2-D QUS methods were used to compute spectral and backscatter coefficient (BSC) parametric maps on the selected ROIs. Next, the computed 2-D parameters were mapped to a Cartesian 3-D space, interpolated, and rendered to provide a transparent color-coded visualization of the entire breast tumour. Such 3-D visualization can potentially be used for further analysis of the breast tumours in terms of their size and extension. Moreover, the 3-D volumetric scans can be used for tissue characterization and the categorization of breast tumours as benign or malignant by quantifying the computed parametric maps over the whole tumour volume.
Optical Sensors and Methods for Underwater 3D Reconstruction.
Massot-Campos, Miquel; Oliver-Codina, Gabriel
2015-01-01
This paper presents a survey on optical sensors and methods for 3D reconstruction in underwater environments. The techniques to obtain range data have been listed and explained, together with the different sensor hardware that makes them possible. The literature has been reviewed, and a classification has been proposed for the existing solutions. New developments, commercial solutions and previous reviews in this topic have also been gathered and considered. PMID:26694389
Optical Sensors and Methods for Underwater 3D Reconstruction
Massot-Campos, Miquel; Oliver-Codina, Gabriel
2015-01-01
This paper presents a survey on optical sensors and methods for 3D reconstruction in underwater environments. The techniques to obtain range data have been listed and explained, together with the different sensor hardware that makes them possible. The literature has been reviewed, and a classification has been proposed for the existing solutions. New developments, commercial solutions and previous reviews in this topic have also been gathered and considered. PMID:26694389
Discrete Method of Images for 3D Radio Propagation Modeling
NASA Astrophysics Data System (ADS)
Novak, Roman
2016-09-01
Discretization by rasterization is introduced into the method of images (MI) in the context of 3D deterministic radio propagation modeling as a way to exploit spatial coherence of electromagnetic propagation for fine-grained parallelism. Traditional algebraic treatment of bounding regions and surfaces is replaced by computer graphics rendering of 3D reflections and double refractions while building the image tree. The visibility of reception points and surfaces is also resolved by shader programs. The proposed rasterization is shown to be of comparable run time to that of the fundamentally parallel shooting and bouncing rays. The rasterization does not affect the signal evaluation backtracking step, thus preserving its advantage over the brute force ray-tracing methods in terms of accuracy. Moreover, the rendering resolution may be scaled back for a given level of scenario detail with only marginal impact on the image tree size. This allows selection of scene optimized execution parameters for faster execution, giving the method a competitive edge. The proposed variant of MI can be run on any GPU that supports real-time 3D graphics.
Parallel 3D Mortar Element Method for Adaptive Nonconforming Meshes
NASA Technical Reports Server (NTRS)
Feng, Huiyu; Mavriplis, Catherine; VanderWijngaart, Rob; Biswas, Rupak
2004-01-01
High order methods are frequently used in computational simulation for their high accuracy. An efficient way to avoid unnecessary computation in smooth regions of the solution is to use adaptive meshes which employ fine grids only in areas where they are needed. Nonconforming spectral elements allow the grid to be flexibly adjusted to satisfy the computational accuracy requirements. The method is suitable for computational simulations of unsteady problems with very disparate length scales or unsteady moving features, such as heat transfer, fluid dynamics or flame combustion. In this work, we select the Mark Element Method (MEM) to handle the non-conforming interfaces between elements. A new technique is introduced to efficiently implement MEM in 3-D nonconforming meshes. By introducing an "intermediate mortar", the proposed method decomposes the projection between 3-D elements and mortars into two steps. In each step, projection matrices derived in 2-D are used. The two-step method avoids explicitly forming/deriving large projection matrices for 3-D meshes, and also helps to simplify the implementation. This new technique can be used for both h- and p-type adaptation. This method is applied to an unsteady 3-D moving heat source problem. With our new MEM implementation, mesh adaptation is able to efficiently refine the grid near the heat source and coarsen the grid once the heat source passes. The savings in computational work resulting from the dynamic mesh adaptation is demonstrated by the reduction of the the number of elements used and CPU time spent. MEM and mesh adaptation, respectively, bring irregularity and dynamics to the computer memory access pattern. Hence, they provide a good way to gauge the performance of computer systems when running scientific applications whose memory access patterns are irregular and unpredictable. We select a 3-D moving heat source problem as the Unstructured Adaptive (UA) grid benchmark, a new component of the NAS Parallel
NASA Astrophysics Data System (ADS)
Sharkawi, K.-H.; Abdul-Rahman, A.
2013-09-01
to LoD4. The accuracy and structural complexity of the 3D objects increases with the LoD level where LoD0 is the simplest LoD (2.5D; Digital Terrain Model (DTM) + building or roof print) while LoD4 is the most complex LoD (architectural details with interior structures). Semantic information is one of the main components in CityGML and 3D City Models, and provides important information for any analyses. However, more often than not, the semantic information is not available for the 3D city model due to the unstandardized modelling process. One of the examples is where a building is normally generated as one object (without specific feature layers such as Roof, Ground floor, Level 1, Level 2, Block A, Block B, etc). This research attempts to develop a method to improve the semantic data updating process by segmenting the 3D building into simpler parts which will make it easier for the users to select and update the semantic information. The methodology is implemented for 3D buildings in LoD2 where the buildings are generated without architectural details but with distinct roof structures. This paper also introduces hybrid semantic-geometric 3D segmentation method that deals with hierarchical segmentation of a 3D building based on its semantic value and surface characteristics, fitted by one of the predefined primitives. For future work, the segmentation method will be implemented as part of the change detection module that can detect any changes on the 3D buildings, store and retrieve semantic information of the changed structure, automatically updates the 3D models and visualize the results in a userfriendly graphical user interface (GUI).
System and method for 3D printing of aerogels
Worsley, Marcus A.; Duoss, Eric; Kuntz, Joshua; Spadaccini, Christopher; Zhu, Cheng
2016-03-08
A method of forming an aerogel. The method may involve providing a graphene oxide powder and mixing the graphene oxide powder with a solution to form an ink. A 3D printing technique may be used to write the ink into a catalytic solution that is contained in a fluid containment member to form a wet part. The wet part may then be cured in a sealed container for a predetermined period of time at a predetermined temperature. The cured wet part may then be dried to form a finished aerogel part.
Method and simulation to study 3D crosstalk perception
NASA Astrophysics Data System (ADS)
Khaustova, Dar'ya; Blondé, Laurent; Huynh-Thu, Quan; Vienne, Cyril; Doyen, Didier
2012-03-01
To various degrees, all modern 3DTV displays suffer from crosstalk, which can lead to a decrease of both visual quality and visual comfort, and also affect perception of depth. In the absence of a perfect 3D display technology, crosstalk has to be taken into account when studying perception of 3D stereoscopic content. In order to improve 3D presentation systems and understand how to efficiently eliminate crosstalk, it is necessary to understand its impact on human perception. In this paper, we present a practical method to study the perception of crosstalk. The approach consists of four steps: (1) physical measurements of a 3DTV, (2) building of a crosstalk surface based on those measurements and representing specifically the behavior of that 3TV, (3) manipulation of the crosstalk function and application on reference images to produce test images degraded by crosstalk in various ways, and (4) psychophysical tests. Our approach allows both a realistic representation of the behavior of a 3DTV and the easy manipulation of its resulting crosstalk in order to conduct psycho-visual experiments. Our approach can be used in all studies requiring the understanding of how crosstalk affects perception of stereoscopic content and how it can be corrected efficiently.
Davis, Holly J; Kavanagh, Madeline E; Balan, Tudor; Abell, Chris; Coyne, Anthony G
2016-08-01
The search for new scaffolds to complement current HTS and fragment libraries is an active area of research. The development of novel strategies to synthesise compounds with 3D character in order to expand the diversity of a fragment library was explored. A range of substituted bicyclo[2,2,1]spirooxindoles were synthesised using a Diels-Alder [4+2] cycloaddition reaction. Both diastereoisomers were isolated from the reactions and these 3D fragment scaffolds were screened against the cytochrome P450 enzyme CYP121 from Mycobacterium tuberculosis. A number of hits were identified to bind to CYP121 and were shown to exhibit Type I binding interactions with the heme group. PMID:27287372
Modeling of constructional elements fragmentation:3-D statement and probabilistic approach
NASA Astrophysics Data System (ADS)
Gerasimov, Alexander; Pashkov, Sergey
2011-06-01
The heterogeneity of real materials structure influencing on distribution of material characteristics is one of the factors determining a destruction character. The introduction of the given factor in the equations of mechanics of deformed solid is possible at the use of probabilistic laws of characteristics distribution in the volume of a considered design. The explosive fragmentation of the open and closed shells, thick plate punching by HE charged shell on a normal and at an angle, plate and a shell fragmentation after plate piercing and under HE charge explosion, thin barrier punching on a normal and at an angle, crushing of metal rings dressed on a copper tube, process of high-speed impact of laminated-spaced metallic plates with steel spheres modeling debris of space bodies and artificial objects are considered. The processes are calculated in view of material heterogeneity. To calculate elastoplastic flows and detonation products we used the technique realized on tetrahedral cells and based on Wilkins method for calculation of internal points of a body and on Johnson method for calculation of contact interactions.
Single-camera fixed perspective 360-deg 3D method
NASA Astrophysics Data System (ADS)
Harding, Kevin G.; Fergan, Robert K.
1997-01-01
The use of 3D methods for such applications as feature locations within a wide field-of-view, such as for automated guided vehicles or large assembly work, offers some distinct challenges. The use of stereo viewing has often been the method of choice due to the wide area coverage and hardware simplicity. However, stereo based methods suffer from a loss of spatial position resolution for more distant object as compared to close objects due to the high demagnification needed to cover large fields-of-view. A long depth-of-field in such systems may also degrade the general ability to perform correlations due to poor focus. In addition, stereo looses distance resolution for features nearing the line of the two cameras, typically requiring movement of the cameras. The paper presents a novel method of obtaining 3D scene information as seen from the center of a cylindrical field. The method described uses a single camera with a view that is rotated through 360 degrees by means of a continuously rotating mirror. The viewing systems uses a constant field of view optical system that provides a constant X-Y resolution of features in the scene over depths of several meters. Comparing successive images with the readout from an encoder on the rotating mirror generates all locations of objects within a limited height cylinder. This paper will discuss the sources of errors and typical capabilities of this approach in light of a real-time part location tracking application useful in assembly systems.
3D reconstruction methods of coronal structures by radio observations
NASA Astrophysics Data System (ADS)
Aschwanden, Markus J.; Bastian, T. S.; White, Stephen M.
1992-11-01
The ability to carry out the three dimensional (3D) reconstruction of structures in the solar corona would represent a major advance in the study of the physical properties in active regions and in flares. Methods which allow a geometric reconstruction of quasistationary coronal structures (for example active region loops) or dynamic structures (for example flaring loops) are described: stereoscopy of multi-day imaging observations by the VLA (Very Large Array); tomography of optically thin emission (in radio or soft x-rays); multifrequency band imaging by the VLA; and tracing of magnetic field lines by propagating electron beams.
3D reconstruction methods of coronal structures by radio observations
NASA Technical Reports Server (NTRS)
Aschwanden, Markus J.; Bastian, T. S.; White, Stephen M.
1992-01-01
The ability to carry out the three dimensional (3D) reconstruction of structures in the solar corona would represent a major advance in the study of the physical properties in active regions and in flares. Methods which allow a geometric reconstruction of quasistationary coronal structures (for example active region loops) or dynamic structures (for example flaring loops) are described: stereoscopy of multi-day imaging observations by the VLA (Very Large Array); tomography of optically thin emission (in radio or soft x-rays); multifrequency band imaging by the VLA; and tracing of magnetic field lines by propagating electron beams.
A perceptual preprocess method for 3D-HEVC
NASA Astrophysics Data System (ADS)
Shi, Yawen; Wang, Yongfang; Wang, Yubing
2015-08-01
A perceptual preprocessing method for 3D-HEVC coding is proposed in the paper. Firstly we proposed a new JND model, which accounts for luminance contrast masking effect, spatial masking effect, and temporal masking effect, saliency characteristic as well as depth information. We utilize spectral residual approach to obtain the saliency map and built a visual saliency factor based on saliency map. In order to distinguish the sensitivity of objects in different depth. We segment each texture frame into foreground and background by a automatic threshold selection algorithm using corresponding depth information, and then built a depth weighting factor. A JND modulation factor is built with a linear combined with visual saliency factor and depth weighting factor to adjust the JND threshold. Then, we applied the proposed JND model to 3D-HEVC for residual filtering and distortion coefficient processing. The filtering process is that the residual value will be set to zero if the JND threshold is greater than residual value, or directly subtract the JND threshold from residual value if JND threshold is less than residual value. Experiment results demonstrate that the proposed method can achieve average bit rate reduction of 15.11%, compared to the original coding scheme with HTM12.1, while maintains the same subjective quality.
A method of PSF generation for 3D brightfield deconvolution.
Tadrous, P J
2010-02-01
This paper addresses the problem of 3D deconvolution of through focus widefield microscope datasets (Z-stacks). One of the most difficult stages in brightfield deconvolution is finding the point spread function. A theoretically calculated point spread function (called a 'synthetic PSF' in this paper) requires foreknowledge of many system parameters and still gives only approximate results. A point spread function measured from a sub-resolution bead suffers from low signal-to-noise ratio, compounded in the brightfield setting (by contrast to fluorescence) by absorptive, refractive and dispersal effects. This paper describes a method of point spread function estimation based on measurements of a Z-stack through a thin sample. This Z-stack is deconvolved by an idealized point spread function derived from the same Z-stack to yield a point spread function of high signal-to-noise ratio that is also inherently tailored to the imaging system. The theory is validated by a practical experiment comparing the non-blind 3D deconvolution of the yeast Saccharomyces cerevisiae with the point spread function generated using the method presented in this paper (called the 'extracted PSF') to a synthetic point spread function. Restoration of both high- and low-contrast brightfield structures is achieved with fewer artefacts using the extracted point spread function obtained with this method. Furthermore the deconvolution progresses further (more iterations are allowed before the error function reaches its nadir) with the extracted point spread function compared to the synthetic point spread function indicating that the extracted point spread function is a better fit to the brightfield deconvolution model than the synthetic point spread function. PMID:20096049
The COMET method in 3-D hexagonal geometry
Connolly, K. J.; Rahnema, F.
2012-07-01
The hybrid stochastic-deterministic coarse mesh radiation transport (COMET) method developed at Georgia Tech now solves reactor core problems in 3-D hexagonal geometry. In this paper, the method is used to solve three preliminary test problems designed to challenge the method with steep flux gradients, high leakage, and strong asymmetry and heterogeneity in the core. The test problems are composed of blocks taken from a high temperature test reactor benchmark problem. As the method is still in development, these problems and their results are strictly preliminary. Results are compared to whole core Monte Carlo reference solutions in order to verify the method. Relative errors are on the order of 50 pcm in core eigenvalue, and mean relative error in pin fission density calculations is less than 1% in these difficult test cores. The method requires the one-time pre-computation of a response expansion coefficient library, which may be compiled in a comparable amount of time to a single whole core Monte Carlo calculation. After the library has been computed, COMET may solve any number of core configurations on the order of an hour, representing a significant gain in efficiency over other methods for whole core transport calculations. (authors)
The 3D inelastic analysis methods for hot section components
NASA Technical Reports Server (NTRS)
Dame, L. T.; Mcknight, R. L.
1983-01-01
The objective of this research is to develop an analytical tool capable of economically evaluating the cyclic time dependent plasticity which occurs in hot section engine components in areas of strain concentration resulting from the combination of both mechanical and thermal stresses. The techniques developed must be capable of accommodating large excursions in temperatures with the associated variations in material properties including plasticity and creep. The overall objective of this proposed program is to develop advanced 3-D inelastic structural/stress analysis methods and solution strategies for more accurate and yet more cost effective analysis of combustors, turbine blades, and vanes. The approach will be to develop four different theories, one linear and three higher order with increasing complexities including embedded singularities.
On 3D inelastic analysis methods for hot section components
NASA Technical Reports Server (NTRS)
Mcknight, R. L.; Chen, P. C.; Dame, L. T.; Holt, R. V.; Huang, H.; Hartle, M.; Gellin, S.; Allen, D. H.; Haisler, W. E.
1986-01-01
Accomplishments are described for the 2-year program, to develop advanced 3-D inelastic structural stress analysis methods and solution strategies for more accurate and cost effective analysis of combustors, turbine blades and vanes. The approach was to develop a matrix of formulation elements and constitutive models. Three constitutive models were developed in conjunction with optimized iterating techniques, accelerators, and convergence criteria within a framework of dynamic time incrementing. Three formulations models were developed; an eight-noded mid-surface shell element, a nine-noded mid-surface shell element and a twenty-noded isoparametric solid element. A separate computer program was developed for each combination of constitutive model-formulation model. Each program provides a functional stand alone capability for performing cyclic nonlinear structural analysis. In addition, the analysis capabilities incorporated into each program can be abstracted in subroutine form for incorporation into other codes or to form new combinations.
The 3D inelastic analysis methods for hot section components
NASA Technical Reports Server (NTRS)
Mcknight, R. L.; Maffeo, R. J.; Tipton, M. T.; Weber, G.
1992-01-01
A two-year program to develop advanced 3D inelastic structural stress analysis methods and solution strategies for more accurate and cost effective analysis of combustors, turbine blades, and vanes is described. The approach was to develop a matrix of formulation elements and constitutive models. Three constitutive models were developed in conjunction with optimized iterating techniques, accelerators, and convergence criteria within a framework of dynamic time incrementing. Three formulation models were developed: an eight-noded midsurface shell element; a nine-noded midsurface shell element; and a twenty-noded isoparametric solid element. A separate computer program has been developed for each combination of constitutive model-formulation model. Each program provides a functional stand alone capability for performing cyclic nonlinear structural analysis. In addition, the analysis capabilities incorporated into each program can be abstracted in subroutine form for incorporation into other codes or to form new combinations.
Lattice Boltzmann Method for 3-D Flows with Curved Boundary
NASA Technical Reports Server (NTRS)
Mei, Renwei; Shyy, Wei; Yu, Dazhi; Luo, Li-Shi
2002-01-01
In this work, we investigate two issues that are important to computational efficiency and reliability in fluid dynamics applications of the lattice, Boltzmann equation (LBE): (1) Computational stability and accuracy of different lattice Boltzmann models and (2) the treatment of the boundary conditions on curved solid boundaries and their 3-D implementations. Three athermal 3-D LBE models (D3QI5, D3Ql9, and D3Q27) are studied and compared in terms of efficiency, accuracy, and robustness. The boundary treatment recently developed by Filippova and Hanel and Met et al. in 2-D is extended to and implemented for 3-D. The convergence, stability, and computational efficiency of the 3-D LBE models with the boundary treatment for curved boundaries were tested in simulations of four 3-D flows: (1) Fully developed flows in a square duct, (2) flow in a 3-D lid-driven cavity, (3) fully developed flows in a circular pipe, and (4) a uniform flow over a sphere. We found that while the fifteen-velocity 3-D (D3Ql5) model is more prone to numerical instability and the D3Q27 is more computationally intensive, the 63Q19 model provides a balance between computational reliability and efficiency. Through numerical simulations, we demonstrated that the boundary treatment for 3-D arbitrary curved geometry has second-order accuracy and possesses satisfactory stability characteristics.
3D Wavelet-Based Filter and Method
Moss, William C.; Haase, Sebastian; Sedat, John W.
2008-08-12
A 3D wavelet-based filter for visualizing and locating structural features of a user-specified linear size in 2D or 3D image data. The only input parameter is a characteristic linear size of the feature of interest, and the filter output contains only those regions that are correlated with the characteristic size, thus denoising the image.
Methods for Geometric Data Validation of 3d City Models
NASA Astrophysics Data System (ADS)
Wagner, D.; Alam, N.; Wewetzer, M.; Pries, M.; Coors, V.
2015-12-01
Geometric quality of 3D city models is crucial for data analysis and simulation tasks, which are part of modern applications of the data (e.g. potential heating energy consumption of city quarters, solar potential, etc.). Geometric quality in these contexts is however a different concept as it is for 2D maps. In the latter case, aspects such as positional or temporal accuracy and correctness represent typical quality metrics of the data. They are defined in ISO 19157 and should be mentioned as part of the metadata. 3D data has a far wider range of aspects which influence their quality, plus the idea of quality itself is application dependent. Thus, concepts for definition of quality are needed, including methods to validate these definitions. Quality on this sense means internal validation and detection of inconsistent or wrong geometry according to a predefined set of rules. A useful starting point would be to have correct geometry in accordance with ISO 19107. A valid solid should consist of planar faces which touch their neighbours exclusively in defined corner points and edges. No gaps between them are allowed, and the whole feature must be 2-manifold. In this paper, we present methods to validate common geometric requirements for building geometry. Different checks based on several algorithms have been implemented to validate a set of rules derived from the solid definition mentioned above (e.g. water tightness of the solid or planarity of its polygons), as they were developed for the software tool CityDoctor. The method of each check is specified, with a special focus on the discussion of tolerance values where they are necessary. The checks include polygon level checks to validate the correctness of each polygon, i.e. closeness of the bounding linear ring and planarity. On the solid level, which is only validated if the polygons have passed validation, correct polygon orientation is checked, after self-intersections outside of defined corner points and edges
Two new methods for simulating photolithography development in 3D
Helmsen, J.; Colella, P.; Dorr, M.; Puckett, E.G.
1997-01-30
Two methods are presented for simulating the development of photolithographic profiles during the resist dissolution phase. These algorithms are the volume-of-fluid algorithm, and the steady level-set algorithm. They are compared with the ray-trace, cell, and level-set techniques employed in SAMPLE-3D. The volume-of-fluid algorithm employs an Euclidean Grid with volume fractions. At each time step, the surface is reconstructed by computing an approximation of the tangent plane of the surface in each cell that contains a value between 0 and 1. The geometry constructed in this manner is used to determine flow velocity vectors and the flux across each edge. The material is then advanced by a split advection scheme. The steady Level Set algorithm is an extension of the Iterative Level Set algorithm. The steady Level Set algorithm combines Fast Level Set concepts and a technique for finding zero residual solutions to the ( ) function. The etch time for each cell is calculated in a time ordered manner. Use of heap sorting data structures allows the algorithm to execute extremely quickly. Comparisons of the methods have been performed and results shown.
New method for 3D reconstruction in digital tomosynthesis
NASA Astrophysics Data System (ADS)
Claus, Bernhard E. H.; Eberhard, Jeffrey W.
2002-05-01
Digital tomosynthesis mammography is an advanced x-ray application that can provide detailed 3D information about the imaged breast. We introduce a novel reconstruction method based on simple backprojection, which yields high contrast reconstructions with reduced artifacts at a relatively low computational complexity. The first step in the proposed reconstruction method is a simple backprojection with an order statistics-based operator (e.g., minimum) used for combining the backprojected images into a reconstructed slice. Accordingly, a given pixel value does generally not contribute to all slices. The percentage of slices where a given pixel value does not contribute, as well as the associated reconstructed values, are collected. Using a form of re-projection consistency constraint, one now updates the projection images, and repeats the order statistics backprojection reconstruction step, but now using the enhanced projection images calculated in the first step. In our digital mammography application, this new approach enhances the contrast of structures in the reconstruction, and allows in particular to recover the loss in signal level due to reduced tissue thickness near the skinline, while keeping artifacts to a minimum. We present results obtained with the algorithm for phantom images.
The 3D Shape of the Dendrite by WDT Method
NASA Astrophysics Data System (ADS)
Tang, Chao; Mitobe, Kazutaka; Yoshimura, Noboru
The purpose of this study is use of a three dimension (3D) measuring system that can automatically measure surface condition. We applied the WDT method that is one of the migration acceleration testing methods, to calculate the spatial variation of the electrodes of ion immigration on a glass epoxy printed wiring board. We also investigated the spatial shape and its variation of dendrite after short circuit for the cases of uniform and nonuniform field strength. As a result the phenomenon of immigration peak of separated matter from cathode to anode due to nonuniform was reported.The moving of the peak of the separated matter is supposed to be due to Cu(OH)2's change in accumulation status. Under the nonuniform and uniform situation, the behavior of separated matter will change after occurring short circuit between the electrodes. Therefore in order to avoid the progress of ion immigration, it is necessary to pay attention to the field strength in hardwiring and the curvature so that the field strength of the wiring pattern cannot be very high.
A support-operator method for 3-D rupture dynamics
NASA Astrophysics Data System (ADS)
Ely, Geoffrey P.; Day, Steven M.; Minster, Jean-Bernard
2009-06-01
We present a numerical method to simulate spontaneous shear crack propagation within a heterogeneous, 3-D, viscoelastic medium. Wave motions are computed on a logically rectangular hexahedral mesh, using the generalized finite-difference method of Support Operators (SOM). This approach enables modelling of non-planar surfaces and non-planar fault ruptures. Our implementation, the Support Operator Rupture Dynamics (SORD) code, is highly scalable, enabling large-scale, multiprocessors calculations. The fault surface is modelled by coupled double nodes, where rupture occurs as dictated by the local stress conditions and a frictional failure law. The method successfully performs test problems developed for the Southern California Earthquake Center (SCEC)/U.S. Geological Survey (USGS) dynamic earthquake rupture code validation exercise, showing good agreement with semi-analytical boundary integral method results. We undertake further dynamic rupture tests to quantify numerical errors introduced by shear deformations to the hexahedral mesh. We generate a family of meshes distorted by simple shearing, in the along-strike direction, up to a maximum of 73°. For SCEC/USGS validation problem number 3, grid-induced errors increase with mesh shear angle, with the logarithm of error approximately proportional to angle over the range tested. At 73°, rms misfits are about 10 per cent for peak slip rate, and 0.5 per cent for both rupture time and total slip, indicating that the method (which, up to now, we have applied mainly to near-vertical strike-slip faulting) is also capable of handling geometries appropriate to low-angle surface-rupturing thrust earthquakes. Additionally, we demonstrate non-planar rupture effects, by modifying the test geometry to include, respectively, cylindrical curvature and sharp kinks.
NASA Astrophysics Data System (ADS)
Lague, D.; Brodu, N.; Leroux, J.
2012-12-01
Ground based lidar and photogrammetric techniques are increasingly used to track the evolution of natural surfaces in 3D at an unprecedented resolution and precision. The range of applications encompass many type of natural surfaces with different geometries and roughness characteristics (landslides, cliff erosion, river beds, bank erosion,....). Unravelling surface change in these contexts requires to compare large point clouds in 2D or 3D. The most commonly used method in geomorphology is based on a 2D difference of the gridded point clouds. Yet this is hardly adapted to many 3D natural environments such as rivers (with horizontal beds and vertical banks), while gridding complex rough surfaces is a complex task. On the other hand, tools allowing to perform 3D comparison are scarce and may require to mesh the point clouds which is difficult on rough natural surfaces. Moreover, existing 3D comparison tools do not provide an explicit calculation of confidence intervals that would factor in registration errors, roughness effects and instrument related position uncertainties. To unlock this problem, we developed the first algorithm combining a 3D measurement of surface change directly on point clouds with an estimate of spatially variable confidence intervals (called M3C2). The method has two steps : (1) surface normal estimation and orientation in 3D at a scale consistent with the local roughness ; (2) measurement of mean surface change along the normal direction with explicit calculation of a local confidence interval. Comparison with existing 3D methods based on a closest-point calculation demonstrates the higher precision of the M3C2 method when mm changes needs to be detected. The M3C2 method is also simple to use as it does not require surface meshing or gridding, and is not sensitive to missing data or change in point density. We also present a 3D classification tool (CANUPO) for vegetation removal based on a new geometrical measure: the multi
NASA Astrophysics Data System (ADS)
Commerçon, B.; Hennebelle, P.; Levrier, F.; Launhardt, R.; Henning, Th.
2012-03-01
I will present radiation-magneto-hydrodynamics calculations of low-mass and massive dense core collapse, focusing on the first collapse and the first hydrostatic core (first Larson core) formation. The influence of magnetic field and initial mass on the fragmentation properties will be investigated. In the first part reporting low mass dense core collapse calculations, synthetic observations of spectral energy distributions will be derived, as well as classical observational quantities such as bolometric temperature and luminosity. I will show how the dust continuum can help to target first hydrostatic cores and to state about the nature of VeLLOs. Last, I will present synthetic ALMA observation predictions of first hydrostatic cores which may give an answer, if not definitive, to the fragmentation issue at the early Class 0 stage. In the second part, I will report the results of radiation-magneto-hydrodynamics calculations in the context of high mass star formation, using for the first time a self-consistent model for photon emission (i.e. via thermal emission and in radiative shocks) and with the high resolution necessary to resolve properly magnetic braking effects and radiative shocks on scales <100 AU (Commercon, Hennebelle & Henning ApJL 2011). In this study, we investigate the combined effects of magnetic field, turbulence, and radiative transfer on the early phases of the collapse and the fragmentation of massive dense cores (M=100 M_⊙). We identify a new mechanism that inhibits initial fragmentation of massive dense cores, where magnetic field and radiative transfer interplay. We show that this interplay becomes stronger as the magnetic field strength increases. We speculate that highly magnetized massive dense cores are good candidates for isolated massive star formation, while moderately magnetized massive dense cores are more appropriate to form OB associations or small star clusters. Finally we will also present synthetic observations of these
3-D dynamic rupture simulations by a finite volume method
NASA Astrophysics Data System (ADS)
Benjemaa, M.; Glinsky-Olivier, N.; Cruz-Atienza, V. M.; Virieux, J.
2009-07-01
Dynamic rupture of a 3-D spontaneous crack of arbitrary shape is investigated using a finite volume (FV) approach. The full domain is decomposed in tetrahedra whereas the surface, on which the rupture takes place, is discretized with triangles that are faces of tetrahedra. First of all, the elastodynamic equations are described into a pseudo-conservative form for an easy application of the FV discretization. Explicit boundary conditions are given using criteria based on the conservation of discrete energy through the crack surface. Using a stress-threshold criterion, these conditions specify fluxes through those triangles that have suffered rupture. On these broken surfaces, stress follows a linear slip-weakening law, although other friction laws can be implemented. For The Problem Version 3 of the dynamic-rupture code verification exercise conducted by the SCEC/USGS, numerical solutions on a planar fault exhibit a very high convergence rate and are in good agreement with the reference one provided by a finite difference (FD) technique. For a non-planar fault of parabolic shape, numerical solutions agree satisfactorily well with those obtained with a semi-analytical boundary integral method in terms of shear stress amplitudes, stopping phases arrival times and stress overshoots. Differences between solutions are attributed to the low-order interpolation of the FV approach, whose results are particularly sensitive to the mesh regularity (structured/unstructured). We expect this method, which is well adapted for multiprocessor parallel computing, to be competitive with others for solving large scale dynamic ruptures scenarios of seismic sources in the near future.
Methods For Electronic 3-D Moving Pictures Without Glasses
NASA Astrophysics Data System (ADS)
Collender, Robert B.
1987-06-01
This paper describes implementation approaches in image acquisition and playback for 3-D computer graphics, 3-D television and 3-D theatre movies without special glasses. Projection lamps, spatial light modulators, CRT's and dynamic scanning are all eliminated by the application of an active image array, all static components and a semi-specular screen. The resulting picture shows horizontal parallax with a wide horizontal view field (up to 360 de-grees) giving a holographic appearance in full color with smooth continuous viewing without speckle. Static component systems are compared with dynamic component systems using both linear and circular arrays. Implementation of computer graphic systems are shown that allow complex shaded color images to extend from the viewer's eyes to infinity. Large screen systems visible by hundreds of people are feasible by the use of low f-stops and high gain screens in projection. Screen geometries and special screen properties are shown. Viewing characteristics offer no restrictions in view-position over the entire view-field and have a "look-around" feature for all the categories of computer graphics, television and movies. Standard video cassettes and optical discs can also interface the system to generate a 3-D window viewable without glasses. A prognosis is given for technology application to 3-D pictures without glasses that replicate the daily viewing experience. Super-position of computer graphics on real-world pictures is shown feasible.
Pak, Tannaz; Butler, Ian B.; Geiger, Sebastian; van Dijke, Marinus I. J.; Sorbie, Ken S.
2015-01-01
Using X-ray computed microtomography, we have visualized and quantified the in situ structure of a trapped nonwetting phase (oil) in a highly heterogeneous carbonate rock after injecting a wetting phase (brine) at low and high capillary numbers. We imaged the process of capillary desaturation in 3D and demonstrated its impacts on the trapped nonwetting phase cluster size distribution. We have identified a previously unidentified pore-scale event during capillary desaturation. This pore-scale event, described as droplet fragmentation of the nonwetting phase, occurs in larger pores. It increases volumetric production of the nonwetting phase after capillary trapping and enlarges the fluid−fluid interface, which can enhance mass transfer between the phases. Droplet fragmentation therefore has implications for a range of multiphase flow processes in natural and engineered porous media with complex heterogeneous pore spaces. PMID:25646491
X3D moving grid methods for semiconductor applications
Kuprat, A.; Cartwright, D.; Gammel, J.T.; George, D.; Kendrick, B.; Kilcrease, D.; Trease, H.; Walker, R.
1997-11-01
The Los Alamos 3D grid toolbox handles grid maintenance chores and provides access to a sophisticated set of optimization algorithms for unstructured grids. The application of these tools to semiconductor problems is illustrated in three examples: grain growth, topographic deposition and electrostatics. These examples demonstrate adaptive smoothing, front tracking, and automatic, adaptive refinement/derefinement.
Object-oriented urban 3D spatial data model organization method
NASA Astrophysics Data System (ADS)
Li, Jing-wen; Li, Wen-qing; Lv, Nan; Su, Tao
2015-12-01
This paper combined the 3d data model with object-oriented organization method, put forward the model of 3d data based on object-oriented method, implemented the city 3d model to quickly build logical semantic expression and model, solved the city 3d spatial information representation problem of the same location with multiple property and the same property with multiple locations, designed the space object structure of point, line, polygon, body for city of 3d spatial database, and provided a new thought and method for the city 3d GIS model and organization management.
3D range scan enhancement using image-based methods
NASA Astrophysics Data System (ADS)
Herbort, Steffen; Gerken, Britta; Schugk, Daniel; Wöhler, Christian
2013-10-01
This paper addresses the problem of 3D surface scan refinement, which is desirable due to noise, outliers, and missing measurements being present in the 3D surfaces obtained with a laser scanner. We present a novel algorithm for the fusion of absolute laser scanner depth profiles and photometrically estimated surface normal data, which yields a noise-reduced and highly detailed depth profile with large scale shape robustness. In contrast to other approaches published in the literature, the presented algorithm (1) regards non-Lambertian surfaces, (2) simultaneously computes surface reflectance (i.e. BRDF) parameters required for 3D reconstruction, (3) models pixelwise incident light and viewing directions, and (4) accounts for interreflections. The algorithm as such relies on the minimization of a three-component error term, which penalizes intensity deviations, integrability deviations, and deviations from the known large-scale surface shape. The solution of the error minimization is obtained iteratively based on a calculus of variations. BRDF parameters are estimated by initially reducing and then iteratively refining the optical resolution, which provides the required robust data basis. The 3D reconstruction of concave surface regions affected by interreflections is improved by compensating global illumination in the image data. The algorithm is evaluated based on eight objects with varying albedos and reflectance behaviors (diffuse, specular, metallic). The qualitative evaluation shows a removal of outliers and a strong reduction of noise, while the large scale shape is preserved. Fine surface details Which are previously not contained in the surface scans, are incorporated through using image data. The algorithm is evaluated with respect to its absolute accuracy using two caliper objects of known shape, and based on synthetically generated data. The beneficial effect of interreflection compensation on the reconstruction accuracy is evaluated quantitatively in a
Calibration Methods for a 3D Triangulation Based Camera
NASA Astrophysics Data System (ADS)
Schulz, Ulrike; Böhnke, Kay
A sensor in a camera takes a gray level image (1536 x 512 pixels), which is reflected by a reference body. The reference body is illuminated by a linear laser line. This gray level image can be used for a 3D calibration. The following paper describes how a calibration program calculates the calibration factors. The calibration factors serve to determine the size of an unknown reference body.
Parallel 3-D method of characteristics in MPACT
Kochunas, B.; Dovvnar, T. J.; Liu, Z.
2013-07-01
A new parallel 3-D MOC kernel has been developed and implemented in MPACT which makes use of the modular ray tracing technique to reduce computational requirements and to facilitate parallel decomposition. The parallel model makes use of both distributed and shared memory parallelism which are implemented with the MPI and OpenMP standards, respectively. The kernel is capable of parallel decomposition of problems in space, angle, and by characteristic rays up to 0(104) processors. Initial verification of the parallel 3-D MOC kernel was performed using the Takeda 3-D transport benchmark problems. The eigenvalues computed by MPACT are within the statistical uncertainty of the benchmark reference and agree well with the averages of other participants. The MPACT k{sub eff} differs from the benchmark results for rodded and un-rodded cases by 11 and -40 pcm, respectively. The calculations were performed for various numbers of processors and parallel decompositions up to 15625 processors; all producing the same result at convergence. The parallel efficiency of the worst case was 60%, while very good efficiency (>95%) was observed for cases using 500 processors. The overall run time for the 500 processor case was 231 seconds and 19 seconds for the case with 15625 processors. Ongoing work is focused on developing theoretical performance models and the implementation of acceleration techniques to minimize the number of iterations to converge. (authors)
An automated 3D reconstruction method of UAV images
NASA Astrophysics Data System (ADS)
Liu, Jun; Wang, He; Liu, Xiaoyang; Li, Feng; Sun, Guangtong; Song, Ping
2015-10-01
In this paper a novel fully automated 3D reconstruction approach based on low-altitude unmanned aerial vehicle system (UAVs) images will be presented, which does not require previous camera calibration or any other external prior knowledge. Dense 3D point clouds are generated by integrating orderly feature extraction, image matching, structure from motion (SfM) and multi-view stereo (MVS) algorithms, overcoming many of the cost, time limitations of rigorous photogrammetry techniques. An image topology analysis strategy is introduced to speed up large scene reconstruction by taking advantage of the flight-control data acquired by UAV. Image topology map can significantly reduce the running time of feature matching by limiting the combination of images. A high-resolution digital surface model of the study area is produced base on UAV point clouds by constructing the triangular irregular network. Experimental results show that the proposed approach is robust and feasible for automatic 3D reconstruction of low-altitude UAV images, and has great potential for the acquisition of spatial information at large scales mapping, especially suitable for rapid response and precise modelling in disaster emergency.
A method of multi-view intraoral 3D measurement
NASA Astrophysics Data System (ADS)
Zhao, Huijie; Wang, Zhen; Jiang, Hongzhi; Xu, Yang; Lv, Peijun; Sun, Yunchun
2015-02-01
In dental restoration, its important to achieve a high-accuracy digital impression. Most of the existing intraoral measurement systems can only measure the tooth from a single view. Therfore - if we are wilng to acquire the whole data of a tooth, the scans of the tooth from multi-direction ad the data stitching based on the features of the surface are needed, which increases the measurement duration and influence the measurement accuracy. In this paper, we introduce a fringe-projection based on multi-view intraoral measurement system. It can acquire 3D data of the occlusal surface, the buccal surface and the lingual surface of a tooth synchronously, by using a senor with three mirrors, which aim at the three surfaces respectively and thus expand the measuring area. The constant relationship of the three mirrors is calibrated before measurement and can help stitch the data clouds acquired through different mirrors accurately. Therefore the system can obtain the 3D data of a tooth without the need to measure it from different directions for many times. Experiments proved the availability and reliability of this miniaturized measurement system.
MOM3D method of moments code theory manual
NASA Astrophysics Data System (ADS)
Shaeffer, John F.
1992-03-01
MOM3D is a FORTRAN algorithm that solves Maxwell's equations as expressed via the electric field integral equation for the electromagnetic response of open or closed three dimensional surfaces modeled with triangle patches. Two joined triangles (couples) form the vector current unknowns for the surface. Boundary conditions are for perfectly conducting or resistive surfaces. The impedance matrix represents the fundamental electromagnetic interaction of the body with itself. A variety of electromagnetic analysis options are possible once the impedance matrix is computed including backscatter radar cross section (RCS), bistatic RCS, antenna pattern prediction for user specified body voltage excitation ports, RCS image projection showing RCS scattering center locations, surface currents excited on the body as induced by specified plane wave excitation, and near field computation for the electric field on or near the body.
3D sensitivity of 6-electrode Focused Impedance Method (FIM)
NASA Astrophysics Data System (ADS)
Masum Iquebal, A. H.; Siddique-e Rabbani, K.
2010-04-01
The present work was taken up to have an understanding of the depth sensitivity of the 6 electrode FIM developed by our laboratory earlier, so that it may be applied judiciously for the measurement of organs in 3D, with electrodes on the skin surface. For a fixed electrode geometry sensitivity is expected to depend on the depth, size and conductivity of the target object. With current electrodes 18 cm apart and potential electrodes 5 cm apart, depth sensitivity of spherical conductors, insulators and of pieces of potato of different diameters were measured. The sensitivity dropped sharply with depth gradually leveling off to background, and objects could be sensed down to a depth of about twice their diameters. The sensitivity at a certain depth increases almost linearly with volume for objects with the same conductivity. Thus these results increase confidence in the use of FIM for studying organs at depths of the body.
MOM3D method of moments code theory manual
NASA Technical Reports Server (NTRS)
Shaeffer, John F.
1992-01-01
MOM3D is a FORTRAN algorithm that solves Maxwell's equations as expressed via the electric field integral equation for the electromagnetic response of open or closed three dimensional surfaces modeled with triangle patches. Two joined triangles (couples) form the vector current unknowns for the surface. Boundary conditions are for perfectly conducting or resistive surfaces. The impedance matrix represents the fundamental electromagnetic interaction of the body with itself. A variety of electromagnetic analysis options are possible once the impedance matrix is computed including backscatter radar cross section (RCS), bistatic RCS, antenna pattern prediction for user specified body voltage excitation ports, RCS image projection showing RCS scattering center locations, surface currents excited on the body as induced by specified plane wave excitation, and near field computation for the electric field on or near the body.
Real-time rendering method and performance evaluation of composable 3D lenses for interactive VR.
Borst, Christoph W; Tiesel, Jan-Phillip; Best, Christopher M
2010-01-01
We present and evaluate a new approach for real-time rendering of composable 3D lenses for polygonal scenes. Such lenses, usually called "volumetric lenses," are an extension of 2D Magic Lenses to 3D volumes in which effects are applied to scene elements. Although the composition of 2D lenses is well known, 3D composition was long considered infeasible due to both geometric and semantic complexity. Nonetheless, for a scene with multiple interactive 3D lenses, the problem of intersecting lenses must be considered. Intersecting 3D lenses in meaningful ways supports new interfaces such as hierarchical 3D windows, 3D lenses for managing and composing visualization options, or interactive shader development by direct manipulation of lenses providing component effects. Our 3D volumetric lens approach differs from other approaches and is one of the first to address efficient composition of multiple lenses. It is well-suited to head-tracked VR environments because it requires no view-dependent generation of major data structures, allowing caching and reuse of full or partial results. A Composite Shader Factory module composes shader programs for rendering composite visual styles and geometry of intersection regions. Geometry is handled by Boolean combinations of region tests in fragment shaders, which allows both convex and nonconvex CSG volumes for lens shape. Efficiency is further addressed by a Region Analyzer module and by broad-phase culling. Finally, we consider the handling of order effects for composed 3D lenses. PMID:20224135
Development of 3-D Ice Accretion Measurement Method
NASA Technical Reports Server (NTRS)
Lee, Sam; Broeren, Andy P.; Addy, Harold E., Jr.; Sills, Robert; Pifer, Ellen M.
2012-01-01
A research plan is currently being implemented by NASA to develop and validate the use of a commercial laser scanner to record and archive fully three-dimensional (3-D) ice shapes from an icing wind tunnel. The plan focused specifically upon measuring ice accreted in the NASA Icing Research Tunnel (IRT). The plan was divided into two phases. The first phase was the identification and selection of the laser scanning system and the post-processing software to purchase and develop further. The second phase was the implementation and validation of the selected system through a series of icing and aerodynamic tests. Phase I of the research plan has been completed. It consisted of evaluating several scanning hardware and software systems against an established selection criteria through demonstrations in the IRT. The results of Phase I showed that all of the scanning systems that were evaluated were equally capable of scanning ice shapes. The factors that differentiated the scanners were ease of use and the ability to operate in a wide range of IRT environmental conditions.
Teo, B G; Sarinder, K K S; Lim, L H S
2010-08-01
Three-dimensional (3D) models of the marginal hooks, dorsal and ventral anchors, bars and haptoral reservoirs of a parasite, Sundatrema langkawiense Lim & Gibson, 2009 (Monogenea) were developed using the polygonal modelling method in Autodesk 3ds Max (Version 9) based on two-dimensional (2D) illustrations. Maxscripts were written to rotate the modelled 3D structures. Appropriately orientated 3D haptoral hard-parts were then selected and positioned within the transparent 3D outline of the haptor and grouped together to form a complete 3D haptoral entity. This technique is an inexpensive tool for constructing 3D models from 2D illustrations for 3D visualisation of the spatial relationships between the different structural parts within organisms. PMID:20962723
Method for 3D fibre reconstruction on a microrobotic platform.
Hirvonen, J; Myllys, M; Kallio, P
2016-07-01
Automated handling of a natural fibrous object requires a method for acquiring the three-dimensional geometry of the object, because its dimensions cannot be known beforehand. This paper presents a method for calculating the three-dimensional reconstruction of a paper fibre on a microrobotic platform that contains two microscope cameras. The method is based on detecting curvature changes in the fibre centreline, and using them as the corresponding points between the different views of the images. We test the developed method with four fibre samples and compare the results with the references measured with an X-ray microtomography device. We rotate the samples through 16 different orientations on the platform and calculate the three-dimensional reconstruction to test the repeatability of the algorithm and its sensitivity to the orientation of the sample. We also test the noise sensitivity of the algorithm, and record the mismatch rate of the correspondences provided. We use the iterative closest point algorithm to align the measured three-dimensional reconstructions with the references. The average point-to-point distances between the reconstructed fibre centrelines and the references are 20-30 μm, and the mismatch rate is low. Given the manipulation tolerance, this shows that the method is well suited to automated fibre grasping. This has also been demonstrated with actual grasping experiments. PMID:26695385
New 3-D flow interpolation method on moving ADCP data
NASA Astrophysics Data System (ADS)
Tsubaki, R.; Kawahara, Y.; Muto, Y.; Fujita, I.
2012-05-01
A simple but accurate interpolation procedure for obtaining the three-dimensional distribution of three-component velocity data, from moving acoustic doppler current profiler (ADCP) observation data, is proposed. For understanding actual flow structure within a river with complex bathymetry, the three-dimensional mean velocity field provides a basic picture of the flow. For obtaining the three-dimensional distribution of three-component velocity data, in this work, anisotropic gridding was introduced in order to remove the random component of measured velocity data caused by the turbulence of the flow and measurement error. A continuity correction based on the pressure equation was used to reduce both random and systematic errors. The accuracy of the developed method was evaluated using three-dimensional flow simulation data from a detached-eddy simulation (DES). By using the procedure developed, the complex flow structure surrounding the spur dikes section in the Uji River was successfully visualized and explored. The proposed method shows superiorities in both accuracy and consistency for the interpolated velocity field, as compared to the kriging and inverse-distance weighted (IDW) methods.
Filtering method for 3D laser scanning point cloud
NASA Astrophysics Data System (ADS)
Liu, Da; Wang, Li; Hao, Yuncai; Zhang, Jun
2015-10-01
In recent years, with the rapid development of the hardware and software of the three-dimensional model acquisition, three-dimensional laser scanning technology is utilized in various aspects, especially in space exploration. The point cloud filter is very important before using the data. In the paper, considering both the processing quality and computing speed, an improved mean-shift point cloud filter method is proposed. Firstly, by analyze the relevance of the normal vector between the upcoming processing point and the near points, the iterative neighborhood of the mean-shift is selected dynamically, then the high frequency noise is constrained. Secondly, considering the normal vector of the processing point, the normal vector is updated. Finally, updated position is calculated for each point, then each point is moved in the normal vector according to the updated position. The experimental results show that the large features are retained, at the same time, the small sharp features are also existed for different size and shape of objects, so the target feature information is protected precisely. The computational complexity of the proposed method is not high, it can bring high precision results with fast speed, so it is very suitable for space application. It can also be utilized in civil, such as large object measurement, industrial measurement, car navigation etc. In the future, filter with the help of point strength will be further exploited.
Jacob's Interpretation Method Revisited: Accounting for 3-D Spatial Heterogeneity
NASA Astrophysics Data System (ADS)
Sanchez-Vila, X.; Riva, M.; Guadagnini, A.; Carrera, J.
2005-12-01
Traditional approaches to hydraulic test interpretation provide typically individual aquifer parameters, such as hydraulic conductivity (K) and storativity (S) values. The values obtained somehow incorporate some averaging values of aquifer heterogeneity, while the averaging functions are a direct consequence of the method of analysis employed. In recent years most work, casted in a stochastic framework, focused on the relationship between pumping rate and ensemble mean or variance of drawdown, thus having to pre-specify the parameters characterizing the underlying random spatial function. On the contrary, we contend that additional highly relevant information about heterogeneity can be obtained by looking to the spatial distribution of drawdown in individual realizations of the heterogeneous K field, without the need for invoking ergodic arguments. We present an analysis of the spatial distribution of time-dependent drawdown in a tridimensional aquifer produced by constant rate pumping in a fully penetrating well. The aquifer is considered of infinite extension in the x, y directions, and we assume no-flow boundaries in the aquifer top and bottom. The observation point is a fully penetrating piezometer. We consider an unknown spatial distribution of K(x,y,z), and using a perturbation expansion up to second order, we look at the late-time behavior of drawdown at any given observation vertical line. We conclude that: (1) at any given observation line the late-time behavior of drawdown would display a straight line in a drawdown versus log time plot, thus allowing the use of Jacob's method for test interpretation; (2) the slope of the straight line is the same for each observation line, thus providing a global average of K(x,y,z) through the aquifer; (3) the intercept point of the line in the same plot depends on location and is related to connectivity issues between the pumping and observation locations; (4) the intercept value is a weighted function of the local
Multi-crosswell profile 3D imaging and method
Washbourne, John K.; Rector, III, James W.; Bube, Kenneth P.
2002-01-01
Characterizing the value of a particular property, for example, seismic velocity, of a subsurface region of ground is described. In one aspect, the value of the particular property is represented using at least one continuous analytic function such as a Chebychev polynomial. The seismic data may include data derived from at least one crosswell dataset for the subsurface region of interest and may also include other data. In either instance, data may simultaneously be used from a first crosswell dataset in conjunction with one or more other crosswell datasets and/or with the other data. In another aspect, the value of the property is characterized in three dimensions throughout the region of interest using crosswell and/or other data. In still another aspect, crosswell datasets for highly deviated or horizontal boreholes are inherently useful. The method is performed, in part, by fitting a set of vertically spaced layer boundaries, represented by an analytic function such as a Chebychev polynomial, within and across the region encompassing the boreholes such that a series of layers is defined between the layer boundaries. Initial values of the particular property are then established between the layer boundaries and across the subterranean region using a series of continuous analytic functions. The continuous analytic functions are then adjusted to more closely match the value of the particular property across the subterranean region of ground to determine the value of the particular property for any selected point within the region.
Molecular energies from an incremental fragmentation method
NASA Astrophysics Data System (ADS)
Meitei, Oinam Romesh; Heßelmann, Andreas
2016-02-01
The systematic molecular fragmentation method by Collins and Deev [J. Chem. Phys. 125, 104104 (2006)] has been used to calculate total energies and relative conformational energies for a number of small and extended molecular systems. In contrast to the original approach by Collins, we have tested the accuracy of the fragmentation method by utilising an incremental scheme in which the energies at the lowest level of the fragmentation are calculated on an accurate quantum chemistry level while lower-cost methods are used to correct the low-level energies through a high-level fragmentation. In this work, the fragment energies at the lowest level of fragmentation were calculated using the random-phase approximation (RPA) and two recently developed extensions to the RPA while the incremental corrections at higher levels of the fragmentation were calculated using standard density functional theory (DFT) methods. The complete incremental fragmentation method has been shown to reproduce the supermolecule results with a very good accuracy, almost independent on the molecular type, size, or type of decomposition. The fragmentation method has also been used in conjunction with the DFT-SAPT (symmetry-adapted perturbation theory) method which enables a breakdown of the total nonbonding energy contributions into individual interaction energy terms. Finally, the potential problems of the method connected with the use of capping hydrogen atoms are analysed and two possible solutions are supplied.
A method of 3-D data information storage with virtual holography
NASA Astrophysics Data System (ADS)
Huang, Zhen; Liu, Guodong; Ren, Zhong; Zeng, Lüming
2008-12-01
In this paper, a new method of 3-D data cube based on virtual holographic storage is presented. Firstly, the data information is encoded in the form of 3-D data cube with a certain algorithm, in which the interval along coordinates between every data is d. Using the plane-scanning method, the 3-D cube can be described as a assembly of slices which are parallel planes along the coordinates at an interval of d. The dot on the slice represents a bit. The bright one means "1", while the dark one means "0". Secondly, a hologram of the 3-D cube is obtained by computer with virtual optics technology. All the information of a 3-D cube can be described by a 2-D hologram. At last, the hologram is inputted in the SLM, and recorded in the recording material by intersecting two coherent laser beams. When the 3-D data is exported, a reference light illuminates the hologram, and a CCD is used to get the object image which is a hologram of the 3-D data. Then the 3-D data is computed with virtual optical technology. Compared with 2-D data page storage, the 3-D data cube storage has outstanding performance in larger capacity of data storage and higher security of data.
Evaluation of a new method for stenosis quantification from 3D x-ray angiography images
NASA Astrophysics Data System (ADS)
Betting, Fabienne; Moris, Gilles; Knoplioch, Jerome; Trousset, Yves L.; Sureda, Francisco; Launay, Laurent
2001-05-01
A new method for stenosis quantification from 3D X-ray angiography images has been evaluated on both phantom and clinical data. On phantoms, for the parts larger or equal to 3 mm, the standard deviation of the measurement error has always found to be less or equal to 0.4 mm, and the maximum measurement error less than 0.17 mm. No clear relationship has been observed between the performances of the quantification method and the acquisition FoV. On clinical data, the 3D quantification method proved to be more robust to vessel bifurcations than its 3D equivalent. On a total of 15 clinical cases, the differences between 2D and 3D quantification were always less than 0.7 mm. The conclusion is that stenosis quantification from 3D X-4ay angiography images is an attractive alternative to quantification from 2D X-ray images.
Reconstructing photorealistic 3D models from image sequence using domain decomposition method
NASA Astrophysics Data System (ADS)
Xiong, Hanwei; Pan, Ming; Zhang, Xiangwei
2009-11-01
In the fields of industrial design, artistic design and heritage conservation, physical objects are usually digitalized by reverse engineering through some 3D scanning methods. Structured light and photogrammetry are two main methods to acquire 3D information, and both are expensive. Even if these expensive instruments are used, photorealistic 3D models are seldom available. In this paper, a new method to reconstruction photorealistic 3D models using a single camera is proposed. A square plate glued with coded marks is used to place the objects, and a sequence of about 20 images is taken. From the coded marks, the images are calibrated, and a snake algorithm is used to segment object from the background. A rough 3d model is obtained using shape from silhouettes algorithm. The silhouettes are decomposed into a combination of convex curves, which are used to partition the rough 3d model into some convex mesh patches. For each patch, the multi-view photo consistency constraints and smooth regulations are expressed as a finite element formulation, which can be resolved locally, and the information can be exchanged along the patches boundaries. The rough model is deformed into a fine 3d model through such a domain decomposition finite element method. The textures are assigned to each element mesh, and a photorealistic 3D model is got finally. A toy pig is used to verify the algorithm, and the result is exciting.
The Combined Finite-Discrete Element Method applied to the Study of Rock Fracturing Behavior in 3D
Rougier, Esteban; Bradley, Christopher R.; Broom, Scott T.; Knight, Earl E.; Munjiza, Ante; Sussman, Aviva J.; Swift, Robert P.
2011-01-01
Since its introduction the combined finite-discrete element method (FEM/DEM), has become an excellent tool to address a wide range of problems involving fracturing and fragmentation of solids. Within the context of rock mechanics, the FEM/DEM method has been applied to many complex industrial problems such as block caving, deep mining techniques, rock blasting, seismic waves, packing problems, rock crushing problems, etc. In the real world most of the problems involving fracture and fragmentation of solids are three dimensional problems. With the aim of addressing these problems an improved 2D/3D FEM/DEM capability has been developed at Los Alamos National Laboratory (LANL). These capabilities include state of the art 3D contact detection, contact interaction, constitutive material models, and fracture models. In this paper, Split Hopkinson Pressure Bar (SHPB) Brazilian experiments are simulated using this improved 2D/3D FEM/DEM approach which is implemented in LANL's MUNROU (Munjiza-Rougier) code. The results presented in this work show excellent agreement with both the SHPB experiments and previous 2D numerical simulations performed by other FEM/DEM research groups.
3D-2D registration of cerebral angiograms: a method and evaluation on clinical images.
Mitrovic, Uroš; Špiclin, Žiga; Likar, Boštjan; Pernuš, Franjo
2013-08-01
Endovascular image-guided interventions (EIGI) involve navigation of a catheter through the vasculature followed by application of treatment at the site of anomaly using live 2D projection images for guidance. 3D images acquired prior to EIGI are used to quantify the vascular anomaly and plan the intervention. If fused with the information of live 2D images they can also facilitate navigation and treatment. For this purpose 3D-2D image registration is required. Although several 3D-2D registration methods for EIGI achieve registration accuracy below 1 mm, their clinical application is still limited by insufficient robustness or reliability. In this paper, we propose a 3D-2D registration method based on matching a 3D vasculature model to intensity gradients of live 2D images. To objectively validate 3D-2D registration methods, we acquired a clinical image database of 10 patients undergoing cerebral EIGI and established "gold standard" registrations by aligning fiducial markers in 3D and 2D images. The proposed method had mean registration accuracy below 0.65 mm, which was comparable to tested state-of-the-art methods, and execution time below 1 s. With the highest rate of successful registrations and the highest capture range the proposed method was the most robust and thus a good candidate for application in EIGI. PMID:23649179
Fragmentation methods in laser lithotripsy
NASA Astrophysics Data System (ADS)
Jiang, Zhi X.; Whitehurst, Colin; King, Terence A.
1991-07-01
Following a series of opto-acoustic-mechanical and spectroscopic studies into the basic mechanisms of laser lithotripsy, a plasma-mediated opto-mechanical energy transfer model is presented. Laser energy, first absorbed by the calculi material at the surface and couples into the initiated plasma following ionization of vaporized material, is finally transformed into destructive mechanical energy via a shock wave induced by the impulsive expansion of the resulting plasma. This leads to the fragmentation of the calculi. The laser-plasma energy coupling gives a new definition for the fluence threshold to laser induced breakdown, which agrees with shock wave detection and analysis. A laser pulse shape with initial low intensity and sufficient fluence to vaporize a required amount of target material (lasting a few microsecond(s) ) and terminating in a short, intense pulse of about 1microsecond(s) or less, to couple most of the laser energy into the dense young plasma and so create high pressures, is required to produce optimum effect for laser lithotripsy. An opto-mechanical coupler has been designed which transfers the maximum laser energy into mechanically destructive energy, and successfully fragments various types of urinary and biliary calculi even including those calculi with poor surface absorptance, like pure white cystine. A solid state laser, Ho:YAG (2.1 micrometers wavelength and 150 microsecond(s) pulse width), has also been tested as an alternative to the flashlamp-excited dye laser. The underwater shock wave induced by this laser has been measured and has successfully fragmented calculi with poor absorptance in the visible region.
A Monte Carlo method for 3D thermal infrared radiative transfer
NASA Astrophysics Data System (ADS)
Chen, Y.; Liou, K. N.
2006-09-01
A 3D Monte Carlo model for specific application to the broadband thermal radiative transfer has been developed in which the emissivities for gases and cloud particles are parameterized by using a single cubic element as the building block in 3D space. For spectral integration in the thermal infrared, the correlated k-distribution method has been used for the sorting of gaseous absorption lines in multiple-scattering atmospheres involving 3D clouds. To check the Monte-Carlo simulation, we compare a variety of 1D broadband atmospheric fluxes and heating rates to those computed from the conventional plane-parallel (PP) model and demonstrate excellent agreement between the two. Comparisons of the Monte Carlo results for broadband thermal cooling rates in 3D clouds to those computed from the delta-diffusion approximation for 3D radiative transfer and the independent pixel-by-pixel approximation are subsequently carried out to understand the relative merits of these approaches.
A simple method for producing freestanding 3D microstructures by integrated photomask micromolding
NASA Astrophysics Data System (ADS)
Li, Hui
2015-12-01
Freestanding three-dimensional (3D) microstructures are widely used in micro-electro-mechanical system (MEMS) applications or can function as microdevices themselves. However, microfabrication methods for freestanding 3D microstructures have limitations in shape, size, cost, and mass production, etc. In this work, integrated photomask micromolding is demonstrated, which uses a portable UV light source and chrome glass micromolding to fabricate 3D microstructures without alignment. Specifically, a chrome layer on one side of the glass micromold shields the excess filling SU-8 photoresist from UV exposure and only the SU-8 photoresist in mold cavities is crosslinked. The 3D microstructures produced using this method have very high dimensional accuracy and the profile error is approximately 1.5%. This method can be used with features of virtually any size and shape and can be integrated into highly-parallel micromolding processes and has potential for MEMS applications.
3D modeling method for computer animate based on modified weak structured light method
NASA Astrophysics Data System (ADS)
Xiong, Hanwei; Pan, Ming; Zhang, Xiangwei
2010-11-01
A simple and affordable 3D scanner is designed in this paper. Three-dimensional digital models are playing an increasingly important role in many fields, such as computer animate, industrial design, artistic design and heritage conservation. For many complex shapes, optical measurement systems are indispensable to acquiring the 3D information. In the field of computer animate, such an optical measurement device is too expensive to be widely adopted, and on the other hand, the precision is not as critical a factor in that situation. In this paper, a new cheap 3D measurement system is implemented based on modified weak structured light, using only a video camera, a light source and a straight stick rotating on a fixed axis. For an ordinary weak structured light configuration, one or two reference planes are required, and the shadows on these planes must be tracked in the scanning process, which destroy the convenience of this method. In the modified system, reference planes are unnecessary, and size range of the scanned objects is expanded widely. A new calibration procedure is also realized for the proposed method, and points cloud is obtained by analyzing the shadow strips on the object. A two-stage ICP algorithm is used to merge the points cloud from different viewpoints to get a full description of the object, and after a series of operations, a NURBS surface model is generated in the end. A complex toy bear is used to verify the efficiency of the method, and errors range from 0.7783mm to 1.4326mm comparing with the ground truth measurement.
Accurate compressed look up table method for CGH in 3D holographic display.
Gao, Chuan; Liu, Juan; Li, Xin; Xue, Gaolei; Jia, Jia; Wang, Yongtian
2015-12-28
Computer generated hologram (CGH) should be obtained with high accuracy and high speed in 3D holographic display, and most researches focus on the high speed. In this paper, a simple and effective computation method for CGH is proposed based on Fresnel diffraction theory and look up table. Numerical simulations and optical experiments are performed to demonstrate its feasibility. The proposed method can obtain more accurate reconstructed images with lower memory usage compared with split look up table method and compressed look up table method without sacrificing the computational speed in holograms generation, so it is called accurate compressed look up table method (AC-LUT). It is believed that AC-LUT method is an effective method to calculate the CGH of 3D objects for real-time 3D holographic display where the huge information data is required, and it could provide fast and accurate digital transmission in various dynamic optical fields in the future. PMID:26831987
A multimaterial bioink method for 3D printing tunable, cell-compatible hydrogels.
Rutz, Alexandra L; Hyland, Kelly E; Jakus, Adam E; Burghardt, Wesley R; Shah, Ramille N
2015-03-01
A multimaterial bio-ink method using polyethylene glycol crosslinking is presented for expanding the biomaterial palette required for 3D bioprinting of more mimetic and customizable tissue and organ constructs. Lightly crosslinked, soft hydrogels are produced from precursor solutions of various materials and 3D printed. Rheological and biological characterizations are presented, and the promise of this new bio-ink synthesis strategy is discussed. PMID:25641220
A practical salient region feature based 3D multi-modality registration method for medical images
NASA Astrophysics Data System (ADS)
Hahn, Dieter A.; Wolz, Gabriele; Sun, Yiyong; Hornegger, Joachim; Sauer, Frank; Kuwert, Torsten; Xu, Chenyang
2006-03-01
We present a novel representation of 3D salient region features and its integration into a hybrid rigid-body registration framework. We adopt scale, translation and rotation invariance properties of those intrinsic 3D features to estimate a transform between underlying mono- or multi-modal 3D medical images. Our method combines advantageous aspects of both feature- and intensity-based approaches and consists of three steps: an automatic extraction of a set of 3D salient region features on each image, a robust estimation of correspondences and their sub-pixel accurate refinement with outliers elimination. We propose a region-growing based approach for the extraction of 3D salient region features, a solution to the problem of feature clustering and a reduction of the correspondence search space complexity. Results of the developed algorithm are presented for both mono- and multi-modal intra-patient 3D image pairs (CT, PET and SPECT) that have been acquired for change detection, tumor localization, and time based intra-person studies. The accuracy of the method is clinically evaluated by a medical expert with an approach that measures the distance between a set of selected corresponding points consisting of both anatomical and functional structures or lesion sites. This demonstrates the robustness of the proposed method to image overlap, missing information and artefacts. We conclude by discussing potential medical applications and possibilities for integration into a non-rigid registration framework.
Cylinder Fragmentation Using Gas Gun Methods
NASA Astrophysics Data System (ADS)
Thornhill, Tom; Reinhart, William; Chhabildas, Lalit; Grady, Dennis; Wilson, Leonard
2001-06-01
An experimental technique for investigating fracture and fragmentation characteristics of materials has been developed for use on the gas guns. In this method the candidate material is in the cylindrical form. This technique involves the precision alignment of the candidate cylinder, and symmetric impact of a stationary cylinder plug with the moving projectile from the gun. This test method allows the study of cylinder fragmentation in a laboratory environment under well-controlled loading conditions. In this presentation, results of several experiments on Aermet steel will be presented. The fragmentation toughness of the material can be estimated through knowledge of the material strain-rate and mean fragment size derived from the statistical distribution of the fragments. The values for fragmentation toughness will be compared with those obtained from other experimental methods such as explosives loading or ball on plate impact methods. Future developments and directions in test geometry, test methods and diagnostics will also be reported. This work was supported by the U. S. Department of Energy under contract DE-AC04-94AL85000.
NASA Astrophysics Data System (ADS)
Bai, Lu; Wang, Hongbo; Zhou, Jiangfan; Yang, Rong; Zhang, Hui
2014-11-01
In this paper, the aperture change of the aluminium alloy aerospace structure under real load is researched. Static experiments are carried on which is simulated the load environment of flight course. Compared with the traditional methods, through experiments results, it's proved that 3D digital speckle correlation method has good adaptability and precision on testing aperture change, and it can satisfy measurement on non-contact,real-time 3D deformation or stress concentration. The test results of new method is compared with the traditional method.
High efficient methods of content-based 3D model retrieval
NASA Astrophysics Data System (ADS)
Wu, Yuanhao; Tian, Ling; Li, Chenggang
2013-03-01
Content-based 3D model retrieval is of great help to facilitate the reuse of existing designs and to inspire designers during conceptual design. However, there is still a gap to apply it in industry due to the low time efficiency. This paper presents two new methods with high efficiency to build a Content-based 3D model retrieval system. First, an improvement is made on the "Shape Distribution (D2)" algorithm, and a new algorithm named "Quick D2" is proposed. Four sample 3D mechanical models are used in an experiment to compare the time cost of the two algorithms. The result indicates that the time cost of Quick D2 is much lower than that of D2, while the descriptors extracted by the two algorithms are almost the same. Second, an expandable 3D model repository index method with high performance, namely, RBK index, is presented. On the basis of RBK index, the search space is pruned effectively during the search process, leading to a speed up of the whole system. The factors that influence the values of the key parameters of RBK index are discussed and an experimental method to find the optimal values of the key parameters is given. Finally, "3D Searcher", a content-based 3D model retrieval system is developed. By using the methods proposed, the time cost for the system to respond one query online is reduced by 75% on average. The system has been implemented in a manufacturing enterprise, and practical query examples during a case of the automobile rear axle design are also shown. The research method presented shows a new research perspective and can effectively improve the content-based 3D model retrieval efficiency.
Modified Anderson Method for Accelerating 3D-RISM Calculations Using Graphics Processing Unit.
Maruyama, Yutaka; Hirata, Fumio
2012-09-11
A fast algorithm is proposed to solve the three-dimensional reference interaction site model (3D-RISM) theory on a graphics processing unit (GPU). 3D-RISM theory is a powerful tool for investigating biomolecular processes in solution; however, such calculations are often both memory-intensive and time-consuming. We sought to accelerate these calculations using GPUs, but to work around the problem of limited memory size in GPUs, we modified the less memory-intensive "Anderson method" to give faster convergence to 3D-RISM calculations. Using this method on a Tesla C2070 GPU, we reduced the total computational time by a factor of 8, 1.4 times by the modified Andersen method and 5.7 times by GPU, compared to calculations on an Intel Xeon machine (eight cores, 3.33 GHz) with the conventional method. PMID:26605714
Implementation of algebraic stress models in a general 3-D Navier-Stokes method (PAB3D)
NASA Technical Reports Server (NTRS)
Abdol-Hamid, Khaled S.
1995-01-01
A three-dimensional multiblock Navier-Stokes code, PAB3D, which was developed for propulsion integration and general aerodynamic analysis, has been used extensively by NASA Langley and other organizations to perform both internal (exhaust) and external flow analysis of complex aircraft configurations. This code was designed to solve the simplified Reynolds Averaged Navier-Stokes equations. A two-equation k-epsilon turbulence model has been used with considerable success, especially for attached flows. Accurate predicting of transonic shock wave location and pressure recovery in separated flow regions has been more difficult. Two algebraic Reynolds stress models (ASM) have been recently implemented in the code that greatly improved the code's ability to predict these difficult flow conditions. Good agreement with Direct Numerical Simulation (DNS) for a subsonic flat plate was achieved with ASM's developed by Shih, Zhu, and Lumley and Gatski and Speziale. Good predictions were also achieved at subsonic and transonic Mach numbers for shock location and trailing edge boattail pressure recovery on a single-engine afterbody/nozzle model.
Sippl, Wolfgang
2002-12-01
We have recently reported the development of a 3-D QSAR model for estrogen receptor ligands showing a significant correlation between calculated molecular interaction fields and experimentally measured binding affinity. The ligand alignment obtained from docking simulations was taken as basis for a comparative field analysis applying the GRID/GOLPE program. Using the interaction field derived with a water probe and applying the smart region definition (SRD) variable selection procedure, a significant and robust model was obtained (q(2)(LOO)=0.921, SDEP=0.345). To further analyze the robustness and the predictivity of the established model several recently developed estrogen receptor ligands were selected as external test set. An excellent agreement between predicted and experimental binding data was obtained indicated by an external SDEP of 0.531. Two other traditionally used prediction techniques were applied in order to check the performance of the receptor-based 3-D QSAR procedure. The interaction energies calculated on the basis of receptor-ligand complexes were correlated with experimentally observed affinities. Also ligand-based 3-D QSAR models were generated using program FlexS. The interaction energy-based model, as well as the ligand-based 3-D QSAR models yielded models with lower predictivity. The comparison with the interaction energy-based model and with the ligand-based 3-D QSAR models, respectively, indicates that the combination of receptor-based and 3-D QSAR methods is able to improve the quality of prediction. PMID:12413831
Method for making a single-step etch mask for 3D monolithic nanostructures.
Grishina, D A; Harteveld, C A M; Woldering, L A; Vos, W L
2015-12-18
Current nanostructure fabrication by etching is usually limited to planar structures as they are defined by a planar mask. The realization of three-dimensional (3D) nanostructures by etching requires technologies beyond planar masks. We present a method for fabricating a 3D mask that allows one to etch three-dimensional monolithic nanostructures using only CMOS-compatible processes. The mask is written in a hard-mask layer that is deposited on two adjacent inclined surfaces of a Si wafer. By projecting in a single step two different 2D patterns within one 3D mask on the two inclined surfaces, the mutual alignment between the patterns is ensured. Thereby after the mask pattern is defined, the etching of deep pores in two oblique directions yields a three-dimensional structure in Si. As a proof of concept we demonstrate 3D mask fabrication for three-dimensional diamond-like photonic band gap crystals in silicon. The fabricated crystals reveal a broad stop gap in optical reflectivity measurements. We propose how 3D nanostructures with five different Bravais lattices can be realized, namely cubic, tetragonal, orthorhombic, monoclinic and hexagonal, and demonstrate a mask for a 3D hexagonal crystal. We also demonstrate the mask for a diamond-structure crystal with a 3D array of cavities. In general, the 2D patterns on the different surfaces can be completely independently structured and still be in perfect mutual alignment. Indeed, we observe an alignment accuracy of better than 3.0 nm between the 2D mask patterns on the inclined surfaces, which permits one to etch well-defined monolithic 3D nanostructures. PMID:26581317
Efficient fabrication method of nano-grating for 3D holographic display with full parallax views.
Wan, Wenqiang; Qiao, Wen; Huang, Wenbin; Zhu, Ming; Fang, Zongbao; Pu, Donglin; Ye, Yan; Liu, Yanhua; Chen, Linsen
2016-03-21
Without any special glasses, multiview 3D displays based on the diffractive optics can present high resolution, full-parallax 3D images in an ultra-wide viewing angle. The enabling optical component, namely the phase plate, can produce arbitrarily distributed view zones by carefully designing the orientation and the period of each nano-grating pixel. However, such 3D display screen is restricted to a limited size due to the time-consuming fabricating process of nano-gratings on the phase plate. In this paper, we proposed and developed a lithography system that can fabricate the phase plate efficiently. Here we made two phase plates with full nano-grating pixel coverage at a speed of 20 mm^{2}/mins, a 500 fold increment in the efficiency when compared to the method of E-beam lithography. One 2.5-inch phase plate generated 9-view 3D images with horizontal-parallax, while the other 6-inch phase plate produced 64-view 3D images with full-parallax. The angular divergence in horizontal axis and vertical axis was 1.5 degrees, and 1.25 degrees, respectively, slightly larger than the simulated value of 1.2 degrees by Finite Difference Time Domain (FDTD). The intensity variation was less than 10% for each viewpoint, in consistency with the simulation results. On top of each phase plate, a high-resolution binary masking pattern containing amplitude information of all viewing zone was well aligned. We achieved a resolution of 400 pixels/inch and a viewing angle of 40 degrees for 9-view 3D images with horizontal parallax. In another prototype, the resolution of each view was 160 pixels/inch and the view angle was 50 degrees for 64-view 3D images with full parallax. As demonstrated in the experiments, the homemade lithography system provided the key fabricating technology for multiview 3D holographic display. PMID:27136814
Method for making a single-step etch mask for 3D monolithic nanostructures
NASA Astrophysics Data System (ADS)
Grishina, D. A.; Harteveld, C. A. M.; Woldering, L. A.; Vos, W. L.
2015-12-01
Current nanostructure fabrication by etching is usually limited to planar structures as they are defined by a planar mask. The realization of three-dimensional (3D) nanostructures by etching requires technologies beyond planar masks. We present a method for fabricating a 3D mask that allows one to etch three-dimensional monolithic nanostructures using only CMOS-compatible processes. The mask is written in a hard-mask layer that is deposited on two adjacent inclined surfaces of a Si wafer. By projecting in a single step two different 2D patterns within one 3D mask on the two inclined surfaces, the mutual alignment between the patterns is ensured. Thereby after the mask pattern is defined, the etching of deep pores in two oblique directions yields a three-dimensional structure in Si. As a proof of concept we demonstrate 3D mask fabrication for three-dimensional diamond-like photonic band gap crystals in silicon. The fabricated crystals reveal a broad stop gap in optical reflectivity measurements. We propose how 3D nanostructures with five different Bravais lattices can be realized, namely cubic, tetragonal, orthorhombic, monoclinic and hexagonal, and demonstrate a mask for a 3D hexagonal crystal. We also demonstrate the mask for a diamond-structure crystal with a 3D array of cavities. In general, the 2D patterns on the different surfaces can be completely independently structured and still be in perfect mutual alignment. Indeed, we observe an alignment accuracy of better than 3.0 nm between the 2D mask patterns on the inclined surfaces, which permits one to etch well-defined monolithic 3D nanostructures.
2D-3D hybrid stabilized finite element method for tsunami runup simulations
NASA Astrophysics Data System (ADS)
Takase, S.; Moriguchi, S.; Terada, K.; Kato, J.; Kyoya, T.; Kashiyama, K.; Kotani, T.
2016-09-01
This paper presents a two-dimensional (2D)-three-dimensional (3D) hybrid stabilized finite element method that enables us to predict a propagation process of tsunami generated in a hypocentral region, which ranges from offshore propagation to runup to urban areas, with high accuracy and relatively low computational costs. To be more specific, the 2D shallow water equation is employed to simulate the propagation of offshore waves, while the 3D Navier-Stokes equation is employed for the runup in urban areas. The stabilized finite element method is utilized for numerical simulations for both of the 2D and 3D domains that are independently discretized with unstructured meshes. The multi-point constraint and transmission methods are applied to satisfy the continuity of flow velocities and pressures at the interface between the resulting 2D and 3D meshes, since neither their spatial dimensions nor node arrangements are consistent. Numerical examples are presented to demonstrate the performance of the proposed hybrid method to simulate tsunami behavior, including offshore propagation and runup to urban areas, with substantially lower computation costs in comparison with full 3D computations.
2D-3D hybrid stabilized finite element method for tsunami runup simulations
NASA Astrophysics Data System (ADS)
Takase, S.; Moriguchi, S.; Terada, K.; Kato, J.; Kyoya, T.; Kashiyama, K.; Kotani, T.
2016-05-01
This paper presents a two-dimensional (2D)-three-dimensional (3D) hybrid stabilized finite element method that enables us to predict a propagation process of tsunami generated in a hypocentral region, which ranges from offshore propagation to runup to urban areas, with high accuracy and relatively low computational costs. To be more specific, the 2D shallow water equation is employed to simulate the propagation of offshore waves, while the 3D Navier-Stokes equation is employed for the runup in urban areas. The stabilized finite element method is utilized for numerical simulations for both of the 2D and 3D domains that are independently discretized with unstructured meshes. The multi-point constraint and transmission methods are applied to satisfy the continuity of flow velocities and pressures at the interface between the resulting 2D and 3D meshes, since neither their spatial dimensions nor node arrangements are consistent. Numerical examples are presented to demonstrate the performance of the proposed hybrid method to simulate tsunami behavior, including offshore propagation and runup to urban areas, with substantially lower computation costs in comparison with full 3D computations.
NASA Technical Reports Server (NTRS)
Zhao, W.; Newman, J. C., Jr.; Sutton, M. A.; Wu, X. R.; Shivakumar, K. N.
1995-01-01
Stress intensity factors for quarter-elliptical corner cracks emanating from a circular hole are determined using a 3-D weight function method combined with a 3-D finite element method. The 3-D finite element method is used to analyze uncracked configuration and provide stress distribution in the region where crack is to occur. Using this stress distribution as input, the 3-D weight function method is used to determine stress intensity factors. Three different loading conditions, i.e. remote tension, remote bending and wedge loading, are considered for a wide range in geometrical parameters. The significance in using 3-D uncracked stress distribution and the difference between single and double corner cracks are studied. Typical crack opening displacements are also provided. Comparisons are made with solutions available in the literature.
Detecting and estimating errors in 3D restoration methods using analog models.
NASA Astrophysics Data System (ADS)
José Ramón, Ma; Pueyo, Emilio L.; Briz, José Luis
2015-04-01
Some geological scenarios may be important for a number of socio-economic reasons, such as water or energy resources, but the available underground information is often limited, scarce and heterogeneous. A truly 3D reconstruction, which is still necessary during the decision-making process, may have important social and economic implications. For this reason, restoration methods were developed. By honoring some geometric or mechanical laws, they help build a reliable image of the subsurface. Pioneer methods were firstly applied in 2D (balanced and restored cross-sections) during the sixties and seventies. Later on, and due to the improvements of computational capabilities, they were extended to 3D. Currently, there are some academic and commercial restoration solutions; Unfold by the Université de Grenoble, Move by Midland Valley Exploration, Kine3D (on gOcad code) by Paradigm, Dynel3D by igeoss-Schlumberger. We have developed our own restoration method, Pmag3Drest (IGME-Universidad de Zaragoza), which is designed to tackle complex geometrical scenarios using paleomagnetic vectors as a pseudo-3D indicator of deformation. However, all these methods have limitations based on the assumptions they need to establish. For this reason, detecting and estimating uncertainty in 3D restoration methods is of key importance to trust the reconstructions. Checking the reliability and the internal consistency of every method, as well as to compare the results among restoration tools, is a critical issue never tackled so far because of the impossibility to test out the results in Nature. To overcome this problem we have developed a technique using analog models. We built complex geometric models inspired in real cases of superposed and/or conical folding at laboratory scale. The stratigraphic volumes were modeled using EVA sheets (ethylene vinyl acetate). Their rheology (tensile and tear strength, elongation, density etc) and thickness can be chosen among a large number of values
A simple method for the production of anti-C3d monoclonal antibody.
Cruz, Carlos; León, Graciela
2007-12-01
Production of monoclonal antibodies to C3d usually involves the purification of protein. Our method does not require C3 purification; it relies on attachment of C3b to mouse erythrocytes by activation of alternative pathways and further conversion in C3d. We prepared human complement-coated mouse red cells and sensitized mice of the same strain with our own schedule of immunization and applied the classical methods to obtain a mouse monoclonal antibody. We obtained a clone called BMS-11 which produces a monoclonal antibody of IgM class, to C3d with a title of 1:500000. The monoclonal antibody obtained has shown that it is suitable for use as an antiglobulin reagent. PMID:18158789
Simulations of Coulomb systems with slab geometry using an efficient 3D Ewald summation method
NASA Astrophysics Data System (ADS)
dos Santos, Alexandre P.; Girotto, Matheus; Levin, Yan
2016-04-01
We present a new approach to efficiently simulate electrolytes confined between infinite charged walls using a 3d Ewald summation method. The optimal performance is achieved by separating the electrostatic potential produced by the charged walls from the electrostatic potential of electrolyte. The electric field produced by the 3d periodic images of the walls is constant inside the simulation cell, with the field produced by the transverse images of the charged plates canceling out. The non-neutral confined electrolyte in an external potential can be simulated using 3d Ewald summation with a suitable renormalization of the electrostatic energy, to remove a divergence, and a correction that accounts for the conditional convergence of the resulting lattice sum. The new algorithm is at least an order of magnitude more rapid than the usual simulation methods for the slab geometry and can be further sped up by adopting a particle-particle particle-mesh approach.
A Multiscale Constraints Method Localization of 3D Facial Feature Points
Li, Hong-an; Zhang, Yongxin; Li, Zhanli; Li, Huilin
2015-01-01
It is an important task to locate facial feature points due to the widespread application of 3D human face models in medical fields. In this paper, we propose a 3D facial feature point localization method that combines the relative angle histograms with multiscale constraints. Firstly, the relative angle histogram of each vertex in a 3D point distribution model is calculated; then the cluster set of the facial feature points is determined using the cluster algorithm. Finally, the feature points are located precisely according to multiscale integral features. The experimental results show that the feature point localization accuracy of this algorithm is better than that of the localization method using the relative angle histograms. PMID:26539244
A Review of Failure Analysis Methods for Advanced 3D Microelectronic Packages
NASA Astrophysics Data System (ADS)
Li, Yan; Srinath, Purushotham Kaushik Muthur; Goyal, Deepak
2016-01-01
Advanced three dimensional (3D) packaging is a key enabler in driving form factor reduction, performance benefits, and package cost reduction, especially in the fast paced mobility and ultraportable consumer electronics segments. The high level of functional integration and the complex package architecture pose a significant challenge for conventional fault isolation (FI) and failure analysis (FA) methods. Innovative FI/FA tools and techniques are required to tackle the technical and throughput challenges. In this paper, the applications of FI and FA techniques such as Electro Optic Terahertz Pulse Reflectometry, 3D x-ray computed tomography, lock-in thermography, and novel physical sample preparation methods to 3D packages with package on package and stacked die with through silicon via configurations are reviewed, along with the key FI and FA challenges.
Uncovering the true nature of deformation microstructures using 3D analysis methods
NASA Astrophysics Data System (ADS)
Ferry, M.; Quadir, M. Z.; Afrin, N.; Xu, W.; Loeb, A.; Soe, B.; McMahon, C.; George, C.; Bassman, L.
2015-08-01
Three-dimensional electron backscatter diffraction (3D EBSD) has emerged as a powerful technique for generating 3D crystallographic information in reasonably large volumes of a microstructure. The technique uses a focused ion beam (FIB) as a high precision serial sectioning device for generating consecutive ion milled surfaces of a material, with each milled surface subsequently mapped by EBSD. The successive EBSD maps are combined using a suitable post-processing method to generate a crystallographic volume of the microstructure. The first part of this paper shows the usefulness of 3D EBSD for understanding the origin of various structural features associated with the plastic deformation of metals. The second part describes a new method for automatically identifying the various types of low and high angle boundaries found in deformed and annealed metals, particularly those associated with grains exhibiting subtle and gradual variations in orientation. We have adapted a 2D image segmentation technique, fast multiscale clustering, to 3D EBSD data using a novel variance function to accommodate quaternion data. This adaptation is capable of segmenting based on subtle and gradual variation as well as on sharp boundaries within the data. We demonstrate the excellent capabilities of this technique with application to 3D EBSD data sets generated from a range of cold rolled and annealed metals described in the paper.
Device and methods for "gold standard" registration of clinical 3D and 2D cerebral angiograms
NASA Astrophysics Data System (ADS)
Madan, Hennadii; Likar, Boštjan; Pernuš, Franjo; Å piclin, Žiga
2015-03-01
Translation of any novel and existing 3D-2D image registration methods into clinical image-guidance systems is limited due to lack of their objective validation on clinical image datasets. The main reason is that, besides the calibration of the 2D imaging system, a reference or "gold standard" registration is very difficult to obtain on clinical image datasets. In the context of cerebral endovascular image-guided interventions (EIGIs), we present a calibration device in the form of a headband with integrated fiducial markers and, secondly, propose an automated pipeline comprising 3D and 2D image processing, analysis and annotation steps, the result of which is a retrospective calibration of the 2D imaging system and an optimal, i.e., "gold standard" registration of 3D and 2D images. The device and methods were used to create the "gold standard" on 15 datasets of 3D and 2D cerebral angiograms, whereas each dataset was acquired on a patient undergoing EIGI for either aneurysm coiling or embolization of arteriovenous malformation. The use of the device integrated seamlessly in the clinical workflow of EIGI. While the automated pipeline eliminated all manual input or interactive image processing, analysis or annotation. In this way, the time to obtain the "gold standard" was reduced from 30 to less than one minute and the "gold standard" of 3D-2D registration on all 15 datasets of cerebral angiograms was obtained with a sub-0.1 mm accuracy.
Methods of constructing a 3D geological model from scatter data
Horsman, J.; Bethel, W.
1995-04-01
Most geoscience applications, such as assessment of an oil reservoir or hazardous waste site, require geological characterization of the site. Geological characterization involves analysis of spatial distributions of lithology, porosity, etc. Because of the complexity of the spatial relationships, the authors find that a 3-D model of geology is better suited for integration of many different types of data and provides a better representation of a site than a 2-D one. A 3-D model of geology is constructed from sample data obtained from field measurements, which are usually scattered. To create a volume model from scattered data, interpolation between points is required. The interpolation can be computed using one of several computational algorithms. Alternatively, a manual method may be employed, in which an interactive graphics device is used to input by hand the information that lies between the data points. For example, a mouse can be used to draw lines connecting data points with equal values. The combination of these two methods presents yet another approach. In this study, the authors will compare selected methods of 3-D geological modeling, They used a flow-based, modular visualization environment (AVS) to construct the geological models computationally. Within this system, they used three modules, scat{_}3d, trivar and scatter{_}to{_}ucd, as examples of computational methods. They compare these methods to the combined manual and computational approach. Because there are no tools readily available in AVS for this type of construction, they used a geological modeling system to demonstrate this method.
Flatbed-type 3D display systems using integral imaging method
NASA Astrophysics Data System (ADS)
Hirayama, Yuzo; Nagatani, Hiroyuki; Saishu, Tatsuo; Fukushima, Rieko; Taira, Kazuki
2006-10-01
We have developed prototypes of flatbed-type autostereoscopic display systems using one-dimensional integral imaging method. The integral imaging system reproduces light beams similar of those produced by a real object. Our display architecture is suitable for flatbed configurations because it has a large margin for viewing distance and angle and has continuous motion parallax. We have applied our technology to 15.4-inch displays. We realized horizontal resolution of 480 with 12 parallaxes due to adoption of mosaic pixel arrangement of the display panel. It allows viewers to see high quality autostereoscopic images. Viewing the display from angle allows the viewer to experience 3-D images that stand out several centimeters from the surface of the display. Mixed reality of virtual 3-D objects and real objects are also realized on a flatbed display. In seeking reproduction of natural 3-D images on the flatbed display, we developed proprietary software. The fast playback of the CG movie contents and real-time interaction are realized with the aid of a graphics card. Realization of the safety 3-D images to the human beings is very important. Therefore, we have measured the effects on the visual function and evaluated the biological effects. For example, the accommodation and convergence were measured at the same time. The various biological effects are also measured before and after the task of watching 3-D images. We have found that our displays show better results than those to a conventional stereoscopic display. The new technology opens up new areas of application for 3-D displays, including arcade games, e-learning, simulations of buildings and landscapes, and even 3-D menus in restaurants.
NASA Astrophysics Data System (ADS)
Crockett, Ethan Van
The need for clinically intuitive metrics for patient-specific quality assurance in radiation therapy has been well-documented (Zhen, Nelms et al. 2011). A novel transform method has shown to be effective at converting full-density 3D dose measurements made in a phantom to dose values in the patient geometry, enabling comparisons using clinically intuitive metrics such as dose-volume histograms (Oldham et al. 2011). This work investigates the transform method and compares its calculated dose-volume histograms (DVHs) to DVH values calculated by a Delta4 QA device (Scandidos), marking the first comparison of a true 3D system to a semi-3D device using clinical metrics. Measurements were made using Presage 3D dosimeters, which were readout by an in-house optical-CT scanner. Three patient cases were chosen for the study: one head-and-neck VMAT treatment and two spine IMRT treatments. The transform method showed good agreement with the planned dose values for all three cases. Furthermore, the transformed DVHs adhered to the planned dose with more accuracy than the Delta4 DVHs. The similarity between the Delta4 DVHs and the transformed DVHs, however, was greater for one of the spine cases than it was for the head-and-neck case, implying that the accuracy of the Delta4 Anatomy software may vary from one treatment site to another. Overall, the transform method, which incorporates data from full-density 3D dose measurements, provides clinically intuitive results that are more accurate and consistent than the corresponding results from a semi-3D Delta 4 system.
High fidelity digital inline holographic method for 3D flow measurements.
Toloui, Mostafa; Hong, Jiarong
2015-10-19
Among all the 3D optical flow diagnostic techniques, digital inline holographic particle tracking velocimetry (DIH-PTV) provides the highest spatial resolution with low cost, simple and compact optical setups. Despite these advantages, DIH-PTV suffers from major limitations including poor longitudinal resolution, human intervention (i.e. requirement for manually determined tuning parameters during tracer field reconstruction and extraction), limited tracer concentration, and expensive computations. These limitations prevent this technique from being widely used for high resolution 3D flow measurements. In this study, we present a novel holographic particle extraction method with the goal of overcoming all the major limitations of DIH-PTV. The proposed method consists of multiple steps involving 3D deconvolution, automatic signal-to-noise ratio enhancement and thresholding, and inverse iterative particle extraction. The entire method is implemented using GPU-based algorithm to increase the computational speed significantly. Validated with synthetic particle holograms, the proposed method can achieve particle extraction rate above 95% with fake particles less than 3% and maximum position error below 1.6 particle diameter for holograms with particle concentration above 3000 particles/mm^{3}. The applicability of the proposed method for DIH-PTV has been further validated using the experiment of laminar flow in a microchannel and the synthetic tracer flow fields generated using a DNS turbulent channel flow database. Such improvements will substantially enhance the implementation of DIH-PTV for 3D flow measurements and enable the potential commercialization of this technique. PMID:26480377
A 3-D aerodynamic method for the analysis of isolated horizontal-axis wind turbines
Ammara, I.; Masson, C.; Paraschivoiu, I.
1997-12-31
In most existing performance-analysis methods, wind turbines are considered isolated so that interference effects caused by other rotors or by the site topography are neglected. The main objective of this paper is to propose a practical 3-D method suitable for the study of these effects, in order to optimize the arrangement and the positioning of Horizontal-Axis Wind Turbines (HAWTs) in a wind farm. In the proposed methodology, the flow field around isolated HAWTs is predicted by solving the 3-D, time-averaged, steady-state, incompressible, Navier-Stokes equations in which the turbines are represented by distributions of momentum sources. The resulting governing equations are solved using a Control-Volume Finite Element Method (CVFEM). The fundamental aspects related to the development of a practical 3-D method are discussed in this paper, with an emphasis on some of the challenges that arose during its implementation. The current implementation is limited to the analysis of isolated HAWTs. Preliminary results have indicated that, the proposed 3-D method reaches the same level of accuracy, in terms of performance predictions, that the previously developed 2-D axisymmetric model and the well-known momentum-strip theory, while still using reasonable computers resources. It can be considered as a useful tool for the design of HAWTs. Its main advantages, however, are its intrinsic capacity to predict the details of the flow in the wake, and its capabilities of modelling arbitrary wind-turbine arrangements and including ground effects.
Accident or homicide--virtual crime scene reconstruction using 3D methods.
Buck, Ursula; Naether, Silvio; Räss, Beat; Jackowski, Christian; Thali, Michael J
2013-02-10
The analysis and reconstruction of forensically relevant events, such as traffic accidents, criminal assaults and homicides are based on external and internal morphological findings of the injured or deceased person. For this approach high-tech methods are gaining increasing importance in forensic investigations. The non-contact optical 3D digitising system GOM ATOS is applied as a suitable tool for whole body surface and wound documentation and analysis in order to identify injury-causing instruments and to reconstruct the course of event. In addition to the surface documentation, cross-sectional imaging methods deliver medical internal findings of the body. These 3D data are fused into a whole body model of the deceased. Additional to the findings of the bodies, the injury inflicting instruments and incident scene is documented in 3D. The 3D data of the incident scene, generated by 3D laser scanning and photogrammetry, is also included into the reconstruction. Two cases illustrate the methods. In the fist case a man was shot in his bedroom and the main question was, if the offender shot the man intentionally or accidentally, as he declared. In the second case a woman was hit by a car, driving backwards into a garage. It was unclear if the driver drove backwards once or twice, which would indicate that he willingly injured and killed the woman. With this work, we demonstrate how 3D documentation, data merging and animation enable to answer reconstructive questions regarding the dynamic development of patterned injuries, and how this leads to a real data based reconstruction of the course of event. PMID:22727689
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.
NASA Astrophysics Data System (ADS)
Toschi, I.; Capra, A.; De Luca, L.; Beraldin, J.-A.; Cournoyer, L.
2014-05-01
This paper discusses a methodology to evaluate the accuracy of recently developed image-based 3D modelling techniques. So far, the emergence of these novel methods has not been supported by the definition of an internationally recognized standard which is fundamental for user confidence and market growth. In order to provide an element of reflection and solution to the different communities involved in 3D imaging, a promising approach is presented in this paper for the assessment of both metric quality and limitations of an open-source suite of tools (Apero/MicMac), developed for the extraction of dense 3D point clouds from a set of unordered 2D images. The proposed procedural workflow is performed within a metrological context, through inter-comparisons with "reference" data acquired with two hemispherical laser scanners, one total station, and one laser tracker. The methodology is applied to two case studies, designed in order to analyse the software performances in dealing with both outdoor and environmentally controlled conditions, i.e. the main entrance of Cathédrale de la Major (Marseille, France) and a custom-made scene located at National Research Council of Canada 3D imaging Metrology Laboratory (Ottawa). Comparative data and accuracy evidence produced for both tests allow the study of some key factors affecting 3D model accuracy.
Estimating the complexity of 3D structural models using machine learning methods
NASA Astrophysics Data System (ADS)
Mejía-Herrera, Pablo; Kakurina, Maria; Royer, Jean-Jacques
2016-04-01
Quantifying the complexity of 3D geological structural models can play a major role in natural resources exploration surveys, for predicting environmental hazards or for forecasting fossil resources. This paper proposes a structural complexity index which can be used to help in defining the degree of effort necessary to build a 3D model for a given degree of confidence, and also to identify locations where addition efforts are required to meet a given acceptable risk of uncertainty. In this work, it is considered that the structural complexity index can be estimated using machine learning methods on raw geo-data. More precisely, the metrics for measuring the complexity can be approximated as the difficulty degree associated to the prediction of the geological objects distribution calculated based on partial information on the actual structural distribution of materials. The proposed methodology is tested on a set of 3D synthetic structural models for which the degree of effort during their building is assessed using various parameters (such as number of faults, number of part in a surface object, number of borders, ...), the rank of geological elements contained in each model, and, finally, their level of deformation (folding and faulting). The results show how the estimated complexity in a 3D model can be approximated by the quantity of partial data necessaries to simulated at a given precision the actual 3D model without error using machine learning algorithms.
NASA Astrophysics Data System (ADS)
Bhavanam, Sharada
The aim of this thesis is to numerically evaluate the mixed-mode Stress Intensity Factors (SIFs) of complex 3D structural geometries with arbitrary 3D cracks using the Symmetric Galerkin Boundary Element Method-Finite Element Method (SGBEM-FEM) Alternating Method. Various structural geometries with different loading scenarios and crack configurations were examined in this thesis to understand the behavior and trends of the mixed-mode SIFs as well as the fatigue life for these complex structural geometries. Although some 3D structures have empirical and numerical solutions that are readily available in the open literature, some do not; therefore this thesis presents the results of fracture and fatigue analyses of these 3D complex structures using the SGBEM-FEM Alternating Method to serve as reference for future studies. Furthermore, there are advantages of using the SGBEM-FEM Alternating Method compared to traditional FEM methods. For example, the fatigue-crack-growth and fatigue life can be better estimated for a structure because different fatigue models (i.e. Walker, Paris, and NASGRO) can be used within the same framework of the SGBEM-FEM Alternating Method. The FEM (un-cracked structure)/BEM(crack model) meshes are modeled independently, which speeds up the computation process and reduces the cost of human labor. A simple coarse mesh can be used for all fracture and fatigue analyses of complex structures. In this thesis, simple coarse meshes were used for 3D complex structures, which were below 5000 elements as compared to traditional FEM, which require meshes where the elements range on the order of ˜250,000 to ˜106 and sometimes even more than that.
2D and 3D Method of Characteristic Tools for Complex Nozzle Development
NASA Technical Reports Server (NTRS)
Rice, Tharen
2003-01-01
This report details the development of a 2D and 3D Method of Characteristic (MOC) tool for the design of complex nozzle geometries. These tools are GUI driven and can be run on most Windows-based platforms. The report provides a user's manual for these tools as well as explains the mathematical algorithms used in the MOC solutions.
Simulation of surface tension in 2D and 3D with smoothed particle hydrodynamics method
NASA Astrophysics Data System (ADS)
Zhang, Mingyu
2010-09-01
The methods for simulating surface tension with smoothed particle hydrodynamics (SPH) method in two dimensions and three dimensions are developed. In 2D surface tension model, the SPH particle on the boundary in 2D is detected dynamically according to the algorithm developed by Dilts [G.A. Dilts, Moving least-squares particle hydrodynamics II: conservation and boundaries, International Journal for Numerical Methods in Engineering 48 (2000) 1503-1524]. The boundary curve in 2D is reconstructed locally with Lagrangian interpolation polynomial. In 3D surface tension model, the SPH particle on the boundary in 3D is detected dynamically according to the algorithm developed by Haque and Dilts [A. Haque, G.A. Dilts, Three-dimensional boundary detection for particle methods, Journal of Computational Physics 226 (2007) 1710-1730]. The boundary surface in 3D is reconstructed locally with moving least squares (MLS) method. By transforming the coordinate system, it is guaranteed that the interface function is one-valued in the local coordinate system. The normal vector and curvature of the boundary surface are calculated according to the reconstructed boundary surface and then surface tension force can be calculated. Surface tension force acts only on the boundary particle. Density correction is applied to the boundary particle in order to remove the boundary inconsistency. The surface tension models in 2D and 3D have been applied to benchmark tests for surface tension. The ability of the current method applying to the simulation of surface tension in 2D and 3D is proved.
3D Spectral Element Method Simulations Of The Seismic Response of Caracas (Venezuela) Basin
NASA Astrophysics Data System (ADS)
Delavaud, E.; Vilotte, J.; Festa, G.; Cupillard, P.
2007-12-01
We present here 3D numerical simulations of the response of the Caracas (Venezuela) valley up to 5 Hz for different scenarios of plane wave excitation based on the regional seismicity. Attention is focused on the effects of the 3D basin geometry and of the adjacent regional topography. The simulations are performed using Spectral Element method (SEM) together with an unstructured hexahedral mesh discretization and perfectly matched layers (PML). These simulations show 3D amplification phenomena associated with complex wave reflexion, diffraction and focalisation patterns linked to the geometry of the basin. Time and frequency analysis reveal some interesting features both in terms of amplification and energy residence in the basin. The low frequency amplification pattern is mainly controlled by the early response of the basin to the incident plane wave while the high frequency amplification patterns result mainly from late arrivals where complex 3D wave diffraction phenomena are dominating and the memory of the initial excitation is lost. Interestingly enough, it is shown that H/V method correctly predict the low frequency amplification pattern when apply to the late part of the recorded seismograms. The complex high frequency amplification pattern is shown to be associated with surface wave generation at, and propagation from, sharp edges of the basin. Importance of 3D phenomena is assessed by comparison with simple 2D simulations. Significant differences in terms of time of residence, energy and amplification levels point out the interest of complete 3D modeling. In conclusions some of the limitations associated with the use of unstructured hexahedral meshes will be adressed. Despite the use of unstructured meshing tool, modeling the geometry of geological basins remain a complex and time consuming task. Possible extensions using more elaborate techniques like non conforming domain decomposition will be also discussed in conclusion.
A new method of 3D scene recognition from still images
NASA Astrophysics Data System (ADS)
Zheng, Li-ming; Wang, Xing-song
2014-04-01
Most methods of monocular visual three dimensional (3D) scene recognition involve supervised machine learning. However, these methods often rely on prior knowledge. Specifically, they learn the image scene as part of a training dataset. For this reason, when the sampling equipment or scene is changed, monocular visual 3D scene recognition may fail. To cope with this problem, a new method of unsupervised learning for monocular visual 3D scene recognition is here proposed. First, the image is made using superpixel segmentation based on the CIELAB color space values L, a, and b and on the coordinate values x and y of pixels, forming a superpixel image with a specific density. Second, a spectral clustering algorithm based on the superpixels' color characteristics and neighboring relationships was used to reduce the dimensions of the superpixel image. Third, the fuzzy distribution density functions representing sky, ground, and façade are multiplied with the segment pixels, where the expectations of these segments are obtained. A preliminary classification of sky, ground, and façade is generated in this way. Fourth, the most accurate classification images of sky, ground, and façade were extracted through the tier-1 wavelet sampling and Manhattan direction feature. Finally, a depth perception map is generated based on the pinhole imaging model and the linear perspective information of ground surface. Here, 400 images of Make3D Image data from the Cornell University website were used to test the algorithm. The experimental results showed that this unsupervised learning method provides a more effective monocular visual 3D scene recognition model than other methods.
NASA Astrophysics Data System (ADS)
Wen, D. H.; Jiang, Y. S.; Zhang, Y. Z.; Gao, Q.
2014-03-01
The theoretical and experimental investigations on the polarization imagery system of speckle statistical characteristics and speckle removing method are researched. A method to obtain two images encoded by polarization degree with a single measurement process is proposed. A theoretical model for polarization imagery system on Müller matrix is proposed. According to modern charge coupled device (CCD) imaging characteristics, speckles are divided into two kinds, namely small speckle and big speckle. Based on this model, a speckle reduction algorithm based on a dual-tree complex wavelet transform (DTCWT) and blockmatching 3D filter (BM3D) is proposed (DTBM3D). Original laser image data transformed by logarithmic compression is decomposed by DTCWT into approximation and detail subbands. Bilateral filtering is applied to the approximation subbands, and a suited BM3D filter is applied to the detail subbands. The despeckling results show that contrast improvement index and edge preserve index outperform those of traditional methods. The researches have important reference value in research of speckle noise level and removing speckle noise.
Gap-filling methods for 3D PlanTIS data.
Loukiala, A; Tuna, U; Beer, S; Jahnke, S; Ruotsalainen, U
2010-10-21
The range of positron emitters and their labeled compounds have led to high-resolution PET scanners becoming widely used, not only in clinical and pre-clinical studies but also in plant studies. A high-resolution PET scanner, plant tomographic imaging system (PlanTIS), was designed to study metabolic and physiological functions of plants noninvasively. The gantry of the PlanTIS scanner has detector-free regions. Even when the gantry of the PlanTIS is rotated during the scan, these regions result in missing sinogram bins in the acquired data. Missing data need to be estimated prior to the analytical image reconstructions in order to avoid artifacts in the final reconstructed images. In this study, we propose three gap-filling methods for estimation of the unique gaps existing in the 3D PlanTIS sinogram data. The 3D sinogram data were gap-filled either by linear interpolation in the transaxial planes or by the bicubic interpolation method (proposed for the ECAT high-resolution research tomograph) in the transradial planes or by the inpainting method in the transangular planes. Each gap-filling method independently compensates for slices in one of three orthogonal sinogram planes (transaxial, transradial and transangular planes). A 3D numerical Shepp-Logan phantom and the NEMA image quality phantom were used to evaluate the methods. The gap-filled sinograms were reconstructed using the analytical 3D reprojection (3DRP) method. The NEMA phantom sinograms were also reconstructed by the iterative reconstruction method, ordered subsets maximum a posteriori one step late (OSMAPOSL), to compare the results of gap filling followed by 3DRP with the results of OSMAPOSL reconstruction without gap filling. The three methods were evaluated quantitatively (by mean square error and coefficients of variation) over the selected regions of the 3D numerical Shepp-Logan phantom at eight different Poisson noise levels. Moreover, the NEMA phantom scan data were used in visual assessments
Gap-filling methods for 3D PlanTIS data
NASA Astrophysics Data System (ADS)
Loukiala, A.; Tuna, U.; Beer, S.; Jahnke, S.; Ruotsalainen, U.
2010-10-01
The range of positron emitters and their labeled compounds have led to high-resolution PET scanners becoming widely used, not only in clinical and pre-clinical studies but also in plant studies. A high-resolution PET scanner, plant tomographic imaging system (PlanTIS), was designed to study metabolic and physiological functions of plants noninvasively. The gantry of the PlanTIS scanner has detector-free regions. Even when the gantry of the PlanTIS is rotated during the scan, these regions result in missing sinogram bins in the acquired data. Missing data need to be estimated prior to the analytical image reconstructions in order to avoid artifacts in the final reconstructed images. In this study, we propose three gap-filling methods for estimation of the unique gaps existing in the 3D PlanTIS sinogram data. The 3D sinogram data were gap-filled either by linear interpolation in the transaxial planes or by the bicubic interpolation method (proposed for the ECAT high-resolution research tomograph) in the transradial planes or by the inpainting method in the transangular planes. Each gap-filling method independently compensates for slices in one of three orthogonal sinogram planes (transaxial, transradial and transangular planes). A 3D numerical Shepp-Logan phantom and the NEMA image quality phantom were used to evaluate the methods. The gap-filled sinograms were reconstructed using the analytical 3D reprojection (3DRP) method. The NEMA phantom sinograms were also reconstructed by the iterative reconstruction method, ordered subsets maximum a posteriori one step late (OSMAPOSL), to compare the results of gap filling followed by 3DRP with the results of OSMAPOSL reconstruction without gap filling. The three methods were evaluated quantitatively (by mean square error and coefficients of variation) over the selected regions of the 3D numerical Shepp-Logan phantom at eight different Poisson noise levels. Moreover, the NEMA phantom scan data were used in visual assessments
Efficient calculation method for realistic deep 3D scene hologram using orthographic projection
NASA Astrophysics Data System (ADS)
Igarashi, Shunsuke; Nakamura, Tomoya; Matsushima, Kyoji; Yamaguchi, Masahiro
2016-03-01
We propose a fast calculation method to synthesize a computer-generated hologram (CGH) of realistic deep three-dimensional (3D) scene. In our previous study, we have proposed a calculation method of CGH for reproducing such scene called ray-sampling-plane (RSP) method, in which light-ray information of a scene is converted to wavefront, and the wavefront is numerically propagated based on diffraction theory. In this paper, we introduce orthographic projection to the RSP method for accelerating calculation time. By numerical experiments, we verified the accelerated calculation with the ratio of 28-times compared to the conventional RSP method. The calculated CGH was fabricated by the printing system using laser lithography and demonstrated deep 3D image reconstruction in 52mm×52mm with realistic appearance effect such as gloss and translucent effect.
Computational methods for constructing protein structure models from 3D electron microscopy maps
Esquivel-Rodríguez, Juan; Kihara, Daisuke
2013-01-01
Protein structure determination by cryo-electron microscopy (EM) has made significant progress in the past decades. Resolutions of EM maps have been improving as evidenced by recently reported structures that are solved at high resolutions close to 3 Å. Computational methods play a key role in interpreting EM data. Among many computational procedures applied to an EM map to obtain protein structure information, in this article we focus on reviewing computational methods that model protein three-dimensional (3D) structures from a 3D EM density map that is constructed from two-dimensional (2D) maps. The computational methods we discuss range from de novo methods, which identify structural elements in an EM map, to structure fitting methods, where known high resolution structures are fit into a low-resolution EM map. A list of available computational tools is also provided. PMID:23796504
Charged-particle Gun Design with 3D Finite-element Methods
NASA Astrophysics Data System (ADS)
Humphries, Stanley
2002-04-01
The DARHT second-axis injector poses a major challenge for computer simulation. The relativistic electrons are subject to strong beam-generated electric and magnetic forces. The beam and applied fields are fully three-dimensional. Furthermore, accurate field calculations at surfaces are critical to model Child-law emission. Although several 2D relativistic beam codes are available, there is presently no 3D tool that can address all important processes in the DARHT injector. As a result, we created the OmniTrak 3D finite-element code suite. This talk gives a basic tutorial on finite-element methods with emphasis on electron gun design via the ray-tracing technique. Four main areas are covered: 1) the mesh as a tool to organize space, 2) transformation of the Poisson equation through the minimum residual principle, 3) orbit tracking in a complex environment and 4) handling self-consistent beam-generated fields. The components of a volume mesh (elements, nodes and facets) are reviewed. We consider motivations for choosing a 3D mesh style: structured versus unstructured, tetrahedrons versus hexahedrons. We discuss methods for taking volume integrals over arbitrary hexahedrons through normal coordinates and shape functions, leading to the fundamental field equations. The special problems of 3D magnetic field solutions and the advantages of the reduced potential method are outlined. Accurate field interpolations for orbit calculations require fast identification of occupied elements. A method for fast element identification that also yields the orbit penetration point on the element surface is described. The final topics are the assignment of charge and current to meshes from calculated orbits and techniques for space-charge-limited emission from multiple arbitrary 3D surfaces.
Importance of a 3D forward modeling tool for surface wave analysis methods
NASA Astrophysics Data System (ADS)
Pageot, Damien; Le Feuvre, Mathieu; Donatienne, Leparoux; Philippe, Côte; Yann, Capdeville
2016-04-01
Since a few years, seismic surface waves analysis methods (SWM) have been widely developed and tested in the context of subsurface characterization and have demonstrated their effectiveness for sounding and monitoring purposes, e.g., high-resolution tomography of the principal geological units of California or real time monitoring of the Piton de la Fournaise volcano. Historically, these methods are mostly developed under the assumption of semi-infinite 1D layered medium without topography. The forward modeling is generally based on Thomson-Haskell matrix based modeling algorithm and the inversion is driven by Monte-Carlo sampling. Given their efficiency, SWM have been transfered to several scale of which civil engineering structures in order to, e.g., determine the so-called V s30 parameter or assess other critical constructional parameters in pavement engineering. However, at this scale, many structures may often exhibit 3D surface variations which drastically limit the efficiency of SWM application. Indeed, even in the case of an homogeneous structure, 3D geometry can bias the dispersion diagram of Rayleigh waves up to obtain discontinuous phase velocity curves which drastically impact the 1D mean velocity model obtained from dispersion inversion. Taking advantages of high-performance computing center accessibility and wave propagation modeling algorithm development, it is now possible to consider the use of a 3D elastic forward modeling algorithm instead of Thomson-Haskell method in the SWM inversion process. We use a parallelized 3D elastic modeling code based on the spectral element method which allows to obtain accurate synthetic data with very low numerical dispersion and a reasonable numerical cost. In this study, we choose dike embankments as an illustrative example. We first show that their longitudinal geometry may have a significant effect on dispersion diagrams of Rayleigh waves. Then, we demonstrate the necessity of 3D elastic modeling as a forward
A new 3D tracking method exploiting the capabilities of digital holography in microscopy
NASA Astrophysics Data System (ADS)
Miccio, L.; Memmolo, P.; Merola, F.; Fusco, S.; Embrione, V.; Netti, P. A.; Ferraro, P.
2013-04-01
A method for 3D tracking has been developed exploiting Digital Holographic Microscopy (DHM) features. In the framework of self-consistent platform for manipulation and measurement of biological specimen we use DHM for quantitative and completely label free analysis of specimen with low amplitude contrast. Tracking capability extend the potentiality of DHM allowing to monitor the motion of appropriate probes and correlate it with sample properties. Complete 3D tracking has been obtained for the probes avoiding the issue of amplitude refocusing in traditional tracking processing. Our technique belongs to the video tracking methods that, conversely from Quadrant Photo-Diode method, opens the possibility to track multiples probes. All the common used video tracking algorithms are based on the numerical analysis of amplitude images in the focus plane and the shift of the maxima in the image plane are measured after the application of an appropriate threshold. Our approach for video tracking uses different theoretical basis. A set of interferograms is recorded and the complex wavefields are managed numerically to obtain three dimensional displacements of the probes. The procedure works properly on an higher number of probes and independently from their size. This method overcomes the traditional video tracking issues as the inability to measure the axial movement and the choice of suitable threshold mask. The novel configuration allows 3D tracking of micro-particles and simultaneously can furnish Quantitative Phase-contrast maps of tracked micro-objects by interference microscopy, without changing the configuration. In this paper, we show a new concept for a compact interferometric microscope that can ensure the multifunctionality, accomplishing accurate 3D tracking and quantitative phase-contrast analysis. Experimental results are presented and discussed for in vitro cells. Through a very simple and compact optical arrangement we show how two different functionalities
Flexible 3D reconstruction method based on phase-matching in multi-sensor system.
Wu, Qingyang; Zhang, Baichun; Huang, Jinhui; Wu, Zejun; Zeng, Zeng
2016-04-01
Considering the measuring range limitation of a single sensor system, multi-sensor system has become essential in obtaining complete image information of the object in the field of 3D image reconstruction. However, for the traditional multi-sensors worked independently in its system, there was some point in calibrating each sensor system separately. And the calibration between all single sensor systems was complicated and required a long time. In this paper, we present a flexible 3D reconstruction method based on phase-matching in multi-sensor system. While calibrating each sensor, it realizes the data registration of multi-sensor system in a unified coordinate system simultaneously. After all sensors are calibrated, the whole 3D image data directly exist in the unified coordinate system, and there is no need to calibrate the positions between sensors any more. Experimental results prove that the method is simple in operation, accurate in measurement, and fast in 3D image reconstruction. PMID:27137020
A correction method of color projection fringes in 3D contour measurement
NASA Astrophysics Data System (ADS)
Song, Li-mei; Li, Zong-yan; Chen, Chang-man; Xi, Jiang-tao; Guo, Qing-hua; Li, Xiao-jie
2015-07-01
In the three-dimensional (3D) contour measurement, the phase shift profilometry (PSP) method is the most widely used one. However, the measurement speed of PSP is very low because of the multiple projections. In order to improve the measurement speed, color grating stripes are used for measurement in this paper. During the measurement, only one color sinusoidal fringe is projected on the measured object. Therefore, the measurement speed is greatly improved. Since there is coupling or interference phenomenon between the adjacent color grating stripes, a color correction method is used to improve the measurement results. A method for correcting nonlinear error of measurement system is proposed in this paper, and the sinusoidal property of acquired image after correction is better than that before correction. Experimental results show that with these correction methods, the measurement errors can be reduced. Therefore, it can support a good foundation for the high-precision 3D reconstruction.
A support-operator method for viscoelastic wave modelling in 3-D heterogeneous media
NASA Astrophysics Data System (ADS)
Ely, Geoffrey P.; Day, Steven M.; Minster, Jean-Bernard
2008-01-01
We apply the method of support operators (SOM) to solve the 3-D, viscoelastic equations of motion for use in earthquake simulations. SOM is a generalized finite-difference method that can utilize meshes of arbitrary structure and incorporate irregular geometry. Our implementation uses a 3-D, logically rectangular, hexahedral mesh. Calculations are second-order in space and time. A correction term is employed for suppression of spurious zero-energy modes (hourglass oscillations). We develop a free surface boundary condition, and an absorbing boundary condition using the method of perfectly matched layers (PML). Numerical tests using a layered material model in a highly deformed mesh show good agreement with the frequency-wavenumber method, for resolutions greater than 10 nodes per wavelength. We also test a vertically incident P wave on a semi-circular canyon, for which results match boundary integral solutions at resolutions greater that 20 nodes per wavelength. We also demonstrate excellent parallel scalability of our code.
An improved 3D shape context registration method for non-rigid surface registration
NASA Astrophysics Data System (ADS)
Xiao, Di; Zahra, David; Bourgeat, Pierrick; Berghofer, Paula; Acosta Tamayo, Oscar; Wimberley, Catriona; Gregoire, Marie-Claude; Salvado, Olivier
2010-03-01
3D shape context is a method to define matching points between similar shapes as a pre-processing step to non-rigid registration. The main limitation of the approach is point mismatching, which includes long geodesic distance mismatch and neighbors crossing mismatch. In this paper, we propose a topological structure verification method to correct the long geodesic distance mismatch and a correspondence field smoothing method to correct the neighbors crossing mismatch. A robust 3D shape context model is proposed and further combined with thin-plate spline model for non-rigid surface registration. The method was tested on phantoms and rat hind limb skeletons from micro CT images. The results from experiments on mouse hind limb skeletons indicate that the approach is robust.
Robust method for extracting the pulmonary vascular trees from 3D MDCT images
NASA Astrophysics Data System (ADS)
Taeprasartsit, Pinyo; Higgins, William E.
2011-03-01
Segmentation of pulmonary blood vessels from three-dimensional (3D) multi-detector CT (MDCT) images is important for pulmonary applications. This work presents a method for extracting the vascular trees of the pulmonary arteries and veins, applicable to both contrast-enhanced and unenhanced 3D MDCT image data. The method finds 2D elliptical cross-sections and evaluates agreement of these cross-sections in consecutive slices to find likely cross-sections. It next employs morphological multiscale analysis to separate vessels from adjoining airway walls. The method then tracks the center of the likely cross-sections to connect them to the pulmonary vessels in the mediastinum and forms connected vascular trees spanning both lungs. A ground-truth study indicates that the method was able to detect on the order of 98% of the vessel branches having diameter >= 3.0 mm. The extracted vascular trees can be utilized for the guidance of safe bronchoscopic biopsy.
Numerical solution of 3-D magnetotelluric using vector finite element method
NASA Astrophysics Data System (ADS)
Prihantoro, Rudy; Sutarno, Doddy; Nurhasan
2015-09-01
Magnetotelluric (MT) is a passive electromagnetic (EM) method which measure natural variations of electric and magnetic vector fields at the Earth surface to map subsurface electrical conductivity/resistivity structure. In this study, we obtained numerical solution of three-dimensional (3-D) MT using vector finite element method by solving second order Maxwell differential equation describing diffusion of plane wave through the conductive earth. Rather than the nodes of the element, the edges of the element is used as a vector basis to overcome the occurrence of nonphysical solutions that usually faced by scalar (node based) finite element method. Electric vector fields formulation was used and the resulting system of equation was solved using direct solution method to obtain the electric vector field distribution throughout the earth resistivity model structure. The resulting MT response functions was verified with 1-D layered Earth and 3-D2 COMMEMI outcropping structure. Good agreement is achieved for both structure models.
Small pitch fringe projection method with multiple linear fiber arrays for 3D shape measurement
NASA Astrophysics Data System (ADS)
Hayashi, Takumi; Fujigaki, Motoharu; Murata, Yorinobu
2014-07-01
3-D shape measurement systems by contactless method are required in the quality inspections of metal molds and electronic parts in industrial fields. A grating projection method with phase-shifting method has advantages of high precision and high speed. Recently, the size of a BGA (ball grid array) becomes smaller. So the pitch of a grating pattern projected onto the specimen should be smaller. In conventional method, fringe pattern is projected using an imaging lens. The focal depth becomes smaller in the case of reduced projection. It is therefore difficult to project a grating pattern with small pitch onto an object with large incident angles. Authors recently proposed a light source stepping method using a linear LED device. It is easy to shrink the projected grating pitch with a lens because this projection method does not use an imaging lens. The pitch of the projected grating depends on the width of the light source. There is a limit to shrink the projected grating pitch according to the size of the LED chip. In this paper, a small pitch fringe projection method with multiple linear fiber arrays for 3D shape measurement is proposed. The width of the fiber array is 30μm. It is one digit smaller than the width of the LED chip. The experimental result of 3-D shape measurement with small pitch projection with large incident angles is shown.
NASA Astrophysics Data System (ADS)
Hosseininaveh Ahmadabadian, Ali; Robson, Stuart; Boehm, Jan; Shortis, Mark
2013-04-01
Multi-View Stereo (MVS) as a low cost technique for precise 3D reconstruction can be a rival for laser scanners if the scale of the model is resolved. A fusion of stereo imaging equipment with photogrammetric bundle adjustment and MVS methods, known as photogrammetric MVS, can generate correctly scaled 3D models without using any known object distances. Although a huge number of stereo images (e.g. 200 high resolution images from a small object) captured of the object contains redundant data that allows detailed and accurate 3D reconstruction, the capture and processing time is increased when a vast amount of high resolution images are employed. Moreover, some parts of the object are often missing due to the lack of coverage of all areas. These problems demand a logical selection of the most suitable stereo camera views from the large image dataset. This paper presents a method for clustering and choosing optimal stereo or optionally single images from a large image dataset. The approach focusses on the two key steps of image clustering and iterative image selection. The method is developed within a software application called Imaging Network Designer (IND) and tested by the 3D recording of a gearbox and three metric reference objects. A comparison is made between IND and CMVS, which is a free package for selecting vantage images. The final 3D models obtained from the IND and CMVS approaches are compared with datasets generated with an MMDx Nikon Laser scanner. Results demonstrate that IND can provide a better image selection for MVS than CMVS in terms of surface coordinate uncertainty and completeness.
Development of direct-inverse 3-D methods for applied aerodynamic design and analysis
NASA Technical Reports Server (NTRS)
Carlson, Leland A.
1988-01-01
Several inverse methods have been compared and initial results indicate that differences in results are primarily due to coordinate systems and fuselage representations and not to design procedures. Further, results from a direct-inverse method that includes 3-D wing boundary layer effects, wake curvature, and wake displacement are presented. These results show that boundary layer displacements must be included in the design process for accurate results.
Caldwell, John J; Davies, Thomas G; Donald, Alastair; McHardy, Tatiana; Rowlands, Martin G; Aherne, G Wynne; Hunter, Lisa K; Taylor, Kevin; Ruddle, Ruth; Raynaud, Florence I; Verdonk, Marcel; Workman, Paul; Garrett, Michelle D; Collins, Ian
2008-04-10
Fragment-based screening identified 7-azaindole as a protein kinase B inhibitor scaffold. Fragment elaboration using iterative crystallography of inhibitor-PKA-PKB chimera complexes efficiently guided improvements in the potency and selectivity of the compounds, resulting in the identification of nanomolar 6-(piperidin-1-yl)purine, 4-(piperidin-1-yl)-7-azaindole, and 4-(piperidin-1-yl)pyrrolo[2,3- d]pyrimidine inhibitors of PKBbeta with antiproliferative activity and showing pathway inhibition in cells. A divergence in the binding mode was seen between 4-aminomethylpiperidine and 4-aminopiperidine containing molecules. Selectivity for PKB vs PKA was observed with 4-aminopiperidine derivatives, and the most PKB-selective inhibitor (30-fold) showed significantly different bound conformations between PKA and PKA-PKB chimera. PMID:18345609
Performance and sensitivity evaluation of 3D spot detection methods in confocal microscopy.
Štěpka, Karel; Matula, Pavel; Matula, Petr; Wörz, Stefan; Rohr, Karl; Kozubek, Michal
2015-08-01
Reliable 3D detection of diffraction-limited spots in fluorescence microscopy images is an important task in subcellular observation. Generally, fluorescence microscopy images are heavily degraded by noise and non-specifically stained background, making reliable detection a challenging task. In this work, we have studied the performance and parameter sensitivity of eight recent methods for 3D spot detection. The study is based on both 3D synthetic image data and 3D real confocal microscopy images. The synthetic images were generated using a simulator modeling the complete imaging setup, including the optical path as well as the image acquisition process. We studied the detection performance and parameter sensitivity under different noise levels and under the influence of uneven background signal. To evaluate the parameter sensitivity, we propose a novel measure based on the gradient magnitude of the F1 score. We measured the success rate of the individual methods for different types of the image data and found that the type of image degradation is an important factor. Using the F1 score and the newly proposed sensitivity measure, we found that the parameter sensitivity is not necessarily proportional to the success rate of a method. This also provided an explanation why the best performing method for synthetic data was outperformed by other methods when applied to the real microscopy images. On the basis of the results obtained, we conclude with the recommendation of the HDome method for data with relatively low variations in quality, or the Sorokin method for image sets in which the quality varies more. We also provide alternative recommendations for high-quality images, and for situations in which detailed parameter tuning might be deemed expensive. PMID:26033916
NASA Astrophysics Data System (ADS)
Fargier, Yannick; Dore, Ludovic; Antoine, Raphael; Palma Lopes, Sérgio; Fauchard, Cyrille
2016-04-01
The extraction of subsurface materials is a key element for the economy of a nation. However, natural degradation of underground quarries is a major issue from an economic and public safety point of view. Consequently, the quarries stakeholders require relevant tools to define hazards associated to these structures. Safety assessment methods of underground quarries are recent and mainly based on rock physical properties. This kind of method leads to a certain homogeneity assumption of pillar internal properties that can cause an underestimation of the risk. Electrical Resistivity Imaging (ERI) is a widely used method that possesses two advantages to overcome this limitation. The first is to provide a qualitative understanding for the detection and monitoring of anomalies in the pillar body (e.g. faults). The second is to provide a quantitative description of the electrical resistivity distribution inside the pillar. This quantitative description can be interpreted with constitutive laws to help decision support (water content decreases the mechanical resistance of a chalk). However, conventional 2D and 3D Imaging techniques are usually applied to flat surface surveys or to surfaces with moderate topography. A 3D inversion of more complex media (case of the pillar) requires a full consideration of the geometry that was never taken into account before. The Photogrammetric technique presents a cost effective solution to obtain an accurate description of the external geometry of a complex media. However, this method has never been fully coupled with a geophysical method to enhance/improve the inversion process. Consequently we developed a complete procedure showing that photogrammetric and ERI tools can be efficiently combined to assess a complex 3D structure. This procedure includes in a first part a photogrammetric survey, a processing stage with an open source software and a post-processing stage finalizing a 3D surface model. The second part necessitates the
Finite volume and finite element methods applied to 3D laminar and turbulent channel flows
Louda, Petr; Příhoda, Jaromír; Sváček, Petr; Kozel, Karel
2014-12-10
The work deals with numerical simulations of incompressible flow in channels with rectangular cross section. The rectangular cross section itself leads to development of various secondary flow patterns, where accuracy of simulation is influenced by numerical viscosity of the scheme and by turbulence modeling. In this work some developments of stabilized finite element method are presented. Its results are compared with those of an implicit finite volume method also described, in laminar and turbulent flows. It is shown that numerical viscosity can cause errors of same magnitude as different turbulence models. The finite volume method is also applied to 3D turbulent flow around backward facing step and good agreement with 3D experimental results is obtained.
An Application of the Method of Arbitrary Lines to 3D Elastic Stress Analysis
NASA Astrophysics Data System (ADS)
Kaminishi, Ken; Ando, Ryuma
The MAL (Method of Arbitrary Lines) is a technique of reducing a partial differential equation to a system of ordinary differential equations. It is known that relevant use of this procedure yields high accuracy in some problems of two-dimensional elasticity and elastoplasticity. Since the basic concept of MAL is simple and based on generality, it is expected that many problems in other fields will be effectively solvable by this method. In this study, we consider the application of MAL to 3D (three-dimensional) elasticity analysis. We first give a MAL formulation of 3D elasticity problems, and demonstrate its effectiveness and accuracy for a typical problem. The reported numerical results are compared with the exact solution or that of the finite element method (FEM).
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.
Analysis of method of 3D shape reconstruction using scanning deflectometry
NASA Astrophysics Data System (ADS)
Novák, Jiří; Novák, Pavel; Mikš, Antonín.
2013-04-01
This work presents a scanning deflectometric approach to solving a 3D surface reconstruction problem, which is based on measurements of a surface gradient of optically smooth surfaces. It is shown that a description of this problem leads to a nonlinear partial differential equation (PDE) of the first order, from which the surface shape can be reconstructed numerically. The method for effective finding of the solution of this differential equation is proposed, which is based on the transform of the problem of PDE solving to the optimization problem. We describe different types of surface description for the shape reconstruction and a numerical simulation of the presented method is performed. The reconstruction process is analyzed by computer simulations and presented on examples. The performed analysis confirms a robustness of the reconstruction method and a good possibility for measurements and reconstruction of the 3D shape of specular surfaces.
Intrathoracic tumour motion estimation from CT imaging using the 3D optical flow method
NASA Astrophysics Data System (ADS)
Guerrero, Thomas; Zhang, Geoffrey; Huang, Tzung-Chi; Lin, Kang-Ping
2004-09-01
The purpose of this work was to develop and validate an automated method for intrathoracic tumour motion estimation from breath-hold computed tomography (BH CT) imaging using the three-dimensional optical flow method (3D OFM). A modified 3D OFM algorithm provided 3D displacement vectors for each voxel which were used to map tumour voxels on expiration BH CT onto inspiration BH CT images. A thoracic phantom and simulated expiration/inspiration BH CT pairs were used for validation. The 3D OFM was applied to the measured inspiration and expiration BH CT images from one lung cancer and one oesophageal cancer patient. The resulting displacements were plotted in histogram format and analysed to provide insight regarding the tumour motion. The phantom tumour displacement was measured as 1.20 and 2.40 cm with full-width at tenth maximum (FWTM) for the distribution of displacement estimates of 0.008 and 0.006 cm, respectively. The maximum error of any single voxel's motion estimate was 1.1 mm along the z-dimension or approximately one-third of the z-dimension voxel size. The simulated BH CT pairs revealed an rms error of less than 0.25 mm. The displacement of the oesophageal tumours was nonuniform and up to 1.4 cm, this was a new finding. A lung tumour maximum displacement of 2.4 cm was found in the case evaluated. In conclusion, 3D OFM provided an accurate estimation of intrathoracic tumour motion, with estimated errors less than the voxel dimension in a simulated motion phantom study. Surprisingly, oesophageal tumour motion was large and nonuniform, with greatest motion occurring at the gastro-oesophageal junction. Presented at The IASTED Second International Conference on Biomedical Engineering (BioMED 2004), Innsbruck, Austria, 16-18 February 2004.
NASA Astrophysics Data System (ADS)
Uehara, Shin-ichi; Ujike, Hiroyasu; Hamagishi, Goro; Taira, Kazuki; Koike, Takafumi; Kato, Chiaki; Nomura, Toshio; Horikoshi, Tsutomu; Mashitani, Ken; Yuuki, Akimasa; Izumi, Kuniaki; Hisatake, Yuzo; Watanabe, Naoko; Umezu, Naoaki; Nakano, Yoshihiko
2010-02-01
We are engaged in international standardization activities for 3D displays. We consider that for a sound development of 3D displays' market, the standards should be based on not only mechanism of 3D displays, but also human factors for stereopsis. However, we think that there is no common understanding on what the 3D display should be and that the situation makes developing the standards difficult. In this paper, to understand the mechanism and human factors, we focus on a double image, which occurs in some conditions on an autostereoscopic display. Although the double image is generally considered as an unwanted effect, we consider that whether the double image is unwanted or not depends on the situation and that there are some allowable double images. We tried to classify the double images into the unwanted and the allowable in terms of the display mechanism and visual ergonomics for stereopsis. The issues associated with the double image are closely related to performance characteristics for the autostereoscopic display. We also propose performance characteristics, measurement and analysis methods to represent interocular crosstalk and motion parallax.
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
Novel high speed method using gray level vector modulation for 3D shape measurement
NASA Astrophysics Data System (ADS)
Lin, Gui-Wen; Li, Dong; Tian, Jin-Dong
2014-11-01
Binocular Vision Technique is widely used in three-dimensional (3-D) measurement. Matching of pictures captured from two cameras is the most critical and difficult step in 3-D shape reconstruction. The method combines codedstructured light and spatial phase is usually adopted. However, being time consuming in matching, this method could not meet the requirements of real-time 3-D vision. In order to satisfy the high speed characteristic of real-time measurement, a novel method using gray level vector modulation is introduced. Combining binary code with gray coding principle, new coding patterns using gray level vector method is designed and projected onto the object surface. Each pixel corresponds to the designed sequence of gray values as a feature vector. The unique gray level vector is then dimensionally reduced to a resulting value which could be used as characteristic information for binocular matching. Experimental results further demonstrated the correctness and feasibility of the proposed method with fewer component patterns and less computational time.
Earthquake source tensor inversion with the gCAP method and 3D Green's functions
NASA Astrophysics Data System (ADS)
Zheng, J.; Ben-Zion, Y.; Zhu, L.; Ross, Z.
2013-12-01
We develop and apply a method to invert earthquake seismograms for source properties using a general tensor representation and 3D Green's functions. The method employs (i) a general representation of earthquake potency/moment tensors with double couple (DC), compensated linear vector dipole (CLVD), and isotropic (ISO) components, and (ii) a corresponding generalized CAP (gCap) scheme where the continuous wave trains are broken into Pnl and surface waves (Zhu & Ben-Zion, 2013). For comparison, we also use the waveform inversion method of Zheng & Chen (2012) and Ammon et al. (1998). Sets of 3D Green's functions are calculated on a grid of 1 km3 using the 3-D community velocity model CVM-4 (Kohler et al. 2003). A bootstrap technique is adopted to establish robustness of the inversion results using the gCap method (Ross & Ben-Zion, 2013). Synthetic tests with 1-D and 3-D waveform calculations show that the source tensor inversion procedure is reasonably reliable and robust. As initial application, the method is used to investigate source properties of the March 11, 2013, Mw=4.7 earthquake on the San Jacinto fault using recordings of ~45 stations up to ~0.2Hz. Both the best fitting and most probable solutions include ISO component of ~1% and CLVD component of ~0%. The obtained ISO component, while small, is found to be a non-negligible positive value that can have significant implications for the physics of the failure process. Work on using higher frequency data for this and other earthquakes is in progress.
A Quality Assurance Method that Utilizes 3D Dosimetry and Facilitates Clinical Interpretation
Oldham, Mark; Thomas, Andrew; O'Daniel, Jennifer; Juang, Titania; Ibbott, Geoffrey; Adamovics, John; Kirkpatrick, John P.
2012-10-01
Purpose: To demonstrate a new three-dimensional (3D) quality assurance (QA) method that provides comprehensive dosimetry verification and facilitates evaluation of the clinical significance of QA data acquired in a phantom. Also to apply the method to investigate the dosimetric efficacy of base-of-skull (BOS) intensity-modulated radiotherapy (IMRT) treatment. Methods and Materials: Two types of IMRT QA verification plans were created for 6 patients who received BOS IMRT. The first plan enabled conventional 2D planar IMRT QA using the Varian portal dosimetry system. The second plan enabled 3D verification using an anthropomorphic head phantom. In the latter, the 3D dose distribution was measured using the DLOS/Presage dosimetry system (DLOS = Duke Large-field-of-view Optical-CT System, Presage Heuris Pharma, Skillman, NJ), which yielded isotropic 2-mm data throughout the treated volume. In a novel step, measured 3D dose distributions were transformed back to the patient's CT to enable calculation of dose-volume histograms (DVH) and dose overlays. Measured and planned patient DVHs were compared to investigate clinical significance. Results: Close agreement between measured and calculated dose distributions was observed for all 6 cases. For gamma criteria of 3%, 2 mm, the mean passing rate for portal dosimetry was 96.8% (range, 92.0%-98.9%), compared to 94.9% (range, 90.1%-98.9%) for 3D. There was no clear correlation between 2D and 3D passing rates. Planned and measured dose distributions were evaluated on the patient's anatomy, using DVH and dose overlays. Minor deviations were detected, and the clinical significance of these are presented and discussed. Conclusions: Two advantages accrue to the methods presented here. First, treatment accuracy is evaluated throughout the whole treated volume, yielding comprehensive verification. Second, the clinical significance of any deviations can be assessed through the generation of DVH curves and dose overlays on the patient
Enhanced Rgb-D Mapping Method for Detailed 3d Modeling of Large Indoor Environments
NASA Astrophysics Data System (ADS)
Tang, Shengjun; Zhu, Qing; Chen, Wu; Darwish, Walid; Wu, Bo; Hu, Han; Chen, Min
2016-06-01
RGB-D sensors are novel sensing systems that capture RGB images along with pixel-wise depth information. Although they are widely used in various applications, RGB-D sensors have significant drawbacks with respect to 3D dense mapping of indoor environments. First, they only allow a measurement range with a limited distance (e.g., within 3 m) and a limited field of view. Second, the error of the depth measurement increases with increasing distance to the sensor. In this paper, we propose an enhanced RGB-D mapping method for detailed 3D modeling of large indoor environments by combining RGB image-based modeling and depth-based modeling. The scale ambiguity problem during the pose estimation with RGB image sequences can be resolved by integrating the information from the depth and visual information provided by the proposed system. A robust rigid-transformation recovery method is developed to register the RGB image-based and depth-based 3D models together. The proposed method is examined with two datasets collected in indoor environments for which the experimental results demonstrate the feasibility and robustness of the proposed method
Equivalent Body Force Finite Elements Method and 3-D Earth Model Applied In 2004 Sumatra Earthquake
NASA Astrophysics Data System (ADS)
Qu, W.; Cheng, H.; Shi, Y.
2015-12-01
The 26 December 2004 Sumatra-Andaman earthquake with moment magnitude (Mw) of 9.1 to 9.3 is the first great earthquake recorded by digital broadband, high-dynamic-range seismometers and global positioning system (GPS) equipment, which recorded many high-quality geophysical data sets. The spherical curvature is not negligible in far field especially for large event and the real Earth is laterally inhomogeneity and the analytical results still are difficult to explain the geodetic measurements. We use equivalent body force finite elements method Zhang et al. (2015) and mesh the whole earth, to compute global co-seismic displacements using four fault slip models of the 2004 Sumatra earthquake provided by different authors. Comparisons of calculated co-seismic displacements and GPS show that the confidences are well in near field for four models, and the confidences are according to different models. In the whole four models, the Chlieh model (Chlieh et al., 2007) is the best as this slip model not only accord well with near field data but also far field data. And then we use the best slip model, Chlieh model to explore influence of three dimensional lateral earth structure on both layered spherically symmetric (PREM) and real 3-D heterogeneous earth model (Crust 1.0 model and GyPSuM). Results show that the effects of 3-D heterogeneous earth model are not negligible and decrease concomitantly with increasing distance from the epicenter. The relative effects of 3-D crust model are 23% and 40% for horizontal and vertical displacements, respectively. The effects of the 3-D mantle model are much smaller than that of 3-D crust model but with wider impacting area.
Parlea, Lorena G; Sweeney, Blake A; Hosseini-Asanjan, Maryam; Zirbel, Craig L; Leontis, Neocles B
2016-07-01
RNA 3D motifs occupy places in structured RNA molecules that correspond to the hairpin, internal and multi-helix junction "loops" of their secondary structure representations. As many as 40% of the nucleotides of an RNA molecule can belong to these structural elements, which are distinct from the regular double helical regions formed by contiguous AU, GC, and GU Watson-Crick basepairs. With the large number of atomic- or near atomic-resolution 3D structures appearing in a steady stream in the PDB/NDB structure databases, the automated identification, extraction, comparison, clustering and visualization of these structural elements presents an opportunity to enhance RNA science. Three broad applications are: (1) identification of modular, autonomous structural units for RNA nanotechnology, nanobiology and synthetic biology applications; (2) bioinformatic analysis to improve RNA 3D structure prediction from sequence; and (3) creation of searchable databases for exploring the binding specificities, structural flexibility, and dynamics of these RNA elements. In this contribution, we review methods developed for computational extraction of hairpin and internal loop motifs from a non-redundant set of high-quality RNA 3D structures. We provide a statistical summary of the extracted hairpin and internal loop motifs in the most recent version of the RNA 3D Motif Atlas. We also explore the reliability and accuracy of the extraction process by examining its performance in clustering recurrent motifs from homologous ribosomal RNA (rRNA) structures. We conclude with a summary of remaining challenges, especially with regard to extraction of multi-helix junction motifs. PMID:27125735
Methods for Measuring the Orientation and Rotation Rate of 3D-printed Particles in Turbulence.
Cole, Brendan C; Marcus, Guy G; Parsa, Shima; Kramel, Stefan; Ni, Rui; Voth, Greg A
2016-01-01
Experimental methods are presented for measuring the rotational and translational motion of anisotropic particles in turbulent fluid flows. 3D printing technology is used to fabricate particles with slender arms connected at a common center. Shapes explored are crosses (two perpendicular rods), jacks (three perpendicular rods), triads (three rods in triangular planar symmetry), and tetrads (four arms in tetrahedral symmetry). Methods for producing on the order of 10,000 fluorescently dyed particles are described. Time-resolved measurements of their orientation and solid-body rotation rate are obtained from four synchronized videos of their motion in a turbulent flow between oscillating grids with Rλ = 91. In this relatively low-Reynolds number flow, the advected particles are small enough that they approximate ellipsoidal tracer particles. We present results of time-resolved 3D trajectories of position and orientation of the particles as well as measurements of their rotation rates. PMID:27404898
Optoranger: A 3D pattern matching method for bin picking applications
NASA Astrophysics Data System (ADS)
Sansoni, Giovanna; Bellandi, Paolo; Leoni, Fabio; Docchio, Franco
2014-03-01
This paper presents a new method, based on 3D vision, for the recognition of free-form objects in the presence of clutters and occlusions, ideal for robotic bin picking tasks. The method can be considered as a compromise between complexity and effectiveness. A 3D point cloud representing the scene is generated by a triangulation-based scanning system, where a fast camera acquires a blade projected by a laser source. Image segmentation is based on 2D images, and on the estimation of the distances between point pairs, to search for empty areas. Object recognition is performed using commercial software libraries integrated with custom-developed segmentation algorithms, and a database of model clouds created by means of the same scanning system.
A novel 3D constellation-masked method for physical security in hierarchical OFDMA system.
Zhang, Lijia; Liu, Bo; Xin, Xiangjun; Liu, Deming
2013-07-01
This paper proposes a novel 3D constellation-masked method to ensure the physical security in hierarchical optical orthogonal frequency division multiplexing access (OFDMA) system. The 3D constellation masking is executed on the two levels of hierarchical modulation and among different OFDM subcarriers, which is realized by the masking vectors. The Lorenz chaotic model is adopted for the generation of masking vectors in the proposed scheme. A 9.85 Gb/s encrypted hierarchical QAM OFDM signal is successfully demonstrated in the experiment. The performance of illegal optical network unit (ONU) with different masking vectors is also investigated. The proposed method is demonstrated to be secure and efficient against the commonly known attacks in the experiment. PMID:23842348
Gras, Christophe; Smith, Nikaïa; Sengmanivong, Lucie; Gandini, Mariana; Kubelka, Claire Fernandes; Herbeuval, Jean-Philippe
2013-01-31
The apoptotic ligand TNF-related apoptosis ligand (TRAIL) is expressed on the membrane of immune cells during HIV infection. The intracellular stockade of TRAIL in human primary CD4(+) T cells is not known. Here we investigated whether primary CD4(+) T cells expressed TRAIL in their intracellular compartment and whether TRAIL is relocalized on the plasma membrane under HIV activation. We found that TRAIL protein was stocked in intracellular compartment in non activated CD4(+) T cells and that the total level of TRAIL protein was not increased under HIV-1 stimulation. However, TRAIL was massively relocalized on plasma membrane when cells were cultured with HIV. Using three dimensional (3D) microscopy we localized TRAIL protein in human T cells and developed a new method to visualize plasma membrane without the need of a membrane marker. This method used the 3D interactive surface plot and bright light acquired images. PMID:23085529
Perfetti, Christopher M; Rearden, Bradley T
2014-01-01
This work introduces a new approach for calculating sensitivity coefficients for generalized neutronic responses to nuclear data uncertainties using continuous-energy Monte Carlo methods. The approach presented in this paper, known as the GEAR-MC method, allows for the calculation of generalized sensitivity coefficients for multiple responses in a single Monte Carlo calculation with no nuclear data perturbations or knowledge of nuclear covariance data. The theory behind the GEAR-MC method is presented here, and proof of principle is demonstrated by using the GEAR-MC method to calculate sensitivity coefficients for responses in several 3D, continuous-energy Monte Carlo applications.
Efficient solution on solving 3D Maxwell equations using stable semi-implicit splitting method
NASA Astrophysics Data System (ADS)
Cen, Wei; Gu, Ning
2016-05-01
In this paper, we propose an efficient solution on solving 3-dimensional (3D) time-domain Maxwell equations using the semi-implicit Crank-Nicholson (CN) method for time domain discretization with advantage of unconditional time stability. By applying the idea of fractional steps method (FSM) to the CN scheme, the proposed method provides a much simpler and efficient implementation than a direct implementation of the CN scheme. Compared with the alternating-direction implicit (ADI) method and explicit finite-difference time-domain approach (FDTD), it significantly saves the computational resource like memory and CPU time while remains similar numerical accuracy.
Reconstruction of 3D structure using stochastic methods: morphology and transport properties
NASA Astrophysics Data System (ADS)
Karsanina, Marina; Gerke, Kirill; Čapek, Pavel; Vasilyev, Roman; Korost, Dmitry; Skvortsova, Elena
2013-04-01
One of the main factors defining numerous flow phenomena in rocks, soils and other porous media, including fluid and solute movements, is pore structure, e.g., pore sizes and their connectivity. Numerous numerical methods were developed to quantify single and multi-phase flow in such media on microscale. Among most popular ones are: 1) a wide range of finite difference/element/volume solutions of Navier-Stokes equations and its simplifications; 2) lattice-Boltzmann method; 3) pore-network models, among others. Each method has some advantages and shortcomings, so that different research teams usually utilize more than one, depending on the study case. Recent progress in 3D imaging of internal structure, e.g., X-ray tomography, FIB-SEM and confocal microscopy, made it possible to obtain digitized input pore parameters for such models, however, a trade-off between resolution and sample size is usually unavoidable. There are situations then only standard two-dimensional information of porous structure is known due to tomography high cost or resolution limitations. However, physical modeling on microscale requires 3D information. There are three main approaches to reconstruct (using 2D cut(s) or some other limited information/properties) porous media: 1) statistical methods (correlation functions and simulated annealing, multi-point statistics, entropy methods), 2) sequential methods (sphere or other granular packs) and 3) morphological methods. Stochastic reconstructions using correlation functions possess some important advantage - they provide a statistical description of the structure, which is known to have relationships with all physical properties. In addition, this method is more flexible for other applications to characterize porous media. Taking different 3D scans of natural and artificial porous materials (sandstones, soils, shales, ceramics) we choose some 2D cut/s as sources of input correlation functions. Based on different types of correlation functions
Analysis of the 3D acoustic cloaking problems using optimization method
NASA Astrophysics Data System (ADS)
Alekseev, G. V.; Spivak, Yu E.
2016-06-01
Control problems for the 3D model of acoustic scattering which describes scattering acoustic waves by a permeable obstacle with the form of a spherical layer are considered. These problems arise while developing the design technologies of acoustic cloaking devices using the wave flow method. The solvability of direct and control problems for the acoustic scattering model under study is proved. The sufficient conditions which provide local uniqueness and stability of optimal solutions are established.
A Method for 3D Histopathology Reconstruction Supporting Mouse Microvasculature Analysis.
Xu, Yiwen; Pickering, J Geoffrey; Nong, Zengxuan; Gibson, Eli; Arpino, John-Michael; Yin, Hao; Ward, Aaron D
2015-01-01
Structural abnormalities of the microvasculature can impair perfusion and function. Conventional histology provides good spatial resolution with which to evaluate the microvascular structure but affords no 3-dimensional information; this limitation could lead to misinterpretations of the complex microvessel network in health and disease. The objective of this study was to develop and evaluate an accurate, fully automated 3D histology reconstruction method to visualize the arterioles and venules within the mouse hind-limb. Sections of the tibialis anterior muscle from C57BL/J6 mice (both normal and subjected to femoral artery excision) were reconstructed using pairwise rigid and affine registrations of 5 µm-thick, paraffin-embedded serial sections digitized at 0.25 µm/pixel. Low-resolution intensity-based rigid registration was used to initialize the nucleus landmark-based registration, and conventional high-resolution intensity-based registration method. The affine nucleus landmark-based registration was developed in this work and was compared to the conventional affine high-resolution intensity-based registration method. Target registration errors were measured between adjacent tissue sections (pairwise error), as well as with respect to a 3D reference reconstruction (accumulated error, to capture propagation of error through the stack of sections). Accumulated error measures were lower (p < 0.01) for the nucleus landmark technique and superior vasculature continuity was observed. These findings indicate that registration based on automatic extraction and correspondence of small, homologous landmarks may support accurate 3D histology reconstruction. This technique avoids the otherwise problematic "banana-into-cylinder" effect observed using conventional methods that optimize the pairwise alignment of salient structures, forcing them to be section-orthogonal. This approach will provide a valuable tool for high-accuracy 3D histology tissue reconstructions for
Shao, Yan-Lin Faltinsen, Odd M.
2014-10-01
We propose a new efficient and accurate numerical method based on harmonic polynomials to solve boundary value problems governed by 3D Laplace equation. The computational domain is discretized by overlapping cells. Within each cell, the velocity potential is represented by the linear superposition of a complete set of harmonic polynomials, which are the elementary solutions of Laplace equation. By its definition, the method is named as Harmonic Polynomial Cell (HPC) method. The characteristics of the accuracy and efficiency of the HPC method are demonstrated by studying analytical cases. Comparisons will be made with some other existing boundary element based methods, e.g. Quadratic Boundary Element Method (QBEM) and the Fast Multipole Accelerated QBEM (FMA-QBEM) and a fourth order Finite Difference Method (FDM). To demonstrate the applications of the method, it is applied to some studies relevant for marine hydrodynamics. Sloshing in 3D rectangular tanks, a fully-nonlinear numerical wave tank, fully-nonlinear wave focusing on a semi-circular shoal, and the nonlinear wave diffraction of a bottom-mounted cylinder in regular waves are studied. The comparisons with the experimental results and other numerical results are all in satisfactory agreement, indicating that the present HPC method is a promising method in solving potential-flow problems. The underlying procedure of the HPC method could also be useful in other fields than marine hydrodynamics involved with solving Laplace equation.
A new combined prior based reconstruction method for compressed sensing in 3D ultrasound imaging
NASA Astrophysics Data System (ADS)
Uddin, Muhammad S.; Islam, Rafiqul; Tahtali, Murat; Lambert, Andrew J.; Pickering, Mark R.
2015-03-01
Ultrasound (US) imaging is one of the most popular medical imaging modalities, with 3D US imaging gaining popularity recently due to its considerable advantages over 2D US imaging. However, as it is limited by long acquisition times and the huge amount of data processing it requires, methods for reducing these factors have attracted considerable research interest. Compressed sensing (CS) is one of the best candidates for accelerating the acquisition rate and reducing the data processing time without degrading image quality. However, CS is prone to introduce noise-like artefacts due to random under-sampling. To address this issue, we propose a combined prior-based reconstruction method for 3D US imaging. A Laplacian mixture model (LMM) constraint in the wavelet domain is combined with a total variation (TV) constraint to create a new regularization regularization prior. An experimental evaluation conducted to validate our method using synthetic 3D US images shows that it performs better than other approaches in terms of both qualitative and quantitative measures.
NASA Astrophysics Data System (ADS)
Grigoraş, I.-R.; Covăsnianu, A.; Pleşu, G.; Benedict, B.
2009-04-01
The paper describes an experiment which took place in Iasi town, Romania, consisted in two different topographical survey techniques applied for one and the same objective placed in a block within the city (western part) - a thermal power station. The purpose was to compare those methods and to determine which one is proper to be used in this domain in terms of fastness, optimization and speed of data processing. First technique applied for our survey was the classical one, with a total station. Using the CAD technique, we obtained a final product (a dwg file) and a list of coordinates (a text file). The second method, which we focused our attention more, was the measurement with a very precise 3D laser scanstation, also very suitable in archeology. The data obtained were processed with special software. Result was a 3D model of the thermal power plant composed of measurable cloud point data. Finally, analyzing the advantages and disadvantages of each method, we came to the conclusion that the 3D laser scanning which we used matches well the application, in this case civil engineering, but the future of accepting and implementing this technique is in the hands of Romanian authorities.
NASA Astrophysics Data System (ADS)
Celli, Jonathan P.; Rizvi, Imran; Blanden, Adam R.; Abu-Yousif, Adnan O.; Spring, Bryan Q.; Hasan, Tayyaba
2011-02-01
Three-dimensional in vitro tumor models have emerged as powerful research tools in cancer biology, though the vast potential of these systems as high-throughput, biologically relevant reporters of treatment response has yet to be adequately explored. Here, building on previous studies, we demonstrate the utility of using 3D models for ovarian and pancreatic cancers in conjunction with quantitative image processing to reveal aspects of growth behavior and treatment response that would not be evident without either modeling or quantitative analysis component. In this report we specifically focus on recent improvements in the imaging component of this integrative research platform and emphasize analysis to establish reproducible growth properties in 3D tumor arrays, a key consideration in establishing the utility of this platform as a reliable reporter of therapeutic response. Building on previous studies using automated segmentation of low magnification image fields containing large numbers of nodules to study size dependent treatment effects, we introduce an improvement to this method using multiresolution decomposition to remove gradient background from transmitted light images for more reliable feature identification. This approach facilitates the development of a new treatment response metric, disruption fraction (Dfrac), which quantifies dose dependent distribution shifts from nodular fragmentation induced by cytotoxic therapies. Using this approach we show that PDT treatment is associated with significant dose-dependent increases in Dfrac, while this is not observed with carboplatin treatment. The ability to quantify this response to therapy could play a key role in design of combination regimens involving these two modalities.
Woie, Leik; Måløy, Frode; Eftestøl, Trygve; Engan, Kjersti; Edvardsen, Thor; Kvaløy, Jan Terje; Ørn, Stein
2014-02-01
Current methods for the estimation of infarct size by late-enhanced cardiac magnetic imaging are based upon 2D analysis that first determines the size of the infarction in each slice, and thereafter adds the infarct sizes from each slice to generate a volume. We present a novel, automatic 3D method that estimates infarct size by a simultaneous analysis of all pixels from all slices. In a population of 54 patients with ischemic scars, the infarct size estimated by the automatic 3D method was compared with four established 2D methods. The new 3D method defined scar as the sum of all pixels with signal intensity (SI) ≥35 % of max SI from the complete myocardium, border zone: SI 35-50 % of max SI and core as SI ≥50 % of max SI. The 3D method yielded smaller infarct size (-2.8 ± 2.3 %) and core size (-3.0 ± 1.7 %) than the 2D method most similar to ours. There was no difference in the size of the border zone (0.2 ± 1.4 %). The 3D method demonstrated stronger correlations between scar size and left ventricular (LV) remodelling parameters (LV ejection fraction: r = -0.71, p < 0.0005, LV end-diastolic index: r = 0.54, p < 0.0005, and LV end-systolic index: r = 0.59, p < 0.0005) compared with conventional 2D methods. Infarct size estimation by our novel 3D automatic method is without the need for manual demarcation of the scar; it is less time-consuming and has a stronger correlation with remodelling parameters compared with existing methods. PMID:24249515
A Novel 2D-to-3D Video Conversion Method Using Time-Coherent Depth Maps
Yin, Shouyi; Dong, Hao; Jiang, Guangli; Liu, Leibo; Wei, Shaojun
2015-01-01
In this paper, we propose a novel 2D-to-3D video conversion method for 3D entertainment applications. 3D entertainment is getting more and more popular and can be found in many contexts, such as TV and home gaming equipment. 3D image sensors are a new method to produce stereoscopic video content conveniently and at a low cost, and can thus meet the urgent demand for 3D videos in the 3D entertaiment market. Generally, 2D image sensor and 2D-to-3D conversion chip can compose a 3D image sensor. Our study presents a novel 2D-to-3D video conversion algorithm which can be adopted in a 3D image sensor. In our algorithm, a depth map is generated by combining global depth gradient and local depth refinement for each frame of 2D video input. Global depth gradient is computed according to image type while local depth refinement is related to color information. As input 2D video content consists of a number of video shots, the proposed algorithm reuses the global depth gradient of frames within the same video shot to generate time-coherent depth maps. The experimental results prove that this novel method can adapt to different image types, reduce computational complexity and improve the temporal smoothness of generated 3D video. PMID:26131674
A Novel 2D-to-3D Video Conversion Method Using Time-Coherent Depth Maps.
Yin, Shouyi; Dong, Hao; Jiang, Guangli; Liu, Leibo; Wei, Shaojun
2015-01-01
In this paper, we propose a novel 2D-to-3D video conversion method for 3D entertainment applications. 3D entertainment is getting more and more popular and can be found in many contexts, such as TV and home gaming equipment. 3D image sensors are a new method to produce stereoscopic video content conveniently and at a low cost, and can thus meet the urgent demand for 3D videos in the 3D entertaiment market. Generally, 2D image sensor and 2D-to-3D conversion chip can compose a 3D image sensor. Our study presents a novel 2D-to-3D video conversion algorithm which can be adopted in a 3D image sensor. In our algorithm, a depth map is generated by combining global depth gradient and local depth refinement for each frame of 2D video input. Global depth gradient is computed according to image type while local depth refinement is related to color information. As input 2D video content consists of a number of video shots, the proposed algorithm reuses the global depth gradient of frames within the same video shot to generate time-coherent depth maps. The experimental results prove that this novel method can adapt to different image types, reduce computational complexity and improve the temporal smoothness of generated 3D video. PMID:26131674
A 3D front tracking method on a CPU/GPU system
Bo, Wurigen; Grove, John
2011-01-21
We describe the method to port a sequential 3D interface tracking code to a GPU with CUDA. The interface is represented as a triangular mesh. Interface geometry properties and point propagation are performed on a GPU. Interface mesh adaptation is performed on a CPU. The convergence of the method is assessed from the test problems with given velocity fields. Performance results show overall speedups from 11 to 14 for the test problems under mesh refinement. We also briefly describe our ongoing work to couple the interface tracking method with a hydro solver.
NASA Technical Reports Server (NTRS)
Nakazawa, S.
1988-01-01
This annual status report presents the results of work performed during the fourth year of the 3-D Inelastic Analysis Methods for Hot Section Components program (NASA Contract NAS3-23697). The objective of the program is to produce a series of new computer codes permitting more accurate and efficient 3-D analysis of selected hot section components, i.e., combustor liners, turbine blades and turbine vanes. The computer codes embody a progression of math models and are streamlined to take advantage of geometrical features, loading conditions, and forms of material response that distinguish each group of selected components. Volume 1 of this report discusses the special finite element models developed during the fourth year of the contract.
OPTIMIZATION OF 3-D IMAGE-GUIDED NEAR INFRARED SPECTROSCOPY USING BOUNDARY ELEMENT METHOD
Srinivasan, Subhadra; Carpenter, Colin; Pogue, Brian W.; Paulsen, Keith D.
2010-01-01
Multimodality imaging systems combining optical techniques with MRI/CT provide high-resolution functional characterization of tissue by imaging molecular and vascular biomarkers. To optimize these hybrid systems for clinical use, faster and automatable algorithms are required for 3-D imaging. Towards this end, a boundary element model was used to incorporate tissue boundaries from MRI/CT into image formation process. This method uses surface rendering to describe light propagation in 3-D using diffusion equation. Parallel computing provided speedup of up to 54% in time of computation. Simulations showed that location of NIRS probe was crucial for quantitatively accurate estimation of tumor response. A change of up to 61% was seen between cycles 1 and 3 in monitoring tissue response to neoadjuvant chemotherapy. PMID:20523751
3D In Vitro Model for Breast Cancer Research Using Magnetic Levitation and Bioprinting Method.
Leonard, Fransisca; Godin, Biana
2016-01-01
Tumor microenvironment composition and architecture are known as a major factor in orchestrating the tumor growth and its response to various therapies. In this context, in vivo studies are necessary to evaluate the responses. However, while tumor cells can be of human origin, tumor microenvironment in the in vivo models is host-based. On the other hand, in vitro studies in a flat monoculture of tumor cells (the most frequently used in vitro tumor model) are unable to recapitulate the complexity of tumor microenvironment. Three-dimensional (3D) in vitro cell cultures of tumor cells have been proven to be an important experimental tool in understanding mechanisms of tumor growth, response to therapeutics, and transport of nutrients/drugs. We have recently described a novel tool to create 3D co-cultures of tumor cells and cells in the tumor microenvironment. Our method utilizes magnetic manipulation/levitation of the specific ratios of tumor cells and cells in the tumor microenvironment (from human or animal origin) aiding in the formation of tumor spheres with defined cellular composition and density, as quickly as within 24 h. This chapter describes the experimental protocols developed to model the 3D structure of the cancer environment using the above method. PMID:26820961
3D measurement method based on combined temporal encoding structured light
NASA Astrophysics Data System (ADS)
Yu, Xiaoyang; Wang, Yang; Yu, Shuang; Cheng, Hao; Sun, Xiaoming; Yu, Shuchun; Chen, Deyun
2013-10-01
Three-dimensional (3D) vision measurement technology based on encoding structured light plays an important role and has become the main development trend in the field of 3D non-contact measurement. However, how to synthetically improve measurement speed, accuracy and sampling density is still a difficult problem. Thus in the present work, a novel 3D measurement method based on temporal encoding structured light by combining trapezoidal phase-shifting pattern and cyclic code pattern is proposed. Due to trapezoidal phase-shifting has the advantages of high sampling density and high-speed, the proposed method can maintain these advantages by using cyclic code to expand the range of trapezoidal phase-shifting. In addition, the correction scheme is designed to solve the problem of cycle dislocation. Finally, simulation experimental platform is built with 3ds max and MATLAB. Experimental analyses and results show that, the maximal error is less than 3 mm in the range from 400 mm to 1100 mm, cycle dislocation correction has a good effect.
A 3D neurovascular bundles segmentation method based on MR-TRUS deformable registration
NASA Astrophysics Data System (ADS)
Yang, Xiaofeng; Rossi, Peter; Jani, Ashesh B.; Mao, Hui; Ogunleye, Tomi; Curran, Walter J.; Liu, Tian
2015-03-01
In this paper, we propose a 3D neurovascular bundles (NVB) segmentation method for ultrasound (US) image by integrating MR and transrectal ultrasound (TRUS) images through MR-TRUS deformable registration. First, 3D NVB was contoured by a physician in MR images, and the 3D MRdefined NVB was then transformed into US images using a MR-TRUS registration method, which models the prostate tissue as an elastic material, and jointly estimates the boundary deformation and the volumetric deformations under the elastic constraint. This technique was validated with a clinical study of 6 patients undergoing radiation therapy (RT) treatment for prostate cancer. The accuracy of our approach was assessed through the locations of landmarks, as well as previous ultrasound Doppler images of patients. MR-TRUS registration was successfully performed for all patients. The mean displacement of the landmarks between the post-registration MR and TRUS images was less than 2 mm, and the average NVB volume Dice Overlap Coefficient was over 89%. This NVB segmentation technique could be a useful tool as we try to spare the NVB in prostate RT, monitor NVB response to RT, and potentially improve post-RT potency outcomes.
Method for 3D noncontact measurements of cut trees package area
NASA Astrophysics Data System (ADS)
Knyaz, Vladimir A.; Vizilter, Yuri V.
2001-02-01
Progress in imaging sensors and computers create the background for numerous 3D imaging application for wide variety of manufacturing activity. Many demands for automated precise measurements are in wood branch of industry. One of them is the accurate volume definition for cut trees carried on the truck. The key point for volume estimation is determination of the front area of the cut tree package. To eliminate slow and inaccurate manual measurements being now in practice the experimental system for automated non-contact wood measurements is developed. The system includes two non-metric CCD video cameras, PC as central processing unit, frame grabbers and original software for image processing and 3D measurements. The proposed method of measurement is based on capturing the stereo pair of front of trees package and performing the image orthotranformation into the front plane. This technique allows to process transformed image for circle shapes recognition and calculating their area. The metric characteristics of the system are provided by special camera calibration procedure. The paper presents the developed method of 3D measurements, describes the hardware used for image acquisition and the software realized the developed algorithms, gives the productivity and precision characteristics of the system.
2D and 3D visualization methods of endoscopic panoramic bladder images
NASA Astrophysics Data System (ADS)
Behrens, Alexander; Heisterklaus, Iris; Müller, Yannick; Stehle, Thomas; Gross, Sebastian; Aach, Til
2011-03-01
While several mosaicking algorithms have been developed to compose endoscopic images of the internal urinary bladder wall into panoramic images, the quantitative evaluation of these output images in terms of geometrical distortions have often not been discussed. However, the visualization of the distortion level is highly desired for an objective image-based medical diagnosis. Thus, we present in this paper a method to create quality maps from the characteristics of transformation parameters, which were applied to the endoscopic images during the registration process of the mosaicking algorithm. For a global first view impression, the quality maps are laid over the panoramic image and highlight image regions in pseudo-colors according to their local distortions. This illustration supports then surgeons to identify geometrically distorted structures easily in the panoramic image, which allow more objective medical interpretations of tumor tissue in shape and size. Aside from introducing quality maps in 2-D, we also discuss a visualization method to map panoramic images onto a 3-D spherical bladder model. Reference points are manually selected by the surgeon in the panoramic image and the 3-D model. Then the panoramic image is mapped by the Hammer-Aitoff equal-area projection onto the 3-D surface using texture mapping. Finally the textured bladder model can be freely moved in a virtual environment for inspection. Using a two-hemisphere bladder representation, references between panoramic image regions and their corresponding space coordinates within the bladder model are reconstructed. This additional spatial 3-D information thus assists the surgeon in navigation, documentation, as well as surgical planning.
NASA Astrophysics Data System (ADS)
Takahashi, M.; Kawabata, Y.; Washitani, T.; Tanaka, S.; Maeda, S.; Mimotogi, S.
2014-03-01
In progress of lithography technologies, the importance of Mask3D analysis has been emphasized because the influence of mask topography effects is not avoidable to be increased explosively. An electromagnetic filed simulation method, such as FDTD, RCWA and FEM, is applied to analyze those complicated phenomena. We have investigated Constrained Interpolation Profile (CIP) method, which is one of the Method of Characteristics (MoC), for Mask3D analysis in optical lithography. CIP method can reproduce the phase of propagating waves with less numerical error by using high order polynomial function. The restrictions of grid distance are relaxed with spatial grid. Therefore this method reduces the number of grid points in complex structure. In this paper, we study the feasibility of CIP scheme applying a non-uniform and spatial-interpolated grid to practical mask patterns. The number of grid points might be increased in complex layout and topological structure since these structures require a dense grid to remain the fidelity of each design. We propose a spatial interpolation method based on CIP method same as time-domain interpolation to reduce the number of grid points to be computed. The simulation results of two meshing methods with spatial interpolation are shown.
A Retrospective Research for 1999 Chi-Chi, Taiwan Earthquake by 3-D PI Method
NASA Astrophysics Data System (ADS)
Li, H.; Chen, C.; Tiampo, K.; Rundle, J.; Klein, W.
2007-12-01
The PI (Pattern Informatics) method was proposed by Tiampo et al., 2002, for the identification of future seismicity in California after 1999. A plausible result was published in Proc. Natl. Acad. Sci. U.S.A. 99, suppl., 2002. Chen et al. modified the calculation of the original PI method to de-emphasize the effect of current events, and applied the modified PI method to make a retrospective analysis for the 1999 Chi-Chi earthquake in Taiwan. In the case study of the Chi-Chi sequence, the main shock and most of the aftershocks were located in and around the ß¡±hot spot ßÆØ region produced by the modified PI method. Tiampo et al. (2007) applied the Thirumalai- Mountain metric to three earthquake catalogs, California, Spain and eastern Canada, which belong to different tectonic environments. Under specific spatial-temporal resolutions, effectively ergodic behaviors of seismic rate all exist in these regions. Ongoing research for Taiwan suggests that, once depth factor is considered in seismic event distribution, a similar effectively ergodicity also exists in the seismicity data. It motivates us to improve the original PI method to a 3-D version on order to consider the depth effect in a very condensed, high seismicity region. In this study, we used 3-D PI method to make a retrospective forecast of the 1999 M=7.3 Chi-Chi earthquake. The CWB (Central Weather Bureau) catalogue was used. An effectively ergodic period chosen from inverse TM metric-time plot was used as the forecast period. The main shock and several large aftershocks, which magnitudes are ß « 6.0, are well located in or near hotspots in this 3-D PI forecast. In a relative operating characteristic test (Jolliffe and Stephenson, 2003), the performance of PI forecast is also better than relative intensity (RI).
A brain-computer interface method combined with eye tracking for 3D interaction.
Lee, Eui Chul; Woo, Jin Cheol; Kim, Jong Hwa; Whang, Mincheol; Park, Kang Ryoung
2010-07-15
With the recent increase in the number of three-dimensional (3D) applications, the need for interfaces to these applications has increased. Although the eye tracking method has been widely used as an interaction interface for hand-disabled persons, this approach cannot be used for depth directional navigation. To solve this problem, we propose a new brain computer interface (BCI) method in which the BCI and eye tracking are combined to analyze depth navigation, including selection and two-dimensional (2D) gaze direction, respectively. The proposed method is novel in the following five ways compared to previous works. First, a device to measure both the gaze direction and an electroencephalogram (EEG) pattern is proposed with the sensors needed to measure the EEG attached to a head-mounted eye tracking device. Second, the reliability of the BCI interface is verified by demonstrating that there is no difference between the real and the imaginary movements for the same work in terms of the EEG power spectrum. Third, depth control for the 3D interaction interface is implemented by an imaginary arm reaching movement. Fourth, a selection method is implemented by an imaginary hand grabbing movement. Finally, for the independent operation of gazing and the BCI, a mode selection method is proposed that measures a user's concentration by analyzing the pupil accommodation speed, which is not affected by the operation of gazing and the BCI. According to experimental results, we confirmed the feasibility of the proposed 3D interaction method using eye tracking and a BCI. PMID:20580646
A cut cell method for the 3D simulation of Crookes radiometer
Dechriste, Guillaume; Mieussens, Luc
2014-12-09
Devices involved in engineering applications, such as vacuum pumps or MEMS, may be made of several moving parts. This raise the issue of the simulation of rarefied gas flow around moving boundaries. We propose a simple process, known as cut cell method, to treat the motion of a solid body in the framework of the deterministic solving of a kinetic equation. Up to our knowledge, this is the first time that this approach has been used for this kind of simulations. The method is illustrated by the 2D and 3D simulations of a Crookes radiometer.
Shot noise limit of the optical 3D measurement methods for smooth surfaces
NASA Astrophysics Data System (ADS)
Pavliček, Pavel; Pech, Miroslav
2016-03-01
The measurement uncertainty of optical 3D measurement methods for smooth surfaces caused by shot noise is investigated. The shot noise is a fundamental property of the quantum nature of light. If all noise sources are eliminated, the shot noise represents the ultimate limit of the measurement uncertainty. The measurement uncertainty is calculated for several simple model methods. The analysis shows that the measurement uncertainty depends on the wavelength of used light, the number of photons used for the measurement, and on a factor that is connected with the geometric arrangement of the measurement setup.
GPU accelerated simulations of 3D deterministic particle transport using discrete ordinates method
NASA Astrophysics Data System (ADS)
Gong, Chunye; Liu, Jie; Chi, Lihua; Huang, Haowei; Fang, Jingyue; Gong, Zhenghu
2011-07-01
Graphics Processing Unit (GPU), originally developed for real-time, high-definition 3D graphics in computer games, now provides great faculty in solving scientific applications. The basis of particle transport simulation is the time-dependent, multi-group, inhomogeneous Boltzmann transport equation. The numerical solution to the Boltzmann equation involves the discrete ordinates ( Sn) method and the procedure of source iteration. In this paper, we present a GPU accelerated simulation of one energy group time-independent deterministic discrete ordinates particle transport in 3D Cartesian geometry (Sweep3D). The performance of the GPU simulations are reported with the simulations of vacuum boundary condition. The discussion of the relative advantages and disadvantages of the GPU implementation, the simulation on multi GPUs, the programming effort and code portability are also reported. The results show that the overall performance speedup of one NVIDIA Tesla M2050 GPU ranges from 2.56 compared with one Intel Xeon X5670 chip to 8.14 compared with one Intel Core Q6600 chip for no flux fixup. The simulation with flux fixup on one M2050 is 1.23 times faster than on one X5670.
3D modelling of the electromagnetic response of geophysical targets using the FDTD method
Debroux, P.S.
1996-05-01
A publicly available and maintained electromagnetic finite-difference time domain (FDTD) code has been applied to the forward modelling of the response of 1D, 2D and 3D geophysical targets to a vertical magnetic dipole excitation. The FDTD method is used to analyze target responses in the 1 MHz to 100 MHz range, where either conduction or displacement currents may have the controlling role. The response of the geophysical target to the excitation is presented as changes in the magnetic field ellipticity. The results of the FDTD code compare favorably with previously published integral equation solutions of the response of 1D targets, and FDTD models calculated with different finite-difference cell sizes are compared to find the effect of model discretization on the solution. The discretization errors, calculated as absolute error in ellipticity, are presented for the different ground geometry models considered, and are, for the most part, below 10% of the integral equation solutions. Finally, the FDTD code is used to calculate the magnetic ellipticity response of a 2D survey and a 3D sounding of complicated geophysical targets. The response of these 2D and 3D targets are too complicated to be verified with integral equation solutions, but show the proper low- and high-frequency responses.
GPU accelerated simulations of 3D deterministic particle transport using discrete ordinates method
Gong Chunye; Liu Jie; Chi Lihua; Huang Haowei; Fang Jingyue; Gong Zhenghu
2011-07-01
Graphics Processing Unit (GPU), originally developed for real-time, high-definition 3D graphics in computer games, now provides great faculty in solving scientific applications. The basis of particle transport simulation is the time-dependent, multi-group, inhomogeneous Boltzmann transport equation. The numerical solution to the Boltzmann equation involves the discrete ordinates (S{sub n}) method and the procedure of source iteration. In this paper, we present a GPU accelerated simulation of one energy group time-independent deterministic discrete ordinates particle transport in 3D Cartesian geometry (Sweep3D). The performance of the GPU simulations are reported with the simulations of vacuum boundary condition. The discussion of the relative advantages and disadvantages of the GPU implementation, the simulation on multi GPUs, the programming effort and code portability are also reported. The results show that the overall performance speedup of one NVIDIA Tesla M2050 GPU ranges from 2.56 compared with one Intel Xeon X5670 chip to 8.14 compared with one Intel Core Q6600 chip for no flux fixup. The simulation with flux fixup on one M2050 is 1.23 times faster than on one X5670.
NASA Astrophysics Data System (ADS)
Wang, Minjie; Shen, Shuwei; Yang, Jie; Dong, Erbao; Xu, Ronald
2014-03-01
The performance of biomedical optical imaging devices heavily relies on appropriate calibration. However, many of existing calibration phantoms for biomedical optical devices are based on homogenous materials without considering the multi-layer heterogeneous structures observed in biological tissue. Using such a phantom for optical calibration may result in measurement bias. To overcome this problem, we propose a 3D printing method for freeform fabrication of tissue simulating phantoms with multilayer heterogeneous structure. The phantom simulates not only the morphologic characteristics of biological tissue but also absorption and scattering properties. The printing system is based on a 3D motion platform with coordinated control of the DC motors. A special jet nozzle is designed to mix base, scattering, and absorption materials at different ratios. 3D tissue structures are fabricated through layer-by-layer printing with selective deposition of phantom materials of different ingredients. Different mixed ratios of base, scattering and absorption materials have been tested in order to optimize the printing outcome. A spectrometer and a tissue spectrophotometer are used for characterizing phantom absorption and scattering properties. The goal of this project is to fabricate skin tissue simulating phantoms as a traceable standard for the calibration of biomedical optical spectral devices.
A new method to combine 3D reconstruction volumes for multiple parallel circular cone beam orbits
Baek, Jongduk; Pelc, Norbert J.
2010-01-01
Purpose: This article presents a new reconstruction method for 3D imaging using a multiple 360° circular orbit cone beam CT system, specifically a way to combine 3D volumes reconstructed with each orbit. The main goal is to improve the noise performance in the combined image while avoiding cone beam artifacts. Methods: The cone beam projection data of each orbit are reconstructed using the FDK algorithm. When at least a portion of the total volume can be reconstructed by more than one source, the proposed combination method combines these overlap regions using weighted averaging in frequency space. The local exactness and the noise performance of the combination method were tested with computer simulations of a Defrise phantom, a FORBILD head phantom, and uniform noise in the raw data. Results: A noiseless simulation showed that the local exactness of the reconstructed volume from the source with the smallest tilt angle was preserved in the combined image. A noise simulation demonstrated that the combination method improved the noise performance compared to a single orbit reconstruction. Conclusions: In CT systems which have overlap volumes that can be reconstructed with data from more than one orbit and in which the spatial frequency content of each reconstruction can be calculated, the proposed method offers improved noise performance while keeping the local exactness of data from the source with the smallest tilt angle. PMID:21089770
NASA Astrophysics Data System (ADS)
Chen, Ping-Feng; Krim, Hamid
2008-02-01
In this paper, we propose using two methods to determine the canonical views of 3D objects: minimum description length (MDL) criterion and compressive sensing method. MDL criterion searches for the description length that achieves the balance between model accuracy and parsimony. It takes the form of the sum of a likelihood and a penalizing term, where the likelihood is in favor of model accuracy such that more views assists the description of an object, while the second term penalizes lengthy description to prevent overfitting of the model. In order to devise the likelihood term, we propose a model to represent a 3D object as the weighted sum of multiple range images, which is used in the second method to determine the canonical views as well. In compressive sensing method, an intelligent way of parsimoniously sampling an object is presented. We make direct inference from Donoho1 and Candes'2 work, and adapt it to our model. Each range image is viewed as a projection, or a sample, of a 3D model, and by using compressive sensing theory, we are able to reconstruct the object with an overwhelming probability by scarcely sensing the object in a random manner. Compressive sensing is different from traditional compressing method in the sense that the former compress things in the sampling stage while the later collects a large number of samples and then compressing mechanism is carried out thereafter. Compressive sensing scheme is particularly useful when the number of sensors are limited or the sampling machinery cost much resource or time.
A faster method for 3D/2D medical image registration--a simulation study.
Birkfellner, Wolfgang; Wirth, Joachim; Burgstaller, Wolfgang; Baumann, Bernard; Staedele, Harald; Hammer, Beat; Gellrich, Niels Claudius; Jacob, Augustinus Ludwig; Regazzoni, Pietro; Messmer, Peter
2003-08-21
3D/2D patient-to-computed-tomography (CT) registration is a method to determine a transformation that maps two coordinate systems by comparing a projection image rendered from CT to a real projection image. Iterative variation of the CT's position between rendering steps finally leads to exact registration. Applications include exact patient positioning in radiation therapy, calibration of surgical robots, and pose estimation in computer-aided surgery. One of the problems associated with 3D/2D registration is the fact that finding a registration includes solving a minimization problem in six degrees of freedom (dof) in motion. This results in considerable time requirements since for each iteration step at least one volume rendering has to be computed. We show that by choosing an appropriate world coordinate system and by applying a 2D/2D registration method in each iteration step, the number of iterations can be grossly reduced from n6 to n5. Here, n is the number of discrete variations around a given coordinate. Depending on the configuration of the optimization algorithm, this reduces the total number of iterations necessary to at least 1/3 of it's original value. The method was implemented and extensively tested on simulated x-ray images of a tibia, a pelvis and a skull base. When using one projective image and a discrete full parameter space search for solving the optimization problem, average accuracy was found to be 1.0 +/- 0.6(degrees) and 4.1 +/- 1.9 (mm) for a registration in six parameters, and 1.0 +/- 0.7(degrees) and 4.2 +/- 1.6 (mm) when using the 5 + 1 dof method described in this paper. Time requirements were reduced by a factor 3.1. We conclude that this hardware-independent optimization of 3D/2D registration is a step towards increasing the acceptance of this promising method for a wide number of clinical applications. PMID:12974581
A method for the evaluation of thousands of automated 3D stem cell segmentations.
Bajcsy, P; Simon, M; Florczyk, S J; Simon, C G; Juba, D; Brady, M C
2015-12-01
There is no segmentation method that performs perfectly with any dataset in comparison to human segmentation. Evaluation procedures for segmentation algorithms become critical for their selection. The problems associated with segmentation performance evaluations and visual verification of segmentation results are exaggerated when dealing with thousands of three-dimensional (3D) image volumes because of the amount of computation and manual inputs needed. We address the problem of evaluating 3D segmentation performance when segmentation is applied to thousands of confocal microscopy images (z-stacks). Our approach is to incorporate experimental imaging and geometrical criteria, and map them into computationally efficient segmentation algorithms that can be applied to a very large number of z-stacks. This is an alternative approach to considering existing segmentation methods and evaluating most state-of-the-art algorithms. We designed a methodology for 3D segmentation performance characterization that consists of design, evaluation and verification steps. The characterization integrates manual inputs from projected surrogate 'ground truth' of statistically representative samples and from visual inspection into the evaluation. The novelty of the methodology lies in (1) designing candidate segmentation algorithms by mapping imaging and geometrical criteria into algorithmic steps, and constructing plausible segmentation algorithms with respect to the order of algorithmic steps and their parameters, (2) evaluating segmentation accuracy using samples drawn from probability distribution estimates of candidate segmentations and (3) minimizing human labour needed to create surrogate 'truth' by approximating z-stack segmentations with 2D contours from three orthogonal z-stack projections and by developing visual verification tools. We demonstrate the methodology by applying it to a dataset of 1253 mesenchymal stem cells. The cells reside on 10 different types of biomaterial
Color decomposition method for multiprimary display using 3D-LUT in linearized LAB space
NASA Astrophysics Data System (ADS)
Kang, Dong-Woo; Kim, Yun-Tae; Cho, Yang-Ho; Park, Kee-Hyon; Choe, Wonhee; Ha, Yeong-Ho
2005-01-01
This paper proposes a color decomposition method for a multi-primary display (MPD) using a 3-dimensional look-up-table (3D-LUT) in linearized LAB space. The proposed method decomposes the conventional three primary colors into multi-primary control values for a display device under the constraints of tristimulus matching. To reproduce images on an MPD, the color signals are estimated from a device-independent color space, such as CIEXYZ and CIELAB. In this paper, linearized LAB space is used due to its linearity and additivity in color conversion. First, the proposed method constructs a 3-D LUT containing gamut boundary information to calculate the color signals for the MPD in linearized LAB space. For the image reproduction, standard RGB or CIEXYZ is transformed to linearized LAB, then the hue and chroma are computed with reference to the 3D-LUT. In linearized LAB space, the color signals for a gamut boundary point are calculated to have the same lightness and hue as the input point. Also, the color signals for a point on the gray axis are calculated to have the same lightness as the input point. Based on the gamut boundary points and input point, the color signals for the input point are then obtained using the chroma ratio divided by the chroma of the gamut boundary point. In particular, for a change of hue, the neighboring boundary points are also employed. As a result, the proposed method guarantees color signal continuity and computational efficiency, and requires less memory.
Color decomposition method for multiprimary display using 3D-LUT in linearized LAB space
NASA Astrophysics Data System (ADS)
Kang, Dong-Woo; Kim, Yun-Tae; Cho, Yang-Ho; Park, Kee-Hyon; Choe, Wonhee; Ha, Yeong-Ho
2004-12-01
This paper proposes a color decomposition method for a multi-primary display (MPD) using a 3-dimensional look-up-table (3D-LUT) in linearized LAB space. The proposed method decomposes the conventional three primary colors into multi-primary control values for a display device under the constraints of tristimulus matching. To reproduce images on an MPD, the color signals are estimated from a device-independent color space, such as CIEXYZ and CIELAB. In this paper, linearized LAB space is used due to its linearity and additivity in color conversion. First, the proposed method constructs a 3-D LUT containing gamut boundary information to calculate the color signals for the MPD in linearized LAB space. For the image reproduction, standard RGB or CIEXYZ is transformed to linearized LAB, then the hue and chroma are computed with reference to the 3D-LUT. In linearized LAB space, the color signals for a gamut boundary point are calculated to have the same lightness and hue as the input point. Also, the color signals for a point on the gray axis are calculated to have the same lightness as the input point. Based on the gamut boundary points and input point, the color signals for the input point are then obtained using the chroma ratio divided by the chroma of the gamut boundary point. In particular, for a change of hue, the neighboring boundary points are also employed. As a result, the proposed method guarantees color signal continuity and computational efficiency, and requires less memory.
The image adaptive method for solder paste 3D measurement system
NASA Astrophysics Data System (ADS)
Xiaohui, Li; Changku, Sun; Peng, Wang
2015-03-01
The extensive application of Surface Mount Technology (SMT) requires various measurement methods to evaluate the circuit board. The solder paste 3D measurement system utilizing laser light projecting on the printed circuit board (PCB) surface is one of the critical methods. The local oversaturation, arising from the non-consistent reflectivity of the PCB surface, will lead to inaccurate measurement. The paper reports a novel optical image adaptive method of remedying the local oversaturation for solder paste measurement. The liquid crystal on silicon (LCoS) and image sensor (CCD or CMOS) are combined as the high dynamic range image (HDRI) acquisition system. The significant characteristic of the new method is that the image after adjustment is captured by specially designed HDRI acquisition system programmed by the LCoS mask. The formation of the LCoS mask, depending on a HDRI combined with the image fusion algorithm, is based on separating the laser light from the local oversaturated region. Experimental results demonstrate that the method can significantly improve the accuracy for the solder paste 3D measurement system with local oversaturation.
A 3D finite element ALE method using an approximate Riemann solution
Chiravalle, V. P.; Morgan, N. R.
2016-08-09
Arbitrary Lagrangian–Eulerian finite volume methods that solve a multidimensional Riemann-like problem at the cell center in a staggered grid hydrodynamic (SGH) arrangement have been proposed. This research proposes a new 3D finite element arbitrary Lagrangian–Eulerian SGH method that incorporates a multidimensional Riemann-like problem. Here, two different Riemann jump relations are investigated. A new limiting method that greatly improves the accuracy of the SGH method on isentropic flows is investigated. A remap method that improves upon a well-known mesh relaxation and remapping technique in order to ensure total energy conservation during the remap is also presented. Numerical details and test problemmore » results are presented.« less
A novel 3D absorption correction method for quantitative EDX-STEM tomography.
Burdet, Pierre; Saghi, Z; Filippin, A N; Borrás, A; Midgley, P A
2016-01-01
This paper presents a novel 3D method to correct for absorption in energy dispersive X-ray (EDX) microanalysis of heterogeneous samples of unknown structure and composition. By using STEM-based tomography coupled with EDX, an initial 3D reconstruction is used to extract the location of generated X-rays as well as the X-ray path through the sample to the surface. The absorption correction needed to retrieve the generated X-ray intensity is then calculated voxel-by-voxel estimating the different compositions encountered by the X-ray. The method is applied to a core/shell nanowire containing carbon and oxygen, two elements generating highly absorbed low energy X-rays. Absorption is shown to cause major reconstruction artefacts, in the form of an incomplete recovery of the oxide and an erroneous presence of carbon in the shell. By applying the correction method, these artefacts are greatly reduced. The accuracy of the method is assessed using reference X-ray lines with low absorption. PMID:26484792
A new 3D shape measurement method using digital fringe projection technique
NASA Astrophysics Data System (ADS)
Zhang, Jiarui; Zhang, Yingjie; Yu, Mingrang; Xiang, Dehu
2015-10-01
This paper proposes a novel optical three-dimensional (3D) measurement method using the traditional space-time stereo system. In the proposed method, the projector not only shoots fringe pattern onto the measurement object to achieve precise matching, but also plays a vital role in the 3D information calculation. With the combination of two cameras and a projector, two digital fringe projection (DFP) measurement systems and one traditional space-time stereo measurement system can be obtained. In another word, the measurand will be measured three times simultaneously, which results in three independent point clouds of the same region of the object to be measured. So it is necessary to register these three sets of points for obtaining one final data set. The iterative closest points (ICP) method, which is known as the most popular registration approach, is sensitive to the initial estimation of the transformation between the two sets of points to be matched. Thus, a robust rough registration, which is introduced from Natasha, is useful for ICP to realize accurate registration. After registration, a scattered point set with redundant and errors, which are caused by overlapping, is obtained. Then some local surfaces are constructed for those overlapping regions using the moving least squares (MLS) method, and the points extracted from those surfaces are used to replace the points of the overlapping regions. Finally, a simplified, precise point cloud can be obtained.
Computation of an Underexpanded 3-D Rectangular Jet by the CE/SE Method
NASA Technical Reports Server (NTRS)
Loh, Ching Y.; Himansu, Ananda; Wang, Xiao Y.; Jorgenson, Philip C. E.
2000-01-01
Recently, an unstructured three-dimensional space-time conservation element and solution element (CE/SE) Euler solver was developed. Now it is also developed for parallel computation using METIS for domain decomposition and MPI (message passing interface). The method is employed here to numerically study the near-field of a typical 3-D rectangular under-expanded jet. For the computed case-a jet with Mach number Mj = 1.6. with a very modest grid of 1.7 million tetrahedrons, the flow features such as the shock-cell structures and the axis switching, are in good qualitative agreement with experimental results.
Dynamic Analysis of 2D Electromagnetic Resonant Optical Scanner Using 3D Finite Element Method
NASA Astrophysics Data System (ADS)
Hirata, Katsuhiro; Hong, Sara; Maeda, Kengo
The optical scanner is a scanning device in which a laser beam is reflected by a mirror that can be rotated or oscillated. In this paper, we propose a new 2D electromagnetic resonant optical scanner that employs electromagnets and leaf springs. Torque characteristics and resonance characteristics of the scanner are analyzed using the 3D finite element method. The validity of the analysis is shown by comparing the characteristics inferred from the analysis with the characteristics of the prototype. Further, 2D resonance is investigated by introducing a superimposed-frequency current in a single coil.
On 3-D inelastic analysis methods for hot section components (base program)
NASA Technical Reports Server (NTRS)
Wilson, R. B.; Bak, M. J.; Nakazawa, S.; Banerjee, P. K.
1986-01-01
A 3-D Inelastic Analysis Method program is described. This program consists of a series of new computer codes embodying a progression of mathematical models (mechanics of materials, special finite element, boundary element) for streamlined analysis of: (1) combustor liners, (2) turbine blades, and (3) turbine vanes. These models address the effects of high temperatures and thermal/mechanical loadings on the local (stress/strain)and global (dynamics, buckling) structural behavior of the three selected components. Three computer codes, referred to as MOMM (Mechanics of Materials Model), MHOST (Marc-Hot Section Technology), and BEST (Boundary Element Stress Technology), have been developed and are briefly described in this report.
NASA Technical Reports Server (NTRS)
Wilson, R. B.; Banerjee, P. K.
1987-01-01
This Annual Status Report presents the results of work performed during the third year of the 3-D Inelastic Analysis Methods for Hot Sections Components program (NASA Contract NAS3-23697). The objective of the program is to produce a series of computer codes that permit more accurate and efficient three-dimensional analyses of selected hot section components, i.e., combustor liners, turbine blades, and turbine vanes. The computer codes embody a progression of mathematical models and are streamlined to take advantage of geometrical features, loading conditions, and forms of material response that distinguish each group of selected components.
A hybrid method for the computation of quasi-3D seismograms.
NASA Astrophysics Data System (ADS)
Masson, Yder; Romanowicz, Barbara
2013-04-01
The development of powerful computer clusters and efficient numerical computation methods, such as the Spectral Element Method (SEM) made possible the computation of seismic wave propagation in a heterogeneous 3D earth. However, the cost of theses computations is still problematic for global scale tomography that requires hundreds of such simulations. Part of the ongoing research effort is dedicated to the development of faster modeling methods based on the spectral element method. Capdeville et al. (2002) proposed to couple SEM simulations with normal modes calculation (C-SEM). Nissen-Meyer et al. (2007) used 2D SEM simulations to compute 3D seismograms in a 1D earth model. Thanks to these developments, and for the first time, Lekic et al. (2011) developed a 3D global model of the upper mantle using SEM simulations. At the local and continental scale, adjoint tomography that is using a lot of SEM simulation can be implemented on current computers (Tape, Liu et al. 2009). Due to their smaller size, these models offer higher resolution. They provide us with images of the crust and the upper part of the mantle. In an attempt to teleport such local adjoint tomographic inversions into the deep earth, we are developing a hybrid method where SEM computation are limited to a region of interest within the earth. That region can have an arbitrary shape and size. Outside this region, the seismic wavefield is extrapolated to obtain synthetic data at the Earth's surface. A key feature of the method is the use of a time reversal mirror to inject the wavefield induced by distant seismic source into the region of interest (Robertsson and Chapman 2000). We compute synthetic seismograms as follow: Inside the region of interest, we are using regional spectral element software RegSEM to compute wave propagation in 3D. Outside this region, the wavefield is extrapolated to the surface by convolution with the Green's functions from the mirror to the seismic stations. For now, these
Tadrous, Paul Joseph
2012-01-01
Anatomy has advanced using 3-dimensional (3D) studies at macroscopic (e.g., dissection, injection moulding of vessels, radiology) and microscopic (e.g., serial section reconstruction with light and electron microscopy) levels. This paper presents the first results in human cells of a new method of subcellular 3D brightfield microscopy. Unlike traditional 3D deconvolution and confocal techniques, this method is suitable for general application to brightfield microscopy. Unlike brightfield serial sectioning it has subcellular resolution. Results are presented of the 3D structure of chromatin in the interphase nucleus of two human cell types, hepatocyte and plasma cell. I show how the freedom to examine these structures in 3D allows greater morphological discrimination between and within cell types and the 3D structural basis for the classical “clock-face” motif of the plasma cell nucleus is revealed. Potential for further applications discussed. PMID:22567315
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.
ROI-preserving 3D video compression method utilizing depth information
NASA Astrophysics Data System (ADS)
Ti, Chunli; Xu, Guodong; Guan, Yudong; Teng, Yidan
2015-09-01
Efficiently transmitting the extra information of three dimensional (3D) video is becoming a key issue of the development of 3DTV. 2D plus depth format not only occupies the smaller bandwidth and is compatible transmission under the condition of the existing channel, but also can provide technique support for advanced 3D video compression in some extend. This paper proposes an ROI-preserving compression scheme to further improve the visual quality at a limited bit rate. According to the connection between the focus of Human Visual System (HVS) and depth information, region of interest (ROI) can be automatically selected via depth map progressing. The main improvement from common method is that a meanshift based segmentation is executed to the depth map before foreground ROI selection to keep the integrity of scene. Besides, the sensitive areas along the edges are also protected. The Spatio-temporal filtering adapting to H.264 is used to the non-ROI of both 2D video and depth map before compression. Experiments indicate that, the ROI extracted by this method is more undamaged and according with subjective feeling, and the proposed method can keep the key high-frequency information more effectively while the bit rate is reduced.
3D Multi-spectral Image-guided Near-infrared Spectroscopy using Boundary Element Method
Srinivasan, Subhadra; Pogue, Brian W.; Paulsen, Keith D.
2010-01-01
Image guided (IG) Near-Infrared spectroscopy (NIRS) has the ability to provide high-resolution metabolic and vascular characterization of tissue, with clinical applications in diagnosis of breast cancer. This method is specific to multimodality imaging where tissue boundaries obtained from alternate modalities such as MRI/CT, are used for NIRS recovery. IG-NIRS is severely limited in 3D by challenges such as volumetric meshing of arbitrary anatomical shapes and computational burden encountered by existing models which use finite element method (FEM). We present an efficient and feasible alternative to FEM using boundary element method (BEM). The main advantage is the use of surface discretization which is reliable and more easily generated than volume grids in 3D and enables automation for large number of clinical data-sets. The BEM has been implemented for the diffusion equation to model light propagation in tissue. Image reconstruction based on BEM has been tested in a multi-threading environment using four processors which provides 60% improvement in computational time compared to a single processor. Spectral priors have been implemented in this framework and applied to a three-region problem with mean error of 6% in recovery of NIRS parameters. PMID:21179380
3D modeling for solving forward model of no-contact fluorescence diffuse optical tomography method
NASA Astrophysics Data System (ADS)
Nouizi, F.; Chabrier, R.; Torregrossa, M.; Poulet, P.
2009-07-01
This paper presents detailed computational aspects of a new 3D modeling for solving the direct problem in a no-contact time-resolved Fluorescent Diffuse Optical Tomography (FDOT) method that rely on near-infrared scattered and fluorescent photons to image the optical properties and distribution of fluorescent probes in small laboratory animals. An optical scanner allowing performing in-vivo measurements in no-contact scheme was built in our laboratory and is presented. We use the three-dimensional Finite Element Method (FEM) to solve the coupled diffusion equations of excitation and fluorescence photons in highly scattering objects. The computed results allowed yielding photon density maps and the temporal profiles of photons on the surface of the small animal. Our 3D modeling of propagation of photons in the void space between the surface of the object and the detectors allows calculating the quantity of photons reaching the optodes. Simulations were carried-out on two test objects: a resin cylinder and a mouse phantom. The results demonstrate the potential applications of the method to pre-clinical imaging.
Calculation by the finite element method of 3-D turbulent flow in a centrifugal pump
NASA Astrophysics Data System (ADS)
Combes, J. F.
1992-02-01
In order to solve industrial flow problems in complex geometries, a finite element code, N3S, was developed. It allows the computation of a wide variety of 2-D or 3-D unsteady incompressible flows, by solving the Reynolds averaged Navier-Stokes equations together with a k-epsilon turbulence model. Some recent developments of this code concern turbomachinery flows, where one has to take into account periodic boundary conditions, as well as Coriolis and centrifugal forces. The numerical treatment is based on a fractional step method: at each time step, an advection step is solved successively by means of a characteristic method; a diffusion step for the scalar terms; and finally, a Generalized Stokes Problem by using a preconditioned Uzawa algorithm. The space discretization uses a standard Galerkin finite element method with a mixed formulation for the velocity and pressure. An application is presented of this code to the flow inside a centrifugal pump which was extensively tested on several air and water test rigs, and for which many quasi-3-D or Euler calculations were reported. The present N3S calculation is made on a finite element mesh comprising about 28000 tetrahedrons and 43000 nodes.
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
NASA Astrophysics Data System (ADS)
Kochi, Nobuo; Ito, Tadayuki; Kitamura, Kazuo; Kaneko, Syun'ichi
The three dimensional measurement & modeling system with digital cameras on PC is now making progress and its need and hope is increasingly felt in terrestrial (close-range) photogrammetry for such sectors as cultural heritage preservation, architecture, civil engineering, manufacturing, measurement etc. Therefore, we have developed a system to improve the accuracy of stereo-matching, which is the very core of 3D measurement. As for stereo-matching method, in order to minimize the mismatching and to be robust in geometric distortions, occlusion, as well as brightness change, we invented Coarse-to-Fine Strategy Method by integrating OCM (Orientation Code Matching) with LSM (Least Squares Matching). Thus this system could attain the accuracy of 0.26mm, when we experimented on a mannequin. And when we actually experimented on the archeological ruins in Greece and Turkey, the accuracy was within the range of 1cm, compared with their blue-print plan. Besides, formally workers used to take at least 1.5 month for this kind of survey operation with the existing method, but now workers need only 3 or 4 days. Thus, its practicality and efficiency was confirmed. This paper demonstrates our new system of 3D measurement and stereo-matching with some concrete examples as its practical application.
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.
A least-squares finite element method for 3D incompressible Navier-Stokes equations
NASA Technical Reports Server (NTRS)
Jiang, Bo-Nan; Lin, T. L.; Hou, Lin-Jun; Povinelli, Louis A.
1993-01-01
The least-squares finite element method (LSFEM) based on the velocity-pressure-vorticity formulation is applied to three-dimensional steady incompressible Navier-Stokes problems. This method can accommodate equal-order interpolations, and results in symmetric, positive definite algebraic system. An additional compatibility equation, i.e., the divergence of vorticity vector should be zero, is included to make the first-order system elliptic. The Newton's method is employed to linearize the partial differential equations, the LSFEM is used to obtain discretized equations, and the system of algebraic equations is solved using the Jacobi preconditioned conjugate gradient method which avoids formation of either element or global matrices (matrix-free) to achieve high efficiency. The flow in a half of 3D cubic cavity is calculated at Re = 100, 400, and 1,000 with 50 x 52 x 25 trilinear elements. The Taylor-Gortler-like vortices are observed at Re = 1,000.
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
Multilevel local refinement and multigrid methods for 3-D turbulent flow
Liao, C.; Liu, C.; Sung, C.H.; Huang, T.T.
1996-12-31
A numerical approach based on multigrid, multilevel local refinement, and preconditioning methods for solving incompressible Reynolds-averaged Navier-Stokes equations is presented. 3-D turbulent flow around an underwater vehicle is computed. 3 multigrid levels and 2 local refinement grid levels are used. The global grid is 24 x 8 x 12. The first patch is 40 x 16 x 20 and the second patch is 72 x 32 x 36. 4th order artificial dissipation are used for numerical stability. The conservative artificial compressibility method are used for further improvement of convergence. To improve the accuracy of coarse/fine grid interface of local refinement, flux interpolation method for refined grid boundary is used. The numerical results are in good agreement with experimental data. The local refinement can improve the prediction accuracy significantly. The flux interpolation method for local refinement can keep conservation for a composite grid, therefore further modify the prediction accuracy.
Accurate, finite-volume methods for 3D MHD on unstructured Lagrangian meshes
Barnes, D.C.; Rousculp, C.L.
1998-10-01
Previous 2D methods for magnetohydrodynamics (MHD) have contributed both to development of core code capability and to physics applications relevant to AGEX pulsed-power experiments. This strategy is being extended to 3D by development of a modular extension of an ASCI code. Extension to 3D not only increases complexity by problem size, but also introduces new physics, such as magnetic helicity transport. The authors have developed a method which incorporates all known conservation properties into the difference scheme on a Lagrangian unstructured mesh. Because the method does not depend on the mesh structure, mesh refinement is possible during a calculation to prevent the well known problem of mesh tangling. Arbitrary polyhedral cells are decomposed into tetrahedrons. The action of the magnetic vector potential, A {center_dot} {delta}l, is centered on the edges of this extended mesh. For ideal flow, this maintains {del} {center_dot} B = 0 to round-off error. Vertex forces are derived by the variation of magnetic energy with respect to vertex positions, F = {minus}{partial_derivative}W{sub B}/{partial_derivative}r. This assures symmetry as well as magnetic flux, momentum, and energy conservation. The method is local so that parallelization by domain decomposition is natural for large meshes. In addition, a simple, ideal-gas, finite pressure term has been included. The resistive diffusion part is calculated using the support operator method, to obtain an energy conservative, symmetric method on an arbitrary mesh. Implicit time difference equations are solved by preconditioned, conjugate gradient methods. Results of convergence tests are presented. Initial results of an annular Z-pinch implosion problem illustrate the application of these methods to multi-material problems.
A novel window based method for approximating the Hausdorff in 3D range imagery.
Koch, Mark William
2004-10-01
Matching a set of 3D points to another set of 3D points is an important part of any 3D object recognition system. The Hausdorff distance is known for it robustness in the face of obscuration, clutter, and noise. We show how to approximate the 3D Hausdorff fraction with linear time complexity and quadratic space complexity. We empirically demonstrate that the approximation is very good when compared to actual Hausdorff distances.
A comparative study for 2D and 3D computer-aided diagnosis methods for solitary pulmonary nodules.
Yeh, Chinson; Wang, Jen-Feng; Wu, Ming-Ting; Yen, Chen-Wen; Nagurka, Mark L; Lin, Chen-Liang
2008-06-01
Many computer-aided diagnosis (CAD) methods, including 2D and 3D approaches, have been proposed for solitary pulmonary nodules (SPNs). However, the detection and diagnosis of SPNs remain challenging in many clinical circumstances. One goal of this work is to investigate the relative diagnostic accuracy of 2D and 3D methods. An additional goal is to develop a two-stage approach that combines the simplicity of 2D and the accuracy of 3D methods. The experimental results show statistically significant differences between the diagnostic accuracy of 2D and 3D methods. The results also show that with a very minor drop in diagnostic performance the two-stage approach can significantly reduce the number of nodules needed to be processed by the 3D method, streamlining the computational demand. PMID:18313899
NASA Technical Reports Server (NTRS)
Zhao, W.; Newman, J. C., Jr.; Sutton, M. A.; Shivakumar, K. N.; Wu, X. R.
1995-01-01
Parallel with the work in Part-1, stress intensity factors for semi-elliptical surface cracks emanating from a circular hole are determined. The 3-D weight function method with the 3D finite element solutions for the uncracked stress distribution as in Part-1 is used for the analysis. Two different loading conditions, i.e. remote tension and wedge loading, are considered for a wide range in geometrical parameters. Both single and double surface cracks are studied and compared with other solutions available in the literature. Typical crack opening displacements are also provided.
Determination of 3D Surface Roughness Parameters by Cross-Section Method
NASA Astrophysics Data System (ADS)
Rudzitis, J.; Krizbergs, J.; Kumermanis, M.; Mozga, N.; Ancans, A.; Leitans, A.
2014-04-01
Currently, in the production engineering the surface roughness parameters are estimated in three dimensions, however, the equipment for these measurements is rather expensive and not always available. In many cases to buy such equipment is not economically justified. Therefore, the 3D surface roughness parameters are usually determined from the well-known 2D profile ones using the existing 2D equipment. This could be done best using the cross-section (or profile) method, especially in the case of nanoroughness estimation, with calculation of the mean values for the roughness height, spacing, and shape. This method - though mainly meant for irregular rough surfaces - can also be used for other types of rough surfaces. Particular emphasis is here given to the correlation between the surface cross-section (profile) parameters and 3D parameters as well as to the choice of the number of cross-cuttings and their orientation on the surface. Mūsdienu ražošanā ir nepieciešams novērtēt virsmas raupjuma parametrus trijās dimensijās, tomēr, aprīkojums šādu mērījumu veikšanai ir ļoti dārgs un ne vienmēr pieejams. Tādēļ bieži rodas nepieciešamība noteikt 3D virsmas raupjuma parametrus pēc labi zināmajiem profila (2D) parametriem, izmantojot eksistējošo 2D mērīšanas aprīkojumu. Labākais risinājums šai problēmai ir izmantot 3D raupjuma parametru noteikšanai šķēlumu jeb profilu metodi. Metode uzrāda labus rezultātus arī novērtējot nanoraupjumu. Iespējams aprēķināt sekojošu virsmas raupjuma mikrotopogrāfisko parametru vidējās vērtības: raupjuma augstumu; soļu parametrus un formu. Metode ir paredzēta izmantošanai virsmām ar neregulāru raksturu, bet var tikt pielāgota arī citu tipu virsmām.
A novel method for vaginal cylinder treatment planning: a seamless transition to 3D brachytherapy
Wu, Vincent; Wang, Zhou; Patil, Sachin
2012-01-01
Purpose Standard treatment plan libraries are often used to ensure a quick turn-around time for vaginal cylinder treatments. Recently there is increasing interest in transitioning from conventional 2D radiograph based brachytherapy to 3D image based brachytherapy, which has resulted in a substantial increase in treatment planning time and decrease in patient through-put. We describe a novel technique that significantly reduces the treatment planning time for CT-based vaginal cylinder brachytherapy. Material and methods Oncentra MasterPlan TPS allows multiple sets of data points to be classified as applicator points which has been harnessed in this method. The method relies on two hard anchor points: the first dwell position in a catheter and an applicator configuration specific dwell position as the plan origin and a soft anchor point beyond the last active dwell position to define the axis of the catheter. The spatial location of various data points on the applicator's surface and at 5 mm depth are stored in an Excel file that can easily be transferred into a patient CT data set using window operations and then used for treatment planning. The remainder of the treatment planning process remains unaffected. Results The treatment plans generated on the Oncentra MasterPlan TPS using this novel method yielded results comparable to those generated on the Plato TPS using a standard treatment plan library in terms of treatment times, dwell weights and dwell times for a given optimization method and normalization points. Less than 2% difference was noticed between the treatment times generated between both systems. Using the above method, the entire planning process, including CT importing, catheter reconstruction, multiple data point definition, optimization and dose prescription, can be completed in ~5–10 minutes. Conclusion The proposed method allows a smooth and efficient transition to 3D CT based vaginal cylinder brachytherapy planning. PMID:23349650
Spectral triangulation: a 3D method for locating single-walled carbon nanotubes in vivo.
Lin, Ching-Wei; Bachilo, Sergei M; Vu, Michael; Beckingham, Kathleen M; Bruce Weisman, R
2016-05-21
Nanomaterials with luminescence in the short-wave infrared (SWIR) region are of special interest for biological research and medical diagnostics because of favorable tissue transparency and low autofluorescence backgrounds in that region. Single-walled carbon nanotubes (SWCNTs) show well-known sharp SWIR spectral signatures and therefore have potential for noninvasive detection and imaging of cancer tumours, when linked to selective targeting agents such as antibodies. However, such applications face the challenge of sensitively detecting and localizing the source of SWIR emission from inside tissues. A new method, called spectral triangulation, is presented for three dimensional (3D) localization using sparse optical measurements made at the specimen surface. Structurally unsorted SWCNT samples emitting over a range of wavelengths are excited inside tissue phantoms by an LED matrix. The resulting SWIR emission is sampled at points on the surface by a scanning fibre optic probe leading to an InGaAs spectrometer or a spectrally filtered InGaAs avalanche photodiode detector. Because of water absorption, attenuation of the SWCNT fluorescence in tissues is strongly wavelength-dependent. We therefore gauge the SWCNT-probe distance by analysing differential changes in the measured SWCNT emission spectra. SWCNT fluorescence can be clearly detected through at least 20 mm of tissue phantom, and the 3D locations of embedded SWCNT test samples are found with sub-millimeter accuracy at depths up to 10 mm. Our method can also distinguish and locate two embedded SWCNT sources at distinct positions. PMID:27140495
A method of 3D object recognition and localization in a cloud of points
NASA Astrophysics Data System (ADS)
Bielicki, Jerzy; Sitnik, Robert
2013-12-01
The proposed method given in this article is prepared for analysis of data in the form of cloud of points directly from 3D measurements. It is designed for use in the end-user applications that can directly be integrated with 3D scanning software. The method utilizes locally calculated feature vectors (FVs) in point cloud data. Recognition is based on comparison of the analyzed scene with reference object library. A global descriptor in the form of a set of spatially distributed FVs is created for each reference model. During the detection process, correlation of subsets of reference FVs with FVs calculated in the scene is computed. Features utilized in the algorithm are based on parameters, which qualitatively estimate mean and Gaussian curvatures. Replacement of differentiation with averaging in the curvatures estimation makes the algorithm more resistant to discontinuities and poor quality of the input data. Utilization of the FV subsets allows to detect partially occluded and cluttered objects in the scene, while additional spatial information maintains false positive rate at a reasonably low level.
Liu, Miao; Yang, Shourui; Wang, Zhangying; Huang, Shujun; Liu, Yue; Niu, Zhenqi; Zhang, Xiaoxuan; Zhu, Jigui; Zhang, Zonghua
2016-05-30
Augmented reality system can be applied to provide precise guidance for various kinds of manual works. The adaptability and guiding accuracy of such systems are decided by the computational model and the corresponding calibration method. In this paper, a novel type of augmented reality guiding system and the corresponding designing scheme are proposed. Guided by external positioning equipment, the proposed system can achieve high relative indication accuracy in a large working space. Meanwhile, the proposed system is realized with a digital projector and the general back projection model is derived with geometry relationship between digitized 3D model and the projector in free space. The corresponding calibration method is also designed for the proposed system to obtain the parameters of projector. To validate the proposed back projection model, the coordinate data collected by a 3D positioning equipment is used to calculate and optimize the extrinsic parameters. The final projecting indication accuracy of the system is verified with subpixel pattern projecting technique. PMID:27410124
A Parallelized 3D Particle-In-Cell Method With Magnetostatic Field Solver And Its Applications
NASA Astrophysics Data System (ADS)
Hsu, Kuo-Hsien; Chen, Yen-Sen; Wu, Men-Zan Bill; Wu, Jong-Shinn
2008-10-01
A parallelized 3D self-consistent electrostatic particle-in-cell finite element (PIC-FEM) code using an unstructured tetrahedral mesh was developed. For simulating some applications with external permanent magnet set, the distribution of the magnetostatic field usually also need to be considered and determined accurately. In this paper, we will firstly present the development of a 3D magnetostatic field solver with an unstructured mesh for the flexibility of modeling objects with complex geometry. The vector Poisson equation for magnetostatic field is formulated using the Galerkin nodal finite element method and the resulting matrix is solved by parallel conjugate gradient method. A parallel adaptive mesh refinement module is coupled to this solver for better resolution. Completed solver is then verified by simulating a permanent magnet array with results comparable to previous experimental observations and simulations. By taking the advantage of the same unstructured grid format of this solver, the developed PIC-FEM code could directly and easily read the magnetostatic field for particle simulation. In the upcoming conference, magnetron is simulated and presented for demonstrating the capability of this code.
Sutradhar, Alok; Park, Jaejong; Carrau, Diana; Nguyen, Tam H; Miller, Michael J; Paulino, Glaucio H
2016-07-01
Large craniofacial defects require efficient bone replacements which should not only provide good aesthetics but also possess stable structural function. The proposed work uses a novel multiresolution topology optimization method to achieve the task. Using a compliance minimization objective, patient-specific bone replacement shapes can be designed for different clinical cases that ensure revival of efficient load transfer mechanisms in the mid-face. In this work, four clinical cases are introduced and their respective patient-specific designs are obtained using the proposed method. The optimized designs are then virtually inserted into the defect to visually inspect the viability of the design . Further, once the design is verified by the reconstructive surgeon, prototypes are fabricated using a 3D printer for validation. The robustness of the designs are mechanically tested by subjecting them to a physiological loading condition which mimics the masticatory activity. The full-field strain result through 3D image correlation and the finite element analysis implies that the solution can survive the maximum mastication of 120 lb. Also, the designs have the potential to restore the buttress system and provide the structural integrity. Using the topology optimization framework in designing the bone replacement shapes would deliver surgeons new alternatives for rather complicated mid-face reconstruction. PMID:26660897
Microelectro discharge machining: an innovative method for the fabrication of 3D microdevices
NASA Astrophysics Data System (ADS)
Lesche, Claudia; Krah, Thomas; Büttgenbach, Stephanus
2011-06-01
This paper reports on the potential of microelectro discharge machining (μEDM) as an innovative method for the fabrication of 3D microdevices. To demonstrate the wide capabilities of μEDM two different high-potential 3D microsystems - a microfluidic device for the dispersion of nanoparticles and a star probe for microcoordinate metrology - are presented. For the fabrication of these microdevices a μEDM-milling machine with integrated microwire electro discharge grinding (μWEDG) module is utilized. To gain optimized process conditions as well as a high surface quality an adequate adaption of the single erosion parameters such as energy, pulse frequency and spark gap has to be carried out and are discussed below. The dispersion micromodule is used for pharmaceutical screening applications in a high pressure range up to 2000 bar. At the channel bottom a surface roughness of Ra = 80 nm is achieved. In case of the star probe it is possible to produce shaft and sphere out of one piece. The fabricated stylus elements have sphere diameters of 40-200 μm. For both applications μEDM offers a flexible, precise, effective and cost-efficient fabrication method for the machining of hard and resistant materials.
A coordinate-free method for the analysis of 3D facial change
NASA Astrophysics Data System (ADS)
Mao, Zhili; Siebert, Jan Paul; Cockshott, W. Paul; Ayoub, Ashraf Farouk
2004-05-01
Euclidean Distance Matrix Analysis (EDMA) is widely held as the most important coordinate-free method by which to analyze landmarks. It has been used extensively in the field of medical anthropometry and has already produced many useful results. Unfortunately this method renders little information regarding the surface on which these points are located and accordingly is inadequate for the 3D analysis of surface anatomy. Here we shall present a new inverse surface flatness metric, the ratio between the Geodesic and the Euclidean inter-landmark distances. Because this metric also only reflects one aspect of three-dimensional shape, i.e. surface flatness, we have combined it with the Euclidean distance to investigate 3D facial change. The goal of this investigation is to be able to analyze three-dimensional facial change in terms of bilateral symmetry as encoded both by surface flatness and by geometric configuration. Our initial study, based on 25 models of surgically managed children (unilateral cleft lip repair) and 40 models of control children at the age of 2 years, indicates that the faces of the surgically managed group were found to be significantly less symmetric than those of the control group in terms of surface flatness, geometric configuration and overall symmetry.
Improvement of advanced nodal method used in 3D core design system
Rauck, S.; Dall'Osso, A.
2006-07-01
This paper deals with AREVA NP progress in the modelling of neutronic phenomena, evaluated through 3D determinist core codes and using 2-group diffusion theory. Our report highlights the advantages of taking into account the assembly environment in the process used for the building of the 2-group collapsed neutronic parameters, such as cross sections or discontinuity factors. The interest of the present method, developed in order to account for the impact of the environment on the above mentioned parameters, resides (i) in the very definition of a global correlation between collapsed neutronic data calculated in an infinite medium and those calculated in a 3D-geometry, and (ii) in the use of a re-homogenization method. Using this approach, computations match better with actual measurements on control rod worth. They also present smaller differences on pin by pin power values compared to the ones computed with another code considered as a reference since it relies on multigroup transport theory. (authors)
Gietzelt, Matthias; Wolf, Klaus-Hendrik; Marschollek, Michael; Haux, Reinhold
2013-07-01
Calibration of accelerometers can be reduced to 3D-ellipsoid fitting problems. Changing extrinsic factors like temperature, pressure or humidity, as well as intrinsic factors like the battery status, demand to calibrate the measurements permanently. Thus, there is a need for fast calibration algorithms, e.g. for online analyses. The primary aim of this paper is to propose a non-iterative calibration algorithm for accelerometers with the focus on minimal execution time and low memory consumption. The secondary aim is to benchmark existing calibration algorithms based on 3D-ellipsoid fitting methods. We compared the algorithms regarding the calibration quality and the execution time as well as the number of quasi-static measurements needed for a stable calibration. As evaluation criterion for the calibration, both the norm of calibrated real-life measurements during inactivity and simulation data was used. The algorithms showed a high calibration quality, but the execution time differed significantly. The calibration method proposed in this paper showed the shortest execution time and a very good performance regarding the number of measurements needed to produce stable results. Furthermore, this algorithm was successfully implemented on a sensor node and calibrates the measured data on-the-fly while continuously storing the measured data to a microSD-card. PMID:23566707
Advanced methods for 3-D inelastic structural analysis for hot engine structures
NASA Technical Reports Server (NTRS)
Chamis, C. C.
1989-01-01
Three-dimensional Inelastic Analysis Methods are described. These methods were incorporated into a series of new computer codes embodying a progression of mathematical models (mechanics of materials, specialty finite element, boundary element) for streamlined analysis of hot engine structures such as: (1) combustor liners, (2) turbine blades, and (3) turbine vanes. These models address the effects of high temperatures and thermal/mechanical loadings on the local (stress/strain) and global (displacements, frequencies, amplitudes, buckling) structural behavior of the three respective components. The methods and the three computer codes, referred to as MOMM (Mechanics Of Materials Model), MHOST (MARC-Hot Section Technology), and BEST3D (Boundary Element Stress Technology), have been developed and are briefly described.
A multi-thread scheduling method for 3D CT image reconstruction using multi-GPU.
Zhu, Yining; Zhao, Yunsong; Zhao, Xing
2012-01-01
As a whole process, we present a concept that the complete reconstruction of CT image should include the computation part on GPUs and the data storage part on hard disks. From this point of view, we propose a Multi-Thread Scheduling (MTS) method to implement the 3D CT image reconstruction such as using FDK algorithm, to trade off the computing and storage time. In this method we use Multi-Threads to control GPUs and a separate thread to accomplish data storage, so that we make the calculation and data storage simultaneously. In addition, we use the 4-channel texture to maintain symmetrical projection data in CUDA framework, which can reduce the calculation time significantly. Numerical experiment shows that the time for the whole process with our method is almost the same as the data storage time. PMID:22635174
Block-Iterative Methods for 3D Constant-Coefficient Stencils on GPUs and Multicore CPUs
Philip, Bobby; Wang, Zhen; Berrill, Mark A
2014-06-01
Block iterative methods are extremely important as smoothers for multigrid methods, as preconditioners for Krylov methods, and as solvers for diagonally dominant linear systems. Developing robust and efficient smoother algorithms suitable for current and evolving GPU and multicore CPU systems is a significant challenge. We address this issue in the case of constant-coefficient stencils arising in the solution of elliptic partial differential equations on structured 3D uniform and adaptively refined block structured grids. Robust, highly parallel implementations of block Jacobi and chaotic block Gauss-Seidel algorithms with exact inversion of the blocks are developed using different parallelization techniques. Experimental results for NVIDIA Fermi/Kepler GPUs and AMD multicore systems are presented.
3D shape reconstruction of medical images using a perspective shape-from-shading method
NASA Astrophysics Data System (ADS)
Yang, Lei; Han, Jiu-qiang
2008-06-01
A 3D shape reconstruction approach for medical images using a shape-from-shading (SFS) method was proposed in this paper. A new reflectance map equation of medical images was analyzed with the assumption that the Lambertian reflectance surface was irradiated by a point light source located at the light center and the image was formed under perspective projection. The corresponding static Hamilton-Jacobi (H-J) equation of the reflectance map equation was established. So the shape-from-shading problem turned into solving the viscosity solution of the static H-J equation. Then with the conception of a viscosity vanishing approximation, the Lax-Friedrichs fast sweeping numerical method was used to compute the viscosity solution of the H-J equation and a new iterative SFS algorithm was gained. Finally, experiments on both synthetic images and real medical images were performed to illustrate the efficiency of the proposed SFS method.
A Tool-Free Calibration Method for Turntable-Based 3D Scanning Systems.
Pang, Xufang; Lau, Rynson W H; Song, Zhan; Li, Yangyan; He, Shengfeng
2016-01-01
Turntable-based 3D scanners are popular but require calibration of the turntable axis. Existing methods for turntable calibration typically make use of specially designed tools, such as a chessboard or criterion sphere, which users must manually install and dismount. In this article, the authors propose an automatic method to calibrate the turntable axis without any calibration tools. Given a scan sequence of the input object, they first recover the initial rotation axis from an automatic registration step. Then they apply an iterative procedure to obtain the optimized turntable axis. This iterative procedure alternates between two steps: refining the initial pose of the input scans and approximating the rotation matrix. The performance of the proposed method was evaluated on a structured light-based scanning system. PMID:25137724
Review on applications of 3D inverse design method for pump
NASA Astrophysics Data System (ADS)
Yin, Junlian; Wang, Dezhong
2014-05-01
The 3D inverse design method, which methodology is far superior to the conventional design method that based on geometrical description, is gradually applied in pump blade design. However, no complete description about the method is outlined. Also, there are no general rules available to set the two important input parameters, blade loading distribution and stacking condition. In this sense, the basic theory and the mechanism why the design method can suppress the formation of secondary flow are summarized. And also, several typical pump design cases with different specific speeds ranging from centrifugal pump to axial pump are surveyed. The results indicates that, for centrifugal pump and mixed pump or turbine, the ratio of blade loading on the hub to that on the shroud is more than unit in the fore part of the blade, whereas in the aft part, the ratio is decreased to satisfy the same wrap angle for hub and shroud. And the choice of blade loading type depends on the balancing of efficiency and cavitation. If the cavitation is more weighted, the better choice is aft-loaded, otherwise, the fore-loaded or mid-loaded is preferable to improve the efficiency. The stacking condition, which is an auxiliary to suppress the secondary flow, can have great effect on the jet-wake outflow and the operation range for pump. Ultimately, how to link the design method to modern optimization techniques is illustrated. With the know-how design methodology and the know-how systematic optimization approach, the application of optimization design is promising for engineering. This paper summarizes the 3D inverse design method systematically.
An unstaggered constrained transport method for the 3D ideal magnetohydrodynamic equations
NASA Astrophysics Data System (ADS)
Helzel, Christiane; Rossmanith, James A.; Taetz, Bertram
2011-05-01
Numerical methods for solving the ideal magnetohydrodynamic (MHD) equations in more than one space dimension must either confront the challenge of controlling errors in the discrete divergence of the magnetic field, or else be faced with nonlinear numerical instabilities. One approach for controlling the discrete divergence is through a so-called constrained transport method, which is based on first predicting a magnetic field through a standard finite volume solver, and then correcting this field through the appropriate use of a magnetic vector potential. In this work we develop a constrained transport method for the 3D ideal MHD equations that is based on a high-resolution wave propagation scheme. Our proposed scheme is the 3D extension of the 2D scheme developed by Rossmanith [J.A. Rossmanith, An unstaggered, high-resolution constrained transport method for magnetohydrodynamic flows, SIAM J. Sci. Comput. 28 (2006) 1766], and is based on the high-resolution wave propagation method of Langseth and LeVeque [J.O. Langseth, R.J. LeVeque, A wave propagation method for threedimensional hyperbolic conservation laws, J. Comput. Phys. 165 (2000) 126]. In particular, in our extension we take great care to maintain the three most important properties of the 2D scheme: (1) all quantities, including all components of the magnetic field and magnetic potential, are treated as cell-centered; (2) we develop a high-resolution wave propagation scheme for evolving the magnetic potential; and (3) we develop a wave limiting approach that is applied during the vector potential evolution, which controls unphysical oscillations in the magnetic field. One of the key numerical difficulties that is novel to 3D is that the transport equation that must be solved for the magnetic vector potential is only weakly hyperbolic. In presenting our numerical algorithm we describe how to numerically handle this problem of weak hyperbolicity, as well as how to choose an appropriate gauge condition. The
Quantitative validation of the 3D SAR profile of hyperthermia applicators using the gamma method.
de Bruijne, Maarten; Samaras, Theodoros; Chavannes, Nicolas; van Rhoon, Gerard C
2007-06-01
For quality assurance of hyperthermia treatment planning systems, quantitative validation of the electromagnetic model of an applicator is essential. The objective of this study was to validate a finite-difference time-domain (FDTD) model implementation of the Lucite cone applicator (LCA) for superficial hyperthermia. The validation involved (i) the assessment of the match between the predicted and measured 3D specific absorption rate (SAR) distribution, and (ii) the assessment of the ratio between model power and real-world power. The 3D SAR distribution of seven LCAs was scanned in a phantom bath using the DASY4 dosimetric measurement system. The same set-up was modelled in SEMCAD X. The match between the predicted and the measured SAR distribution was quantified with the gamma method, which combines distance-to-agreement and dose difference criteria. Good quantitative agreement was observed: more than 95% of the measurement points met the acceptance criteria 2 mm/2% for all applicators. The ratio between measured and predicted power absorption ranged from 0.75 to 0.92 (mean 0.85). This study shows that quantitative validation of hyperthermia applicator models is feasible and is worth considering as a part of hyperthermia quality assurance procedures. PMID:17505090
Women's Preferences for Penis Size: A New Research Method Using Selection among 3D Models
Park, Jaymie; Leung, Shannon
2015-01-01
Women’s preferences for penis size may affect men’s comfort with their own bodies and may have implications for sexual health. Studies of women’s penis size preferences typically have relied on their abstract ratings or selecting amongst 2D, flaccid images. This study used haptic stimuli to allow assessment of women’s size recall accuracy for the first time, as well as examine their preferences for erect penis sizes in different relationship contexts. Women (N = 75) selected amongst 33, 3D models. Women recalled model size accurately using this method, although they made more errors with respect to penis length than circumference. Women preferred a penis of slightly larger circumference and length for one-time (length = 6.4 inches/16.3 cm, circumference = 5.0 inches/12.7 cm) versus long-term (length = 6.3 inches/16.0 cm, circumference = 4.8 inches/12.2 cm) sexual partners. These first estimates of erect penis size preferences using 3D models suggest women accurately recall size and prefer penises only slightly larger than average. PMID:26332467
Method for accurate sizing of pulmonary vessels from 3D medical images
NASA Astrophysics Data System (ADS)
O'Dell, Walter G.
2015-03-01
Detailed characterization of vascular anatomy, in particular the quantification of changes in the distribution of vessel sizes and of vascular pruning, is essential for the diagnosis and management of a variety of pulmonary vascular diseases and for the care of cancer survivors who have received radiation to the thorax. Clinical estimates of vessel radii are typically based on setting a pixel intensity threshold and counting how many "On" pixels are present across the vessel cross-section. A more objective approach introduced recently involves fitting the image with a library of spherical Gaussian filters and utilizing the size of the best matching filter as the estimate of vessel diameter. However, both these approaches have significant accuracy limitations including mis-match between a Gaussian intensity distribution and that of real vessels. Here we introduce and demonstrate a novel approach for accurate vessel sizing using 3D appearance models of a tubular structure along a curvilinear trajectory in 3D space. The vessel branch trajectories are represented with cubic Hermite splines and the tubular branch surfaces represented as a finite element surface mesh. An iterative parameter adjustment scheme is employed to optimally match the appearance models to a patient's chest X-ray computed tomography (CT) scan to generate estimates for branch radii and trajectories with subpixel resolution. The method is demonstrated on pulmonary vasculature in an adult human CT scan, and on 2D simulated test cases.
NASA Astrophysics Data System (ADS)
Lin, Zhili; Li, Xiaoyan; Zhao, Kuixia; Chen, Xudong; Chen, Mingyu; Pu, Jixiong
2016-06-01
For an inertial confinement fusion (ICF) system, the light intensity distribution in the hohlraum is key to the initial plasma excitation and later laser-plasma interaction process. Based on the concept of coordinate transformation of spatial points and vector, we present a robust method with a detailed procedure that makes the calculation of the three dimensional (3D) light intensity distribution in hohlraum easily. The method is intuitive but powerful enough to solve the complex cases of random number of laser beams with arbitrary polarization states and incidence angles. Its application is exemplified in the Shenguang III Facility (SG-III) that verifies its effectiveness and it is useful for guiding the design of hohlraum structure parameter.
Segmentation of Brain MRI Using SOM-FCM-Based Method and 3D Statistical Descriptors
Ortiz, Andrés; Palacio, Antonio A.; Górriz, Juan M.; Ramírez, Javier; Salas-González, Diego
2013-01-01
Current medical imaging systems provide excellent spatial resolution, high tissue contrast, and up to 65535 intensity levels. Thus, image processing techniques which aim to exploit the information contained in the images are necessary for using these images in computer-aided diagnosis (CAD) systems. Image segmentation may be defined as the process of parcelling the image to delimit different neuroanatomical tissues present on the brain. In this paper we propose a segmentation technique using 3D statistical features extracted from the volume image. In addition, the presented method is based on unsupervised vector quantization and fuzzy clustering techniques and does not use any a priori information. The resulting fuzzy segmentation method addresses the problem of partial volume effect (PVE) and has been assessed using real brain images from the Internet Brain Image Repository (IBSR). PMID:23762192
A new method using orthogonal two-frequency grating in online 3D measurement
NASA Astrophysics Data System (ADS)
Peng, Kuang; Cao, Yiping; Wu, Yingchun; Lu, Mingteng
2016-09-01
In online 3D measurement, a new method using orthogonal two-frequency grating based on Phase Measuring Profilometry(PMP) is proposed. The modulation of the entire measured object is used to match pixels and this proposed method successfully resolves the contradiction of the demand for different frequency fringes between the extraction of the modulation information and the phase unwrapping. The high-frequency fringe is used to catch the better modulation patterns for pixel matching, and the low-frequency fringe is used to calculate the phase distribution and avoid phase unwrapping error. In addition, to extract the better modulation patterns for pixel matching, the flat filtering window replaces the circular filtering window to avoid the spectrum aliasing phenomenon. The simulations and experiments show its feasibility.
3D-DIP-Chip: a microarray-based method to measure genomic DNA damage
Powell, James Rees; Bennett, Mark Richard; Evans, Katie Ellen; Yu, Shirong; Webster, Richard Michael; Waters, Raymond; Skinner, Nigel; Reed, Simon Huw
2015-01-01
Genotoxins cause DNA damage, which can result in genomic instability. The genetic changes induced have far-reaching consequences, often leading to diseases such as cancer. A wide range of genotoxins exists, including radiations and chemicals found naturally in the environment, and in man-made forms created by human activity across a variety of industries. Genomic technologies offer the possibility of unravelling the mechanisms of genotoxicity, including the repair of genetic damage, enhancing our ability to develop, test and safely use existing and novel materials. We have developed 3D-DIP-Chip, a microarray-based method to measure the prevalence of genomic genotoxin-induced DNA damage. We demonstrate the measurement of both physical and chemical induced DNA damage spectra, integrating the analysis of these with the associated changes in histone acetylation induced in the epigenome. We discuss the application of the method in the context of basic and translational sciences. PMID:25609656
A projection method to extract biological membrane models from 3D material models.
Roohbakhshan, Farshad; Duong, Thang X; Sauer, Roger A
2016-05-01
This paper presents a projection method for deriving membrane models from the corresponding three-dimensional material models. As a particular example the anisotropic Holzapfel-Gasser-Ogden model is considered. The projection procedure is based on the kinematical and constitutive assumptions of a general membrane theory, considering the membrane to be a general two-dimensional manifold. By assuming zero transverse stress, the Lagrange multiplier associated with the incompressibility constraint can be eliminated from the formulation. The resulting nonlinear model is discretized and linearized within the finite element method. Several numerical examples are shown, considering quadratic Lagrange and NURBS finite elements. These show that the proposed model is in very good agreement with analytical solutions and with full 3D finite element computations. PMID:26455810
A 3D Frictional Segment-to-Segment Contact Method for Large Deformations and Quadratic Elements
Puso, M; Laursen, T; Solberg, J
2004-04-01
Node-on-segment contact is the most common form of contact used today but has many deficiencies ranging from potential locking to non-smooth behavior with large sliding. Furthermore, node-on-segment approaches are not at all applicable to higher order discretizations (e.g. quadratic elements). In a previous work, [3, 4] we developed a segment-to-segment contact approach for eight node hexahedral elements based on the mortar method that was applicable to large deformation mechanics. The approach proved extremely robust since it eliminated the over-constraint that caused 'locking' and provided smooth force variations in large sliding. Here, we extend this previous approach to treat frictional contact problems. In addition, the method is extended to 3D quadratic tetrahedrals and hexahedrals. The proposed approach is then applied to several challenging frictional contact problems that demonstrate its effectiveness.
3D Pharmacophore, hierarchical methods, and 5-HT4 receptor binding data.
Varin, Thibault; Saettel, Nicolas; Villain, Jonathan; Lesnard, Aurelien; Dauphin, François; Bureau, Ronan; Rault, Sylvain
2008-10-01
5-Hydroxytryptamine subtype-4 (5-HT(4)) receptors have stimulated considerable interest amongst scientists and clinicians owing to their importance in neurophysiology and potential as therapeutic targets. A comparative analysis of hierarchical methods applied to data from one thousand 5-HT(4) receptor-ligand binding interactions was carried out. The chemical structures were described as chemical and pharmacophore fingerprints. The definitions of indices, related to the quality of the hierarchies in being able to distinguish between active and inactive compounds, revealed two interesting hierarchies with the Unity (1 active cluster) and pharmacophore fingerprints (4 active clusters). The results of this study also showed the importance of correct choice of metrics as well as the effectiveness of a new alternative of the Ward clustering algorithm named Energy (Minimum E-Distance method). In parallel, the relationship between these classifications and a previously defined 3D 5-HT(4) antagonist pharmacophore was established. PMID:18821249