Iterative coherent mapping of cardiac electrophysiological (EP) activation including reentry effects
A method includes receiving an input mesh representation of a cardiac chamber, a set of measured locations on a wall tissue of the cardiac chamber, and a respective set of local activation times (LATs) measured at the locations. The input mesh is re-meshed into a regular mesh including regularized polygons. The set of measured locations and respective LATs is data fitted to the regularized polygons. Respective LAT values are iteratively calculated for the regularized polygons, so as to obtain a cyclic EP activation wave solution over the regular mesh that take account of reentry of an EP wave. An electroanatomical map including the cyclic EP activation wave overlaid on the regular mesh is presented.
1. A method, comprising:
receiving an input mesh representation of a cardiac chamber, a set of measured locations on a wall tissue of the cardiac chamber, and a respective set of local activation times (LATs) measured at the locations;
re-meshing the input mesh into a regular mesh comprising regularized polygons;
data fitting the set of measured locations and respective LATs to the regularized polygons;
iteratively calculating for the regularized polygons respective LAT values, so as to obtain a cyclic EP activation wave solution over the regular mesh that take account of reentry of an EP wave, wherein iteratively calculating the LAT values comprises iteratively solving a set of complex valued equations that are configured to describe the reentry of an EP wave by including a reentry cycle length value; and
presenting an electroanatomical map comprising the cyclic EP activation wave overlaid on the regular mesh.
2. The method according to claim 1 , wherein iteratively calculating the LAT values comprises iteratively solving a set of three linear equations for the LAT values, slowness values and propagation of the EP wave.
3. The method according to claim 1 , wherein iteratively calculating the LAT values comprises reassigning to each regularized polygon a slowness vector recalculated using a proximity weight multiplier.
4. A system, comprising:
an interface, which is configured to receive an input mesh representation of a cardiac chamber, a set of measured locations on a wall tissue of the cardiac chamber and a respective set of local activation times (LAT) measured at the locations; and
a processor, which is configured to:
re-mesh the input mesh into a regular mesh comprising regularized polygons; data fit the set of measured locations and respective LATs to the regularized polygons;
iteratively calculate for the regularized polygons respective LAT values, so as to obtain a cyclic EP activation wave solution over the regular mesh that take account of reentry of an EP wave, wherein the processor is configured to iteratively calculate the LAT values by iteratively solving a set of complex valued equations that are configured to describe the reentry of an EP wave by including a reentry cycle length value; and
present an electroanatomical map comprising the cyclic EP activation wave overlaid on the regular mesh.
5. The system according to claim 4 , wherein the processor is configured to iteratively calculate the LAT values by iteratively solving a set of three linear equations for the LAT values, slowness values and propagation of the EP wave.
6. The system according to claim 5 , wherein the processor is configured to iteratively calculate the LAT values by reassigning to each regularized polygon a slowness vector recalculated using a proximity weight multiplier.