A new coarse-to-fine strategy was proposed for nonrigid registration of computed tomography(CT) and magnetic resonance(MR) images of a liver.This hierarchical framework consisted of an affine transformation and a B-sp...A new coarse-to-fine strategy was proposed for nonrigid registration of computed tomography(CT) and magnetic resonance(MR) images of a liver.This hierarchical framework consisted of an affine transformation and a B-splines free-form deformation(FFD).The affine transformation performed a rough registration targeting the mismatch between the CT and MR images.The B-splines FFD transformation performed a finer registration by correcting local motion deformation.In the registration algorithm,the normalized mutual information(NMI) was used as similarity measure,and the limited memory Broyden-Fletcher- Goldfarb-Shannon(L-BFGS) optimization method was applied for optimization process.The algorithm was applied to the fully automated registration of liver CT and MR images in three subjects.The results demonstrate that the proposed method not only significantly improves the registration accuracy but also reduces the running time,which is effective and efficient for nonrigid registration.展开更多
An intensity-based non-rigid registration algorithm is discussed, which uses Gaussian smoothing to constrain the transformation to be smooth, and thus preserves the topology of images. In view of the insufficiency of ...An intensity-based non-rigid registration algorithm is discussed, which uses Gaussian smoothing to constrain the transformation to be smooth, and thus preserves the topology of images. In view of the insufficiency of the uniform Gaussian filtering of the deformation field, an automatic and accurate non-rigid image registration method based on B-splines approximation is proposed. The regularization strategy is adopted by using multi-level B-splines approximation to regularize the displacement fields in a coarse-to-fine manner. Moreover, it assigns the different weights to the estimated displacements according to their reliabilities. In this way, the level of regularity can be adapted locally. Experiments were performed on both synthetic and real medical images of brain, and the results show that the proposed method improves the registration accuracy and robustness.展开更多
Methods of digital human modeling have been developed and utilized to reflect human shape features.However,most of published works focused on dynamic visualization or fashion design,instead of high-accuracy modeling,w...Methods of digital human modeling have been developed and utilized to reflect human shape features.However,most of published works focused on dynamic visualization or fashion design,instead of high-accuracy modeling,which was strongly demanded by medical or rehabilitation scenarios.Prior to a high-accuracy modeling of human legs based on non-uniform rational B-splines(NURBS),the method of extracting the required quasi-grid network of feature points for human legs is presented in this work.Given the 3 D scanned human body,the leg is firstly segmented and put in standardized position.Then re-sampling of the leg is conducted via a set of equidistant cross sections.Through analysis of leg circumferences and circumferential curvature,the characteristic sections of the leg as well as the characteristic points on the sections are then identified according to the human anatomy and shape features.The obtained collection can be arranged to form a grid of data points for knots calculation and high-accuracy shape reconstruction in future work.展开更多
Polynomial splines have played an important role in image processing, medical imaging and wavelet theory. Exponential splines which are of more general concept have been recently investigated.We focus on cardinal expo...Polynomial splines have played an important role in image processing, medical imaging and wavelet theory. Exponential splines which are of more general concept have been recently investigated.We focus on cardinal exponential splines and develop a method to implement the exponential B-splines which form a Riesz basis of the space of cardinal exponential splines with finite energy.展开更多
A time integration algorithm for structural dynamic analysis is proposed by uniform cubic B-spline functions. The proposed algorithm is successfully used to solve the dynamic response of a single degree of freedom （S...A time integration algorithm for structural dynamic analysis is proposed by uniform cubic B-spline functions. The proposed algorithm is successfully used to solve the dynamic response of a single degree of freedom （SDOF） system, and then is generalized for a multiple-degree of freedom （MDOF） system. Stability analysis shows that, with an adjustable algorithmic parameter, the proposed method can achieve both conditional and unconditional stabilities. Validity of the method is shown with four numerical simulations. Comparison between the proposed method and other methods shows that the proposed method possesses high computation accuracy and desirable computation efficiency.展开更多
This paper describes an approximating solution, based on Lagrange interpolation and spline functions, to treat functional integral equations of Fredholm type and Volterra type. This method extended to functional integ...This paper describes an approximating solution, based on Lagrange interpolation and spline functions, to treat functional integral equations of Fredholm type and Volterra type. This method extended to functional integral and integro-differential equations. For showing efficiency of the method we give some numerical examples.展开更多
A new kind of spline with variable frequencies, called ωB-spline, is presented. It not only unifies B-splines, trigonometric and hyperbolic polynomial B-splines, but also produces more new types of splines, ωB-splin...A new kind of spline with variable frequencies, called ωB-spline, is presented. It not only unifies B-splines, trigonometric and hyperbolic polynomial B-splines, but also produces more new types of splines, ωB-spline bases are defined in the space spanned by {coso） t, sino）t, ], t, ..., t^n, ...} with the sequence of frequencies m where n is an arbitrary nonnegative integer, ωB-splines persist all desirable properties of B-splines. Furthermore, they have some special properties advantageous for modeling free form curves and surfaces.展开更多
is not completely clear which elements constitute the frame sets of the B-splines currently,but some considerable results have been obtained.In this paper,firstly,the background of frame set is introduced.Secondly,the...is not completely clear which elements constitute the frame sets of the B-splines currently,but some considerable results have been obtained.In this paper,firstly,the background of frame set is introduced.Secondly,the main progress of the frame sets of the B-splines in the past more than twenty years are reviewed,and particularly the progress for the frame set of the 2 order Bspline and the frame set of the 3 order B-spline are explained,respectively.展开更多
Triangular B-splines are powerful and flexible in modeling a broader class of geometric objects defined over arbitrary, non-rectangular domains. Despite their great potential and advantages in theory, practical techni...Triangular B-splines are powerful and flexible in modeling a broader class of geometric objects defined over arbitrary, non-rectangular domains. Despite their great potential and advantages in theory, practical techniques and computational tools with triangular B-splines are less-developed. This is mainly because users have to handle a large number of irregularly distributed control points over arbitrary triangulation. In this paper, an automatic and efficient method is proposed to generate visually pleasing, high-quality triangular B-splines of arbitrary topology. The experimental results on several real datasets show that triangular B-splines are powerful and effective in both theory and practice.展开更多
Abstract Generalized B-splines have been employed as geometric modeling and numerical simu- lation tools for isogeometric analysis （IGA for short）. However, the previous models used in IGA, such as trigonometric gen...Abstract Generalized B-splines have been employed as geometric modeling and numerical simu- lation tools for isogeometric analysis （IGA for short）. However, the previous models used in IGA, such as trigonometric generalized B-splines or hyperbolic generalized B-splines, are not the unified mathematical representation of conics and polynomial parametric curves/surfaces. In this paper, a unified approach to construct the generalized non-uniform B-splines over the space spanned by {α（t）,β（t）,ξ（t）, η（t）, 1, t,……. , tn-4} is proposed, and the corresponding isogeometric analysis framework for PDE solving is also studied. Compared with the NURBS-IGA method, the proposed frameworks have several advantages such as high accuracy, easy-to-compute derivatives and integrals due to the non-rational form. Furthermore, with the proposed spline models, isogeometric analysis can be performed on the computational domain bounded by transcendental curves/surfaces, such as the involute of circle, the helix/helicoid, the catenary/catenoid and the cycloid. Several numerical examples for isogeometrie heat conduction problems are presented to show the effectiveness of the proposed methods.展开更多
Regression splines are often used for fitting nonparametric functions, and they work especially well for additivity models. In this paper, we consider two simple tests of additivity: an adaptation of Tukey’s one degr...Regression splines are often used for fitting nonparametric functions, and they work especially well for additivity models. In this paper, we consider two simple tests of additivity: an adaptation of Tukey’s one degree of freedom test and a nonparametric version of Rao’s score test. While the Tukey-type test can detect most forms of the local non-additivity at the parametric rate of O(n-1/2), the score test is consistent for all alternative at a nonparametric rate. The asymptotic distribution of these test statistics is derived under both the null and local alternative hypotheses. A simulation study is conducted to compare their finite-sample performances with some existing kernel-based tests. The score test is found to have a good overall performance.展开更多
A numerical computation method based on B-splines is presented for thehydrodynamic interaction forces between two ships. The B-spline functions are adopted to approximatethe fully three-dimensional ship hull geometry ...A numerical computation method based on B-splines is presented for thehydrodynamic interaction forces between two ships. The B-spline functions are adopted to approximatethe fully three-dimensional ship hull geometry and unknown velocity potential in the fluid domain.The results and analysis are given in detail for the hydrodynamic interaction forces between twoships moving on parallel courses. All the computations show that the numerical results are in goodagreement with some experimental or other theoretical results.展开更多
This work puts forward an explicit isogeometric topology optimization(ITO)method using moving morphable components(MMC),which takes the suitably graded truncated hierarchical B-Spline based isogeometric analysis as th...This work puts forward an explicit isogeometric topology optimization(ITO)method using moving morphable components(MMC),which takes the suitably graded truncated hierarchical B-Spline based isogeometric analysis as the solver of physical unknown(SGTHB-ITO-MMC).By applying properly basis graded constraints to the hierarchical mesh of truncated hierarchical B-splines(THB),the convergence and robustness of the SGTHB-ITOMMC are simultaneously improved and the tiny holes occurred in optimized structure are eliminated,due to the improved accuracy around the explicit structural boundaries.Moreover,an efficient computational method is developed for the topological description functions(TDF)ofMMC under the admissible hierarchicalmesh,which consists of reducing the dimensionality strategy for design space and the locally computing strategy for hierarchical mesh.We apply the above SGTHB-ITO-MMC with improved efficiency to a series of 2D and 3Dcompliance design problems.The numerical results show that the proposed SGTHB-ITO-MMC method outperforms the traditional THB-ITO-MMCmethod in terms of convergence rate and efficiency.Therefore,the proposed SGTHB-ITO-MMC is an effective way of solving topology optimization(TO)problems.展开更多
基金Project(61240010)supported by the National Natural Science Foundation of ChinaProject(20070007070)supported by Specialized Research Fund for the Doctoral Program of Higher Education of China
文摘A new coarse-to-fine strategy was proposed for nonrigid registration of computed tomography(CT) and magnetic resonance(MR) images of a liver.This hierarchical framework consisted of an affine transformation and a B-splines free-form deformation(FFD).The affine transformation performed a rough registration targeting the mismatch between the CT and MR images.The B-splines FFD transformation performed a finer registration by correcting local motion deformation.In the registration algorithm,the normalized mutual information(NMI) was used as similarity measure,and the limited memory Broyden-Fletcher- Goldfarb-Shannon(L-BFGS) optimization method was applied for optimization process.The algorithm was applied to the fully automated registration of liver CT and MR images in three subjects.The results demonstrate that the proposed method not only significantly improves the registration accuracy but also reduces the running time,which is effective and efficient for nonrigid registration.
基金Supported by National Natural Science Foundation of China (No60373061)Joint Programof National Natural Science Foundation of ChinaGeneral Administration of Civil Aviation of China (No60672168)
文摘An intensity-based non-rigid registration algorithm is discussed, which uses Gaussian smoothing to constrain the transformation to be smooth, and thus preserves the topology of images. In view of the insufficiency of the uniform Gaussian filtering of the deformation field, an automatic and accurate non-rigid image registration method based on B-splines approximation is proposed. The regularization strategy is adopted by using multi-level B-splines approximation to regularize the displacement fields in a coarse-to-fine manner. Moreover, it assigns the different weights to the estimated displacements according to their reliabilities. In this way, the level of regularity can be adapted locally. Experiments were performed on both synthetic and real medical images of brain, and the results show that the proposed method improves the registration accuracy and robustness.
基金National Natural Science Foundation of China(Nos.12002085 and 51603039)Shanghai Pujiang Program,China(No.19PC002)+1 种基金Fundamental Research Funds for the Central Universities,China(No.2232019D3-58)Initial Research Funds for Young Teachers of Donghua University,China(No.104-07-0053088)。
文摘Methods of digital human modeling have been developed and utilized to reflect human shape features.However,most of published works focused on dynamic visualization or fashion design,instead of high-accuracy modeling,which was strongly demanded by medical or rehabilitation scenarios.Prior to a high-accuracy modeling of human legs based on non-uniform rational B-splines(NURBS),the method of extracting the required quasi-grid network of feature points for human legs is presented in this work.Given the 3 D scanned human body,the leg is firstly segmented and put in standardized position.Then re-sampling of the leg is conducted via a set of equidistant cross sections.Through analysis of leg circumferences and circumferential curvature,the characteristic sections of the leg as well as the characteristic points on the sections are then identified according to the human anatomy and shape features.The obtained collection can be arranged to form a grid of data points for knots calculation and high-accuracy shape reconstruction in future work.
文摘Polynomial splines have played an important role in image processing, medical imaging and wavelet theory. Exponential splines which are of more general concept have been recently investigated.We focus on cardinal exponential splines and develop a method to implement the exponential B-splines which form a Riesz basis of the space of cardinal exponential splines with finite energy.
基金Project supported by the National Natural Science Foundation of China(Nos.11602004 and11602081)the Fundamental Research Funds for the Central Universities(No.531107040934)
文摘A time integration algorithm for structural dynamic analysis is proposed by uniform cubic B-spline functions. The proposed algorithm is successfully used to solve the dynamic response of a single degree of freedom （SDOF） system, and then is generalized for a multiple-degree of freedom （MDOF） system. Stability analysis shows that, with an adjustable algorithmic parameter, the proposed method can achieve both conditional and unconditional stabilities. Validity of the method is shown with four numerical simulations. Comparison between the proposed method and other methods shows that the proposed method possesses high computation accuracy and desirable computation efficiency.
文摘This paper describes an approximating solution, based on Lagrange interpolation and spline functions, to treat functional integral equations of Fredholm type and Volterra type. This method extended to functional integral and integro-differential equations. For showing efficiency of the method we give some numerical examples.
基金the National Natural Science Foundation of China(Grant No.60773179)Foundation of State Key Basic Research 973 Development Programming Item of China(Grant No.G2004CB318000)
文摘A new kind of spline with variable frequencies, called ωB-spline, is presented. It not only unifies B-splines, trigonometric and hyperbolic polynomial B-splines, but also produces more new types of splines, ωB-spline bases are defined in the space spanned by {coso） t, sino）t, ], t, ..., t^n, ...} with the sequence of frequencies m where n is an arbitrary nonnegative integer, ωB-splines persist all desirable properties of B-splines. Furthermore, they have some special properties advantageous for modeling free form curves and surfaces.
基金supported in part by the National Natural Science Foundation of China(Grant No.61471410).
文摘is not completely clear which elements constitute the frame sets of the B-splines currently,but some considerable results have been obtained.In this paper,firstly,the background of frame set is introduced.Secondly,the main progress of the frame sets of the B-splines in the past more than twenty years are reviewed,and particularly the progress for the frame set of the 2 order Bspline and the frame set of the 3 order B-spline are explained,respectively.
文摘Triangular B-splines are powerful and flexible in modeling a broader class of geometric objects defined over arbitrary, non-rectangular domains. Despite their great potential and advantages in theory, practical techniques and computational tools with triangular B-splines are less-developed. This is mainly because users have to handle a large number of irregularly distributed control points over arbitrary triangulation. In this paper, an automatic and efficient method is proposed to generate visually pleasing, high-quality triangular B-splines of arbitrary topology. The experimental results on several real datasets show that triangular B-splines are powerful and effective in both theory and practice.
基金supported by Zhejiang Provincial Natural Science Foundation of China under Grant No.LR16F020003the National Nature Science Foundation of China under Grant Nos.61472111,61602138+1 种基金the Open Project Program of the State Key Lab of CAD&CG(A1703)Zhejiang University
文摘Abstract Generalized B-splines have been employed as geometric modeling and numerical simu- lation tools for isogeometric analysis （IGA for short）. However, the previous models used in IGA, such as trigonometric generalized B-splines or hyperbolic generalized B-splines, are not the unified mathematical representation of conics and polynomial parametric curves/surfaces. In this paper, a unified approach to construct the generalized non-uniform B-splines over the space spanned by {α（t）,β（t）,ξ（t）, η（t）, 1, t,……. , tn-4} is proposed, and the corresponding isogeometric analysis framework for PDE solving is also studied. Compared with the NURBS-IGA method, the proposed frameworks have several advantages such as high accuracy, easy-to-compute derivatives and integrals due to the non-rational form. Furthermore, with the proposed spline models, isogeometric analysis can be performed on the computational domain bounded by transcendental curves/surfaces, such as the involute of circle, the helix/helicoid, the catenary/catenoid and the cycloid. Several numerical examples for isogeometrie heat conduction problems are presented to show the effectiveness of the proposed methods.
基金This work was partially supported by the National Natural Science Foundation of China (Grant No. 10231030)the Excellent Young Teacher Foundation of Education Ministry of China and University of Illinois Campus Research Board and by NSF Award SBR-9617278 and DMS-0102411
文摘Regression splines are often used for fitting nonparametric functions, and they work especially well for additivity models. In this paper, we consider two simple tests of additivity: an adaptation of Tukey’s one degree of freedom test and a nonparametric version of Rao’s score test. While the Tukey-type test can detect most forms of the local non-additivity at the parametric rate of O(n-1/2), the score test is consistent for all alternative at a nonparametric rate. The asymptotic distribution of these test statistics is derived under both the null and local alternative hypotheses. A simulation study is conducted to compare their finite-sample performances with some existing kernel-based tests. The score test is found to have a good overall performance.
文摘A numerical computation method based on B-splines is presented for thehydrodynamic interaction forces between two ships. The B-spline functions are adopted to approximatethe fully three-dimensional ship hull geometry and unknown velocity potential in the fluid domain.The results and analysis are given in detail for the hydrodynamic interaction forces between twoships moving on parallel courses. All the computations show that the numerical results are in goodagreement with some experimental or other theoretical results.
基金supported by the National Key R&D Program of China (2020YFB1708300)the Project funded by the China Postdoctoral Science Foundation (2021M701310).
文摘This work puts forward an explicit isogeometric topology optimization(ITO)method using moving morphable components(MMC),which takes the suitably graded truncated hierarchical B-Spline based isogeometric analysis as the solver of physical unknown(SGTHB-ITO-MMC).By applying properly basis graded constraints to the hierarchical mesh of truncated hierarchical B-splines(THB),the convergence and robustness of the SGTHB-ITOMMC are simultaneously improved and the tiny holes occurred in optimized structure are eliminated,due to the improved accuracy around the explicit structural boundaries.Moreover,an efficient computational method is developed for the topological description functions(TDF)ofMMC under the admissible hierarchicalmesh,which consists of reducing the dimensionality strategy for design space and the locally computing strategy for hierarchical mesh.We apply the above SGTHB-ITO-MMC with improved efficiency to a series of 2D and 3Dcompliance design problems.The numerical results show that the proposed SGTHB-ITO-MMC method outperforms the traditional THB-ITO-MMCmethod in terms of convergence rate and efficiency.Therefore,the proposed SGTHB-ITO-MMC is an effective way of solving topology optimization(TO)problems.
