Biomedical subjects

P Cerveri

Publications and source records attributed to P Cerveri.

16 recordsLinked to original sources

PubMed2005

CT-3D rotational angiography automatic registration: a sensitivity analysis.

Preprocessing, binning and dataset subsampling are investigated with regard to simultaneous maximisation of the speed, accuracy and robustness of CT-3D rotational angiography (3DRA) registration. Clinical diagnosis and treatment can both take advantage of this integration, because 3DRA allows the shape of vessel structures to be evaluated three-dimensionally with respect to standard 2D projective angiography. The method for optimising preprocessing, binning and subsampling consisted of independent variation of the corresponding parameters to maximise robustness and speed while maintaining subvoxel accuracy; the latter was computed as the sum of the mean squared errors initially present in the registrations with the errors relative to both binning and subsampling. The results suggest the choice of 256 bins, steps between 14 mm (coarse optimisation) and 2.5 mm (fine optimisation) and bone segmentation by threshold, for binning, subsampling and preprocessing, respectively. The application of this parameter set-up to 50 CT-3DRA registrations resulted in a saving, on average, of 40% of the time with respect to the method previously used, while registration error was maintained within 2 mm (1.97 mm, 90% confidence interval) and robustness was increased, so that no manual initial realignment was needed in 48 registrations. Validation by the registration of images acquired for a head phantom showed subvoxel residual errors. In conclusion, the proposed procedure can be considered a satisfactory strategy to optimise CT-3DRA registration.

Angiography↗

PubMed2005

Derivation of centers and axes of rotation for wrist and fingers in a hand kinematic model: methods and reliability results.

In the field of 3D reconstruction of human motion from video, model-based techniques have been proposed to increase the estimation accuracy and the degree of automation. The feasibility of this approach is strictly connected with the adopted biomechanical model. Particularly, the representation of the kinematic chain and the assessment of the corresponding parameters play a relevant role for the success of the motion assessment. In this paper, the focus is on the determination of the kinematic parameters of a general hand skeleton model using surface measurements. A novel method that integrates nonrigid sphere fitting and evolutionary optimization is proposed to estimate the centers and the functional axes of rotation of the skeletal joints. The reliability of the technique is tested using real movement data and simulated motions with known ground truth 3D measurement noise and different ranges of motion (RoM). With respect to standard nonrigid sphere fitting techniques, the proposed method performs 10-50% better in the best condition (very low noise and wide RoM) and over 100% better with physiological artifacts and RoM. Repeatability in the range of a couple of millimeters, on the localization of the centers of rotation, and in the range of one degree, on the axis directions is obtained from real data experiments.

Adult↗

PubMed2005

Preliminary study on the use of nonrigid registration for thoraco-abdominal radiosurgery.

The inclusion of organ deformation and movement in radiosurgery treatment planning is of increasing importance as research and clinical applications begin to take into consideration the effects of physiological processes, like breathing, on the shape and position of lesions. In this scenario, the challenge is to localize the target in toto (not only by means of marker sampling) and to calculate the dose distribution as the sum of all the contributions from the positions assumed by the target during the respiratory cycle. The aim of this work is to investigate the use of nonrigid registration for target tracking and dynamic treatment planning, i.e., treatment planning based not on one single CT scan but on multiple CT scans representative of the respiration. Twenty patients were CT scanned at end-inhale and end-exhale. An expert radiation oncologist identified the PTV in both examinations. The two CT data sets per patient were nonrigidly registered using a free-form deformation algorithm based on B-splines. The optimized objective function consisted of a weighted sum of a similarity criterion (Mutual Information) and a regularization factor which constrains the transformation to be locally rigid. Once the transformation was obtained and the registration validated, its parameters were applied to the target only. Finally, the deformed target was compared to the PTV delineated by the radiation oncologist in the other study. The results of this procedure show an agreement between the center of mass as well as volume of the target identified automatically by deformable registration and manually by the radiation oncologist. Moreover, obtained displacements were in agreement with body structure constraints and considerations usually accepted in radiation therapy practice. No significant influence of initial target volume on displacements was found. In conclusion, the proposed method seems to offer the possibility of using nonrigid registrations in radiosurgery treatment planning, even if more cases need to be investigated in order to give a statistical consistency to parameter setup and proposed considerations.

Abdominal Neoplasms↗

PubMed2004

Kinematical models to reduce the effect of skin artifacts on marker-based human motion estimation.

The estimation of the skeletal motion obtained from marker-based motion capture systems is known to be affected by significant bias caused by skin movement artifacts, which affects joint center and rotation axis estimation. Among different techniques proposed in the literature, that based on rigid body model, still the most used by commercial motion capture systems, can smooth only part of the above effects without eliminating their main components. In order to sensibly improve the accuracy of the motion estimation, a novel technique, named local motion estimation (LME), is proposed. This rests on a recently described approach that, using virtual humans and extended Kalman filters, estimates the kinematical variables directly from 2D measurements without requiring the 3D marker reconstruction. In this paper, we show how such method can be extended to include the computation of the local marker displacement due to skin artifacts. The 3D marker coordinates, expressed in the corresponding local reference coordinate frames, are inserted into the state vector of the filter and their dynamics is automatically estimated, with adequate accuracy, without assuming any particular deformation function. Simulated experiments of lower limb motion, involving systematic mislocations (5, 10, 20 mm) and random errors of the marker coordinates and joint center locations (+/-5, +/-10, +/-15 mm), have shown that artifact motion can be substantially decoupled from the global skeletal motion with an effective increase of the accuracy wrt standard techniques. In particular, the comparison between the nominal kinematical variables and the one recovered from markers attached to the skin surface proved LME to be sensibly superior (50% in the worse condition) to the methods imposing marker-bone rigidity. In conclusion, while requiring further validation on real movement data, we argue that the proposed method can constitute an appropriate approach toward the improvement of the human motion estimation.

Biomechanical Phenomena↗

PubMed2004

Evolutionary optimization for robust hierarchical computation of the rotation centres of kinematic chains from reduced ranges of motion the lower spine case.

A novel technique based on evolutionary optimization is proposed here to compute the average rotation centres (RCs) of ball joints linked into kinematic chains using 3D trajectories of the markers attached to the external surface of the corresponding articulated structures. The chain is hierarchically solved by iteratively minimizing the variance of the marker distances from the actual RC through an evolutional strategy method (ESM) from proximal to distal joints. In particular, the technique is compared to the non-rigid sphere-fitting method, recently proposed in literature and implemented through a closed-form solution (CFS), in conditions of random and systematic noise superimposed to the marker coordinates. Results from simulated motions showed that, in case of small range of motion (5 degrees , 10 degrees ) the performance of CFS is really unreliable whereas ESM provided satisfactory accuracy. Error propagation along the kinematic chain was found to be negligible. Also in the case of systematic errors, ESM provides an accuracy that is sensibly better than that of the CFS. As a case study, ESM was applied to the in vivo computation of the RCs of the vertebrae in the lower spine region using a specific marker protocol. A set of spine movements by a normal adult male, recorded by an optoelectronic motion capture system, were processed with the developed method. The variability of the estimated average RCs was small (few millimeters) in agreement with the literature data from cadaveric studies and X-ray imaging.

Algorithms↗

PubMed2004

Non-invasive approach towards the in vivo estimation of 3D inter-vertebral movements: methods and preliminary results.

A kinematical model of the lower spine was designed and used to obtain a robust estimation of the vertebral rotations during torso movements from skin-surface markers recorded by video-cameras. Markers were placed in correspondence of the anatomical landmarks of the pelvic bone and vertebral spinous and transverse processes, and acquired during flexion, lateral bending and axial motions. In the model calibration stage, a motion-based approach was used to compute the rotation axes and centres of the functional segmental units. Markers were mirrored into virtual points located on the model surface, expressed in the local reference system of coordinates. The spine motion assessment was solved into the domain of extended Kalman filters: at each frame of the acquisition, the model pose was updated by minimizing the distances between the measured 2D marker projections on the cameras and the corresponding back-projections of virtual points located on the model surface. The novelty of the proposed technique rests on the fact that the varying location of the rotation centres of the functional segmental units can be tracked directly during motion computation. In addition, we show how the effects of skin artefacts on orientation data can be taken into account. As a result, the kinematical estimation of simulated motions shows that orientation artefacts were reduced by a factor of at least 50%. Preliminary experiments on real motion confirmed the reliability of the proposed method with results in agreement with classical studies in literature.

Algorithms↗

PubMed2003

GAAS: gene array analyzer software for management, analysis and visualization of gene expression data.

SUMMARY: GAAS, Gene Array Analyzer Software supports multi-user efficient management and suitable analyses of large amounts of gene expression data across replicated experiments. Its management framework handles input data generated by different technologies. A multi-user environment allows each user to store his/her own data visualization scheme, analysis parameters used, values and formats of the output data. The analysis engine performs: background and spot quality evaluation, data normalization, differential gene expression analyses in single and multiple replica experiments. Results of expression profiles can be interactively navigated through graphical interfaces and stored into output databases.

Algorithms↗

PubMed2003

Real-time human motion estimation using biomechanical models and non-linear state-space filters.

In the field of sports biomechanics and rehabilitation engineering, the possibility of computing, in real time, the angular displacements and derivatives of human joints, from a video of motion sequences, represents an appealing goal. In particular, applications of biofeedback protocols in rehabilitation can benefit from this capability. The focus of the investigation was concerned with the application of biomechanical models, comprising of a kinematic chain and surface envelopes, and state-space filters, to the computation, in real time and with high accuracy, of the angular data and derivatives. By minimising the distances, measured with TV cameras, between the 2D marker projections and the corresponding back-projected markers located on the mannequin, the configuration of the biomechanical model was automatically updated. The use of state-space estimation allowed the computation of smooth derivatives of the orientation data. Owing to the non-linearity of the functions involved, the derivatives of the observation model were obtained through a multidimensional extension of Stirling's interpolation formula. Proper algorithms were developed to cope with the model calibration, initialisation and data labelling. Extensive experiments on real and simulated motions proved the reliability (maximum angular error less than 1 degree, maximum point reconstruction less than 1 mm) of the developed system, which is robust to false matching caused by marker occlusions. Moreover, orientation artifacts due to skin motion can be reduced by a factor of 50%.

Algorithms↗

PubMed2003

Hierarchical radial basis function networks and local polynomial un-warping for X-ray image intensifier distortion correction: a comparison with global techniques.

Global polynomial (GP) methods have been widely used to correct geometric image distortion of small-size (up to 30 cm) X-ray image intensifiers (XRIIs). This work confirms that this kind of approach is suitable for 40 cm XRIIs (now increasingly used). Nonetheless, two local methods, namely 3rd-order local un-warping polynomials (LUPs) and hierarchical radial basis function (HRBF) networks are proposed as alternative solutions. Extensive experimental tests were carried out to compare these methods with classical low-order local polynomial and GP techniques, in terms of residual error (RMSE) measured at points not used for parameter estimation. Simulations showed that the LUP and HRBF methods had accuracies comparable with that attained using GP methods. In detail, the LUP method (0.353 microm) performed worse than HRBF (0.348 microm) only for small grid spacing (15 x 15 control points); the accuracy of both HRBF (0.157 microm) and LUP (0.160 microm) methods was little affected by local distortions (30 x 30 control points); weak local distortions made the GP method poorer (0.320 microm). Tests on real data showed that LUP and HRBF had accuracies comparable with that of GP for both 30 cm (GP: 0.238 microm; LUP: 0.240 microm; HRBF: 0.238 microm) and 40 cm (GP: 0.164 microm; LUP: 0.164 microm; HRBF: 0.164 microm) XRIIs. The LUP-based distortion correction was implemented in real time for image correction in digital tomography applications.

Algorithms↗

PubMed2003

Robust recovery of human motion from video using Kalman filters and virtual humans.

In sport science, as in clinical gait analysis, optoelectronic motion capture systems based on passive markers are widely used to recover human movement. By processing the corresponding image points, as recorded by multiple cameras, the human kinematics is resolved through multistage processing involving spatial reconstruction, trajectory tracking, joint angle determination, and derivative computation. Key problems with this approach are that marker data can be indistinct, occluded or missing from certain cameras, that phantom markers may be present, and that both 3D reconstruction and tracking may fail. In this paper, we present a novel technique, based on state space filters, that directly estimates the kinematical variables of a virtual mannequin (biomechanical model) from 2D measurements, that is, without requiring 3D reconstruction and tracking. Using Kalman filters, the configuration of the model in terms of joint angles, first and second order derivatives is automatically updated in order to minimize the distances, as measured on TV-cameras, between the 2D measured markers placed on the subject and the corresponding back-projected virtual markers located on the model. The Jacobian and Hessian matrices of the nonlinear observation function are computed through a multidimensional extension of Stirling's interpolation formula. Extensive experiments on simulated and real data confirmed the reliability of the developed system that is robust against false matching and severe marker occlusions. In addition, we show how the proposed technique can be extended to account for skin artifacts and model inaccuracy.

Humans↗

PubMed2002

Distortion correction for x-ray image intensifiers: local unwarping polynomials and RBF neural networks.

In this paper we present two novel techniques, namely a local unwarping polynomial (LUP) and a hierarchical radial basis function (HRBF) network, to correct geometric distortions in XRII images. The two techniques have been implemented and compared, in terms of residual error measured at control and intermediate points, with local and global methods reported in the previous literature. In particular, LUP rests on a locally optimized 3rd degree polynomial applied within each quadrilateral cell on the rectilinear calibration grid of points. HRBF, based on a feed-forward neural network paradigm, is constituted by a set of hierarchical layers at increasing cut-off frequency, each characterized by a set of Gaussian functions. Extensive experiments have been performed both on simulated and real data. In simulation, we tested the effect of pincushion, sigmoidal and local distortions, along with the number of calibration points. Provided that a sufficient number of cells of the calibration grid is available, the obtained accuracy for both LUP and HRBF is comparable to or better than that of global polynomial technique. Tests on real data, carried out by using two different (12 in. and 16 in.) XRIIs, showed that the global polynomial accuracy (0.16+/-0.08 pixels) is slightly worse than that of LUP (0.07+/-0.05 pixels) and HRBF (0.08+/-0.04 pixels). The effects of the discontinuity at the border of the local areas and the decreased accuracy at intermediate points, typical of local techniques, have been proved to be smoothed for both LUP and HRBF.

Algorithms↗

PubMed2001

Symbolic representation of anatomical knowledge: concept classification and development strategies.

In this paper a novel approach to anatomy knowledge representation is described. The focus of the present research has been on the development of a representational framework where the conceptual level has been implemented by using hierarchical and nonhierarchical conceptual networks. This has allowed handling the demand for multiple views of anatomy (systemic and topographical views). The terminological level of the knowledge representation has been implemented by using a compositional strategy which has avoided the explicit storage of the terms used to express composite concepts. Hierarchical relations and composite concept representations have required supervision of both the inheritance and concept reconstruction. For this purpose heuristic knowledge has been stored in terms of consistency rules in the knowledge base. As proof of the capability of this system, we show how the knowledge base has been used to provide symbolic access to spatial information consisting of a reduced set of images from the Visible Human Dataset.

Anatomy, Cross-Sectional↗

PubMed2001

A fast method for calibrating video-based motion analysers using only a rigid bar.

Video-camera systems are widely used in biomechanics and clinical fields to measure the 3D kinematic measurements of human motion. To be used, they need to be calibrated, that is the parameters which geometrically define the cameras have to be determined. It is shown here how this can be achieved by surveying a rigid bar in motion inside the working volume, and in a very short time: less than 15 s on a Pentium III. The exterior parameters are estimated through the coplanarity constraint, the camera focal lengths through the properties of epipolar geometry and the principal points with a fast evolutionary optimisation which guarantees convergence when the initial principal points cannot be adequately estimated. The method has been widely tested on simulated and real data. Results show that its accuracy is comparable with that obtained using methods based on points of known 3D coordinates (DLT): 0.37 mm RMS error over a volume with a diagonal approximately 1.5 m. A preferential absolute reference system is obtained from the same bar motion data and is used to guide an intelligent decimation of the data. Finally, the role that the principal points play in achieving a high accuracy, which is questioned in the computer vision domain, is assessed through simulations.

Biomechanical Phenomena↗

PubMed2000

Remote access to anatomical information: an integration between semantic knowledge and visual data.

A novel internet-based application is presented which provides access to anatomy knowledge through symbolic modality expressed by keywords taken from controlled or non-controlled terminology. The system is based on a database where anatomical concepts have been organized into a hierarchical framework. Along with term queries that allow retrieving concepts containing or exactly matching the used keyword, the system also provides semantic access to anatomical information. Queries can be setup, which retrieve concepts relying to a particular meaning and sharing a particular relationship. Moreover, the application has the capability to refine the search of the terms by querying the UMLS knowledge server. Anatomical image data have been integrated by using Visible Human Dataset. A set of these images has been indexed according to our anatomical classification and is used inside the application. The system has been implemented through Java client-server technology and works within standard Internet browsers.

Anatomy↗

PubMed1998

Complete calibration of a stereo photogrammetric system through control points of unknown coordinates.

This paper presents a new method for calibrating a video 3D stereo-photogrammetric system. The external parameters and the focal lengths of the cameras are determined from the epipolar constraint and the principal points are computed through the minimisation of a cost function carried out through an evolutionary optimisation. The method has been made more robust with a deterministic annealing procedure of the search region amplitude. Calibration is carried out by moving a rigid bar, carrying two markers on its extremities, inside the working volume. The distance between the two markers is the only measure required. Tests on real data are reported which show that the obtained accuracy is comparable to the one achieved calibrating with control points of known 3D coordinates.

Calibration↗

PubMed1997

Statistical comparison of DLT versus ILSSC in the calibration of a photogrammetric stereo-system.

This paper compares the DLT and ILSSC approaches in the geometrical calibration of a photogrammetric stereo-system in terms of accuracy and speed. To come up with an unbiased quantitative evaluation of the accuracy of the algorithms, the concept of reliable estimate has been introduced: the statistical distribution of the accuracy is assessed over different calibration experiments performed with the same data but with different noise distribution and different test sets. Results show that in the simulations where the only error on the two-dimensional points was Gaussian, zero mean, and on real data which were corrected for distortions through polynomial or linear interpolation, the accuracy of the two methods was quite similar. DLT showed more accurate than ILSSC on simulated data with residual distortion errors and on real data which were not corrected for distortions. As far as speed is concerned, a fast triangulation algorithm is associated with ILSSC while the simultaneous solution of two pairs of DLT equations is associated to DLT. The first algorithm is much faster, requiring 113 flops per point versus 259 of DLT; the fast triangulation with DLT parameters does not achieve the same accuracy on the reconstructed three-dimensional position. Taken all together the results suggest that ILSSC can be theoretically considered the best approach to three-dimensional reconstruction, provided that distortions are corrected in advance. The statistical evaluation of the accuracy allows a fair judgement of the performances of the algorithms to be obtained, unbiased by particular distributions of measurement errors and test points.

Algorithms↗