Research Article  Open Access
Jing Guo, Yi Liu, Jin Wang, Chao Zeng, Jie Huang, Chao Liu, "The Design of Compact RoboticAssisted Needle Position System with MPCBased Remote Control", Complexity, vol. 2020, Article ID 5406084, 13 pages, 2020. https://doi.org/10.1155/2020/5406084
The Design of Compact RoboticAssisted Needle Position System with MPCBased Remote Control
Abstract
This article introduces the design and control performance of a lightweight, flexible, 4degreeoffreedom (DOF) parallel robot for percutaneous biopsy guided by computed tomography (CT). At present, the CT guidance method allows surgeons to quickly locate the lesion area; however, it is necessary to manually adjust the position of the puncture needle for insertion. In this paper, a threedimensional assisted method is used to infer the control input required to reach the target point through the kinematic model of the robot. A Kalman filter is designed to estimate model parameters and obtain a more accurate model. To further improve the control performance of the robot system, a modelbased control method—the model predictive control (MPC) controller—is used to increase the accuracy of the needle position in the developed robot system. In this way, medical efficiency is improved while reducing the burden on the surgeon.
1. Introduction
The surgeries have been advanced in the past decades thanks to novel surgical techniques [1], such as laparoscopy surgeries, in which a surgeon can operate a surgery through tiny holes on a human body for removing pathological tissue. However, the handhold tools still make the operation hard to be performed due to limited ergonomics and accuracy. Meanwhile, the surgeons may feel fatigued when operating for long time during the surgery. In order to improve the quality of surgical operation, the robot assistant technique has been introduced into the operation room which has eased the burden of the surgeon during the surgery [2–4]. Nowadays, the robotic systems have been developed and used for various surgeries such as biopsy, brachytherapy, and tumor ablation [5–7].
Before performing a minimally invasive surgical (MIS) task, the coordinates of the robot are usually registered to the medical image coordinate system, so that the robot is controlled by analysis of the lesion’s position for further operations. According to the characteristics of MIS, the pose of the robot endeffector before needle insertion during surgery becomes very important. The needle positioning system of the surgical assistant robot can be independent of the insertion action [8–10]. The surgeon can guide the robot’s needle to the insertion position of the lesion on the basis of the medical image through image guidance [11–16] and then manually insert the needle into the lesion area. Medical imaging, acting as a realtime feedback tool for the needle pose, is crucial in the process of needle positioning and control. In current needleguided systems, CT [11–13, 17, 18] is often used. Today’s commercially available minimally invasive robotic systems (MIRS) maintain and expand the flexibility of the surgeon’s hand [19–22]. NavarroAlarcon et al. [12] developed a 3DOF needle driver for biopsy. The positioning and insertion of the needle are achieved through three interfaces, the first two of which are aligned to the target, and the third is inserted. Koseki et al. [23] established a cooperative manipulation structure which uses the optical linear encoder to measure the needle position. However, manual remote operation control requires multiple scans using CT to obtain accurate targets. The robotassisted positioning approach can locally record the workspace and medical images of the robot system to reduce the risk of cancer induced by repeated CT scans [24]. This approach poses challenges to the compatibility and size of the robot system.
In the development of a surgical robot, another major problem is the network system in the remote operating system, which is caused by long distance or wireless link [23–27]. Sampling and delay are involved in this problem, which is a challenge for remote controlling of the robot. The solution to the delay problem in remote is discussed in [28–31]. However, the local optimizer on the slave side of the robot will also cause a delay [32, 33]. This part of the delay will also affect the remote operating system following the feedback signal. Dong et al. [34] modified standard functions to improve the response speed. Norizuki and Uchimura [35] embedded the MPC controller into the remote operating system to reduce the impact of the optimizer’s delay in the system.
In order to improve the quality of the surgery, the accuracy of the needle pose becomes crucial. The needle pose is affected by the robot control system and medical images. However, factors such as friction in the design of the robot structure and noise in the feedback of medical images are inevitable, which makes it difficult to accurately control the needle positioning [36, 37]. Thus, the closedloop control method becomes a good solution [38, 39]. In this paper, firstly, as shown in Figure 1, a novel lightweighted puncture robotic system is provided. A surgeon controls the medical robot through the master in the remote system. The medical robot can provide feedback of the endeffector pose and force to the surgeon. The size of the designed robot is high, long, and wide. In addition, the robot is controlled by four motors, which leads to the coupling of robot control. Due to the characteristics of the robot structure, specific constraints are proposed to decouple the kinematics. Finally, we develop a related needle positioning system on this robotic platform, which uses an MPC strategy to achieve the robotic local closedloop control of the needle tip position. This approach not only solves the problem of the time delay (TD) of the needle positioning system in local, but also reduces the cost and time of personnel training required for the operation of the robotic needle positioning system.
The structure of the paper is organized as follows: in Section 2, the overall robotic system design is introduced, including the derivation of constrained kinematic decoupling approach. Section 3 introduces model parameter estimation. Section 4 describes MPC control method. In Section 5, the experiments and the analysis of the experimental results are introduced. We conclude this article in Section 6.
In our previous work [40], we have designed a remote system which contains a master and a slave robot (see Figure 1). A surgeon can control the slave robot through the master side, and during a task, the desired position of the endeffector can be marked and recorded. The endeffector position of the slave robot is updated via medical imaging or other sensors, and such position feedback is sent to master. There are many ways to solve the remote communication scheme between the master and slave robots [28–31]. However, when the slave robot uses the optimized control method to improve the control effect, a delay will occur in the slave robot, which will affect the use of the remote operating system. Thus, the local controller design of the slave robot needs to be considered.
Because the robotic control system is a multiinput multioutput (MIMO) system, the controller selected for the slave robot should be able to adapt to the MIMO system. The proposed slave robot control system can be seen in Figure 2(a). The surgeon uses CT to mark the target lesion area or trajectory in the robot workspace. The surgical robot obtains the control trajectory which serves as the reference input for the MPC controller through kinematic model. The MPC controller optimizes the current control input to the surgical robot according to the reference trajectory and the reference state of the surgical robot. Furthermore, in order to improve the control performance, the statespace model of the robot required by the MPC controller should be accurate. In this paper, according to the structure of the surgical robot, the kinematics model is divided into a linear part and a nonlinear part. The linear part uses the model parameter estimator to identify the model parameters, and the nonlinear part is linearized using the Taylor formula. As shown in Figure 2(b), the model parameter estimator of the linear part is composed of an order judger and a parameter estimator. The order judger is used to judge the order of the linear system. The Kalman filter is selected as the parameter estimator to estimate the parameters of the linear system.
2. Kinematic Modeling
In this section, the modeling process of surgical robot kinematics will be described. The robotic kinematic is analyzed for the establishment of model parameter estimator and MPC controller. A suitable kinematics model is constructed by analyzing the mechanical structure of the parallel puncture needle robot [40]. Figure 3 shows the threedimensional mechanical structure and sectional structure of the surgical robot. The blocks, i.e.,  and , can be moved by controlling the angle of four motors. When the positions of the two blocks change, the position and insertion direction of the puncture needle will change accordingly. Among this structure, the vertical distance between the two blocks and the length of the puncture needle below the  are both constant.
In order to obtain a suitable kinematics model, the coordinate relationship as shown in Figure 4 is established. The fixed world coordinate system is at the initial center position of . is the coordinate system of , which is used to describe the moving distance of . Similarly, the value of the coordinate system is the moving distance of . Coordinate system indicates the position of the tip of the puncture needle.
2.1. Forward Kinematic
The purpose of forward kinematics is to obtain the position of the tip of the puncture needle , , and , and the forward kinematics of the surgical robot can be obtained using geometric theorem. As shown in Figure 4, firstly, the rotation of the motor will drive the block. The moving distance of the two blocks , , , and will change the angle of the puncture needle .where is the pitch of the screw connected to the block; is the rotation angle of the motor, and the relationship between the rotation angle of other motors and its block coordinates is similar; is the angle generated by two blocks on the  plane; next, for coordinate system , the change in the position of the tip of the puncture needle , , and can be described by , , and the constant length of the needle :
Finally, the position of the puncture needle’s tip , , and relative to the fixed world coordinate system can be obtained as follows:
2.2. Inverse Kinematic
The surgical robot is a parallel structure; hence, there is a coupling relationship between the control inputs. To address this problem, constraints are first defined to decouple the kinematic model. From the point of view of the positioning range of the puncture needle, it is better to have a larger range that the puncture needle can reach in a limited space. Thus, the first constraint is that the two blocks should move in opposite directions. On the other hand, from the view of the lesion area marked by surgeon, it can be seen that the position of the tip of the puncture needle in the  plane is of more significance. In order to obtain a feasible solution for the motor control input, the second constraint is defined. That is, when  moves to the limit position and still cannot reach the lesion area,  is then moved. The inverse kinematics model is constructed using geometric theorem based on the constraints:andwhere is the largest distance that the tip of the puncture needle can reach when only  is in motion; and are the lesion’s target point coordinates marked by the surgeon. However, no complete feasible solution can be obtained yet. Then, the solution of the coordinates of  within the limit is discussed:where is the angle between the block and puncture needle:
So far, a feasible solution within the limit from equations (10) and (12) is obtained:
Finally, the feasible solution in the other case will be analyzed to obtain the complete feasible solution. When the lesion area exceeds the limit , the coordinate values of  are fixed. Therefore, according to the angle between the block and the puncture needle, the relationship between and can be acquired:
The solution of coordinate value of  can be acquired by equations (15) and (16):where , , , and are known parameters as follows:
Then, the feasible solution that exceeds the limit is obtained by equations (11) and (17):
Now, the complete solution of the kinematic model of the surgical robot has been obtained by equations (14) and (19). The corresponding motor control angle can be obtained through the lesion area marked by the surgeon.
3. Model Parameter Estimation
After the desired control input is obtained, it cannot, in fact, directly act on the robotic controller. In order to make the statespace model of the robot closer to the actual model, the transmission process of the surgical robot will be divided into two parts for modeling. The first part is the linear system model of the block movement, and the second part is the nonlinear system model of the puncture needle movement. The determinant ratio is used to estimate the order of the linear system model, and a Kalman filter is established to estimate the linear system model parameters. Taylor formula is used to linearize the nonlinear system model into the controller.
3.1. Linear Part
Firstly, the linear system model is a MIMO system; however, each input corresponds to only one output. Thus, the MIMO system is equivalent to four singleinput singleoutput (SISO) systems. The order and parameters of the SISO system are considered:where is the Gaussian white noise with mean zero, is the input signal, and is the output signal. The input and output signals at each moment are known, but parameters are unknown. Then, input signal and output signal are used to construct determinants for comparison to obtain the ratio :where is the estimated order and is a matrix composed of input and output and has the following relationship:where is the sampling length of the signal. gradually increases from one. If only has a significant growth compared to , is the estimated order.
Next, a Kalman filter is designed to dynamically identify the system’s unknown parameters and . The Kalman filter is usually used to optimally estimate the state value of the statespace equation, which is composed of the unknown parameters:where both and are noise matrix composed of Gaussian white noise with zero mean. is the process noise drive matrix. is the state matrix at the th sampling moment. is the state transition matrix at the th sampling moment. is the output matrix at the th sampling moment. And we assume that the variances of the two noises are and . Then, the update equation of Kalman filter can be derived based on Euclidean theorem:where is the optimal estimate at the  sampling moment, is the error covariance of the current sampling time, and is the gain of the Kalman filter, which is used to update the values of and . In this paper, the estimated state matrix with unknown parameters is constructed:
Then the system’s state equation is obtained as follows:where is composed of . We use the observed input and output data to construct the matrix and state equation:
Finally, state equation (28) can be updated by equation (25) to obtain the actual system parameters. The SISO systems of four motors are combined to obtain the MIMO system of the linear part:where the linear part’s MIMO system is expressed by rewriting the state equations (32) and (33); consists of samples of the output before the th time step; , , , and are composed of each motor model parameter matrix ; similarly, other parameters of the MIMO system can also be expressed by equations (30) and (31).
3.2. Nonlinear Part
The nonlinear part of the system mainly represents the relationship between the rotation of the puncture needle and the coordinates of the block. The relationship between them through equations (2)–(9) can be obtained as follows:where is their nonlinear mapping. Then, a firstorder expansion of Taylor’s formula is performed for to obtain the following linear relationship:where is the target reference values obtained by the inverse kinematics model, is an infinitesimal quantity, and is the linearized Jacobian matrix. Finally, the complete input and output linear system model of the surgical robot can be obtained through equations (33) and (35) as follows:where
4. Model Predictive Control
Model predictive control (MPC) is now a wellknown modelbased control strategy, and it has gone through a period of development [41–43]. Taking the statespace model as an example, MPC uses the state of the model at the current time and the target state to estimate the control sequences at the next time under constraints. Compared with proportionalintegralderivative (PID) controller, MPC can predict the future behavior of the system. The optimal future control sequences are estimated by the known data. These sequences are evaluated at each sample time by the optimizer with constraints.
MPC consists of three components [44]: (1) a loss function, (2) constraints in the form of equality and inequality, and (3) initial conditions. The purpose of the loss function is to minimize the distance between the target state and the state at the current time [45], thereby minimizing the energy required by the control system. The equality constraints consist of the system’s state equations. The limits on the state variables are considered as the inequality constraints. If all constraints are linear, the problem forms a convex set. Meanwhile, if the control system is also convex, MPC can be solved as a convex optimization problem. Similarly, the MPC can be solved as a quadratic problem if the control system is a quadratic problem.
MPC uses the system’s model to predict the future control sequences [46]. As shown in Figure 5, the optimizer in MPC is to get the best control input. The linear timeinvariant system will be considered in the discrete time domain as shown in equation (37). However, in equation (37), the state and input variables are added by two matrices to obtain the final output. To comprehensively consider the state and input variables, equations (37) and (38) are rewritten as follows:where
Equations (42) and (43) are statespace equations that consider both state and input variables. We mainly discuss the optimization problem of equation (42), namely, the state space. This equation can make predictions about the future state, that is, the prediction step required for optimization can be obtained:wherewhere is the predicted step size of the system and is the control step size of the system. Next, the reference trajectory is used to construct the following loss function:where is the weight matrix. Considering that the control input cannot change too much within a period of time, the following loss function is constructed:where is the weight matrix. And the objective function to be optimized can be obtained by equations (53) and (54), which can be seen as a quadratic programming problem:
Finally, the constraints of the control variable of the loss function to construct a standard form of quadratic programming will be set:where and are the lower and upper bounds of the variable . Then, the best control sequences can be obtained by equation (45).
5. Simulation Setup and Results
This section mainly describes the simulated experimental settings and experimental results about this article. This article uses MATLAB to build a remote subsystem, which contains the MPC controller and parameter identification system. Under the function of the parameter identification system, the MPC controller obtains the robotic kinematics model realistically and outputs the optimal control sequence. In the first place, as shown in Figure 2, there is actually a TD when the MPC controller is applied to the system; thus, in the experiment, a delay system will be considered for comparison. We will, in addition, divide it into two experiments: (1) order and parameter identification of the system’s linear part and (2) performance of MPC controller and PID controller in delay system.
5.1. Model Identification Experiment
Above all, only the model identification of the linear part of the system is considered. And the order needs to be identified through offline data. The system’s order can be obtained by equation (21). In this paper, the system’s order , the pitch of the screw , and the largest movement distance of the block . Next, the Kalman filter is designed to identify the system parameters. The state equation of the system is rewritten to identify the parameters of the system, and the result is represented by in equation (28). The results of parameter identification are shown in Figure 6.
It can be seen from Figure 6 that the three parameters can stabilize at a certain value , , and within the step size. However, it is not known whether a system composed of stable parameters can approximate the real system. Thus, as shown in Figure 7, the output of the actual system is compared with the output of the identified model. The figure shows the target curve, the output curve of the identification model, and the error curve between them. The value of the error curve floats above and below zero.
5.2. MPC Performance
First, the model of the linear part of the control system can be obtained by equation (28), and the parameters of the model are obtained in the previous experiment, which contains the effective length of the needle and the distance between blocks . And the linearized model of the nonlinear part of the control system can be obtained by equation (35). Then, the two models are combined by equation (37) to obtain a complete model of the control system. In addition, in order to comprehensively consider the control input and system state, the statespace equation of the control system is rewritten by equations (42) and (43). Finally, the prediction equation (49) required by the MPC controller is obtained by equation (42), which is used for rolling optimization in the MPC controller. The constraints of quadratic programming and are constraints on the difference values of control sequence. The optimizer of the MPC controller outputs the optimal control sequence through equation (55). In this paper, traditional PID controller is designed to compare with the MPC controller. In fact, as the controller takes time to run, there is a delay in the output of the controller, which is equivalent to the control delay system. Thus, we set different TDs in the control system. We set the weight matrix and in the MPC controller as the identity matrix and the parameters of the PID controller , , and .
Figure 8 shows the axis motion trajectory of , which illustrates that with increase in the delay time, the response time of the system also increases and the desired trajectory can be tracked by the MPC controller. As the delay increases, the response speed of the system also becomes slower, and the output of the MPC controller will produce overshoot. However, the PID controller does not perform well in the case of TD. The output of the PID controller diverges when TD = , making the system unstable. When TD = , the PID controller makes the system gradually stable, but the setting time and the overshoot exceed the constraints of the surgical robot.
Then, the output of the position of the tip of the puncture needle will be analyzed. As shown in Figure 9, in the openloop state without controller, the output value of the system has a large fluctuation around the desire trajectory. The closedloop system formed by the MPC controller can track the desired trajectory with TD, which will produce overshoot in the delay system. However, the closedloop system formed by the PID controller is sensitive to TD. When there is a large TD in the system, it will produce a larger overshoot than the MPC controller and even make the system unstable.
Finally, in order to more clearly show the advantages of the MPC controller, a differential operation is performed on the state value of the sine wave shown in Figure 8 to obtain the distance from the target value shown in Figure 10. In the figure, the output of the MPC controller applied to the delay system produces phase shifts at different distances relative to the expected value. However, when the PID controller acts on a nondelayed system, phase shift occurs. As the delay increases, even the output diverges. The above experimental results show that the effect of MPC controller is better than that of the PID controller. The error variance is shown in Table 1, the desired error variance is zero, and the closer to the ideal error variance, the more ideal the system’s outputs are. The error variance of the MPC controller in the TD system is smaller than the error variance generated by the PID controller. And when , the error variance generated by the MPC controller is even smaller than the error variance generated by the nondelay system with noise. From the result, we can get that for the parallel mechanism robot involved in this article, MPC has a superior control effect because MPC can control the MIMO system.

6. Conclusion and Future Work
We introduced a surgical robotic puncture needle positioning platform that can improve the efficiency of biopsy. Kinematic analysis of the platform illustrates that the puncture needle can work in the working space through multiple angles that can be obtained by decoupling the robotic kinematic. In this paper, two constraints are defined for decoupling the kinematic, which ensure the robotic largest workspace. The designed MPC controller makes the control flow of the positioning device form a closed loop. And the state equation of the MPC controller needed is divided into two parts: linear part and nonlinear part. In order to guarantee the simulation of robotic system closer to the reality, the Kalman filter system parameters identification approach is used to obtain the system’s parameters in the linear part. In the nonlinear part, Taylor formula is used to linearize the system, and the linearization system only contains a Jacobian matrix. Next, the two parts of the system are combined into a complete robot system state equation. In simulation, the closedloop control significantly improves the stability of robot needle positioning. And the performance of MPC controller is better than that of the PID controller. Finally, the system is open source in Github (https://github.com/tKsome/MPCPunctureRobot) to reduce the amount of development time for other researchers involved in the control of the medical robotic MIMO system. The system is flexible and may be useful for most imagebased intervention procedures. In future work, we will test the performance of the proposed approach on the realwork robotic platform. Furthermore, a force control strategy will be proposed such that a force feedback can be provided to the surgeon through remote operation.
Data Availability
The data included in this study are available upon proper request by contact with the corresponding author.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the National Natural Science Foundation of China with grant no. 61803103 and Science and Technology Program of Guangdong with grant nos. 2019B010140002 and 2019B010142001.
References
 P. Kulkarni, S. Sikander, P. Biswas, S. Frawley, and S.E. Song, “Review of robotic needle guide systems for percutaneous intervention,” Annals of Biomedical Engineering, vol. 47, p. 7, 2019. View at: Publisher Site  Google Scholar
 B. Schell, K. Eichler, M. Mack et al., “Roboterunterstützte punktion in einem hochfeldkernspintomografen—erste klinische ergebnisse,” RöFo—Fortschritte auf dem Gebiet der Röntgenstrahlen und der Bildgebenden VERFAHREN, vol. 184, no. 1, pp. 42–47, 2012. View at: Publisher Site  Google Scholar
 N. Hungr, I. Bricault, P. Cinquin, and C. Fouard, “Design and validation of a CT and MRIguided robot for percutaneous needle procedures,” IEEE Transactions on Robotics, vol. 32, no. 4, pp. 973–987, 2016. View at: Publisher Site  Google Scholar
 N. Simaan, R. M. Yasin, and L. Wang, “Medical technologies and challenges of robot assisted minimally invasive intervention and diagnostics,” 2018, https://arxiv.org/abs/1807.03731. View at: Google Scholar
 Y. Chen, S. Xu, A. Squires et al., “MRIguided robotically assisted focal laser ablation of the prostate using canine cadavers,” IEEE Transactions on Biomedical Engineering, vol. 65, no. 7, pp. 1434–1442, 2017. View at: Google Scholar
 S. Lim, K. Sharma, P. Li et al., “Robotically assisted long bone biopsy under MRI: cadaver study results,” International Journal of ComputerAssisted Radiology and Surgery, vol. 14, no. 1, pp. 147–156, 2018. View at: Publisher Site  Google Scholar
 M. Giannakou, C. Yiallouras, G. Menikou, C. Ioannides, and C. Damianou, “MRIguided frameless biopsy robotic system with the inclusion of unfocused ultrasound transducer for brain cancer ablation,” The International Journal of Medical Robotics and Computer Assisted Surgery :MRCAS, vol. 15, p. e1951, 2019. View at: Publisher Site  Google Scholar
 B. Fallahi, C. Rossa, R. Sloboda, N. Usmani, and M. Tavakoli, “Partial estimation of needle tip orientation in generalized coordinates in ultrasound imageguided needle insertion,” in Proceedings of the 2016 IEEE International Conference on Advanced Intelligent Mechatronics (AIM), pp. 1604–1609, Banff, Canada, 2016. View at: Google Scholar
 Y. Kurose, Y. M. Baek, Y. Kamei et al., “Preliminary study of needle tracking in a microsurgical robotic system for automated operations,” in Proceedings of the 2013 13th International Conference on Control, Automation and Systems (ICCAS 2013), pp. 627–630, Gwangju, Republic of Korea, 2013. View at: Google Scholar
 A. Melzer, B. Gutmann, T. Remmele et al., “Innomotion for percutaneous imageguided interventions,” IEEE Engineering in Medicine and Biology Magazine, vol. 27, no. 3, pp. 66–73, 2008. View at: Publisher Site  Google Scholar
 H. J. Won, N. Kim, G. Kim, J. B. Seo, and H. Kim, “Validation of a CTguided intervention robot for biopsy and radiofrequency ablation: experimental study with an abdominal phantom,” Diagnostic and Interventional Radiology, vol. 23, pp. 233–237, 2017. View at: Google Scholar
 D. NavarroAlarcon, S. Singh, T. Zhang et al., “Developing a compact robotic needle driver for MRIguided breast biopsy in tight environments,” IEEE Robotics and Automation Letters, vol. 2, no. 3, pp. 1648–1655, 2017. View at: Publisher Site  Google Scholar
 S. Park, J.G. Kim, K.W. Lim et al., “A magnetic resonance imageguided breast needle intervention robot system: overview and design considerations,” International Journal of Computer Assisted Radiology and Surgery, vol. 12, pp. 1319–1331, 2017. View at: Google Scholar
 T. Zhang, D. NavarroAlarcon, K. W. Ng, M. K. Chow, Y. Liu, and H. L. Chung, “A novel palmshape breast deformation robot for MRIguided biopsy,” in Proceedings of the 2016 IEEE International Conference on Robotics and Biomimetics (ROBIO), pp. 527–532, Qingdao, China, December 2016. View at: Google Scholar
 Y. Chen, A. Squires, R. Seifabadi et al., “Robotic system for MRIguided focal laser ablation in the prostate,” IEEE/ASME Transactions on Mechatronics, vol. 22, no. 1, pp. 107–114, 2017. View at: Publisher Site  Google Scholar
 Z. Cao, T. Liao, W. Song, Z. Chen, and C. Li, “Detecting the shuttlecock for a badminton robot: a yolo based approach,” Expert Systems with Applications, vol. 164, Article ID 113833, 2020. View at: Publisher Site  Google Scholar
 W. Wang, A. N. Viswanathan, A. L. Damato et al., “Evaluation of an active magnetic resonance tracking system for interstitial brachytherapy,” Medical Physics, vol. 42, no. 12, p. 7114, 2015. View at: Publisher Site  Google Scholar
 Y. Chen, W. Wang, E. J. Schmidt et al., “Design and fabrication of MRtracked metallic stylet for gynecologic brachytherapy,” IEEE/ASME Transactions on Mechatronics, vol. 21, no. 2, pp. 956–962, 2016. View at: Publisher Site  Google Scholar
 J. Ladoiye, D. Necsulescu, and J. Sasiadek, “Force control of surgical robot with time delay using model predictive control,” in Proceedings of the 15th International Conference on Informatics in Control, Automation and Robotics, Porto, Portugal, 2018. View at: Google Scholar
 A. Squires, J. Oshinski, J. Lamanna, and Z. T. H. Tse, “Spinotemplate: a platform for MRIguided spinal cord injections,” Journal of Medical Robotics Research, vol. 1, no. 1, Article ID 1640006, 2016. View at: Publisher Site  Google Scholar
 L. Zhang, Z. Chen, W. Cui et al., “Wifibased indoor robot positioning using deep fuzzy forests,” IEEE Internet of Things Journal, 2020. View at: Publisher Site  Google Scholar
 W. Chen, T. Liao, Z. Li et al., “Using FTOC to track shuttlecock for the badminton robot,” Neurocomputing, vol. 334, pp. 182–196, 2019. View at: Publisher Site  Google Scholar
 Y. Koseki, D. De Lorenzo, K. Chinzei, and A. Okamura, “Coaxial needle insertion assistant for epidural puncture,” in Proceedings of the 2011 IEEE/RSJ International Conference on Intelligent Robots and Systems, San Francisco, CA, USA, 2011. View at: Google Scholar
 J. Kettenbach and G. Kronreif, “Robotic systems for percutaneous needleguided interventions,” Minimally Invasive Therapy & Allied Technologies, vol. 24, 2014. View at: Publisher Site  Google Scholar
 Y. Tipsuwan and M.Y. Chow, “Control methodologies in networked control systems,” Control Engineering Practice, vol. 11, no. 10, pp. 1099–1111, 2003. View at: Publisher Site  Google Scholar
 S. Zampieri, “Trends in networked control systems,” IFAC Proceedings Volumes (IFACPapersOnline), vol. 17, 2008. View at: Google Scholar
 X.M. Zhang, Q.L. Han, X. Ge et al., “Networked control systems: a survey of trends and techniques,” IEEE/CAA Journal of Automatica Sinica, vol. 7, no. 1, pp. 1–17, 2020. View at: Publisher Site  Google Scholar
 D. Uygun, “Lowpower and costeffective speed control card for remotelycontrolled dc motors,” Electronics World, vol. 121, pp. 16–19, 2015. View at: Google Scholar
 T. Savi and M. Radonji, “One solution for remote control through internet application,” in Proceedings of the 2016 Information Technologies Conference, Žabljak, Montenegro, 2016. View at: Google Scholar
 H. Su, C. Yang, G. Ferrigno, and E. De Momi, “Improved humanrobot collaborative control of redundant robot for teleoperated minimally invasive surgery,” IEEE Robotics and Automation Letters, vol. 4, no. 2, pp. 1447–1453, 2019. View at: Publisher Site  Google Scholar
 H. Su, Y. Schmirander, S. ValderramaHincapie et al., “Asymmetric bimanual control of dualarm serial manipulator for robotassisted minimally invasive surgeries,” Sensors and Materials, vol. 32, p. 1223, 2020. View at: Google Scholar
 C. Liu, W. H. Chen, and J. Andrews, “Model predictive control for autonomous helicopters with computational delay,” in Proceedings of the UKACC International Conference on Control, Coventry, UK, 2010. View at: Google Scholar
 P. Choomkasien, P. Chomphooyod, and D. Banjerdpongchai, “Design of model predictive control for industrial process with input time delay,” in Proceedings of the 2017 17th International Conference on Control, Automation and Systems (ICCAS), Jeju, Republic of Korea, October 2017. View at: Google Scholar
 N. Dong, Y. Feng, X. Han, and A. Wu, “An improved modelfree adaptive predictive control algorithm for nonlinear systems with large time delay,” in Proceedings of the 2018 IEEE 7th Data Driven Control and Learning Systems Conference (DDCLS), Enshi, China, May 2018. View at: Google Scholar
 H. Norizuki and Y. Uchimura, “Design of model predictive control considering time delay in remote excavation work,” in Proceedings of the 2018 IEEE 15th International Workshop on Advanced Motion Control (AMC), Tokyo, Japan, March 2018. View at: Google Scholar
 C. Cheng, S. Jiang, J. Liu, and J. Lou, “Design and kinematics analysis of a novel MRcompatible robot for needle insertion,” in Advances in Reconfigurable Mechanisms and Robots I, J. S. Dai, M. Zoppi, and X. Kong, Eds., pp. 443–453, Springer, London, UK, 2012. View at: Publisher Site  Google Scholar
 S. Jiang, F. Sun, W. Feng, L. F. Hofman, and Y. Yu, “Analysis of a novel highprecision 5degrees of freedom magnetic resonance imagingcompatible surgery robot for needleinsertion prostate brachytherapy,” Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, vol. 228, no. 5, pp. 865–876, 2014. View at: Publisher Site  Google Scholar
 E. G. Christoforou, I. Seimenis, E. Andreou, E. Eracleous, and N. V. Tsekos, “A novel, generalpurpose, MRcompatible, manually actuated robotic manipulation system for minimally invasive interventions under direct MRI guidance,” The International Journal of Medical Robotics and Computer Assisted Surgery, vol. 10, no. 1, pp. 22–34, 2014. View at: Publisher Site  Google Scholar
 S. Guo, Y. Wang, N. Xiao, Y. Zhao, Y. Zeng, and J. Wu, “A novel masterslave robotic system with close loop control for vascular interventional surgery,” in Proceedings of the 2018 IEEE International Conference on Mechatronics and Automation (ICMA), Tokyo, Japan, 2018. View at: Google Scholar
 J. Guo, Y. Liu, Q. Qiu et al., “A novel robotic guidance system with eye gaze tracking control for needle based interventions,” IEEE Transactions on Cognitive and Developmental Systems, p. 1, 2019. View at: Publisher Site  Google Scholar
 R. Kalman, “A new approach to linear filtering and prediction problems,” Journal of Basic Engineering, vol. 82D, pp. 35–45, 1959. View at: Google Scholar
 Y. Shang and T. Lu, “Optimal tracking control of mini unmanned helicopter based iterative approximating,” in Proceedings of the 2019 IEEE 3rd Information Technology, Networking, Electronic and Automation Control Conference (ITNEC), Chengdu, China, 2019. View at: Google Scholar
 R. Tanaka and T. Koga, “An approach to linear active disturbance rejection controller design with a linear quadratic regulator for a nonminimum phase system,” in Proceedings of the 2019 Chinese Control Conference (CCC), Gangzhou, China, 2019. View at: Google Scholar
 M. Dominici and R. Cortesäo, “Model predictive control architectures with force feedback for roboticassisted beating heart surgery,” in Proceedings of the 2014 IEEE International Conference on Robotics and Automation (ICRA), Hong Kong, China, 2014. View at: Google Scholar
 R. Abe, T. Satoh, N. Saito, J. Nagase, and N. Saga, “Disturbance observerbased model predictive control using timevarying constraints,” in Proceedings of the 2015 World Congress on Industrial Control Systems Security (WCICSS), London, UK, 2015. View at: Google Scholar
 J. Jerez, E. Kerrigan, and G. Constantinides, “A condensed and sparse QP formulation for predictive control,” in Proceedings of the 2011 50th IEEE Conference on Decision and Control and European Control Conference, Orlando, FL, USA, 2011. View at: Google Scholar
Copyright
Copyright © 2020 Jing Guo 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.