Translator Disclaimer
1 July 2007 Three-dimensional panoramic imaging of cardiac arrhythmias in rabbit heart
Author Affiliations +
Cardiac fluorescent optical imaging provides the unique opportunity to investigate the dynamics of propagating electrical waves during ventricular arrhythmias and the termination of arrhythmias by strong electric shocks. Panoramic imaging systems using charge-coupled device (CCD) cameras as the photodetector have been developed to overcome the inability to monitor electrical activity from the entire cardiac surface. Photodiode arrays (PDAs) are known to have higher temporal resolution and signal quality, but lower spatial resolution compared to CCD cameras. We construct a panoramic imaging system with three PDAs and image Langendorff perfused rabbit hearts (n=18) during normal sinus rhythm, epicardial pacing, and arrhythmias. The recorded spatiotemporal dynamics of electrical activity is texture mapped onto a reconstructed 3-D geometrical heart model specific to each heart studied. The PDA-based system provides sufficient spatial resolution (1.72 mm without interpolation) for the study of wavefront propagation in the rabbit heart. The reconstructed 3-D electrical activity provides us with a powerful tool to investigate the fundamental mechanisms of arrhythmia maintenance and termination.



Characterization of the spread of electrical activity is essential for understanding the mechanisms responsible for normal cardiac rhythm, arrhythmias, and antiarrhythmia therapies. Cardiac optical mapping, in which myocardial electrical activity is simultaneously recorded from hundreds or thousands of sites, has made great strides in furthering our understanding of the initiation, maintenance, and termination of arrhythmias. In optical mapping of transmembrane potential, heart tissue is stained with a voltage-sensitive dye and illuminated with an excitation light source.1, 2 The resulting emission fluorescence is proportional to the transmembrane potential. In contrast to electrode mapping techniques, optical mapping has the ability to faithfully reproduce transmembrane action potential morphology while being optically isolated from the overwhelming electric field applied during defibrillation shocks.1 Therefore, the optical mapping technique is a powerful tool for elucidating the exact physiological mechanisms of cardiac arrhythmias and defibrillation.

The monocular principle is predominantly used for cardiac optical mapping. The mapped region is limited to the field of view of the optical sensor. In cardiac arrhythmias, single or multiple coexisting reentrant wavefronts have been observed in numerous studies.3, 4, 5 The core of the reentrant arrhythmia can be highly unstable and meander across the epicardium.6 Thus, a study using a monocular imaging system cannot collect the full information during arrhythmia if the core of reentry leaves the field of view or if a reentry core is beyond the field of view. This limitation strongly motivated efforts to build panoramic imaging systems that reveal the electrical activity on the entire ventricular epicardium. The first implementations of the panoramic technique7, 8, 9, 10 introduced a cost-efficient method in which the investigators optically mapped electrical activity using a charge-coupled device (CCD) camera and a panoramic mirror arrangement to obtain the full ventricular epicardial view. After image registration, the electrical data was texture mapped onto the reconstructed heart geometry. Kay, Amison, and Rogers11 extended this idea to a panoramic optical mapping system capable of imaging large hearts, which used two CCD cameras and one mirror to obtain four views of the heart.

Optical imaging of the intact heart is usually performed with CCD cameras and photodiode array (PDA) detectors. CCD technology could potentially offer a significant advantage of higher spatial resolution due to the large number of pixels on a CCD sensor. However, the rate of data acquisition is usually substantially lower than that achieved with a PDA system. The rate can be increased by pixel binning, yet this defeats the major advantage of CCD technology, since binning effectively reduces the spatial resolution. This limitation is particularly acute when studying the shock response of defibrillation, since the duration of defibrillation waveforms is only a few milliseconds. Moreover, the transmembrane potential response to strong electric shocks can be very fast (less than 1ms )12 and contains much higher frequency components compared to a normal propagated response, reinforcing the need for fast sampling rates.

In addition to higher temporal resolution, PDAs also provide higher quality signals as compared to that of a CCD camera. One of the problems during studies of the biophysical mechanisms of defibrillation is the fact that the imaging sensor can be partially obstructed by the shock electrodes in Langendorff perfused heart experiments that mimic external defibrillation. As a result, the signal-to-noise ratio (SNR) of recordings obtained from the obstructed area are lower than that of recordings from unobstructed areas. Furthermore, when imaging diseased hearts such as a heart with a healed myocardial infarction, the optical signals from the unhealthy regions are also lower as compared to those from normal tissue. In addition to both of these factors, optical signals recorded during an arrhythmia can have extremely low amplitudes near the reentry core, regardless of if the tissue is healthy or diseased. All of these difficulties reinforce the requirement of high signal quality, which may not be achieved with a CCD-based system.

Defibrillation has been studied extensively in various in-vivo and in-vitro heart models. However, many findings have been limited, since these studies used functionally and structurally normal heart models, whereas a large percentage of patients who receive defibrillation therapy are actually suffering from coronary diseases such as ischemia and myocardial infarction. Thus far, vulnerability and defibrillation have not been widely studied by optical mapping at the whole heart level under these disease conditions. In this study, we developed a PDA-based 3-D fast fluorescence panoramic imaging (FFPI) system with high temporal resolution and signal quality, making this system well suited to study the mechanisms of defibrillation in the diseased heart. This system operates at a 5000framessec sampling rate and has 768pixels in total, providing us with the unique opportunity to visualize the electrical activity on the epicardial surface of the rabbit heart.


Materials and Methods


Isolated Rabbit Heart

The protocol was approved by the Institutional Animal Care and Use Committee at Washington University. The hearts of New Zealand white rabbits (n=18) were imaged in this study. These animals consisted of healthy control rabbits as well as those with diseased hearts, including healed myocardial infarction and hypertrophic cardiomyopathy.13 The rabbits were injected intravenously with sodium pentobarbital (50mgkg) and with 2000 U heparin. The hearts were quickly removed, placed on a Langendorff apparatus, and perfused with oxygenated modified Tyrode’s solution as previously described.12

The hearts were stained by a gradual injection of 50μL of stock solution (1.25mgmL) of the voltage-sensitive dye di-4-ANEPPS (Molecular Probes, Eugene, Oregon) diluted in dimethylsulfoxide (DMSO; Fisher Scientific, Fair Lawn, New Jersey), delivered by a micropump over 5min . The excitation-contraction uncoupler 2,3-butanedione monoxime (BDM, 15mM ; Fisher Scientific, Fair Lawn, New Jersey) was added to the perfusate to suppress motion artifacts in the optical recordings.


Fast Fluorescence Panoramic Imaging System

As shown in Figs. 1a and 1b , the hearts were positioned in a hexagonal perfusion chamber filled with Tyrode’s solution. Three photodiode arrays (PDAs model C4675-103, Hamamatsu, Bridgewater, New Jersey) were spaced 120deg apart and directed toward the center of the chamber. The other three faces of the chamber were used for illumination by three commercially available light emitting diode (LED) arrays (Luxeon Flood 18-up, Quadica Developments, Calgary, Canada). To mimic the electrode configuration of external defibrillation, two stainless-steel mesh electrodes were placed into the solution chamber in an orientation perpendicular to the projection of PDA-1. The hearts were oriented so that the right ventricle faced PDA-1 and the left ventricle faced the mesh electrode distant from PDA-1. The perfusion cannula was connected to a rotation stage. A digital camera (Sony DSC-S70) was used to take images of the heart ( 640×480pixel resolution) at 10-deg increments as the heart was rotated a full 360deg . These images were used to reconstruct the 3-D heart geometry. As shown in Fig. 1c, for each individual PDA, the fluorescence emitted from the heart was filtered using an emission filter (> 610nm) , and collected by a PDA with built-in first-stage preamplifiers. The outputs of each PDA were fed into a custom-made 256-channel second-stage amplifier (Innovative Technology, Brooksville, Florida) and then recorded by a data acquisition system (National Instruments, Austin, Texas) at 5000framessec with 16-bit resolution. Instrumentation channels recorded the shock field strength, electrocardiogram, shock voltage, pacing stimuli, and defibrillation triggers, which were saved for off-line data analysis.

Fig. 1

(A) Schematic diagram of PDA-based panoramic imaging system. (B) Panoramic imaging system setup. (C) Schematic diagram of optical mapping components.



Camera Model

A commonly accepted camera calibration model, the affine distortion model, states that a 3-D point M(M=[X,Y,Z]T) is projected with a perspective projection onto an image plane on a 2-D point m(m=[μ,v,1]T) based on the projection equation:

Eq. 1

where fμ and fv are the focal distances expressed in units of horizontal and vertical pixels. The principal coordinates μ0 and v0 correspond to the image of the optical center. The skew factor γ encodes the angle between the x and y pixel axes, allowing the camera model to handle nonsquare pixels. The digital camera we used to document the heart geometry has square pixels so that the skew factor γ was known to be 0. A more detailed description of this camera model can be found at the following website:ḇdoc/htmls/parameters.html We used this model for heart surface reconstruction and texture mapping, as described next.


Heart Surface Reconstruction

The left panel of Fig. 2a shows the flowchart of the algorithm for reconstructing the heart surface. The occluding contours algorithm14 has been used in a previous panoramic imaging study9 as the principle method to reconstruct the heart surface geometry. The essence of the occluding contours method is to iteratively shave a virtual 3-D cube by silhouette edges to obtain the volume of an object inside the cube. Kay, Amison, and Rogers11 incorporated an adaptive octree mesh refinement algorithm into the occluding contours method to reduce computational load and memory requirements. We implemented these algorithms as follows.

  • 1. The digital camera was positioned at a fixed distance from the solution chamber, while its optical axis was aligned to intersect with the axis of rotation of the rotation stage. After we solved the intrinsic parameters of the camera model, we took 36 images (up to 0.12-mmpixel resolution) of the heart while the rotation stage rotated around its axis for a full 360-deg revolution with a 10-deg rotation step.

  • 2. The heart boundary in these images was extracted by a combined image processing procedure, including intensity adjustment, intensity thresholding, image opening, and image closing. After heart boundary detection, we created silhouettes (36 in total) for these images by setting the pixels on the heart to a silhouette value (α) of 1, and the rest of pixels to a silhouette value of 0.

  • 3. A virtual 3-D cube (40×40×40mm3) just large enough to contain the heart was created. The cube was initially divided into eight voxels.

  • 4. The voxel vertices were projected to the camera imaging plane based on Eq. 1 to determine their silhouette value (α) using bilinear interpolation of the silhouette image at the corresponding rotation angle. The value of α of a single vertex was clamped to 0 for the remainder of the analysis whenever α was found to be equal to 0, thereby indicating that the vertex was outside the heart. We then rotated the cube 10deg around the axis of rotation and computed α of the vertices from the next corresponding silhouette image. This procedure was repeated for all the silhouette images. Voxels that had all eight corners outside the heart volume (α<0.95) and voxels that had all eight corners inside the heart volume (α0.95) were excluded from further analysis.

  • 5. For the remaining voxels, each of them was further evenly divided into another eight voxels, and step 4 was repeated. The entire process was repeated until the desired resolution was achieved (0.2mm) . The centroid of these voxels formed a set of scattered points approaching the heart surface.

Fig. 2

(a) Flowchart of the algorithms for reconstructing the heart surface and texture mapping the epicardial action potentials. (b) Visible angle calculation.


We reconstructed the heart surface by defining an analytic function E that maps the abstract, topological sphere S2 into R3 . The function E was built by blending together multiple, overlapping polynomial functions, called a chart, each of which maps a subset of the sphere to R3 . The result is a rational polynomial embedding of S2 15 into R3 with guaranteed continuity Ck , i.e., the first k derivatives are defined and smooth. The data points found in the previous step were placed on the sphere, respecting their local connectivity. Each chart was then fit to a connected subset of these points to locally approximate the surface. Once the surface is reconstructed, it can be tessellated at any desired resolution with near-equilateral triangles. The charts also provide a local parameterization for every point on the surface, suitable for use in texture mapping and representing other data on the surface.

More specifically, we began by finding a local tangent neighborhood for each data point. We first estimated a surface normal16 if one was not given. The tangent neighborhood consisted of the four (or more) nearest data points, which together formed an enclosing ring around the given data point, when projected onto the tangent plane. We next grouped the points into connected subsets. A chart was created by choosing a seed point, then taking all of the data points within a geodesic distance r from the seed point, as measured in the tangent neighborhood graph. We used a greedy algorithm to choose the chart seed points by choosing an uncovered point that is close to the ideal distance from one or more existing chart centers. The goal was to place the chart centers at a distance 2r-gr from each other, where g0.3 is the desired chart overlap. This tends to produce a hexagonal tiling. Note that data points may (and will) appear in multiple charts. For this dataset, we set r to be 110 of the height of the heart.

Once the chart seeds and groupings were identified, they (and the data points) were mapped to an abstract representation of the sphere, preserving local neighborhood information for both.15 This guarantees that the resulting surface will have spherical topology. At this stage, any gaps due to missing data were filled in by adding additional charts. Each chart was then fit to its corresponding data and blended into the final function using a Ck blend function, which is 1 in the middle of the chart and decays to zero by the boundary. Each individual polynomial function approximates its data within a given epsilon (0.1 of the average distance between points); the approximation error of the entire surface E is less than, or equal to, the individual function’s error.


Texture Mapping of Optically Recorded Action Potential

The next step, shown in the right panel of Fig. 2a, was texture mapping the optically recorded action potential onto the reconstructed heart surface mesh. We developed a robust algorithm to assign such data to each element in the surface mesh.

  • 1. Register mesh to PDA projections. Shown in Fig. 2b, the reconstructed heart surface mesh was registered to PDA projections by calculating the visible angle (θ) between the outward normal vector (n) of each mesh cell and each normalized PDA projection vector (p) . Then we defined the projection angle (Φ) as:

    Eq. 2

    Mesh cells visible from a particular view have projection angles greater than 0deg for that view. The view with the maximum projection angle at 90deg for a mesh cell provides the best vantage point for viewing the surface of the mesh cell.

  • 2. Single- or dual-projection texture mapping. The texture mapping procedure depends on the registration of mesh cells. For each mesh cell registered to a single PDA projection, its centroid was back projected onto the corresponding PDA imaging plane to determine the fluorescence using bilinear interpolation. For each mesh cell registered to two PDA projections, the same procedure was performed twice but for different PDA projections. We then computed the weighted average fluorescence from the two raw fluorescence signals as the fluorescence of the mesh cell based on the following equation:

    Eq. 3

    where F1 , F2 , and Φ1 , Φ2 are the fluorescence signal and projection angle of the two registered PDAs.

  • 3. Fluorescent signals were scaled to mV, assuming that a normal resting potential of 85mV and action potential amplitude of 100mV were present at all of the mesh cells.

  • 4. We performed 3-D visualization in Matlab (The Mathworks, Incorporated, Natick, Massachusetts) using the Matlab multifaceted patches function (patch.m).


Experimental Protocol and Data Analysis

A bipolar Ag–AgCl pacing electrode with 1-mm interelectrode distance was placed at the anterior epicardium. We first recorded sinus rhythm, then the heart was paced at 300-ms basic cycle length by 2-ms stimuli and the electrical activity during this epicardial pacing was recorded. After the pacing stimuli, a test shock was delivered to the heart through the mesh electrodes spaced 100mm apart inside the solution chamber. Shocks were delivered using a custom-made defibrillator, which consisted of five capacitors ( 3100μF each) and a triggering circuit controlled by a TTL pulse from the computer. Arrhythmias were introduced by either burst pacing from the bipolar pacing electrode or a T-wave shock. Sustained arrhythmias were recorded, and then an extra shock was delivered to restore the normal rhythm.

The SNR presented in this study is peak-to-peak SNR. We selected a single beat of sinus rhythm, and the peak-to-peak amplitude of the noise was computed during the action potential phase-0 (diastole), and the peak-to-peak amplitude of the florescent action potential was computed during action potential systole.



White-light sources such as tungsten halogen lamps and mercury arc lamps have been widely used in early optical mapping systems in combination with narrow bandpass filters to select the desired excitation wavelength and spectral bandwidth. Light emitting diodes (LEDs) provide an attractive option for excitation light sources,17, 18 since LEDs are significantly less expensive than white light sources. In this study, we compared fluorescence recordings using two excitation light sources: a 250-W tungsten halogen lamp and Luxeon Flood LEDs. Figure 3 shows the amplitude changes of an optically recorded action potential from a single site [in Fig. 3a, cyan color] under different LED illumination configurations as well as with the filtered tungsten halogen lamp. Figure 3b shows the SNR (mean±std) within the dashed rectangle in Fig. 3a for each illumination configuration. These results show that the use of Luxeon LEDs as the excitation light source produced fluorescence signals with a higher SNR under most illumination configurations compared to the tungsten halogen lamp.

Fig. 3

Comparison of optically recorded transmembrane action potential signals using two different types of light sources ( 250-W tungsten-halogen lamp and LED matrix). (A) Optically recorded transmembrane action potentials superimposed onto the anterior epicardial heart surface. The solid rectangle represents the field of view of the PDA. The thickened cyan color trace shows the site from where signals were illustrated in (B). (B) Amplitude changes of optical action potential from a single site [shown in (A), cyan color] under different illumination configurations: 1. when the anterior epicardium was best illuminated by two LEDs, 2. when the heart was panoramically illuminated by three LEDs located around the heart with 120deg between any two of them [shown in Figs. 1a and 1b], 3. when the heart was illuminated by bandpass filtered (520±45nm) light from a 250-W tungsten-halogen lamp, and 4. when the PDA was looking through a shock mesh electrode inside the solution chamber ( 15mm away from the heart) while the heart was panoramically illuminated. Shown at the right is the signal-to-noise ratio (SNR, mean±std ) within the dashed rectangle area in (A).


We then compared fluorescence recordings recorded by a PDA and a 128×128pixel CCD camera (CA-D1-0128T-STDL, Dalsa, Waterloo, Ontario, Canada). For this analysis, the heart was illuminated by an LED matrix. As indicated in Fig. 4a , the signals within the white rectangle were used for SNR comparison. The PDA system collected data at a 5000framessec sample rate, whereas the CCD system used 466framessec . As evident from this figure, PDA imaging yielded not only high temporal resolution, but also a significantly higher SNR compared to this particular CCD camera. Even with 8×8 binning of the CCD data, the PDA provided a much larger SNR while providing comparable spatial resolution.

Fig. 4

Comparison of transmembrane action potentials recorded by a PDA and a CCD camera. (A) Heart was illuminated by an LED matrix. Signals within the white rectangle were used for SNR comparison shown in (B), right side. The cyan rectangle shows the site and equivalent pixel size of the signal illustrated in (B), top trace. The green rectangle shows the site and equivalent pixel size of the signal illustrated in (B), middle and bottom trace. (B) Optically recorded action potential from a single site. Top: 2×2 binning of a raw CCD recording; the blue line is the signal filtered by a simple FIR low-pass filter. Middle. 8×8 binning of the raw CCD recording, the dark green line is the signal filtered by the same low-pass filter. Bottom: signal recorded by PDA from the same location as the CCD 8×8 binning. Shown at the right is the SNR (mean±std) within the solid white rectangle in (A). The PDA system collected data at 5000-framessec sample rate; the CCD system used 466-framessec sample rate.


One potential advantage of using multiple cameras to visualize an object is that it may be possible to improve the SNR in the areas visible by multiple sensors. To confirm this, all the heart surface mesh cells were registered. Each single mesh cell can be registered in one of the following values: 1 to 3 (only visible to PDA1/2/3), 4 (visible to PDA 1 and 2), 5 (visible to PDA 2 and 3), 6 (visible to PDA 3 and 1), or 7 (not visible to any PDAs). Figure 5a shows the mesh cells registration, values 1 to 6 are represented by different gray colors, and value 7 is represented by black. For all the dual-registered mesh cells with projection angles Φ1 and Φ2 , we computed two types of SNRs: SNRO and SNRN . SNRO is the SNR of the fluorescence signal after a weighted averaging process [see Eq. 3], SNRN is the SNR of the fluorescence signal from a single PDA that has the larger projection angle. Figure 5b demonstrates that SNRO is larger than SNRN in a majority of the dual-registered mesh cells, which indicates that we can improve the SNR by averaging the two fluorescence signals recorded by different PDAs. Also, when the absolute difference between Φ1 and Φ2 (ΔΦ,ΔΦ=Φ1Φ2) increased, the fluorescence signal from the PDA with the larger projection angle dominates the SNR calculation, thus the difference between SNRO and SNRN becomes smaller [shown in Fig. 5b].

Fig. 5

(a) Mesh cells projection registration. (b) Improvement of SNR by averaging fluorescence signals from two PDAs (see text for details).


Figure 6a shows an example of the reconstructed rabbit heart surface and epicardial action potential texture mapping as visualized from the three PDA views. In the left column, the heart surface geometry is shown, and texture mapping of the electrical activity present during epicardial pacing and shock-induced ventricular tachycardia are shown in the middle and right columns, respectively. Individual optical signals from locations a through f are shown in Fig. 6b. The heart was first ventricularly paced (t=t1) , then a shock from two mesh electrodes was delivered at the plateau of the action potentials, which induced a sustained ventricular tachycardia (t=t2) .

Fig. 6

Reconstructed heart surface and epicardial action potential texture mapping. (A) Left: reconstructed heart surface visualized from projections of PDA arrays. Middle: epicardial action potential texture mapping during epicardial pacing ( p in PDA-2 projection is the pacing site). Right: epicardial action potential texture mapping during shock-induced ventricular tachycardia. Signals from locations a to f are shown in (B). (B) Action potentials from locations a to f [shown in (A), right column]. The heart was first ventricularly paced ( t=t1 is shown in (A), middle column), then a shock from two mesh electrodes was delivered at the plateau of the action potentials, which induced a sustained ventricular tachycardia [ t=t2 is shown in (A), right column].




In this study, we developed a novel PDA-based fast fluorescence panoramic imaging (FFPI) system operated at 5000framessec with 768pixels in total. The FFPI system provides high quality fluorescence signals from a majority of the epicardium of the Langendorff perfused rabbit heart.

Previously, two CCD-based panoramic imaging systems have been developed:9, 11 CCD technology has a significant advantage of higher spatial resolution due to the large number of pixels on a CCD sensor. Bray, Lin, and Wikswo10 have demonstrated sufficient spatial resolution of a CCD-based panoramic imaging system in the study of 3-D cardiac electrodynamic behavior. Another CCD-based system developed by Kay, Amison, and Rogers11 demonstrated sufficient spatial resolution ( 1.7mm average spatial resolution before image processing) for the study of ventricular fibrillation (VF) in large heart models. However, PDAs are more commonly used in studies where high temporal resolution and high dynamic range are needed, such as the study of defibrillation. Our results demonstrate that the PDA system reported in this study can achieve a higher SNR at an approximately 10 times faster sample rate compared to the CCD camera tested. However, other commercially available CCD cameras may provide higher SNRs and faster sample rates than the tested camera.

The spatial resolution of optical mapping is dependent on the surface area mapped and the number of available pixels. In our studies, each PDA (16×16pixels) imaged 760mm2 . This area contained a small portion of the atrial epicardium and most of the ventricular epicardium. Three PDAs provided 768pixels in total. Approximately 570 of those pixels contained data. On average, each pixel mapped an area of 2.9mm2 , providing an average spatial resolution of 1.72mm before the application of bilinear interpolation. This spatial resolution is very similar to the spatial resolution achieved in the previous CCD-based system11 in large heart models, and is high enough for the study of wavefront propagation during arrhythmias.19 However, because the spatial resolution of this panoramic imaging system is limited by the number of photodiodes in each PDA, this system cannot be directly used on large hearts. Improvements in complementary metal-oxide semiconductor (CMOS) technology have produced a family of novel image sensors with high speed image acquisition while retaining the quantum efficiency of CCD. CMOS cameras are more costly than both CCD and PDA cameras. However, due to these clear advantages, CMOS cameras should become competitive soon.

During defibrillation, the time constant of the membrane response to a shock depends on the shock strength and refractory stage of the tissue when the shock is delivered. At the early plateau of the action potential, the fastest time constant can be less than 1ms when the applied shock is strong enough to create electroporation.20, 21 Another fast membrane potential change occurs when a shock is applied during diastole. Using a 5000-framessec sampling rate, Sharifov and Fast22 observed fast activation at about 0.6to0.7ms when an intermediate strength shock was delivered during diastole. Therefore, the fastest frequency component could be as high as approximately 2000Hz . According to the Nyquist sampling theory, 4000framessec is needed to accurately record these signals. Although even higher sampling frequencies are desirable, decreased signal quality at higher sampling rates is a tradeoff. Therefore, we used 5000framessec , despite the fact that our FFPI system can be operated as high as 10,000framessec .

One of the major advantages of our FFPI system is its ability to record high SNR signals even when the PDAs are partially obstructed by the mesh electrodes used to deliver external defibrillation shocks, making this system well suited for the study of diseased hearts, which can have very low amplitude optical signals. According to our experiments, a clear attenuation effect was observed when the mesh electrode was positioned within 1cm of the nearest heart surface. The attenuation effect rapidly decreased as we moved the mesh electrode away from the heart. At distances larger than 2.5cm , we could not see any difference in the morphology of the recorded action potentials with and without the mesh electrode, except in the amplitude of the signal [shown in Fig. 3b.]

In our experiments, each rotation of the heart results in a slight swing. Thus, it takes several seconds to let the heart stabilize before image capture. It takes approximately 9to10min to complete the full 360-deg rotation procedure at a 10-deg step size. Initially, 5-deg steps were used. The difference between these two step sizes has a minimum effect on the geometric reconstruction, primarily because the curvature of the ventricles is very smooth. Therefore, we selected a 10-deg step size to expedite the procedure. However, a finer step size is probably needed to accurately reconstruct a more complex anatomical structure.

Many light sources have been used in optical mapping systems, including lasers,23 DC-powered tungsten-halogen lamps,24 and most recently, light-emitting diodes.17, 18, 25 All of these light sources have their own unique properties and limitations. In this system, we used commercially available LED arrays, the Luxeon Flood, as the excitation light source. The Luxeon Flood is constructed by 18 Luxeon emitters (green, typical wavelength=530nm , spectral half width=35nm ) mounted on a rectangular PCB to deliver the most light output in the smallest possible space. Compared with the traditional illumination method of a tungsten-halogen lamp accompanied with a bandpass green filter and dichroic mirror, the Luxeon Flood is much more cost effective, ranging from $100 to $200 per LED array, compared to several thousand dollars for the light source. Another advantage is that each emitter has a 110-deg viewing angle. Thus, the Luxeon Flood provides uniform illumination at a distance larger than 45mm . For the FFPI system, we connected three Luxeon Floods in parallel and powered them with a constant voltage power supply at 20V and 2.1A ( 0.7A per flood). Although we have achieved satisfactory signal quality [see Figs. 1a and 1b], it is possible to further improve the florescence signal quality by increasing the excitation light intensity [as indicated in Fig. 3b as blue and red signals]. This can be done by increasing the number of floods, increasing the number of emitters on each flood, or by increasing the driving current (up to 1.05A ).


Study Limitations

In this study, we did not directly address the volume change in the Langendorff preparation. An advantage of CCD-based systems is that any change in volume of the heart over the course of the experiment can be examined directly through the CCD camera during the geometric reconstruction phase, as well as throughout the optical data collection phase. Using this method, Bray, Lin, and Wikswo did not see a significant change in volume throughout the course of their experiments.9 Kay, Amison, and Rogers demonstrated that the heart volume increases rapidly within the first 40min in Langendorff perfused swine hearts after exposure to DAM.11 In our FFPI system, we cannot directly examine volume changes from the PDA. Therefore, to minimize the overall effects of volume changes on the geometric reconstruction and texture mapping procedures, the heart was given a longer time (at least 30min ) to stabilize on the Langendorff apparatus before heart rotation and image acquisition began. In addition to contributing to volume changes of the Langendorff perfused heart, BDM also has an effect on a variety of ion channels and may alter the action potential duration in a number of species.26, 27 Therefore, the effects of BDM need to be taken into consideration for an appropriately designed experiment. However, we have recently identified a new excitation-contraction uncoupler, blebbistatin, which may solve this problem.28

In this study, we allowed at least 15sec for the LED light sources to reach a steady state before data acquisition. However, we did not directly measure the time to steady-state spectrum, intensity, or noise of the LED. These characteristics need to be examined in a future study.

Finally, cardiac electrical activity is essentially a 3-D phenomenon, in particular during complex arrhythmias. Results of this study are limited due to the typical epicaridal penetration depth of the optical mapping technique.


This work was supported by National Institutes of Health (NIH) grants (HL67322, HL074283) and National Science Foundation grant 049856.



I. R. Efimov, V. P. Nikolski, and G. Salama, “Optical imaging of the heart,” Circ. Res., 95 21 –33 (2004). 0009-7330 Google Scholar


M. Morad and G. Salama, “Optical probes of membrane potential in heart muscle,” J. Physiol. (London), 292 267 –295 (1979). 0022-3751 Google Scholar


F. X. Witkowski, L. J. Leon, P. A. Penkoske, W. R. Giles, M. L. Spano, W. L. Ditto, and A. T. Winfree, “Spatiotemporal evolution of ventricular fibrillation [see comments],” Nature (London), 392 78 –82 (1998). 0028-0836 Google Scholar


A. Garfinkel, Y. H. Kim, O. Voroshilovsky, Z. Qu, J. R. Kil, M. H. Lee, H. S. Karagueuzian, J. N. Weiss, and P. S. Chen, “Preventing ventricular fibrillation by flattening cardiac restitution,” Proc. Natl. Acad. Sci. U.S.A., 97 6061 –6066 (2000). 0027-8424 Google Scholar


A. V. Zaitsev, P. K. Guha, F. Sarmast, A. Kolli, O. Berenfeld, A. M. Pertsov, J. R. de Groot, R. Coronel, and J. Jalife, “Wavebreak formation during ventricular fibrillation in the isolated, regionally ischemic pig heart,” Circ. Res., 92 546 –553 (2003). 0009-7330 Google Scholar


J. M. Davidenko, A. V. Pertsov, R. Salomonsz, W. Baxter, and J. Jalife, “Stationary and drifting spiral waves of excitation in isolated cardiac muscle,” Nature (London), 355 349 –351 (1992). 0028-0836 Google Scholar


M. A. Bray, S. F. Lin, and J. P. Wikswo, “Panoramic epifluorescent visualization of cardiac action potential activity,” Proc. SPIE, 3658 99 –107 (1999). 0277-786X Google Scholar


S. F. Lin and J. P. Wikswo, “Panoramic optical imaging of electrical propagation in isolated heart,” J. Biomed. Opt., 4 (2), 200 –207 (1999). 1083-3668 Google Scholar


M. A. Bray, S. F. Lin, J. P. Wikswo Jr., “Three-dimensional surface reconstruction and fluorescent visualization of cardiac activation,” IEEE Trans. Biomed. Eng., 47 1382 –1391 (2000). 0018-9294 Google Scholar


M. A. Bray, S. F. Lin, and J. Wikswo, “Three-dimensional visualization of phase singularities on the isolated rabbit heart,” J. Cardiovasc. Electrophysiol., 13 1311 (2002). 1045-3873 Google Scholar


M. W. Kay, P. M. Amison, and J. M. Rogers, “Three-dimensional surface reconstruction and panoramic optical mapping of large hearts,” IEEE Trans. Biomed. Eng., 51 1219 –1229 (2004). 0018-9294 Google Scholar


I. R. Efimov, Y. N. Cheng, M. Biermann, D. R. Van Wagoner, T. Mazgalev, and P. J. Tchou, “Transmembrane voltage changes produced by real and virtual electrodes during monophasic defibrillation shock delivered by an implantable electrode,” J. Cardiovasc. Electrophysiol., 8 1031 –1045 (1997). 1045-3873 Google Scholar


A. J. Marian, Y. Wu, D. S. Lim, M. McCluggage, K. Youker, Q. T. Yu, R. Brugada, F. DeMayo, M. Quinones, and R. Roberts, “A transgenic rabbit model for human hypertrophic cardiomyopathy,” J. Clin. Invest., 104 1683 –1692 (1999). 0021-9738 Google Scholar


W. Niem, “Robust and fast modeling of 3D natural objects from multiple views,” Proc. SPIE, 2182 388 –397 (1994). 0277-786X Google Scholar


C. Grimm, “Spherical manifolds for adaptive resolution surface modeling,” Graphite, (2005) Google Scholar


N. J. Mitra, A. Nguyen, and L. Guibas, “Estimating surface normals in noisy point cloud data,” Int. J. Comput. Geom. Appl., 14 261 –276 (2004). 0218-1959 Google Scholar


I. Kodama, I. Sakuma, N. Shibata, S. B. Knisley, R. Niwa, and H. Honjo, “Regional differences in arrhythmogenic aftereffects of high intensity DC stimulation in the ventricles,” Pacing Clin. Electrophysiol., 23 807 –817 (2000). 0147-8389 Google Scholar


E. Entcheva, Y. Kostov, E. Tchernev, and L. Tung, “Fluorescence imaging of electrical activity in cardiac cells using an all-solid-state system,” IEEE Trans. Biomed. Eng., 51 333 –341 (2004). 0018-9294 Google Scholar


P. V. Bayly, E. E. Johnson, S. F. Idriss, R. E. Ideker, and W. M. Smith, “Efficient electrode spacing for examining spatial organization during ventricular fibrillation,” IEEE Trans. Biomed. Eng., 40 1060 –1066 (1993). 0018-9294 Google Scholar


D. K. Cheng, L. Tung, and E. A. Sobie, “Nonuniform responses of transmembrane potential during electric field stimulation of single cardiac cells,” Am. J. Physiol., 277 H351 –H362 (1999). 0002-9513 Google Scholar


V. P. Nikolski, A. T. Sambelashvili, V. I. Krinsky, and I. R. Efimov, “Effects of electroporation on optically recorded transmembrane potential responses to high-intensity electrical shocks,” Am. J. Physiol. Heart Circ. Physiol., 286 H412 –H418 (2004). 0363-6135 Google Scholar


O. F. Sharifov and V. G. Fast, “Optical mapping of transmural activation induced by electrical shocks in isolated left ventricular wall wedge preparations,” J. Cardiovasc. Electrophysiol., 14 1215 –1222 (2003). 1045-3873 Google Scholar


S. Dillon and M. Morad, “A new laser scanning system for measuring action potential propagation in the heart,” Science, 214 453 –456 (1981). 0036-8075 Google Scholar


G. Salama and M. Morad, “Merocyanine 540 as an optical probe of transmembrane electrical activity in the heart,” Science, 191 485 –487 (1976). 0036-8075 Google Scholar


M. Amino, M. Yamazaki, H. Nakagawa, H. Honjo, Y. Okuno, K. Yoshioka, T. Tanabe, K. Yasui, J. K. Lee, M. Horiba, K. Kamiya, and I. Kodama, “Combined effects of nifekalant and lidocaine on the spiral-type re-entry in a perfused 2-dimensional layer of rabbit ventricular myocardium,” Jpn. Circ. J., 69 576 –584 (2005). 0047-1828 Google Scholar


Y. C. Liu, C. Cabo, R. Salomonsz, M. Delmar, J. Davidenko, and J. Jalife, “Effects of diacetyl monoxime on the electrical properties of sheep and guinea pig ventricular muscle,” Cardiovasc. Res., 27 1991 –1997 (1993). 0008-6363 Google Scholar


A. Coulombe, I. A. Lefevre, E. Deroubaix, D. Thuringer, and E. Coraboeuf, “Effect of 2, 3-butanedione 2-monoxime on slow inward and transient outward currents in rat ventricular myocytes,” J. Mol. Cell. Cardiol., 22 921 –932 (1990). 0022-2828 Google Scholar


V. V. Fedorov, I. Lozinsky, E. A. Sosunov, E. P. Anyukhovsky, M. R. Rosen, C. W. Balke, and I. R. Efimov, “Application of blebbistatin as an excitation-contraction uncoupler for electrophysiological study of rat and rabbit hearts,” Heart Rhythm, 4 619 –626 (2007). 1547-5271 Google Scholar
©(2007) Society of Photo-Optical Instrumentation Engineers (SPIE)
Fujian Qu, Crystal M. Ripplinger, Vladimir Nikolski, Cindy Grimm, and Igor R. Efimov "Three-dimensional panoramic imaging of cardiac arrhythmias in rabbit heart," Journal of Biomedical Optics 12(4), 044019 (1 July 2007).
Published: 1 July 2007

Back to Top