Comparison of methods for curvature estimation from volume fractions

This paper evaluates and compares the accuracy and robustness of curvature estimation methods for three-dimensional interfaces represented implicitly by discrete volume fractions on a Cartesian mesh. The height function (HF) method is compared to three paraboloid fitting methods: fitting to the piecewise linear interface reconstruction centroids (PC), fitting to the piecewise linear interface reconstruction volumetrically (PV), and volumetrically fitting (VF) the paraboloid directly to the volume fraction field. The numerical studies presented in this work find that while the curvature error from the VF method converges with second-order accuracy as with the HF method for static interfaces represented by exact volume fractions, the PV method best balances low curvature errors with low computational cost for dynamic interfaces when the interface reconstruction and advection are coupled to a two-phase Navier-Stokes solver

Direct estimation of kinetic parametric images for dynamic PET.

Dynamic positron emission tomography (PET) can monitor spatiotemporal distribution of radiotracer in vivo. The spatiotemporal information can be used to estimate parametric images of radiotracer kinetics that are of physiological and biochemical interests. Direct estimation of parametric images from raw projection data allows accurate noise modeling and has been shown to offer better image quality than conventional indirect methods, which reconstruct a sequence of PET images first and then perform tracer kinetic modeling pixel-by-pixel. Direct reconstruction of parametric images has gained increasing interests with the advances in computing hardware. Many direct reconstruction algorithms have been developed for different kinetic models. In this paper we review the recent progress in the development of direct reconstruction algorithms for parametric image estimation. Algorithms for linear and nonlinear kinetic models are described and their properties are discussed

Surface generating on the basis of experimental data

The problem of surface generating on the basis of experimental data is presented in this lecture. Special attention is given to the implementation of moving ordinary least squares and moving total least squares. Some results done in the Institute for Applied Mathematics in Osijek are mentioned which were published in the last several years

A Minimalist Approach to Type-Agnostic Detection of Quadrics in Point Clouds

This paper proposes a segmentation-free, automatic and efficient procedure to detect general geometric quadric forms in point clouds, where clutter and occlusions are inevitable. Our everyday world is dominated by man-made objects which are designed using 3D primitives (such as planes, cones, spheres, cylinders, etc.). These objects are also omnipresent in industrial environments. This gives rise to the possibility of abstracting 3D scenes through primitives, thereby positions these geometric forms as an integral part of perception and high level 3D scene understanding. As opposed to state-of-the-art, where a tailored algorithm treats each primitive type separately, we propose to encapsulate all types in a single robust detection procedure. At the center of our approach lies a closed form 3D quadric fit, operating in both primal & dual spaces and requiring as low as 4 oriented-points. Around this fit, we design a novel, local null-space voting strategy to reduce the 4-point case to 3. Voting is coupled with the famous RANSAC and makes our algorithm orders of magnitude faster than its conventional counterparts. This is the first method capable of performing a generic cross-type multi-object primitive detection in difficult scenes. Results on synthetic and real datasets support the validity of our method.Comment: Accepted for publication at CVPR 201

Restoration of the cantilever bowing distortion in Atomic Force Microscopy

Due to the mechanics of the Atomic Force Microscope (AFM), there is a curvature distortion (bowing effect) present in the acquired images. At present, flattening such images requires human intervention to manually segment object data from the background, which is time consuming and highly inaccurate. In this paper, an automated algorithm to flatten lines from AFM images is presented. The proposed method classifies the data into objects and background, and fits convex lines in an iterative fashion. Results on real images from DNA wrapped carbon nanotubes (DNACNTs) and synthetic experiments are presented, demonstrating the effectiveness of the proposed algorithm in increasing the resolution of the surface topography. In addition a link between the flattening problem and MRI inhomogeneity (shading) is given and the proposed method is compared to an entropy based MRI inhomogeniety correction method

BENDING THE DOMING EFFECT IN STRUCTURE FROM MOTION RECONSTRUCTIONS THROUGH BUNDLE ADJUSTMENT

Structure from Motion techniques provides low-cost and flexible methods that can be adopted in arial surveying to collect topographic data with accurate results. Nevertheless, the so-called "doming effect", due to unfortunate acquisition conditions or unreliable modeling of radial distortion, has been recognized as a critical issue that disrupts the quality of the attained 3D reconstruction. In this paper we propose a novel method, that works effectively in the presence of a nearly flat soil, to tackle a posteriori the doming effect: an automatic ground detection method is used to capture the doming deformation flawing the reconstruction, which in turn is wrapped to the correct geometry by iteratively enforcing a planarity constraint through a Bundle Adjustment framework. Experiments on real word datasets demonstrate promising results

The Three-Dimensional Circumstellar Environment of SN 1987A

We present the detailed construction and analysis of the most complete map to date of the circumstellar environment around SN 1987A, using ground and space-based imaging from the past 16 years. PSF-matched difference-imaging analyses of data from 1988 through 1997 reveal material between 1 and 28 ly from the SN. Careful analyses allows the reconstruction of the probable circumstellar environment, revealing a richly-structured bipolar nebula. An outer, double-lobed ``Peanut,'' which is believed to be the contact discontinuity between red supergiant and main sequence winds, is a prolate shell extending 28 ly along the poles and 11 ly near the equator. Napoleon's Hat, previously believed to be an independent structure, is the waist of this Peanut, which is pinched to a radius of 6 ly. Interior to this is a cylindrical hourglass, 1 ly in radius and 4 ly long, which connects to the Peanut by a thick equatorial disk. The nebulae are inclined 41\degr south and 8\degr east of the line of sight, slightly elliptical in cross section, and marginally offset west of the SN. From the hourglass to the large, bipolar lobes, echo fluxes suggest that the gas density drops from 1--3 cm^{-3} to >0.03 cm^{-3}, while the maximum dust-grain size increases from ~0.2 micron to 2 micron, and the Si:C dust ratio decreases. The nebulae have a total mass of ~1.7 Msun. The geometry of the three rings is studied, suggesting the northern and southern rings are located 1.3 and 1.0 ly from the SN, while the equatorial ring is elliptical (b/a < 0.98), and spatially offset in the same direction as the hourglass.Comment: Accepted for publication in the ApJ Supplements. 38 pages in apjemulate format, with 52 figure

Imaging Fabry-Perot Spectroscopy of NGC 5775: Kinematics of the Diffuse Ionized Gas Halo

We present imaging Fabry-Perot observations of Halpha emission in the nearly edge-on spiral galaxy NGC 5775. We have derived a rotation curve and a radial density profile along the major axis by examining position-velocity (PV) diagrams from the Fabry-Perot data cube as well as a CO 2-1 data cube from the literature. PV diagrams constructed parallel to the major axis are used to examine changes in azimuthal velocity as a function of height above the midplane. The results of this analysis reveal the presence of a vertical gradient in azimuthal velocity. The magnitude of this gradient is approximately 1 km/s/arcsec, or about 8 km/s/kpc, though a higher value of the gradient may be appropriate in localized regions of the halo. The evidence for an azimuthal velocity gradient is much stronger for the approaching half of the galaxy, although earlier slit spectra are consistent with a gradient on both sides. There is evidence for an outward radial redistribution of gas in the halo. The form of the rotation curve may also change with height, but this is not certain. We compare these results with those of an entirely ballistic model of a disk-halo flow. The model predicts a vertical gradient in azimuthal velocity which is shallower than the observed gradient, indicating that an additional mechanism is required to further slow the rotation speeds in the halo.Comment: 18 pages, 18 figures. Uses emulateapj.cls. Accepted for publication in Ap

Proper motions of the HH1 jet

We describe a new method for determining proper motions of extended objects, and a pipeline developed for the application of this method. We then apply this method to an analysis of four epochs of [S~II] HST images of the HH~1 jet (covering a period of $\sim 20$~yr). We determine the proper motions of the knots along the jet, and make a reconstruction of the past ejection velocity time-variability (assuming ballistic knot motions). This reconstruction shows an "acceleration" of the ejection velocities of the jet knots, with higher velocities at more recent times. This acceleration will result in an eventual merging of the knots in $\sim 450$~yr and at a distance of $\sim 80"$ from the outflow source, close to the present-day position of HH~1.Comment: 12 pages, 8 figure