Adaptive Parametrization of Multivariate B-splines for Image Registration
Hansen, Michael Sass; Glocker, Benjamin; Navab, Nassir;
2008-01-01
We present an adaptive parametrization scheme for dynamic mesh refinement in the application of parametric image registration. The scheme is based on a refinement measure ensuring that the control points give an efficient representation of the warp fields, in terms of minimizing the registration...... cost function. In the current work we introduce multivariate B-splines as a novel alternative to the widely used tensor B-splines enabling us to make efficient use of the derived measure.The multivariate B-splines of order n are Cn- 1 smooth and are based on Delaunay configurations of arbitrary 2D or 3...... reside on a regular grid. In contrast, by efficient non- constrained placement of the knots, the multivariate B- splines are shown to give a good representation of inho- mogeneous objects in natural settings. The wide applicability of the method is illustrated through its application on medical data and...
宋红; 李佳佳; 王树良; 马婧婷
2014-01-01
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.
Robust Nonrigid Multimodal Image Registration using Local Frequency Maps*
Jian, Bing; Vemuri, Baba C.; Marroquin, José L.
2008-01-01
Automatic multi-modal image registration is central to numerous tasks in medical imaging today and has a vast range of applications e.g., image guidance, atlas construction, etc. In this paper, we present a novel multi-modal 3D non-rigid registration algorithm where in 3D images to be registered are represented by their corresponding local frequency maps efficiently computed using the Riesz transform as opposed to the popularly used Gabor filters. The non-rigid registration between these local frequency maps is formulated in a statistically robust framework involving the minimization of the integral squared error a.k.a. L2E (L2 error). This error is expressed as the squared difference between the true density of the residual (which is the squared difference between the non-rigidly transformed reference and the target local frequency representations) and a Gaussian or mixture of Gaussians density approximation of the same. The non-rigid transformation is expressed in a B-spline basis to achieve the desired smoothness in the transformation as well as computational efficiency. The key contributions of this work are (i) the use of Riesz transform to achieve better efficiency in computing the local frequency representation in comparison to Gabor filter-based approaches, (ii) new mathematical model for local-frequency based non-rigid registration, (iii) analytic computation of the gradient of the robust non-rigid registration cost function to achieve efficient and accurate registration. The proposed non-rigid L2E-based registration is a significant extension of research reported in literature to date. We present experimental results for registering several real data sets with synthetic and real non-rigid misalignments. PMID:17354721
Nonrigid registration of myocardial perfusion MRI
Ólafsdóttir, Hildur
2005-01-01
This paper describes a fully automatic registration of 10 multi-slice myocardial perfusion magnetic resonance image sequences. The registration of these sequences is crucial for the clinical interpretation, which currently is subjected to manual labour. The approach used in this study is a nonrigid...... registration algorithm based on free-form deformations due to Rueckert et al. Inspection of difference images from the wash-out part of the perfusion sequences indicates that a good registration accuracy is obtained....
Nonrigid registration of myocardial perfusion MRI
Ólafsdóttir, Hildur
2005-01-01
This paper describes a fully automatic registration of 10 multi-slice myocardial perfusion magnetic resonance image sequences. The registration of these sequences is crucial for the clinical interpretation, which currently is subjected to manual labour. The approach used in this study is a nonrigid registration algorithm based on free-form deformations due to Rueckert et al. Inspection of difference images from the wash-out part of the perfusion sequences indicates that a good registration ac...
Efficient nonrigid registration using ranked order statistics
Tennakoon, R.B.; Bab-Hadiashar, A.; de Bruijne, Marleen;
2013-01-01
research. In this paper we propose a fast and accurate non-rigid registration method for intra-modality volumetric images. Our approach exploits the information provided by an order statistics based segmentation method, to find the important regions for registration and use an appropriate sampling scheme...
Non-Rigid Registration meets Surface Reconstruction
Rouhani, Mohammad; Boyer, Edmond; Angel D. Sappa
2014-01-01
Non rigid registration is an important task in computer vision with many applications in shape and motion modeling. A fundamental step of the registration is the data association between the source and the target sets. Such association proves difficult in practice, due to the discrete nature of the information and its corruption by various types of noise, e.g. outliers and missing data. In this paper we investigate the benefit of the implicit representations for the non-rigid registration of ...
Ji, Songbai; Fan, Xiaoyao; Roberts, David W.; Hartov, Alex; Paulsen, Keith D.
2011-03-01
Compensating for brain shift as surgery progresses is important to ensure sufficient accuracy in patient-to-image registration in the operating room (OR) for reliable neuronavigation. Ultrasound has emerged as an important and practical imaging technique for brain shift compensation either by itself or through computational modeling that estimates whole-brain deformation. Using volumetric true 3D ultrasound (3DUS), it is possible to nonrigidly (e.g., based on B-splines) register two temporally different 3DUS images directly to generate feature displacement maps for data assimilation in the biomechanical model. Because of a large amount of data and number of degrees-of-freedom (DOFs) involved, however, a significant computational cost may be required that can adversely influence the clinical feasibility of the technique for efficiently generating model-updated MR (uMR) in the OR. This paper parametrically investigates three B-splines registration parameters and their influence on the computational cost and registration accuracy: number of grid nodes along each direction, floating image volume down-sampling rate, and number of iterations. A simulated rigid body displacement field was employed as a ground-truth against which the accuracy of displacements generated from the B-splines nonrigid registration was compared. A set of optimal parameters was then determined empirically that result in a registration computational cost of less than 1 min and a sub-millimetric accuracy in displacement measurement. These resulting parameters were further applied to a clinical surgery case to demonstrate their practical use. Our results indicate that the optimal set of parameters result in sufficient accuracy and computational efficiency in model computation, which is important for future application of the overall biomechanical modeling to generate uMR for image-guidance in the OR.
Deformable image registration (DIR) is fundamental technique for adaptive radiotherapy and image-guided radiotherapy. However, further improvement of DIR is still needed. We evaluated the accuracy of B-spline transformation-based DIR implemented in elastix. This registration package is largely based on the Insight Segmentation and Registration Toolkit (ITK), and several new functions were implemented to achieve high DIR accuracy. The purpose of this study was to clarify whether new functions implemented in elastix are useful for improving DIR accuracy. Thoracic 4D computed tomography images of ten patients with esophageal or lung cancer were studied. Datasets for these patients were provided by DIR-lab (dir-lab.com) and included a coordinate list of anatomical landmarks that had been manually identified. DIR between peak-inhale and peak-exhale images was performed with four types of parameter settings. The first one represents original ITK (Parameter 1). The second employs the new function of elastix (Parameter 2), and the third was created to verify whether new functions improve DIR accuracy while keeping computational time (Parameter 3). The last one partially employs a new function (Parameter 4). Registration errors for these parameter settings were calculated using the manually determined landmark pairs. 3D registration errors with standard deviation over all cases were 1.78 (1.57), 1.28 (1.10), 1.44 (1.09) and 1.36 (1.35) mm for Parameter 1, 2, 3 and 4, respectively, indicating that the new functions are useful for improving DIR accuracy, even while maintaining the computational time, and this B-spline-based DIR could be used clinically to achieve high-accuracy adaptive radiotherapy. (author)
To explore the effects of computed tomography (CT) image characteristics and B-spline knot spacing (BKS) on the spatial accuracy of a B-spline deformable image registration (DIR) in the head-and-neck geometry. The effect of image feature content, image contrast, noise, and BKS on the spatial accuracy of a B-spline DIR was studied. Phantom images were created with varying feature content and varying contrast-to-noise ratio (CNR), and deformed using a known smooth B-spline deformation. Subsequently, the deformed images were repeatedly registered with the original images using different BKSs. The quality of the DIR was expressed as the mean residual displacement (MRD) between the known imposed deformation and the result of the B-spline DIR. Finally, for three patients, head-and-neck planning CT scans were deformed with a realistic deformation field derived from a rescan CT of the same patient, resulting in a simulated deformed image and an a-priori known deformation field. Hence, a B-spline DIR was performed between the simulated image and the planning CT at different BKSs. Similar to the phantom cases, the DIR accuracy was evaluated by means of MRD. In total, 162 phantom registrations were performed with varying CNR and BKSs. MRD-values < 1.0 mm were observed with a BKS between 10–20 mm for image contrast ≥ ± 250 HU and noise < ± 200 HU. Decreasing the image feature content resulted in increased MRD-values at all BKSs. Using BKS = 15 mm for the three clinical cases resulted in an average MRD < 1.0 mm. For synthetically generated phantoms and three real CT cases the highest DIR accuracy was obtained for a BKS between 10–20 mm. The accuracy decreased with decreasing image feature content, decreasing image contrast, and higher noise levels. Our results indicate that DIR accuracy in clinical CT images (typical noise levels < ± 100 HU) will not be effected by the amount of image noise
Non-rigid consistent registration of 2D image sequences
We present a novel algorithm for the registration of 2D image sequences that combines the principles of multiresolution B-spline-based elastic registration and those of bidirectional consistent registration. In our method, consecutive triples of images are iteratively registered to gradually extend the information through the set of images of the entire sequence. The intermediate results are reused for the registration of the following triple. We choose to interpolate the images and model the deformation fields using B-spline multiresolution pyramids. Novel boundary conditions are introduced to better characterize the deformations at the boundaries. In the experimental section, we quantitatively show that our method recovers from barrel/pincushion and fish-eye deformations with subpixel error. Moreover, it is more robust against outliers-occasional strong noise and large rotations-than the state-of-the-art methods. Finally, we show that our method can be used to realign series of histological serial sections, which are often heavily distorted due to folding and tearing of the tissues.
Non-rigid image registration using bone growth model
Bro-Nielsen, Morten; Gramkow, Claus; Kreiborg, Sven;
1997-01-01
Non-rigid registration has traditionally used physical models like elasticity and fluids. These models are very seldom valid models of the difference between the registered images. This paper presents a non-rigid registration algorithm, which uses a model of bone growth as a model of the change...... between time sequence images of the human mandible. By being able to register the images, this paper at the same time contributes to the validation of the growth model, which is based on the currently available medical theories and knowledge...
The goal of this study is to develop a matching algorithm that can handle large geometric changes in x-ray computed tomography (CT)-derived lung geometry occurring during deep breath maneuvers. These geometric relationships are further utilized to build a dynamic lung airway model for computational fluid dynamics (CFD) studies of pulmonary air flow. The proposed algorithm is based on a cubic B-spline-based hybrid registration framework that incorporates anatomic landmark information with intensity patterns. A sequence of invertible B-splines is composed in a multiresolution framework to ensure local invertibility of the large deformation transformation and a physiologically meaningful similarity measure is adopted to compensate for changes in voxel intensity due to inflation. Registrations are performed using the proposed approach to match six pairs of 3D CT human lung datasets. Results show that the proposed approach has the ability to match the intensity pattern and the anatomical landmarks, and ensure local invertibility for large deformation transformations. Statistical results also show that the proposed hybrid approach yields significantly improved results as compared with approaches using either landmarks or intensity alone.
Yin Youbing; Lin, Ching-Long [Department of Mechanical and Industrial Engineering, University of Iowa, Iowa City, IA 52242 (United States); Hoffman, Eric A [Department of Radiology, University of Iowa, Iowa City, IA 52242 (United States); Ding Kai; Reinhardt, Joseph M, E-mail: ching-long-lin@uiowa.edu [Department of Biomedical Engineering, University of Iowa, Iowa City, IA 52242 (United States)
2011-01-07
The goal of this study is to develop a matching algorithm that can handle large geometric changes in x-ray computed tomography (CT)-derived lung geometry occurring during deep breath maneuvers. These geometric relationships are further utilized to build a dynamic lung airway model for computational fluid dynamics (CFD) studies of pulmonary air flow. The proposed algorithm is based on a cubic B-spline-based hybrid registration framework that incorporates anatomic landmark information with intensity patterns. A sequence of invertible B-splines is composed in a multiresolution framework to ensure local invertibility of the large deformation transformation and a physiologically meaningful similarity measure is adopted to compensate for changes in voxel intensity due to inflation. Registrations are performed using the proposed approach to match six pairs of 3D CT human lung datasets. Results show that the proposed approach has the ability to match the intensity pattern and the anatomical landmarks, and ensure local invertibility for large deformation transformations. Statistical results also show that the proposed hybrid approach yields significantly improved results as compared with approaches using either landmarks or intensity alone.
Non-rigid registration by geometry-constrained diffusion
Andresen, Per Rønsholt; Nielsen, Mads
1999-01-01
Assume that only partial knowledge about a non-rigid registration is given so that certain point, curves, or surfaces in one 3D image map to certain points, curves, or surfaces in another 3D image. We are facing the aperture problem because along the curves and surfaces, point correspondences are...
Non-rigid registration by geometry-constrained diffusion
Andresen, Per Rønsholt; Nielsen, Mads
2001-01-01
Assume that only partial knowledge about a non-rigid registration is given: certain points, curves, or surfaces in one 3D image are known to map to certain points, curves, or surfaces in another 3D image. In trying to identify the non-rigid registration field, we face a generalized aperture problem...... by a geometry-constrained diffusion, which in a precise sense yields the simplest displacement field. The point registration obtained may be used for segmentation, growth modeling, shape analysis, or kinematic interpolation. The algorithm applies to geometrical objects of any dimensionality. We may...... thus keep any number of fiducial points, curves, and/or surfaces fixed while finding the simplest registration. Examples of inferred point correspondences in a synthetic example and a longitudinal growth study of the human mandible are given....
van der Put, R. W.; Kerkhof, E. M.; Raaymakers, B. W.; Jürgenliemk-Schulz, I. M.; Lagendijk, J. J. W.
2009-12-01
External beam radiation treatment for patients with cervical cancer is hindered by the relatively large motion of the target volume. A hybrid MRI-accelerator system makes it possible to acquire online MR images during treatment in order to correct for motion and deformation. To fully benefit from such a system, online delineation of the target volumes is necessary. The aim of this study is to investigate the accuracy of rigid, non-rigid and semi-automatic registrations of MR images for interfractional contour propagation in patients with cervical cancer. Registration using mutual information was performed on both bony anatomy and soft tissue. A B-spline transform was used for the non-rigid method. Semi-automatic registration was implemented with a point set registration algorithm on a small set of manual landmarks. Online registration was simulated by application of each method to four weekly MRI scans for each of 33 cervical cancer patients. Evaluation was performed by distance analysis with respect to manual delineations. The results show that soft-tissue registration significantly (P registration based on bony anatomy. A combination of user-assisted and non-rigid registration provides the best results with a median error of 3.2 mm (1.4-9.9 mm) compared to 5.9 mm (1.7-19.7 mm) with bone registration (P registration (P = 0.01). In a clinical setting, the benefit may be further increased when outliers can be removed by visual inspection of the online images. We conclude that for external beam radiation treatment of cervical cancer, online MRI imaging will allow target localization based on soft tissue visualization, which provides a significantly higher accuracy than localization based on bony anatomy. The use of limited user input to guide the registration increases overall accuracy. Additional non-rigid registration further reduces the propagation error and negates errors caused by small observer variations.
Learning Nonrigid Deformations for Constrained Multi-modal Image Registration
Onofrey, John A.; Staib, Lawrence H.; Papademetris, Xenophon
2013-01-01
We present a new strategy to constrain nonrigid registrations of multi-modal images using a low-dimensional statistical deformation model and test this in registering pre-operative and post-operative images from epilepsy patients. For those patients who may undergo surgical resection for treatment, the current gold-standard to identify regions of seizure involves craniotomy and implantation of intracranial electrodes. To guide surgical resection, surgeons utilize pre-op anat...
Non-rigid multi-modal registration on the GPU
Vetter, Christoph; Guetter, Christoph; Xu, Chenyang; Westermann, Rüdiger
2007-03-01
Non-rigid multi-modal registration of images/volumes is becoming increasingly necessary in many medical settings. While efficient registration algorithms have been published, the speed of the solutions is a problem in clinical applications. Harnessing the computational power of graphics processing unit (GPU) for general purpose computations has become increasingly popular in order to speed up algorithms further, but the algorithms have to be adapted to the data-parallel, streaming model of the GPU. This paper describes the implementation of a non-rigid, multi-modal registration using mutual information and the Kullback-Leibler divergence between observed and learned joint intensity distributions. The entire registration process is implemented on the GPU, including a GPU-friendly computation of two-dimensional histograms using vertex texture fetches as well as an implementation of recursive Gaussian filtering on the GPU. Since the computation is performed on the GPU, interactive visualization of the registration process can be done without bus transfer between main memory and video memory. This allows the user to observe the registration process and to evaluate the result more easily. Two hybrid approaches distributing the computation between the GPU and CPU are discussed. The first approach uses the CPU for lower resolutions and the GPU for higher resolutions, the second approach uses the GPU to compute a first approximation to the registration that is used as starting point for registration on the CPU using double-precision. The results of the CPU implementation are compared to the different approaches using the GPU regarding speed as well as image quality. The GPU performs up to 5 times faster per iteration than the CPU implementation.
Brownian Warps for Non-Rigid Registration
Nielsen, Mads; Johansen, Peter; Jackson, Andrew D.;
2008-01-01
A Brownian motion model in the group of diffeomorphisms has been introduced as inducing a least committed prior on warps. This prior is source-destination symmetric, fulfills a natural semi-group property for warps, and with probability 1 creates invertible warps. Using this as a least committed ...... images, and show that the obtained warps are also in practice source-destination symmetric and in an example on X-ray spine registration provides extrapolations from landmark point superior to those of spline solutions. Udgivelsesdato: July......A Brownian motion model in the group of diffeomorphisms has been introduced as inducing a least committed prior on warps. This prior is source-destination symmetric, fulfills a natural semi-group property for warps, and with probability 1 creates invertible warps. Using this as a least committed...... prior, we formulate a Partial Differential Equation for obtaining the maximally likely warp given matching constraints derived from the images. We solve for the free boundary conditions, and the bias toward smaller areas in the finite domain setting. Furthermore, we demonstrate the technique on 2D...
External beam radiation treatment for patients with cervical cancer is hindered by the relatively large motion of the target volume. A hybrid MRI-accelerator system makes it possible to acquire online MR images during treatment in order to correct for motion and deformation. To fully benefit from such a system, online delineation of the target volumes is necessary. The aim of this study is to investigate the accuracy of rigid, non-rigid and semi-automatic registrations of MR images for interfractional contour propagation in patients with cervical cancer. Registration using mutual information was performed on both bony anatomy and soft tissue. A B-spline transform was used for the non-rigid method. Semi-automatic registration was implemented with a point set registration algorithm on a small set of manual landmarks. Online registration was simulated by application of each method to four weekly MRI scans for each of 33 cervical cancer patients. Evaluation was performed by distance analysis with respect to manual delineations. The results show that soft-tissue registration significantly (P < 0.001) improves the accuracy of contour propagation compared to registration based on bony anatomy. A combination of user-assisted and non-rigid registration provides the best results with a median error of 3.2 mm (1.4-9.9 mm) compared to 5.9 mm (1.7-19.7 mm) with bone registration (P < 0.001) and 3.4 mm (1.3-19.1 mm) with non-rigid registration (P = 0.01). In a clinical setting, the benefit may be further increased when outliers can be removed by visual inspection of the online images. We conclude that for external beam radiation treatment of cervical cancer, online MRI imaging will allow target localization based on soft tissue visualization, which provides a significantly higher accuracy than localization based on bony anatomy. The use of limited user input to guide the registration increases overall accuracy. Additional non-rigid registration further reduces the propagation
An improved 3D shape context registration method for non-rigid surface registration
Xiao, Di; Zahra, David; Bourgeat, Pierrick; Berghofer, Paula; Acosta Tamayo, Oscar; Wimberley, Catriona; Gregoire, Marie-Claude; Salvado, Olivier
2010-03-01
3D shape context is a method to define matching points between similar shapes as a pre-processing step to non-rigid registration. The main limitation of the approach is point mismatching, which includes long geodesic distance mismatch and neighbors crossing mismatch. In this paper, we propose a topological structure verification method to correct the long geodesic distance mismatch and a correspondence field smoothing method to correct the neighbors crossing mismatch. A robust 3D shape context model is proposed and further combined with thin-plate spline model for non-rigid surface registration. The method was tested on phantoms and rat hind limb skeletons from micro CT images. The results from experiments on mouse hind limb skeletons indicate that the approach is robust.
Onofrey, John A; Staib, Lawrence H; Papademetris, Xenophon
2016-01-01
This paper describes a framework for learning a statistical model of non-rigid deformations induced by interventional procedures. We make use of this learned model to perform constrained non-rigid registration of pre-procedural and post-procedural imaging. We demonstrate results applying this framework to non-rigidly register post-surgical computed tomography (CT) brain images to pre-surgical magnetic resonance images (MRIs) of epilepsy patients who had intra-cranial electroencephalography electrodes surgically implanted. Deformations caused by this surgical procedure, imaging artifacts caused by the electrodes, and the use of multi-modal imaging data make non-rigid registration challenging. Our results show that the use of our proposed framework to constrain the non-rigid registration process results in significantly improved and more robust registration performance compared to using standard rigid and non-rigid registration methods. PMID:26900569
John A. Onofrey
2016-01-01
Full Text Available This paper describes a framework for learning a statistical model of non-rigid deformations induced by interventional procedures. We make use of this learned model to perform constrained non-rigid registration of pre-procedural and post-procedural imaging. We demonstrate results applying this framework to non-rigidly register post-surgical computed tomography (CT brain images to pre-surgical magnetic resonance images (MRIs of epilepsy patients who had intra-cranial electroencephalography electrodes surgically implanted. Deformations caused by this surgical procedure, imaging artifacts caused by the electrodes, and the use of multi-modal imaging data make non-rigid registration challenging. Our results show that the use of our proposed framework to constrain the non-rigid registration process results in significantly improved and more robust registration performance compared to using standard rigid and non-rigid registration methods.
Nonrigid registration of carotid ultrasound and MR images using a "twisting and bending" model
Nanayakkara, Nuwan D.; Chiu, Bernard; Samani, Abbas; Spence, J. David; Parraga, Grace; Samarabandu, Jagath; Fenster, Aaron
2008-03-01
Atherosclerosis at the carotid bifurcation resulting in cerebral emboli is a major cause of ischemic stroke. Most strokes associated with carotid atherosclerosis can be prevented by lifestyle/dietary changes and pharmacological treatments if identified early by monitoring carotid plaque changes. Plaque composition information from magnetic resonance (MR) carotid images and dynamic characteristics information from 3D ultrasound (US) are necessary for developing and validating US imaging tools to identify vulnerable carotid plaques. Combining these images requires nonrigid registration to correct the non-linear miss-alignments caused by relative twisting and bending in the neck due to different head positions during the two image acquisitions sessions. The high degree of freedom and large number of parameters associated with existing nonrigid image registration methods causes several problems including unnatural plaque morphology alteration, computational complexity, and low reliability. Our approach was to model the normal movement of the neck using a "twisting and bending model" with only six parameters for nonrigid registration. We evaluated our registration technique using intra-subject in-vivo 3D US and 3D MR carotid images acquired on the same day. We calculated the Mean Registration Error (MRE) between the segmented vessel surfaces in the target image and the registered image using a distance-based error metric after applying our "twisting bending model" based nonrigid registration algorithm. We achieved an average registration error of 1.33+/-0.41mm using our nonrigid registration technique. Visual inspection of segmented vessel surfaces also showed a substantial improvement of alignment with our non-rigid registration technique.
MUTUAL INFORMATION BASED 3D NON-RIGID REGISTRATION OF CT/MR ABDOMEN IMAGES
无
2001-01-01
A mutual information based 3D non-rigid registration approach was proposed for the registration of deformable CT/MR body abdomen images. The Parzen Windows Density Estimation (PWDE) method is adopted to calculate the mutual information between the two modals of CT and MRI abdomen images. By maximizing MI between the CT and MR volume images, the overlapping part of them reaches the biggest, which means that the two body images of CT and MR matches best to each other. Visible Human Project (VHP) Male abdomen CT and MRI Data are used as experimental data sets. The experimental results indicate that this approach of non-rigid 3D registration of CT/MR body abdominal images can be achieved effectively and automatically, without any prior processing procedures such as segmentation and feature extraction, but has a main drawback of very long computation time. Key words: medical image registration; multi-modality; mutual information; non-rigid; Parzen window density estimation
Nonrigid Image Registration for Head and Neck Cancer Radiotherapy Treatment Planning With PET/CT
Purpose: Head and neck radiotherapy planning with positron emission tomography/computed tomography (PET/CT) requires the images to be reliably registered with treatment planning CT. Acquiring PET/CT in treatment position is problematic, and in practice for some patients it may be beneficial to use diagnostic PET/CT for radiotherapy planning. Therefore, the aim of this study was first to quantify the image registration accuracy of PET/CT to radiotherapy CT and, second, to assess whether PET/CT acquired in diagnostic position can be registered to planning CT. Methods and Materials: Positron emission tomography/CT acquired in diagnostic and treatment position for five patients with head and neck cancer was registered to radiotherapy planning CT using both rigid and nonrigid image registration. The root mean squared error for each method was calculated from a set of anatomic landmarks marked by four independent observers. Results: Nonrigid and rigid registration errors for treatment position PET/CT to planning CT were 2.77 ± 0.80 mm and 4.96 ± 2.38 mm, respectively, p = 0.001. Applying the nonrigid registration to diagnostic position PET/CT produced a more accurate match to the planning CT than rigid registration of treatment position PET/CT (3.20 ± 1.22 mm and 4.96 ± 2.38 mm, respectively, p = 0.012). Conclusions: Nonrigid registration provides a more accurate registration of head and neck PET/CT to treatment planning CT than rigid registration. In addition, nonrigid registration of PET/CT acquired with patients in a standardized, diagnostic position can provide images registered to planning CT with greater accuracy than a rigid registration of PET/CT images acquired in treatment position. This may allow greater flexibility in the timing of PET/CT for head and neck cancer patients due to undergo radiotherapy
Free-form deformation based non-rigid registration on breast cancer MR imaging
Zhang, Liangbin; Suo, Shiteng; Lu, Xuesong; Li, Yuehua; Chen, Li; Zhang, Su
2013-07-01
High-Intensity Focused Ultrasound treatment combined with magnetic resonance technology (MRI-guided HIFU, MRgHIFU) can protect the thermal ablation without harming the surrounding tissue by using MRI for target positioning, where image registration plays an important role in the implementation of precise treatment. In this paper, we apply three-dimension free-form deformation non-rigid registration on treatment plan amendments and tracking of breast cancer. Free-form deformation based and demons based non-rigid registration are respectively employed on breast cancer MR imaging required at different times before and after for comparison. The results of the experiments show that the registration performed on the breast tumor image data with slight and larger deformation is effective, and the mutual information of the ROI increased from 1.49 before registration to 1.53.
Non-rigid registration of volumetric images using ranked order statistics
Tennakoon, Ruwan; Bab-Hadiashar, Alireza; Cao, Zhenwei; de Bruijne, Marleen
2014-01-01
Non-rigid image registration techniques using intensity based similarity measures are widely used in medical imaging applications. Due to high computational complexities of these techniques, particularly for volumetric images, finding appropriate registration methods to both reduce the computation...... the important regions for registration and use an appropriate sampling scheme to target those areas and reduce the registration computation time. A unique advantage of the proposed method is its ability to identify the point of diminishing returns and stop the registration process. Our experiments on...
Robust non-rigid point set registration using student's-t mixture model.
Zhiyong Zhou
Full Text Available The Student's-t mixture model, which is heavily tailed and more robust than the Gaussian mixture model, has recently received great attention on image processing. In this paper, we propose a robust non-rigid point set registration algorithm using the Student's-t mixture model. Specifically, first, we consider the alignment of two point sets as a probability density estimation problem and treat one point set as Student's-t mixture model centroids. Then, we fit the Student's-t mixture model centroids to the other point set which is treated as data. Finally, we get the closed-form solutions of registration parameters, leading to a computationally efficient registration algorithm. The proposed algorithm is especially effective for addressing the non-rigid point set registration problem when significant amounts of noise and outliers are present. Moreover, less registration parameters have to be set manually for our algorithm compared to the popular coherent points drift (CPD algorithm. We have compared our algorithm with other state-of-the-art registration algorithms on both 2D and 3D data with noise and outliers, where our non-rigid registration algorithm showed accurate results and outperformed the other algorithms.
Onofrey, John A.; Staib, Lawrence H.; Xenophon Papademetris
2016-01-01
This paper describes a framework for learning a statistical model of non-rigid deformations induced by interventional procedures. We make use of this learned model to perform constrained non-rigid registration of pre-procedural and post-procedural imaging. We demonstrate results applying this framework to non-rigidly register post-surgical computed tomography (CT) brain images to pre-surgical magnetic resonance images (MRIs) of epilepsy patients who had intra-cranial electroencephalography el...
Evaluation and validation methods for intersubject nonrigid 3D image registration of the human brain
Guo, Ting; Starreveld, Yves P.; Peters, Terry M.
2005-04-01
This work presents methodologies for assessing the accuracy of non-rigid intersubject registration algorithms from both qualitative and quantitative perspectives. The first method was based on a set of 43 anatomical landmarks. MRI brain images of 12 subjects were non-rigidly registered to the standard MRI dataset. The "gold-standard" coordinates of the 43 landmarks in the target were estimated by averaging their coordinates after 6 tagging sessions. The Euclidean distance between each landmark of a subject after warping to the reference space and the homologous "gold-standard" landmark on the reference image was considered as the registration error. Another method based on visual inspection software displaying the spatial change of colour-coded spheres, before and after warping, was also developed to evaluate the performance of the non-rigid warping algorithms within the homogeneous regions in the deep-brain. Our methods were exemplified by assessing and comparing the accuracy of two intersubject non-rigid registration approaches, AtamaiWarp and ANIMAL algorithms. From the first method, the average registration error was 1.04mm +/- 0.65mm for AtamaiWarp, and 1.59mm +/- 1.47mm for ANIMAL. With maximum registration errors of 2.78mm and 3.90mm respectively, AtamaiWarp and ANIMAL located 58% and 35% landmarks respectively with registration errors less than 1mm. A paired t-test showed that the differences in registration error between AtamaiWarp and ANIMAL were significant (P ANIMAL, also provides more accurate results. From the second method, both algorithms treated the interior of homogeneous regions in an appropriate manner.
Xiao, Di; Zahra, David; Bourgeat, Pierrick; Berghofer, Paula; Acosta Tamayo, Oscar; Wimberley, Catriona; Gregoire, Marie C.; Salvado, Olivier
2011-03-01
Small animal image registration is challenging because of its joint structure, and posture and position difference in each acquisition without a standard scan protocol. In this paper, we face the issue of mouse whole-body skeleton registration from CT images. A novel method is developed for analyzing mouse hind-limb and fore-limb postures based on geodesic path descriptor and then registering the major skeletons and fore limb skeletons initially by thin-plate spline (TPS) transform based on the obtained geodesic paths and their enhanced correspondence fields. A target landmark correction method is proposed for improving the registration accuracy of the improved 3D shape context non-rigid registration method we previously proposed. A novel non-rigid registration framework, combining the skeleton posture analysis, geodesic path based initial alignment and 3D shape context model, is proposed for mouse whole-body skeleton registration. The performance of the proposed methods and framework was tested on 12 pairs of mouse whole-body skeletons. The experimental results demonstrated the flexibility, stability and accuracy of the proposed framework for automatic mouse whole body skeleton registration.
Deng, Zhipeng; Lei, Lin; Zhou, Shilin
2015-10-01
Automatic image registration is a vital yet challenging task, particularly for non-rigid deformation images which are more complicated and common in remote sensing images, such as distorted UAV (unmanned aerial vehicle) images or scanning imaging images caused by flutter. Traditional non-rigid image registration methods are based on the correctly matched corresponding landmarks, which usually needs artificial markers. It is a rather challenging task to locate the accurate position of the points and get accurate homonymy point sets. In this paper, we proposed an automatic non-rigid image registration algorithm which mainly consists of three steps: To begin with, we introduce an automatic feature point extraction method based on non-linear scale space and uniform distribution strategy to extract the points which are uniform distributed along the edge of the image. Next, we propose a hybrid point matching algorithm using DaLI (Deformation and Light Invariant) descriptor and local affine invariant geometric constraint based on triangulation which is constructed by K-nearest neighbor algorithm. Based on the accurate homonymy point sets, the two images are registrated by the model of TPS (Thin Plate Spline). Our method is demonstrated by three deliberately designed experiments. The first two experiments are designed to evaluate the distribution of point set and the correctly matching rate on synthetic data and real data respectively. The last experiment is designed on the non-rigid deformation remote sensing images and the three experimental results demonstrate the accuracy, robustness, and efficiency of the proposed algorithm compared with other traditional methods.
3D nonrigid medical image registration using a new information theoretic measure
Li, Bicao; Yang, Guanyu; Coatrieux, Jean Louis; Li, Baosheng; Shu, Huazhong
2015-11-01
This work presents a novel method for the nonrigid registration of medical images based on the Arimoto entropy, a generalization of the Shannon entropy. The proposed method employed the Jensen-Arimoto divergence measure as a similarity metric to measure the statistical dependence between medical images. Free-form deformations were adopted as the transformation model and the Parzen window estimation was applied to compute the probability distributions. A penalty term is incorporated into the objective function to smooth the nonrigid transformation. The goal of registration is to optimize an objective function consisting of a dissimilarity term and a penalty term, which would be minimal when two deformed images are perfectly aligned using the limited memory BFGS optimization method, and thus to get the optimal geometric transformation. To validate the performance of the proposed method, experiments on both simulated 3D brain MR images and real 3D thoracic CT data sets were designed and performed on the open source elastix package. For the simulated experiments, the registration errors of 3D brain MR images with various magnitudes of known deformations and different levels of noise were measured. For the real data tests, four data sets of 4D thoracic CT from four patients were selected to assess the registration performance of the method, including ten 3D CT images for each 4D CT data covering an entire respiration cycle. These results were compared with the normalized cross correlation and the mutual information methods and show a slight but true improvement in registration accuracy.
3D nonrigid medical image registration using a new information theoretic measure
This work presents a novel method for the nonrigid registration of medical images based on the Arimoto entropy, a generalization of the Shannon entropy. The proposed method employed the Jensen–Arimoto divergence measure as a similarity metric to measure the statistical dependence between medical images. Free-form deformations were adopted as the transformation model and the Parzen window estimation was applied to compute the probability distributions. A penalty term is incorporated into the objective function to smooth the nonrigid transformation. The goal of registration is to optimize an objective function consisting of a dissimilarity term and a penalty term, which would be minimal when two deformed images are perfectly aligned using the limited memory BFGS optimization method, and thus to get the optimal geometric transformation. To validate the performance of the proposed method, experiments on both simulated 3D brain MR images and real 3D thoracic CT data sets were designed and performed on the open source elastix package. For the simulated experiments, the registration errors of 3D brain MR images with various magnitudes of known deformations and different levels of noise were measured. For the real data tests, four data sets of 4D thoracic CT from four patients were selected to assess the registration performance of the method, including ten 3D CT images for each 4D CT data covering an entire respiration cycle. These results were compared with the normalized cross correlation and the mutual information methods and show a slight but true improvement in registration accuracy. (paper)
B-spline techniques for volatility modeling
Corlay, Sylvain
2013-01-01
This paper is devoted to the application of B-splines to volatility modeling, specifically the calibration of the leverage function in stochastic local volatility models and the parameterization of an arbitrage-free implied volatility surface calibrated to sparse option data. We use an extension of classical B-splines obtained by including basis functions with infinite support. We first come back to the application of shape-constrained B-splines to the estimation of conditional expectations, ...
Direct estimation of nonrigid registrations with image-based self-occlusion reasoning.
Gay-Bellile, Vincent; Bartoli, Adrien; Sayd, Patrick
2010-01-01
The registration problem for images of a deforming surface has been well studied. External occlusions are usually well handled. In 2D image-based registration, self-occlusions are more challenging. Consequently, the surface is usually assumed to be only slightly self-occluding. This paper is about image-based nonrigid registration with self-occlusion reasoning. A specific framework explicitly modeling self-occlusions is proposed. It is combined with an intensity-based, "direct" data term for registration. Self-occlusions are detected as shrinkage areas in the 2D warp. Experimental results on several challenging data sets show that our approach successfully registers images with self-occlusions while effectively detecting the self-occluded regions. PMID:19926901
Kjer, Hans Martin; Fagertun, Jens; Vera, Sergio;
2016-01-01
Better understanding of the anatomical variability of the human cochlear is important for the design and function of Cochlear Implants. Proper non-rigid alignment of high-resolution cochlear μCT data is a challenge for the typical cubic B-spline registration model. In this paper we study one way ...
A kidney deformation model for use in non-rigid registration during image-guided surgery
Ong, Rowena E.; Herrell, S. Duke, III; Miga, Michael I.; Galloway, Robert L., Jr.
2008-03-01
In order to facilitate the removal of tumors during partial nephrectomies, an image-guided surgery system may be useful. This system would require a registration of the physical kidney to a pre-operative image volume; however, it is unclear whether a rigid registration would be sufficient. One possible source of non-rigid deformation is the clamping of the renal artery during surgery and the subsequent loss of pressure as the kidney is punctured and blood loss occurs. To explore this issue, a model of kidney deformation due to loss of perfusion and pressure was developed based on Biot's consolidation model. The model was tested on two resected porcine kidneys in which the renal artery and vein were clamped. CT image volumes of the kidney were obtained before and after the deformation caused unclamping, and fiducial markers embedded on the kidney surface allowed the deformation to be tracked. The accuracy of the kidney model was accessed by calculating the model error at the fiducial locations and using image similarity measures. Preliminary results indicate that the model may be useful in a non-rigid registration scheme; however, further refinements to the model may be necessary to better simulate the deformation due to loss of perfusion and pressure.
MUTUAL INFORMATION BASED 3D NON-RIGID REGISTRATION OF CT/MR ABDOMEN IMAGES
HU; Hai-bo(
2001-01-01
［1］Maintz J B, Viergever M A. A survey of medical image registration[J]. Medical Image Analysis, 1998, 3(1):1～37.［2］Collignon A. Automated multi-modality image registration based on information theory[J]. Computational Imaging and vision, 1995, 3:263～274.［3］Eberl S, Braun M. Intra-and inter-modality registration of functional and anatomical clinical images[A]. Pham B, et al. eds. New Approaches in Medical Image Analysis, SPIE 3747[C].[s.l.]:[s.n.], 1999. 102～114.［4］Lau Y H, Braun M, Hutton B F. Non-rigid 3D image registration using regionally constrained matching and the correlation ratio[A]. Pernus F, et al.eds. Biomedical Image Registration, Proc Int Workshop[C]. Bled, Slovenia, 1999. 137～148.［5］Wells Ⅲ W M, Viola P, Atsumi H, et al. Multi-modal volume registration by maximization of mutual information[J]. Medical Image Analysis, 1996, 1(1):35～51.［6］Feldmar J, Ayache N. Rigid, affine and locally affine registration of free-form surfaces[J]. Int J of Computer Vision, 1996, 23(3):97～104.
Purpose: In oxygen-enhanced magnetic resonance imaging of the lung (O2-MRI), motion artifacts related to breathing hamper the quality of the parametric O2-maps. In this study, fully automatic non-rigid image registration was assessed as a post-processing method to improve the quality of O2-MRI. Materials and methods: Twenty healthy volunteers were investigated on a 1.5 T MR system. O2-MRI was obtained in four coronal sections using an IR-HASTE sequence with TE/TI of 12/1200 ms. Each section was repeatedly imaged during oxygen and room-air ventilation. Spatial differences among the images were corrected by fully automatic non-rigid registration. Signal variability, relative enhancement ratio between oxygen and room air images, and spatial heterogeneity of lung enhancement were assessed before and after image registration. Results: Motion artifacts were corrected in 5–10 s. Non-rigid registration reduced signal variability of the source images and heterogeneity of the O2-maps by 1.1 ± 0.2% and 11.2 ± 2.9%, respectively (p < 0.0001). Registration did not influence O2 relative enhancement ratio (p = 0.06). Conclusion: Fully automatic non-rigid image registration improves the quality of multislice oxygen-enhanced MRI of the lung
Non-rigid registration and KLT filter to improve SNR and CNR in GRE-EPI myocardial perfusion imaging
Mihai, Georgeta; Ding, Yu; Xue, Hui; Chung, Yiu-Cho; Rajagopalan, Sanjay; Guehring, Jens; Simonetti, Orlando P.
2012-01-01
The purpose of the study was to evaluate the effect of motion compensation by non-rigid registration combined with the Karhunen-Loeve Transform (KLT) filter on the signal to noise (SNR) and contrast-to-noise ratio (CNR) of hybrid gradient-echo echoplanar (GRE-EPI) first-pass myocardial perfusion imaging. Twenty one consecutive first-pass adenosine stress perfusion MR data sets interpreted positive for ischemia or infarction were processed by non-rigid Registration followed by KLT filtering. S...
Fuzzy B-Spline Surface Modeling
Rozaimi Zakaria
2014-01-01
Full Text Available This paper discusses the construction of a fuzzy B-spline surface model. The construction of this model is based on fuzzy set theory which is based on fuzzy number and fuzzy relation concepts. The proposed theories and concepts define the uncertainty data sets which represent fuzzy data/control points allowing the uncertainties data points modeling which can be visualized and analyzed. The fuzzification and defuzzification processes were also defined in detail in order to obtain the fuzzy B-spline surface crisp model. Final section shows an application of fuzzy B-spline surface modeling for terrain modeling which shows its usability in handling uncertain data.
A deformation model for non-rigid registration of the kidney
Ong, Rowena E.; Glisson, Courtenay L.; Herrell, S. Duke; Miga, Michael I.; Galloway, Robert
2009-02-01
The development of an image-guided renal surgery system may aid tumor resection during partial nephrectomies. This system would require the registration of pre-operative kidney CT or MR scans to the physical kidney; however, the amount of non-rigid deformation occurring during surgery and whether it can be corrected for in an image-guided system is unknown. One possible source of non-rigid deformation is a change in pressure within the kidney: during surgery, clamping of the renal artery and vein results in a loss of perfusion, such that the subsequent cutting of the kidney and fluid outflow may cause a decrease in intrarenal pressure. In this work, we attempt to characterize the deformation due to cutting of the kidney and subsequent changes in intrarenal pressure. To accomplish this, we perfused a resected porcine kidney at a physiologically realistic pressure, clamped the renal vessels, and cut the kidney using a tracked scalpel. The resulting deformation was tracked in a CT scanner using 15-20 glass bead fiducials attached to the kidney surface. A modified form of Biot's consolidation model was used to simulate the deformation, and the accuracy was assessed by calculating the target registration error and image similarity.
Technical improvements in planning and dose delivery and in verification of patient positioning have substantially widened the therapeutic window for radiation treatment of cancer. However, changes in patient anatomy during the treatment limit the exploitation of these new techniques. To further improve radiation treatments, anatomical changes need to be modeled and accounted for. Nonrigid registration can be used for this purpose. This article describes the design, the implementation, and the validation of a new framework for nonrigid registration for radiotherapy applications. The core of this framework is an improved version of the thin plate spline robust point matching (TPS-RPM) algorithm. The TPS-RPM algorithm estimates a global correspondence and a transformation between the points that represent organs of interest belonging to two image sets. However, the algorithm does not allow for the inclusion of prior knowledge on the correspondence of subset of points, and therefore, it can lead to inconsistent anatomical solutions. In this article TPS-RPM was improved by employing a novel correspondence filter that supports simultaneous registration of multiple structures. The improved method allows for coherent organ registration and for the inclusion of user-defined landmarks, lines, and surfaces inside and outside of structures of interest. A procedure to generate control points from segmented organs is described. The framework parameters r and λ, which control the number of points and the nonrigidness of the transformation, respectively, were optimized for three sites with different degrees of deformation (head and neck, prostate, and cervix) using two cases per site. For the head and neck cases, the salivary glands were manually contoured on CT scans, for the prostate cases the prostate and the vesicles, and for the cervix cases the cervix uterus, the bladder, and the rectum. The transformation error obtained using the best set of parameters was below 1 mm for
Non-Rigid Registration of Liver CT Images for CT-Guided Ablation of Liver Tumors.
Luu, Ha Manh; Klink, Camiel; Niessen, Wiro; Moelker, Adriaan; Walsum, Theo van
2016-01-01
CT-guided percutaneous ablation for liver cancer treatment is a relevant technique for patients not eligible for surgery and with tumors that are inconspicuous on US imaging. The lack of real-time imaging and the use of a limited amount of CT contrast agent make targeting the tumor with the needle challenging. In this study, we evaluate a registration framework that allows the integration of diagnostic pre-operative contrast enhanced CT images and intra-operative non-contrast enhanced CT images to improve image guidance in the intervention. The liver and tumor are segmented in the pre-operative contrast enhanced CT images. Next, the contrast enhanced image is registered to the intra-operative CT images in a two-stage approach. First, the contrast-enhanced diagnostic image is non-rigidly registered to a non-contrast enhanced image that is conventionally acquired at the start of the intervention. In case the initial registration is not sufficiently accurate, a refinement step is applied using non-rigid registration method with a local rigidity term. In the second stage, the intra-operative CT-images that are used to check the needle position, which often consist of only a few slices, are registered rigidly to the intra-operative image that was acquired at the start of the intervention. Subsequently, the diagnostic image is registered to the current intra-operative image, using both transformations, this allows the visualization of the tumor region extracted from pre-operative data in the intra-operative CT images containing needle. The method is evaluated on imaging data of 19 patients at the Erasmus MC. Quantitative evaluation is performed using the Dice metric, mean surface distance of the liver border and corresponding landmarks in the diagnostic and the intra-operative images. The registration of the diagnostic CT image to the initial intra-operative CT image did not require a refinement step in 13 cases. For those cases, the resulting registration had a Dice
Comparing nonrigid registration techniques for motion corrected MR prostate diffusion imaging
Purpose: T2-weighted magnetic resonance imaging (MRI) is commonly used for anatomical visualization in the pelvis area, such as the prostate, with high soft-tissue contrast. MRI can also provide functional information such as diffusion-weighted imaging (DWI) which depicts the molecular diffusion processes in biological tissues. The combination of anatomical and functional imaging techniques is widely used in oncology, e.g., for prostate cancer diagnosis and staging. However, acquisition-specific distortions as well as physiological motion lead to misalignments between T2 and DWI and consequently to a reduced diagnostic value. Image registration algorithms are commonly employed to correct for such misalignment. Methods: The authors compare the performance of five state-of-the-art nonrigid image registration techniques for accurate image fusion of DWI with T2. Results: Image data of 20 prostate patients with cancerous lesions or cysts were acquired. All registration algorithms were validated using intensity-based as well as landmark-based techniques. Conclusions: The authors’ results show that the “fast elastic image registration” provides most accurate results with a target registration error of 1.07 ± 0.41 mm at minimum execution times of 11 ± 1 s
Comparing nonrigid registration techniques for motion corrected MR prostate diffusion imaging
Buerger, C., E-mail: christian.buerger@philips.com; Sénégas, J.; Kabus, S.; Carolus, H.; Schulz, H.; Renisch, S. [Philips Research Hamburg, Hamburg 22335 (Germany); Agarwal, H. [Philips Research North America, Briarcliff Manor, New York 10510 and Molecular Imaging Program, NCI, National Institute of Health, Bethesda, Maryland 20892 (United States); Turkbey, B.; Choyke, P. L. [Molecular Imaging Program, NCI, National Institute of Health, Bethesda, Maryland 20892 (United States)
2015-01-15
Purpose: T{sub 2}-weighted magnetic resonance imaging (MRI) is commonly used for anatomical visualization in the pelvis area, such as the prostate, with high soft-tissue contrast. MRI can also provide functional information such as diffusion-weighted imaging (DWI) which depicts the molecular diffusion processes in biological tissues. The combination of anatomical and functional imaging techniques is widely used in oncology, e.g., for prostate cancer diagnosis and staging. However, acquisition-specific distortions as well as physiological motion lead to misalignments between T{sub 2} and DWI and consequently to a reduced diagnostic value. Image registration algorithms are commonly employed to correct for such misalignment. Methods: The authors compare the performance of five state-of-the-art nonrigid image registration techniques for accurate image fusion of DWI with T{sub 2}. Results: Image data of 20 prostate patients with cancerous lesions or cysts were acquired. All registration algorithms were validated using intensity-based as well as landmark-based techniques. Conclusions: The authors’ results show that the “fast elastic image registration” provides most accurate results with a target registration error of 1.07 ± 0.41 mm at minimum execution times of 11 ± 1 s.
van Velden, Floris H. P.; Ida A Nissen; Wendy Hayes; Velasquez, Linda M; Hoekstra, Otto S.; Ronald Boellaard
2014-01-01
OBJECTIVES: Reusing baseline volumes of interest (VOI) by applying non-rigid and to some extent (local) rigid image registration showed good test-retest variability similar to delineating VOI on both scans individually. The aim of the present study was to compare response assessments and classifications based on various types of image registration with those based on (semi)-automatic tumour delineation. METHODS: Baseline (n = 13), early (n = 12) and late (n = 9) response (after one and three ...
Non-rigid registration of CT-MRI images for visualizing vertebra and nerves
The operation of the orthopedics department under the endoscope is an operation that approaches the diseased part to cut only a part of the vertebra unlike a fixed art of the spine with a past fastening device, and eases the pressure of nerve. It aims at the preoperative planning. Lumbar vertebra and vertebrae cervicales taken picture of by CT and Nerve that can be taken picture by MRI visible making simultaneously is requested. But taking picture of CT and MRI are difficult to keep patient's posture to be the same. Problem that gap is caused by solid registration. This present study proposes the method of the non-rigid registration that paid attention to shape of the body in both of CT and MRI taken picture. And generate mesh by the characteristic point based on the curvature in outside of the body. We defined evaluation function the distance between the mesh and shape preservation. And registration is achieved by solving the least square problem. This presentation reports proposal technique about the medical image data. (author)
Validation of non-rigid point-set registration methods using a porcine bladder pelvic phantom
Zakariaee, Roja; Hamarneh, Ghassan; Brown, Colin J.; Spadinger, Ingrid
2016-01-01
The problem of accurate dose accumulation in fractionated radiotherapy treatment for highly deformable organs, such as bladder, has garnered increasing interest over the past few years. However, more research is required in order to find a robust and efficient solution and to increase the accuracy over the current methods. The purpose of this study was to evaluate the feasibility and accuracy of utilizing non-rigid (affine or deformable) point-set registration in accumulating dose in bladder of different sizes and shapes. A pelvic phantom was built to house an ex vivo porcine bladder with fiducial landmarks adhered onto its surface. Four different volume fillings of the bladder were used (90, 180, 360 and 480 cc). The performance of MATLAB implementations of five different methods were compared, in aligning the bladder contour point-sets. The approaches evaluated were coherent point drift (CPD), gaussian mixture model, shape context, thin-plate spline robust point matching (TPS-RPM) and finite iterative closest point (ICP-finite). The evaluation metrics included registration runtime, target registration error (TRE), root-mean-square error (RMS) and Hausdorff distance (HD). The reference (source) dataset was alternated through all four points-sets, in order to study the effect of reference volume on the registration outcomes. While all deformable algorithms provided reasonable registration results, CPD provided the best TRE values (6.4 mm), and TPS-RPM yielded the best mean RMS and HD values (1.4 and 6.8 mm, respectively). ICP-finite was the fastest technique and TPS-RPM, the slowest.
Validation of non-rigid point-set registration methods using a porcine bladder pelvic phantom.
Zakariaee, Roja; Hamarneh, Ghassan; Brown, Colin J; Spadinger, Ingrid
2016-01-21
The problem of accurate dose accumulation in fractionated radiotherapy treatment for highly deformable organs, such as bladder, has garnered increasing interest over the past few years. However, more research is required in order to find a robust and efficient solution and to increase the accuracy over the current methods. The purpose of this study was to evaluate the feasibility and accuracy of utilizing non-rigid (affine or deformable) point-set registration in accumulating dose in bladder of different sizes and shapes. A pelvic phantom was built to house an ex vivo porcine bladder with fiducial landmarks adhered onto its surface. Four different volume fillings of the bladder were used (90, 180, 360 and 480 cc). The performance of MATLAB implementations of five different methods were compared, in aligning the bladder contour point-sets. The approaches evaluated were coherent point drift (CPD), gaussian mixture model, shape context, thin-plate spline robust point matching (TPS-RPM) and finite iterative closest point (ICP-finite). The evaluation metrics included registration runtime, target registration error (TRE), root-mean-square error (RMS) and Hausdorff distance (HD). The reference (source) dataset was alternated through all four points-sets, in order to study the effect of reference volume on the registration outcomes. While all deformable algorithms provided reasonable registration results, CPD provided the best TRE values (6.4 mm), and TPS-RPM yielded the best mean RMS and HD values (1.4 and 6.8 mm, respectively). ICP-finite was the fastest technique and TPS-RPM, the slowest. PMID:26740511
Lung segmentation in chest radiographs using anatomical atlases with nonrigid registration.
Candemir, Sema; Jaeger, Stefan; Palaniappan, Kannappan; Musco, Jonathan P; Singh, Rahul K; Zhiyun Xue; Karargyris, Alexandros; Antani, Sameer; Thoma, George; McDonald, Clement J
2014-02-01
The National Library of Medicine (NLM) is developing a digital chest X-ray (CXR) screening system for deployment in resource constrained communities and developing countries worldwide with a focus on early detection of tuberculosis. A critical component in the computer-aided diagnosis of digital CXRs is the automatic detection of the lung regions. In this paper, we present a nonrigid registration-driven robust lung segmentation method using image retrieval-based patient specific adaptive lung models that detects lung boundaries, surpassing state-of-the-art performance. The method consists of three main stages: 1) a content-based image retrieval approach for identifying training images (with masks) most similar to the patient CXR using a partial Radon transform and Bhattacharyya shape similarity measure, 2) creating the initial patient-specific anatomical model of lung shape using SIFT-flow for deformable registration of training masks to the patient CXR, and 3) extracting refined lung boundaries using a graph cuts optimization approach with a customized energy function. Our average accuracy of 95.4% on the public JSRT database is the highest among published results. A similar degree of accuracy of 94.1% and 91.7% on two new CXR datasets from Montgomery County, MD, USA, and India, respectively, demonstrates the robustness of our lung segmentation approach. PMID:24239990
Lei, Peng; Dandekar, Omkar; Widlus, David; Shekhar, Raj
2009-01-01
This study evaluates the accuracy of augmenting initial intraprocedural computed tomography (CT) during radiofrequency ablation (RFA) of hepatic metastases with preprocedural positron emission tomography (PET) through a hardware-accelerated implementation of an automatic nonrigid PET–CT registration algorithm. The feasibility of augmenting intraprocedural CT with preprocedural PET to improve localization of CT-invisible but PET-positive tumors with images from actual RFA was explored. Preproc...
Quadrotor system identification using the multivariate multiplex b-spline
Visser, T.; De Visser, C.C.; Van Kampen, E.J.
2015-01-01
A novel method for aircraft system identification is presented that is based on a new multivariate spline type; the multivariate multiplex B-spline. The multivariate multiplex B-spline is a generalization of the recently introduced tensor-simplex B-spline. Multivariate multiplex splines obtain simil
DIRECT MANIPULATION OF B-SPLINE SURFACES
Wang Zhiguo; Zhou Laishui; Wang Xiaoping
2005-01-01
Engineering design and geometric modeling often require the ability to modify the shape of parametric curves and surfaces so that their shape satisfy some given geometric constraints, including point, normal vector, curve and surface. Two approaches are presented to directly manipulate the shape of B-spline surface. The former is based on the least-square, whereas the latter is based on minimizing the bending energy of surface. For each method, since unified and explicit formulae are derived to compute new control points of modified surface, these methods are simple, fast and applicable for CAD systems. Algebraic technique is used to simplify the computation of B-spline composition and multiplication. Comparisons and examples are also given.
Mihai, Georgeta; Ding, Yu; Xue, Hui; Chung, Yiu-Cho; Rajagopalan, Sanjay; Guehring, Jens; Simonetti, Orlando P
2012-12-01
The purpose of the study was to evaluate the effect of motion compensation by non-rigid registration combined with the Karhunen-Loeve Transform (KLT) filter on the signal to noise (SNR) and contrast-to-noise ratio (CNR) of hybrid gradient-echo echoplanar (GRE-EPI) first-pass myocardial perfusion imaging. Twenty one consecutive first-pass adenosine stress perfusion MR data sets interpreted positive for ischemia or infarction were processed by non-rigid Registration followed by KLT filtering. SNR and CNR were measured in abnormal and normal myocardium in unfiltered and KLT filtered images following non-rigid registration to compensate for respiratory and other motions. Image artifacts introduced by filtering in registered and nonregistered images were evaluated by two observers. There was a statistically significant increase in both SNR and CNR between normal and abnormal myocardium with KLT filtering (mean SNR increased by 62.18% ± 21.05% and mean CNR increased by 58.84% ± 18.06%; p = 0.01). Motion correction prior to KLT filtering reduced significantly the occurrence of filter induced artifacts (KLT only-artifacts in 42 out of 55 image series vs. registered plus KLT-artifacts in 3 out of 55 image series). In conclusion the combination of non- rigid registration and KLT filtering was shown to increase the SNR and CNR of GRE-EPI perfusion images. Subjective evaluation of image artifacts revealed that prior motion compensation significantly reduced the artifacts introduced by the KLT filtering process. PMID:23936584
Purpose: This study introduces a probabilistic nonrigid registration method for use in image-guided prostate brachytherapy. Intraoperative imaging for prostate procedures, usually transrectal ultrasound (TRUS), is typically inferior to diagnostic-quality imaging of the pelvis such as endorectal magnetic resonance imaging (MRI). MR images contain superior detail of the prostate boundaries and provide substructure features not otherwise visible. Previous efforts to register diagnostic prostate images with the intraoperative coordinate system have been deterministic and did not offer a measure of the registration uncertainty. The authors developed a Bayesian registration method to estimate the posterior distribution on deformations and provide a case-specific measure of the associated registration uncertainty. Methods: The authors adapted a biomechanical-based probabilistic nonrigid method to register diagnostic to intraoperative images by aligning a physician's segmentations of the prostate in the two images. The posterior distribution was characterized with a Markov Chain Monte Carlo method; the maximum a posteriori deformation and the associated uncertainty were estimated from the collection of deformation samples drawn from the posterior distribution. The authors validated the registration method using a dataset created from ten patients with MRI-guided prostate biopsies who had both diagnostic and intraprocedural 3 Tesla MRI scans. The accuracy and precision of the estimated posterior distribution on deformations were evaluated from two predictive distance distributions: between the deformed central zone-peripheral zone (CZ-PZ) interface and the physician-labeled interface, and based on physician-defined landmarks. Geometric margins on the registration of the prostate's peripheral zone were determined from the posterior predictive distance to the CZ-PZ interface separately for the base, mid-gland, and apical regions of the prostate. Results: The authors observed
Scalable low-complexity B-spline discretewavelet transform architecture
Martina, Maurizio; Masera, Guido; Piccinini, Gianluca
2010-01-01
A scalable discrete wavelet transform architecture based on the B-spline factorisation is presented. In particular, it is shown that several wavelet filters of practical interest have a common structure in the distributed part of their B-spline factorisation. This common structure is effectively exploited to achieve scalability and to save multipliers compared with a direct polyphase B-spline implementation. Since the proposed solution is more robust to coefficient quantisation than direct po...
Faired MISO B-Spline Fuzzy Systems and Its Applications
Tan Yanhua
2013-01-01
Full Text Available We construct two classes of faired MISO B-spline fuzzy systems using the fairing method in computer-aided geometric design (CAGD for reducing adverse effects of the inexact data. Towards this goal, we generalize the faring method to high-dimension cases so that the faring method only for SISO and DISO B-spline fuzzy systems is extended to fair the MISO ones. Then the problem to construct a faired MISO B-spline fuzzy systems is transformed into solving an optimization problem with a strictly convex quadratic objective function and the unique optimal solution vector is taken as linear combination coefficients of the basis functions for a certain B-spline fuzzy system to obtain a faired MISO B-spline fuzzy system. Furthermore, we design variable universe adaptive fuzzy controllers by B-spline fuzzy systems and faired B-spline fuzzy systems to stabilize the double inverted pendulum. The simulation results show that the controllers by faired B-spline fuzzy systems perform better than those by B-spline fuzzy systems, especially when the data for fuzzy systems are inexact.
Berkels, Benjamin; Deserno, Thomas; Ehrlich, Eva E.; Fritz, Ulrike B.; Sirazitdinova, Ekaterina; Tatano, Rosalia
2016-03-01
Quantitative light-induced fluorescence (QLF) is widely used to assess the damage of a tooth due to decalcification. In digital photographs, decalcification appears as white spot lesions, i.e. white spots on the tooth surface. We propose a novel multimodal registration approach for the matching of digital photographs and QLF images of decalcified teeth. The registration is based on the idea of contour-to-pixel matching. Here, the curve, which represents the shape of the tooth, is extracted from the QLF image using a contour segmentation by binarization and morphological processing. This curve is aligned to the photo with a non-rigid variational registration approach. Thus, the registration problem is formulated as minimization problem with an objective function that consists of a data term and a regularizer for the deformation. To construct the data term, the photo is pointwise classified into tooth and non-tooth regions. Then, the signed distance function of the tooth region allows to measure the mismatch between curve and photo. As regularizer a higher order, linear elastic prior is used. The resulting minimization problem is solved numerically using bilinear Finite Elements for the spatial discretization and the Gauss-Newton algorithm. The evaluation is based on 150 image pairs, where an average of 5 teeth have been captured from 32 subjects. All registrations have been confirmed correctly by a dental expert. The contour-to-pixel methods can directly be used in 3D for surface-to-voxel tasks.
Image Registration of Cochlear µCT Data Using Heat Distribution Similarity
Kjer, Hans Martin; Vera, Sergio; Fagertun, Jens;
2015-01-01
Better understanding of the anatomical variability of the human cochlear is important for the design and function of Cochlear Implants. Good non-rigid alignment of high-resolution cochlear μCT data is a challenging task. In this paper we study the use of heat distribution similarity between samples...... as an anatomical registration prior. We set-up and present our heat distribution model for the cochlea and utilize it in a typical cubic B-spline registration model. Evaluation and comparison is done against a corresponding normal registration of binary segmentations....
Exponential B-splines and the partition of unity property
Christensen, Ole; Massopust, Peter
2012-01-01
We provide an explicit formula for a large class of exponential B-splines. Also, we characterize the cases where the integer-translates of an exponential B-spline form a partition of unity up to a multiplicative constant. As an application of this result we construct explicitly given pairs of dual...
Lord, Nicholas; Ho, Jeffrey; Vemuri, Baba; Eisenschenk, Stephan
2006-01-01
International audience In clinical applications where structural asymmetries between homologous shapes have been correlated with pathology, the questions of definition and quantification of 'asymmetry' arise naturally. When not only the degree but the position of deformity is thought relevant, asymmetry localization must also be addressed. Asymmetries between paired shapes can and have been formulated in the literature in terms of (nonrigid) diffeomorphisms between the shapes. For the infi...
B-spline Collocation with Domain Decomposition Method
A global B-spline collocation method has been previously developed and successfully implemented by the present authors for solving elliptic partial differential equations in arbitrary complex domains. However, the global B-spline approximation, which is simply reduced to Bezier approximation of any degree p with C0 continuity, has led to the use of B-spline basis of high order in order to achieve high accuracy. The need for B-spline bases of high order in the global method would be more prominent in domains of large dimension. For the increased collocation points, it may also lead to the ill-conditioning problem. In this study, overlapping domain decomposition of multiplicative Schwarz algorithm is combined with the global method. Our objective is two-fold that improving the accuracy with the combination technique, and also investigating influence of the combination technique to the employed B-spline basis orders with respect to the obtained accuracy. It was shown that the combination method produced higher accuracy with the B-spline basis of much lower order than that needed in implementation of the initial method. Hence, the approximation stability of the B-spline collocation method was also increased.
In this paper we describe a method to non-rigidly co-register a 2D slice sequence from real-time 3D echocardiography with a 2D cardiovascular MR image sequence. This is challenging because the imaging modalities have different spatial and temporal resolution. Non-rigid registration is required for accurate alignment due to imprecision of cardiac gating and natural motion variations between cardiac cycles. In our approach the deformation field between the imaging modalities is decoupled into temporal and spatial components. First, temporal alignment is performed to establish temporal correspondence between a real-time 3D echocardiography frame and a cardiovascular MR frame. Spatial alignment is then performed using an adaptive non-rigid registration algorithm based on local phase mutual information on each temporally aligned image pair. Experiments on seven volunteer datasets are reported. Evaluation of registration errors based on expert-identified landmarks shows that the spatio-temporal registration algorithm gives a mean registration error of 3.56 ± 0.49 and 3.54 ± 0.27 mm for the short and long axis sequences, respectively.
A Comparison of FFD-based Nonrigid Registration and AAMs Applied to Myocardial Perfusion MRI
Ólafsdóttir, Hildur; Stegmann, Mikkel Bille; Ersbøll, Bjarne Kjær;
2006-01-01
-slice myocardial perfusion images. The images are acquired by magnetic resonance imaging, from infarct patients. A registration of these sequences is crucial for clinical practice, which currently is subjected to manual labor. In the paper, the pros and cons of the two registration approaches are discussed and...
Knot Insertion Algorithms for ECT B-spline Curves
SONG Huan-huan; TANG Yue-hong; LI Yu-juan
2013-01-01
Knot insertion algorithm is one of the most important technologies of B-spline method. By inserting a knot the local prop-erties of B-spline curve and the control flexibility of its shape can be further improved, also the segmentation of the curve can be re-alized. ECT spline curve is drew by the multi-knots spline curve with associated matrix in ECT spline space;Muehlbach G and Tang Y and many others have deduced the existence and uniqueness of the ECT spline function and developed many of its important properties .This paper mainly focuses on the knot insertion algorithm of ECT B-spline curve.It is the widest popularization of B-spline Behm algorithm and theory. Inspired by the Behm algorithm, in the ECT spline space, structure of generalized Pólya poly-nomials and generalized de Boor Fix dual functional, expressing new control points which are inserted after the knot by linear com-bination of original control vertex the single knot, and there are two cases, one is the single knot, the other is the double knot. Then finally comes the insertion algorithm of ECT spline curve knot. By application of the knot insertion algorithm, this paper also gives out the knot insertion algorithm of four order geometric continuous piecewise polynomial B-spline and algebraic trigonometric spline B-spline, which is consistent with previous results.
Lin Muqing; Chen, Jeon-Hor; Bahri, Shadfar; Nalcioglu, Orhan; Su Minying [Tu and Yuen Center for Functional Onco-Imaging, Department of Radiological Sciences, University of California, Irvine, CA (United States); Mehta, Rita S [Department of Medicine, University of California, Irvine, CA (United States); Chan Siwa, E-mail: msu@uci.edu [Department of Radiology, Taichung Veterans General Hospital, Taichung, Taiwan (China)
2011-09-21
Breast MRI acquires many images from the breast, and computer-aided algorithms and display tools are often used to assist the radiologist's interpretation. Women with lifetime risk greater than 20% of developing breast cancer are recommended to receive annual screening MRI, but the current breast MRI computer-aided-diagnosis systems do not provide the necessary function for comparison of images acquired at different times. The purpose of this work was to develop registration methods for evaluating the spatial change pattern of fibroglandular tissue between two breast MRI scans of the same woman taken at different times. The registration method is based on rigid alignment followed by a non-rigid Demons algorithm. The method was tested on three different subjects who had different degrees of changes in the fibroglandular tissue, including two patients who showed different spatial shrinkage patterns after receiving neoadjuvant chemotherapy before surgery, and one control case from a normal volunteer. Based on the transformation matrix, the collapse of multiple voxels on the baseline images to one voxel on the follow-up images is used to calculate the shrinkage factor. Conversely, based on the reverse transformation matrix the expansion factor can be calculated. The shrinkage/expansion factor, the deformation magnitude and direction, as well as the Jacobian determinate at each location can be displayed in a 3D rendering view to show the spatial changes between two MRI scans. These different parameters show consistent results and can be used for quantitative evaluation of the spatial change patterns. The presented registration method can be further developed into a clinical tool for evaluating therapy-induced changes and for early diagnosis of breast cancer in screening MRI.
Breast MRI acquires many images from the breast, and computer-aided algorithms and display tools are often used to assist the radiologist's interpretation. Women with lifetime risk greater than 20% of developing breast cancer are recommended to receive annual screening MRI, but the current breast MRI computer-aided-diagnosis systems do not provide the necessary function for comparison of images acquired at different times. The purpose of this work was to develop registration methods for evaluating the spatial change pattern of fibroglandular tissue between two breast MRI scans of the same woman taken at different times. The registration method is based on rigid alignment followed by a non-rigid Demons algorithm. The method was tested on three different subjects who had different degrees of changes in the fibroglandular tissue, including two patients who showed different spatial shrinkage patterns after receiving neoadjuvant chemotherapy before surgery, and one control case from a normal volunteer. Based on the transformation matrix, the collapse of multiple voxels on the baseline images to one voxel on the follow-up images is used to calculate the shrinkage factor. Conversely, based on the reverse transformation matrix the expansion factor can be calculated. The shrinkage/expansion factor, the deformation magnitude and direction, as well as the Jacobian determinate at each location can be displayed in a 3D rendering view to show the spatial changes between two MRI scans. These different parameters show consistent results and can be used for quantitative evaluation of the spatial change patterns. The presented registration method can be further developed into a clinical tool for evaluating therapy-induced changes and for early diagnosis of breast cancer in screening MRI.
Validation of TMJ osteoarthritis synthetic defect database via non-rigid registration
Paniagua, Beatriz; Pera, Juliette; Budin, Francois; Gomes, Liliane; Styner, Martin; Lucia, Cevidanes; Nguyen, Tung
2015-03-01
Temporomandibular joint (TMJ) disorders are a group of conditions that cause pain and dysfunction in the jaw joint and the muscles controlling jaw movement. However, diagnosis and treatment of these conditions remain controversial. To date, there is no single sign, symptom, or test that can clearly diagnose early stages of osteoarthritis (OA). Instead, the diagnosis is based on a consideration of several factors, including radiological evaluation. The current radiological diagnosis scores of TMJ pathology are subject to misdiagnosis. We believe these scores are limited by the acquisition procedures, such as oblique cuts of the CT and head positioning errors, and can lead to incorrect diagnoses of flattening of the head of the condyle, formation of osteophytes, or condylar pitting. This study consists of creating and validating a methodological framework to simulate defects in CBCT scans of known location and size, in order to create synthetic TMJ OA database. User-generated defects were created using a non-rigid deformation protocol in CBCT. All segmentation evaluation, surface distances and linear distances from the user-generated to the simulated defects showed our methodological framework to be very precise and within a voxel (0.5 mm) of magnitude. A TMJ OA synthetic database will be created next, and evaluated by expert radiologists, and this will serve to evaluate how sensitive the current radiological diagnosis tools are.
Purpose: Treatments like radiotherapy and focused ultrasound in the abdomen require accurate motion tracking, in order to optimize dosage delivery to the target and minimize damage to critical structures and healthy tissues around the target. 4D ultrasound is a promising modality for motion tracking during such treatments. In this study, the authors evaluate the accuracy of motion tracking in the liver based on deformable registration of 4D ultrasound images. Methods: The offline analysis was performed using a nonrigid registration algorithm that was specifically designed for motion estimation from dynamic imaging data. The method registers the entire 4D image data sequence in a groupwise optimization fashion, thus avoiding a bias toward a specifically chosen reference time point. Three healthy volunteers were scanned over several breathing cycles (12 s) from three different positions and angles on the abdomen; a total of nine 4D scans for the three volunteers. Well-defined anatomic landmarks were manually annotated in all 96 time frames for assessment of the automatic algorithm. The error of the automatic motion estimation method was compared with interobserver variability. The authors also performed experiments to investigate the influence of parameters defining the deformation field flexibility and evaluated how well the method performed with a lower temporal resolution in order to establish the minimum frame rate required for accurate motion estimation. Results: The registration method estimated liver motion with an error of 1 mm (75% percentile over all datasets), which was lower than the interobserver variability of 1.4 mm. The results were only slightly dependent on the degrees of freedom of the deformation model. The registration error increased to 2.8 mm with an eight times lower temporal resolution. Conclusions: The authors conclude that the methodology was able to accurately track the motion of the liver in the 4D ultrasound data. The authors believe
Conventional radiotherapy is planned using free-breathing computed tomography (CT), ignoring the motion and deformation of the anatomy from respiration. New breath-hold-synchronized, gated, and four-dimensional (4D) CT acquisition strategies are enabling radiotherapy planning utilizing a set of CT scans belonging to different phases of the breathing cycle. Such 4D treatment planning relies on the availability of tumor and organ contours in all phases. The current practice of manual segmentation is impractical for 4D CT, because it is time consuming and tedious. A viable solution is registration-based segmentation, through which contours provided by an expert for a particular phase are propagated to all other phases while accounting for phase-to-phase motion and anatomical deformation. Deformable image registration is central to this task, and a free-form deformation-based nonrigid image registration algorithm will be presented. Compared with the original algorithm, this version uses novel, computationally simpler geometric constraints to preserve the topology of the dense control-point grid used to represent free-form deformation and prevent tissue fold-over. Using mean squared difference as an image similarity criterion, the inhale phase is registered to the exhale phase of lung CT scans of five patients and of characteristically low-contrast abdominal CT scans of four patients. In addition, using expert contours for the inhale phase, the corresponding contours were automatically generated for the exhale phase. The accuracy of the segmentation (and hence deformable image registration) was judged by comparing automatically segmented contours with expert contours traced directly in the exhale phase scan using three metrics: volume overlap index, root mean square distance, and Hausdorff distance. The accuracy of the segmentation (in terms of radial distance mismatch) was approximately 2 mm in the thorax and 3 mm in the abdomen, which compares favorably to the
Vijayan, Sinara, E-mail: sinara.vijayan@ntnu.no [Norwegian University of Science and Technology, 7491 Trondheim (Norway); Klein, Stefan [Norwegian University of Science and Technology, 7491 Trondheim, Norway and Biomedical Imaging Group Rotterdam, Department of Medical Informatics and Radiology, Erasmus MC, 3000 CA Rotterdam (Netherlands); Hofstad, Erlend Fagertun; Langø, Thomas [SINTEF, Department Medical Technology, 7465 Trondheim (Norway); Lindseth, Frank [Norwegian University of Science and Technology, 7491 Trondheim, Norway and SINTEF, Department Medical Technology, 7465 Trondheim (Norway); Ystgaard, Brynjulf [Department of Surgery, St. Olavs Hospital, 7030 Trondheim (Norway)
2014-08-15
Purpose: Treatments like radiotherapy and focused ultrasound in the abdomen require accurate motion tracking, in order to optimize dosage delivery to the target and minimize damage to critical structures and healthy tissues around the target. 4D ultrasound is a promising modality for motion tracking during such treatments. In this study, the authors evaluate the accuracy of motion tracking in the liver based on deformable registration of 4D ultrasound images. Methods: The offline analysis was performed using a nonrigid registration algorithm that was specifically designed for motion estimation from dynamic imaging data. The method registers the entire 4D image data sequence in a groupwise optimization fashion, thus avoiding a bias toward a specifically chosen reference time point. Three healthy volunteers were scanned over several breathing cycles (12 s) from three different positions and angles on the abdomen; a total of nine 4D scans for the three volunteers. Well-defined anatomic landmarks were manually annotated in all 96 time frames for assessment of the automatic algorithm. The error of the automatic motion estimation method was compared with interobserver variability. The authors also performed experiments to investigate the influence of parameters defining the deformation field flexibility and evaluated how well the method performed with a lower temporal resolution in order to establish the minimum frame rate required for accurate motion estimation. Results: The registration method estimated liver motion with an error of 1 mm (75% percentile over all datasets), which was lower than the interobserver variability of 1.4 mm. The results were only slightly dependent on the degrees of freedom of the deformation model. The registration error increased to 2.8 mm with an eight times lower temporal resolution. Conclusions: The authors conclude that the methodology was able to accurately track the motion of the liver in the 4D ultrasound data. The authors believe
Fingerprint Representation Methods Based on B-Spline Functions
Ruan Ke; Xia De-lin; Yan Pu-liu
2004-01-01
The global characteristics of a fingerprint image such as the ridge shape and ridge topology are often ignored in most automatic fingerprint verification system. In this paper, a new representative method based on B-Spline curve is proposed to address this problem. The resultant B-Spline curves can represent the global characteristics completely and the curves are analyzable and precise. An algorithm is also proposed to extract the curves from the fingerprint image. In addition to preserve the most information of the fingerprint image, the knot-points number of the B-Spline curve is reduced to minimum in this algorithm. At the same time, the influence of the fingerprint image noise is discussed. In the end, an example is given to demonstrate the effectiveness of the representation method.
Non-rigid registration of CLSM images of physical sections with discontinuous deformations
Michálek, Jan; Čapek, Martin; Kubínová, Lucie
Urbino : Societa Italiana Scienze Micriscopiche, 2011. s. 201-202 R&D Projects: GA ČR(CZ) GA102/08/0691; GA ČR(CZ) GA304/09/0733; GA MŠk(CZ) LC06063; GA MŠk(CZ) ME09010; GA MŠk(CZ) MEB0810092 Institutional research plan: CEZ:AV0Z50110509 Keywords : confocal microscopy * image registration * multilabel optimization Subject RIV: JC - Computer Hardware ; Software
Nonrigid Registration of CLSM Images of Physical Sections with Discontinuous Deformations
Michálek, Jan; Čapek, Martin; Kubínová, Lucie
2011-01-01
Roč. 17, č. 6 (2011), s. 923-936. ISSN 1431-9276 R&D Projects: GA ČR(CZ) GA102/08/0691; GA ČR(CZ) GA304/09/0733; GA MŠk(CZ) LC06063; GA MŠk(CZ) ME09010; GA MŠk(CZ) MEB0810092 Institutional research plan: CEZ:AV0Z50110509 Keywords : confocal microscopy * image registration * target image * reference image * total variation * data term * smoothness term * L1-TV function * E-L1-TV function * graph cuts * multilabel optimization * local entropy Subject RIV: JC - Computer Hardware ; Software Impact factor: 3.007, year: 2011
Nonrigid registration-based coronary artery motion correction for cardiac computed tomography
Bhagalia, Roshni; Pack, Jed D.; Miller, James V.; Iatrou, Maria [GE Global Research, Niskayuna, New York 12309 (United States); GE Healthcare, Hoboken, New Jersey 07030 (United States)
2012-07-15
Purpose: X-ray computed tomography angiography (CTA) is the modality of choice to noninvasively monitor and diagnose heart disease with coronary artery health and stenosis detection being of particular interest. Reliable, clinically relevant coronary artery imaging mandates high spatiotemporal resolution. However, advances in intrinsic scanner spatial resolution (CT scanners are available which combine nearly 900 detector columns with focal spot oversampling) can be tempered by motion blurring, particularly in patients with unstable heartbeats. As a result, recently numerous methods have been devised to improve coronary CTA imaging. Solutions involving hardware, multisector algorithms, or {beta}-blockers are limited by cost, oversimplifying assumptions about cardiac motion, and populations showing contraindications to drugs, respectively. This work introduces an inexpensive algorithmic solution that retrospectively improves the temporal resolution of coronary CTA without significantly affecting spatial resolution. Methods: Given the goal of ruling out coronary stenosis, the method focuses on 'deblurring' the coronary arteries. The approach makes no assumptions about cardiac motion, can be used on exams acquired at high heart rates (even over 75 beats/min), and draws on a fast and accurate three-dimensional (3D) nonrigid bidirectional labeled point matching approach to estimate the trajectories of the coronary arteries during image acquisition. Motion compensation is achieved by employing a 3D warping of a series of partial reconstructions based on the estimated motion fields. Each of these partial reconstructions is created from data acquired over a short time interval. For brevity, the algorithm 'Subphasic Warp and Add' (SWA) reconstruction. Results: The performance of the new motion estimation-compensation approach was evaluated by a systematic observer study conducted using nine human cardiac CTA exams acquired over a range of average heart
Color management with a hammer: the B-spline fitter
Bell, Ian E.; Liu, Bonny H. P.
2003-01-01
To paraphrase Abraham Maslow: If the only tool you have is a hammer, every problem looks like a nail. We have a B-spline fitter customized for 3D color data, and many problems in color management can be solved with this tool. Whereas color devices were once modeled with extensive measurement, look-up tables and trilinear interpolation, recent improvements in hardware have made B-spline models an affordable alternative. Such device characterizations require fewer color measurements than piecewise linear models, and have uses beyond simple interpolation. A B-spline fitter, for example, can act as a filter to remove noise from measurements, leaving a model with guaranteed smoothness. Inversion of the device model can then be carried out consistently and efficiently, as the spline model is well behaved and its derivatives easily computed. Spline-based algorithms also exist for gamut mapping, the composition of maps, and the extrapolation of a gamut. Trilinear interpolation---a degree-one spline---can still be used after nonlinear spline smoothing for high-speed evaluation with robust convergence. Using data from several color devices, this paper examines the use of B-splines as a generic tool for modeling devices and mapping one gamut to another, and concludes with applications to high-dimensional and spectral data.
Rexilius Jan
2005-03-01
Full Text Available Abstract Background The visual combination of different modalities is essential for many medical imaging applications in the field of Computer-Assisted medical Diagnosis (CAD to enhance the clinical information content. Clinically, incontinence is a diagnosis with high clinical prevalence and morbidity rate. The search for a method to identify risk patients and to control the success of operations is still a challenging task. The conjunction of magnetic resonance (MR and 3D ultrasound (US image data sets could lead to a new clinical visual representation of the morphology as we show with corresponding data sets of the female anal canal with this paper. Methods We present a feasibility study for a non-rigid registration technique based on a biomechanical model for MR and US image data sets of the female anal canal as a base for a new innovative clinical visual representation. Results It is shown in this case study that the internal and external sphincter region could be registered elastically and the registration partially corrects the compression induced by the ultrasound transducer, so the MR data set showing the native anatomy is used as a frame for the US data set showing the same region with higher resolution but distorted by the transducer Conclusion The morphology is of special interest in the assessment of anal incontinence and the non-rigid registration of normal clinical MR and US image data sets is a new field of the adaptation of this method incorporating the advantages of both technologies.
RECONSTRUCTION OF LAYER DATA WITH DEFORMABLE B-SPLINES
Cheng Siyuan; Zhang Xiangwei; Xiong Hanwei
2005-01-01
A new B-spline surface reconstruction method from layer data based on deformable model is presented. An initial deformable surface, which is represented as a closed cylinder, is firstly given. The surface is subject to internal forces describing its implicit smoothness property and external forces attracting it toward the layer data points. And then finite element method is adopted to solve its energy minimization problem, which results a bicubic closed B-spline surface with C2 continuity. The proposed method can provide a smoothness and accurate surface model directly from the layer data, without the need to fit cross-sectional curves and make them compatible. The feasibility of the proposed method is verified by the experimental results.
Tatsugami, Fuminari; Higaki, Toru; Nakamura, Yuko; Yamagami, Takuji; Date, Shuji; Awai, Kazuo [Hiroshima University, Department of Diagnostic Radiology, Minami-ku, Hiroshima (Japan); Fujioka, Chikako; Kiguchi, Masao [Hiroshima University, Department of Radiology, Minami-ku, Hiroshima (Japan); Kihara, Yasuki [Hiroshima University, Department of Cardiovascular Medicine, Minami-ku, Hiroshima (Japan)
2015-01-15
To investigate the feasibility of a newly developed noise reduction technique at coronary CT angiography (CTA) that uses multi-phase data-averaging and non-rigid image registration. Sixty-five patients underwent coronary CTA with prospective ECG-triggering. The range of the phase window was set at 70-80 % of the R-R interval. First, three sets of consecutive volume data at 70 %, 75 % and 80 % of the R-R interval were prepared. Second, we applied non-rigid registration to align the 70 % and 80 % images to the 75 % image. Finally, we performed weighted averaging of the three images and generated a de-noised image. The image noise and contrast-to-noise ratio (CNR) in the proximal coronary arteries between the conventional 75 % and the de-noised images were compared. Two radiologists evaluated the image quality using a 5-point scale (1, poor; 5, excellent). On de-noised images, mean image noise was significantly lower than on conventional 75 % images (18.3 HU ± 2.6 vs. 23.0 HU ± 3.3, P < 0.01) and the CNR was significantly higher (P < 0.01). The mean image quality score for conventional 75 % and de-noised images was 3.9 and 4.4, respectively (P < 0.01). Our method reduces image noise and improves image quality at coronary CTA. (orig.)
4D ultrasound and 3D MRI registration of beating heart
To realize intra-cardiac surgery without cardio-pulmonary bypass, a medical imaging technique with both high image quality and data acquisition rate that is fast enough to follow heart beat movements is required. In this research, we proposed a method that utilized the image quality of MRI and the speed of ultrasound. We developed a 4D image reconstruction method using image registration of 3D MRI and 4D ultrasound images. The registration method consists of rigid registration between 3D MRI and 3D ultrasound with the same heart beat phase, and non-rigid registration between 3D ultrasound images from different heart beat phases. Non-rigid registration was performed with B-spline based registration using variable spring model. In phantom experiment using balloon phantom, registration accuracy was less than 2 mm for total heart volume variation range of 10%. We applied our registration method on 3D MRI and 4D ultrasound images of a volunteer's beating heart data and confirmed through visual observation that heart beat pattern was well reproduced. (orig.)
Research on Non-rigid Registration in Uterine Fibroids MR Imaging%子宫肌瘤MRI影像非刚性配准研究
所世腾; 张良宾; 李跃华; 陈亚珠; 张素
2012-01-01
高强度聚焦超声(High-Intensity Focused Ultrasound,HIFU)治疗和磁共振技术结合(MRl-guided HIFU,MRgHIFU)采用MRI进行目标定位、治疗规划和能量沉积的闭环控制以保障热消融不伤及周围组织,其中图像配准是校正定位误差,实施精确治疗的重要环节.针对三维非刚性配准方法,在子宫肌瘤的治疗计划修正和跟踪方面的应用进行研究.针对前后两个不同时段采集的子宫肌瘤MR影像,分别进行基于自由形变模型(Free-From Deformation,FFD)算法和Demons算法的非刚性配准对比.实验结果表明,该基于FFD的非刚性配准算法,对于形变较小的子宫肌瘤真实数据能够取得较为合适的配准效果,重叠区域的互相关系数(CC)从配准前的0.59提高到配准后的0.74.%High-Intensity Focused Ultrasound (H1FU) combined with magnetic resonance achieves thermal ablation without hurting the surrounding tissue through a closed-loop control process of target positioning and treatment planning. Image registration is a vital important section in correction of positioning error and implement of precise treatment. Research of non-rigid registration in amendments and tracing of treatment plan was carried out, in which series of uterine fibroids MR images acquired in two different times are employed as experimental data. Then non-rigid registration based on Free-Form Deformation and Demons were performed successively for a comparison of the results. The experiments shows that the FFD based registration achieved preferable results with tiny deformation, and the correlation coefficient (CC) of overlapping region increased from 0.59 to 0.74 after registration.
Hodograph computation and bound estimation for rational B-spline curves
无
2007-01-01
It is necessary to compute the derivative and estimate the bound of rational B-spline curves in design system, which has not been studied to date. To improve the function of computer aided design (CAD) system, and to enhance the efficiency of different algorithms of rational B-spline curves, the representation of scaled hodograph and bound of derivative magnitude of uniform planar rational B-spline curves are derived by applying Dir function, which indicates the direction of Cartesian vector between homogeneous points, discrete B-spline theory and the formula of translating the product into a summation of B-spline functions. As an application of the result above,upper bound of parametric distance between any two points in a uniform planar rational B-spline curve is further presented.
Non-rigid registration of small animal skeletons from micro-CT using 3D shape context
Xiao, Di; Bourgeat, Pierrick; Fripp, Jurgen; Acosta Tamayo, Oscar; Gregoire, Marie Claude; Salvado, Olivier
2009-02-01
Small animal registration is an important step for molecular image analysis. Skeleton registration from whole-body or only partial micro Computerized Tomography (CT) image is often performed to match individual rats to atlases and templates, for example to identify organs in positron emission tomography (PET). In this paper, we extend the shape context matching technique for 3D surface registration and apply it for rat hind limb skeleton registration from CT images. Using the proposed method, after standard affine iterative closest point (ICP) registration, correspondences between the 3D points from sour and target objects were robustly found and used to deform the limb skeleton surface with thin-plate-spline (TPS). Experiments are described using phantoms and actual rat hind limb skeletons. On animals, mean square errors were decreased by the proposed registration compared to that of its initial alignment. Visually, skeletons were successfully registered even in cases of very different animal poses.
Grgic, Aleksandar; Ballek, Elena; Moca, Norbert; Schaefer, Andrea; Kirsch, Carl-Martin; Hellwig, Dirk [Saarland University Medical Center, Department of Nuclear Medicine, Homburg/Saar (Germany); Fleckenstein, Jochen; Kremp, Stephanie; Ruebe, Christian [Saarland University Medical Center, Department of Radio-Oncology, Homburg/Saar (Germany); Kuhnigk, Jan-Martin [Fraunhofer, MEVIS, Bremen (Germany)
2011-05-15
Assessment of the metabolically active tumour tissue by FDG PET is evolving for use in the diagnosis of non-small-cell lung cancer (NSCLC), in the planning of radiotherapy, and in follow-up and response evaluation. For exact evaluation accurate registration of PET and CT data is required. The registration process is usually based on rigid algorithms; however, nonrigid algorithms are increasingly being used. The influence of the registration method on FDG PET-based standardized uptake value (SUVmax) and metabolic tumour volume (MTV) definition has not yet been evaluated. We compared intra- and interindividual differences in SUV and MTV between rigid- and nonrigid-registered PET and CT acquired during different breathing manoeuvres. The study group comprised 28 radiotherapy candidates with histologically proven NSCLC who underwent FDG PET acquisition and three CT acquisitions (expiration - EXP, inspiration - INS, mid-breath-hold - MID). All scans were registered with both a rigid (R) and a nonrigid (NR) procedure resulting in six fused datasets: R-INS, R-EXP, R-MID, NR-INS, NR-EXP and NR-MID. For the delineation of MTVs a contrast-oriented contouring algorithm developed in-house was used. To accelerate the delineation a semiautomatic software prototype was utilized. Tumour mean SUVmax did not differ for R and NR registration (R 17.5 {+-} 7, NR 17.4 {+-} 7; p=0.2). The mean MTV was higher by 3 {+-} 12 ml (p=0.02) in the NR group than in the R group, as was the mean tumour diameter (by 0.1 {+-} 0.2 cm; p<0.01). With respect to the three different breathing manoeuvres, there were no differences in MTV in the R group (p > 0.7). In intraindividual comparison there were no significant differences in MTVs concerning the registration pairs R-EXP (68 {+-} 88 ml) vs. NR-EXP (69 {+-} 85 ml) und R-MID (68 {+-} 86 ml) vs. NR-MID (69 {+-} 83 ml) (both p > 0.4). However, the MTVs were larger after NR registration during inspiration (R-INS 68 {+-} 82 vs. NR-INS 78 {+-} 93 ml; p=0
Image edges detection through B-Spline filters
B-Spline signal processing was used to detect the edges of a digital image. This technique is based upon processing the image in the Spline transform domain, instead of doing so in the space domain (classical processing). The transformation to the Spline transform domain means finding out the real coefficients that makes it possible to interpolate the grey levels of the original image, with a B-Spline polynomial. There exist basically two methods of carrying out this interpolation, which produces the existence of two different Spline transforms: an exact interpolation of the grey values (direct Spline transform), and an approximated interpolation (smoothing Spline transform). The latter results in a higher smoothness of the gray distribution function defined by the Spline transform coefficients, and is carried out with the aim of obtaining an edge detection algorithm which higher immunity to noise. Finally the transformed image was processed in order to detect the edges of the original image (the gradient method was used), and the results of the three methods (classical, direct Spline transform and smoothing Spline transform) were compared. The results were that, as expected, the smoothing Spline transform technique produced a detection algorithm more immune to external noise. On the other hand the direct Spline transform technique, emphasizes more the edges, even more than the classical method. As far as the consuming time is concerned, the classical method is clearly the fastest one, and may be applied whenever the presence of noise is not important, and whenever edges with high detail are not required in the final image. (author). 9 refs., 17 figs., 1 tab
B-spline collocation methods for numerical solutions of the Burgers' equation
İdris Dağ; Dursun Irk; Ali Şahin
2005-01-01
Both time- and space-splitted Burgers' equations are solved numerically. Cubic B-spline collocation method is applied to the time-splitted Burgers' equation. Quadratic B-spline collocation method is used to get numerical solution of the space-splitted Burgers' equation. The results of both schemes are compared for some test problems.
Carvalho, Diego D. B.; Akkus, Zeynettin; Bosch, Johan G.; van den Oord, Stijn C. H.; Niessen, Wiro J.; Klein, Stefan
2014-03-01
In this work, we investigate nonrigid motion compensation in simultaneously acquired (side-by-side) B-mode ultrasound (BMUS) and contrast enhanced ultrasound (CEUS) image sequences of the carotid artery. These images are acquired to study the presence of intraplaque neovascularization (IPN), which is a marker of plaque vulnerability. IPN quantification is visualized by performing the maximum intensity projection (MIP) on the CEUS image sequence over time. As carotid images contain considerable motion, accurate global nonrigid motion compensation (GNMC) is required prior to the MIP. Moreover, we demonstrate that an improved lumen and plaque differentiation can be obtained by averaging the motion compensated BMUS images over time. We propose to use a previously published 2D+t nonrigid registration method, which is based on minimization of pixel intensity variance over time, using a spatially and temporally smooth B-spline deformation model. The validation compares displacements of plaque points with manual trackings by 3 experts in 11 carotids. The average (+/- standard deviation) root mean square error (RMSE) was 99+/-74μm for longitudinal and 47+/-18μm for radial displacements. These results were comparable with the interobserver variability, and with results of a local rigid registration technique based on speckle tracking, which estimates motion in a single point, whereas our approach applies motion compensation to the entire image. In conclusion, we evaluated that the GNMC technique produces reliable results. Since this technique tracks global deformations, it can aid in the quantification of IPN and the delineation of lumen and plaque contours.
Micropolar Fluids Using B-spline Divergence Conforming Spaces
Sarmiento, Adel
2014-06-06
We discretized the two-dimensional linear momentum, microrotation, energy and mass conservation equations from micropolar fluids theory, with the finite element method, creating divergence conforming spaces based on B-spline basis functions to obtain pointwise divergence free solutions [8]. Weak boundary conditions were imposed using Nitsche\\'s method for tangential conditions, while normal conditions were imposed strongly. Once the exact mass conservation was provided by the divergence free formulation, we focused on evaluating the differences between micropolar fluids and conventional fluids, to show the advantages of using the micropolar fluid model to capture the features of complex fluids. A square and an arc heat driven cavities were solved as test cases. A variation of the parameters of the model, along with the variation of Rayleigh number were performed for a better understanding of the system. The divergence free formulation was used to guarantee an accurate solution of the flow. This formulation was implemented using the framework PetIGA as a basis, using its parallel stuctures to achieve high scalability. The results of the square heat driven cavity test case are in good agreement with those reported earlier.
Pheiffer, Thomas S.; Ou, Jao J.; Miga, Michael I.
2010-02-01
Modality-independent elastography (MIE) is a method of elastography that reconstructs the elastic properties of tissue using images acquired under different loading conditions and a biomechanical model. Boundary conditions are a critical input to the algorithm, and are often determined by time-consuming point correspondence methods requiring manual user input. Unfortunately, generation of accurate boundary conditions for the biomechanical model is often difficult due to the challenge of accurately matching points between the source and target surfaces and consequently necessitates the use of large numbers of fiducial markers. This study presents a novel method of automatically generating boundary conditions by non-rigidly registering two image sets with a Demons diffusion-based registration algorithm. The use of this method was successfully performed in silico using magnetic resonance and X-ray computed tomography image data with known boundary conditions. These preliminary results have produced boundary conditions with accuracy of up to 80% compared to the known conditions. Finally, these boundary conditions were utilized within a 3D MIE reconstruction to determine an elasticity contrast ratio between tumor and normal tissue. Preliminary results show a reasonable characterization of the material properties on this first attempt and a significant improvement in the automation level and viability of the method.
Adaptive B-spline volume representation of measured BRDF data for photorealistic rendering
Hyungjun Park
2015-01-01
Full Text Available Measured bidirectional reflectance distribution function (BRDF data have been used to represent complex interaction between lights and surface materials for photorealistic rendering. However, their massive size makes it hard to adopt them in practical rendering applications. In this paper, we propose an adaptive method for B-spline volume representation of measured BRDF data. It basically performs approximate B-spline volume lofting, which decomposes the problem into three sub-problems of multiple B-spline curve fitting along u-, v-, and w-parametric directions. Especially, it makes the efficient use of knots in the multiple B-spline curve fitting and thereby accomplishes adaptive knot placement along each parametric direction of a resulting B-spline volume. The proposed method is quite useful to realize efficient data reduction while smoothing out the noises and keeping the overall features of BRDF data well. By applying the B-spline volume models of real materials for rendering, we show that the B-spline volume models are effective in preserving the features of material appearance and are suitable for representing BRDF data.
Adaptive non-uniform B-spline dictionaries on a compact interval
Rebollo-Neira, Laura
2009-01-01
Non-uniform B-spline dictionaries on a compact interval are discussed. For each given partition, dictionaries of B-spline functions for the corresponding spline space are constructed. It is asserted that, by dividing the given partition into subpartitions and joining together the bases for the concomitant subspaces, slightly redundant dictionaries of B-splines functions are obtained. Such dictionaries are proved to span the spline space associated to the given partition. The proposed construction is shown to be potentially useful for the purpose of sparse signal representation. With that goal in mind, spline spaces specially adapted to produce a sparse representation of a given signal are considered.
Numerical Method Using Cubic Trigonometric B-Spline Technique for Nonclassical Diffusion Problems
Muhammad Abbas; Ahmad Abd. Majid; Ahmad Izani Md Ismail; Abdur Rashid
2014-01-01
A new two-time level implicit technique based on cubic trigonometric B-spline is proposed for the approximate solution of a nonclassical diffusion problem with nonlocal boundary constraints. The standard finite difference approach is applied to discretize the time derivative while cubic trigonometric B-spline is utilized as an interpolating function in the space dimension. The technique is shown to be unconditionally stable using the von Neumann method. Several numerical examples are discusse...
Adaptive B-spline volume representation of measured BRDF data for photorealistic rendering
Hyungjun Park; Joo-Haeng Lee
2015-01-01
Measured bidirectional reflectance distribution function (BRDF) data have been used to represent complex interaction between lights and surface materials for photorealistic rendering. However, their massive size makes it hard to adopt them in practical rendering applications. In this paper, we propose an adaptive method for B-spline volume representation of measured BRDF data. It basically performs approximate B-spline volume lofting, which decomposes the problem into three sub-problems of mu...
B-SPLINE-BASED SVM MODEL AND ITS APPLICATIONS TO OIL WATER-FLOODED STATUS IDENTIFICATION
Shang Fuhua; Zhao Tiejun; Yi Xiongying
2007-01-01
A method of B-spline transform for signal feature extraction is developed. With the B-spline,the log-signal space is mapped into the vector space. An efficient algorithm based on Support Vector Machine (SVM) to automatically identify the water-flooded status of oil-saturated stratum is described.The experiments show that this algorithm can improve the performances for the identification and the generalization in the case of a limited set of samples.
BSR: B-spline atomic R-matrix codes
Zatsarinny, Oleg
2006-02-01
BSR is a general program to calculate atomic continuum processes using the B-spline R-matrix method, including electron-atom and electron-ion scattering, and radiative processes such as bound-bound transitions, photoionization and polarizabilities. The calculations can be performed in LS-coupling or in an intermediate-coupling scheme by including terms of the Breit-Pauli Hamiltonian. New version program summaryTitle of program: BSR Catalogue identifier: ADWY Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADWY Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland Computers on which the program has been tested: Microway Beowulf cluster; Compaq Beowulf cluster; DEC Alpha workstation; DELL PC Operating systems under which the new version has been tested: UNIX, Windows XP Programming language used: FORTRAN 95 Memory required to execute with typical data: Typically 256-512 Mwords. Since all the principal dimensions are allocatable, the available memory defines the maximum complexity of the problem No. of bits in a word: 8 No. of processors used: 1 Has the code been vectorized or parallelized?: no No. of lines in distributed program, including test data, etc.: 69 943 No. of bytes in distributed program, including test data, etc.: 746 450 Peripherals used: scratch disk store; permanent disk store Distribution format: tar.gz Nature of physical problem: This program uses the R-matrix method to calculate electron-atom and electron-ion collision processes, with options to calculate radiative data, photoionization, etc. The calculations can be performed in LS-coupling or in an intermediate-coupling scheme, with options to include Breit-Pauli terms in the Hamiltonian. Method of solution: The R-matrix method is used [P.G. Burke, K.A. Berrington, Atomic and Molecular Processes: An R-Matrix Approach, IOP Publishing, Bristol, 1993; P.G. Burke, W.D. Robb, Adv. At. Mol. Phys. 11 (1975) 143; K.A. Berrington, W.B. Eissner, P.H. Norrington, Comput
Sacharow, Alexei
2011-12-01
The unprecedented success of the iterative closest point (ICP) method for registration in geometry processing and related fields can be attributed to its efficiency, robustness, and wide spectrum of applications. Its use is however quite limited as soon as the objects to be registered arise from each other by a transformation significantly different from a Euclidean motion. We present a novel variant of ICP, tailored for the specific needs of production engineering, which registers a triangle mesh with a second surface model of arbitrary digital representation. Our method inherits most of ICP\\'s practical advantages but is capable of detecting medium-strength bendings i.e. isometric deformations. Initially, the algorithm assigns to all vertices in the source their closest point on the target mesh and then iteratively establishes isometry, a process which, very similar to ICP, requires intermediate re-projections. A NURBS-based technique for applying the resulting deformation to arbitrary instances of the source geometry, other than the very mesh used for correspondence estimation, is described before we present numerical results on synthetic and real data to underline the viability of our approach in comparison with others. © 2011 Elsevier Ltd. All rights reserved.
Yang, Xiaofeng, E-mail: xyang43@emory.edu; Rossi, Peter; Ogunleye, Tomi; Marcus, David M.; Jani, Ashesh B.; Curran, Walter J.; Liu, Tian [Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, Georgia 30322 (United States); Mao, Hui [Department of Radiology and Imaging Sciences, Emory University, Atlanta, Georgia 30322 (United States)
2014-11-01
Purpose: The technological advances in real-time ultrasound image guidance for high-dose-rate (HDR) prostate brachytherapy have placed this treatment modality at the forefront of innovation in cancer radiotherapy. Prostate HDR treatment often involves placing the HDR catheters (needles) into the prostate gland under the transrectal ultrasound (TRUS) guidance, then generating a radiation treatment plan based on CT prostate images, and subsequently delivering high dose of radiation through these catheters. The main challenge for this HDR procedure is to accurately segment the prostate volume in the CT images for the radiation treatment planning. In this study, the authors propose a novel approach that integrates the prostate volume from 3D TRUS images into the treatment planning CT images to provide an accurate prostate delineation for prostate HDR treatment. Methods: The authors’ approach requires acquisition of 3D TRUS prostate images in the operating room right after the HDR catheters are inserted, which takes 1–3 min. These TRUS images are used to create prostate contours. The HDR catheters are reconstructed from the intraoperative TRUS and postoperative CT images, and subsequently used as landmarks for the TRUS–CT image fusion. After TRUS–CT fusion, the TRUS-based prostate volume is deformed to the CT images for treatment planning. This method was first validated with a prostate-phantom study. In addition, a pilot study of ten patients undergoing HDR prostate brachytherapy was conducted to test its clinical feasibility. The accuracy of their approach was assessed through the locations of three implanted fiducial (gold) markers, as well as T2-weighted MR prostate images of patients. Results: For the phantom study, the target registration error (TRE) of gold-markers was 0.41 ± 0.11 mm. For the ten patients, the TRE of gold markers was 1.18 ± 0.26 mm; the prostate volume difference between the authors’ approach and the MRI-based volume was 7.28% ± 0
Calculation of energy levels and wavefunctions of hydrogen molecular ion using B-splines function
Energy levels and wavefunctions of the ground state and the first excited state of hydrogen molecular ion are calculated by solving stationary Schrodinger equation with B-splines functions. By adopting nuclear positions as knots of B-splines basis, high accuracy of energy levels of the ground state and the first excited state for hydrogen molecular ion can be reached even for the larger internuclear separations, and our ι dependent radial wavefunctions of the ground state are in a good agreement with those computed from GAUSSIAN chemistry software. (authors)
A weighted extended B-spline solver for bending and buckling of stiffened plates
Verschaeve, Joris C G
2015-01-01
The weighted extended B-spline method [Hoellig (2003)] is applied to bending and buckling problems of plates with different shapes and stiffener arrangements. The discrete equations are obtained from the energy contributions of the different components constituting the system by means of the Rayleigh-Ritz approach. The pre-buckling or plane stress is computed by means of Airy's stress function. A boundary data extension algorithm for the weighted extended B-spline method is derived in order to solve for inhomogeneous Dirichlet boundary conditions. A series of benchmark tests is performed touching various aspects influencing the accuracy of the method.
Cardiac motion tracking with multilevel B-splines and SinMod from tagged MRI
Wang, Hui; Amini, Amir A.
2011-03-01
Cardiac motion analysis can play an important role in cardiac disease diagnosis. Tagged magnetic resonance imaging (MRI) has the ability to directly and non-invasively alter tissue magnetization and produce tags on the deforming tissue. This paper proposes an approach to analysis of tagged MR images using a multilevel B-splines fitting model incorporating phase information. The novelty of the proposed technique is that phase information is extracted from SinMod.1 By using real tag intersections extracted directly from tagged MR image data and virtual tag intersections extracted from phase information, both considered to be scattered data, multilevel B-spline fitting can result in accurate displacement motion fields. The B-spline approximation which also serves to remove noise in the displacement measurements is performed without specifying control point locations explicitly and is very fast. Dense virtual tag intersections based on SinMod were created and incorporated into the multilevel B-spline fitting process. Experimental results on simulated data from the 13- parameter kinematic model of Arts et al.2 and in vivo canine data demonstrate further improvement in accuracy and effectiveness of the proposed method.
Calculations of Electron Structure of Endohedrally Confined Helium Atom with B-Spline Type Functions
QIAO HaoXue; SHI TingYun; LI BaiWen
2002-01-01
The B-spline basis set method is used to study the properties of helium confined endohedrally at thegeometrical centre of a fullerene. The boundary conditions of the wavefunctions can be simply satisfied with thismethod. From our results, the phenomenon of "mirror collapse" is found in the case of confining helium. The interestingbehaviors of confining helium are also discussed.
A Mixed Basis Density Functional Approach for Low Dimensional Systems with B-splines
Ren, Chung-Yuan; Chang, Yia-Chung
2014-01-01
A mixed basis approach based on density functional theory is employed for low dimensional systems. The basis functions are taken to be plane waves for the periodic direction multiplied by B-spline polynomials in the non-periodic direction. B-splines have the following advantages:(1) the associated matrix elements are sparse, (2) B-splines possess a superior treatment of derivatives, (3) B-splines are not associated with atomic positions when the geometry structure is optimized, making the geometry optimization easy to implement. With this mixed basis set we can directly calculate the total energy of the system instead of using the conventional supercell model with a slab sandwiched between vacuum regions. A generalized Lanczos-Krylov iterative method is implemented for the diagonalization of the Hamiltonian matrix. To demonstrate the present approach, we apply it to study the C(001)-(2x1) surface with the norm-conserving pseudopotential, the n-type delta-doped graphene, and graphene nanoribbon with Vanderbilt...
Fast Parallel Image Registration on CPU and GPU for Diagnostic Classification of Alzheimer's Disease
Denis P Shamonin
2014-01-01
Full Text Available Nonrigid image registration is an important, but time-consuming taskin medical image analysis. In typical neuroimaging studies, multipleimage registrations are performed, i.e. for atlas-based segmentationor template construction. Faster image registration routines wouldtherefore be beneficial.In this paper we explore acceleration of the image registrationpackage elastix by a combination of several techniques: iparallelization on the CPU, to speed up the cost function derivativecalculation; ii parallelization on the GPU building on andextending the OpenCL framework from ITKv4, to speed up the Gaussianpyramid computation and the image resampling step; iii exploitationof certain properties of the B-spline transformation model; ivfurther software optimizations.The accelerated registration tool is employed in a study ondiagnostic classification of Alzheimer's disease and cognitivelynormal controls based on T1-weighted MRI. We selected 299participants from the publicly available Alzheimer's DiseaseNeuroimaging Initiative database. Classification is performed with asupport vector machine based on gray matter volumes as a marker foratrophy. We evaluated two types of strategies (voxel-wise andregion-wise that heavily rely on nonrigid image registration.Parallelization and optimization resulted in an acceleration factorof 4-5x on an 8-core machine. Using OpenCL a speedup factor of ~2was realized for computation of the Gaussian pyramids, and 15-60 forthe resampling step, for larger images. The voxel-wise and theregion-wise classification methods had an area under thereceiver operator characteristic curve of 88% and 90%,respectively, both for standard and accelerated registration.We conclude that the image registration package elastix wassubstantially accelerated, with nearly identical results to thenon-optimized version. The new functionality will become availablein the next release of elastix as open source under the BSD license.
Non-Parametric Bayesian Registration (NParBR) of Body Tumors in DCE-MRI Data.
Pilutti, David; Strumia, Maddalena; Buchert, Martin; Hadjidemetriou, Stathis
2016-04-01
The identification of tumors in the internal organs of chest, abdomen, and pelvis anatomic regions can be performed with the analysis of Dynamic Contrast Enhanced Magnetic Resonance Imaging (DCE-MRI) data. The contrast agent is accumulated differently by pathologic and healthy tissues and that results in a temporally varying contrast in an image series. The internal organs are also subject to potentially extensive movements mainly due to breathing, heart beat, and peristalsis. This contributes to making the analysis of DCE-MRI datasets challenging as well as time consuming. To address this problem we propose a novel pairwise non-rigid registration method with a Non-Parametric Bayesian Registration (NParBR) formulation. The NParBR method uses a Bayesian formulation that assumes a model for the effect of the distortion on the joint intensity statistics, a non-parametric prior for the restored statistics, and also applies a spatial regularization for the estimated registration with Gaussian filtering. A minimally biased intra-dataset atlas is computed for each dataset and used as reference for the registration of the time series. The time series registration method has been tested with 20 datasets of liver, lungs, intestines, and prostate. It has been compared to the B-Splines and to the SyN methods with results that demonstrate that the proposed method improves both accuracy and efficiency. PMID:26672032
周志勇; 李莉华; 郑健; 蒯多杰; 胡粟; 张涛
2014-01-01
基于高斯混合模型(Gaussian mixture model, GMM)的点集非刚性配准算法易受重尾点和异常点影响,提出含局部空间约束的t 分布混合模型的点集非刚性配准算法。通过期望最大化(Expectation maximization, EM)框架将高斯混合模型推广为t分布混合模型；把Dirichlet 分布作为浮动点的先验权重,并构造含局部空间约束性质的Dirichlet 分布参数。使用EM算法获得配准参数的闭合解；计算浮动点的自由度,改变其概率密度分布,避免异常点水平估计误差。实验表明,本文提出的配准算法具有配准误差小、鲁棒性好、抗干扰能力强等优点。%A robust non-rigid registration framework using the student0s-t mixture model with spatial constraints is proposed in this paper. The Gaussian mixture model which is vulnerable to outliers and data longer than normal tails is a special case of the student0s-t mixture model in theory. The Dirichlet distribution is used as a prior distribution to reduce the impact of outliers. The Dirichlet parameter set with spatial constrains is structured to incorporate the spatial information into the decision process. The closed form solution of the parameter set of the student0s-t mixture model is solved by re-parameterizing the student0s-t mixture model in the expectation maximization (EM) algorithm. The degree of freedom of each moving point is calculated to change the probability density to reduce the registration error. It can also avoid estimating the outlier level of data sets that may bring additional error. The experiments showed that this non-rigid registration algorithm has features of high-accuracy and good robustness compared to other point set registration approaches.
Geometry Modeling of Ship Hull Based on Non-uniform B-spline
WANG Hu; ZOU Zao-jian
2008-01-01
In order to generate the three-dimensional (3-D) hull surface accurately and smoothly, a mixed method which is made up of non-uniform B-spline together with an iterative procedure was developed. By using the iterative method the data points on each section curve are calculated and the generalized waterlines and transverse section curves are determined. Then using the non-uniform B-spline expression, the control vertex net of the hull is calculated based on the generalized waterlines and section curves. A ship with tunnel stern was taken as test case. The numerical results prove that the proposed approach for geometry modeling of 3-D ship hull surface is accurate and effective.
B-Spline Active Contour with Handling of Topology Changes for Fast Video Segmentation
Frederic Precioso
2002-06-01
Full Text Available This paper deals with video segmentation for MPEG-4 and MPEG-7 applications. Region-based active contour is a powerful technique for segmentation. However most of these methods are implemented using level sets. Although level-set methods provide accurate segmentation, they suffer from large computational cost. We propose to use a regular B-spline parametric method to provide a fast and accurate segmentation. Our B-spline interpolation is based on a fixed number of points 2j depending on the level of the desired details. Through this spatial multiresolution approach, the computational cost of the segmentation is reduced. We introduce a length penalty. This results in improving both smoothness and accuracy. Then we show some experiments on real-video sequences.
Ship hull plate processing surface fairing with constraints based on B-spline
无
2005-01-01
The problem of ship hull plate processing surface fairing with constraints based on B-spline is solved in this paper. The algorithm for B-spline curve fairing with constraints is one of the most common methods in plane curve fairing. The algorithm can be applied to global and local curve fairing. It can constrain the perturbation range of the control points and the shape variation of the curve, and get a better fairing result in plane curves. In this paper, a new fairing algorithm with constraints for curves and surfaces in space is presented. Then this method is applied to the experiments of ship hull plate processing surface. Finally numerical results are obtained to show the efficiency of this method.
Investigation of confined hydrogen atom in spherical cavity, using B-splines basis set
M Barezi
2011-03-01
Full Text Available Studying confined quantum systems (CQS is very important in nano technology. One of the basic CQS is a hydrogen atom confined in spherical cavity. In this article, eigenenergies and eigenfunctions of hydrogen atom in spherical cavity are calculated, using linear variational method. B-splines are used as basis functions, which can easily construct the trial wave functions with appropriate boundary conditions. The main characteristics of B-spline are its high localization and its flexibility. Besides, these functions have numerical stability and are able to spend high volume of calculation with good accuracy. The energy levels as function of cavity radius are analyzed. To check the validity and efficiency of the proposed method, extensive convergence test of eigenenergies in different cavity sizes has been carried out.
Left ventricular motion reconstruction with a prolate spheroidal B-spline model
Tagged cardiac magnetic resonance (MR) imaging can non-invasively image deformation of the left ventricular (LV) wall. Three-dimensional (3D) analysis of tag data requires fitting a deformation model to tag lines in the image data. In this paper, we present a 3D myocardial displacement and strain reconstruction method based on a B-spline deformation model defined in prolate spheroidal coordinates, which more closely matches the shape of the LV wall than existing Cartesian or cylindrical coordinate models. The prolate spheroidal B-spline (PSB) deformation model also enforces smoothness across and can compute strain at the apex. The PSB reconstruction algorithm was evaluated on a previously published data set to allow head-to-head comparison of the PSB model with existing LV deformation reconstruction methods. We conclude that the PSB method can accurately reconstruct deformation and strain in the LV wall from tagged MR images and has several advantages relative to existing techniques
The application of cubic B-spline collocation method in impact force identification
Qiao, Baijie; Chen, Xuefeng; Xue, Xiaofeng; Luo, Xinjie; Liu, Ruonan
2015-12-01
The accurate real-time characterization of impact event is vital during the life-time of a mechanical product. However, the identified impact force may seriously diverge from the real one due to the unknown noise contaminating the measured data, as well as the ill-conditioned system matrix. In this paper, a regularized cubic B-spline collocation (CBSC) method is developed for identifying the impact force time history, which overcomes the deficiency of the ill-posed problem. The cubic B-spline function by controlling the mesh size of the collocation point has the profile of a typical impact event. The unknown impact force is approximated by a set of translated cubic B-spline functions and then the original governing equation of force identification is reduced to find the coefficient of the basis function at each collocation point. Moreover, a modified regularization parameter selection criterion derived from the generalized cross validation (GCV) criterion for the truncated singular value decomposition (TSVD) is introduced for the CBSC method to determine the optimum regularization number of cubic B-spline functions. In the numerical simulation of a two degrees-of-freedom (DOF) system, the regularized CBSC method is validated under different noise levels and frequency bands of exciting forces. Finally, an impact experiment is performed on a clamped-free shell structure to confirm the performance of the regularized CBSC method. Experimental results demonstrate that the peak relative errors of impact forces based on the regularized CBSC method are below 8%, while those based on the TSVD method are approximately 30%.
An efficient active B-spline/nurbs model for virtual sculpting
Moore, Patricia
2013-01-01
This thesis presents an Efficient Active B-Spline/Nurbs Model for Virtual Sculpting. In spite of the on-going rapid development of computer graphics and computer-aided design tools, 3D graphics designers still rely on non-intuitive modelling procedures for the creation and manipulation of freeform virtual content. The ’Virtual Sculpting' paradigm is a well-established mechanism for shielding designers from the complex mathematics that underpin freeform shape design. The premise is to emulate ...
Electron scattering from krypton: High-resolution experiments and B-spline R-matrix calculations
In a joint experimental and theoretical effort, we carried out a detailed study of e-Kr collisions. For elastic scattering and excitation of the 4p55s states, we present total and angle-differential cross sections over the entire angular range (0°-180°) for a number of energies, as well as energy scans for selected angles. The experimental results are in very satisfactory agreement with predictions from fully relativistic Dirac B-spline R-matrix models.
The Numerical Approach to the Fisher's Equation via Trigonometric Cubic B-spline Collocation Method
Ersoy, Ozlem; Dag, Idris
2016-01-01
In this study, we set up a numerical technique to get approximate solutions of Fisher's equation which is one of the most important model equation in population biology. We integrate the equation fully by using combination of the trigonometric cubic B-spline functions for space variable and Crank-Nicolson for the time integration. Numerical results have been presented to show the accuracy of the current algorithm. We have seen that the proposed technique is a good alternative to some existing...
A mixed basis density functional approach for one-dimensional systems with B-splines
Ren, Chung-Yuan; Chang, Yia-Chung; Hsue, Chen-Shiung
2016-05-01
A mixed basis approach based on density functional theory is extended to one-dimensional (1D) systems. The basis functions here are taken to be the localized B-splines for the two finite non-periodic dimensions and the plane waves for the third periodic direction. This approach will significantly reduce the number of the basis and therefore is computationally efficient for the diagonalization of the Kohn-Sham Hamiltonian. For 1D systems, B-spline polynomials are particularly useful and efficient in two-dimensional spatial integrations involved in the calculations because of their absolute localization. Moreover, B-splines are not associated with atomic positions when the geometry structure is optimized, making the geometry optimization easy to implement. With such a basis set we can directly calculate the total energy of the isolated system instead of using the conventional supercell model with artificial vacuum regions among the replicas along the two non-periodic directions. The spurious Coulomb interaction between the charged defect and its repeated images by the supercell approach for charged systems can also be avoided. A rigorous formalism for the long-range Coulomb potential of both neutral and charged 1D systems under the mixed basis scheme will be derived. To test the present method, we apply it to study the infinite carbon-dimer chain, graphene nanoribbon, carbon nanotube and positively-charged carbon-dimer chain. The resulting electronic structures are presented and discussed in detail.
Design of Low-Pass Digital Differentiators Based on B-splines
Zijun He
2014-07-01
Full Text Available This paper describes a new method for designing low-pass differentiators that could be widely suitable for low-frequency signals with different sampling rates. The method is based on the differential property of convolution and the derivatives of B-spline bias functions. The first order differentiator is just constructed by the first derivative of the B-spline of degree 5 or 4. A high (>2 order low-pass differentiator is constructed by cascading two low order differentiators, of which the coefficients are obtained from the nth derivative of a B-spline of degree n+2 expanded by factor a. In this paper, the properties of the proposed differentiators were presented. In addition, we gave the examples of designing the first to sixth order differentiators, and several simulations, including the effects of the factor a on the results and the anti-noise capability of the proposed differentiators. These properties analysis and simulations indicate that the proposed differentiator can be applied to a wide range of low-frequency signals, and the trade-off between noise- reduction and signal preservation can be made by selecting the maximum allowable value of a.
Geetha M
2012-03-01
Full Text Available Sign language is the most natural way of expression for the deaf community. The urge to support the integration of deaf people into the hearing society made the automatic sign language recognition, an area of interest for the researchers. Indian Sign Language (ISL is a visual-spatial language which provides linguistic information using hands, arms, facial expressions, and head/body postures. In this paper we propose a novel vision-based recognition of Indian Sign Language Alphabets and Numerals using B-Spline Approximation. Gestures of ISL alphabets are complex since it involves the gestures of both the hands together. Our algorithm approximates the boundary extracted from the Region of Interest, to a B-Spline curve by taking the Maximum Curvature Points (MCPs as the Control points. Then the B-Spline curve is subjected to iterations for smoothening resulting in the extraction of Key Maximum Curvature points (KMCPs, which are the key contributors of the gesture shape. Hence a translation & scale invariant feature vector is obtained from the spatial locations of the KMCPs in the 8 Octant Regions of the 2D Space which isgiven for classification.
A B-Spline-Based Colocation Method to Approximate the Solutions to the Equations of Fluid Dynamics
Johnson, Richard Wayne; Landon, Mark Dee
1999-07-01
The potential of a B-spline collocation method for numerically solving the equations of fluid dynamics is discussed. It is known that B-splines can resolve curves with drastically fewer data than can their standard shape function counterparts. This feature promises to allow much faster numerical simulations of fluid flow than standard finite volume/finite element methods without sacrificing accuracy. An example channel flow problem is solved using the method.
A B-Spline-Based Colocation Method to Approximate the Solutions to the Equations of Fluid Dynamics
M. D. Landon; R. W. Johnson
1999-07-01
The potential of a B-spline collocation method for numerically solving the equations of fluid dynamics is discussed. It is known that B-splines can resolve complex curves with drastically fewer data than can their standard shape function counterparts. This feature promises to allow much faster numerical simulations of fluid flow than standard finite volume/finite element methods without sacrificing accuracy. An example channel flow problem is solved using the method.
Shazalina Mat Zin; Ahmad Abd. Majid; Ahmad Izani Md. Ismail; Muhammad Abbas
2014-01-01
The generalized nonlinear Klien-Gordon equation is important in quantum mechanics and related fields. In this paper, a semi-implicit approach based on hybrid cubic B-spline is presented for the approximate solution of the nonlinear Klien-Gordon equation. The usual finite difference approach is used to discretize the time derivative while hybrid cubic B-spline is applied as an interpolating function in the space dimension. The results of applications to several test problems indicate good agre...
Data assimilation using Bayesian filters and B-spline geological models
This paper proposes a new approach to problems of data assimilation, also known as history matching, of oilfield production data by adjustment of the location and sharpness of patterns of geological facies. Traditionally, this problem has been addressed using gradient based approaches with a level set parameterization of the geology. Gradient-based methods are robust, but computationally demanding with real-world reservoir problems and insufficient for reservoir management uncertainty assessment. Recently, the ensemble filter approach has been used to tackle this problem because of its high efficiency from the standpoint of implementation, computational cost, and performance. Incorporation of level set parameterization in this approach could further deal with the lack of differentiability with respect to facies type, but its practical implementation is based on some assumptions that are not easily satisfied in real problems. In this work, we propose to describe the geometry of the permeability field using B-spline curves. This transforms history matching of the discrete facies type to the estimation of continuous B-spline control points. As filtering scheme, we use the ensemble square-root filter (EnSRF). The efficacy of the EnSRF with the B-spline parameterization is investigated through three numerical experiments, in which the reservoir contains a curved channel, a disconnected channel or a 2-dimensional closed feature. It is found that the application of the proposed method to the problem of adjusting facies edges to match production data is relatively straightforward and provides statistical estimates of the distribution of geological facies and of the state of the reservoir.
Data assimilation using Bayesian filters and B-spline geological models
Duan, Lian
2011-04-01
This paper proposes a new approach to problems of data assimilation, also known as history matching, of oilfield production data by adjustment of the location and sharpness of patterns of geological facies. Traditionally, this problem has been addressed using gradient based approaches with a level set parameterization of the geology. Gradient-based methods are robust, but computationally demanding with real-world reservoir problems and insufficient for reservoir management uncertainty assessment. Recently, the ensemble filter approach has been used to tackle this problem because of its high efficiency from the standpoint of implementation, computational cost, and performance. Incorporation of level set parameterization in this approach could further deal with the lack of differentiability with respect to facies type, but its practical implementation is based on some assumptions that are not easily satisfied in real problems. In this work, we propose to describe the geometry of the permeability field using B-spline curves. This transforms history matching of the discrete facies type to the estimation of continuous B-spline control points. As filtering scheme, we use the ensemble square-root filter (EnSRF). The efficacy of the EnSRF with the B-spline parameterization is investigated through three numerical experiments, in which the reservoir contains a curved channel, a disconnected channel or a 2-dimensional closed feature. It is found that the application of the proposed method to the problem of adjusting facies edges to match production data is relatively straightforward and provides statistical estimates of the distribution of geological facies and of the state of the reservoir.
A cubic B-spline Galerkin approach for the numerical simulation of the GEW equation
S. Battal Gazi Karakoç
2016-02-01
Full Text Available The generalized equal width (GEW wave equation is solved numerically by using lumped Galerkin approach with cubic B-spline functions. The proposed numerical scheme is tested by applying two test problems including single solitary wave and interaction of two solitary waves. In order to determine the performance of the algorithm, the error norms L2 and L∞ and the invariants I1, I2 and I3 are calculated. For the linear stability analysis of the numerical algorithm, von Neumann approach is used. As a result, the obtained findings show that the presented numerical scheme is preferable to some recent numerical methods.
B-Spline Filtering for Automatic Detection of Calcification Lesions in Mammograms
Bueno, G.; Sánchez, S.; Ruiz, M.
2006-10-01
Breast cancer continues to be an important health problem between women population. Early detection is the only way to improve breast cancer prognosis and significantly reduce women mortality. It is by using CAD systems that radiologist can improve their ability to detect, and classify lesions in mammograms. In this study the usefulness of using B-spline based on a gradient scheme and compared to wavelet and adaptative filtering has been investigated for calcification lesion detection and as part of CAD systems. The technique has been applied to different density tissues. A qualitative validation shows the success of the method.
A numerical solution of the Burgers' equation using septic B-splines
Ramadan, Mohamed A. [Department of Mathematics, Faculty of Science, Menoufia University, Shiben El-Koom (Egypt)] e-mail: mramadan@mailer.eun.eg; El-Danaf, Talaat S. [Department of Mathematics, Faculty of Science, Menoufia University, Shiben El-Koom (Egypt); Abd Alaal, Faisal E.I. [Department of Mathematics, Faculty of Science, Menoufia University, Shiben El-Koom (Egypt)
2005-11-01
In this paper, numerical solutions of the nonlinear Burgers' equation are obtained by a method based on collocation of septic B-splines over finite elements. Applying the Von-Neumann stability analysis, the proposed method is shown to be unconditionally stable. Numerical solutions of the modified Burgers' equation are also obtained by making a simple change of the suggested numerical scheme for the Burgers' equation. The accuracy of the presented method is demonstrated by two test problems. The numerical results are found to be in good agreement with the exact solutions.
Energy Spectra of the Confined Atoms Obtained by Using B-Splines
SHI Ting-Yun; BAO Cheng-Guang; LI Bai-Wen
2001-01-01
We have calculated the energy spectra of one- and two-electron atoms (ions) centered in an impenetrable spherical box by variational method with B-splines as basis functions. Accurate results are obtained for both large and small radii of confinement. The critical box radius of confined hydrogen atom is also calculated to show the usefulness of our method. A partial energy degeneracy in confined hydrogen atom is found when the radius of spherical box is equal to the distance at which a node of single-node wavefunctions of free hydrogen atom is located.
Simulating the focusing of light onto 1D nanostructures with a B-spline modal method
Bouchon, P.; Chevalier, P.; Héron, S.; Pardo, F.; Pelouard, J.-L.; Haïdar, R.
2015-03-01
Focusing the light onto nanostructures thanks to spherical lenses is a first step to enhance the field, and is widely used in applications, in particular for enhancing non-linear effects like the second harmonic generation. Nonetheless, the electromagnetic response of such nanostructures, which have subwavelength patterns, to a focused beam can not be described by the simple ray tracing formalism. Here, we present a method to compute the response to a focused beam, based on the B-spline modal method. The simulation of a gaussian focused beam is obtained thanks to a truncated decomposition on plane waves computed on a single period, which limits the computation burden.
Complex wavenumber Fourier analysis of the B-spline based finite element method
Kolman, Radek; Plešek, Jiří; Okrouhlík, Miloslav
2014-01-01
Roč. 51, č. 2 (2014), s. 348-359. ISSN 0165-2125 R&D Projects: GA ČR(CZ) GAP101/11/0288; GA ČR(CZ) GAP101/12/2315; GA ČR GPP101/10/P376; GA ČR GA101/09/1630 Institutional support: RVO:61388998 Keywords : elastic wave propagation * dispersion errors * B-spline * finite element method * isogeometric analysis Subject RIV: JR - Other Machinery Impact factor: 1.513, year: 2014 http://www.sciencedirect.com/science/article/pii/S0165212513001479
Numerical solution of elastic wave propagation problems by B-spline finite element method
Kolman, Radek; Plešek, Jiří; Okrouhlík, Miloslav; Gabriel, Dušan
Aviero : Universidade de Aveiro, 2012, s. 1-10. ISBN 978-972-99784-2-5. [ECCOMAS Young Investigators Conference. Aveiro (PT), 24.04.2012-27.04.2012] R&D Projects: GA ČR GPP101/10/P376; GA ČR GA101/09/1630; GA ČR(CZ) GAP101/11/0288; GA ČR(CZ) GAP101/12/2315 Institutional research plan: CEZ:AV0Z20760514 Keywords : elastic wave propagation * B-spline finite element method * spurious oscillations Subject RIV: BI - Acoustics
Non-Stationary Hydrologic Frequency Analysis using B-Splines Quantile Regression
Nasri, B.; St-Hilaire, A.; Bouezmarni, T.; Ouarda, T.
2015-12-01
Hydrologic frequency analysis is commonly used by engineers and hydrologists to provide the basic information on planning, design and management of hydraulic structures and water resources system under the assumption of stationarity. However, with increasing evidence of changing climate, it is possible that the assumption of stationarity would no longer be valid and the results of conventional analysis would become questionable. In this study, we consider a framework for frequency analysis of extreme flows based on B-Splines quantile regression, which allows to model non-stationary data that have a dependence on covariates. Such covariates may have linear or nonlinear dependence. A Markov Chain Monte Carlo (MCMC) algorithm is used to estimate quantiles and their posterior distributions. A coefficient of determination for quantiles regression is proposed to evaluate the estimation of the proposed model for each quantile level. The method is applied on annual maximum and minimum streamflow records in Ontario, Canada. Climate indices are considered to describe the non-stationarity in these variables and to estimate the quantiles in this case. The results show large differences between the non-stationary quantiles and their stationary equivalents for annual maximum and minimum discharge with high annual non-exceedance probabilities. Keywords: Quantile regression, B-Splines functions, MCMC, Streamflow, Climate indices, non-stationarity.
DBSR_HF: A B-spline Dirac-Hartree-Fock program
Zatsarinny, Oleg; Froese Fischer, Charlotte
2016-05-01
A B-spline version of a general Dirac-Hartree-Fock program is described. The usual differential equations are replaced by a set of generalized eigenvalue problems of the form (Ha -εa B) Pa = 0, where Ha and B are the Hamiltonian and overlap matrices, respectively, and Pa is the two-component relativistic orbit in the B-spline basis. A default universal grid allows for flexible adjustment to different nuclear models. When two orthogonal orbitals are both varied, the energy must also be stationary with respect to orthonormal transformations. At such a stationary point the off-diagonal Lagrange multipliers may be eliminated through projection operators. The self-consistent field procedure exhibits excellent convergence. Several atomic states can be considered simultaneously, including some configuration-interaction calculations. The program provides several options for the treatment of Breit interaction and QED corrections. The information about atoms up to Z = 104 is stored by the program. Along with a simple interface through command-line arguments, this information allows the user to run the program with minimal initial preparations.
A spectral/B-spline method for the Navier-Stokes equations in unbounded domains
The numerical method presented in this paper aims at solving the incompressible Navier-Stokes equations in unbounded domains. The problem is formulated in cylindrical coordinates and the method is based on a Galerkin approximation scheme that makes use of vector expansions that exactly satisfy the continuity constraint. More specifically, the divergence-free basis vector functions are constructed with Fourier expansions in the θ and z directions while mapped B-splines are used in the semi-infinite radial direction. Special care has been taken to account for the particular analytical behaviors at both end points r=0 and r→∞. A modal reduction algorithm has also been implemented in the azimuthal direction, allowing for a relaxation of the CFL constraint on the timestep size and a possibly significant reduction of the number of DOF. The time marching is carried out using a mixed quasi-third order scheme. Besides the advantages of a divergence-free formulation and a quasi-spectral convergence, the local character of the B-splines allows for a great flexibility in node positioning while keeping narrow bandwidth matrices. Numerical tests show that the present method compares advantageously with other similar methodologies using purely global expansions
A critical challenges in urban aeras is slums. In fact, they are considered a source of crime and disease due to poor-quality housing, unsanitary conditions, poor infrastructures and occupancy security. The poor in the dense urban slums are the most vulnerable to infection due to (i) inadequate and restricted access to safety, drinking water and sufficient quantities of water for personal hygiene; (ii) the lack of removal and treatment of excreta; and (iii) the lack of removal of solid waste. This study aims to investigate the capability of ENVISAT ASAR satellite and Google Earth data for three-dimensional (3-D) slum urban reconstruction in developed countries such as Egypt. The main objective of this work is to utilize some 3-D automatic detection algorithm for urban slum in ENVISAT ASAR and Google Erath images were acquired in Cairo, Egypt using Fuzzy B-spline algorithm. The results show that the fuzzy algorithm is the best indicator for chaotic urban slum as it can discriminate between them from its surrounding environment. The combination of Fuzzy and B-spline then used to reconstruct 3-D of urban slum. The results show that urban slums, road network, and infrastructures are perfectly discriminated. It can therefore be concluded that the fuzzy algorithm is an appropriate algorithm for chaotic urban slum automatic detection in ENVSIAT ASAR and Google Earth data
Marghany, Maged
2014-06-01
A critical challenges in urban aeras is slums. In fact, they are considered a source of crime and disease due to poor-quality housing, unsanitary conditions, poor infrastructures and occupancy security. The poor in the dense urban slums are the most vulnerable to infection due to (i) inadequate and restricted access to safety, drinking water and sufficient quantities of water for personal hygiene; (ii) the lack of removal and treatment of excreta; and (iii) the lack of removal of solid waste. This study aims to investigate the capability of ENVISAT ASAR satellite and Google Earth data for three-dimensional (3-D) slum urban reconstruction in developed countries such as Egypt. The main objective of this work is to utilize some 3-D automatic detection algorithm for urban slum in ENVISAT ASAR and Google Erath images were acquired in Cairo, Egypt using Fuzzy B-spline algorithm. The results show that the fuzzy algorithm is the best indicator for chaotic urban slum as it can discriminate between them from its surrounding environment. The combination of Fuzzy and B-spline then used to reconstruct 3-D of urban slum. The results show that urban slums, road network, and infrastructures are perfectly discriminated. It can therefore be concluded that the fuzzy algorithm is an appropriate algorithm for chaotic urban slum automatic detection in ENVSIAT ASAR and Google Earth data.
B-splines as a Tool to Solve Constraints in Non-Hydrostatic Forecast Model
Subias, Alvaro
2016-01-01
Finite elements has been proven to be an useful tool to discretize the vertical coordinate in the hydrostatic forecast models allowing to define model variables in full levels so that no staggering is needed. In the non-hydrostatic case a constraint in the vertical operators appears (called C1) that does not allow to reduce the set of semi-implicit linear equations to a single equation in one variable as in the analytic case. Recently vertical finite elements based in B-splines have been used with an iterative method to relax the C1 constraint. In this paper we want to develop properly some representations of vertical operators in terms of B-splines in order to keep the C1-constraint. An invertibility relation between integral and derivative operators between vertical velocity and vertical divergence is also presented. The final scope of this paper is to provide a theoretical framework of development of finite element vertical operators to be implemented in the nh-Harmonie model
Railroad inspection based on ACFM employing a non-uniform B-spline approach
Chacón Muñoz, J. M.; García Márquez, F. P.; Papaelias, M.
2013-11-01
The stresses sustained by rails have increased in recent years due to the use of higher train speeds and heavier axle loads. For this reason surface and near-surface defects generate by Rolling Contact Fatigue (RCF) have become particularly significant as they can cause unexpected structural failure of the rail, resulting in severe derailments. The accident that took place in Hatfield, UK (2000), is an example of a derailment caused by the structural failure of a rail section due to RCF. Early detection of RCF rail defects is therefore of paramount importance to the rail industry. The performance of existing ultrasonic and magnetic flux leakage techniques in detecting rail surface-breaking defects, such as head checks and gauge corner cracking, is inadequate during high-speed inspection, while eddy current sensors suffer from lift-off effects. The results obtained through rail inspection experiments under simulated conditions using Alternating Current Field Measurement (ACFM) probes, suggest that this technique can be applied for the accurate and reliable detection of surface-breaking defects at high inspection speeds. This paper presents the B-Spline approach used for the accurate filtering the noise of the raw ACFM signal obtained during high speed tests to improve the reliability of the measurements. A non-uniform B-spline approximation is employed to calculate the exact positions and the dimensions of the defects. This method generates a smooth approximation similar to the ACFM dataset points related to the rail surface-breaking defect.
D. Reclik
2008-08-01
Full Text Available Purpose: The main reason of this paper was to prepare the system, which tests the use of elastic band for smoothing the collision-free trajectory. The aided robot off-line programming system is based on NURBS and B-Spline curves. Because there is a lot of information in references about using elastic band algorithm, authors decided to compare these two methods. The most important criterion in robotics is having the smoothest possible robot trajectory, so as a standard there the NURBS curves (C2 smooth class were used.Design/methodology/approach: Pascal language compiler was used for research. All algorithms were coded in this programming language and compiled. Results were set in Microsoft Excel worksheet.Findings: Results show that calculations, which were made with B-Spline method, have taken less time than calculations based on elastic band curves. Moreover, the elastic band method gave the smoothest curves but only in geometrical sense, which is less important (the first and second derivate are not continuous, which is the most important issue in presented case. That is why it was found that using the B-Spline algorithm is a better solution, because it takes less time and gives better quality results.Research limitations/implications: The MS Windows application was created, which generates smooth curves (in geometrical sense by marking the interpolation base points which are calculated by the collision-free movement planner. This application generates curves by using both presented methods - B-Spline and elastic band. Both of these curves were compared in regard of standard deviation and variance of B-Spline and elastic band.Practical implications: Because the elastic band algorithm takes a lot of time (three times longer than B-Spline it is not used in the final application. The authors used B-Spline method to make smoother and optimized trajectory in application for off-line collision-free robot programming.Originality/value: This is a new
Respiratory motion in emission tomography leads to reduced image quality. Developed correction methodology has been concentrating on the use of respiratory synchronized acquisitions leading to gated frames. Such frames, however, are of low signal-to-noise ratio as a result of containing reduced statistics. In this work, we describe the implementation of an elastic transformation within a list-mode-based reconstruction for the correction of respiratory motion over the thorax, allowing the use of all data available throughout a respiratory motion average acquisition. The developed algorithm was evaluated using datasets of the NCAT phantom generated at different points throughout the respiratory cycle. List-mode-data-based PET-simulated frames were subsequently produced by combining the NCAT datasets with Monte Carlo simulation. A non-rigid registration algorithm based on B-spline basis functions was employed to derive transformation parameters accounting for the respiratory motion using the NCAT dynamic CT images. The displacement matrices derived were subsequently applied during the image reconstruction of the original emission list mode data. Two different implementations for the incorporation of the elastic transformations within the one-pass list mode EM (OPL-EM) algorithm were developed and evaluated. The corrected images were compared with those produced using an affine transformation of list mode data prior to reconstruction, as well as with uncorrected respiratory motion average images. Results demonstrate that although both correction techniques considered lead to significant improvements in accounting for respiratory motion artefacts in the lung fields, the elastic-transformation-based correction leads to a more uniform improvement across the lungs for different lesion sizes and locations
Gu, Renliang, E-mail: Venliang@iastate.edu, E-mail: ald@iastate.edu; Dogandžić, Aleksandar, E-mail: Venliang@iastate.edu, E-mail: ald@iastate.edu [Iowa State University, Center for Nondestructive Evaluation, 1915 Scholl Road, Ames, IA 50011 (United States)
2015-03-31
We develop a sparse image reconstruction method for polychromatic computed tomography (CT) measurements under the blind scenario where the material of the inspected object and the incident energy spectrum are unknown. To obtain a parsimonious measurement model parameterization, we first rewrite the measurement equation using our mass-attenuation parameterization, which has the Laplace integral form. The unknown mass-attenuation spectrum is expanded into basis functions using a B-spline basis of order one. We develop a block coordinate-descent algorithm for constrained minimization of a penalized negative log-likelihood function, where constraints and penalty terms ensure nonnegativity of the spline coefficients and sparsity of the density map image in the wavelet domain. This algorithm alternates between a Nesterov’s proximal-gradient step for estimating the density map image and an active-set step for estimating the incident spectrum parameters. Numerical simulations demonstrate the performance of the proposed scheme.
One Fairing Method of Cubic B-spline Curves Based on Weighted Progressive Iterative Approximation
ZHANG Li; YANG Yan; LI Yuan-yuan; TAN Jie-qing
2014-01-01
A new method to the problem of fairing planar cubic B-spline curves is introduced in this paper. The method is based on weighted progressive iterative approximation (WPIA for short) and consists of following steps:finding the bad point which needs to fair, deleting the bad point, re-inserting a new data point to keep the structure of the curve and applying WPIA method with the new set of the data points to obtain the faired curve. The new set of the data points is formed by the rest of the original data points and the new inserted point. The method can be used for shape design and data processing. Numerical examples are provided to demonstrate the effectiveness of the method.
Ionospheric scintillation modeling for high- and mid-latitude using B-spline technique
Priyadarshi, S.
2015-09-01
Ionospheric scintillation is a significant component of space-weather studies and serves as an estimate for the level of perturbation in the satellite radio wave signal caused due to small-scale ionospheric irregularities. B-spline functions are used on the GPS ground based data collected during the year 2007-2012 for modeling high- and mid-latitude ionospheric scintillation. Proposed model is for Hornsund, Svalbard and Warsaw, Poland. The input data used in this model were recorded by GSV 4004b receivers. For validation, results of this model are compared with the observation and other existing models. Physical behavior of the ionospheric scintillation during different seasons and geomagnetic conditions are discussed well. Model is found in good coherence with the ionospheric scintillation theory as well as to the accepted scintillation mechanism for high- and mid-latitude.
We develop a sparse image reconstruction method for polychromatic computed tomography (CT) measurements under the blind scenario where the material of the inspected object and the incident energy spectrum are unknown. To obtain a parsimonious measurement model parameterization, we first rewrite the measurement equation using our mass-attenuation parameterization, which has the Laplace integral form. The unknown mass-attenuation spectrum is expanded into basis functions using a B-spline basis of order one. We develop a block coordinate-descent algorithm for constrained minimization of a penalized negative log-likelihood function, where constraints and penalty terms ensure nonnegativity of the spline coefficients and sparsity of the density map image in the wavelet domain. This algorithm alternates between a Nesterov’s proximal-gradient step for estimating the density map image and an active-set step for estimating the incident spectrum parameters. Numerical simulations demonstrate the performance of the proposed scheme
B-spline parameterization of spatial response in a monolithic scintillation camera
Solovov, V; Chepel, V; Domingos, V; Martins, R
2016-01-01
A framework for parameterization of the light response functions (LRFs) in a scintillation camera was developed. It is based on approximation of the measured or simulated photosensor response with weighted sums of uniform cubic B-splines or their tensor products. The LRFs represented in this way are smooth, computationally inexpensive to evaluate and require much less memory than non-parametric alternatives. The parameters are found in a straightforward way by the linear least squares method. The use of linear fit makes the fitting process stable and predictable enough to be used in non-supervised mode. Several techniques that allow to reduce the storage and processing power requirements were developed. A software library for fitting simulated and measured light response with spline functions was developed and integrated into an open source software package ANTS2 designed for simulation and data processing for Anger camera-type detectors.
3D Profile Filter Algorithm Based on Parallel Generalized B-spline Approximating Gaussian
REN Zhiying; GAO Chenghui; SHEN Ding
2015-01-01
Currently, the approximation methods of the Gaussian filter by some other spline filters have been developed. However, these methods are only suitable for the study of one-dimensional filtering, when these methods are used for three-dimensional filtering, it is found that a rounding error and quantization error would be passed to the next in every part. In this paper, a new and high-precision implementation approach for Gaussian filter is described, which is suitable for three-dimensional reference filtering. Based on the theory of generalized B-spline function and the variational principle, the transmission characteristics of a digital filter can be changed through the sensitivity of the parameters (t1, t2), and which can also reduce the rounding error and quantization error by the filter in a parallel form instead of the cascade form. Finally, the approximation filter of Gaussian filter is obtained. In order to verify the feasibility of the new algorithm, the reference extraction of the conventional methods are also used and compared. The experiments are conducted on the measured optical surface, and the results show that the total calculation by the new algorithm only requires 0.07 s for 480´480 data points;the amplitude deviation between the reference of the parallel form filter and the Gaussian filter is smaller;the new method is closer to the characteristic of the Gaussian filter through the analysis of three-dimensional roughness parameters, comparing with the cascade generalized B-spline approximating Gaussian. So the new algorithm is also efficient and accurate for the implementation of Gaussian filter in the application of surface roughness measurement.
We report a fully ab initio implementation of exterior complex scaling in B-splines to evaluate total, singly and triply differential cross sections in double photoionization problems. Results for He and H2 double photoionization are presented and compared with experiment
For elastic scattering and electron impact excitation of the 4p5s states in Kr, we present independently normalized, absolute angle-differential cross sections over the entire angular range (0° − 180°). Excellent agreement is obtained between the present experimental data and theoretical predictions from a fully relativistic B-spline R-matrix (close-coupling) model.
Lu, Guolan; Halig, Luma; Wang, Dongsheng; Chen, Zhuo G.; Fei, Baowei
2014-03-01
The determination of tumor margins during surgical resection remains a challenging task. A complete removal of malignant tissue and conservation of healthy tissue is important for the preservation of organ function, patient satisfaction, and quality of life. Visual inspection and palpation is not sufficient for discriminating between malignant and normal tissue types. Hyperspectral imaging (HSI) technology has the potential to noninvasively delineate surgical tumor margin and can be used as an intra-operative visual aid tool. Since histological images provide the ground truth of cancer margins, it is necessary to warp the cancer regions in ex vivo histological images back to in vivo hyperspectral images in order to validate the tumor margins detected by HSI and to optimize the imaging parameters. In this paper, principal component analysis (PCA) is utilized to extract the principle component bands of the HSI images, which is then used to register HSI images with the corresponding histological image. Affine registration is chosen to model the global transformation. A B-spline free form deformation (FFD) method is used to model the local non-rigid deformation. Registration experiment was performed on animal hyperspectral and histological images. Experimental results from animals demonstrated the feasibility of the hyperspectral imaging method for cancer margin detection.
Accurate modeling of system response and scatter distribution is crucial for image reconstruction in emission tomography. Monte Carlo simulations are very well suited to calculate these quantities. However, Monte Carlo simulations are also slow and many simulated counts are needed to provide a sufficiently exact estimate of the detection probabilities. In order to overcome these problems, we propose to split the simulation into two parts, the detection system and the object to be imaged (the patient). A so-called 'virtual boundary' that separates these two parts is introduced. Within the patient, particles are simulated conventionally. Whenever a photon reaches the virtual boundary, its detection probability is calculated analytically by evaluating a multi-dimensional B-spline that depends on the photon position, direction and energy. The unknown B-spline knot values that define this B-spline are fixed by a prior 'pre-' simulation that needs to be run once for each scanner type. After this pre-simulation, the B-spline model can be used in any subsequent simulation with different patients. We show that this approach yields accurate results when simulating the Biograph 16 HiREZ PET scanner with Geant4 Application for Emission Tomography (GATE). The execution time is reduced by a factor of about 22 x (scanner with voxelized phantom) to 30 x (empty scanner) with respect to conventional GATE simulations of same statistical uncertainty. The pre-simulation and calculation of the B-spline knots values could be performed within half a day on a medium-sized cluster.
Hardy, David J; Wolff, Matthew A; Xia, Jianlin; Schulten, Klaus; Skeel, Robert D
2016-03-21
The multilevel summation method for calculating electrostatic interactions in molecular dynamics simulations constructs an approximation to a pairwise interaction kernel and its gradient, which can be evaluated at a cost that scales linearly with the number of atoms. The method smoothly splits the kernel into a sum of partial kernels of increasing range and decreasing variability with the longer-range parts interpolated from grids of increasing coarseness. Multilevel summation is especially appropriate in the context of dynamics and minimization, because it can produce continuous gradients. This article explores the use of B-splines to increase the accuracy of the multilevel summation method (for nonperiodic boundaries) without incurring additional computation other than a preprocessing step (whose cost also scales linearly). To obtain accurate results efficiently involves technical difficulties, which are overcome by a novel preprocessing algorithm. Numerical experiments demonstrate that the resulting method offers substantial improvements in accuracy and that its performance is competitive with an implementation of the fast multipole method in general and markedly better for Hamiltonian formulations of molecular dynamics. The improvement is great enough to establish multilevel summation as a serious contender for calculating pairwise interactions in molecular dynamics simulations. In particular, the method appears to be uniquely capable for molecular dynamics in two situations, nonperiodic boundary conditions and massively parallel computation, where the fast Fourier transform employed in the particle-mesh Ewald method falls short. PMID:27004867
CT segmentation of dental shapes by anatomy-driven reformation imaging and B-spline modelling.
Barone, S; Paoli, A; Razionale, A V
2016-06-01
Dedicated imaging methods are among the most important tools of modern computer-aided medical applications. In the last few years, cone beam computed tomography (CBCT) has gained popularity in digital dentistry for 3D imaging of jawbones and teeth. However, the anatomy of a maxillofacial region complicates the assessment of tooth geometry and anatomical location when using standard orthogonal views of the CT data set. In particular, a tooth is defined by a sub-region, which cannot be easily separated from surrounding tissues by only considering pixel grey-intensity values. For this reason, an image enhancement is usually necessary in order to properly segment tooth geometries. In this paper, an anatomy-driven methodology to reconstruct individual 3D tooth anatomies by processing CBCT data is presented. The main concept is to generate a small set of multi-planar reformation images along significant views for each target tooth, driven by the individual anatomical geometry of a specific patient. The reformation images greatly enhance the clearness of the target tooth contours. A set of meaningful 2D tooth contours is extracted and used to automatically model the overall 3D tooth shape through a B-spline representation. The effectiveness of the methodology has been verified by comparing some anatomy-driven reconstructions of anterior and premolar teeth with those obtained by using standard tooth segmentation tools. Copyright © 2015 John Wiley & Sons, Ltd. PMID:26418417
Analysis of Gene Coexpression by B-Spline Based CoD Estimation
Li Huai
2007-01-01
Full Text Available The gene coexpression study has emerged as a novel holistic approach for microarray data analysis. Different indices have been used in exploring coexpression relationship, but each is associated with certain pitfalls. The Pearson's correlation coefficient, for example, is not capable of uncovering nonlinear pattern and directionality of coexpression. Mutual information can detect nonlinearity but fails to show directionality. The coefficient of determination (CoD is unique in exploring different patterns of gene coexpression, but so far only applied to discrete data and the conversion of continuous microarray data to the discrete format could lead to information loss. Here, we proposed an effective algorithm, CoexPro, for gene coexpression analysis. The new algorithm is based on B-spline approximation of coexpression between a pair of genes, followed by CoD estimation. The algorithm was justified by simulation studies and by functional semantic similarity analysis. The proposed algorithm is capable of uncovering both linear and a specific class of nonlinear relationships from continuous microarray data. It can also provide suggestions for possible directionality of coexpression to the researchers. The new algorithm presents a novel model for gene coexpression and will be a valuable tool for a variety of gene expression and network studies. The application of the algorithm was demonstrated by an analysis on ligand-receptor coexpression in cancerous and noncancerous cells. The software implementing the algorithm is available upon request to the authors.
A baseline correction algorithm for Raman spectroscopy by adaptive knots B-spline
Wang, Xin; Fan, Xian-guang; Xu, Ying-jie; Wang, Xiu-fen; He, Hao; Zuo, Yong
2015-11-01
The Raman spectroscopy technique is a powerful and non-invasive technique for molecular fingerprint detection which has been widely used in many areas, such as food safety, drug safety, and environmental testing. But Raman signals can be easily corrupted by a fluorescent background, therefore we presented a baseline correction algorithm to suppress the fluorescent background in this paper. In this algorithm, the background of the Raman signal was suppressed by fitting a curve called a baseline using a cyclic approximation method. Instead of the traditional polynomial fitting, we used the B-spline as the fitting algorithm due to its advantages of low-order and smoothness, which can avoid under-fitting and over-fitting effectively. In addition, we also presented an automatic adaptive knot generation method to replace traditional uniform knots. This algorithm can obtain the desired performance for most Raman spectra with varying baselines without any user input or preprocessing step. In the simulation, three kinds of fluorescent background lines were introduced to test the effectiveness of the proposed method. We showed that two real Raman spectra (parathion-methyl and colza oil) can be detected and their baselines were also corrected by the proposed method.
Hardy, David J.; Wolff, Matthew A.; Xia, Jianlin; Schulten, Klaus; Skeel, Robert D.
2016-03-01
The multilevel summation method for calculating electrostatic interactions in molecular dynamics simulations constructs an approximation to a pairwise interaction kernel and its gradient, which can be evaluated at a cost that scales linearly with the number of atoms. The method smoothly splits the kernel into a sum of partial kernels of increasing range and decreasing variability with the longer-range parts interpolated from grids of increasing coarseness. Multilevel summation is especially appropriate in the context of dynamics and minimization, because it can produce continuous gradients. This article explores the use of B-splines to increase the accuracy of the multilevel summation method (for nonperiodic boundaries) without incurring additional computation other than a preprocessing step (whose cost also scales linearly). To obtain accurate results efficiently involves technical difficulties, which are overcome by a novel preprocessing algorithm. Numerical experiments demonstrate that the resulting method offers substantial improvements in accuracy and that its performance is competitive with an implementation of the fast multipole method in general and markedly better for Hamiltonian formulations of molecular dynamics. The improvement is great enough to establish multilevel summation as a serious contender for calculating pairwise interactions in molecular dynamics simulations. In particular, the method appears to be uniquely capable for molecular dynamics in two situations, nonperiodic boundary conditions and massively parallel computation, where the fast Fourier transform employed in the particle-mesh Ewald method falls short.
Merging quantum-chemistry with B-splines to describe molecular photoionization
Argenti, L.; Marante, C.; Klinker, M.; Corral, I.; Gonzalez, J.; Martin, F.
2016-05-01
Theoretical description of observables in attosecond pump-probe experiments requires a good representation of the system's ionization continuum. For polyelectronic atoms and molecules, however, this is still a challenge, due to the complicated short-range structure of correlated electronic wavefunctions. Whereas quantum chemistry packages (QCP) implementing sophisticated methods to compute bound electronic molecular states are well established, comparable tools for the continuum are not widely available yet. To tackle this problem, we have developed a new approach that, by means of a hybrid Gaussian-B-spline basis, interfaces existing QCPs with close-coupling scattering methods. To illustrate the viability of this approach, we report results for the multichannel ionization of the helium atom and of the hydrogen molecule that are in excellent agreement with existing accurate benchmarks. These findings, together with the flexibility of QCPs, make of this approach a good candidate for the theoretical study of the ionization of poly-electronic systems. FP7/ERC Grant XCHEM 290853.
Based on the multiquadric trigonometric B-spline quasi-interpolant, this paper proposes a meshless scheme for some partial differential equations whose solutions are periodic with respect to the spatial variable. This scheme takes into account the periodicity of the analytic solution by using derivatives of a periodic quasi-interpolant (multiquadric trigonometric B-spline quasi-interpolant) to approximate the spatial derivatives of the equations. Thus, it overcomes the difficulties of the previous schemes based on quasi-interpolation (requiring some additional boundary conditions and yielding unwanted high-order discontinuous points at the boundaries in the spatial domain). Moreover, the scheme also overcomes the difficulty of the meshless collocation methods (i.e., yielding a notorious ill-conditioned linear system of equations for large collocation points). The numerical examples that are presented at the end of the paper show that the scheme provides excellent approximations to the analytic solutions. (general)
Edge detection of molten pool and weld line for CO2 welding based on B-spline wavelet
薛家祥; 贾林; 李海宝; 张丽玲
2004-01-01
Due to the disturbances of spatters, dusts and strong arc light, it is difficult to detect the molten pool edge and the weld line location in CO2 welding processes. The median filtering and self-multiplication was employed to preprocess the image of the CO2 welding in order to detect effectively the edge of molten pool and the location of weld line. The B-spline wavelet algorithm has been investigated, the influence of different scales and thresholds on the results of the edge detection have been compared and analyzed. The experimental results show that better performance to extract the edge of the molten pool and the location of weld line can be obtained by using the B-spline wavelet transform. The proposed edge detection approach can be further applied to the control of molten depth and the seam tracking.
T Nikbakht
2012-12-01
Full Text Available Effects of quantum size and potential shape on the spectra of an electron and a hydrogenic-donor at the center of a permeable spherical cavity have been calculated, using linear variational method. B-splines have been used as basis functions. By extensive convergence tests and comparing with other results given in the literature, the validity and efficiency of the method were confirmed.
Free Form Deformation–Based Image Registration Improves Accuracy of Traction Force Microscopy
Jorge-Peñas, Alvaro; Izquierdo-Alvarez, Alicia; Aguilar-Cuenca, Rocio; Vicente-Manzanares, Miguel; Garcia-Aznar, José Manuel; Van Oosterwyck, Hans; de-Juan-Pardo, Elena M.; Ortiz-de-Solorzano, Carlos; Muñoz-Barrutia, Arrate
2015-01-01
Traction Force Microscopy (TFM) is a widespread method used to recover cellular tractions from the deformation that they cause in their surrounding substrate. Particle Image Velocimetry (PIV) is commonly used to quantify the substrate’s deformations, due to its simplicity and efficiency. However, PIV relies on a block-matching scheme that easily underestimates the deformations. This is especially relevant in the case of large, locally non-uniform deformations as those usually found in the vicinity of a cell’s adhesions to the substrate. To overcome these limitations, we formulate the calculation of the deformation of the substrate in TFM as a non-rigid image registration process that warps the image of the unstressed material to match the image of the stressed one. In particular, we propose to use a B-spline -based Free Form Deformation (FFD) algorithm that uses a connected deformable mesh to model a wide range of flexible deformations caused by cellular tractions. Our FFD approach is validated in 3D fields using synthetic (simulated) data as well as with experimental data obtained using isolated endothelial cells lying on a deformable, polyacrylamide substrate. Our results show that FFD outperforms PIV providing a deformation field that allows a better recovery of the magnitude and orientation of tractions. Together, these results demonstrate the added value of the FFD algorithm for improving the accuracy of traction recovery. PMID:26641883
Efficient Constrained Local Model Fitting for Non-Rigid Face Alignment
Lucey, Simon; Wang, Yang; Cox, Mark; Sridharan, Sridha; Cohn, Jeffery F.
2009-01-01
Active appearance models (AAMs) have demonstrated great utility when being employed for non-rigid face alignment/tracking. The “simultaneous” algorithm for fitting an AAM achieves good non-rigid face registration performance, but has poor real time performance (2-3 fps). The “project-out” algorithm for fitting an AAM achieves faster than real time performance (> 200 fps) but suffers from poor generic alignment performance. In this paper we introduce an extension to a discriminative method for...
Kolman, Radek; Kopačka, Ján; Plešek, Jiří; Okrouhlík, Miloslav; Gabriel, Dušan
Stockholm : KTH Mechanics : The Royal Institute of Technology, 2010 - (Eriksson, A.; Tibert, G.), s. 255-258 ISSN 0348-467X. [Nordic Seminar on Computational Mechanics /23./. Stockholm (SE), 21.10.2010-22.10.2010] R&D Projects: GA ČR(CZ) GPP101/10/P376; GA ČR(CZ) GA101/09/1630; GA ČR(CZ) GA101/07/1471 Institutional research plan: CEZ:AV0Z20760514 Keywords : B-spline * based finite element method * dispersion analysis * wave propagation Subject RIV: BI - Acoustics http://www.mech.kth.se
The 3D tomography reconstruction has been a profitable alternative in the analysis of the FCC-type- riser (Fluid Catalytic Cracking), for appropriately keeping track of the sectional catalyst concentration distribution in the process of oil refining. The method of tomography reconstruction proposed by M. Azzi and colleagues (1991) uses a relatively small amount of trajectories (from 3 to 5) and projections (from 5 to 7) of gamma rays, a desirable feature in the industrial process tomography. Compared to more popular methods, such as the FBP (Filtered Back Projection), which demands a much higher amount of gamma rays projections, the method by Azzi et al. is more appropriate for the industrial process, where the physical limitations and the cost of the process require more economical arrangements. The use of few projections and trajectories facilitates the diagnosis in the flow dynamical process. This article proposes an improvement in the basis functions introduced by Azzi et al., through the use of quadratic B-splines functions. The use of B-splines functions makes possible a smoother surface reconstruction of the density distribution, since the functions are continuous and smooth. This work describes how the modeling can be done. (author)
Joong-Hyun Rhim; Doo-Yeoun Cho; Kyu-Yeul Lee; Tae-Wan Kim
2006-01-01
We propose a method that automatically generates discrete bicubic G1 continuous B-spline surfaces that interpolate the curve network of a ship hullform. First, the curves in the network are classified into two types: boundary curves and "reference curves". The boundary curves correspond to a set of rectangular (or triangular) topological type that can be represented with tensor-product (or degenerate) B-spline surface patches. Next, in the interior of the patches,surface fitting points and cross boundary derivatives are estimated from the reference curves by constructing "virtual" isoparametric curves. Finally, a discrete G1 continuous B-spline surface is generated by a surface fitting algorithm. Several smooth ship hullform surfaces generated from curve networks corresponding to actual ship hullforms demonstrate the quality of the method.
I Nyoman Budiantara
2006-01-01
Full Text Available Regression analysis is constructed for capturing the influences of independent variables to dependent ones. It can be done by looking at the relationship between those variables. This task of approximating the mean function can be done essentially in two ways. The quiet often use parametric approach is to assume that the mean curve has some prespecified functional forms. Alternatively, nonparametric approach, .i.e., without reference to a specific form, is used when there is no information of the regression function form (Haerdle, 1990. Therefore nonparametric approach has more flexibilities than the parametric one. The aim of this research is to find the best fit model that captures relationship between admission test score to the GPA. This particular data was taken from the Department of Design Communication and Visual, Petra Christian University, Surabaya for year 1999. Those two approaches were used here. In the parametric approach, we use simple linear, quadric cubic regression, and in the nonparametric ones, we use B-Spline and Multivariate Adaptive Regression Splines (MARS. Overall, the best model was chosen based on the maximum determinant coefficient. However, for MARS, the best model was chosen based on the GCV, minimum MSE, maximum determinant coefficient. Abstract in Bahasa Indonesia : Analisa regresi digunakan untuk melihat pengaruh variabel independen terhadap variabel dependent dengan terlebih dulu melihat pola hubungan variabel tersebut. Hal ini dapat dilakukan dengan melalui dua pendekatan. Pendekatan yang paling umum dan seringkali digunakan adalah pendekatan parametrik. Pendekatan parametrik mengasumsikan bentuk model sudah ditentukan. Apabila tidak ada informasi apapun tentang bentuk dari fungsi regresi, maka pendekatan yang digunakan adalah pendekatan nonparametrik. (Haerdle, 1990. Karena pendekatan tidak tergantung pada asumsi bentuk kurva tertentu, sehingga memberikan fleksibelitas yang lebih besar. Tujuan penelitian ini
Pizarro Bondia, Carlos; Arasa Marti, Jose; de Lasarte, Marta; Pujol Ramo, Jaume; Arjona Carbonell, Mª Montserrat; Vilaseca Ricart, Meritxell
2008-01-01
La principal motivación de este trabajo es la búsqueda de una única expresión matemática que permita reproducir distribuciones espectrales de forma general. Para ello se consideran polinomios B-Spline rotacionales de segundo orden como expresión matemática base para dicha reproducción. El objetivo fundamental de este trabajo es, por tanto, la determinación de los coeficientes de los polinomios B-Spline que permitan reproducir distribuciones espectrales, así como la evaluación de la exactit...
Duddu, Ravindra
2011-10-05
We present a numerical formulation aimed at modeling the nonlinear response of elastic materials using large deformation continuum mechanics in three dimensions. This finite element formulation is based on the Eulerian description of motion and the transport of the deformation gradient. When modeling a nearly incompressible solid, the transport of the deformation gradient is decomposed into its isochoric part and the Jacobian determinant as independent fields. A homogeneous isotropic hyperelastic solid is assumed and B-splines-based finite elements are used for the spatial discretization. A variational multiscale residual-based approach is employed to stabilize the transport equations. The performance of the scheme is explored for both compressible and nearly incompressible applications. The numerical results are in good agreement with theory illustrating the viability of the computational scheme. © 2011 John Wiley & Sons, Ltd.
HUI Ping; SHI Ting-Yun; BAO Cheng-Guang
2003-01-01
Exciton energies as a function of radii of quantum dots in the range of 5-35 A are calculated based on effective mass approximation model with the B-spline technique and compared with experimental and other theoretical data for the CdS dots. This method leads to accurate and fast convergent exciton energy, which are in good agreement with experimental data in the whole confinement regime. The effect of penetration of wave function from the inside to the outside of the dots and the effect of dielectric constants are taken into account. The magnitudes of dynamical parameters are discussed. It is found that the different materials surrounding the CdS quantum dot affect not only the potential energy and Coulomb interaction energy of the system, but also the effective masses. The comparison shows that the effective mass approximation model can describe very well the quantum size effects observed experimentally on the exciton ground state energy.
Chevalier, Paul; Bouchon, Patrick; Pardo, Fabrice; Haïdar, Riad
2014-08-01
Focusing light onto nanostructures thanks to spherical lenses is a first step in enhancing the field and is widely used in applications. Nonetheless, the electromagnetic response of such nanostructures, which have subwavelength patterns, to a focused beam cannot be described by the simple ray tracing formalism. Here, we present a method for computing the response to a focused beam, based on the B-spline modal method adapted to nanostructures in conical mounting. The eigenmodes are computed in each layer for both polarizations and are then combined for the computation of scattering matrices. The simulation of a Gaussian focused beam is obtained thanks to a truncated decomposition into plane waves computed on a single period, which limits the computation burden. PMID:25121523
Correlation studies for B-spline modeled F2 Chapman parameters obtained from FORMOSAT-3/COSMIC data
Limberger, M.; Liang, W.; Schmidt, M.; Dettmering, D.; Hernández-Pajares, M.; Hugentobler, U.
2014-12-01
The determination of ionospheric key quantities such as the maximum electron density of the F2 layer NmF2, the corresponding F2 peak height hmF2 and the F2 scale height HF2 are of high relevance in 4-D ionosphere modeling to provide information on the vertical structure of the electron density (Ne). The Ne distribution with respect to height can, for instance, be modeled by the commonly accepted F2 Chapman layer. An adequate and observation driven description of the vertical Ne variation can be obtained from electron density profiles (EDPs) derived by ionospheric radio occultation measurements between GPS and low Earth orbiter (LEO) satellites. For these purposes, the six FORMOSAT-3/COSMIC (F3/C) satellites provide an excellent opportunity to collect EDPs that cover most of the ionospheric region, in particular the F2 layer. For the contents of this paper, F3/C EDPs have been exploited to determine NmF2, hmF2 and HF2 within a regional modeling approach. As mathematical base functions, endpoint-interpolating polynomial B-splines are considered to model the key parameters with respect to longitude, latitude and time. The description of deterministic processes and the verification of this modeling approach have been published previously in Limberger et al. (2013), whereas this paper should be considered as an extension dealing with related correlation studies, a topic to which less attention has been paid in the literature. Relations between the B-spline series coefficients regarding specific key parameters as well as dependencies between the three F2 Chapman key parameters are in the main focus. Dependencies are interpreted from the post-derived correlation matrices as a result of (1) a simulated scenario without data gaps by taking dense, homogenously distributed profiles into account and (2) two real data scenarios on 1 July 2008 and 1 July 2012 including sparsely, inhomogeneously distributed F3/C EDPs. Moderate correlations between hmF2 and HF2 as well as inverse
The self-controlling feedback control method requires an external periodic oscillator with special design, which is technically challenging. This paper proposes a chaos control method based on time series non-uniform rational B-splines (SNURBS for short) signal feedback. It first builds the chaos phase diagram or chaotic attractor with the sampled chaotic time series and any target orbit can then be explicitly chosen according to the actual demand. Second, we use the discrete timing sequence selected from the specific target orbit to build the corresponding external SNURBS chaos periodic signal, whose difference from the system current output is used as the feedback control signal. Finally, by properly adjusting the feedback weight, we can quickly lead the system to an expected status. We demonstrate both the effectiveness and efficiency of our method by applying it to two classic chaotic systems, i.e., the Van der Pol oscillator and the Lorenz chaotic system. Further, our experimental results show that compared with delayed feedback control, our method takes less time to obtain the target point or periodic orbit (from the starting point) and that its parameters can be fine-tuned more easily
Shao, Chenxi; Xue, Yong; Fang, Fang; Bai, Fangzhou; Yin, Peifeng; Wang, Binghong
2015-07-01
The self-controlling feedback control method requires an external periodic oscillator with special design, which is technically challenging. This paper proposes a chaos control method based on time series non-uniform rational B-splines (SNURBS for short) signal feedback. It first builds the chaos phase diagram or chaotic attractor with the sampled chaotic time series and any target orbit can then be explicitly chosen according to the actual demand. Second, we use the discrete timing sequence selected from the specific target orbit to build the corresponding external SNURBS chaos periodic signal, whose difference from the system current output is used as the feedback control signal. Finally, by properly adjusting the feedback weight, we can quickly lead the system to an expected status. We demonstrate both the effectiveness and efficiency of our method by applying it to two classic chaotic systems, i.e., the Van der Pol oscillator and the Lorenz chaotic system. Further, our experimental results show that compared with delayed feedback control, our method takes less time to obtain the target point or periodic orbit (from the starting point) and that its parameters can be fine-tuned more easily.
Xiao, Xun; Geyer, Veikko F; Bowne-Anderson, Hugo; Howard, Jonathon; Sbalzarini, Ivo F
2016-08-01
Biological filaments, such as actin filaments, microtubules, and cilia, are often imaged using different light-microscopy techniques. Reconstructing the filament curve from the acquired images constitutes the filament segmentation problem. Since filaments have lower dimensionality than the image itself, there is an inherent trade-off between tracing the filament with sub-pixel accuracy and avoiding noise artifacts. Here, we present a globally optimal filament segmentation method based on B-spline vector level-sets and a generalized linear model for the pixel intensity statistics. We show that the resulting optimization problem is convex and can hence be solved with global optimality. We introduce a simple and efficient algorithm to compute such optimal filament segmentations, and provide an open-source implementation as an ImageJ/Fiji plugin. We further derive an information-theoretic lower bound on the filament segmentation error, quantifying how well an algorithm could possibly do given the information in the image. We show that our algorithm asymptotically reaches this bound in the spline coefficients. We validate our method in comprehensive benchmarks, compare with other methods, and show applications from fluorescence, phase-contrast, and dark-field microscopy. PMID:27104582
Côrtes, A.M.A.
2015-02-20
The recently introduced divergence-conforming B-spline discretizations allow the construction of smooth discrete velocity–pressure pairs for viscous incompressible flows that are at the same time inf-sup stable and pointwise divergence-free. When applied to discretized Stokes equations, these spaces generate a symmetric and indefinite saddle-point linear system. Krylov subspace methods are usually the most efficient procedures to solve such systems. One of such methods, for symmetric systems, is the Minimum Residual Method (MINRES). However, the efficiency and robustness of Krylov subspace methods is closely tied to appropriate preconditioning strategies. For the discrete Stokes system, in particular, block-diagonal strategies provide efficient preconditioners. In this article, we compare the performance of block-diagonal preconditioners for several block choices. We verify how the eigenvalue clustering promoted by the preconditioning strategies affects MINRES convergence. We also compare the number of iterations and wall-clock timings. We conclude that among the building blocks we tested, the strategy with relaxed inner conjugate gradients preconditioned with incomplete Cholesky provided the best results.
Verónica S. Martínez
2015-12-01
Full Text Available Metabolic flux analysis (MFA is widely used to estimate intracellular fluxes. Conventional MFA, however, is limited to continuous cultures and the mid-exponential growth phase of batch cultures. Dynamic MFA (DMFA has emerged to characterize time-resolved metabolic fluxes for the entire culture period. Here, the linear DMFA approach was extended using B-spline fitting (B-DMFA to estimate mass balanced fluxes. Smoother fits were achieved using reduced number of knots and parameters. Additionally, computation time was greatly reduced using a new heuristic algorithm for knot placement. B-DMFA revealed that Chinese hamster ovary cells shifted from 37 °C to 32 °C maintained a constant IgG volume-specific productivity, whereas the productivity for the controls peaked during mid-exponential growth phase and declined afterward. The observed 42% increase in product titer at 32 °C was explained by a prolonged cell growth with high cell viability, a larger cell volume and a more stable volume-specific productivity.
Filtering of Airborne LiDAR Point Cloud Based on Variable Radius Circle and B-spline Fitting
ZHENG Jitao
2015-12-01
Full Text Available Filtering is the key to acquire digital terrain model from airborne LiDAR point cloud. In this paper, a new LiDAR point cloud filtering method is proposed. First, the scanning lines formed with point sequence are obtained throuth scanning the point cloud along the same direction. Then a circle with the variable radius rolling over the bottom of these scanning lines, the purpose is to acquire the points on the ground surface and delete the points on the objects at the same time. The next step is interval sampling from the scanning lines. On this basis, after fitting terrain surface with uniform B-spline surface, every point is projected to the fitting surface and calculate its height. According to compare the real height and its projection height to judge every point is on the terrain surface or not. The experiments show that filtering precision of the algorithm proposed in this paper is improved 1 to 5 times of the traditional methods, it can be used for the city, mountains and forest, and the time complexity of the algorithm is O(n.
Zhiwei Pan
2016-05-01
Full Text Available Global look-up table strategy proposed recently has been proven to be an efficient method to accelerate the interpolation, which is the most time-consuming part in the iterative sub-pixel digital image correlation (DIC algorithms. In this paper, a global look-up table strategy with cubic B-spline interpolation is developed for the DIC method based on the inverse compositional Gauss–Newton (IC-GN algorithm. The performance of this strategy, including accuracy, precision, and computation efficiency, is evaluated through a theoretical and experimental study, using the one with widely employed bicubic interpolation as a benchmark. The global look-up table strategy with cubic B-spline interpolation improves significantly the accuracy of the IC-GN algorithm-based DIC method compared with the one using the bicubic interpolation, at a trivial price of computation efficiency.
Image registration in intra-oral radiography
Leung, CC; Yiu, KL; Zee, KY; Tsui, WK
2005-01-01
Image registration is one of the image processing methods which is widely used in computer vision, pattern recognition, and medical imaging. In digital subtraction radiography, image registration is one of the important prerequisites to match the reference and subsequent images. In this paper, we propose an automatic non-rigid registration method namely curvature-based registration that relies on a curvature based penalizing term and its application on dental radiography. The regularizing ter...
In a joint experimental and theoretical effort, we carried out a detailed study of electron impact excitation of the 4p5 5s states of Kr. We present angle-differential cross sections over the entire angular range (00-1800) for a number of energies in the near-threshold region, as well as energy scans for selected angles. The experimental results are in very satisfactory agreement with predictions from a fully relativistic Dirac B-spline R-matrix model.
This research concerns with the development of a linear three-dimensional numerical model in a quantum environment. We use the semi inverse variational method together with B-spline bases to extract the structures of bound states of the Schrödinger equation. The model performances are demonstrated for the Coulomb type problem. From realistic examples, some state configurations are presented to illustrate the effectiveness and the exactitude of the proposed method. (general)
Soares Filho, Djalma Manoel; Eckhardt, Wietze; Braganca, Ricardo Silva Nunes de [PETROBRAS, Rio de Janeiro, RJ (Brazil). Centro de Pesquisas. Setor de Geofisica e Estratigrafia
1997-07-01
A gauss-Newton type method for seismic crosswell traveltime tomography based on B-splines products linear combinations for velocity field representation is presented. Knot points and order are established through Dierckx approximation technique according to the initial guess model parameters. The Frechet derivatives are contour integrations along the rays and the regularization factor is estimated using the concepts of condition number and Backus-Gilbert spread factor for normalized resolution matrices. A numerical example demonstrates the robustness of the method. (author)
DI TOMMASO, DEVIS
2006-01-01
Lo scopo della presente tesi è l'applicazione e lo sviluppo del metodo LCAO B-spline DFT, basato sulla Combinazione Lineare degli Orbitali Atomici (LCAO), sulla Teoria del Funzionale Densità (DFT) e sull'impiego delle funzioni di base Bspline, per lo studio della fotoionizzazione molecolare. Nel Capitolo 1 viene considerata una breve introduzione al processo di fotoionizzazione e alla teoria della fotoionizzazione molecolare. Nel Capitolo 2 vengono presentati i metodi computazionali utilizzat...
Dong Chun-wang; Ye Yang; Zhang Jian-qiang; Zhu hong-kai; Liu Fei
2014-01-01
In order to ifnd an effective method of detecting thrips defect on green-peel citrus, a defect segmentation method was developed using a single threshold value based on combination of characteristic wavelengths principal component analysis (PCA) and B-spline lighting correction method in this study. At ifrst, four characteristic wavelengths (523, 587, 700 and 768 nm) were obtained using PCA of Vis-NIR (visible and near-infrared) bands and analysis of weighting coefifcients;secondarily, PCA was performed using characteristic wavelengths and the second principal component (PC2) was selected to classify images;then, B-spline lighting correction method was proposed to overcome the inlfuence of lighting non-uniform on citrus when thrips defect was segmented;ifnally, thrips defect on citrus was extracted by global threshold segmentation and morphological image processing. The experimental results show that thrips defect in citrus can be detected with an accuracy of 96.5%by characteristic wavelengths PCA and B-spline lighting correction method. This study shows that thrips defect on green-peel citrus can be effectively identiifed using hyperspectral imaging technology.
Stochastic Image Registration with User Constraints
Kolesov, Ivan; Lee, JeHoon; Vela, Patricio; Tannenbaum, Allen
2013-01-01
Constrained registration is an active area of research and is the focus of this work. This note describes a non-rigid image registration framework for incorporating landmark constraints. Points that must remain stationary are selected, the user chooses the spatial extent of the inputs, and an automatic step computes the deformable registration, respecting the constraints. Parametrization of the deformation field is by an additive composition of a similarity transformation and a set of Gaussia...
Mass preserving image registration for lung CT
Gorbunova, Vladlena; Sporring, Jon; Lo, Pechin;
2012-01-01
This paper presents a mass preserving image registration algorithm for lung CT images. To account for the local change in lung tissue intensity during the breathing cycle, a tissue appearance model based on the principle of preservation of total lung mass is proposed. This model is incorporated...... into a standard image registration framework with a composition of a global affine and several free-form B-Spline transformations with increasing grid resolution. The proposed mass preserving registration method is compared to registration using the sum of squared intensity differences as a similarity...... function on four groups of data: 44 pairs of longitudinal inspiratory chest CT scans with small difference in lung volume; 44 pairs of longitudinal inspiratory chest CT scans with large difference in lung volume; 16 pairs of expiratory and inspiratory CT scans; and 5 pairs of images extracted at end exhale...
Analysis of deformable image registration accuracy using computational modeling
Computer aided modeling of anatomic deformation, allowing various techniques and protocols in radiation therapy to be systematically verified and studied, has become increasingly attractive. In this study the potential issues in deformable image registration (DIR) were analyzed based on two numerical phantoms: One, a synthesized, low intensity gradient prostate image, and the other a lung patient's CT image data set. Each phantom was modeled with region-specific material parameters with its deformation solved using a finite element method. The resultant displacements were used to construct a benchmark to quantify the displacement errors of the Demons and B-Spline-based registrations. The results show that the accuracy of these registration algorithms depends on the chosen parameters, the selection of which is closely associated with the intensity gradients of the underlying images. For the Demons algorithm, both single resolution (SR) and multiresolution (MR) registrations required approximately 300 iterations to reach an accuracy of 1.4 mm mean error in the lung patient's CT image (and 0.7 mm mean error averaged in the lung only). For the low gradient prostate phantom, these algorithms (both SR and MR) required at least 1600 iterations to reduce their mean errors to 2 mm. For the B-Spline algorithms, best performance (mean errors of 1.9 mm for SR and 1.6 mm for MR, respectively) on the low gradient prostate was achieved using five grid nodes in each direction. Adding more grid nodes resulted in larger errors. For the lung patient's CT data set, the B-Spline registrations required ten grid nodes in each direction for highest accuracy (1.4 mm for SR and 1.5 mm for MR). The numbers of iterations or grid nodes required for optimal registrations depended on the intensity gradients of the underlying images. In summary, the performance of the Demons and B-Spline registrations have been quantitatively evaluated using numerical phantoms. The results show that parameter
and body motion phantoms with a varying extent and character for each individual. In addition to these phantoms, we recorded the position of markers placed on the chest of the volunteers for the body motion studies, which could be used as external motion measurement. Using these phantoms and external motion data, investigators will be able to test their motion correction approaches for realistic motion obtained from different individuals. The non-uniform rational B-spline data and the parameter files for these phantoms are freely available for downloading and can be used with the XCAT license. (paper)
Jednorozměrná disperzní analýza B-Spline varianty metody konečných prvků
Kolman, Radek; Plešek, Jiří; Okrouhlík, Miloslav; Gabriel, Dušan
Praha : ČVUT, Fakulta strojní,Ústav mechaniky,biomechaniky a mechatroniky, 2010 - (Španiel, M.; Horák, Z.), s. 28-35 ISBN 978-80-01-04670-8. [Výpočty konstrukcí metodou konečných prvků. Praha (CZ), 25.11.2010] R&D Projects: GA ČR(CZ) GPP101/10/P376; GA ČR(CZ) GA101/07/1471; GA ČR(CZ) GA101/09/1630 Institutional research plan: CEZ:AV0Z20760514 Keywords : B-spline finite element method * dispersion analysis * wave propagation Subject RIV: BI - Acoustics
Modeling and animation of 3D characters with ball B-splines%基于球B样条的3D人物角色建模与动画
徐欣; 武仲科; 周明全; 骆岩林
2009-01-01
提出一种新的人物角色的快速建模和动画的方法.通过使用球B样条实现以骨架为基础自动构建角色模型以及数据驱动的角色动画生成方法.给定人物角色骨架后,依据人体测量学根据骨骼长度计算出关节点处的半径.通过插值关节点以及半径生成球B样条曲线曲面表示的角色模型.根据所获取的运动数据(来自运动捕获),通过更新控制顶点来更新角色模型,从而实现角色的实时动画.%This paper proposed a new approach for modeling and animation of characters. Through using ball B-splines, an automatic character building model based on skeleton was realized and a character animation generation method driven by data was obtained. When the skeleton of a character was given, these radii on joints were obtained according to the length of bones in anthropometry. Then ball B-spline curves and surfaces were generated to represent characters through interpolating joint points and radii on them. The modeling method is rapid. Furthermore, by using motion data to update control points, real-time animation is achieved.
基于B样条梯度的融合图像算法%Gradient Image Fusion Algorithm Based on B-Spline
刘军华; 雷超阳
2012-01-01
图像融合在目标识别、计算机视觉、医学、遥感图像处理等领域应用广泛．利用B样条梯度融合图像算法来解决图像融合，将融合的图像分成3×3予块，计算每个子块的B样务梯度．把待融合图像对应的子块梯度模值累加求和，计算每个子块梯度模值的权值并与每个子块相乘并求和，得到融合后的图像子块．模拟实验结果表明该方法是有效的，优于传统的主分量分析、小波变换和Laplacian金字塔等方法．%Image fusion is widely used in target recognition, computer vision, medicine, remote sensing, im- age processing and many other fields. This paper proposes an image fusion method based on B-spline gra- dient. It first divides the source image into small blocks of 3 × 3 size, calculates separately the B-spline gradient for each sub-block,adds up the gradient modulus values,calculates the weight value of each sub- block gradient modulus value,multiplies each sub-block and respective weight value and obtains the fu- sion image sub-block. Simulation experiments and analysis show that this method is effective and has a better performance than traditional methods such as principal component analysis, wavelet transforma- tion and the Laplacian pyramid.
Fast fluid registration of medical images
Bro-Nielsen, Morten; Gramkow, Claus
This paper offers a new fast algorithm for non-rigid viscous fluid registration of medical images that is at least an order of magnitude faster than the previous method by (Christensen et al., 1994). The core algorithm in the fluid registration method is based on a linear elastic deformation of the...... velocity field of the fluid. Using the linearity of this deformation we derive a convolution filter which we use in a scale-space framework. We also demonstrate that the `demon'-based registration method of (Thirion, 1996) can be seen as an approximation to the fluid registration method and point to...
Evaluation of various deformable image registration algorithms for thoracic images
We evaluated the accuracy of one commercially available and three publicly available deformable image registration (DIR) algorithms for thoracic four-dimensional (4D) computed tomography (CT) images. Five patients with esophagus cancer were studied. Datasets of the five patients were provided by DIR-lab (dir-lab.com) and consisted of thoracic 4D CT images and a coordinate list of anatomical landmarks that had been manually identified. Expert landmark correspondence was used for evaluating DIR spatial accuracy. First, the manually measured displacement vector field (mDVF) was obtained from the coordinate list of anatomical landmarks. Then the automatically calculated displacement vector field (aDVF) was calculated by using the following four DIR algorithms: B-spine implemented in Velocity AI (Velocity Medical, Atlanta, GA, USA), free-form deformation (FFD), Horn–Schunk optical flow (OF) and Demons in DIRART of MATLAB software. Registration error is defined as the difference between mDVF and aDVF. The mean 3D registration errors were 2.7 ± 0.8 mm for B-spline, 3.6 ± 1.0 mm for FFD, 2.4 ± 0.9 mm for OF and 2.4 ± 1.2 mm for Demons. The results showed that reasonable accuracy was achieved in B-spline, OF and Demons, and that these algorithms have the potential to be used for 4D dose calculation, automatic image segmentation and 4D CT ventilation imaging in patients with thoracic cancer. However, for all algorithms, the accuracy might be improved by using the optimized parameter setting. Furthermore, for B-spline in Velocity AI, the 3D registration error was small with displacements of less than ∼10 mm, indicating that this software may be useful in this range of displacements
Purpose: Three deformable image registration (DIR) algorithms are utilized to perform deformable dose accumulation for head and neck tomotherapy treatment, and the differences of the accumulated doses are evaluated. Methods: Daily MVCT data for 10 patients with pathologically proven nasopharyngeal cancers were analyzed. The data were acquired using tomotherapy (TomoTherapy, Accuray) at the PLA General Hospital. The prescription dose to the primary target was 70Gy in 33 fractions.Three DIR methods (B-spline, Diffeomorphic Demons and MIMvista) were used to propagate parotid structures from planning CTs to the daily CTs and accumulate fractionated dose on the planning CTs. The mean accumulated doses of parotids were quantitatively compared and the uncertainties of the propagated parotid contours were evaluated using Dice similarity index (DSI). Results: The planned mean dose of the ipsilateral parotids (32.42±3.13Gy) was slightly higher than those of the contralateral parotids (31.38±3.19Gy)in 10 patients. The difference between the accumulated mean doses of the ipsilateral parotids in the B-spline, Demons and MIMvista deformation algorithms (36.40±5.78Gy, 34.08±6.72Gy and 33.72±2.63Gy ) were statistically significant (B-spline vs Demons, P<0.0001, B-spline vs MIMvista, p =0.002). And The difference between those of the contralateral parotids in the B-spline, Demons and MIMvista deformation algorithms (34.08±4.82Gy, 32.42±4.80Gy and 33.92±4.65Gy ) were also significant (B-spline vs Demons, p =0.009, B-spline vs MIMvista, p =0.074). For the DSI analysis, the scores of B-spline, Demons and MIMvista DIRs were 0.90, 0.89 and 0.76. Conclusion: Shrinkage of parotid volumes results in the dose increase to the parotid glands in adaptive head and neck radiotherapy. The accumulated doses of parotids show significant difference using the different DIR algorithms between kVCT and MVCT. Therefore, the volume-based criterion (i.e. DSI) as a quantitative evaluation of
Xu, S [Key Laboratory of Particle & Radiation Imaging (Tsinghua University), Ministry of Education, Beijing, 100084 China (China); Chinese PLA General Hospital, Beijing, 100853 China (China); Liu, B [Image processing center, Beihang University, Beijing, 100191 China (China)
2015-06-15
Purpose: Three deformable image registration (DIR) algorithms are utilized to perform deformable dose accumulation for head and neck tomotherapy treatment, and the differences of the accumulated doses are evaluated. Methods: Daily MVCT data for 10 patients with pathologically proven nasopharyngeal cancers were analyzed. The data were acquired using tomotherapy (TomoTherapy, Accuray) at the PLA General Hospital. The prescription dose to the primary target was 70Gy in 33 fractions.Three DIR methods (B-spline, Diffeomorphic Demons and MIMvista) were used to propagate parotid structures from planning CTs to the daily CTs and accumulate fractionated dose on the planning CTs. The mean accumulated doses of parotids were quantitatively compared and the uncertainties of the propagated parotid contours were evaluated using Dice similarity index (DSI). Results: The planned mean dose of the ipsilateral parotids (32.42±3.13Gy) was slightly higher than those of the contralateral parotids (31.38±3.19Gy)in 10 patients. The difference between the accumulated mean doses of the ipsilateral parotids in the B-spline, Demons and MIMvista deformation algorithms (36.40±5.78Gy, 34.08±6.72Gy and 33.72±2.63Gy ) were statistically significant (B-spline vs Demons, P<0.0001, B-spline vs MIMvista, p =0.002). And The difference between those of the contralateral parotids in the B-spline, Demons and MIMvista deformation algorithms (34.08±4.82Gy, 32.42±4.80Gy and 33.92±4.65Gy ) were also significant (B-spline vs Demons, p =0.009, B-spline vs MIMvista, p =0.074). For the DSI analysis, the scores of B-spline, Demons and MIMvista DIRs were 0.90, 0.89 and 0.76. Conclusion: Shrinkage of parotid volumes results in the dose increase to the parotid glands in adaptive head and neck radiotherapy. The accumulated doses of parotids show significant difference using the different DIR algorithms between kVCT and MVCT. Therefore, the volume-based criterion (i.e. DSI) as a quantitative evaluation of
Efficient Constrained Local Model Fitting for Non-Rigid Face Alignment
Wang, Yang; Cox, Mark; Sridharan, Sridha; Cohn, Jeffery F.
2009-01-01
Active appearance models (AAMs) have demonstrated great utility when being employed for non-rigid face alignment/tracking. The “simultaneous” algorithm for fitting an AAM achieves good non-rigid face registration performance, but has poor real time performance (2-3 fps). The “project-out” algorithm for fitting an AAM achieves faster than real time performance (> 200 fps) but suffers from poor generic alignment performance. In this paper we introduce an extension to a discriminative method for non-rigid face registration/tracking referred to as a constrained local model (CLM). Our proposed method is able to achieve superior performance to the “simultaneous” AAM algorithm along with real time fitting speeds (35 fps). We improve upon the canonical CLM formulation, to gain this performance, in a number of ways by employing: (i) linear SVMs as patch-experts, (ii) a simplified optimization criteria, and (iii) a composite rather than additive warp update step. Most notably, our simplified optimization criteria for fitting the CLM divides the problem of finding a single complex registration/warp displacement into that of finding N simple warp displacements. From these N simple warp displacements, a single complex warp displacement is estimated using a weighted least-squares constraint. Another major advantage of this simplified optimization lends from its ability to be parallelized, a step which we also theoretically explore in this paper. We refer to our approach for fitting the CLM as the “exhaustive local search” (ELS) algorithm. Experiments were conducted on the CMU Multi-PIE database. PMID:20046797
Spectral Distance Distributions for Non-rigid Objects
CAO Wei-guo; LI Hai-yang; LI Shi-rui; LIU Yu-jie; LI Hua
2013-01-01
Non-rigid shape deformation without tearing or stretching is called isometry. There are many difficulties to research non-rigid shape in Euclidean space. Therefore, non-rigid shapes are firstly embedded into a none-Euclidean space. Spectral space is chosen in this paper. Then three descriptors are proposed based on three spectral distances. The existence of zero-eigenvalue has negative effects on computation of spectral distance. Therefore the spectral distance should be computed from the first non-zero-eigenvalue. Experiments show that spectral distance distributions are very effective to describe the non-rigid shapes.
Purpose: The actual dose delivered to critical organs will differ from the simulated dose because of interfractional organ motion and deformation. Here, we developed a method to estimate the rectal dose in prostate intensity modulated radiation therapy with consideration to interfractional organ motion using daily megavoltage cone-beam computed tomography (MVCBCT). Methods and Materials: Under exemption status from our institutional review board, we retrospectively reviewed 231 series of MVCBCT of 8 patients with prostate cancer. On both planning CT (pCT) and MVCBCT images, the rectal contours were delineated and the CT value within the contours was replaced by the mean CT value within the pelvis, with the addition of 100 Hounsfield units. MVCBCT images were rigidly registered to pCT and then nonrigidly registered using B-Spline deformable image registration (DIR) with Velocity AI software. The concordance between the rectal contours on MVCBCT and pCT was evaluated using the Dice similarity coefficient (DSC). The dose distributions normalized for 1 fraction were also deformed and summed to estimate the actual total dose. Results: The DSC of all treatment fractions of 8 patients was improved from 0.75±0.04 (mean ±SD) to 0.90 ±0.02 by DIR. Six patients showed a decrease of the generalized equivalent uniform dose (gEUD) from total dose compared with treatment plans. Although the rectal volume of each treatment fraction did not show any correlation with the change in gEUD (R2=0.18±0.13), the displacement of the center of gravity of rectal contours in the anterior-posterior (AP) direction showed an intermediate relationship (R2=0.61±0.16). Conclusion: We developed a method for evaluation of rectal dose using DIR and MVCBCT images and showed the necessity of DIR for the evaluation of total dose. Displacement of the rectum in the AP direction showed a greater effect on the change in rectal dose compared with the rectal volume
Sandra P Mateus
2010-01-01
Full Text Available Dentro de las técnicas existentes de Inteligencia Artificial, se escogieron y adaptaron dos Redes Neuronales Artificiales (RNA para realizar el ajuste de uno de los elementos que definen una B-Spline Racional No Uniforme (NURBS y con ello obtener un modelado adecuado de la NURBS. Los elementos escogidos fueron los puntos de control. Las RNA utilizadas son las de Función de Base Radial y las de Kohonen o Mapas Auto-organizativos. Con base en el análisis de resultados y la caracterización de las RNA, la Función de Base Radial tuvo un desempeño más adecuado y óptimo para un número elevado de datos, lo cual es una desventaja de los Mapas Auto-organizativos. En este modelo se tiene que realizar procesos extras para determinar la neurona ganadora y realizar el reajuste de los pesos.In the existing techniques of Artificial Intelligence, two Artificial Neural Networks (ANN were selected and adapted to fit one of the elements that define a Non-Uniform Rational B-Spline (NURBS and thus obtaining an appropriate modeling of the NURBS. The selected elements were the checkpoints. The ANN used were the Radial Basis Function and the Kohonen model or Self-Organizing Maps. Based on the analysis of the results and characterization of the ANN the Radial Basis Function had a more appropriate and optimum performance for a large number of data, which is a disadvantage of the Self-Organizing Maps. In this model, additional processes must be done to determine the winning neuron and the weights must be refitted.
Liu, Hong; Yan, Meng; Song, Enmin; Wang, Jie; Wang, Qian; Jin, Renchao; Jin, Lianghai; Hung, Chih-Cheng
2016-05-01
Myocardial motion estimation of tagged cardiac magnetic resonance (TCMR) images is of great significance in clinical diagnosis and the treatment of heart disease. Currently, the harmonic phase analysis method (HARP) and the local sine-wave modeling method (SinMod) have been proven as two state-of-the-art motion estimation methods for TCMR images, since they can directly obtain the inter-frame motion displacement vector field (MDVF) with high accuracy and fast speed. By comparison, SinMod has better performance over HARP in terms of displacement detection, noise and artifacts reduction. However, the SinMod method has some drawbacks: 1) it is unable to estimate local displacements larger than half of the tag spacing; 2) it has observable errors in tracking of tag motion; and 3) the estimated MDVF usually has large local errors. To overcome these problems, we present a novel motion estimation method in this study. The proposed method tracks the motion of tags and then estimates the dense MDVF by using the interpolation. In this new method, a parameter estimation procedure for global motion is applied to match tag intersections between different frames, ensuring specific kinds of large displacements being correctly estimated. In addition, a strategy of tag motion constraints is applied to eliminate most of errors produced by inter-frame tracking of tags and the multi-level b-splines approximation algorithm is utilized, so as to enhance the local continuity and accuracy of the final MDVF. In the estimation of the motion displacement, our proposed method can obtain a more accurate MDVF compared with the SinMod method and our method can overcome the drawbacks of the SinMod method. However, the motion estimation accuracy of our method depends on the accuracy of tag lines detection and our method has a higher time complexity. PMID:26712656
Symmetric and Transitive Registration of Image Sequences
Oskar Škrinjar
2008-01-01
Full Text Available This paper presents a method for constructing symmetric and transitive algorithms for registration of image sequences from image registration algorithms that do not have these two properties. The method is applicable to both rigid and nonrigid registration and it can be used with linear or periodic image sequences. The symmetry and transitivity properties are satisfied exactly (up to the machine precision, that is, they always hold regardless of the image type, quality, and the registration algorithm as long as the computed transformations are invertable. These two properties are especially important in motion tracking applications since physically incorrect deformations might be obtained if the registration algorithm is not symmetric and transitive. The method was tested on two sequences of cardiac magnetic resonance images using two different nonrigid image registration algorithms. It was demonstrated that the transitivity and symmetry errors of the symmetric and transitive modification of the algorithms could be made arbitrary small when the computed transformations are invertable, whereas the corresponding errors for the nonmodified algorithms were on the order of the pixel size. Furthermore, the symmetric and transitive modification of the algorithms had higher registration accuracy than the nonmodified algorithms for both image sequences.
Objective: To study the morphology of normal liver and tumors by breathing motion of hepatocellular carcinoma patients, through comparing the modified demons algorithm and FFD algorithm based on B-spline, and combing four-dimensional computed tomography (4DCT). Methods: The 4DCT images of 8 HCC patients were segmented into 10-series which were named CT0, CT10…CT80, CT90 according to the respiratory phases, CT0 and CT50 are defined to be end-inhale and end-exhale respectively. CT50 was chosen as the reference image. We used the modified demons algorithm and FFD algorithm based on B-spline to deform the images. Linear interpolation was used in both mode 1 and mode 2. The normalized mutual information (NMI), Hausdorff distance (dH) and registration speed were used to verify the registration performance. Results: The average NMI for the end-inhale and end-exhale images of 8 HCC patients after demons registration in mode 1 improved 4.75% with FFD algorithm based on B-spline (P = 0.002). And the difference of dH after demons reduced 15.2% comparing with FFD model algorithm (P = 0.02). In addition, demons algorithm has the absolute advantage in registration speed (P = 0.036). Conclusions: The breathing movement for deformation of normal liver and tumor targets is significant. These two algorithms can achieve the registration of 4DCT images and the modified demons registration can deform 4DCT images effectively. (authors)
Park, Seyoun; Robinson, Adam; Quon, Harry; Kiess, Ana P.; Shen, Colette; Wong, John; Plishker, William; Shekhar, Raj; Lee, Junghoon
2016-03-01
In this paper, we propose a CT-CBCT registration method to accurately predict the tumor volume change based on daily cone-beam CTs (CBCTs) during radiotherapy. CBCT is commonly used to reduce patient setup error during radiotherapy, but its poor image quality impedes accurate monitoring of anatomical changes. Although physician's contours drawn on the planning CT can be automatically propagated to daily CBCTs by deformable image registration (DIR), artifacts in CBCT often cause undesirable errors. To improve the accuracy of the registration-based segmentation, we developed a DIR method that iteratively corrects CBCT intensities by local histogram matching. Three popular DIR algorithms (B-spline, demons, and optical flow) with the intensity correction were implemented on a graphics processing unit for efficient computation. We evaluated their performances on six head and neck (HN) cancer cases. For each case, four trained scientists manually contoured the nodal gross tumor volume (GTV) on the planning CT and every other fraction CBCTs to which the propagated GTV contours by DIR were compared. The performance was also compared with commercial image registration software based on conventional mutual information (MI), VelocityAI (Varian Medical Systems Inc.). The volume differences (mean±std in cc) between the average of the manual segmentations and automatic segmentations are 3.70+/-2.30 (B-spline), 1.25+/-1.78 (demons), 0.93+/-1.14 (optical flow), and 4.39+/-3.86 (VelocityAI). The proposed method significantly reduced the estimation error by 9% (B-spline), 38% (demons), and 51% (optical flow) over the results using VelocityAI. Although demonstrated only on HN nodal GTVs, the results imply that the proposed method can produce improved segmentation of other critical structures over conventional methods.
王玉珠; 杨丹; 张小洪
2007-01-01
This paper presents a multi-scale curvature product corner detection technique in the framework of B-spline curvature scale space. A scale product function is derived from the curvature product of the contour at different scales. Corners are constructed as the local maxima by thresholding the curvature product results across several scales. Through scale product, the localization accuracy and detection performance can be notably improved in terms of CNN criteria. Experiments also demonstrate that proposed method shows robustness to high frequency details and provides promising detection results.
Börger, Klaus; Schmidt, Michael; Dettmering, Denise; Limberger, Marco; Erdogan, Eren; Seitz, Florian; Brandert, Sylvia; Görres, Barbara; Kersten, Wilhelm; Bothmer, Volker; Hinrichs, Johannes; Venzmer, Malte; Mrotzek, Niclas
2016-04-01
Today, the observations of space geodetic techniques are usually available with a rather low latency which applies to space missions observing the solar terrestrial environment, too. Therefore, we can use all these measurements in near real-time to compute and to provide ionosphere information, e.g. the vertical total electron content (VTEC). GSSAC and BGIC support a project aiming at a service for providing ionosphere information. This project is called OPTIMAP, meaning "Operational Tool for Ionosphere Mapping and Prediction"; the scientific work is mainly done by the German Geodetic Research Institute of the Technical University Munich (DGFI-TUM) and the Institute for Astrophysics of the University of Goettingen (IAG). The OPTIMAP strategy for providing ionosphere target quantities of high quality, such as VTEC or the electron density, includes mathematical approaches and tools allowing for the model adaptation to the real observational scenario as a significant improvement w.r.t. the traditional well-established methods. For example, OPTIMAP combines different observation types such as GNSS (GPS, GLONASS), Satellite Altimetry (Jason-2), DORIS as well as radio-occultation measurements (FORMOSAT#3/COSMIC). All these observations run into a Kalman-filter to compute global ionosphere maps, i.e. VTEC, for the current instant of time and as a forecast for a couple of subsequent days. Mathematically, the global VTEC is set up as a series expansion in terms of two-dimensional basis functions defined as tensor products of trigonometric B-splines for longitude and polynomial B-splines for latitude. Compared to the classical spherical harmonics, B-splines have a localizing character and, therefore, can handle an inhomogeneous data distribution properly. Finally, B-splines enable a so-called multi-resolution-representation (MRR) enabling the combination of global and regional modelling approaches. In addition to the geodetic measurements, Sun observations are pre
A Variational Approach to Video Registration with Subspace Constraints
Garg, Ravi; Roussos, Anastasios; Agapito, Lourdes
2013-01-01
This paper addresses the problem of non-rigid video registration, or the computation of optical flow from a reference frame to each of the subsequent images in a sequence, when the camera views deformable objects. We exploit the high correlation between 2D trajectories of different points on the same non-rigid surface by assuming that the displacement of any point throughout the sequence can be expressed in a compact way as a linear combination of a low-rank motion basis. This subspace constr...
Medical Image Registration and Surgery Simulation
Bro-Nielsen, Morten
1996-01-01
surgery simulation Real-time deformable models, using finite element models of linear elasticity, have been developed for surgery simulation. The time consumption of the finite element method is reduced dramaticly, by the use of condensation techniques, explicit inversion of the stiffness matrix, and the...... implements the linear elasticity operator. Using the filter results in a speedup of at least an order of magnitude. Use of convolution hardware is expected to improve the performance even more. Mandibular growth for time registration of mandibles Non-rigid registration using a physically valid model of bone...
The objective of this work is to present the capabilities of the NUMERICS web platform for evaluation of the performance of image registration algorithms. The NUMERICS platform is a web accessible tool which provides access to dedicated numerical algorithms for registration and comparison of medical images (http://numerics.phys.uni-sofia.bg). The platform allows comparison of noisy medical images by means of different types of image comparison algorithms, which are based on statistical tests for outliers. The platform also allows 2D image registration with different techniques like Elastic Thin-Plate Spline registration, registration based on rigid transformations, affine transformations, as well as non-rigid image registration based on Mobius transformations. In this work we demonstrate how the platform can be used as a tool for evaluation of the quality of the image registration process. We demonstrate performance evaluation of a deformable image registration technique based on Mobius transformations. The transformations are applied with appropriate cost functions like: Mutual information, Correlation coefficient, Sum of Squared Differences. The accent is on the results provided by the platform to the user and their interpretation in the context of the performance evaluation of 2D image registration. The NUMERICS image registration and image comparison platform provides detailed statistical information about submitted image registration jobs and can be used to perform quantitative evaluation of the performance of different image registration techniques. (authors)
Registration of different imaging modalities such as CT, MRI, functional MRI (fMRI), positron (PET) and single photon (SPECT) emission tomography is used in many clinical applications. Determining the quality of any automatic registration procedure has been a challenging part because no gold standard is available to evaluate the registration. In this note we present a method, called the 'multiple sub-volume registration' (MSR) method, for assessing the consistency of a rigid registration. This is done by registering sub-images of one data set on the other data set, performing a crude non-rigid registration. By analysing the deviations (local deformations) of the sub-volume registrations from the full registration we get a measure of the consistency of the rigid registration. Registration of 15 data sets which include CT, MR and PET images for brain, head and neck, cervix, prostate and lung was performed utilizing a rigid body registration with normalized mutual information as the similarity measure. The resulting registrations were classified as good or bad by visual inspection. The resulting registrations were also classified using our MSR method. The results of our MSR method agree with the classification obtained from visual inspection for all cases (p < 0.02 based on ANOVA of the good and bad groups). The proposed method is independent of the registration algorithm and similarity measure. It can be used for multi-modality image data sets and different anatomic sites of the patient. (note)
Disperzní analýza B-spline metody konečných prvků v jednorozměrném šíření elastických vln
Kolman, Radek
Prague : Institute of Thermomechanics AS CR, v. v. i, 2011 - (Kolman, R.; Linkeová, I.; Okrouhlík, M.; Pařík, P.). s. 1-3 ISBN 978-80-87012-30-7. [Splines and IsoGeometric Analysis. 09.02.2011, Praha] R&D Projects: GA ČR(CZ) GPP101/10/P376; GA ČR(CZ) GA101/09/1630; GA ČR(CZ) GA101/07/1471 Institutional research plan: CEZ:AV0Z20760514 Keywords : finite element method * B-spline * dispersion errors Subject RIV: BI - Acoustics http://www.it.cas.cz/cs/siga2011
Evaluation of Registration Methods on Thoracic CT
Murphy, K.; van Ginneken, B.; Reinhardt, J.;
2011-01-01
method and the evaluation is independent, using the same criteria for all participants. All results are published on the EMPIRE10 website (http://empire10.isi.uu.nl). The challenge remains ongoing and open to new participants. Full results from 24 algorithms have been published at the time of writing......EMPIRE10 (Evaluation of Methods for Pulmonary Image REgistration 2010) is a public platform for fair and meaningful comparison of registration algorithms which are applied to a database of intra-patient thoracic CT image pairs. Evaluation of non-rigid registration techniques is a non trivial task....... This article details the organisation of the challenge, the data and evaluation methods and the outcome of the initial launch with 20 algorithms. The gain in knowledge and future work are discussed....
Registration of structurally dissimilar images in MRI-based brachytherapy
A serious challenge in image registration is the accurate alignment of two images in which a certain structure is present in only one of the two. Such topological changes are problematic for conventional non-rigid registration algorithms. We propose to incorporate in a conventional free-form registration framework a geometrical penalty term that minimizes the volume of the missing structure in one image. We demonstrate our method on cervical MR images for brachytherapy. The intrapatient registration problem involves one image in which a therapy applicator is present and one in which it is not. By including the penalty term, a substantial improvement in the surface distance to the gold standard anatomical position and the residual volume of the applicator void are obtained. Registration of neighboring structures, i.e. the rectum and the bladder is generally improved as well, albeit to a lesser degree. (paper)
Motion Geometric Active Contours: Tracking Nonrigid Objects in Clutter Background
Cen Feng (岑峰); Qi Feihu
2003-01-01
MGAC (Motion Geometric Active Contours), a new variational framework of geometric active contours to track multiple nonrigid moving objects in the clutter background in image sequences is presented. This framework, incorporating with the motion edge information, consists of motion detection and tracking stages. At the motion detection stage, the motion edge map provides an approximate edge map of the moving objects. Then, a tracking stage, merely using the static edge information, is considered to improve the motion detection result. Force field regularization method is used to extend the capture range of the edge attraction force field in both stages. Experiments demonstrate that the proposed framework is valid for tracking multiple nonrigid objects in the clutter background.
Nonrigid spherical real analytic hypersurfaces in C^2
Merker, Joel
2009-01-01
A Levi nondegenerate real analytic hypersurface M of C^2 represented in local coordinates (z, w) in C^2 by a complex defining equation of the form w = Theta (z, \\bar z, \\bar w) which satisfies an appropriate reality condition, is spherical if and only if its complex graphing function Theta satisfies an explicitly written sixth-order polynomial complex partial differential equation. In the rigid case (known before), this system simplifies considerably, but in the general nonrigid case, its com...
GS Department
2010-01-01
Following a reorganization in Building 55, please note that the Registration Service is now organised as follows : Ground floor: access cards (76903). 1st floor : registration of external firms’ personnel (76611 / 76622); car access stickers (76633); biometric registration (79710). Opening hours: 07-30 to 16-00 non-stop. GS-SEM Group General Infrastructure Services Department
Learning-Based Tracking of Complex Non-Rigid Motion
Qiang Wang; Hai-Zhou Ai; Guang-You Xu
2004-01-01
This paper describes a novel method for tracking complex non-rigid motions by learning the intrinsic object structure.The approach builds on and extends the studies on non-linear dimensionality reduction for object representation,object dynamics modeling and particle filter style tracking.First,the dimensionality reduction and density estimation algorithm is derived for unsupervised learning of object intrinsic representation,and the obtained non-rigid part of object state reduces even to 2-3 dimensions.Secondly the dynamical model is derived and trained based on this intrinsic representation.Thirdly the learned intrinsic object structure is integrated into a particle filter style tracker.It is shown that this intrinsic object representation has some interesting properties and based on which the newly derived dynamical model makes particle filter style tracker more robust and reliable.Extensive experiments are done on the tracking of challenging non-rigid motions such as fish twisting with selfocclusion,large inter-frame lip motion and facial expressions with global head rotation.Quantitative results are given to make comparisons between the newly proposed tracker and the existing tracker.The proposed method also has the potential to solve other type of tracking problems.
Montin, E; Cutri, E; Spadola, G; Testori, A; Pennati, G; Mainardi, L
2015-08-01
Several studies report the mechanical properties of skin tissues but their values largely depend on the measurement method. Therefore, we investigated the feasibility of recognizing the cellular constituents mechanical properties of pigmented skin by Confocal Laser Scanner Microscopy (CLSM). With this purpose, an healthy volunteer was examined in three areas nearby a pigmented skin lesion in two configurations: deforming and non deforming the nevus. The tissue displacement of the nevus was then assessed by means of deformable registration of the images in these two configurations. There are several registration strategy able to overcome this task, among them, we proposed two methods with different deformation models: a Free Form Deformation (FFD) model based on b-spline and a second one based on Demons Registration Algorithm (DRA). These two strategies need the definition of several parameters in order to obtain optimal registration performances. Thus, we tuned these parameters by means of simulated data and evaluated their registration abilities on the real in vivo CLSM acquisitions in the two configurations. The results showed that the registration using DRA had a better performance in comparison to the FFD one, in particular in two out of the three areas the DRA performance was significantly better than the FFD one. The registration procedure highlighted deformation differences among the chosen areas. PMID:26737734
Longitudinal brain MRI analysis with uncertain registration.
Simpson, Ivor J A; Woolrich, Mark W; Groves, Adrian R; Schnabel, Julia A
2011-01-01
In this paper we propose a novel approach for incorporating measures of spatial uncertainty, which are derived from non-rigid registration, into spatially normalised statistics. Current approaches to spatially normalised statistical analysis use point-estimates of the registration parameters. This is limiting as the registration will rarely be completely accurate, and therefore data smoothing is often used to compensate for the uncertainty of the mapping. We derive localised measurements of spatial uncertainty from a probabilistic registration framework, which provides a principled approach to image smoothing. We evaluate our method using longitudinal deformation features from a set of MR brain images acquired from the Alzheimer's Disease Neuroimaging Initiative. These images are spatially normalised using our probabilistic registration algorithm. The spatially normalised longitudinal features are adaptively smoothed according to the registration uncertainty. The proposed adaptive smoothing shows improved classification results, (84% correct Alzheimer's Disease vs. controls), over either not smoothing (79.6%), or using a Gaussian filter with sigma = 2mm (78.8%). PMID:21995084
Amir Pasha Mahmoudzadeh
2013-01-01
Full Text Available Interpolation has become a default operation in image processing and medical imaging and is one of the important factors in the success of an intensity-based registration method. Interpolation is needed if the fractional unit of motion is not matched and located on the high resolution (HR grid. The purpose of this work is to present a systematic evaluation of eight standard interpolation techniques (trilinear, nearest neighbor, cubic Lagrangian, quintic Lagrangian, hepatic Lagrangian, windowed Sinc, B-spline 3rd order, and B-spline 4th order and to compare the effect of cost functions (least squares (LS, normalized mutual information (NMI, normalized cross correlation (NCC, and correlation ratio (CR for optimized automatic image registration (OAIR on 3D spoiled gradient recalled (SPGR magnetic resonance images (MRI of the brain acquired using a 3T GE MR scanner. Subsampling was performed in the axial, sagittal, and coronal directions to emulate three low resolution datasets. Afterwards, the low resolution datasets were upsampled using different interpolation methods, and they were then compared to the high resolution data. The mean squared error, peak signal to noise, joint entropy, and cost functions were computed for quantitative assessment of the method. Magnetic resonance image scans and joint histogram were used for qualitative assessment of the method.
Han, Lianghao; Hipwell, John H; Eiben, Björn; Barratt, Dean; Modat, Marc; Ourselin, Sebastien; Hawkes, David J
2014-03-01
Preoperative diagnostic magnetic resonance (MR) breast images can provide good contrast between different tissues and 3-D information about suspicious tissues. Aligning preoperative diagnostic MR images with a patient in the theatre during breast conserving surgery could assist surgeons in achieving the complete excision of cancer with sufficient margins. Typically, preoperative diagnostic MR breast images of a patient are obtained in the prone position, while surgery is performed in the supine position. The significant shape change of breasts between these two positions due to gravity loading, external forces and related constraints makes the alignment task extremely difficult. Our previous studies have shown that either nonrigid intensity-based image registration or biomechanical modelling alone are limited in their ability to capture such a large deformation. To tackle this problem, we proposed in this paper a nonlinear biomechanical model-based image registration method with a simultaneous optimization procedure for both the material parameters of breast tissues and the direction of the gravitational force. First, finite element (FE) based biomechanical modelling is used to estimate a physically plausible deformation of the pectoral muscle and the major deformation of breast tissues due to gravity loading. Then, nonrigid intensity-based image registration is employed to recover the remaining deformation that FE analyses do not capture due to the simplifications and approximations of biomechanical models and the uncertainties of external forces and constraints. We assess the registration performance of the proposed method using the target registration error of skin fiducial markers and the Dice similarity coefficient (DSC) of fibroglandular tissues. The registration results on prone and supine MR image pairs are compared with those from two alternative nonrigid registration methods for five breasts. Overall, the proposed algorithm achieved the best registration
Purpose: Measuring changes in lung perfusion resulting from radiation therapy dose requires registration of the functional imaging to the radiation therapy treatment planning scan. This study investigates registration accuracy and utility for positron emission tomography (PET)/computed tomography (CT) perfusion imaging in radiation therapy for non–small cell lung cancer. Methods: 68Ga 4-dimensional PET/CT ventilation-perfusion imaging was performed before, during, and after radiation therapy for 5 patients. Rigid registration and deformable image registration (DIR) using B-splines and Demons algorithms was performed with the CT data to obtain a deformation map between the functional images and planning CT. Contour propagation accuracy and correspondence of anatomic features were used to assess registration accuracy. Wilcoxon signed-rank test was used to determine statistical significance. Changes in lung perfusion resulting from radiation therapy dose were calculated for each registration method for each patient and averaged over all patients. Results: With B-splines/Demons DIR, median distance to agreement between lung contours reduced modestly by 0.9/1.1 mm, 1.3/1.6 mm, and 1.3/1.6 mm for pretreatment, midtreatment, and posttreatment (P<.01 for all), and median Dice score between lung contours improved by 0.04/0.04, 0.05/0.05, and 0.05/0.05 for pretreatment, midtreatment, and posttreatment (P<.001 for all). Distance between anatomic features reduced with DIR by median 2.5 mm and 2.8 for pretreatment and midtreatment time points, respectively (P=.001) and 1.4 mm for posttreatment (P>.2). Poorer posttreatment results were likely caused by posttreatment pneumonitis and tumor regression. Up to 80% standardized uptake value loss in perfusion scans was observed. There was limited change in the loss in lung perfusion between registration methods; however, Demons resulted in larger interpatient variation compared with rigid and B-splines registration. Conclusions
Hardcastle, Nicholas, E-mail: nick.hardcastle@gmail.com [Department of Physical Sciences, Peter MacCallum Cancer Centre, East Melbourne (Australia); Centre for Medical Radiation Physics, University of Wollongong, Wollongong (Australia); Hofman, Michael S. [Molecular Imaging, Centre for Cancer Imaging, Peter MacCallum Cancer Centre, Melbourne (Australia); Hicks, Rodney J. [Molecular Imaging, Centre for Cancer Imaging, Peter MacCallum Cancer Centre, Melbourne (Australia); Department of Medicine, University of Melbourne, Melbourne (Australia); Callahan, Jason [Molecular Imaging, Centre for Cancer Imaging, Peter MacCallum Cancer Centre, Melbourne (Australia); Kron, Tomas [Department of Medical Imaging and Radiation Sciences, Monash University, Clayton (Australia); The Sir Peter MacCallum Department of Oncology, Melbourne University, Victoria (Australia); MacManus, Michael P.; Ball, David L. [Division of Radiation Oncology and Cancer Imaging, Peter MacCallum Cancer Centre, East Melbourne (Australia); The Sir Peter MacCallum Department of Oncology, University of Melbourne, Melbourne (Australia); Jackson, Price [Department of Physical Sciences, Peter MacCallum Cancer Centre, East Melbourne (Australia); Siva, Shankar [Division of Radiation Oncology and Cancer Imaging, Peter MacCallum Cancer Centre, East Melbourne (Australia)
2015-09-01
Purpose: Measuring changes in lung perfusion resulting from radiation therapy dose requires registration of the functional imaging to the radiation therapy treatment planning scan. This study investigates registration accuracy and utility for positron emission tomography (PET)/computed tomography (CT) perfusion imaging in radiation therapy for non–small cell lung cancer. Methods: {sup 68}Ga 4-dimensional PET/CT ventilation-perfusion imaging was performed before, during, and after radiation therapy for 5 patients. Rigid registration and deformable image registration (DIR) using B-splines and Demons algorithms was performed with the CT data to obtain a deformation map between the functional images and planning CT. Contour propagation accuracy and correspondence of anatomic features were used to assess registration accuracy. Wilcoxon signed-rank test was used to determine statistical significance. Changes in lung perfusion resulting from radiation therapy dose were calculated for each registration method for each patient and averaged over all patients. Results: With B-splines/Demons DIR, median distance to agreement between lung contours reduced modestly by 0.9/1.1 mm, 1.3/1.6 mm, and 1.3/1.6 mm for pretreatment, midtreatment, and posttreatment (P<.01 for all), and median Dice score between lung contours improved by 0.04/0.04, 0.05/0.05, and 0.05/0.05 for pretreatment, midtreatment, and posttreatment (P<.001 for all). Distance between anatomic features reduced with DIR by median 2.5 mm and 2.8 for pretreatment and midtreatment time points, respectively (P=.001) and 1.4 mm for posttreatment (P>.2). Poorer posttreatment results were likely caused by posttreatment pneumonitis and tumor regression. Up to 80% standardized uptake value loss in perfusion scans was observed. There was limited change in the loss in lung perfusion between registration methods; however, Demons resulted in larger interpatient variation compared with rigid and B-splines registration
Non-Rigid Object Tracking by Anisotropic Kernel Mean Shift
无
2007-01-01
Mean shift, an iterative procedure that shifts each data point to the average of data points in its neighborhood, has been applied to object tracker. However, the traditional mean shift tracker by isotropic kernel often loses the object with the changing object structure in video sequences, especially when the object structure varies fast. This paper proposes a non-rigid object tracker by anisotropic kernel mean shift in which the shape, scale, and orientation of the kernels adapt to the changing object structure. The experimental results show that the new tracker is self-adaptive and approximately twice faster than the traditional tracker, which ensures the robustness and real time of tracking.
Full Non-Rigid Group and Symmetry of Dimethyltrichlorophosphorus
ASHRAFI; AliReza
2005-01-01
In this work, a simple method is described, by means of which it is possible to calculate character tables for the symmetry group of molecules consisting of a number of NH3 groups attached to a rigid framework. The full non-rigid group (f-NRG) of dimethyltrichlorophosphorus with the symmetry group D3h was studied. It has been proven that it is a group of order 216 with 27 conjugacy classes and its character table computed. Finally, the Permutation-lnversion group of this molecule was calculated.
Purpose: To generate virtual phantoms with clinically relevant deformation and use them to objectively evaluate geometric and dosimetric uncertainties of deformable image registration (DIR) algorithms. Methods: Ten lung cancer patients undergoing adaptive 3DCRT planning were selected. For each patient, a pair of planning CT (pCT) and replanning CT (rCT) were used as the basis for virtual phantom generation. Manually adjusted meshes were created for selected ROIs (e.g. PTV, lungs, spinal cord, esophagus, and heart) on pCT and rCT. The mesh vertices were input into a thin-plate spline algorithm to generate a reference displacement vector field (DVF). The reference DVF was used to deform pCT to generate a simulated replanning CT (srCT) that was closely matched to rCT. Three DIR algorithms (Demons, B-Spline, and intensity-based) were applied to these ten virtual phantoms. The images, ROIs, and doses were mapped from pCT to srCT using the DVFs computed by these three DIRs and compared to those mapped using the reference DVF. Results: The average Dice coefficients for selected ROIs were from 0.85 to 0.96 for Demons, from 0.86 to 0.97 for intensity-based, and from 0.76 to 0.95 for B-Spline. The average Hausdorff distances for selected ROIs were from 2.2 to 5.4 mm for Demons, from 2.3 to 6.8 mm for intensity-based, and from 2.4 to 11.4 mm for B-Spline. The average absolute dose errors for selected ROIs were from 0.2 to 0.6 Gy for Demons, from 0.1 to 0.5 Gy for intensity-based, and from 0.5 to 1.5 Gy for B-Spline. Conclusion: Virtual phantoms were modeled after patients with lung cancer and were clinically relevant for adaptive radiotherapy treatment replanning. Virtual phantoms with known DVFs serve as references and can provide a fair comparison when evaluating different DIRs. Demons and intensity-based DIRs were shown to have smaller geometric and dosimetric uncertainties than B-Spline. Z Shen: None; K Bzdusek: an employee of Philips Healthcare; J Greskovich: None; P Xia
Purpose: To evaluate geometric and dosimetric uncertainties of CT-CBCT deformable image registration (DIR) algorithms using digital phantoms generated from real patients. Methods: We selected ten H&N cancer patients with adaptive IMRT. For each patient, a planning CT (CT1), a replanning CT (CT2), and a pretreatment CBCT (CBCT1) were used as the basis for digital phantom creation. Manually adjusted meshes were created for selected ROIs (e.g. PTVs, brainstem, spinal cord, mandible, and parotids) on CT1 and CT2. The mesh vertices were input into a thin-plate spline algorithm to generate a reference displacement vector field (DVF). The reference DVF was applied to CBCT1 to create a simulated mid-treatment CBCT (CBCT2). The CT-CBCT digital phantom consisted of CT1 and CBCT2, which were linked by the reference DVF. Three DIR algorithms (Demons, B-Spline, and intensity-based) were applied to these ten digital phantoms. The images, ROIs, and volumetric doses were mapped from CT1 to CBCT2 using the DVFs computed by these three DIRs and compared to those mapped using the reference DVF. Results: The average Dice coefficients for selected ROIs were from 0.83 to 0.94 for Demons, from 0.82 to 0.95 for B-Spline, and from 0.67 to 0.89 for intensity-based DIR. The average Hausdorff distances for selected ROIs were from 2.4 to 6.2 mm for Demons, from 1.8 to 5.9 mm for B-Spline, and from 2.8 to 11.2 mm for intensity-based DIR. The average absolute dose errors for selected ROIs were from 0.7 to 2.1 Gy for Demons, from 0.7 to 2.9 Gy for B- Spline, and from 1.3 to 4.5 Gy for intensity-based DIR. Conclusion: Using clinically realistic CT-CBCT digital phantoms, Demons and B-Spline were shown to have similar geometric and dosimetric uncertainties while intensity-based DIR had the worst uncertainties. CT-CBCT DIR has the potential to provide accurate CBCT-based dose verification for H&N adaptive radiotherapy. Z Shen: None; K Bzdusek: an employee of Philips Healthcare; S Koyfman: None; P Xia
Deformable Registration of Digital Images
管伟光; 解林; 等
1998-01-01
is paper proposes a novel elastic model and presents a deformable registration method based on the model.The method registers images without the need to extract reatures from the images,and therefore works directly on grey-level images.A new similarity metric is given on which the formation of external forces is based.The registration method,taking the coarse-to-fine strategy,constructs external forces in larger scales for the first few iterations to rely more on global evidence,and ther in smaller scales for later iterations to allow local refinements.The stiffness of the elastic body decreases as the process proceeds.To make it widely applicable,the method is not restricted to any type of transformation.The variations between images are thought as general free-form deformations.Because the elastic model designed is linearized,it can be solved very efficiently with high accuracy.The method has been successfully tested on MRI images.It will certainly find other uses such as matching time-varying sequences of pictures for motion analysis,fitting templates into images for non-rigid object recognition,matching stereo images for shape recovery,etc.
Exact surface registration of retinal surfaces from 3-D optical coherence tomography images.
Lee, Sieun; Lebed, Evgeniy; Sarunic, Marinko V; Beg, Mirza Faisal
2015-02-01
Nonrigid registration of optical coherence tomography (OCT) images is an important problem in studying eye diseases, evaluating the effect of pharmaceuticals in treating vision loss, and performing group-wise cross-sectional analysis. High dimensional nonrigid registration algorithms required for cross-sectional and longitudinal analysis are still being developed for accurate registration of OCT image volumes, with the speckle noise in images presenting a challenge for registration. Development of algorithms for segmentation of OCT images to generate surface models of retinal layers has advanced considerably and several algorithms are now available that can segment retinal OCT images into constituent retinal surfaces. Important morphometric measurements can be extracted if accurate surface registration algorithm for registering retinal surfaces onto corresponding template surfaces were available. In this paper, we present a novel method to perform multiple and simultaneous retinal surface registration, targeted to registering surfaces extracted from ocular volumetric OCT images. This enables a point-to-point correspondence (homology) between template and subject surfaces, allowing for a direct, vertex-wise comparison of morphometric measurements across subject groups. We demonstrate that this approach can be used to localize and analyze regional changes in choroidal and nerve fiber layer thickness among healthy and glaucomatous subjects, allowing for cross-sectional population wise analysis. We also demonstrate the method's ability to track longitudinal changes in optic nerve head morphometry, allowing for within-individual tracking of morphometric changes. This method can also, in the future, be used as a precursor to 3-D OCT image registration to better initialize nonrigid image registration algorithms closer to the desired solution. PMID:25312906
Sequential Non-Rigid Structure from Motion Using Physical Priors.
Agudo, Antonio; Moreno-Noguer, Francesc; Calvo, Begona; Montiel, Jose M Martinez
2016-05-01
We propose a new approach to simultaneously recover camera pose and 3D shape of non-rigid and potentially extensible surfaces from a monocular image sequence. For this purpose, we make use of the Extended Kalman Filter based Simultaneous Localization And Mapping (EKF-SLAM) formulation, a Bayesian optimization framework traditionally used in mobile robotics for estimating camera pose and reconstructing rigid scenarios. In order to extend the problem to a deformable domain we represent the object's surface mechanics by means of Navier's equations, which are solved using a Finite Element Method (FEM). With these main ingredients, we can further model the material's stretching, allowing us to go a step further than most of current techniques, typically constrained to surfaces undergoing isometric deformations. We extensively validate our approach in both real and synthetic experiments, and demonstrate its advantages with respect to competing methods. More specifically, we show that besides simultaneously retrieving camera pose and non-rigid shape, our approach is adequate for both isometric and extensible surfaces, does not require neither batch processing all the frames nor tracking points over the whole sequence and runs at several frames per second. PMID:27046840
Testing for additivity with B-splines
2007-01-01
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.
Testing for additivity with B-splines
Heng-jian CUI; Xu-ming HE; Li LIU
2007-01-01
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 kernelbased tests. The score test is found to have a good overall performance.
Registration accuracy of 153Gd transmission images of the brain
The aim of the study was to determine the accuracy of non-rigid nine-parameter image registrations based on 153Gd transmission computed tomography (TCT) images as compared with those based on 99mTc-ethyl cysteinate dimer (ECD) images and to assess whether normalised mutual information (NMI) or count difference (CD) should be used. TCT and ECD data were acquired in 25 randomly selected patients. Emission images were registered to an ECD template with a CD cost function. The same registration parameters were applied to the transmission images to create a TCT template. All TCT images were registered to the TCT template and the same registration parameters were applied to the ECD images. The procedure was repeated with NMI as cost function. Accuracy of both ECD-based and TCT-based registrations was assessed by comparing the normalisation parameter values and regional activities in the spatially normalised ECD images, using a mixed-model analysis of variance (ANOVA). Scheffe post hoc tests were performed. No significant differences were found between ECD/CD, ECD/NMI and TCT/CD, suggesting that ECD registration can be done with either CD or NMI, and that TCT registration using CD is equally as accurate as ECD registration. The accuracy of TCT registration with NMI was lower, with discrepancies occurring in the frontal inferior region and the cerebellum. The analysis of normalisation parameters indicated that z-scaling is underestimated and yz-rotation overestimated with TCT/NMI registration. We conclude that ECD registrations with CD or NMI are as accurate as TCT registrations with CD and that TCT registrations with NMI should be avoided. (orig.)
Non-rigid connector in fixed partial dentures with pier abutment: An enigma simplified
Prince Kumar
2012-01-01
Full Text Available In the fixed partial denture (FPD prosthesis with rigid connectors, forces of mastication are transmitted to the terminal retainers, which make the middle abutment act as a class I lever fulcrum. This techno-clinical cycle often causes failure of the fixed partial dentures. To overcome this dilemma, a non-rigid connector may be incorporated on the distal aspect of the middle (pier abutment. The non-rigid connector counterbalances the effects of these forces of leverage. This clinical report presents the use of a non-rigid connector in a long-span, five-unit FPD, replacing two missing teeth with an intermediate pier abutment.
Higher Order Kernels and Locally Affine LDDMM Registration
Sommer, Stefan; Darkner, Sune; Pennec, Xavier
2011-01-01
To achieve sparse description that allows intuitive analysis, we aim to represent deformation with a basis containing interpretable elements, and we wish to use elements that have the description capacity to represent the deformation compactly. We accomplish this by introducing higher order kernels in the LDDMM registration framework. The kernels allow local description of affine transformations and subsequent compact description of non-translational movement and of the entire non-rigid deformation. This is obtained with a representation that contains directly interpretable information from both mathematical and modeling perspectives. We develop the mathematical construction behind the higher order kernels, we show the implications for sparse image registration and deformation description, and we provide examples of how the capacity of the kernels enables registration with a very low number of parameters. The capacity and interpretability of the kernels lead to natural modeling of articulated movement, and th...
Supplementary active stabilization of nonrigid gravity gradient satellites
Keat, J. E.
1972-01-01
The use of active control for stability augmentation of passive gravity gradient satellites is investigated. The reaction jet method of control is the main interest. Satellite nonrigidity is emphasized. The reduction in the Hamiltonian H is used as a control criteria. The velocities, relative to local vertical, of the jets along their force axes are shown to be of fundamental significance. A basic control scheme which satisfies the H reduction criteria is developed. Each jet is fired when its velocity becomes appropriately large. The jet is de-energized when velocity reaches zero. Firing constraints to preclude orbit alteration may be needed. Control is continued until H has been minimized. This control policy is investigated using impulse and rectangular pulse models of the jet outputs.
Automatic Detection of Wild-type Mouse Cranial Sutures
Ólafsdóttir, Hildur; Darvann, Tron Andre; Hermann, Nuno V.;
, automatic detection of the cranial sutures becomes important. We have previously built a craniofacial, wild-type mouse atlas from a set of 10 Micro CT scans using a B-spline-based nonrigid registration method by Rueckert et al. Subsequently, all volumes were registered nonrigidly to the atlas. Using these...... observer traced the sutures on each of the mouse volumes as well. The observer outperforms the automatic approach by approximately 0.1 mm. All mice have similar errors while the suture error plots reveal that suture 1 and 2 are cumbersome, both for the observer and the automatic approach. These sutures can...
3D ultrasound-CT registration of the liver using combined landmark-intensity information
An important issue in computer-assisted surgery of the liver is a fast and reliable transfer of preoperative resection plans to the intraoperative situation. One problem is to match the planning data, derived from preoperative CT or MR images, with 3D ultrasound images of the liver, acquired during surgery. As the liver deforms significantly in the intraoperative situation non-rigid registration is necessary. This is a particularly challenging task because pre- and intraoperative image data stem from different modalities and ultrasound images are generally very noisy. One way to overcome these problems is to incorporate prior knowledge into the registration process. We propose a method of combining anatomical landmark information with a fast non-parametric intensity registration approach. Mathematically, this leads to a constrained optimization problem. As distance measure we use the normalized gradient field which allows for multimodal image registration. A qualitative and quantitative validation on clinical liver data sets of three different patients has been performed. We used the distance of dense corresponding points on vessel center lines for quantitative validation. The combined landmark and intensity approach improves the mean and percentage of point distances above 3 mm compared to rigid and thin-plate spline registration based only on landmarks. The proposed algorithm offers the possibility to incorporate additional a priori knowledge - in terms of few landmarks - provided by a human expert into a non-rigid registration process. (orig.)
Non-rigid connector in fixed partial dentures with pier abutment: An enigma simplified
Prince Kumar; Vishal Singh; Roshni Goel; Harkanwal P Singh
2012-01-01
In the fixed partial denture (FPD) prosthesis with rigid connectors, forces of mastication are transmitted to the terminal retainers, which make the middle abutment act as a class I lever fulcrum. This techno-clinical cycle often causes failure of the fixed partial dentures. To overcome this dilemma, a non-rigid connector may be incorporated on the distal aspect of the middle (pier) abutment. The non-rigid connector counterbalances the effects of these forces of leverage. This clinical report...
Goshtasby, A Ardeshir
2012-01-01
This book presents a thorough and detailed guide to image registration, outlining the principles and reviewing state-of-the-art tools and methods. The book begins by identifying the components of a general image registration system, and then describes the design of each component using various image analysis tools. The text reviews a vast array of tools and methods, not only describing the principles behind each tool and method, but also measuring and comparing their performances using synthetic and real data. Features: discusses similarity/dissimilarity measures, point detectors, feature extr
Objects: In functional stereotactic neurosurgery, one of the cornerstones upon which the success and the operating time depends is an accurate targeting. The subthalamic nucleus (STN) is the usual target involved when applying deep brain stimulation for Parkinson's disease (PD). Unfortunately, STN is usually not clearly visible in common medical imaging modalities, which justifies the use of atlas-based segmentation techniques to infer the STN location. Materials and methods: Eight bilaterally implanted PD patients were included in this study. A three-dimensional T1-weighted sequence and inversion recovery T2-weighted coronal slices were acquired pre-operatively. We propose a methodology for the construction of a ground truth of the STN location and a scheme that allows both, to perform a comparison between different non-rigid registration algorithms and to evaluate their usability to locate the STN automatically. Results: The intra-expert variability in identifying the STN location is 1.06±0.61 mm while the best non-rigid registration method gives an error of 1.80±0.62 mm. On the other hand, statistical tests show that an affine registration with only 12 degrees of freedom is not enough for this application. Conclusions: Using our validation-evaluation scheme, we demonstrate that automatic STN localization is possible and accurate with non-rigid registration algorithms. (orig.)
Mammogram CAD, hybrid registration and iconic analysis
Boucher, A.; Cloppet, F.; Vincent, N.
2013-03-01
This paper aims to develop a computer aided diagnosis (CAD) based on a two-step methodology to register and analyze pairs of temporal mammograms. The concept of "medical file", including all the previous medical information on a patient, enables joint analysis of different acquisitions taken at different times, and the detection of significant modifications. The developed registration method aims to superimpose at best the different anatomical structures of the breast. The registration is designed in order to get rid of deformation undergone by the acquisition process while preserving those due to breast changes indicative of malignancy. In order to reach this goal, a referent image is computed from control points based on anatomical features that are extracted automatically. Then the second image of the couple is realigned on the referent image, using a coarse-to-fine approach according to expert knowledge that allows both rigid and non-rigid transforms. The joint analysis detects the evolution between two images representing the same scene. In order to achieve this, it is important to know the registration error limits in order to adapt the observation scale. The approach used in this paper is based on an image sparse representation. Decomposed in regular patterns, the images are analyzed under a new angle. The evolution detection problem has many practical applications, especially in medical images. The CAD is evaluated using recall and precision of differences in mammograms.
Image registration with auto-mapped control volumes
inhale and exhale phases of a lung 4D CT. Algorithm convergence was confirmed by starting the registration calculations from a large number of initial transformation parameters. An accuracy of ∼2 mm was achieved for both deformable and rigid registration. The proposed image registration method greatly reduces the complexity involved in the determination of homologous control points and allows us to minimize the subjectivity and uncertainty associated with the current manual interactive approach. Patient studies have indicated that the two-step registration technique is fast, reliable, and provides a valuable tool to facilitate both rigid and nonrigid image registrations
The purpose of this work is to validate an in-house deformable image registration (DIR) algorithm for adaptive radiotherapy for head and neck patients. We aim to use the registrations to estimate the ''dose of the day'' and assess the need to replan. NiftyReg is an open-source implementation of the B-splines deformable registration algorithm, developed in our institution. We registered a planning CT to a CBCT acquired midway through treatment for 5 HN patients that required replanning. We investigated 16 different parameter settings that previously showed promising results. To assess the registrations, structures delineated in the CT were warped and compared with contours manually drawn by the same clinical expert on the CBCT. This structure set contained vertebral bodies and soft tissue. Dice similarity coefficient (DSC), overlap index (OI), centroid position and distance between structures' surfaces were calculated for every registration, and a set of parameters that produces good results for all datasets was found. We achieve a median value of 0.845 in DSC, 0.889 in OI, error smaller than 2 mm in centroid position and over 90% of the warped surface pixels are distanced less than 2 mm of the manually drawn ones. By using appropriate DIR parameters, we are able to register the planning geometry (pCT) to the daily geometry (CBCT).
A Condition Number for Non-Rigid Shape Matching
Ovsjanikov, Maks
2011-08-01
© 2011 The Author(s). Despite the large amount of work devoted in recent years to the problem of non-rigid shape matching, practical methods that can successfully be used for arbitrary pairs of shapes remain elusive. In this paper, we study the hardness of the problem of shape matching, and introduce the notion of the shape condition number, which captures the intuition that some shapes are inherently more difficult to match against than others. In particular, we make a connection between the symmetry of a given shape and the stability of any method used to match it while optimizing a given distortion measure. We analyze two commonly used classes of methods in deformable shape matching, and show that the stability of both types of techniques can be captured by the appropriate notion of a condition number. We also provide a practical way to estimate the shape condition number and show how it can be used to guide the selection of landmark correspondences between shapes. Thus we shed some light on the reasons why general shape matching remains difficult and provide a way to detect and mitigate such difficulties in practice.
Dealing with difficult deformations: Construction of a knowledge-based deformation atlas
Thorup, Signe Strann; Darvann, T.A.; Hermann, N.V.; Larsen, P.; Olafsdottir, H.; Paulsen, Rasmus Reinhold; Kane, A.A.; Govier, D.; Lo, L.J.; Kreiborg, S.; Larsen, Rasmus
2010-01-01
(atlases) for the pre- and post-surgical populations, respectively, were automatically constructed by non-rigid registration. An expert placed corresponding landmarks in the cleft area in the two atlases; this provided prior information used to build a knowledge-based deformation atlas. We model the change...... from pre- to post-surgery using thin-plate spline warping. The registration results are convincing and represent a first move towards an automatic registration method for dealing with difficult deformations due to this type of surgery. New or breakthrough work to be presented: The method provides a...... needed. We have previously demonstrated that non-rigid registration using B-splines is able to provide automated determination of point correspondences in populations of infants without cleft lip. However, this type of registration fails when applied to the task of determining the complex deformation...
Elastic image registration via rigid object motion induced deformation
Zheng, Xiaofen; Udupa, Jayaram K.; Hirsch, Bruce E.
2011-03-01
In this paper, we estimate the deformations induced on soft tissues by the rigid independent movements of hard objects and create an admixture of rigid and elastic adaptive image registration transformations. By automatically segmenting and independently estimating the movement of rigid objects in 3D images, we can maintain rigidity in bones and hard tissues while appropriately deforming soft tissues. We tested our algorithms on 20 pairs of 3D MRI datasets pertaining to a kinematic study of the flexibility of the ankle complex of normal feet as well as ankles affected by abnormalities in foot architecture and ligament injuries. The results show that elastic image registration via rigid object-induced deformation outperforms purely rigid and purely nonrigid approaches.
Automatic Mrf-Based Registration of High Resolution Satellite Video Data
Platias, C.; Vakalopoulou, M.; Karantzalos, K.
2016-06-01
In this paper we propose a deformable registration framework for high resolution satellite video data able to automatically and accurately co-register satellite video frames and/or register them to a reference map/image. The proposed approach performs non-rigid registration, formulates a Markov Random Fields (MRF) model, while efficient linear programming is employed for reaching the lowest potential of the cost function. The developed approach has been applied and validated on satellite video sequences from Skybox Imaging and compared with a rigid, descriptor-based registration method. Regarding the computational performance, both the MRF-based and the descriptor-based methods were quite efficient, with the first one converging in some minutes and the second in some seconds. Regarding the registration accuracy the proposed MRF-based method significantly outperformed the descriptor-based one in all the performing experiments.
Real-time Animation Technique for a Kind of Non-rigid Objects
无
2003-01-01
A real-time animation technique for a kind of non-rigid objects, flexible and thin objects, is proposed, which can update with stability the state of n mass points of the mass-spring (MS) model with time complexity of O(n). The new implicit numerical integration technique of the authors, which is based on a simple approximation of the linear system, has great advantages over the existing implicit integration methods. Moreover, experiment shows that the new technique is highly efficient in animating a kind of non-rigid objects, and suitable for the draping module of the 3D garment CAD system.
Augmenting Images of Non-Rigid Scenes Using Point and Curve Correspondences
Bartoli, Adrien; Von Tunzelmann, Eugénie; Zisserman, Andrew
2004-01-01
Our goal is to augment images of non-rigid scenes coming from single-camera footage. We do not assume any a priori information about the scene being viewed, such as for example a parameterized 3D model or the motion of the camera. One possible solution is to use non-rigid factorization of points, from which a dense interpolating function modeled by a thin-plane spline can be computed. However, in many cases, point correspondences fail to capture precisely all the deformations occurring in the...