Computational and Mathematical Methods in Medicine
Volume 2018 (2018), Article ID 4534579, 7 pages
https://doi.org/10.1155/2018/4534579

Research Article

Method for Calculating the Bending Angle of Puncture Needle in Preoperative Planning for Transjugular Intrahepatic Portal Systemic Shunt (TIPS)

Xiaoli Zhu,¹ Zhao Ran,² Wanci Li,¹ Wansheng Wang,¹ Kangshun Zhu,³ Wensou Huang,³ and Xin Gao²

¹Invasive Technology Department, The First Affiliated Hospital of Soochow University, No. 899, Pinghai Road, Suzhou, Jiangsu 215006, China
²Department of Medical Imaging, Suzhou Institute of Biomedical Engineering and Technology, Chinese Academy of Sciences, No. 88, Keling Road, Suzhou, Jiangsu 215163, China
³Department of Minimally Invasive Interventional Radiology, The Second Affiliated Hospital of Guangzhou Medical University, No. 250, Changgang East Road, Guangzhou, Guangdong 510260, China

Correspondence should be addressed to Xin Gao; xingaosam@yahoo.com

Received 4 January 2018; Revised 7 March 2018; Accepted 11 April 2018; Published 30 May 2018

Academic Editor: Hiro Yoshida

Copyright © 2018 Xiaoli Zhu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Transjugular Intrahepatic Portal Systemic Shunt is a comprehensive interventional therapy for portal hypertension. During this intervention, puncturing from hepatic vein into portal vein is a difficult step. Selecting puncture needle with a proper bending angle is vital to accurate puncture. Thus, this prospective study provides a method to calculate the angle of the puncture needle using preinterventional contrast-enhanced CT imaging. According to the geometrical characteristics of puncture needle, Bezier curve equation was adopted to describe its bending part. By testing whether each point in a specific region satisfied the equation set of Bezier curves, the possible position of needle tip was obtained. Then, the bending angle of puncture needle was obtained by calculating curvature. The method was evaluated in 13 patients from 2 centers showing now a success rate of 100% and a duration of the procedure of 141 and 161 minutes. The method based on Bezier curve equation for calculating a proper bending angle of puncture needle was proven to be effective. And the clinical study is preliminary and additional work for clinical evaluation is necessary.

1. Introduction

Portal hypertension (PH) refers to a group of clinical syndromes characterized by abnormal changes of hemodynamics in portal vein system. It is also a major complication of cirrhosis and can even lead to death. Transjugular Intrahepatic Portal Systemic Shunt (TIPS) is a kind of comprehensive interventional therapy for PH [1, 2]. TIPS adopts radiation technology and uses special equipment (puncture needle with an bending angle, stent, etc.) to build a shunt between inferior vena cava and the branch of portal vein, so part of the bloodstream in portal vein can flow into inferior vena cava. Finally, portal vein pressure is decreased and severe complications caused by PH can be controlled and even prevented [3, 4].

Currently, the frequently used puncture tools are RUPS100 or RTPS100 puncture kits manufactured by Cook Company (USA). In the puncture kits, the puncture needle has only one type of bending angle. However, the causes of PH vary with patients, and the structure of intrahepatic vascular system, especially the spatial relationship between hepatic vein and portal vein, also varies with patients. Thus, the puncture needle currently used cannot meet the requirement of personalized treatment. Clinical practice indicates that choosing a puncture needle with a proper bending angle is necessary for accurate puncture of the portal vein [5]. Currently, the bending angle of puncture needle is adjusted by doctors according to their experiences. This process lacks accuracy. Once an improper bending angle is chosen, puncture cannot be performed according to the preplanned path and there is even a high risk of abdominal bleeding, resulting in intervention failure.

Therefore, the aim of this preliminary study was to provide a method that enables calculating the bending angle of the puncture needle on the basis of a previous CT imaging. The method was then evaluated in a small pilot study including 13 patients receiving a TIPS. The clinical study is preliminary and additional work for clinical evaluation is necessary.

2. Methods

2.1. Subjects

Between July 2016 and December 2016, clinical image data of 4 subjects from the First Affiliated Hospital of Soochow University (SU hospital) and 9 subjects from the Second Affiliated Hospital of Guangzhou Medical University were collected. These 13 subjects were supposed to undergo TIPS. This study was approved by the local Institutional Review Board. Among the subjects, eight patients were male and five were female. Their average age was 45 ± 4 years. Inclusion criteria include the following: (1) endoscopy confirmed esophageal or gastric varices bleeding, and the medical treatment cannot control bleeding effectively; (2) a history of portal hypertension in upper gastrointestinal bleeding and endoscopic or CT examination showed that esophageal or gastric fundus is still significantly varicose and not suitable for endoscopic therapy. Exclusion criteria included the following: (1) suffering from hepatic vein or intrahepatic portal vein or inferior vena cava occlusion; (2) patients with absolute contraindications to TIPS; (3) unwilling to accept the TIPS intervention patients

2.2. Path Planning

Iterative threshold segmentation combined with morphological operation was adopted to segment liver. Then, double threshold method was used to extract hepatic veins and portal veins in the segmented liver. Subsequently, 3D visualization was performed by using a surface-rendering method. Finally, the visualized image was put into 3D Slicer and an interventional physician was asked to designate two puncture points according to the anatomical spatial relationship between hepatic vein and portal vein. The point at which the needle punctured outside the hepatic vein was denoted as the Initial Point (IP). The point at which the needle punctured into the portal vein was denoted as the Terminal Point (TP). The preferred IP location was on the right hepatic veins, about 2 cm from the IVC. The preferred TP position was on the right branch of portal or on the left branch of portal, about 2 cm from the portal vein bifurcation. According to IP and TP, the puncture path was constructed (Figure 1).

Figure 1: Schematic illustration of puncture path planning.

2.3. Calculation of Bending Angle

In the plane where the puncture path was, puncture needle moved forward from IP to TP. Puncturing into portal vein (branch) was a critical step, after which guide wire was induced into the portal vein (trunk) for subsequent stent implantation. The whole puncture process depended on the bending angle of puncture needle. Only a proper bending angle can ensure the needle to move toward portal vein after puncturing outside hepatic vein and then puncture into the branch of portal vein. Otherwise, puncture process might be performed several times before success, and the risk of liver bleeding might be increased, resulting in intervention failure. As shown in Figure 2, the angle between the tangent of needle body and the tangent of needle tip was defined as the bending angle.

Figure 2: Puncture needle and its bending angle.

2.3.1. Construction of Puncture Plane

CT imaging system had self-defined coordinate system, which was denoted as . According to parameters (voxel position, voxel size, and layer distance) in DICOM file of CT image, the coordinates of voxel in can be determined. The coordinates of preplanned puncture points IP and TP in can also be obtained, which were and , respectively. A successful puncture operation meant that the puncture needle passed through IP and TP as well as a random point (RP) on the puncture path. These three points determined a puncture plane. According to equation of a plane passing through three points, that is,

we could obtain the equation of the puncture plane:

where , , , and were coefficients. Then, the unit normal vector of the puncture plane was , which was denoted as .

To conveniently use equations of curves for calculation, the 3D coordinate system was projected into a 2D coordinate system. As shown in Figure 3, the origin O of was projected onto the puncture plane and was the projection of O. Then, was taken as the origin of the new 2D coordinate system, called . -axis of the 2D coordinate system was the line passing through and IP. -axis was perpendicular to -axis and also in the puncture plane. Thus, the 2D coordinate system was built for the puncture plane:

Figure 3: Projecting the 3D coordinate system into a 2D coordinate system. The gray plane was the puncture plane.

(a) Origin : . Because , . Since was on the puncture plane, it should satisfy the equation of the plane (see (2)). Thus, we had (length of unit vector was 1).

(b) Unit vector of -axis:. was the vector passing through O′ and IP. Coordinates of were obtained in (a) and coordinates of IP were obtained during above-mentioned preplanning process.

(c) Unit vector of -axis: . Unit vector of -axis was the cross product of and , which was perpendicular to the puncture plane.

According to (b) and (c), the basis of was (, ). Any point inside a patient's body in the original 3D coordinate system projected onto the puncture plane had coordinates in the 2D coordinate system . This projection process was denoted as : →.

2.3.2. Curve Fitting

Bezier curve was an important parametric curve in computer graphics. It can use mathematical language to precisely describe any complicated graphs. Bezier curve was developed by Paul de Casteljau in 1959 using de Casteljau algorithm and was published by Pierre Bezier, a French engineer, in 1962 [6]. It had wide application in industrial design. Bezier curve had multiorder expressions. The second-order expression was defined as follows. As shown in Figure 4(a), three points including endpoints P₀ and P₂ and control points P₁ were given. Point Q₀ moved along line segment P₀P₁ and divided this line segment into two parts with a ratio t:1-t (first-order Bezier curve). Another point Q₁ on line segment P₁P₂ was chosen and it divided the line segment P₁P₂ into two parts also with the ratio t:1-t. Finally, a point B on line segment Q₀Q₁ was chosen and it also divided the line segment P₁P₂ into two parts with the ratio t:1-t. When t changed from 0 to 1, Q₀ moved from P₀ to P₁ and Q₁ moved from P₁ to P₂. During this process, the movement trajectory of B was a second-order Bezier curve and had the following equation:

where , , and are the coordinates of P₀, P₁, and P₂.

Figure 4: (a) Two-order Bezier curve. (b) Bezier curve matching puncture needle.

It was found that Bezier curve tool in vector mapping software can well fit the bending part of the needle (Figure 4(b)). It was thus assumed that the bending part of needle can be described by the Bezier curve. In Figure 4(b), the point at which the needle began to bend (bending point) was the starting point P₀ of Bezier curve, and the needle tip was the endpoint P₂ of Bezier curve.

Before intervention, the bending angle of puncture needle should be adjusted manually to meet the requirement of the preplanned path. When the bending point of puncture needle reached IP, the needle could hardly move forward. It was expected that the needle tip can reach TP. However, the currently used puncture needle had only one type of bending angle and length of the bending part was fixed. Moreover, the spatial relationship between IP and TP varied with patients. Therefore, the needle tip might not be able to reach TP. In fact, if the extension line of needle tip passed through TP, accurate puncture can also be achieved. As shown in Figure 5(a), it was assumed that the bending part of puncture needle can always be described by the equation of second-order Bezier curve when the bending angle was manually adjusted. In addition, the position of bending point and tangent at the bending point were assumed not to change, whereas the position of needle tip and tangent at the needle tip changed with the bending angle. The needle was rigid, so the length of its bending part was not supposed to change when the bending angle was adjusted. In the Bezier curve, this assumption conformed to the calculation rules. As shown in Figure 5(b), the change of the Bezier curve for puncture needle was simulated in Mathematica.

Figure 5: (a) Change of the bending part of puncture needle. (b) Change of corresponding Bezier curve in Mathematica.

According to above assumptions, two conditions were always satisfied when the bending angle was adjusted.

(1) The Extension Line of Needle Tip Must Pass through Puncture Point TP. When the bending angle was proper, puncture can be successfully performed according to the preplanned path. In that case, bending point must be where IP was, and needle tip or its extension line must be able to reach TP. Bending point or IP was the starting point P₀ of Bezier curve, and needle tip was the endpoint P₂ of Bezier curve. The position of the needle tip changed with bending angle, but its extension line always passed through puncture point TP. Thus, the movement pattern of the needle tip can be obtained. On this basis, we could obtain the slope of tangent at needle tip according to Bezier curve and further obtain the bending angle.

The slope at any point on the Bezier curve was

where , , and were coordinates of P₀, P₁, and P₂, respectively, of Bezier curve in 2D coordinate system. When the bending angle was adjusted, (i.e., coordinates of IP in ) remained the same, but and changed with bending angle. The slope at needle tip (t = 1, P2) was , where and were unknown. The slope of line segment TP P₂ was , where were coordinates of TP in .

(2) The Length of Bending Part of the Needle Remained the Same. First, we measured the length of bending part of the needle by using a soft thread. The bending part was found to have a length of 4.35 cm.

The bending part can be described by the equation of Bezier curve, so we had

According to integral formula

we can obtain the length of bending part of the needle. The calculation formula was

where , , and were coordinates of P₀, P₁, and P₂, respectively, in 2D coordinate system.

According to the above two conditions, the following equations were obtained:

When the bending angle was adjusted, the position of starting point P₀ did not change, whereas the positions of control point P₁ and endpoint P₂ changed. There were four unknown parameters but only two equations, so we cannot obtain the exact solutions of the equations. Even obtaining the movement trajectory of P₂ (i.e., the equation that satisfied) was difficult. In actual case, the range within which the bending angle changed would not be very large, since the needle was made of rigid material. Thus, we could use all the points in the circle whose center was at the initial position of the needle tip and radius was 2 cm to solve the problem. According to the resolution of the voxel space, all the points in the circle whose center was at the initial position of the needle tip and radius was 4 were tested if they satisfied (5). The step was set at 1. By doing so, we can obtain the position of P₂. The calculation was performed in Matlab and Mathematica.

According to the definition of the bending angle, we had

With the obtained position of , we could obtain .

3. Results

13 patients with portal hypertension underwent TIPS placement in two centers from July 2016 to December 2016. CECT was performed on each subject before TIPS. The CECT data were then used for preoperative planning.

After the original DICOM data were obtained, liver and blood vessels were segmented and were three-dimensionally visualized. Then an interventional physician was asked to designate puncture points for each subject. The coordinates of puncture points in for each subject were shown in Table 1.

Table 1: Coordinates of puncture points in chosen by interventional physician for each subject.

Subject number
IP
TP

1

2

3

4

5

6

7

8

9

10

11

12

13

After calculation, the bending angle of puncture needle for each subject was obtained (Table 2).

Table 2: Calculated bending angle of puncture needle for each subject.

Subject number
Bending angle (degree)

1
30.7

2
41.0

3
37.5

4
36.2

5
39.5

6
40.2

7
66.7

8
52.5

9
54.2

10
36.0

11
36.2

12
45.1

13
41.0

Because of a lack of intraoperative real-time imaging, the positions of actual puncture points IP and TP were a little different from those of the preplanned puncture points. Thus, the actual path was also different from the preplanned one. To evaluate the effect of bending angle on puncture efficiency, we adopted subjective and objective evaluations. The indicators were average puncture times, whether puncture needle deviated from the targeted point even punctured outside liver during intervention, whether there were puncture-related complications after intervention, and so forth. The clinical results were in Table 3.

Table 3: Clinical results of subjects in two research centers.

Research center
Average puncture times
Average procedure time (min)
Puncture outside liver
Puncture-related complications

(n=4)
3
161
1
0

(n=9)
2
142
0
0

refers to SU Hospital.

refers to SYSU Hospital.

Here, the procedure time referred to the time it took to perform the standard operations of TIPS. There were only two subjects whose results were not very good; the procedure time and puncture times of these two subjects were not as good as other patients. One was because the contrast time for CECT was improper, leading to the fact that the hepatic veins cannot be well segmented. Position of IP selected by interventional physician might be very inaccurate. The other was because the patient had undergone hepatectomy; hepatic tissue reconstruction led to difficulty in puncture.

4. Discussion

The bending angle of puncture needle is only one of the parameters determining the success of TIPS. Some parameters associated with the morphology of hepatic and portal venous systems as well as the positional relationship between them may also influence the accuracy and efficiency of puncture from hepatic vein to liver parenchyma and then into portal vein. However, the morphological structure of hepatic vascular system in patients with portal hypertension of cirrhosis is abnormal. This is caused by the atrophy of liver parenchyma and liver ascites after hepatic fibrosis. Thus, the puncture needle bending angle that we calculated is a reference value. Applying our method, probably many punctures will still not be successful because the needle bending may change during advancing the needle and the liver may move and twist during the puncture. Moreover, the procedure does not provide real-time guidance. Therefore, this does not replace sonographic real-time 3D guidance. It is definitely the gold standard. However, our method can be used as an adjunct, which may always improve the technique irrespective of other guiding attempts.

TIPS is a very difficult and highly risky intervention. In fact, it is thought to be one of the most complicated interventions by Peripheral Vascular Interventional Diagnosis and Surgery Hierarchical Directory. The process of puncture needle puncturing outside hepatic vein in a specific region and then precisely puncturing into portal vein branch is a difficult step of TIPS. To solve this problem, Image-guided TIPS intervention [7–9] and path planning method [10–13] have been developed. However, none of them performed calculation of the bending angle of puncture needle. In fact, the bending angle of puncture needle was often subjectively adjusted because of the lack of precise calculation, which might lead to the deviation of puncture needle from the targeted point and the damage to other tissues. In that case, the puncture needle needed to be taken out and reinserted in again, resulting in more puncture times. This paper proposed a method based on Bezier curve equation to precisely calculate a proper bending angle of puncture needle for TIPS. Bezier curve was used to describe the bending part of puncture needle, and Bezier curve equation was used to solve the problem of bending angle. Adjusting puncture needle to obtain a proper bending angle was a critical step during preoperative planning.

A path planning method based on numerical calculation model was proposed, where the critical puncture points can be selected according to preoperative 3D CECT. In addition, a proper bending angle of puncture needle can be calculated based on Bezier curve equation.

Clinical results showed that our method for adjusting the bending angle of puncture needle can well guide the puncture of the portal vein. In this case, the preplanned puncture path and the precisely calculated bending angle of puncture needle can improve puncture accuracy and reduce puncture times, ensuring an effective puncture operation. Knowledge of the bending angle has its own advantage whatever additional guidance is used.

The clinical study is preliminary; the method should be further evaluated intraoperatively and postoperatively on more cases. In the future, we will evaluate the performance of this method by using it in more clinical trials from multiple centers and will propose methods for calculating other parameters related to 3D puncture path.

Data Availability

The data presented in this study are provided by the First Affiliated Hospital of Soochow University and the Second Affiliated Hospital of Guangzhou Medical University. For the acquisition of raw data, please contact the corresponding author Professor Xin Gao at xingaosam@yahoo.com.

Ethical Approval

All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Declaration of Helsinki and its later amendments or comparable ethical standards.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.

Authors’ Contributions

Xiaoli Zhu and Zhao Ran are primary authors and contributed equally to this work.

Acknowledgments

This work was partially supported by Science and Technology Plan Projects of Jiangsu-Society Development Project (BE2017671), Science and Technology Plan Projects of Suzhou Health and Family Planning Commission (LCZX201704), and the National Natural Science Foundation of China (81571772).

References

R. S. Rahimi and D. C. Rockey, “Complications of cirrhosis,” Current Opinion in Gastroenterology, vol. 28, no. 3, pp. 223–229, 2012.

J. Rösch, B. T. Uchida, J. S. Putnam, R. W. Buschman, R. D. Law, and A. L. Hershey, “Experimental intrahepatic portacaval anastomosis: Use of expandable Gianturco stents,” Radiology, vol. 162, no. 2, pp. 481–485, 1987.

G. M. Richter, G. Noeldge, J. C. Palmaz et al., “Transjugular intrahepatic portacaval stent shunt: Preliminary clinical results,” Radiology, vol. 174, no. 3, pp. 1027–1030, 1990.

D. Bettinger, E. Knüppel, W. Euringer et al., “Efficacy and safety of transjugular intrahepatic portosystemic shunt (TIPSS) in 40 patients with hepatocellular carcinoma,” Alimentary Pharmacology and Therapeutics, vol. 41, no. 1, pp. 126–136, 2015.

C. F. Cuijpers, Tips for TIPS [D]; TU Delft, Delft University of Technology, 2015.

Y. Y. Xiong, “Collection of expressions of Bezier curve and their application,” in Proceedings of The Proceedings of The Sixth Geometric Design and Computing Academic Conference, China Society for Industrial and Mathematics: Geometric Design and Computing Committee, p. 4, 2013.

M. Darcy, “Transjugular intrahepatic portosystemic shunt: Techniques for portal localization,” Techniques in Vascular and Interventional Radiology, vol. 3, no. 3, pp. 147–157, 2000.

K. Farsad and J. A. Kaufman, “Novel image guidance techniques for portal vein targeting during transjugular intrahepatic portosystemic shunt creation,” Techniques in Vascular and Interventional Radiology, vol. 19, no. 1, pp. 10–20, 2016.

S.-H. Tang, J.-P. Qin, Q.-F. Shu, and M.-D. Jiang, “The imaging guidance for the portal vein branch puncturing in performing TIPS: recent progress in research,” Journal of Interventional Radiology, vol. 23, no. 7, pp. 640–643, 2014.

R. Adamus, M. Pfister, and R. W. R. Loose, “Enhancing transjugular intrahepatic portosystemic shunt puncture by using three-dimensional path planning based on the back projection of two two-dimensional portographs,” Radiology, vol. 251, no. 2, pp. 543–547, 2009.

J. Tsauo, X. Luo, L. Ye, and X. Li, “Three-dimensional path planning software-assisted transjugular intrahepatic portosystemic shunt: a technical modification,” CardioVascular and Interventional Radiology, vol. 38, no. 3, pp. 742–746, 2015.

J.-P. Qin, S.-H. Tang, M.-D. Jiang et al., “Contrast enhanced computed tomography and reconstruction of hepatic vascular system for transjugular intrahepatic portal systemic shunt puncture path planning,” World Journal of Gastroenterology, vol. 21, no. 32, pp. 9623–9629, 2015.

K. Li, Z. Tang, G.-J. Liu, and S.-X. Zhang, “Three-dimensional reconstruction of paracentesis approach in transjugular intrahepatic portosystemic shunt,” Anatomical Science International, vol. 87, no. 2, pp. 71–79, 2012.