Buscar en
Revista Brasileira de Cardiologia Invasiva (English Edition)
Toda la web
Inicio Revista Brasileira de Cardiologia Invasiva (English Edition) Three-dimensional reconstruction of coronary arteries based on the integration o...
Journal Information
Vol. 23. Issue 2.
Pages 134-138 (April - June 2015)
Vol. 23. Issue 2.
Pages 134-138 (April - June 2015)
Original Article
Open Access
Three-dimensional reconstruction of coronary arteries based on the integration of intravascular ultrasound and conventional angiography
Reconstrução tridimensional de artérias coronárias a partir da integração do ultrassom intracoronário e da angiografia convencional
Cristiano Guedes Bezerraa,b, Gonzalo Daniel Maso Talouc,d, Carlos Alberto Bulantc,d, Breno de Alencar Araripe Falcãoa,b, José Mariani Jr.a,b, Pablo Javier Blancoc,d, Raúl Antonino Feijóoc,d, Pedro Alves Lemos Netoa,b,⿿
Corresponding author

Corresponding author. Avenida Dr. Enéas de Carvalho Aguiar, 44, bloco I, 3o andar, Hemodinâmica, Cerqueira César, CEP: 05403-000, São Paulo, SP, Brasil.
a Instituto do Coração, Hospital das Clínicas, Faculdade de Medicina, Universidade de São Paulo, São Paulo, SP, Brazil
b Hospital Sírio-Libanês, São Paulo, SP, Brazil
c Laboratório Nacional de Computação Científica, Ministério da Ciência, Tecnologia e Inovação, Petrópolis, RJ, Brazil
d Instituto Nacional de Ciências e Tecnologia em Medicina Assistida por Computação Científica, Petrópolis, RJ, Brazil
Article information
Full Text
Download PDF
Figures (1)
Tables (1)
Table 1. Geometric and anatomical descriptors. The variables in these equations are described in the studies by Bogunoviḿ et al.,12 Mut et al.,13 and Meng et al.14

Coronary three-dimensional reconstruction with the combination of intravascular ultrasound and angiography offers advantages over computed tomography angiography of coronary arteries. The authors aimed to present the pilot phase of the validation of a new model of three-dimensional reconstruction of coronary arteries.


This study used angiography and intravascular ultrasound examinations already performed by clinical indication in individuals with known or suspected stable coronary artery disease. Image processing, segmentation, and three-dimensional reconstruction were conducted following specific methodology. For geometrical characterization purposes, tridimensional center lines were obtained.


Three vessels were reconstructed: two left anterior descending arteries and one left circumflex artery. The vessel lumen volume and the overall plaque burden could be easily viewed with three-dimensional reconstruction. The geometric characterization revealed increased absolute values of length, tortuosity, curvature, and torsion, featuring a greater complexity of the center line of the diseased lumen relative to the center line of the external elastic membrane.


This new methodology, which integrates conventional angiography and intravascular ultrasound, has increased the practicality of the reconstructions, with a gain in volumetric accuracy of the vessel and overall visualization of key aspects of atherosclerotic disease, such as plaque remodeling and distribution.

Interventional ultrasonography
Three-dimensional imaging

A reconstrução tridimensional coronária com a combinação do ultrassom intracoronário e da angiografia apresenta vantagens em relação à angiotomografia de coronárias. Objetivamos apresentar a fase piloto de validação de um novo modelo de reconstrução tridimensional de artérias coronárias.


Foram utilizados exames de angiografia e ultrassom intracoronário já realizados por indicação clínica em indivíduos com suspeita ou diagnóstico de doença arterial coronária estável. O processamento das imagens, a segmentação e a reconstrução tridimensional foram realizados seguindo metodologia específica. Para fins de caracterização geométrica, foram obtidas as linhas de centro tridimensionais.


Foram reconstruídos três vasos, sendo duas artérias descendentes anteriores e uma artéria circunflexa. O volume da luz do vaso e a carga de placa global puderam ser visualizados com facilidade com a reconstrução tridimensional. A caracterização geométrica revelou aumento dos valores absolutos do comprimento, tortuosidade, curvatura e torsão, caracterizando uma maior complexidade da linha de centro da luz doente, em relação à linha de centro da membrana elástica externa.


Essa nova metodologia, que integrou angiografia convencional e ultrassom intracoronário, aumentou a praticidade das reconstruções, com ganho em acurácia volumétrica do vaso e visualização global de aspectos-chave da doença aterosclerótica, como remodelamento e distribuição da placa.

Ultrassonografia de intervenção
Imagem tridimensional
Full Text

Some studies use biplane angiography (AX) for three-dimensional reconstruction of the coronary artery lumen.1,2 However, AX is limited to defining the vessel lumen. Intravascular ultrasound (IVUS) is an accurate intravascular imaging technique for plaque quantification and location, since it allows for an accurate assessment from the adventitia to the intima. The combined use of IVUS and AX offers an interesting alternative for coronary three-dimensional reconstruction, outperforming coronary angiotomography (CT) in terms of accuracy. The ability to reconstruct the entire three-dimensional arterial wall allows for exploitation of the spatial characteristics of plaque. It can also be used for extracting geometrical features for studying the development and progression of atherosclerotic plaques, and for computational fluid dynamics models assessing the hemodynamics repercussions of the plaque.3

This study aimed to present the pilot phase of validation of a new three-dimensional reconstruction model.

MethodsAcquisition of images

For the development of this method, AX and IVUS, performed by clinical indication in patients with known or suspected stable coronary artery disease, from two institutions, the Hospital Sírio-Libanês and the Instituto do Coração, Hospital das Clínicas, University of São Paulo Medical School (FMUSP), both located in the city of São Paulo, SP, Brazil. After compilation, the anonymized data were analyzed by the research team in association with the Laboratório Nacional de Computação Científica (LNCC), where they were processed and integrated. The research project that incorporates this manuscript was approved by the Research Ethics Committee of Hospital Sírio-Libanês and of Hospital das Clínicas (FMUSP).

Reconstructions resulting from IVUS-AX were reprocessed, and geometric shapes of the lumen and the external elastic membrane (EEM) were extracted. Subsequently, the center lines and corresponding volumes of the lumen and the plaque were estimated. Finally, a set of geometric descriptors, calculated from the center line, were used to characterize the geometric shapes of the lumen and the EEM, and the difference between center lines was attributed to the presence of atherosclerotic plaque.

For the acquisition of IVUS images, the iLabTM system was used (Boston Scientific Corporation, Natick, USA), which allowed for the acquisition grayscale-scanned ultrasound images. The intracoronary ultrasound catheter Atlantis⿢ SR Pro was used (Boston Scientific Corporation, Natick, USA), consisting of a mechanical ultrasound catheter with a frequency of 40MHz. Automatic IVUS catheter pullbacks were performed in the catheter sheath at the rate of 0.5mm/s, beginning in the mid-distal third towards the artery ostium, with sectional tomographic image acquisition at a rate of 30 frames/sec.

The movement of the IVUS catheter and vessel curvature hampered the volumetric estimation of plaque burden. The movement of the catheter was corrected using only frames acquired in the diastolic phase of the cardiac cycle.4 However, the incorporation of vessel curvature required a location within the catheter space in order to estimate the correct position of the frames in cross-sectional images. Therefore, before the pullback, orthogonal AX was performed in the left anterior oblique and right anterior oblique projections, with cranial and caudal angulation, to estimate the catheter and the sheath spatial location.

Preprocessing and segmentation

To correct the catheter movement due to heartbeat, the frames associated with end-diastolic phase of the cardiac cycle were selected. The mathematical details were provided in a previous publication of this group.5 The choice of the end-diastolic phase was due to the minor displacement of the catheter, ensuring a more accurate spatial location during catheter pullback.

The study with IVUS provides a low signal-to-noise ratio (SNR), making vessel lumen and EEM recognition and, in turn, that of plaque burden, difficult. To reduce noise and preserve the structures of interest, an anisotropic diffusion method was employed.6

The vessel geometry extraction (segmentation) was performed using an adapted method of active contours.7 The energy functions used specifically in the process of minimization of active contours to extract the contour of the vessel lumen and EEM were, respectively,

the variables of these equations are described in articles by Maso Talou,7 Kass et al.,8 and Xu and Prince.9

To begin the process of segmentation, manual initialization of the contours was made in the first frame of the sequence of end-diastolic phases of IVUS. Subsequent frames took into account the contours of the previous frame to initialize the method. In case of a failed frame contour or if not sufficiently accurate, the operator redefined the required contours in that frame and continued the segmentation process.

Three-dimensional reconstruction

The orthogonal AX allowed optimal viewing of the catheter and of bifurcations. Observing the cardiac cycle, two frames of AX in diastolic phase were drawn: one with and the other without contrast. The lack of contrast allowed for a full visualization of the IVUS catheter. In addition to these conditions, the IVUS catheter sheath was segmented from the aortic sinus until the end of the transducer, with a biplane snake strategy.10 During snake initialization, the operator indicated only two points in AX: the location of the transducer and the aortic sinus, generating a straight line between these points. Then, the functional power of snake was optimized to get the segmentation of the catheter as a parametric curve c (s), s ⿿ [0.1].

The next step was to locate the contours of light (SiL) and EEM (SiE) along the trajectory of the catheter c (s). Using the acquisition time of each frame at the end of the diastolic phase of the cardiac cycle (ti), the position of each frame on the trajectory of the catheter was defined as:

where vϿb is the speed of automatic pullback and Fmax is the number of frames along the whole trajectory. The contour was displayed in the perpendicular plane for T, according to Frenet-Serret formula described to c (s).7

Finally, the rotation about the axis defined by c(s) was applied, R(θ), so that the encounter between the contours of IVUS and of the vessel lumen obtained by the contrasted projections of AX was increased. The θo angle, which maximized this encounter, was computed as follows:

where AR and AAX were the shadow of the projected lumen by the vascular geometry rebuilt and the vessel lumen by AX, respectively.

Geometric characterization

The three-dimensional reconstruction process provided the spatial configuration of the lumen and EEM (surface meshes) for each heartbeat time. Using ECG synchronization, the surfaces associated with the frames in the end-diastolic phase of the cardiac cycle were extracted.5 At this point, for each patient, there were two surface meshes to be characterized, one associated with the lumen and the other with EEM.

Both surfaces could be represented by their center lines,11 which are three-dimensional polylines representing the path taken by the artery, with punctual information of the radius. Next, all center lines were characterized using classical descriptors, such as bending and twisting along the polyline. The centerlines of the arteries were used extensively to characterize the vessels⿿ geometry.11⿿14


Fig. 1 depicts the three-dimensional reconstruction of three vessels, the first two corresponding to the left anterior descending artery (LAD), and the third to the left circumflex artery (LCx). Table 1 shows the anatomical and geometric descriptors used in this study.

Figure 1.

From left to right, three-dimensional reconstruction of left anterior descending artery (cases 1 and 2) and left circumflex artery (case 3). The lumen and the external elastic membrane are represented by red and blue center lines, and by pink and white volumetric filling, respectively.

Table 1.

Geometric and anatomical descriptors. The variables in these equations are described in the studies by Bogunoviḿ et al.,12 Mut et al.,13 and Meng et al.14

Geometric characteristics  Case 1: left anterior descendingCase 2: left anterior descendingCase 3: left circumflex
Number of branches  362
Lumen volume, mm3  333.23381.61,153.0
Plaque burden, mm3  304.61424.61496.73
  EEM  Lumen  EEM  Lumen  EEM  Lumen 
Length, mm  45.85  46.43  64.02  64.81  72.00  72.72 
Tortuosity12  0.12  0.13  0.15  0.16  0.15  0.16 
Mean radius, mm  1.88  1.28  1.28  1.21  2.43  1.96 
Aspect ratio13  24.39  36.27  35.43  53.39  29.65  37.19 
Curvature relationship, ÿ10⿿3  8.80  5.54  5.86  3.33  14.7  11.1 
Torsion value12  0.94  0.24  0.57  0.32  1.18  0.55 
Bending energy14  4.48  5.31  7.94  9.32  10.61  11.66 
Torsion energy14  334.4  229.50  779.5  759.46  1,103.13  870.65 
Mean curvature, mm⿿1  0.04  0.04  0.038  0.04  0.04  0.04 
Total curvature, mm⿿1  4.3  5.00  6.47  6.54  7.30 
Mean torsion, mm⿿1  -0.15  -0.14  -0.1  -0.04  -0.17  -0.16 
Total torsion mm⿿1  20.35  23.23  34.33  33.64  38.63  42.58 
Mean combined curvature,13 mm⿿1  0.20  0.22  0.23  0.22  0.25  0.26 
Total combined curvature,13 mm⿿1  22.03  24.78  36.43  36.21  40.94  45.01 

EEM: external elastic membrane.

Each case was filtered using the Oriented Speckle Reducing Anisotropic Diffusion (OSRAD), in order to reduce the noise (typical speckles in IVUS images), while preserving the vessel contours. The following parameters were established: Ͽ=0.7, δ=20, cmax=0.1, and cmin=0.5 for 40 repetitions. Then the IVUS segmentation process was developed, setting the parameters at α=1ÿ10⿿4; β=2; κ=2; κϿ=0.1; η=1.5; and γ=15, for vessel lumen, and α=1ÿ10⿿4; β=160; κ=1; and κϿ=0, for EEM segmentation.

The presence of stenosis resulted in disagreement between center lines of the vessel lumen and EEM (asymmetry). The volume of the vessel lumen and the overall plaque burden were obtained by three-dimensional reconstruction. Assuming little remodeling in healthy vessels, the EEM area was almost equal to that of vessel lumen. Thus, comparison of EEM with lumen characteristics provided a broad spectrum of atherosclerotic disease in the vessel in question. Table 1 shows the increased complexity of the center line of the unhealthy lumen relative to the center line of EEM (e.g., increases in length, tortuosity, total bending, and twisting). The difference of these characteristics stemmed from lesions, but could also be related to the presence of branches (bifurcations).


In comparison with other tests related to three-dimensional reconstruction such as CT, this model has some strengths: better spatial resolution of the lumen vessel, and greater accuracy in the delimitation of the vessel wall and in the details of the type of vascular remodeling. Unlike CT, which provides information of the entire coronary tree, reconstruction by IVUS-AX provides information about the study's target vessel. However, in the vessel in question, this model allows for visualization of coronary atherosclerotic disease in a more accurate way, through the precise location of the lesion and of plaque distribution. The synergistic integration of IVUS with AX allows for increasing the volumetric accuracy of the vessel lumen and the vascular wall, considering the curves and the spatial location of the vessel. This tool has the power to facilitate the recognition of characteristics of atherosclerotic disease, as exposed by IVUS (EEM data, remodeling, plaque burden, etc.), as well as detailing aspects insufficiently described by coronary angiotomography.

The geometric characterization of coronary arteries, as set out in this study, opens the door for future research and studies on the development and progression of atherosclerotic plaques. Due to the small diameter and movement of coronary arteries, it was not possible to directly measure the shear stress. In this context, the three-dimensional reconstruction is essential to calculate variables such as wall shear stress (WSS) or oscillatory shear index (OSI), by performing computational fluid dynamics simulations. Atherosclerotic plaques are related to shear stress, both as cause and as effect. In the development of atherosclerotic disease, positive remodeling prevents the growth of plaque into the vessel lumen. However, when the remodeling fails to offset this growth, the plaque will protrude into the vessel lumen by modifying the shear stress (cause and effect). The biological impact of changes in shear stress on plaque composition and on its vulnerability are not yet fully understood; nonetheless, it is known that the regions of the plaque subjected to high shear stress are more vulnerable and will have thinner and more fragile fibrous layers, predisposing to rupture. Therefore, a robust model of three-dimensional coronary reconstruction allows for the study of the relationship between shear stress, location, and plaque progression, as well as in-stent restenosis.

Moreover, through three-dimensional reconstructions, computational fluid dynamics can be applied to obtain local flow velocity and a functional evaluation of intermediate lesions, through numerical simulation of hemodynamic phenomena (virtual fractional flow reserve).

The classic technique of three-dimensional coronary reconstruction using a fusion of IVUS with biplane AX is known as ANGUS. The lumen contours and EEM are defined by IVUS, and the AX biplane data are used to determine the three-dimensional position of each IVUS frame. This procedure has been validated and is being used by several research groups, including for evaluation of shear stress in human coronary arteries. Compared to the ANGUS method, the method presented in this study exempts biplane AX, a technique not available in Brazilian hospitals. In addition, there was no need for online gating with electrocardiogram, which only retrieves a cardiac cycle phase.4 In fact, offline gating was applied based on image5 retrieval of 15-30 phases/study, allowing, in the future, four-dimensional reconstructions (space and time). Furthermore, this study used biplane snakes10 for the reconstruction of the catheter center line (core line) instead of the MacKay method.15 In this proposed alternative, the interaction with the operator consists of indicating the beginning and end of the pullback (points clearly distinguishable in IVUS images), whereas in the McKay method, over six non-coplanar points should be provided for each of the two images, thus making the initialization a more complex task. Additionally, biplane snakes allow a simple extension to use non-orthogonal AXs.16


The methodology presented for vessel reconstruction allows for greater data extraction compared to the classical reconstruction by CT angiography, as well as an optimization of the reconstruction process, in relation to the ANGUS method. The synergistic integration between angiographic studies and intravascular ultrasound increases the volumetric accuracy of the vessel and the global visualization of the key aspects of the atherosclerotic disease, such as plaque remodeling and distribution.

Funding source

This research was supported by the following Brazilian agencies: Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Fundação de Amparo à Pesquisa do Estado do Rio de Janeiro (FAPERJ), and Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) for the following centers: Laboratório Nacional de Ciências da Computação (LNCC) and Instituto Nacional de Ciências e Tecnologia em Medicina Assistida por Computação Científica (INCT/MACC).

Conflicts of interest

The authors declare no conflicts of interest.

D.L. Parker, D.L. Pope, R. Van Bree, H.W. Marshall.
Three-dimensional reconstruction of moving arterial beds from digital subtraction angiography.
Comput Biomed Res., 20 (1987), pp. 166-168
A. Wahle, E. Wellnhofer, I. Mugaragu, H.U. Saner, H. Oswald, E. Fleck.
Assessment of diffuse coro- nary artery disease by quantitative analysis of coronary morphology based upon 3-D reconstruction from biplane angiograms.
IEEE Trans Med Imaging., 14 (1995), pp. 230-241
R. Krams, J.J. Wentzel, J.A. Oomen, R. Vinke, J.C. Schuurbiers, P.J. de Feyter, et al.
Evaluation of endothelial shear stress and 3D geometry as factors determining the development of ath- erosclerosis and remodeling in human coronary arteries in vivo: com- bining 3D reconstruction from angiography and IVUS (ANGUS) with computational fluid dynamics.
Arterioscler Thromb Vasc Biol., 17 (1997), pp. 2061-2065
N. Bruining, C. von Birgelen, P.J. de Feyter, J. Ligthart, W. Li, P.W. Serruys, et al.
ECG-gated versus non gated three-dimensional intracoronary ultrasound analysis: implications for volumetric measurements.
Cathet Cardiovasc Diagn., 43 (1998), pp. 254-260
G. Maso Talou, I. Larrabide, P. Blanco, C. Bezerra, P. Lemos, R. Feijóo.
Improving Cardiac Phase Extraction in IVUS Studies by Integration of Gating Methods.
IEEE Trans Biomed Eng., (2015 Jun 24),
Epub ahead of print
K. Krissian, C.F. Westin, R. Kikinis, K.G. Vosburgh.
Oriented speckle reducing anisotropic diffusion.
IEEE Trans Image Process., 16 (2007), pp. 1412-1424
G.D. Maso Talou.
Segmentação de imagens IVUS via contornos ativos e reconstrução espaço-temporal dos vasos coronários assistida por angiografias. Dissertação de Mestrado.
Laboratório Nacional de Ciências da Computação, (2013),
M. Kass, A. Witkin, D. Terzopoulos.
Snakes: Active contour models.
Int J Comput Math., 1 (1988), pp. 321-331
C. Xu, J.L. Prince.
Snakes, shapes, and gradient vector flow.
IEEE Trans Image Process., 7 (1998), pp. 359-369
Canero C, Radeva P, Toledo R, Villanueva J, Mauri J.;1; 3D curve reconstruction by biplane snakes. In: Proceeding IEEE of 15th International Conference on Pattern Recognition, ICPR⿿00. Barcelona, Spain: IEEE Computer Society; September 3-8, 2000. v. 4. p. 563-6.
L. Antiga, B. Ene-Iordache, A. Remuzzi.
Computational geometry for patient-specific reconstruction and meshing of blood vessels from MR and CT angiography.
IEEE Trans Med Imaging., 22 (2003), pp. 674-684
H. Bogunoviḿ, J.M. Pozo, R. Cárdenes, M.C. Villa-Uriol, R. Blanc, M. Piotin, et al.
Automated landmarking and geometric characterization of the carotid siphon.
Med Image Anal., 16 (2012), pp. 889-903
F. Mut, S. Wright, G.A. Ascoli, J.R. Cebral.
Morphometric, geographic, and territorial characterization of brain arterial trees.
Int J Numer Method Biomed Eng., 30 (2014), pp. 755-766
S. Meng, S.H. Geyer, F. Costa Lda, M.P. Viana, W.J. Weninger.
Objective characterization of the course of the parasellar internal carotid artery using mathematical tools.
Surg Radiol Anat., 30 (2008), pp. 519-526
S.A. MacKay, M.J. Potel, J.M. Rubin.
Graphics methods for tracking three-dimensional heart wall motion.
Comput Biomed Res., 15 (1982), pp. 455-473
J. Yang, W. Cong, Y. Chen, J. Fan, Y. Liu, Y. Wang.
External force back-projective composition and globally deformable optimization for 3-D coronary artery reconstruction.
Phys Med Biol., 59 (2014), pp. 975-1003

Peer Review under the responsability of Sociedade Brasileira de Hemodinâmica e Cardiologia Intervencionista.

Copyright © 2015. Sociedade Brasileira de Hemodinâmica e Cardiologia Intervencionista
Article options