A mutual information-based non-rigid medical image registration algorithm is presented. An approximate function of Hanning windowed sinc is used as kernel function of partial volume (PV) interpolation to estimate the ...A mutual information-based non-rigid medical image registration algorithm is presented. An approximate function of Hanning windowed sinc is used as kernel function of partial volume (PV) interpolation to estimate the joint histogram, which is the key to calculating the mutual information. And a new method is proposed to compute the gradient of mutual information with respect to the model parameters. The transformation of object is modeled by a free-form deformation (FFD) based on B-splines. The experiments on 3D synthetic and real image data show that the algorithm can converge at the global optimum and restrain the emergency of local extreme.展开更多
In this paper we present a C-1 interpolation scheme on a triangle. The interpolant assumes given values and one order derivatives at the vertices of the triangle. It is made up of partial interpolants blended with cor...In this paper we present a C-1 interpolation scheme on a triangle. The interpolant assumes given values and one order derivatives at the vertices of the triangle. It is made up of partial interpolants blended with corresponding weight functions. Any partial interpolant is a piecewise cubics defined on a split of the triangle, while the weight function is just the respective barycentric coordinate. Hence the interpolant can be regarded as a piecewise quartic. We device a simple algorithm for the evaluation of the interpolant. It's easy to represent the interpolant with B-net method. We also depict the Franke's function and its interpolant, the illustration of which shows good visual effect of the scheme.展开更多
In this paper, we consider the higher divided difference of a composite function f(g(t)) in which g(t) is an s-dimensional vector. By exploiting some properties from mixed partial divided differences and multiva...In this paper, we consider the higher divided difference of a composite function f(g(t)) in which g(t) is an s-dimensional vector. By exploiting some properties from mixed partial divided differences and multivariate Newton interpolation, we generalize the divided difference form of Faà di Bruno's formula with a scalar argument. Moreover, a generalized Faà di Bruno's formula with a vector argument is derived.展开更多
基金Supported bythe National Basic Research Programof China ("973"Program) (No2003CB716103)Key Project of Shanghai Scienceand Technology Committee(No05DZ19509)
文摘A mutual information-based non-rigid medical image registration algorithm is presented. An approximate function of Hanning windowed sinc is used as kernel function of partial volume (PV) interpolation to estimate the joint histogram, which is the key to calculating the mutual information. And a new method is proposed to compute the gradient of mutual information with respect to the model parameters. The transformation of object is modeled by a free-form deformation (FFD) based on B-splines. The experiments on 3D synthetic and real image data show that the algorithm can converge at the global optimum and restrain the emergency of local extreme.
文摘In this paper we present a C-1 interpolation scheme on a triangle. The interpolant assumes given values and one order derivatives at the vertices of the triangle. It is made up of partial interpolants blended with corresponding weight functions. Any partial interpolant is a piecewise cubics defined on a split of the triangle, while the weight function is just the respective barycentric coordinate. Hence the interpolant can be regarded as a piecewise quartic. We device a simple algorithm for the evaluation of the interpolant. It's easy to represent the interpolant with B-net method. We also depict the Franke's function and its interpolant, the illustration of which shows good visual effect of the scheme.
基金Acknowledgments. This work was supported by the National Science Foundation of China (Grant Nos. 10471128, 10731060).
文摘In this paper, we consider the higher divided difference of a composite function f(g(t)) in which g(t) is an s-dimensional vector. By exploiting some properties from mixed partial divided differences and multivariate Newton interpolation, we generalize the divided difference form of Faà di Bruno's formula with a scalar argument. Moreover, a generalized Faà di Bruno's formula with a vector argument is derived.