Realistic Force Reflection in the Spine Biopsy Simulator Dong-Soo Kwon*, Ki-uk Kyung*, Sung Min Kwon**, Jong Beom Ra**, Hyun Wook Park** Heung Sik Kang***, Jianchao Zeng****, and Kevin R Cleary**** * Dept. of Mech. Eng. KAIST (Tel : +82-42-869-3042; Fax : +82-42-869-3210 ; E-mail: kwonds@me.kaist.ac.kr) ** Dept. of Elec. Eng & Computer Science. KAIST (Tel : +82-42-869-5434; E-mail: jbra@ee.kaist.ac.kr) *** Dept. of Diagnostic Radiology, Seoul National University College of Medicine (E-mail: kanghs@radcom.snu.ac.kr) **** ISIS of Georgetown Univ. (Tel : 202-687-8253; Fax : 202-42-784-3479 ; E-mail: cleary@isis.imac.georgetown.edu) Abstract This paper proposes a scheme to produce realistic force reflection in a needle insertion problem. The target system is the spine needle biopsy simulator for tumor inspection by needle insertion. Simulated force is calculated using the 3D human tissue model, the orientation and position of the needle, and this force is produced through the PHANToM TM device. To generate realistic force reflection, the directional force of the needle has been generated by related tissue model, and the other rotational force is generated using a pivot to keep the needle in the initial inserted direction after puncturing the skin. Since the haptic device has limitation for generating high stiffness and large damping, scale downed model and digital filter are used to stabilize the system. 1. Introduction Medical engineering technology has kept a pace with remarkable development of the computer hardware. Such development of medical technology has introduced many changes on the public well fare and on the medical education environment. Especially there is long stride in medical simulators that help acquiring medical skill and knowledge. The major change comparing with the previous simulator is haptic feedback. Watson et al. developed stereoscopic-haptic virtual environments [1]. Karl et al. studied real-time visually and haptically accurate surgical simulation. In this study, they implemented an algorithm to feel the surgical cutting [2]. The application area of force feedback is getting larger and larger. Many basic researches have been also performed. Heimenz and Maaβ experimented to find physical characteristics of different tissues [3,4]. Hopkins analyzed biomechanical characteristics of muscle and bone [5]. Margaret developed a 3 DOF (degree of freedom) joystick that provides force feedback. Kwon studied a 6 DOF haptic controller for telerobot system [6]. Massie developed a PHANToM TM device that has 6 degree of freedom and could make 3 DOF force [9]. Immersion Corporation developed Laparoscopic Impulse Engine [10]. These PHANToM TM and Laparoscopic Impulse Engine are commercialized and used in many medical simulators. In this paper, we propose a method to generate realistic force reflection during needle insertion into a hollow dummy in the spine biopsy simulator. A spine needle biopsy is a useful non-invasive operation to detect and verify a spine tumor. For this operation, a tissue is extracted from the spine of a living body using a hollow needle. However, since there are many critical organs near the spine, the spine needle biopsy operation requires complicated and accurate procedures. Hence, a lot of training is essential for reliable operations. Sunil and Popa measured the force from the needle during a needle biopsy and they tried to generate a same force reflection in the needle biopsy simulation [7,8]. Kevin applied PHANToM TM device to a 2D spine needle biopsy simulator [11]. The rest of the paper is organized as follows. We briefly introduce a surgical simulator for spine needle biopsy in Section 2. In Section 3, we present implementation details for force reflection. Finally, we present result and state our conclusions in Section 4. 2. Spine Needle Biopsy Simulator In this session, we introduce a spine needle biopsy simulator that provides realistic visual and force feedback to a trainee in the PC environment [12]. This system is composed of four parts: a 3D human model, a visual feedback part, a force feedback part, and an evaluation part. 2.1 3D Human Model A 3D human model plays an important role in the spine needle biopsy simulator because most parts of this system use this model. Therefore all kinds of data should be included. The 3D human model has two types of data. The one data type is common data. All voxels have common data such as CT density data, segmentation data, gradient data, and normal vector data. These common data are used by the visual feedback part, but especially segmentation data is used by force feedback part also. The other data type is individual data that depends on a kind of organ. Each organ has special characteristics such as 3D rendering parameters (color, opacity), and force generation parameter (spring constant). The 3D model uses 335 XCT slices having an intraslice resolution of 0.7mm and an inter-slice resolution of 1mm. This data set is given by Seoul National University Hospital. In order to get segmented 3D organ models based on these CT slices, a homemade semiautomatic segmentation tool is used. The target organs are selected as organs that should be needed in the procedure of simulating
Figure 1. The Spine Needle Biopsy Simulator spine needle biopsy. Doctors in the department of diagnostic radiology decide these organs. For the lumbar vertebra region, bone, lung, esophagus, artery, skin, muscle, and fat are selected as the target. For the thoracic vertebra region, vein and kidney are chosen additionally. Therefore total 9 organs are modeled in our 3D human model. 2.2 Visual Feedback Part Visual feedback includes all the GUI (graphic user interface) of the simulation system. Through the GUI of whole procedure, a trainee has a visual experience. In our simulator, there are two types of visual feedback. The one is 2D virtual CT console and the other is 3D view. The 2D virtual CT console box mimics the GUI and functions of a real CT system. In the 2D virtual CT console, a trainee can select a lesion that is a candidate for spine tumor. Then he simulates scanning CT images around the lesion. Many kinds of analysis tool are provided to the trainee in this console, such as finding a distance between two points, measuring an angle between two lines, zooming a specific CT image, etc. The second visual feedback part is 3D view. In this visual feedback part, 3D path planning tool is provided for the needle path from the skin to the lesion. Unlike previous 2D simulators, it supports realistic 3D volume rendering views interactively for 3D path planning. It provides axial, sagittal, coronal and 3D visualization images around the region also. Especially for the 3D visualization image, a trainee selects one image among color volume rendering result, gray volume rendering result, MIP (maximum intensity projection), and summed voxel projection. For the fast 3D volume rendering, a block-based technique is adopted for the efficient handling of large amounts of data and for the easy control of rendering parameters such as opacity [15]. 2.3 Force Feedback Part Force feedback increases the reality in the simulation system. The PHANToM TM haptic device of SensAble Technologies, Inc. is used and the spine biopsy needle replaces a stylus pen in order to implement directional and rotational force without any additional mechanical calculation. A haptic device is used to provide the force feedback to a trainee through the biopsy needle during simulation. For proper force feedback generation, the 3D human model data that contain force feedback parameters of organs near the spine is used and the haptic device is kept operating in 1kHz. A mock operation table that has same height as a real one is constructed and a dummy is laid down on the table. After that, a trainee is supposed to insert the spine biopsy needle into the hollow dummy on the table to feel the force in realistic environment. We explain more details in section 3. 2.4 Simulation This system is implemented by attaching a 3 DOF PHANToM TM device to a PC that has 600MHz Pentium III Dual CPUs and 512Mbyte RAM. A trainee simulates taking CT images of a patient and selects a lesion that is candidate for spine tumor in the 2D virtual CT console. In the 3D view, the 256x256x256 3D human model is rendered so that the trainee could plan the appropriate needle path from the skin to the lesion not to touch the other critical organs. After planning, he operates a needle biopsy virtually using haptic device. When he moves the needle, the haptic device captures the movements of the needle and these captured movements are rendered over the 3D human model. According to the position of the needle, the trainee feels different force reflection. For a 256x256x256 XCT abdomen volume data set, if he moves the spine needle only, it provides an update-rate over 25 Hz for visual-feedback and over 1kHz for force-feedback. When he touches the lesion with the needle tip, this simulation is over. Figure 2. The physical user interface of the developed spine biopsy simulator.
a correction force to keep the needle along the direction of movement( F ); and 3) an ambient force to compensate for C gravity( F ). This section explains how to calculate these G forces. Figure 3. Graphic user interface 2.5 Evaluation Part The evaluation part provides a full performance analysis to the trainee after simulation. This part includes information regarding the final position of needle tip, a list of punctured critical organs, a maximum deviation from the planned path, the number of trials, and the total time spent. 3. Implementation of Force Reflection Figure 5 illustrates force reflection mechanism for a needle biopsy simulator. F = F + F + F (1) R C Equation (1) shows the total force( F ) has three components: 1) the force required to penetrate tissue( F ); 2) R G 3.1 Tissue Penetration Resistive Force This force is simulated from the initial skin puncture until bone is touched. Since the maximum force that can be generated using PHANToM TM is 8.5N, the device has limited ability to simulate stiff bone[7]. In addition, since soft tissue, muscle and fat have nonlinear stiffness and damping properties, experiments and analysis are required for accurate modeling. Maaβ developed a model after measuring their elastic properties using ultrasonic wave[3]. However, the needle biopsy has not only an elastic process, but also a tissue fracture process. Also, the resistive force depends on the needle insertion speed[4]. In our simulator, experimental data were used to develop a tissue model. Sunil showed experimental results for needle punctures and developed an interactive lumbar puncture simulator with tactile feedback[7],[8]. The direction of the tissue penetration resistive force depends on the needle direction. The upper limit of this force is set to 6N because the PHANToM can only generate 8.5N as a peak force. Figure 6 shows Sunil s experimental data which were collected by a needle insertion with constant velocity, 1cm/sec[3]. We modeled each tissue as a spring without damping. Since position resolution at a needle tip is low and damping display of a haptic device is poor[9], damping is not considered. To simulate a realistic force profile, stiffness varies with the depth of inserted needle and the image segmentation value. Figure 4. Evaluation reports Figure 5. Force Reflection Mechanism
Figure 6. Needle Puncturing Resistive Force [7] Figure 8. Simulated Resistive Force the model is similar to the experimental data of Figure 6. Figure 7. Force Generation Algorithm Although the image segmentation classifies tissues as air, skin, muscle, lung, vein, artery, kidney and bone, the simulated force is derived from only air, skin, fat and bone because these are the tissues penetrated by the needle in most spine biopsies. If the biopsy path touches a lung or a kidney, the force falls to zero and a red warning sign appears on the screen. Figure 7 shows the algorithm used to calculate the tissue penetration force. Xr is a reference position to calculate force and this point is located on the needle. And Xp is a its tip position. k is stiffness of the tissue and vel is a speed and d is a depth of the needle insertion. Xr is determined with the inserted needle direction and the tissue that is punctured by needle. Variables, Xd, Xd2 are used to decide the time to change Xr. It is assumed that the skin is punctured when the needle is inserted into some desired depth. Figure 8 shows simulation results of the resistive force. Note the result of 3.2 Correction Force to Simulate Tissue Constraint In real needle biopsy, the motion of the needle is constrainted after insertion by the internal tissue. Moving the needle along the direction of the motion is not difficult, but it is very difficult to move to an orthogonal direction. Since the PHANToM(ver1.5) can only generate haptic translational forces (and not torques), an additional force is required to constrain the needle once it is inserted into the dummy. To solve this problem we make a hole at the surface of a dummy which is used as a pivot point to create rotational forces to constrain the needle along the inserted direction. At first, the direction vector of the needle is computed after the skin puncture. For calculation of the correction force, we assumed the following: The needle is touching only muscle or fat. The skin punctured position is a point and does not change. The rotational angle is much smaller than 1 radian. The first assumption is reasonable because the skin is very thin and there is only muscle and fat between the skin and the bone. The second assumption is trivial. Since the side surface of the needle is not sharp like tip and its motion is constrainted by internal tissues except the movement of inserted direction, the third one is acceptable. Using these assumptions, the correction force is calculated using the torque derived from the force at the needle tip. Figure 5 shows the initial needle direction vector Vi and the current needle direction vector Vc. Vi is derived by multiplying the needle inserted unit vector with the needle length. Vc is derived from the current tip point and the gimbal point of the needle. The correction force Fc can then be derived from the difference of these two vectors. The upper limit of this force is 2N since the haptic device can only generate 8.5N and the limit of the tissue penetration force is 6N. Since the tissues between skin and bone are muscle and fat, the stiffness for
Figure 9. Conversion of Correction Force correction force is determined from muscle and fat and this force is calculated by multiplying the correction vector by a stiffness and a constant. From experimental results, we verified that this correction force helps to keep the needle along the initial inserted direction can confirm that correction help to keep needle in initial inserted direction. Figure 9 shows the simulated force that must be applied by the haptic device. Fc_sim is derived as follows: Figure 10. Puncture Resistive Force over Time ll Fc _ sim = Fc (2) l 2 3.2 Compensation of Stiffness only Model and Stabilization for hard contact to bone Since the spring model is applied to virtual tissue, the needle jumps out when user releases it during insertion. And user even feels the pushed force during pulling needle out. Figure 7 shows the algorithm to simulate the cut feeling after the tissue puncture. The velocity value is used to detect current moving direction, and puling out and inserting is distinguished from this value. Using these, the resistive force falls to zero as the needle is pulled out. For more realistic modeling, it is required viscosity or viscoelastic tissue model. Vibration occurs when the needle contact to the bone. There are two reasons for this. First, mixed volume data with low resolution and abruptly changing stiffness can cause force discontinuities near the bone. This problem has been solved by using a digital filter. Second, the haptic interface has a limited ability to simulate high stiffness. Since the stiffness of the bone is high, it is not typically penetrated by needle puncture alone. As Colgate suggested passivity condition to guarantee stability[13], stiffness is scale downed to satisfy the stability condition. Figure 10 and Figure 11 show the tissue penetration resistive force and inserted depth of the needle over time. The system is stable when it touches the bone. The resistive force falls to zero as the needle is pulled out. 4. Summary The simulated force is calculated from the relationship between volume image data and the orientation and orientation and position of the needle. For more realistic force reflection in the spine needle biopsy simulator, the directional force of the needle has been generated by a tissue model. The rotational force is generated using a pivot to Figure 11. Inserted depth of needle keep the needle in the initial inserted direction after puncturing the skin. Since the haptic device employed has limited ability to generate high stiffness and large damping, a scale downed stiffness model and digital filter are used to stabilize the system. For more realistic simulator, tissue modeling based on biomechanics is required. And the device and controller to generate high stiffness and damping with stability is necessary. Acknowledgement This research was supported by the Korean Ministry of Information and Communication as a international joint project with ISIS Center at Georgetown University Medical Center, USA. REFERENCES [1] Watson K. et al., Development of stereoscopic-haptic virtual environments, 12th IEEE Symposium on Computer-based medical systems, pp. 29-34, 1999. [2] Karl D. Reinig et al., Real-time Visually and Haptically Accurate surgical simulation, http://www.uchsc.edu/sm/chs/research/mmvr4.html. [3] Maaß H.and Kuhnapfel U., Noninvasive Measurement of Elastic Properties of Living Tissue, 13th
International. Congress on Computer Assisted Radiology and Surgery (CARS '99), Paris, France, June 23-26, pp 865-870, 1999. [4] Heimenz L. and Litsky A., Puncture Mechanics for the Insertion of an Epidual Needle, Americal Society of Biomechanics 21th Annual Meeting, 1997. [5] David A. Hopkins, Musculoskeletal Biomechanics Research, http://dahweb.engr.ucdavis.edu/dahweb/ dahsite/dahsite.htm, 2000. [6] Dong-Soo Kwon, Ki Young Woo, and Hyung Suck Cho, Haptic Control of the Master Hand Controller for a Microsurgical Telerobot System, Proc. of the IEEE Int. Conf. on Robotics and Automation, pp. 1722~1727, 1999. [7] Sunil K. Singh and Mikael Bostrom, Design an Interactive lumbar Puncture Simulator with Tactile Feedback, IEEE International Workshop on Robot and Human Communication, 1993. [8] Dan O. Popa et al., Creating Realistic Force Sensations in a Virtual Environment: Experimental System, Fundamental issues and Results, Proc. of the IEEE International Conference on Robotics and Automation, 1998. [9] Thomas Massie, Design of a Force Reflecting Fingertip Stimulator, Bachelor s Thesis, Department of Electrical Engineering and Computer Science, MIT, 1993. [10] Immersion Co., Products and Technologies, http://www.immerse.com/products.html, 2000. [11] Kevin Cleary and Robert Greco, Development and Evaluation of a Spine Biopsy Simulator, Proceeding of Medicine Meets Virtual Reality, vol. 50, pp. 375-376, 1998. [12] J. Kim et al, Development of a Spine Needle Biopsy Simulator with Visual and Force Feedback, Proceeding of Congress on Computer Assisted Radiology and Surgery (CARS) 99, pp. 1052, Paris, France, June 23-26, 1999. [13] Edward J. Colgate and Gerd G. Schenkel, Passivity of a Class of Sampled-Data Systems: Application to Haptic Interface, American Control Conference, Baltimore, pp. 3236-3240, 1994. [14] Paul Gorman, Thomas Krummel et al, A Prototype Haptic Lumbar Puncture Simulator, Proceeding of Medicine Meets Virtual Reality 2000, IOS Press, vol. 70, pp. 106-109, 2000. [15] S. M. Kwon, J. K. Kim, and J. B. Ra, A Block-Based Volume Rendering Scheme for the Spine Biopsy Simulator, Proc. IEEE ICIP 1999, vol. 2, pp. 187-191, Kobe, Japan, 1999.