Hex-dominant mesh generation using 3D constrained triangulation
OWEN,STEVEN J.
2000-05-30
A method for decomposing a volume with a prescribed quadrilateral surface mesh, into a hexahedral-dominated mesh is proposed. With this method, known as Hex-Morphing (H-Morph), an initial tetrahedral mesh is provided. Tetrahedral are transformed and combined starting from the boundary and working towards the interior of the volume. The quadrilateral faces of the hexahedra are treated as internal surfaces, which can be recovered using constrained triangulation techniques. Implementation details of the edge and face recovery process are included. Examples and performance of the H-Morph algorithm are also presented.
3D Laser Triangulation for Plant Phenotyping in Challenging Environments
Kjaer, Katrine Heinsvig; Ottosen, Carl-Otto
2015-01-01
To increase the understanding of how the plant phenotype is formed by genotype and environmental interactions, simple and robust high-throughput plant phenotyping methods should be developed and considered. This would not only broaden the application range of phenotyping in the plant research community, but also increase the ability for researchers to study plants in their natural environments. By studying plants in their natural environment in high temporal resolution, more knowledge on how multiple stresses interact in defining the plant phenotype could lead to a better understanding of the interaction between plant responses and epigenetic regulation. In the present paper, we evaluate a commercial 3D NIR-laser scanner (PlantEye, Phenospex B.V., Herleen, The Netherlands) to track daily changes in plant growth with high precision in challenging environments. Firstly, we demonstrate that the NIR laser beam of the scanner does not affect plant photosynthetic performance. Secondly, we demonstrate that it is possible to estimate phenotypic variation amongst the growth pattern of ten genotypes of Brassica napus L. (rapeseed), using a simple linear correlation between scanned parameters and destructive growth measurements. Our results demonstrate the high potential of 3D laser triangulation for simple measurements of phenotypic variation in challenging environments and in a high temporal resolution. PMID:26066990
Triangulation Based 3D Laser Imaging for Fracture Orientation Analysis
NASA Astrophysics Data System (ADS)
Mah, J.; Claire, S.; Steve, M.
2009-05-01
Laser imaging has recently been identified as a potential tool for rock mass characterization. This contribution focuses on the application of triangulation based, short-range laser imaging to determine fracture orientation and surface texture. This technology measures the distance to the target by triangulating the projected and reflected laser beams, and also records the reflection intensity. In this study, we acquired 3D laser images of rock faces using the Laser Camera System (LCS), a portable instrument developed by Neptec Design Group (Ottawa, Canada). The LCS uses an infrared laser beam and is immune to the lighting conditions. The maximum image resolution is 1024 x 1024 volumetric image elements. Depth resolution is 0.5 mm at 5 m. An above ground field trial was conducted at a blocky road cut with well defined joint sets (Kingston, Ontario). An underground field trial was conducted at the Inco 175 Ore body (Sudbury, Ontario) where images were acquired in the dark and the joint set features were more subtle. At each site, from a distance of 3 m away from the rock face, a grid of six images (approximately 1.6 m by 1.6 m) was acquired at maximum resolution with 20% overlap between adjacent images. This corresponds to a density of 40 image elements per square centimeter. Polyworks, a high density 3D visualization software tool, was used to align and merge the images into a single digital triangular mesh. The conventional method of determining fracture orientations is by manual measurement using a compass. In order to be accepted as a substitute for this method, the LCS should be capable of performing at least to the capabilities of manual measurements. To compare fracture orientation estimates derived from the 3D laser images to manual measurements, 160 inclinometer readings were taken at the above ground site. Three prominent joint sets (strike/dip: 236/09, 321/89, 325/01) were identified by plotting the joint poles on a stereonet. Underground, two main joint
Onomatopoeia characters extraction from comic images using constrained Delaunay triangulation
NASA Astrophysics Data System (ADS)
Liu, Xiangping; Shoji, Kenji; Mori, Hiroshi; Toyama, Fubito
2014-02-01
A method for extracting onomatopoeia characters from comic images was developed based on stroke width feature of characters, since they nearly have a constant stroke width in a number of cases. An image was segmented with a constrained Delaunay triangulation. Connected component grouping was performed based on the triangles generated by the constrained Delaunay triangulation. Stroke width calculation of the connected components was conducted based on the altitude of the triangles generated with the constrained Delaunay triangulation. The experimental results proved the effectiveness of the proposed method.
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.
Computing 2D constrained delaunay triangulation using the GPU.
Qi, Meng; Cao, Thanh-Tung; Tan, Tiow-Seng
2013-05-01
We propose the first graphics processing unit (GPU) solution to compute the 2D constrained Delaunay triangulation (CDT) of a planar straight line graph (PSLG) consisting of points and edges. There are many existing CPU algorithms to solve the CDT problem in computational geometry, yet there has been no prior approach to solve this problem efficiently using the parallel computing power of the GPU. For the special case of the CDT problem where the PSLG consists of just points, which is simply the normal Delaunay triangulation (DT) problem, a hybrid approach using the GPU together with the CPU to partially speed up the computation has already been presented in the literature. Our work, on the other hand, accelerates the entire computation on the GPU. Our implementation using the CUDA programming model on NVIDIA GPUs is numerically robust, and runs up to an order of magnitude faster than the best sequential implementations on the CPU. This result is reflected in our experiment with both randomly generated PSLGs and real-world GIS data having millions of points and edges. PMID:23492377
Performance Analysis of a Low-Cost Triangulation-Based 3d Camera: Microsoft Kinect System
NASA Astrophysics Data System (ADS)
. K. Chow, J. C.; Ang, K. D.; Lichti, D. D.; Teskey, W. F.
2012-07-01
Recent technological advancements have made active imaging sensors popular for 3D modelling and motion tracking. The 3D coordinates of signalised targets are traditionally estimated by matching conjugate points in overlapping images. Current 3D cameras can acquire point clouds at video frame rates from a single exposure station. In the area of 3D cameras, Microsoft and PrimeSense have collaborated and developed an active 3D camera based on the triangulation principle, known as the Kinect system. This off-the-shelf system costs less than 150 USD and has drawn a lot of attention from the robotics, computer vision, and photogrammetry disciplines. In this paper, the prospect of using the Kinect system for precise engineering applications was evaluated. The geometric quality of the Kinect system as a function of the scene (i.e. variation of depth, ambient light conditions, incidence angle, and object reflectivity) and the sensor (i.e. warm-up time and distance averaging) were analysed quantitatively. This system's potential in human body measurements was tested against a laser scanner and 3D range camera. A new calibration model for simultaneously determining the exterior orientation parameters, interior orientation parameters, boresight angles, leverarm, and object space features parameters was developed and the effectiveness of this calibration approach was explored.
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
Spectral triangulation: a 3D method for locating single-walled carbon nanotubes in vivo
NASA Astrophysics Data System (ADS)
Lin, Ching-Wei; Bachilo, Sergei M.; Vu, Michael; Beckingham, Kathleen M.; Bruce Weisman, R.
2016-05-01
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.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
Bound constrained bundle adjustment for reliable 3D reconstruction.
Gong, Yuanzheng; Meng, De; Seibel, Eric J
2015-04-20
Bundle adjustment (BA) is a common estimation algorithm that is widely used in machine vision as the last step in a feature-based three-dimensional (3D) reconstruction algorithm. BA is essentially a non-convex non-linear least-square problem that can simultaneously solve the 3D coordinates of all the feature points describing the scene geometry, as well as the parameters of the camera. The conventional BA takes a parameter either as a fixed value or as an unconstrained variable based on whether the parameter is known or not. In cases where the known parameters are inaccurate but constrained in a range, conventional BA results in an incorrect 3D reconstruction by using these parameters as fixed values. On the other hand, these inaccurate parameters can be treated as unknown variables, but this does not exploit the knowledge of the constraints, and the resulting reconstruction can be erroneous since the BA optimization halts at a dramatically incorrect local minimum due to its non-convexity. In many practical 3D reconstruction applications, unknown variables with range constraints are usually available, such as a measurement with a range of uncertainty or a bounded estimate. Thus to better utilize these pre-known, constrained, but inaccurate parameters, a bound constrained bundle adjustment (BCBA) algorithm is proposed, developed and tested in this study. A scanning fiber endoscope (the camera) is used to capture a sequence of images above a surgery phantom (the object) of known geometry. 3D virtual models are reconstructed based on these images and then compared with the ground truth. The experimental results demonstrate BCBA can achieve a more reliable, rapid, and accurate 3D reconstruction than conventional bundle adjustment. PMID:25969115
Bound constrained bundle adjustment for reliable 3D reconstruction
Gong, Yuanzheng; Meng, De; Seibel, Eric J.
2015-01-01
Bundle adjustment (BA) is a common estimation algorithm that is widely used in machine vision as the last step in a feature-based three-dimensional (3D) reconstruction algorithm. BA is essentially a non-convex non-linear least-square problem that can simultaneously solve the 3D coordinates of all the feature points describing the scene geometry, as well as the parameters of the camera. The conventional BA takes a parameter either as a fixed value or as an unconstrained variable based on whether the parameter is known or not. In cases where the known parameters are inaccurate but constrained in a range, conventional BA results in an incorrect 3D reconstruction by using these parameters as fixed values. On the other hand, these inaccurate parameters can be treated as unknown variables, but this does not exploit the knowledge of the constraints, and the resulting reconstruction can be erroneous since the BA optimization halts at a dramatically incorrect local minimum due to its non-convexity. In many practical 3D reconstruction applications, unknown variables with range constraints are usually available, such as a measurement with a range of uncertainty or a bounded estimate. Thus to better utilize these pre-known, constrained, but inaccurate parameters, a bound constrained bundle adjustment (BCBA) algorithm is proposed, developed and tested in this study. A scanning fiber endoscope (the camera) is used to capture a sequence of images above a surgery phantom (the object) of known geometry. 3D virtual models are reconstructed based on these images and then compared with the ground truth. The experimental results demonstrate BCBA can achieve a more reliable, rapid, and accurate 3D reconstruction than conventional bundle adjustment. PMID:25969115
Improved EEG source localization employing 3D sensing by "Flying Triangulation"
NASA Astrophysics Data System (ADS)
Ettl, Svenja; Rampp, Stefan; Fouladi-Movahed, Sarah; Dalal, Sarang S.; Willomitzer, Florian; Arold, Oliver; Stefan, Hermann; Häusler, Gerd
2013-04-01
With electroencephalography (EEG), a person's brain activity can be monitored over time and sources of activity localized. With this information, brain regions showing pathological activity, such as epileptic spikes, can be delineated. In cases of severe drug-resistant epilepsy, surgical resection of these brain regions may be the only treatment option. This requires a precise localization of the responsible seizure generators. They can be reconstructed from EEG data when the electrode positions are known. The standard method employs a "digitization pen" and has severe drawbacks: It is time consuming, the result is user-dependent, and the patient has to hold still. We present a novel method which overcomes these drawbacks. It is based on the optical "Flying Triangulation" (FlyTri) sensor which allows a motion-robust acquisition of precise 3D data. To compare the two methods, the electrode positions were determined with each method for a real-sized head model with EEG electrodes and their deviation to the ground-truth data calculated. The standard deviation for the current method was 3.39 mm while it was 0.98 mm for the new method. The influence of these results on the final EEG source localization was investigated by simulating EEG data. The digitization pen result deviates substantially from the true source location and time series. In contrast, the FlyTri result agrees with the original information. Our findings suggest that FlyTri might become a valuable tool in the field of medical brain research, because of its improved precision and contactless handling. Future applications might include co-registration of multimodal information.
NASA Astrophysics Data System (ADS)
Willomitzer, Florian; Ettl, Svenja; Arold, Oliver; Häusler, Gerd
2013-05-01
The three-dimensional shape acquisition of objects has become more and more important in the last years. Up to now, there are several well-established methods which already yield impressive results. However, even under quite common conditions like object movement or a complex shaping, most methods become unsatisfying. Thus, the 3D shape acquisition is still a difficult and non-trivial task. We present our measurement principle "Flying Triangulation" which enables a motion-robust 3D acquisition of complex-shaped object surfaces by a freely movable handheld sensor. Since "Flying Triangulation" is scalable, a whole sensor-zoo for different object sizes is presented. Concluding, an overview of current and future fields of investigation is given.
NASA Astrophysics Data System (ADS)
Vauhkonen, J.
2015-03-01
Reconstruction of three-dimensional (3D) forest canopy is described and quantified using airborne laser scanning (ALS) data with densities of 0.6-0.8 points m-2 and field measurements aggregated at resolutions of 400-900 m2. The reconstruction was based on computational geometry, topological connectivity, and numerical optimization. More precisely, triangulations and their filtrations, i.e. ordered sets of simplices belonging to the triangulations, based on the point data were analyzed. Triangulating the ALS point data corresponds to subdividing the underlying space of the points into weighted simplicial complexes with weights quantifying the (empty) space delimited by the points. Reconstructing the canopy volume populated by biomass will thus likely require filtering to exclude that volume from canopy voids. The approaches applied for this purpose were (i) to optimize the degree of filtration with respect to the field measurements, and (ii) to predict this degree by means of analyzing the persistent homology of the obtained triangulations, which is applied for the first time for vegetation point clouds. When derived from optimized filtrations, the total tetrahedral volume had a high degree of determination (R2) with the stem volume considered, both alone (R2=0.65) and together with other predictors (R2=0.78). When derived by analyzing the topological persistence of the point data and without any field input, the R2 were lower, but the predictions still showed a correlation with the field-measured stem volumes. Finally, producing realistic visualizations of a forested landscape using the persistent homology approach is demonstrated.
NASA Astrophysics Data System (ADS)
Nielsen, O. F.; Ploug, C.; Mendoza, J. A.; Martínez, K.
2009-05-01
The need for increaseding accuracy and reduced ambiguities in the inversion results has resulted in focus on the development of more advanced inversion methods of geophysical data. Over the past few years more advanced inversion techniques have been developed to improve the results. Real 3D-inversion is time consuming and therefore often not the best solution in a cost-efficient perspective. This has motivated the development of 3D constrained inversions, where 1D-models are constrained in 3D, also known as a Spatial Constrained Inversion (SCI). Moreover, inversion of several different data types in one inversion has been developed, known as Mutually Constrained Inversion (MCI). In this paper a presentation of a Spatial Mutually Constrained Inversion method (SMCI) is given. This method allows 1D-inversion applied to different geophysical datasets and geological information constrained in 3D. Application of two or more types of geophysical methods in the inversion has proved to reduce the equivalence problem and to increase the resolution in the inversion results. The use of geological information from borehole data or digital geological models can be integrated in the inversion. In the SMCI, a 1D inversion code is used to model soundings that are constrained in three dimensions according to their relative position in space. This solution enhances the accuracy of the inversion and produces distinct layers thicknesses and resistivities. It is very efficient in the mapping of a layered geology but still also capable of mapping layer discontinuities that are, in many cases, related to fracturing and faulting or due to valley fills. Geological information may be included in the inversion directly or used only to form a starting model for the individual soundings in the inversion. In order to show the effectiveness of the method, examples are presented from both synthetic data and real data. The examples include DC-soundings as well as land-based and airborne TEM
NASA Astrophysics Data System (ADS)
Tajbakhsh, Touraj
2010-02-01
A basic concern of computer graphic is the modeling and realistic representation of three-dimensional objects. In this paper we present our reconstruction framework which determines a polygonal surface from a set of dense points such those typically obtained from laser scanners. We deploy the concept of adaptive blobs to achieve a first volumetric representation of the object. In the next step we estimate a coarse surface using the marching cubes method. We propose to deploy a depth-first search segmentation algorithm traversing a graph representation of the obtained polygonal mesh in order to identify all connected components. A so called supervised triangulation maps the coarse surfaces onto the dense point cloud. We optimize the mesh topology using edge exchange operations. For photo-realistic visualization of objects we finally synthesize optimal low-loss textures from available scene captures of different projections. We evaluate our framework on artificial data as well as real sensed data.
Brosed, Francisco Javier; Aguilar, Juan José; Guillomía, David; Santolaria, Jorge
2011-01-01
This article discusses different non contact 3D measuring strategies and presents a model for measuring complex geometry parts, manipulated through a robot arm, using a novel vision system consisting of a laser triangulation sensor and a motorized linear stage. First, the geometric model incorporating an automatic simple module for long term stability improvement will be outlined in the article. The new method used in the automatic module allows the sensor set up, including the motorized linear stage, for the scanning avoiding external measurement devices. In the measurement model the robot is just a positioning of parts with high repeatability. Its position and orientation data are not used for the measurement and therefore it is not directly "coupled" as an active component in the model. The function of the robot is to present the various surfaces of the workpiece along the measurement range of the vision system, which is responsible for the measurement. Thus, the whole system is not affected by the robot own errors following a trajectory, except those due to the lack of static repeatability. For the indirect link between the vision system and the robot, the original model developed needs only one first piece measuring as a "zero" or master piece, known by its accurate measurement using, for example, a Coordinate Measurement Machine. The strategy proposed presents a different approach to traditional laser triangulation systems on board the robot in order to improve the measurement accuracy, and several important cues for self-recalibration are explored using only a master piece. Experimental results are also presented to demonstrate the technique and the final 3D measurement accuracy. PMID:22346569
Brosed, Francisco Javier; Aguilar, Juan José; Guillomía, David; Santolaria, Jorge
2011-01-01
This article discusses different non contact 3D measuring strategies and presents a model for measuring complex geometry parts, manipulated through a robot arm, using a novel vision system consisting of a laser triangulation sensor and a motorized linear stage. First, the geometric model incorporating an automatic simple module for long term stability improvement will be outlined in the article. The new method used in the automatic module allows the sensor set up, including the motorized linear stage, for the scanning avoiding external measurement devices. In the measurement model the robot is just a positioning of parts with high repeatability. Its position and orientation data are not used for the measurement and therefore it is not directly “coupled” as an active component in the model. The function of the robot is to present the various surfaces of the workpiece along the measurement range of the vision system, which is responsible for the measurement. Thus, the whole system is not affected by the robot own errors following a trajectory, except those due to the lack of static repeatability. For the indirect link between the vision system and the robot, the original model developed needs only one first piece measuring as a “zero” or master piece, known by its accurate measurement using, for example, a Coordinate Measurement Machine. The strategy proposed presents a different approach to traditional laser triangulation systems on board the robot in order to improve the measurement accuracy, and several important cues for self-recalibration are explored using only a master piece. Experimental results are also presented to demonstrate the technique and the final 3D measurement accuracy. PMID:22346569
NASA Astrophysics Data System (ADS)
Schneeberger, Raphael; de la Varga, Miguel; Florian Wellmann, J.; Kober, Florian; Berger, Alfons; Herwegh, Marco
2016-04-01
Fluid circulation in crystalline rocks is of key importance when exploring crystalline basement in light of, for example, deep-seated geothermal energy projects or selection of sites for nuclear waste repositories. Due to their enhanced permeability, fluid circulation within crystalline bedrock is mainly controlled by fault zones, which may originate from ductile mylonites but show a strong brittle overprint. In order to better constrain 3D flow paths, a well-founded knowledge on the 3D nature of the fault zone pattern is indispensable. We attempt to constrain the geometry of a complex 3D fault zone pattern in a case study of the Grimsel Test Site (GTS, central Switzerland). The constraints are based on mapping of both the surface as well as the GTS underground tunnel system, offering a unique opportunity to test the 3D model and associated uncertainties. We investigate the effect of increasing geoinformation on the quality and accuracy of the 3D model by using: (i) remote sensing surface data only, (ii) field surface mapping in combination with (i), and (iii) underground data combined with (i) and (ii). This approach allows for defining different steps in 3D geological modelling of a specific area, including a measure of the remaining uncertainty after each step. We obtain a best-estimate model by fitting results between surface and underground data by using a combination of field data and orientation obtained by Delaunay triangulation. We incorporate novel approaches to uncertainty analysis of fault orientations and investigate different fault planes showing the possible variation range of the structures investigated.
NASA Astrophysics Data System (ADS)
Liu, Shuang; Hu, Xiangyun; Xi, Yufei; Liu, Tianyou
2015-03-01
The regular grid discretization is prevalent in the inverse modeling for gravity and magnetic data. However, this subdivision strategy performs lower precision to represent the rugged observation surface. To deal with this problem, we evaluate a non-structured discretization method in which the subsurface with rolling terrain is divided into numbers of Delaunay triangular cells and each mesh has the uniform physical property distributions. The gravity and magnetic anomalies of a complex-shaped anomalous body are represented as the summaries of the single anomaly produced by each triangle field source. When inverting for the potential field data, we specify a minimization objective function composed of data constraints and then use the preconditioned conjugate gradient algorithm to iteratively solve the matrix minimization equations, where the preconditioner is determined by the distances between triangular cells and surface observers. We test our method using synthetic data; all tests return favorable results. In the case studies involving the gravity and magnetic anomalies of the Mengku and Pobei deposits in Xinjiang, northwest China, the inferred magnetite orebodies and ultrabasic rocks distributions are verified by the additional drilling and geological information. The discretization of constrained Delaunay triangulation provides an useful approach of computing and inverting the potential field data on the situations of undulate topography and complicated objects.
Necking and failure of constrained 3D microtissues induced by cellular tension
Wang, Hailong; Svoronos, Alexander A.; Boudou, Thomas; Sakar, Mahmut Selman; Schell, Jacquelyn Youssef; Morgan, Jeffrey R.; Chen, Christopher S.; Shenoy, Vivek B.
2013-01-01
In this paper we report a fundamental morphological instability of constrained 3D microtissues induced by positive chemomechanical feedback between actomyosin-driven contraction and the mechanical stresses arising from the constraints. Using a 3D model for mechanotransduction we find that perturbations in the shape of contractile tissues grow in an unstable manner leading to formation of “necks” that lead to the failure of the tissue by narrowing and subsequent elongation. The magnitude of the instability is shown to be determined by the level of active contractile strain, the stiffness of the extracellular matrix, and the components of the tissue that act in parallel with the active component and the stiffness of the boundaries that constrain the tissue. A phase diagram that demarcates stable and unstable behavior of 3D tissues as a function of these material parameters is derived. The predictions of our model are verified by analyzing the necking and failure of normal human fibroblast tissue constrained in a loop-ended dog-bone geometry and cardiac microtissues constrained between microcantilevers. By analyzing the time evolution of the morphology of the constrained tissues we have quantitatively determined the chemomechanical coupling parameters that characterize the generation of active stresses in these tissues. More generally, the analytical and numerical methods we have developed provide a quantitative framework to study how contractility can influence tissue morphology in complex 3D environments such as morphogenesis and organogenesis. PMID:24324149
Rate-constrained 3D surface estimation from noise-corrupted multiview depth videos.
Sun, Wenxiu; Cheung, Gene; Chou, Philip A; Florencio, Dinei; Zhang, Cha; Au, Oscar C
2014-07-01
Transmitting compactly represented geometry of a dynamic 3D scene from a sender can enable a multitude of imaging functionalities at a receiver, such as synthesis of virtual images at freely chosen viewpoints via depth-image-based rendering. While depth maps—projections of 3D geometry onto 2D image planes at chosen camera viewpoints-can nowadays be readily captured by inexpensive depth sensors, they are often corrupted by non-negligible acquisition noise. Given depth maps need to be denoised and compressed at the encoder for efficient network transmission to the decoder, in this paper, we consider the denoising and compression problems jointly, arguing that doing so will result in a better overall performance than the alternative of solving the two problems separately in two stages. Specifically, we formulate a rate-constrained estimation problem, where given a set of observed noise-corrupted depth maps, the most probable (maximum a posteriori (MAP)) 3D surface is sought within a search space of surfaces with representation size no larger than a prespecified rate constraint. Our rate-constrained MAP solution reduces to the conventional unconstrained MAP 3D surface reconstruction solution if the rate constraint is loose. To solve our posed rate-constrained estimation problem, we propose an iterative algorithm, where in each iteration the structure (object boundaries) and the texture (surfaces within the object boundaries) of the depth maps are optimized alternately. Using the MVC codec for compression of multiview depth video and MPEG free viewpoint video sequences as input, experimental results show that rate-constrained estimated 3D surfaces computed by our algorithm can reduce coding rate of depth maps by up to 32% compared with unconstrained estimated surfaces for the same quality of synthesized virtual views at the decoder. PMID:24876124
3D Geological Model of Nihe ore deposit Constrained by Gravity and Magnetic Modeling
NASA Astrophysics Data System (ADS)
Qi, Guang; Yan, Jiayong; Lv, Qingtan; Zhao, Jinhua
2016-04-01
observed data, and then adjust the model until a satisfactory accuracy of errors is achieved. It is hope that this work can provide reference for similar work in other areas. this study shows that the research of geologic constrained 3D gravity and magnetic modeling has potential value in the aspects of deep mineral exploration and mineral reserves estimation.
Meyer-Hermann, Michael
2008-01-01
Mathematical methods in Biology are of increasing relevance for understanding the control and the dynamics of biological systems with medical relevance. In particular, agent-based methods turn more and more important because of fast increasing computational power which makes even large systems accessible. An overview of different mathematical methods used in Theoretical Biology is provided and a novel agent-based method for cell mechanics based on Delaunay-triangulations and Voronoi-tessellations is explained in more detail: The Delaunay-Object-Dynamics method. It is claimed that the model combines physically realistic cell mechanics with a reasonable computational load. The power of the approach is illustrated with two examples, avascular tumor growth and genesis of lymphoid tissue in a cell-flow equilibrium. PMID:18023735
Constraining 3D Process Sedimentological Models to Geophysical Data Using Image Quilting
NASA Astrophysics Data System (ADS)
Tahmasebi, P.; Da Pra, A.; Pontiggia, M.; Caers, J.
2014-12-01
3D process geological models, whether for carbonate or sedimentological systems, have been proposed for modeling realistic subsurface heterogeneity. The problem with such forward process models is that they are not constrained to any subsurface data whether to wells or geophysical surveys. We propose a new method for realistic geological modeling of complex heterogeneity by hybridizing 3D process modeling of geological deposition with conditioning by means of a novel multiple-point geostatistics (MPS) technique termed image quilting (IQ). Image quilting is a pattern-based techniques that stiches together patterns extracted from training images to generate stochastic realizations that look like the training image. In this paper, we illustrate how 3D process model realizations can be used as training images in image quilting. To constrain the realization to seismic data we first interpret each facies in the geophysical data. These interpretation, while overly smooth and not reflecting finer scale variation are used as auxiliary variables in the generation of the image quilting realizations. To condition to well data, we first perform a kriging of the well data to generate a kriging map and kriging variance. The kriging map is used as additional auxiliary variable while the kriging variance is used as a weight given to the kriging derived auxiliary variable. We present an application to a giant offshore reservoir. Starting from seismic advanced attribute analysis and sedimentological interpretation, we build the 3D sedimentological process based model and use it as non-stationary training image for conditional image quilting.
3D automatic liver segmentation using feature-constrained Mahalanobis distance in CT images.
Salman Al-Shaikhli, Saif Dawood; Yang, Michael Ying; Rosenhahn, Bodo
2016-08-01
Automatic 3D liver segmentation is a fundamental step in the liver disease diagnosis and surgery planning. This paper presents a novel fully automatic algorithm for 3D liver segmentation in clinical 3D computed tomography (CT) images. Based on image features, we propose a new Mahalanobis distance cost function using an active shape model (ASM). We call our method MD-ASM. Unlike the standard active shape model (ST-ASM), the proposed method introduces a new feature-constrained Mahalanobis distance cost function to measure the distance between the generated shape during the iterative step and the mean shape model. The proposed Mahalanobis distance function is learned from a public database of liver segmentation challenge (MICCAI-SLiver07). As a refinement step, we propose the use of a 3D graph-cut segmentation. Foreground and background labels are automatically selected using texture features of the learned Mahalanobis distance. Quantitatively, the proposed method is evaluated using two clinical 3D CT scan databases (MICCAI-SLiver07 and MIDAS). The evaluation of the MICCAI-SLiver07 database is obtained by the challenge organizers using five different metric scores. The experimental results demonstrate the availability of the proposed method by achieving an accurate liver segmentation compared to the state-of-the-art methods. PMID:26501155
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
NASA Astrophysics Data System (ADS)
Tchikaya, Euloge Budet; Chouteau, Michel; Keating, Pierre; Shamsipour, Pejman
2016-02-01
We present an inversion tool for airborne gravity gradient data that yields a 3D density model using stochastic methods i.e. cokriging and conditional simulation. This method uses geostatistical properties of the measured gravity gradient to estimate a 3D density model whose gravity response fits the measured gravity gradient anomaly. Linearity between gravity gradient data and density allows estimation of the model (density) covariance using observed data, i.e. we adjust iteratively the density covariance matrix by fitting experimental and theoretical gravity gradient covariance matrices. Inversion can be constrained by including densities known at some locations. In addition we can explore various reasonable solutions that honour both the estimated density covariance model and the gravity gradient data using geostatistical simulation. The proposed method is first tested with two synthetic datasets generated from a sharp-boundary model and a smooth stochastic model respectively. The results show the method to be capable of retrieving models compatible with the true models; it also allows the integration of complex a priori information. The technique is then applied to gravity gradient survey data collected for the Geological Survey of Canada in the area of McFaulds Lake (Ontario, Canada) using the Falcon airborne gravity system. Unconstrained inversion returns a density model that is geologically plausible and the computed response exactly fits the observed gravity gradient anomaly.
Statistical 3D Prostate Imaging Atlas Construction via Anatomically Constrained Registration
Rusu, Mirabela; Bloch, B. Nicolas; Jaffe, Carl C.; Rofsky, Neil M.; Genega, Elizabeth M.; Feleppa, Ernest; Lenkinski, Robert E.; Madabhushi, Anant
2013-01-01
Statistical imaging atlases allow for integration of information from multiple patient studies collected across different image scales and modalities, such as multi-parametric (MP) MRI and histology, providing population statistics regarding a specific pathology within a single canonical representation. Such atlases are particularly valuable in the identification and validation of meaningful imaging signatures for disease characterization in vivo within a population. Despite the high incidence of prostate cancer, an imaging atlas focused on different anatomic structures of the prostate, i.e. an anatomic atlas, has yet to be constructed. In this work we introduce a novel framework for MRI atlas construction that uses an iterative, anatomically constrained registration (AnCoR) scheme to enable the proper alignment of the prostate (Pr) and central gland (CG) boundaries. Our current implementation uses endorectal, 1.5T or 3T, T2-weighted MRI from 51 patients with biopsy confirmed cancer; however, the prostate atlas is seamlessly extensible to include additional MRI parameters. In our cohort, radical prostatectomy is performed following MP-MR image acquisition; thus ground truth annotations for prostate cancer are available from the histological specimens. Once mapped onto MP-MRI through elastic registration of histological slices to corresponding T2-w MRI slices, the annotations are utilized by the AnCoR framework to characterize the 3D statistical distribution of cancer per anatomic structure. Such distributions are useful for guiding biopsies toward regions of higher cancer likelihood and understanding imaging profiles for disease extent in vivo. We evaluate our approach via the Dice similarity coefficient (DSC) for different anatomic structures (delineated by expert radiologists): Pr, CG and peripheral zone (PZ). The AnCoR-based atlas had a CG DSC of 90.36%, and Pr DSC of 89.37%. Moreover, we evaluated the deviation of anatomic landmarks, the urethra and
A 3-D density model of Greece constrained by gravity and seismic data
NASA Astrophysics Data System (ADS)
Makris, Jannis; Papoulia, Joanna; Yegorova, Tamara
2013-07-01
A 3-D density model of Greece was developed by gravity modelling constrained by 2-D seismic profiles. Densities were defined from seismic velocities using the Nafe & Drake and Birch empirical functions for the sediments, crust and upper mantle. Sediments in the North Aegean are 6 km thick, and are deposited in transtensional basins developing by dextral strike slip motion of the North Anatolian Fault. The Cyclades, central Aegean Sea, are free of sediments. South of Crete, in the Libyan Sea, sediments are approximately 11 km thick. At the western Hellenides sediments of up to 8 km thickness have been accumulated in basins formed by crustal bending and southwestwards thrusting of the Hellenic napes. At a deeper crustal level variations of crustal type and thickness cause density variations explaining large part of the observed gravity field. The North Aegean domain is characterized by a 24-km-thick continental crust, including sediments, whereas the western Cyclades, in central Aegean area, have a slightly thickened crust of 26 km. Crustal thicknesses vary between 16 km in the deep Ionian and Cretan Seas to 40 km in the western Hellenides. In western Crete crust is 30-32 km thick, thinning eastwards to only 26 km. The deep Ionian basin, the Mediterranean Ridge, as well as most of the Libyan Sea are underlain by oceanic crust. In western Turkey the crust thickens from 30 km along the coast to 34 km to the interior. A third deeper level of density variations occurs in the upper mantle. Subduction of the oceanic lithosphere below the Aegean continental domain destabilizes the thermal field, uplifting the isotherms by convection and conduction below the Aegean Sea. Consequently, volume expansion of the upper mantle and lithological changes reduce its density and depress the gravity intensity. This low density-velocity upper mantle extends from the Sporades islands in the North Aegean to the Cretan Sea, occupying the space between the cold subducted Ionian oceanic
Statistical 3D prostate imaging atlas construction via anatomically constrained registration
NASA Astrophysics Data System (ADS)
Rusu, Mirabela; Bloch, B. Nicolas; Jaffe, Carl C.; Rofsky, Neil M.; Genega, Elizabeth M.; Feleppa, Ernest; Lenkinski, Robert E.; Madabhushi, Anant
2013-03-01
Statistical imaging atlases allow for integration of information from multiple patient studies collected across different image scales and modalities, such as multi-parametric (MP) MRI and histology, providing population statistics regarding a specific pathology within a single canonical representation. Such atlases are particularly valuable in the identification and validation of meaningful imaging signatures for disease characterization in vivo within a population. Despite the high incidence of prostate cancer, an imaging atlas focused on different anatomic structures of the prostate, i.e. an anatomic atlas, has yet to be constructed. In this work we introduce a novel framework for MRI atlas construction that uses an iterative, anatomically constrained registration (AnCoR) scheme to enable the proper alignment of the prostate (Pr) and central gland (CG) boundaries. Our current implementation uses endorectal, 1.5T or 3T, T2-weighted MRI from 51 patients with biopsy confirmed cancer; however, the prostate atlas is seamlessly extensible to include additional MRI parameters. In our cohort, radical prostatectomy is performed following MP-MR image acquisition; thus ground truth annotations for prostate cancer are available from the histological specimens. Once mapped onto MP-MRI through elastic registration of histological slices to corresponding T2-w MRI slices, the annotations are utilized by the AnCoR framework to characterize the 3D statistical distribution of cancer per anatomic structure. Such distributions are useful for guiding biopsies toward regions of higher cancer likelihood and understanding imaging profiles for disease extent in vivo. We evaluate our approach via the Dice similarity coefficient (DSC) for different anatomic structures (delineated by expert radiologists): Pr, CG and peripheral zone (PZ). The AnCoR-based atlas had a CG DSC of 90.36%, and Pr DSC of 89.37%. Moreover, we evaluated the deviation of anatomic landmarks, the urethra and
NASA Astrophysics Data System (ADS)
Przybycin, Anna M.; Scheck-Wenderoth, Magdalena; Schneider, Michael
2014-05-01
The North Alpine Foreland Basin is situated in the northern front of the European Alps and extends over parts of France, Switzerland, Germany and Austria. It formed as a wedge shaped depression since the Tertiary in consequence of the Euro - Adriatic continental collision and the Alpine orogeny. The basin is filled with clastic sediments, the Molasse, originating from erosional processes of the Alps and underlain by Mesozoic sedimentary successions and a Paleozoic crystalline crust. For our study we have focused on the German part of the basin. To investigate the deep structure, the isostatic state and the load distribution of this region we have constructed a 3D structural model of the basin and the Alpine area using available depth and thickness maps, regional scale 3D structural models as well as seismic and well data for the sedimentary part. The crust (from the top Paleozoic down to the Moho (Grad et al. 2008)) has been considered as two-parted with a lighter upper crust and a denser lower crust; the partition has been calculated following the approach of isostatic equilibrium of Pratt (1855). By implementing a seismic Lithosphere-Asthenosphere-Boundary (LAB) (Tesauro 2009) the crustal scale model has been extended to the lithospheric-scale. The layer geometry and the assigned bulk densities of this starting model have been constrained by means of 3D gravity modelling (BGI, 2012). Afterwards the 3D load distribution has been calculated using a 3D finite element method. Our results show that the North Alpine Foreland Basin is not isostatically balanced and that the configuration of the crystalline crust strongly controls the gravity field in this area. Furthermore, our results show that the basin area is influenced by varying lateral load differences down to a depth of more than 150 km what allows a first order statement of the required compensating horizontal stress needed to prevent gravitational collapse of the system. BGI (2012). The International
Automatic 3D kidney segmentation based on shape constrained GC-OAAM
NASA Astrophysics Data System (ADS)
Chen, Xinjian; Summers, Ronald M.; Yao, Jianhua
2011-03-01
The kidney can be classified into three main tissue types: renal cortex, renal medulla and renal pelvis (or collecting system). Dysfunction of different renal tissue types may cause different kidney diseases. Therefore, accurate and efficient segmentation of kidney into different tissue types plays a very important role in clinical research. In this paper, we propose an automatic 3D kidney segmentation method which segments the kidney into the three different tissue types: renal cortex, medulla and pelvis. The proposed method synergistically combines active appearance model (AAM), live wire (LW) and graph cut (GC) methods, GC-OAAM for short. Our method consists of two main steps. First, a pseudo 3D segmentation method is employed for kidney initialization in which the segmentation is performed slice-by-slice via a multi-object oriented active appearance model (OAAM) method. An improved iterative model refinement algorithm is proposed for the AAM optimization, which synergistically combines the AAM and LW method. Multi-object strategy is applied to help the object initialization. The 3D model constraints are applied to the initialization result. Second, the object shape information generated from the initialization step is integrated into the GC cost computation. A multi-label GC method is used to segment the kidney into cortex, medulla and pelvis. The proposed method was tested on 19 clinical arterial phase CT data sets. The preliminary results showed the feasibility and efficiency of the proposed method.
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
PDE constrained optimization of electrical defibrillation in a 3D ventricular slice geometry.
Chamakuri, Nagaiah; Kunisch, Karl; Plank, Gernot
2016-04-01
A computational study of an optimal control approach for cardiac defibrillation in a 3D geometry is presented. The cardiac bioelectric activity at the tissue and bath volumes is modeled by the bidomain model equations. The model includes intramural fiber rotation, axially symmetric around the fiber direction, and anisotropic conductivity coefficients, which are extracted from a histological image. The dynamics of the ionic currents are based on the regularized Mitchell-Schaeffer model. The controls enter in the form of electrodes, which are placed at the boundary of the bath volume with the goal of dampening undesired arrhythmias. The numerical optimization is based on Newton techniques. We demonstrated the parallel architecture environment for the computation of potentials on multidomains and for the higher order optimization techniques. Copyright © 2015 John Wiley & Sons, Ltd. PMID:26249168
3D pulse EPR imaging from sparse-view projections via constrained, total variation minimization
NASA Astrophysics Data System (ADS)
Qiao, Zhiwei; Redler, Gage; Epel, Boris; Qian, Yuhua; Halpern, Howard
2015-09-01
Tumors and tumor portions with low oxygen concentrations (pO2) have been shown to be resistant to radiation therapy. As such, radiation therapy efficacy may be enhanced if delivered radiation dose is tailored based on the spatial distribution of pO2 within the tumor. A technique for accurate imaging of tumor oxygenation is critically important to guide radiation treatment that accounts for the effects of local pO2. Electron paramagnetic resonance imaging (EPRI) has been considered one of the leading methods for quantitatively imaging pO2 within tumors in vivo. However, current EPRI techniques require relatively long imaging times. Reducing the number of projection scan considerably reduce the imaging time. Conventional image reconstruction algorithms, such as filtered back projection (FBP), may produce severe artifacts in images reconstructed from sparse-view projections. This can lower the utility of these reconstructed images. In this work, an optimization based image reconstruction algorithm using constrained, total variation (TV) minimization, subject to data consistency, is developed and evaluated. The algorithm was evaluated using simulated phantom, physical phantom and pre-clinical EPRI data. The TV algorithm is compared with FBP using subjective and objective metrics. The results demonstrate the merits of the proposed reconstruction algorithm.
3D pulse EPR imaging from sparse-view projections via constrained, total variation minimization
Qiao, Zhiwei; Redler, Gage; Epel, Boris; Halpern, Howard
2016-01-01
Tumors and tumor portions with low oxygen concentrations (pO2) have been shown to be resistant to radiation therapy. As such, radiation therapy efficacy may be enhanced if delivered radiation dose is tailored based on the spatial distribution of pO2 within the tumor. A technique for accurate imaging of tumor oxygenation is critically important to guide radiation treatment that accounts for the effects of local pO2. Electron paramagnetic resonance imaging (EPRI) has been considered one of the leading methods for quantitatively imaging pO2 within tumors in vivo. However, current EPRI techniques require relatively long imaging times. Reducing the number of projections can considerably reduce the imaging time. Conventional image reconstruction algorithms, such as filtered back projection (FBP), may produce severe artifacts in images reconstructed from sparse-view projections. This can lower the utility of these reconstructed images. In this work, an optimization based image reconstruction algorithm using constrained, total variation (TV) minimization, subject to data consistency, is developed and evaluated. The algorithm was evaluated using simulated phantom, physical phantom and pre-clinical EPRI data. The TV algorithm is compared with FBP using subjective and objective metrics. The results demonstrate the merits of the proposed reconstruction algorithm. PMID:26225440
A 3-D constrained mixture model for mechanically mediated vascular growth and remodeling
Wan, William; Hansen, Laura
2010-01-01
In contrast to the widely applied approach to model soft tissue remodeling employing the concept of volumetric growth, microstructurally motivated models are capable of capturing many of the underlying mechanisms of growth and remodeling; i.e., the production, removal, and remodeling of individual constituents at different rates and to different extents. A 3-dimensional constrained mixture computational framework has been developed for vascular growth and remodeling, considering new, microstructurally motivated kinematics and constitutive equations and new stress and muscle activation mediated evolution equations. Our computational results for alterations in flow and pressure, using reasonable physiological values for rates of constituent growth and turnover, concur with findings in the literature. For example, for flow-induced remodeling, our simulations predict that, although the wall shear stress is restored completely, the circumferential stress is not restored employing realistic physiological rate parameters. Also, our simulations predict different levels of thickening on inner versus outer wall locations, as shown in numerous reports of pressure-induced remodeling. Whereas the simulations are meant to be illustrative, they serve to highlight the experimental data currently lacking to fully quantify mechanically mediated adaptations in the vasculature. PMID:20039091
Jung, Jaehoon; Yoon, Sanghyun; Ju, Sungha; Heo, Joon
2015-01-01
The growing interest and use of indoor mapping is driving a demand for improved data-acquisition facility, efficiency and productivity in the era of the Building Information Model (BIM). The conventional static laser scanning method suffers from some limitations on its operability in complex indoor environments, due to the presence of occlusions. Full scanning of indoor spaces without loss of information requires that surveyors change the scanner position many times, which incurs extra work for registration of each scanned point cloud. Alternatively, a kinematic 3D laser scanning system, proposed herein, uses line-feature-based Simultaneous Localization and Mapping (SLAM) technique for continuous mapping. Moreover, to reduce the uncertainty of line-feature extraction, we incorporated constrained adjustment based on an assumption made with respect to typical indoor environments: that the main structures are formed of parallel or orthogonal line features. The superiority of the proposed constrained adjustment is its reduction for uncertainties of the adjusted lines, leading to successful data association process. In the present study, kinematic scanning with and without constrained adjustment were comparatively evaluated in two test sites, and the results confirmed the effectiveness of the proposed system. The accuracy of the 3D mapping result was additionally evaluated by comparison with the reference points acquired by a total station: the Euclidean average distance error was 0.034 m for the seminar room and 0.043 m for the corridor, which satisfied the error tolerance for point cloud acquisition (0.051 m) according to the guidelines of the General Services Administration for BIM accuracy. PMID:26501292
Jung, Jaehoon; Yoon, Sanghyun; Ju, Sungha; Heo, Joon
2015-01-01
The growing interest and use of indoor mapping is driving a demand for improved data-acquisition facility, efficiency and productivity in the era of the Building Information Model (BIM). The conventional static laser scanning method suffers from some limitations on its operability in complex indoor environments, due to the presence of occlusions. Full scanning of indoor spaces without loss of information requires that surveyors change the scanner position many times, which incurs extra work for registration of each scanned point cloud. Alternatively, a kinematic 3D laser scanning system, proposed herein, uses line-feature-based Simultaneous Localization and Mapping (SLAM) technique for continuous mapping. Moreover, to reduce the uncertainty of line-feature extraction, we incorporated constrained adjustment based on an assumption made with respect to typical indoor environments: that the main structures are formed of parallel or orthogonal line features. The superiority of the proposed constrained adjustment is its reduction for uncertainties of the adjusted lines, leading to successful data association process. In the present study, kinematic scanning with and without constrained adjustment were comparatively evaluated in two test sites, and the results confirmed the effectiveness of the proposed system. The accuracy of the 3D mapping result was additionally evaluated by comparison with the reference points acquired by a total station: the Euclidean average distance error was 0.034 m for the seminar room and 0.043 m for the corridor, which satisfied the error tolerance for point cloud acquisition (0.051 m) according to the guidelines of the General Services Administration for BIM accuracy. PMID:26501292
NASA Technical Reports Server (NTRS)
Madura, T. I.; Gull, T. R.; Owocki, S. P.; Groh, J. H.; Okazaki, A. T.; Russell, C. M. P.
2011-01-01
We present a three-dimensional (3-D) dynamical model for the broad [Fe III] emission observed in Eta Carinae using the Hubble Space Telescope/Space Telescope Imaging Spectrograph (HST/STIS). This model is based on full 3-D Smoothed Particle Hydrodynamics (SPH) simulations of Eta Car's binary colliding winds. Radiative transfer codes are used to generate synthetic spectro-images of [Fe III] emission line structures at various observed orbital phases and STIS slit position angles (PAs). Through a parameter study that varies the orbital inclination i, the PA(theta) that the orbital plane projection of the line-of-sight makes with the apastron side of the semi-major axis, and the PA on the sky of the orbital axis, we are able, for the first time, to tightly constrain the absolute 3-D orientation of the binary orbit. To simultaneously reproduce the blue-shifted emission arcs observed at orbital phase 0.976, STIS slit PA = +38deg, and the temporal variations in emission seen at negative slit PAs, the binary needs to have an i approx. = 130deg to 145deg, Theta approx. = -15deg to +30deg, and an orbital axis projected on the sky at a P A approx. = 302deg to 327deg east of north. This represents a system with an orbital axis that is closely aligned with the inferred polar axis of the Homunculus nebula, in 3-D. The companion star, Eta(sub B), thus orbits clockwise on the sky and is on the observer's side of the system at apastron. This orientation has important implications for theories for the formation of the Homunculus and helps lay the groundwork for orbital modeling to determine the stellar masses.
NASA Astrophysics Data System (ADS)
Madura, T. I.; Gull, T. R.; Owocki, S. P.; Groh, J. H.; Okazaki, A. T.; Russell, C. M. P.
2012-03-01
We present a three-dimensional (3D) dynamical model for the broad [Fe III] emission observed in η Carinae using the Hubble Space Telescope/Space Telescope Imaging Spectrograph (STIS). This model is based on full 3D smoothed particle hydrodynamics simulations of η Car's binary colliding winds. Radiative transfer codes are used to generate synthetic spectroimages of [Fe III] emission-line structures at various observed orbital phases and STIS slit position angles (PAs). Through a parameter study that varies the orbital inclination i, the PA θ that the orbital plane projection of the line of sight makes with the apastron side of the semimajor axis and the PA on the sky of the orbital axis, we are able, for the first time, to tightly constrain the absolute 3D orientation of the binary orbit. To simultaneously reproduce the blueshifted emission arcs observed at orbital phase 0.976, STIS slit PA =+38° and the temporal variations in emission seen at negative slit PAs, the binary needs to have an i≈ 130° to 145°, θ≈-15° to +30° and an orbital axis projected on the sky at a PA ≈ 302° to 327° east of north. This represents a system with an orbital axis that is closely aligned with the inferred polar axis of the Homunculus nebula, in 3D. The companion star, ηB, thus orbits clockwise on the sky and is on the observer's side of the system at apastron. This orientation has important implications for theories for the formation of the Homunculus and helps lay the groundwork for orbital modelling to determine the stellar masses.