Continuous roadmapping in liver TACE procedures using 2D-3D catheter-based registration.

PURPOSE: Fusion of pre/perioperative images and intra-operative images may add relevant information during image-guided procedures. In abdominal procedures, respiratory motion changes the position of organs, and thus accurate image guidance requires a continuous update of the spatial alignment of the (pre/perioperative) information with the organ position during the intervention.
METHODS: In this paper, we propose a method to register in real time perioperative 3D rotational angiography images (3DRA) to intra-operative single-plane 2D fluoroscopic images for improved guidance in TACE interventions. The method uses the shape of 3D vessels extracted from the 3DRA and the 2D catheter shape extracted from fluoroscopy. First, the appropriate 3D vessel is selected from the complete vascular tree using a shape similarity metric. Subsequently, the catheter is registered to this vessel, and the 3DRA is visualized based on the registration results. The method is evaluated on simulated data and clinical data.
RESULTS: The first selected vessel, ranked with the shape similarity metric, is used more than 39 % in the final registration and the second more than 21 %. The median of the closest corresponding points distance between 2D angiography vessels and projected 3D vessels is 4.7-5.4 mm when using the brute force optimizer and 5.2-6.6 mm when using the Powell optimizer.
CONCLUSION: We present a catheter-based registration method to continuously fuse a 3DRA roadmap arterial tree onto 2D fluoroscopic images with an efficient shape similarity.

Introduction

Transcatheter arterial chemoembolization (TACE) is a minimally invasive procedure to treat liver cancer (mostly hepatocellular carcinoma). In this procedure, a catheter is navigated toward a tumor via the femoral and hepatic artery, after which chemotherapeutic agents are injected. Currently, the interventionalist guides the catheter using single-plane 2D X-ray (fluoroscopy), visualizing only the catheter (Fig. 1). Frequently, contrast is injected to visualize the arteries. Computed tomography angiography (CTA) or 3D rotational angiography (3DRA) are used pre/perioperatively to visualize the tumor and feeding arteries. The navigation of the catheter using only 2D fluoroscopy is hampered by the inability to continuously visualize the arterial tree.

Fig. 1

TACE intervention overview (left) and fluoroscopy example (right)

Purpose of our work is to integrate information of the vasculature from pre/perioperative 3D images by fusing it with the intra-operative 2D X-ray images. Such an approach enables a continuous up-to-date roadmap and thus may improve the guidance during the procedure and consequently has the potential to reduce intervention time, radiation dose and contrast agent use.

2D–3D registration for improving image guidance has been studied in cardiac, cranial, abdominal and orthopedic procedures. An overview of 2D–3D registration methods is presented by Markelj et al. [13] and Liao et al. [10]. Following [13], 2D–3D registration methods can be classified as extrinsic, intrinsic and calibration based. Extrinsic methods use markers to register and update the registration [17]. Usually objects visible on X-ray (e.g., small beads) are inserted close to the region of interest before 3D image acquisition. Intrinsic methods rely on anatomical structures such as bones or the vasculature and are generally intensity, gradient, feature based or a combination of them [13]. In abdominal interventions, the vasculature and catheters are mostly the only structures visible on 2D X-ray images that can be used for registration. In cardiac [3, 4, 14, 21, 22], cranial [7, 15, 26, 27] and abdominal [5, 8, 9] interventions, vessel-based registration have been used between pre- or perioperative 3D/4D CTA (computed tomography angiography), MRA (magnetic resonance angiography) or CBCT (cone-beam computed tomography) and 2D DSA (digital subtraction angiography) or 2D fluoroscopies. The X-ray acquisition can either be single plane or biplane. Rigid as well as non-rigid registration approaches for aligning the vessels from 3D/4D pre- or perioperative images with those from DSA or fluoroscopy have been described. These approaches update the 3D vessels position with regard to the C-arm but do not enable a continuous roadmap of the 3D vessels because continuous contrast agent injection during the intervention would be harmful to the patient. Calibration-based methods can be used when the 3D perioperative image and the 2D images are acquired with the same device. For example, if the 3D position of the C-arm is known accurately, it allows alignment of intra-operative 2D X-ray images with perioperative 3D images. Atasoy et al. [2] and Ruijters et al. [23] use C-arm information to update the registration between perioperative 3DRA (or CBCT) and 2D X-ray. This approach has been demonstrated to work accurately in cranial procedures with no head movement. Utilization in abdominal interventions, however, is hampered by the respiratory motion, which invalidates the initial alignment. Atasoy et al. [2] proposed a semiautomatic method to follow one moving region of interest selected by a physician during the intervention (a part of a catheter) and to update the registration with this information. The transformation model contains in-plane translation to correct for shifts caused by respiratory motion. In cardiac interventions, Ma et al. [12] used manual calibration-based methods to achieve an initial alignment and then used features such as diaphragm/heart border, tracheal bifurcation or the catheter to correct for breathing motion. Another method was proposed by Luan et al. [11] for oral cancer treatment. They track the catheter tip with an electromagnetic sensor, reconstruct the catheter path and then register it with a preoperative image. Although tracking the 3D catheter tip tracking is valuable, breathing motion may hamper the reconstruction of the path in, e.g., the abdomen. Unlike most of the other methods, our previous method [1] performs a 2D/3D catheter-based registration using a 3DRA and the complete catheter visible in the 2D X-ray images. It does not require 2D angiographic images nor user interaction for the initial alignment. However, computation times were not interactive, hampering interventional use.

The major contribution of our current work is to propose a method for generating an automatic continuous roadmap during abdominal catheterization using 2D/3D registration with single-plane 2D X-ray images and perioperative 3DRA. This paper is an extension of our previous work: The metric for alignment has been improved, the registration is faster, and the evaluation has been performed on a larger set of data, containing synthetic images, clinical images and additional evaluation metrics.

Methods

The method is based on the registration of a 3D vessel tree with a 2D catheter shape. Therefore, in a preprocessing step, the arterial tree is extracted from the 3DRA image and the catheter shape position is determined from the single-plane fluoroscopic images. The extraction of the vessel tree itself is relatively straightforward for high-contrast 3DRA images. The segmentation of the catheter in the fluoroscopic images, albeit more challenging, has been subject of other studies [6, 16, 18, 19, 25, 28, 29]. These steps are not addressed in this paper.

Given the 3D vascular model and the 2D catheter centerline, the method consists of two steps (Fig. 2). First a shape similarity metric is used to find the vessel centerlines from the 3DRA that are most similar to the 2D catheter shape. Subsequently, a constrained 2D–3D registration is applied to find the corresponding rigid transformation between the 2D catheter and the 2D projections of the best-ranked vessel centerlines. In the following, we first define our coordinate systems and transformations, and then we describe each registration steps.

Fig. 2

Global overview: vessels/catheter extraction and 2D/3D registration

Definitions

We define the following coordinate systems (CS) for our setup in the intervention room (Fig. 3):Accordingly, the following coordinate transformations are defined: and are known for each X-ray image because the geometry and orientation of the C-arm are known. is unknown and is the result of our registration.

Fig. 3

Coordinate systems and transformations of the C-arm space

, denotes the world 3D CS, with the origin at the iso-center of the C-arm, and oriented along the C-arm in its default position

, the detector 3D CS (X-ray image plane)

, the 2D CS of the fluoroscopic image

, 3D CS of the 3DRA

, transformation from to

, cone-beam projection from to

, transformation that aligns the 3DRA space with the patient, from to

With the projection function (in homogeneous coordinates):we have a 3D point in the 3DRA space, , which can be projected on the fluoroscopic image space using the following equation:The catheter centerline extracted from a 2D fluoroscopic image is defined as an ordered set of points:where are 2D points at the center of the catheter in and denotes the tip of the catheter.

The blood vessel tree centerline extracted from the 3DRA is represented as a directed tree:where is the set of 3D points on the centerlines of the vessels in and is the set of directed edges between points. The root of is in the aorta, and the branches of the tree are in the liver.

We define a vessel centerline as an ordered set of points in , from any point along the directed edges to the root (Fig. 4):where in and is the root of .

Fig. 4

Terminology: blood vessel tree, vessel centerline and leaf vessel centerline

Similarly, we define the 2D projection of the 3D vessel centerline :Additionally, we define with , a linearly interpolated version of the projected centerline, with the length of the projected vessel centerline .

Shape-based vessel centerline selection

The registration is performed on a vessel centerline running from a leaf to the root (Fig. 4). Before performing the registration, the vessels that are the most likely to contain the catheter are selected. This selection is based on the shape similarity metric between the catheter and the projected 3D vessel. The metric quantifies the alignment of the tangent vectors of the catheter and the projected vessel. Therefore, it is not sensitive to the distance between centerlines. The underlying assumption is that the orientation of the vasculature changes little between the 3DRA and the fluoroscopy acquisition, which is valid for our application. The shape similarity for a vessel from a point is defined as:where is the length of the 2D catheter, is the identity matrix, and (resp. ) is the tangent of (resp. ) at the position . is the tip of the catheter, and is the possible location of the tip in the tree . with denoting the maximum similarity. As the catheter centerline is represented as a set of points, the integral over S is approximated by summing the dot products over all catheter positions, thereby interpolating the corresponding vessel positions (Fig. 5).

Fig. 5

Discretized sum of dot products between tangents of catheter and the vessel centerline with

In order to select the leaf vessel centerline for which the registration needs to be performed, for each leaf , the maximum similarity over all points in V(l) is determined:

Based on the values of , we selected the leafs with largest for the registration. When several leafs share the same common part with the catheter, only one is kept.

Rigid 2D/3D registration with forward projection

To register the 2D catheter with the vessel centerline, we need to find the rigid transform that yields the best match with the 2D catheter in . We decompose the transformation as follows:where is a rotation matrix with three unknowns (Euler angles, , and ), and is a translation matrix with three unknowns (, , ) with the translations aligned in . A translation along the projection axis in will only have a very minor effect in the projection. We therefore exclude from the registration parameters, leaving us with a five degrees of freedom transformation.

Our registration metric is the sum of distances between points on the catheter and closest points on the projected leaf vessel centerline. The catheter tip thereto is matched to the closest point of the projected vessel , where is a leaf selected thanks to the shape similarity. The distance between the catheter tip and the vessel centerline , given a rigid transformation T, is given by:Each next point of the catheter is similarly matched with a point of the projected vessel. To ensure continuity of the vessel (and simultaneously reducing computation time), the search range is limited to only a few points proximal to the point closest to the previous catheter point. Thus, let be the point matched with , then the distance to the subsequent catheter point is defined as:where are the consecutive points in , starting at , and is determined such that all points in that range are within a distance of .

Given these definitions, the final registration metric of our registration is a weighted sum of these distances (Fig. 6):where is a weight function and is the length of the catheter between and . As the registration accuracy close to the tip is more important than at the proximal part of the catheter, we use a weight to decrease the distance values that are further from the tip. We use a Gaussian with an offset:where is a parameter to control how fast the weight decrease (Fig. 7).

Fig. 6

Registration metric with the first closest distances of the two points of the catheter centerline

Fig. 7

Weight function to give more weight at the catheter tip with and various

This metric has two advantages: first, it is fast because we only look for the closest point in a specific neighborhood; second, by only matching points that are locally connected the continuity of the vessel centerline is respected.

Lastly, the final transformation is the one with the smallest cumulative distance:where represents the 5 degrees of freedom rigid transformation matrix.

Every selected leaf vessel centerline is registered, and the pair (, ) with the optimal similarity is kept.

Experiments

We performed two series of experiments. In the first one, we used clinical data from TACE interventions. As we do not have a ground truth available in these data, we rather evaluate the registration on the alignment of the vessels distal to the catheter tip using angiography and we investigate the effect of varying parameter settings. In the second experiment, registration has been performed on the same clinical data, but with the catheter position simulated. The simulation allows us to have a ground truth for the catheter position. The error in the localization of the registered catheter tip position was used for evaluation.

Data

We retrospectively acquired anonymized data of 19 TACE interventions (Table 1). The 16 first sets were acquired in the Erasmus MC, University Medical Center, Rotterdam, the Netherlands, between 2012 and 2014 in two different intervention rooms with angiographic C-arm systems (Xper Allura, Philips Healthcare, Best, the Netherlands). The last three sets were acquired in the Hôpitaux Universitaires Henri Mondor, Créteil, Paris, France, and the Ospedale di Circolo e Fondazione Macchi, Varese, Italy. For each intervention, we have a set of images consisting of one 3DRA image where the catheter was inside the hepatic artery, a set of fluoroscopic sequences with contrast agent and a set of digital subtraction angiographies (DSA) (Fig. 8). In these sequences, both the catheter and a part of the vasculature distal to the catheter tip are visible (by using the contrast agent). For each sequence, we selected the image with most of the vasculature visible and we manually annotated both the 2D catheter centerline and the 2D vasculature centerlines. The 3D arterial tree from 3DRA is segmented with a semiautomatic method based on thresholding and skeletonization [24].

Table 1

Number of clinical data

Patients	3DRA	Number of angios	Number of DSAs
01	Complete	4	2
02	Complete	2	3
03	Complete	0	2
04	Complete	1	4
05	Complete	3	1
06	Complete	2	2
07	Complete	5	1
08	Complete	2	3
09	Incomplete	4	3
10	Incomplete	2	3
11	Incomplete	3	2
12	Complete	2	2
13	Complete	1	4
14	Complete	0	2
15	Complete	3	2
16	Complete	3	5
17	Complete	1	2
18	Complete	1	0
19	Complete	0	2
Total	–	39	45

Fig. 8

Fluoroscopy with contrast agent (left) and DSA (right)

We divided our data in two different groups depending on the 3DRA acquisition: complete and incomplete acquisition. Incomplete acquisition occurs when the patient’s liver is not aligned with the C-arm rotation iso-center. In that case, the aorta and the hepatic artery are not visible in the 3DRA which hampers the registration.

Implementation

The method described was implemented in C++ and run on a computer with a 3.4 Ghz Intel Core i7. We set , the number of selected leaf vessel centerline to register, to 5.

In order to minimize our metric, we evaluated two different optimizers: a brute force and the Powell optimizer [20]. The brute force is exhaustive and is more likely to find the global minimum, whereas Powell is faster but because of its local search is more likely to converge to local minimum.

Our brute force optimizer has iterations, and for each iteration , the search space is centered at the minimum found in iteration with an interval size whereby coefficient was set to 0.5 (Fig.  9). Each dimension in the search space interval is subdivided in steps. We set the initial intervals of our brute force search to 50 mm for and and 7 for , and . These intervals are sufficiently large to capture breathing motion.

Fig. 9

Brute force optimizer in a 2D space with iterations, reduction coefficient , number of steps and initial interval size

For the Powell optimizer, we use a two-stage approach. We first optimize the in-plane translation and subsequently use that translation to initialize the full 5 degrees of freedom registration.

Clinical data and parameter optimization

In the first experiment, we investigated optimal parameter settings for the method and evaluated how well the resulting registration aligns the vasculature distal to the catheter tip. To determine optimal parameter settings and evaluate the effect of changing parameters, we applied the method with a large set of different settings (, and ) in a leave-one-out cross-validation scheme (determine the optimal parameter values over the set containing all patients except the one on which the evaluation is done). We tested the following settings; : , , , ; : , , , , ; : , , , ,  mm. Because the 2D catheter centerline and 3D blood vessel tree are discretized, we also investigated the effect of using different samplings: , and  mm between each point. Finally, computation time is recorded.

As we do not have a ground truth for the registration, we used the vasculature visible on the angiographies as a reference. Thus, for validating, we compared how well the projected 3DRA matches the arteries visible in the projection images. To this end, we projected the 3DRA vasculature on the 2D image with the registered and we computed the closest corresponding points distance between the projected 3D vasculature and the 2D vasculature. The most relevant region for the roadmapping is the area close to the catheter tip; we therefore only evaluated in a circular region (3 cm radius) around the catheter tip. Before computing the distance, we manually labelled the vessels such that distances are computed between corresponding vessels (Fig. 10). To prevent bias in this assignment, the manual annotation was done without registration, thus only using the initial projection of the 3DRA. To aid in the annotation, the observer could manually register 3DRA and 2D vasculature by changing translations and rotations of the 3DRA. Vessels that cannot be manually adequately linked were not used in the evaluation. Using the labelled corresponding vessels, we computed the closest corresponding points distance for each pair of vessels (excluding distances to endpoints). The distance was computed both from the 2D angio-vessel and from its registered projected 3D vessel pair (Fig. 11). Our evaluation metric for one image is the average of distances over all pairs of vessels.

Fig. 10

2D vasculature from contrast agent (left), 2D projection of 3DRA vasculature after manual registration (center) and manual-paired vessels (right): Same labels have same color vessels

Fig. 11

Closest corresponding points distance between paired vessels

Clinical data with a simulated catheter

When evaluating the registration on patient data, no accurate ground truth registration is available. We therefore also evaluated our method on simulated data with a known ground truth. To stay as close as possible to the reality, our simulation was based on the clinical data (Table 1) where we used the 3D extracted vessel tree in the 3DRA registered space (using ) as well as the fluoroscopic sequences with their projection information. We choose the position of the 3D simulated catheter tip in the 3D vessel tree such that it matches the catheter tip in the fluoroscopic image after a registration of the 3D vessel tree to the fluoroscopic image. Next, we extract a 3D simulated catheter centerline following the 3D vessel path from the tip to the root and project it on the fluoroscopic image using the angles from the fluoroscopic sequences acquired. Those ground truth projections are used to quantify the accuracy of the registration results. To simulate a smooth catheter that may be stretching the vessel and that may be partially outside the vasculature, we smooth the 2D projection centerline with a Gaussian kernel (with a standard deviation ).

We applied random transformations to the 3DRA volume, divided over three sets depending on the magnitude of the transformations and the Gaussian smoothing parameter (Table 2).

Table 2

Randomizations for the simulation experiments

	Slight	Moderate	Large
Translation x (in mm)	[-30,30]	[-30,30]	[-30,30]
Translation y (in mm)	[-20,20]	[-40,-20]∪[20,40]	[-50,-40]∪[40,50]
Translation z (in mm)	[-30,30]	[-30,30]	[-30,30]
Rotation α,β,γ (in ^∘)	[-6,6]	[-6,6]	[-6,6]
Catheter smoothing σ_simu (in mm)	[1,5]	[5,10]	[10,15]

In this experiment, we used the algorithm parameters obtained in the previous leave-one-out cross-validation. We quantified the Euclidean distance between the known projected 3D tip and the registered projected 3D tip. We also computed the longitudinal and orthogonal distances from the point of view of the known tip (Fig. 12).

Fig. 12

Evaluation of the registered tip position with the Euclidean distance, the longitudinal distance and the orthogonal distance

Results

The results of the leave-one-out cross-validation with the two optimizers Powell and brute force, 3 mm sampling and all images are presented in Table 3. The optimal parameter settings are consistent over the leave-one-out experiments. Based on these results, unless noted otherwise, we used 3 mm sampling and with Powell: , ,  mm and with brute force: , ,  mm.

Table 3

Optimal settings for Powell (left) and brute force (right) with 3 mm sampling after a leave-one-out cross-validation

Patients	λ	σ	d_max	Patients	λ	σ	d_max
01,02,03,04,05,06,08,09,10,11,12,13,14,15,16,17,19	0.2	20	40	02,03,04,05,06,07,08,09,11,12,13,14,15,16,17,18	0.1	80	20
07	0.1	40	40	01	0.2	80	10
18	0	40	20	10	0.1	80	30
				19	0.1	40	20

The average paired vessels distance results are summarized in Fig. 13. Results are grouped depending on the 3DRA acquisition. Compared with our previous method [1], using the same sampling, the new method has a median that is smaller. Also the brute force optimizer performs better than Powell and is more robust. With the complete 3DRA set, medians are around 5 mm for brute force and 6 mm for Powell. Varying the sampling density between 1.5 and 6 mm does not clearly affect the accuracy.

Fig. 13

Average distance between paired vessels (in 3 cm radius from the tip) for each image after registration (in mm). Comparisons between Powell, brute force optimizer and our previous method Ambrosini et al. [1]

In Table 4, we present the distribution of the best registered leaf vessel centerline among the ranked and selected leaf vessel centerlines. This table shows that the best registration result is generally obtained with the leaf vessel centerlines that ranked best using the shape-based metric . This demonstrates that the metric can effectively be used to reduce the number of potential vessels to register. The low percentage for the fifth ranked vessel also suggests that the choice of is a good compromise between registration speed and the robustness of the method. The first ranked leaf vessel centerline is also the one giving the best registration for 39–49 % of the images.

Table 4

Distribution of the best registered leaf vessel centerline among the ranked and ordered selected leaf vessel centerlines; with the optimal settings and the complete 3DRA set

	Sampling (in mm)	1st leaf (%)	2nd leaf (%)	3rd leaf (%)	4th leaf (%)	5th leaf (%)
Powell	1.5	40	29	10	8	13
Powell	3	49	24	12	9	6
Powell	6	45	24	10	16	5
Brute force	1.5	46	21	9	16	8
Brute force	3	39	24	8	19	10
Brute force	6	39	22	12	17	10

Figure 14 shows the paired vessels distance after the registration as function of the distance from the catheter tip. The median becomes less accurate after 3 cm between 1 and 5 mm.

Fig. 14

Average distance between paired vessels for all images (in mm) with optimal settings, 3 mm sampling and the complete 3DRA set. Paired vessels are grouped following their distance from the catheter tip (from 0 to 100 mm)

The sampling and the local distance in the metric are the parameters that affect the computation time. We show in Fig. 15 the relation between accuracy and computation time. Brute force is slower than Powell optimization. Our previous method [1] did the registrations with the samplings: 1.5, 3 and 6 mm in 95, 25 and 10 s, which is at least twice as slow as the brute force of our current approach. After 20 mm, does not seem to change the accuracy with brute force. Powell is less stable with both the sampling and .

Fig. 15

Average distance between paired vessels for all images in mm (left) and average time in second (right) according to the neighborhood distance , with optimal settings and the complete 3DRA set

Figure 16 shows examples of correct and incorrect registrations. We note that when there is a small part of the catheter visible on the image, the optimizers are more likely to yield misregistrations because of the lack of information. A correct tip position and distance metric do not imply a perfect match of the vasculature due to the deformation of the liver and the catheter.

Fig. 16

Projection of the 3DRA blood vessel (in green) with the catheter (in black) and the contrast agent (in purple). Initial position (left). Registered position with Powell (middle). Registered position with brute force (right). a The registration is correct. Here the catheter is long enough to give information. b The catheter part is too short. Powell registered with a good distance metric but the result is wrong. Brute force is correct. c The catheter tip position is correct for both optimizers. The vessels and catheter deformation prevent to have a perfect match. d Here the distance metric and the tip is correct with both optimizers but brute force rotates too much. e As a long part of the aorta is missing in the 3DRA, Powell stops in a local minimum while brute force is more exhaustive and reach the global minimum

We visually inspected all registrations from the complete 3DRA set with optimal settings and labelled them as correct, visually close and incorrect for both the registered tip and the registered vessels distal to the tip (Table 5). Visually close implies that the registration is sufficient to know where the catheter is in the 3D vasculature, while incorrect is of no use for the intervention. For each incorrect case, we also report the likely cause of failure (Table 6). The two main reasons of wrong registrations are a too small part of the catheter visible in the 2D X-ray image and large deformation of both the catheter and the vessels distal to it.

Table 5

Visual registration results with optimal settings, 3 mm sampling and the complete 3DRA set

	Powell (%)	Brute force (%)
Selected registered vessel
Correct	84	94
Incorrect	16	6
Match angio/registered vessels
Visually correct	38	35
Visually close	30	52
Incorrect	32	13
Registered tip
Visually correct	59	73
Visually close	20	22
Incorrect	21	5

Table 6

Registration error details among incorrect match with optimal settings, 3 mm sampling and the complete 3DRA set

	Powell (%)	Brute force (%)
Small catheter part visible	40	38
Large vessels and catheter deformation	25	25
Catheter shape not sufficiently distinctive	5	12
Rotate too much to best fit the catheter		25
Powell stops in a local minimum	20
Catheter only in the aorta (missing information)	5
Large part of the aorta is not visible in the 3DRA	5

Figure 17 shows the distance between the real tip in the simulated catheter (without smoothing) and the tip after registration. Figure 18 shows the results with catheter smoothing. Without catheter smoothing, for the brute force optimizer, the median of the Euclidean distance is below 1 mm whereas for Powell the distance is below 3 mm. With catheter smoothing, the registered tip is less accurate and less robust with both Powell and brute force optimization. The longitudinal and orthogonal distance are similar with slight, moderate or large transformation. Overall, the longitudinal distance is slightly more robust than the orthogonal.

Fig. 17

Euclidean distance , longitudinal distance and orthogonal distance between the real tip and the registered one (in mm) with no catheter smoothing, 3 different simulations (Table 2), optimal settings and 3 mm sampling

Fig. 18

Euclidean distance , longitudinal distance and orthogonal distance between the real tip and the registered one (in mm) with catheter smoothing, 3 different simulations (Table 2), optimal settings and 3 mm sampling

Discussion and conclusion

We proposed and evaluated a method that enables a continuous roadmap during abdominal catheterization. The method registers a 3D vessel model obtained from 3DRA imaging using a catheter that is extracted from single-plane 2D X-ray images with or without contrast agent. The method first selects the vessels using shape similarity and then rigidly registers the selected vessels to the catheter.

With the complete 3DRA set and optimal settings, the median of average paired vessels distances of the roadmap distal to the catheter tip and within a radius of 3 cm from the tip is 5.4 mm for the brute force optimizer and 5.2 mm for the Powell optimizer. The first selected vessel during shape similarity is used more than 39 % in the final registration and the second more than 21 %.

We investigated two optimizers for the registration approach: brute force and Powell. In our setup, with 200 ms computation time on average, the registration is real time with the Powell optimizer and a 3- or 6-mm sample interval. Though the brute force optimizer is slower, it could be improved with parallelization and a dedicated implementation. The brute force optimizer tends to be more accurate and robust than the Powell optimizer. Powell is more sensitive to the initial position of the registration (ends up in local minima) as well as the length and distinct shape of the catheter.

The simulation experiments with catheter deformation demonstrate that the registration is robust for both optimizers with slight deformation. They also show that larger deformation leads to less accurate registration. In the simulated data, the longitudinal distance from the tip shows how well the tip is registered along the catheter direction. This distance is more significative than the orthogonal and is slightly more robust.

An important source of error in our experiments was the lack of vessel information in the 3DRA (especially missing the aorta and the hepatic artery), partly caused by the retrospective nature of our study. Optimization of the 3DRA acquisition protocol could remedy this. Another source of registration errors is the lack of information in the 2D X-ray because only very short part of the catheter is visible. This shows the limitation of the method working with no prior knowledge other than the current image. This may be addressed, during the intervention, by slightly increasing the field of view, moving the patient table, or adding more a priori knowledge into the registration such as previous image registration transformations. If we take into account the previous registered transformations and the table motion (which in principle could be obtained from the C-arm system, but is not available in our acquired fluoroscopic images), the registration should have a better initialization and thus use a smaller search space and both Powell and brute force optimizers will perform more robustly while reducing computation time.

Most related 2D/3D registration methods register angiography with CTA or 3DRA. As the complete vasculature is visible on both 2D and 3D images and non-rigid registration is performed, they reach submillimeter accuracies. Our method, dealing only with the catheter visible on the 2D image, has lower accuracy. However, we are interested in improving guidance, and the fusion provides a continuous roadmap of sufficient accuracy to the clinician to reliably estimate the catheter tip position in the 3D vasculature. As far as we know, in abdominal studies, the presented method can be compared only with the method proposed by Atasoy et al. [2]. They evaluate their method with the overlap of the 3DRA vessels onto the catheter. In our case, the overlap is our distance metric so a comparison will be biased toward our method.

During abdominal catheterization, knowledge of the position of the tip in the 3D vasculature is of crucial importance. Table 5 shows a small percentage of incorrect registered tip positions. This implies that if we use the presented fusion method as a roadmap, combining any of the optimizers, the resulting fused visualization is sufficient to guide the interventionists in localizing the tip and identifying the subsequent bifurcations, also in case of slight misalignment.

A robust automatic 2D catheter segmentation is required after initialization to integrate our method into the interventional workflow. The accuracy of the segmentation will influence the registration method. For example, Heibel et al. [6] obtain a median error of real-time automatic catheter tracking 1.5 pixels for abdominal fluoroscopies. Those results are sufficiently accurate for our registration.

Registration studies often lack ground truth for clinical data. In our case, this also prevented us to evaluate the accuracy of the registration method directly. However, both the simulation experiments and the validation with angiographic images demonstrate the good performance of the method.

In our current setup, each registration is independent from previous registrations. During continuous roadmapping, only slight motion should occur between two registrations. In the future, we intend to use previous registration results to further improve the robustness (especially when the visible catheter part is too small to do an accurate independent registration) and to limit the computation time by reducing the space search. A source of registration errors was due to large vessels and catheter deformation. A non-rigid registration to match the catheter deformation could also improve the accuracy close to the catheter tip.

To conclude, we presented a catheter-based registration method to fuse continuously 3DRA roadmap arterial tree onto 2D fluoroscopic images. We evaluated our work with clinical and simulated data demonstrating an efficient shape similarity and a median accuracy, evaluated on close by vessels, of 4.7–6.6 mm and below 4 mm on simulated experiments.