Assessment of screw displacement axis accuracy and repeatability for joint kinematic description using an electromagnetic tracking device

Screw displacement axes (SDAs) have been employed to describe joint kinematics in biomechanical studies. Previous reports have investigated the accuracy of SDAs combining various motion analysis techniques and smoothing procedures. To our knowledge, no study has assessed SDA accuracy describing the relative movement between adjacent bodies with an electromagnetic tracking system. This is important, since in relative motion, neither body is fixed and consequently sensitivity to potential measurement errors from both bodies may be significant. Therefore, this study assessed the accuracy of SDAs for describing relative motion between two moving bodies. We analyzed numerical simulated data, and physical experimental data recorded using a precision jig and electromagnetic tracking device. The numerical simulations demonstrated SDA position accuracy (p=0.04) was superior for single compared to relative body motion, whereas orientation accuracy (p=0.2) was similar. Experimental data showed data-filtering (Butterworth filter) improved SDA position and orientation accuracies for rotation magnitudes smaller or equal to 5.0 degrees, with no effect at larger rotation magnitudes (p<0.05). This suggests that in absence of a filter, SDAs should only be calculated at rotations of greater than 5.0 degrees. For rotation magnitudes of 0.5 degrees (5.0 degrees ) about the SDA, SDA position and orientation error measurements determined from filtered experimental data were 3.75+/-0.30 mm (3.31+/-0.21 mm), and 1.10+/-0.04 degrees (1.04+/-0.03 degrees ), respectively. Experimental accuracy values describing the translation along and rotation about the SDA, were 0.06+/-0.00 mm and 0.09+/-0.01 degrees, respectively. These small errors establish the capability of SDAs to detect small translations, and rotations. In conclusion, application of SDAs should be a useful tool for describing relative motion in joint kinematic studies.     

Evaluation of an electromagnetic position tracking device for measuring in vivo, dynamic joint kinematics

An electromagnetic position tracking device was evaluated to determine its static and dynamic accuracy and reliability for applications related to measuring in vivo joint kinematics. The device detected the position and orientation of small coiled sensors, maintained in an electromagnetic field. System output was measured against known translations or rotations throughout the measurement volume. Average translational errors during static testing were 0.1 +/- 0.04, 0.2 +/- 0.17, and 0.8 +/- 0.81 mm (mean+/-SD) for sensors 50, 300, and 550 mm away from the field generator, respectively. Average rotational errors were 0.4 +/- 0.31 degrees, 0.4 +/- 0.21 degrees, and 0.9 +/- 0.85 degrees (mean +/- SD) for sensors located at the same distances. Since we intended to use this system in an animal walking on a treadmill, we incrementally moved the sensors under various treadmill conditions. The effects of treadmill operation on translational accuracy were found to be negligible. The effects of dynamic motions on sensor-to-sensor distance were also assessed for future data collection in the animal. Sensor-to-sensor distance showed standard deviations of 2.6 mm and a range of 13 mm for the highest frequency tested (0.23 Hz). We conclude that this system is useful for static or slow dynamic motions, but is of limited use for obtaining gait kinematics at higher speeds.     

Repeatability of 3D gait kinematics obtained from an electromagnetic tracking system during treadmill locomotion

The purpose of this paper was to describe a technique that enables three-dimensional (3D) gait kinematics to be obtained using an electromagnetic tracking system, and to report the intra-trial, intra-day/inter-tester and inter-day/intra-tester repeatability of kinematic gait data obtained using this technique. Ten able-bodied adults underwent four gait assessments; the same two testers tested each subject independently on two different days. Gait assessments were conducted on a custom-built long-bed treadmill with no metal components between the rollers. Each gait assessment involved familiarisation to treadmill walking, subject anatomical and functional calibration, and a period of steady-state treadmill walking at a self-selected speed. Following data collection, 3D joint kinematics were calculated using the joint coordinate system approach. 3D joint angle waveforms for 10 left and right strides were extracted and temporally normalised for each trial. Intra-trial, intra-day/inter-tester and inter-day/intra-tester repeatability of the temporally normalised kinematic waveforms were quantified using the coefficient of multiple determination (CMD). CMDs for joint kinematics averaged 0.942 intra-trial, 0.849 intra-day/inter-tester and 0.773 inter-day/intra-tester. In general, sagittal plane kinematics were more repeatable than frontal or transverse plane kinematics, and kinematics at the hip were more repeatable than at the knee or ankle. The level of repeatability of kinematic gait data obtained during treadmill walking using this protocol was equal or superior to that reported previously for overground walking using image-based protocols.     

A method for measuring joint kinematics designed for accurate registration of kinematic data to models constructed from CT data

A method for measuring three-dimensional kinematics that incorporates the direct cross-registration of experimental kinematics with anatomic geometry from Computed Tomography (CT) data has been developed. Plexiglas registration blocks were attached to the bones of interest and the specimen was CT scanned. Computer models of the bone surface were developed from the CT image data. Determination of discrete kinematics was accomplished by digitizing three pre-selected contiguous surfaces of each registration block using a three-dimensional point digitization system. Cross-registration of bone surface models from the CT data was accomplished by identifying the registration block surfaces within the CT images. Kinematics measured during a biomechanical experiment were applied to the computer models of the bone surface. The overall accuracy of the method was shown to be at or below the accuracy of the digitization system used. For this experimental application, the accuracy was better than +/-0.1mm for position and 0.1 degrees for orientation for linkage digitization and better than +/-0.2mm and +/-0.2 degrees for CT digitization. Surface models of the radius and ulna were constructed from CT data, as an example application. Kinematics of the bones were measured for simulated forearm rotation. Screw-displacement axis analysis showed 0.1mm (proximal) translation of the radius (with respect to the ulna) from supination to neutral (85.2 degrees rotation) and 1.4mm (proximal) translation from neutral to pronation (65.3 degrees rotation). The motion of the radius with respect to the ulna was displayed using the surface models. This methodology is a useful tool for the measurement and application of rigid-body kinematics to computer models.     

Calibration of position and angular data from a magnetic tracking device

This paper describes a method for calibrating data from a magnetic tracking device. Position and orientation data were collected in a 1. 6x0.8x1.4m(3) volume using a Polhemus Fastrak((R)) in conjunction with both a long-range and standard transmitter. Position and orientation data were calibrated using a locally linear model based on the position of the measurement. After calibration, the average position and angular errors were less than 1.8cm and 1.2 degrees up to 1.8m from the transmitter for the long-range transmitter. For the standard transmitter, even after calibration, errors increased sharply when the sensor was more than 1.2m from the transmitter. Up to that distance, post-calibration errors were less than 1.2cm and 1. 2 degrees, while up to 1.8m they were below 5cm and 4 degrees. These errors could be further reduced by noise filtering. However, use of the standard transmitter is not recommended at distance greater than 1.2m due to orientation-based effects. It was concluded that for the volume investigated, tracking devices could provide similar three-dimensional accuracy to video systems.     

Linear fuzzy gene network models obtained from microarray data by exhaustive search

Background  

Recent technological advances in high-throughput data collection allow for experimental study of increasingly complex systems on the scale of the whole cellular genome and proteome. Gene network models are needed to interpret the resulting large and complex data sets. Rationally designed perturbations (e.g., gene knock-outs) can be used to iteratively refine hypothetical models, suggesting an approach for high-throughput biological system analysis. We introduce an approach to gene network modeling based on a scalable linear variant of fuzzy logic: a framework with greater resolution than Boolean logic models, but which, while still semi-quantitative, does not require the precise parameter measurement needed for chemical kinetics-based modeling.     

A portable system for collecting anatomical joint angles during stair ascent: a comparison with an optical tracking device

Background

Assessments of stair climbing in real-life situations using an optical tracking system are lacking, as it is difficult to adapt the system for use in and around full flights of stairs. Alternatively, a portable system that consists of inertial measurement units (IMUs) can be used to collect anatomical joint angles during stair ascent. The purpose of this study was to compare the anatomical joint angles obtained by IMUs to those calculated from position data of an optical tracking device.

Methods

Anatomical joint angles of the thigh, knee and ankle, obtained using IMUs and an optical tracking device, were compared for fourteen healthy subjects. Joint kinematics obtained with the two measurement devices were evaluated by calculating the root mean square error (RMSE) and by calculating a two-tailed Pearson product-moment correlation coefficient (r) between the two signals.

Results

Strong mean correlations (range 0.93 to 0.99) were found for the angles between the two measurement devices, as well as an average root mean square error (RMSE) of 4 degrees over all the joint angles, showing that the IMUs are a satisfactory system for measuring anatomical joint angles.

Conclusion

These highly portable body-worn inertial sensors can be used by clinicians and researchers alike, to accurately collect data during stair climbing in complex real-life situations.

     

Direct comparison of kinematic data collected using an electromagnetic tracking system versus a digital optical system

The purpose of this study was to quantify the dynamic accuracy of kinematics measured by a digital optical motion analysis system in a gait analysis laboratory (capture volume approximately 20m(3)) compared to a standard range direct-current electromagnetic (EM) tracking device (capture volume approximately 1m(3)). This is a subset of a larger effort to establish an appropriate marker set for the optical system to quantify upperlimb kinematics simultaneously with gait, in comparison to previous studies of isolated upperlimb movements that have employed EM tracking devices. Rigid clusters of spherical reflective markers and EM sensors were attached to a mechanical articulator that mimicked three-dimensional joint rotations, similar to the elbow. As the articulator was moved through known ranges of motion (i.e. gold standard), kinematic data were collected simultaneously using both tracking systems. Both systems were tended to underestimate the range of motion; however, the application of post hoc smoothing and least-squares correction algorithms reduced these effects. When smoothing and correction algorithms were used, the magnitude of the mean difference between the gold standard and either the EM or optical system did not exceed 2 degrees for any of the compound motions performed. This level of agreement suggests that the measurements obtained from either system are clinically comparable, provided appropriate smoothing and correction algorithms are employed.     

Factors influencing accuracy of screw displacement axis detection with a D.C.-based electromagnetic tracking system.

Recent technical improvements and cost reductions in electromagnetic motion tracking systems invite their application to motion axis determination in the surgical setting. After evaluation of the accuracy of a state-of-the-art D.C. electromagnetic tracking system, which generates complete three-dimensional kinematic outputs from just a single receiver, we calculated screw displacement axes (SDA's) from its source data. The accuracy of SDA determination from such source data was evaluated for various rotational increment sizes around a revolute joint. A novel smoothing procedure, customized for this type of source data, was developed, enabling SDA detection from incremental rotations of less than 1 deg, at an accuracy appropriate for intra-operative measurement of human joint motion. Examples of SDA determination are given for motion tracking of a ball joint and of the elbow articulation.     

Mammal tracking and nest location in Brazilian forest with an improved spool-and-line device

Spool-and-line mammal tracking devices, with a range of up to 2300 m, were constructed and used to trace Brazilian forest mammals. A retrieval rate of 64.6% was achieved with
263 animals of 16 species. Detailed observations of nests or refuges and behaviour were possible, which contributed to an understanding of the ecological segregation of the species. It was demonstrated that spool-and-line tracking was a valuable, neglected procedure, with some advantages over other methods, that could become an integral part of routine investigations of mammal ecology.     

Evaluation of Adair binding constants from experimental data obtained with an Imai cell. Statistical weighting functions

One group of laboratories uses an 'even-weighted', or unweighted, nonlinear least-squares method for the analysis of experimental oxygen binding data obtained with an Imai oxygenation apparatus. Another group uses an 'end-weighted' nonlinear least-squares analysis. With end weighting each observation is assigned a statistical weight which is proportional to Y(1-Y), where Y is the fractional saturation. In this work we discuss statistical weighting functions as applied to the Imai oxygenation apparatus and then determined what are the best weighting factors for an actual series of published experiments. Based on these calculations, it is concluded that the 'best' weighting for the Imai oxygenation apparatus is a very small amount of end weighting. Furthermore, the amount of end weighting is so small that even weighting, or no weighting at all, is also appropriate.     

Improvement of design of a surgical interface using an eye tracking device

Background

Surgical interfaces are used for helping surgeons in interpretation and quantification of the patient information, and for the presentation of an integrated workflow where all available data are combined to enable optimal treatments. Human factors research provides a systematic approach to design user interfaces with safety, accuracy, satisfaction and comfort. One of the human factors research called user-centered design approach is used to develop a surgical interface for kidney tumor cryoablation. An eye tracking device is used to obtain the best configuration of the developed surgical interface.

Methods

Surgical interface for kidney tumor cryoablation has been developed considering the four phases of user-centered design approach, which are analysis, design, implementation and deployment. Possible configurations of the surgical interface, which comprise various combinations of menu-based command controls, visual display of multi-modal medical images, 2D and 3D models of the surgical environment, graphical or tabulated information, visual alerts, etc., has been developed. Experiments of a simulated cryoablation of a tumor task have been performed with surgeons to evaluate the proposed surgical interface. Fixation durations and number of fixations at informative regions of the surgical interface have been analyzed, and these data are used to modify the surgical interface.

Results

Eye movement data has shown that participants concentrated their attention on informative regions more when the number of displayed Computer Tomography (CT) images has been reduced. Additionally, the time required to complete the kidney tumor cryoablation task by the participants had been decreased with the reduced number of CT images. Furthermore, the fixation durations obtained after the revision of the surgical interface are very close to what is observed in visual search and natural scene perception studies suggesting more efficient and comfortable interaction with the surgical interface. The National Aeronautics and Space Administration Task Load Index (NASA-TLX) and Short Post-Assessment Situational Awareness (SPASA) questionnaire results have shown that overall mental workload of surgeons related with surgical interface has been low as it has been aimed, and overall situational awareness scores of surgeons have been considerably high.

Conclusions

This preliminary study highlights the improvement of a developed surgical interface using eye tracking technology to obtain the best SI configuration. The results presented here reveal that visual surgical interface design prepared according to eye movement characteristics may lead to improved usability.

     

Bootstrap estimates of the statistical accuracy of thresholds obtained from psychometric functions

Estimates of the accuracy of a threshold obtained from a psychometric function are often based on asymptotic theory. When the number of trails is small, however, these estimates may be untrustworthy. A computer program is described that uses a move reliable bootstrap approach to obtaining estimates of the standard deviation and confidence limits of a threshold and of the slope and spread of the psychometric function for any criterion level of performance.     

Characterizing RNA ensembles from NMR data with kinematic models

Functional mechanisms of biomolecules often manifest themselves precisely in transient conformational substates. Researchers have long sought to structurally characterize dynamic processes in non-coding RNA, combining experimental data with computer algorithms. However, adequate exploration of conformational space for these highly dynamic molecules, starting from static crystal structures, remains challenging. Here, we report a new conformational sampling procedure, KGSrna, which can efficiently probe the native ensemble of RNA molecules in solution. We found that KGSrna ensembles accurately represent the conformational landscapes of 3D RNA encoded by NMR proton chemical shifts. KGSrna resolves motionally averaged NMR data into structural contributions; when coupled with residual dipolar coupling data, a KGSrna ensemble revealed a previously uncharacterized transient excited state of the HIV-1 trans-activation response element stem–loop. Ensemble-based interpretations of averaged data can aid in formulating and testing dynamic, motion-based hypotheses of functional mechanisms in RNAs with broad implications for RNA engineering and therapeutic intervention.     

Accurate state estimation from uncertain data and models: an application of data assimilation to mathematical models of human brain tumors

Background

Data assimilation refers to methods for updating the state vector (initial condition) of a complex spatiotemporal model (such as a numerical weather model) by combining new observations with one or more prior forecasts. We consider the potential feasibility of this approach for making short-term (60-day) forecasts of the growth and spread of a malignant brain cancer (glioblastoma multiforme) in individual patient cases, where the observations are synthetic magnetic resonance images of a hypothetical tumor.

Results

We apply a modern state estimation algorithm (the Local Ensemble Transform Kalman Filter), previously developed for numerical weather prediction, to two different mathematical models of glioblastoma, taking into account likely errors in model parameters and measurement uncertainties in magnetic resonance imaging. The filter can accurately shadow the growth of a representative synthetic tumor for 360 days (six 60-day forecast/update cycles) in the presence of a moderate degree of systematic model error and measurement noise.

Conclusions

The mathematical methodology described here may prove useful for other modeling efforts in biology and oncology. An accurate forecast system for glioblastoma may prove useful in clinical settings for treatment planning and patient counseling.

Reviewers

This article was reviewed by Anthony Almudevar, Tomas Radivoyevitch, and Kristin Swanson (nominated by Georg Luebeck).     

The accuracy of measuring glenohumeral motion with a surface humeral cuff

Conclusions about normal and pathologic shoulder motion are frequently made from studies using skin surface markers, yet accuracy of such sensors representing humeral motion is not well known. Nineteen subjects were investigated with flock of birds electromagnetic sensors attached to transcortical pins placed into the scapula and humerus, and a thermoplastic cuff secured on the arm. Subjects completed two repetitions of raising and lowering the arm in the sagittal, scapular and coronal planes, as well as shoulder internal and external rotation with the elbow at the side and abducted to 90°. Humeral motion was recorded simultaneously from surface and bone fixed sensors. The average magnitude of error was calculated for the surface and bone fixed measurements throughout the range of motion. ANOVA tested for differences across angles of elevation, raising and lowering, and differences in body mass index. For all five motions tested, the plane of elevation rotation average absolute error ranged from 0-2°, while the humeral elevation rotation average error ranged from 0-4°. The axial rotation average absolute error was much greater, ranging from 5° during elevation motions to approaching 30° at maximum excursion of internal/external rotation motions. Average absolute error was greater in subjects with body mass index greater than 25. Surface sensors are an accurate way of measuring humeral elevation rotations and plane of elevation rotations. Conversely, there is a large amount of average error for axial rotations when using a humeral cuff to measure glenohumeral internal/external rotation as the primary motion.     

The accuracy of ribosomal RNA comparative structure models

The determination of the 16S and 23S rRNA secondary structure models was initiated shortly after the first complete 16S and 23S rRNA sequences were determined in the late 1970s. The structures that are common to all 16S rRNAs and all 23S rRNAs were determined using comparative methods from the analysis of thousands of rRNA sequences. Twenty-plus years later, the 16S and 23S rRNA comparative structure models have been evaluated against the recently determined high-resolution crystal structures of the 30S and 50S ribosomal subunits. Nearly all of the predicted covariation-based base pairs, including the regular base pairs and helices, and the irregular base pairs and tertiary interactions, were present in the 30S and 50S crystal structures.     

Determining dual Euler angles of the ankle complex in vivo using "flock of birds" electromagnetic tracking device

The dual Euler angles method has been proposed as an alternative approach to describe the general spatial human joint motion. In this study, the dual Euler angles method was applied to study the three-dimensional motion of the ankle complex. The methodology for obtaining dual Euler angles of the ankle complex was developed by using a "Flock of Birds" electromagnetic tracking device. The repeatability of the methodology was studied based on the intertester and intratester variability analysis. Finally kinematic coupling characteristics of the ankle complex during dorsiflexion-plantarflexion, eversion-inversion, and abduction-adduction were analyzed according to the parameters of the dual Euler angles.     

Calibration of the "Flock of Birds" electromagnetic tracking device and its application in shoulder motion studies.

In this paper the applicability in terms of measurement accuracy of the "Flock of Birds" six D.O.F. electromagnetic tracking device in shoulder research is investigated. Position measurements in a workspace of approximately 1 m3 were performed using a stylus. The andom error at the stylus tip appeared to be 1.86, 1.98 and 2.54 mm for x-, y- and z-coordinate, respectively. The error caused by distortion of the magnetic field by metal in the concrete of especially the floor was 20.8, 22.2 and 20.4 mm for the x-, y- and z-coordinate, respectively. Calibration and leaving out the measurements closest to the floor lowered this error to 2.07, 2.38 and 2.35 mm. Orientation errors of the shoulder bones evolving from the measurement inaccuracy were estimated from repeated measurements of shoulder bony landmarks of ten subjects by means of the stylus. These errors were generally below 2 degrees. This is lower than found for the same measurements using a spatial linkage digitizer. It is concluded that the "Flock of Birds" is a useful tool for shoulder kinematic studies.     

Correlation of compartment pressure data from an intraoperative sensing device with postoperative fluoroscopic kinematic results in TKA patients

Fluoroscopy has recently been used to analyze postoperative kinematics in total knee arthroplasty (TKA). These analyses have reported varying results even in patients with similar implant design. In addition, patterns of wear in retrieved tibial polyethylene inserts of similar design have been found to vary substantially. These findings suggest that surgical technique, especially soft tissue balancing, may play a role in postoperative kinematics and implant failure. Accurate soft-tissue balancing is hypothesized to result in similar pressures within the medial and lateral compartments of the knee. However, a method of easily measuring these pressures at TKA has not been developed. In the present study, 32 patients were implanted with a mobile-bearing LCS TKA utilizing the balanced gap technique. An electronic pressure sensor, developed specifically to record pressure magnitude and distribution in the medial and lateral compartments, was incorporated into the implant trials. The knee was then passively taken through a range of motion while pressure data was recorded via computer. Postoperatively, 16 patients underwent active fluoroscopic kinematic analysis to assess for condylar liftoff and femorotibial translation. We found that abnormal compartment pressures and distributions as recorded by the intraoperative pressure sensor were correlated with inappropriate or paradoxical postoperative kinematics. In addition, subjects having similar pressures in both compartments throughout a range of motion did not experience condylar liftoff values greater than 1.0 mm. These data suggest that surgical technique influences the magnitude and distribution of forces at the articulation, postoperative kinematics, and likely, implant longevity.