The accurate calculation of the skeleton of an object is a problem not satisfactorily solved by existing approaches. Most algorithms require a significant amount of user interaction and use a voxel grid to compute discrete and often coarse approximations of this representation of the data. We present a novel, automatic algorithm for computing subvoxel precise skeletons of volumetric data based on subvoxel precise distance fields. Most voxel based centerline and skeleton algorithms start with a binary mask and end with a list of voxels that define the centerline or skeleton. Even though subsequent smoothing may be applied, the results are inherently discrete. Our skeletonization algorithm uses as input a subvoxel precise distance field and employs a number of fast marching method propagations to extract the skeleton at subvoxel precision. We present the skeletons of various three-dimensional (3D) data sets and digital phantom models as validations of our algorithm.

The influence of head tissue conductivity on magnetoencephalography (MEG) was investigated by comparing the normal component of the magnetic field calculated at 61 detectors and the localization accuracy of realistic head finite element method (FEM) models using dipolar sources and containing altered scalp, skull, cerebrospinal fluid, gray, and white matter conductivities to the results obtained using a FEM realistic head model with the same dipolar sources but containing published baseline conductivity values. In the models containing altered conductivity values, the tissue conductivity values were varied, one at a time, between 10% and 200% of their baseline values, and then varied simultaneously. Although changes in conductivity values for a single tissue layer often altered the calculated magnetic field and source localization accuracy only slightly, varying multiple conductivity layers simultaneously caused significant discrepancies in calculated results. The conductivity of scalp, and to a lesser extent that of white and gray matter, appears especially influential in determining the magnetic field. Comparing the results obtained from models containing the baseline conductivity values to the results obtained using other published conductivity values suggests that inaccuracies can occur depending upon which tissue conductivity values are employed. We show the importance of accurate head tissue conductivities for MEG source localization in human brain, especially for deep dipole sources or when an accuracy greater than 1.4 cm is needed.

Most image processing and visualization applications allow users to configure computation parameters and manipulate the resulting visualizations. SCIRun, VolView, MeVisLab, and the Medical Interaction Toolkit (MITK) are four image processing and visualization frameworks that were built for these purposes. All frameworks are freely available and all allow the use of the ITK C++ library. In this paper, the benefits and limitations of each visualization framework are presented to aid both application developers and users in the decision of which framework may be best to use for their application. The analysis is based on more than 50 evaluation criteria, functionalities, and example applications. We report implementation times for various steps in the creation of a reference application in each of the compared frameworks. The data-flow programming frameworks, SCIRun and MeVisLab, were determined to be best for developing application prototypes, while VolView was advantageous for nonautomatic end-user applications based on existing ITK functionalities, and MITK was preferable for automated end-user applications that might include new ITK classes specifically designed for the application.

Volume currents are important for the accurate calculation of magnetoencephalographic (MEG) forward or inverse simulations in realistic head models. We verify the accuracy of our finite element method implementation for MEG simulations by comparing its results for spheres containing dipoles to those obtained from the analytic solution. We then use this finite element method to show that, in an inhomogeneous, nonspherical realistic head model, the magnetic field normal to the MEG detector due to volume currents often has a magnitude on the same order or greater than the magnitude of the normal component of the primary magnetic field from the dipole. We also demonstrate the disparity in forward solutions between a model that employs spheres, one that uses the realistic head and primary currents alone, and a realistic head model that incorporates both primary and volume currents. In forward and inverse MEG simulations using the inhomogeneous realistic model, the results obtained from calculations containing volume currents are more accurate than those derived without considering volume currents.