Postprocessing of 3-D current density reconstruction results with equivalent ellipsoids

A method of postprocessing and visualizing three-dimensional vector fields, such as current density reconstruction results, is presented. This method is based on equivalent ellipsoids fitted to the vector fields. The technique has been tested with simulated data and current density reconstructions based on bioelectromagnetic data obtained from a physical thorax phantom. Three different approaches based on: (1) longest distance; (2) dominant direction; and (3) principal component analysis, for fitting the equivalent ellipsoids are proposed. Multiple foci in vector fields are extracted by multiple ellipsoids which are fitted iteratively. The method enables statistical postprocessing for the sake of comparisons of different source reconstructions algorithms or comparisons of groups of patients or volunteers.