Discover the most talked about and latest scientific content & concepts.

Concept: Splines


This study aims to develop a new optical design method of soft multifocal contact lens (CLs) to obtain uniform optical power in large center-distance zone with optimized Non-Uniform Rational B-spline (NURBS). For the anterior surface profiles of CLs, the NURBS design curves are optimized to match given optical power distributions. Then, the NURBS in the center-distance zones are fitted in the corresponding spherical/aspheric curves for both data points and their centers of curvature to achieve the uniform power. Four cases of soft CLs have been manufactured by casting in shell molds by injection molding and then measured to verify the design specifications. Results of power profiles of these CLs are concord with the given clinical requirements of uniform powers in larger center-distance zone. The developed optical design method has been verified for multifocal CLs design and can be further applied for production of soft multifocal CLs.

Concepts: Optics, Lens, Contact lens, Corrective lens, Focal length, Dioptre, Interpolation, Splines


Non-uniform rational B-spline (NURBS) surface fitting from data points is wildly used in the fields of computer aided design (CAD), medical imaging, cultural relic representation and object-shape detection. Usually, the measured data acquired from coordinate measuring systems is neither gridded nor completely scattered. The distribution of this kind of data is scattered in physical space, but the data points are stored in a way consistent with the order of measurement, so it is named quasi scattered data in this paper. Therefore they can be organized into rows easily but the number of points in each row is random. In order to overcome the difficulty of surface fitting from this kind of data, a new method based on resampling is proposed. It consists of three major steps: (1) NURBS curve fitting for each row, (2) resampling on the fitted curve and (3) surface fitting from the resampled data. Iterative projection optimization scheme is applied in the first and third step to yield advisable parameterization and reduce the time cost of projection. A resampling approach based on parameters, local peaks and contour curvature is proposed to overcome the problems of nodes redundancy and high time consumption in the fitting of this kind of scattered data. Numerical experiments are conducted with both simulation and practical data, and the results show that the proposed method is fast, effective and robust. What’s more, by analyzing the fitting results acquired form data with different degrees of scatterness it can be demonstrated that the error introduced by resampling is negligible and therefore it is feasible.

Concepts: Measurement, Computer-aided design, Systems of measurement, Interpolation, Splines, Multivariate interpolation, Non-uniform rational B-spline


Background Discontinuity of spirometry reference values from childhood into adulthood has been a problem with traditional reference values, thus modern modelling approaches using smoothing spline functions to better depict the transition during growth and ageing have been recently introduced. Following the publication of the new international Global Lung Initiative (GLI2012) reference values also new national Finnish reference values have been calculated using similar GAMLSS-modelling, with spline estimates for mean (Mspline) and standard deviation (Sspline) provided in tables. The aim of this study was to produce polynomial estimates for these spline functions to use in lieu of lookup tables and to assess their validity in the reference population of healthy non-smokers. Methods Linear regression modelling was used to approximate the estimated values for Mspline and Sspline using similar polynomial functions as in the international GLI2012 reference values. Estimated values were compared to original calculations in absolute values, the derived predicted mean and individually calculated z-scores using both values. Results Polynomial functions were estimated for all 10 spirometry variables. The agreement between original lookup table-produced values and polynomial estimates was very good, with no significant differences found. The variation slightly increased in larger predicted volumes, but a range of -0.018 to +0.022 litres of FEV1 representing ± 0.4% of maximum difference in predicted mean. Conclusions Polynomial approximations were very close to the original lookup tables and are recommended for use in clinical practice to facilitate the use of new reference values.

Concepts: Regression analysis, Statistics, Mathematics, Standard deviation, Approximation, Estimation, Continuous function, Splines


The Smoothing Spline ANOVA (SS-ANOVA) requires a specialized construction of basis and penalty terms in order to incorporate prior knowledge about the data to be fitted. Typically, one resorts to the most general approach using tensor product splines. This implies severe constraints on the correlation structure, i.e. the assumption of isotropy of smoothness can not be incorporated in general. This may increase the variance of the spline fit, especially if only a relatively small set of observations are given. In this article, we propose an alternative method that allows to incorporate prior knowledge without the need to construct specialized bases and penalties, allowing the researcher to choose the spline basis and penalty according to the prior knowledge of the observations rather than choosing them according to the analysis to be done. The two approaches are compared with an artificial example and with analyses of fixation durations during reading.

Concepts: Scientific method, Variance, Mathematical analysis, Logic, Penalty, Assumption of Mary, Splines, Spline


The pubic symphysis is frequently used to estimate age-at-death from the adult skeleton. Assessment methods require the visual comparison of the bone morphology against age-informative characteristics that represent a series of phases. Age-at-death is then estimated from the age-range previously associated with the chosen phase. While easily executed, the “morphoscopic” process of feature-scoring and bone-to-phase-matching is known to be subjective. Studies of method and practitioner error demonstrate a need for alternative tools to quantify age-progressive change in the pubic symphysis. This article proposes a more objective, quantitative method that analyzes three-dimensional (3D) surface scans of the pubic symphysis using a thin plate spline algorithm (TPS).

Concepts: Scientific method, Mathematics, Approximation, Estimation, Pubic symphysis, Splines, Thin plate spline


Introduction. It is well known that resolution on a gamma camera varies as a function of distance, scatter and the camera’s characteristics (collimator type, crystal thickness, intrinsic resolution etc). Manufacturers frequently provide only a few pre-calculated resolution values (using a line source in air, 10-15 cm from the collimator surface and without scattering). However, these are typically not obtained in situations resembling a clinical setting. From a diagnostic point of view, it is useful to know the expected resolution of a gamma camera at a given distance from the collimator surface for a particular setting in order to decide whether it is worth scanning patients with “small lesion” or not. When dealing with absolute quantification it is also mandatory to know precisely the expected resolution and its uncertainty in order to make appropriate corrections. Aim. Our aims are: to test a novel mathematical approach, the cubic spline interpolation, for the extraction of the full width at half maximum (FWHM) from the acquisition of a line source (experimental resolution) also considering measurement uncertainty; to compare it with the usually adopted methods such as the gaussian approach; to compare it with the theoretical resolution (analytical resolution) of a gamma camera at different distances; to create a web-based educational program with which to test these theories. Methods. Three mathematical methods (direct calculation, global interpolation using gaussian and local interpolation using splines) for calculating FWHM from a line source (planar scintigraphy) were tested and compared. A NEMA Triple Line Source Phantom was used to obtain static images both in air and with different scattering levels. An advanced, open-source software (MATLAB/Octave and PHP based) was created “ad hoc” to obtain and compare FWHM values and relative uncertainty. Results and Conclusion. Local interpolation using splines proved faster and more reliable than the usually-adopted Gaussian interpolation. The proposed freely available software proved effective in assessing both FWHM and its uncertainty.

Concepts: Mathematics, Measurement, Distance, Spline interpolation, Interpolation, Splines, Polynomial interpolation, Spline


A new horn with high displacement amplification for ultrasonic welding is developed. The profile of the horn is a nonrational B-spline curve with an open uniform knot vector. The ultrasonic actuation of the horn exploits the first longitudinal displacement mode of the horn. The horn is designed by an optimization scheme and finite element analyses. Performances of the proposed horn have been evaluated by experiments. The displacement amplification of the proposed horn is 41.4% and 8.6% higher than that of the traditional catenoidal horn and a Bézier-profile horn, respectively, with the same length and end surface diameters. The developed horn has a lower displacement amplification than the nonuniform rational B-spline profiled horn but a much smoother stress distribution. The developed horn, the catenoidal horn, and the Bézier horn are fabricated and used for ultrasonic welding of lap-shear specimens. The bonding strength of the joints welded by the open uniform nonrational B-spline (OUNBS) horn is the highest among the three horns for the various welding parameters considered. The locations of the failure mode and the distribution of the voids of the specimens are investigated to explain the reason of the high bonding strength achieved by the OUNBS horn.

Concepts: Oxygen, Length, Metric space, Interpolation, Welding, Splines, B-spline, Ultrasonic welding


A set of fluence-to-effective dose conversion coefficients of external exposure to muons were investigated for Chinese hybrid phantom references, which include both male and female. Both polygon meshes and Non-Uniform Rational B-Spline (NURBS) surfaces were used to descried the boundary of the organs and tissues in these phantoms. The 3D-DOCTOR and Rhinoceros software were used to polygonise the colour slice images and generate the NURBS surfaces, respectively. The voxelisation is completed using the BINVOX software and the assembly finished by using MATLAB codes. The voxel resolutions were selected to be 0.22 × 0.22 × 0.22 cm(3) and 0.2 × 0.2 × 0.2 cm(3) for male and female phantoms, respectively. All parts of the final phantoms were matched to their reference organ masses within a tolerance of ±5 %. The conversion coefficients for negative and positive muons were calculated with the FLUKA transport code. There were 21 external monoenergetic beams ranging from 0.01 GeV to 100 TeV in 5 different geometrical conditions of irradiation.

Concepts: Electron, Mass, 3D computer graphics, Joule, The Conversion, Interpolation, Splines, Non-uniform rational B-spline


Dimensionality reduction (DR) has been considered as one of the most significant tools for data analysis. One type of DR algorithms is based on latent variable models (LVM). LVM-based models can handle the preimage problem easily. In this paper we propose a new LVM-based DR model, named thin plate spline latent variable model (TPSLVM). Compared to the well-known Gaussian process latent variable model (GPLVM), our proposed TPSLVM is more powerful especially when the dimensionality of the latent space is low. Also, TPSLVM is robust to shift and rotation. This paper investigates two extensions of TPSLVM, i.e., the back-constrained TPSLVM (BC-TPSLVM) and TPSLVM with dynamics (TPSLVM-DM) as well as their combination BC-TPSLVM-DM. Experimental results show that TPSLVM and its extensions provide better data visualization and more efficient dimensionality reduction compared to PCA, GPLVM, ISOMAP, etc.

Concepts: Dimension, Data analysis, Machine learning, Data mining, Splines, Latent variable model, Latent variable models, Thin plate spline


Branching curves are a technique for modeling curves that change trajectory at a change (branching) point. Currently, the estimation framework is limited to independent data, and smoothing splines are used for estimation. This article aims to extend the branching curve framework to the longitudinal data setting where the branching point varies by subject. If the branching point is modeled as a random effect, then the longitudinal branching curve framework is a semiparametric nonlinear mixed effects model. Given existing issues with using random effects within a smoothing spline, we express the model as a B-spline based semiparametric nonlinear mixed effects model. Simple, clever smoothness constraints are enforced on the B-splines at the change point. The method is applied to Women’s Health data where we model the shape of the labor curve (cervical dilation measured longitudinally) before and after treatment with oxytocin (a labor stimulant). Copyright © 2014 John Wiley & Sons, Ltd.

Concepts: Random effects model, Interpolation, Splines, Smoothing, Spline, Cervical dilation, B-spline, Smoothing spline