Surfing on an uncertain edge: Precision cutting of soft tissue
using torquebased medium classification
Abstract
Precision cutting of softtissue remains a challenging problem in robotics, due to the complex and unpredictable mechanical behaviour of tissue under manipulation. Here, we consider the challenge of cutting along the boundary between two soft mediums, a problem that is made extremely difficult due to visibility constraints, which means that the precise location of the cutting trajectory is typically unknown. This paper introduces a novel strategy to address this task, using a binary medium classifier trained using joint torque measurements, and a closed loop control law that relies on an error signal compactly encoded in the decision boundary of the classifier. We illustrate this on a grapefruit cutting task, successfully modulating a nominal trajectory fit using dynamic movement primitives to follow the boundary between grapefruit pulp and peel using torque based medium classification. Results show that this control strategy is successful in 72 % of attempts in contrast to control using a nominal trajectory, which only succeeds in 50 % of attempts.
I Introduction
The design of effective learning and adaptive control strategies for precision cutting remains an open problem in robotics [1, 2, 3]. This is particularly true for the case where cutting involves moving between two mediums, and where there is uncertainty in the location of these.
Cutting tasks of this form are regularly encountered in surgery, where tumour extraction is guided by well established continuum differences between tumours and normal tissue [4, 5]. This paper is motivated by wide local excision, a surgical procedure that aims to remove a tumour with a clear margin of healthy tissue around it. At present, manipulation tasks like these are inconceivable for autonomous robots, for a variety of reasons. First, the constrained operational space of nontrivial geometry restricts an endeffector’s maneuverability and is severely limited by visibility constraints. These visibility constraints are a particular challenge, and human surgeons often rely strongly on haptic feedback for cutting, using tactile tissue differences to guide procedures instead of vision. In addition, this task is highly variable and uncertain, due to the unpredictable behaviour of deformable tissue and varied tumour shapes. Finally, and most importantly, this contactrich task is characterised by safety constraints imposed on the region of operation and allowable applied forces. It is therefore critical to keep an endeffector inside a desired region while executing the task.
In this study, we move towards addressing the challenge of autonomous cutting with visibility constraints by 1) employing probabilistic inference to identify the boundary between two mediums using torque sensing, and 2) using the medium classification probability as an error signal for online, closedloop movement adaptation. As a testbed, we consider fruit processing, and study the task of scooping a grapefruit segment out of the membrane with a kitchen knife (see Fig. 1). This manipulation task shares several important characteristics with the surgery problem described above, including the complex geometry of the task space, the need for contactrich manipulation in a deformable environment and the presence of two mediums with differing material properties. Precision food processing is itself an industrially useful skill, and the ability to extract fruit portions without damaging food products is particularly valuable.
A key feature of the grapefruit testbed is an implicit task requirement to keep the knife inside the intermediate region between the peel and pulp boundary, such that most of the pulpy segment is extracted without the knife getting stuck in the peel, or too much grapefruit being left on the membrane. Since the exact shape and location of this boundary curve is unknown, it must be inferred during task execution. The key insight for the approach proposed in this work comes from observing grapefruit scooping when executed by human. It is clear that humans do not rely on an accurate geometrical model of the fruit, but instead apply a general scooping movement that is continued until it stops feeling right (i.e. when the knife starts progressively entering the peel). In these cases, the movement is either adjusted, or completely restarted using a different insertion angle. Our hypothesis is thus that cutting occurs using a rough nominal trajectory that is modulated by torque feedback resulting from the differing tissue characteristics of the two mediums being separated.
This work introduces a novel framework to accomplish tasks of this form. Here, we use Dynamic Movement Primitives (DMP) to learn a general scooping motion (the nominal trajectory) using kinesthetic demonstration. However, due to the variations in grapefruit’s shape and its mechanical properties, we show that generalisation of the learned movement primitive is inadequate. We therefore propose a control scheme in which the corrections to the DMP trajectory reflect the probability of knife being inserted into either of the two mediums. In this formulation, the point of highest uncertainty in this belief (probability of 0.5) serves as a proxy for the desired region of operation (i.e. the boundary between the pulp and peel). This probability is estimated at each time step of task execution by classifying torque readings from the joints of robot arm. We use a logistic regression classifier trained to disambiguate the mediums on a dataset of multiple task executions to demonstrate the feasibility of this method.
Ii Related Work
There is a substantial amount of research on the use of the force feedback in robotic manipulation tasks. However, most work is focused on rigid object manipulation (e.g. there is extensive research on the use of force data in the areas of robot door opening [6, 7, 8], grasping [7, 9, 10] and object identification [11, 12]), where task dynamics are relatively well understood and many mature techniques for motion planning and control are readily available. Unfortunately, many of these techniques are not applicable to the manipulation of deformable objects.
A particularly representative deformable object manipulation task involves food cutting, for example fruits or cheese. This process has timevarying nonlinear dynamics that are extremely difficult to describe analytically [13] (although attempts have been made [3, 14, 15]). As a result, learningbased techniques have been proposed to address this challenge. Lenz et. al [16] use deep learning techniques to model food cutting tasks and further use these in a modelpredictive control scheme. Here, robot controls are optimised in real time with respect to the constructed cost function, which penalises the height of the knife and its deviation from a sawing range, thus forcing a cutting movement. This approach was verified on a variety of food objects, such as lemons and potatoes, and showed its ability to adapt to both intraclass and timevarying variations in the physical properties of the objects.
A similar approach to learning the predictive model is described by Tian et. al [17], where the authors demonstrated tactile servoing using high dimensional tactile sensor data. Another use of learning in the latent space is presented in the work of Gemici and Saxena [18], which is concerned with robotic handling of food objects, e.g. grasping or piercing. Here, latent features of objects were learned from force data collected during the manipulation and then used to classify the objects for manipulation planning.
Many manipulation tasks, e.g. scooping, involve nontrivial kinematic trajectories that can be learned from demonstration. In [19], the authors propose a general framework (Dynamic Movement Primitives, or DMPs) for encoding complex movements as a parameterised policy. This framework, when coupled with feedback enables reactive movement adaptations [20, 9]. In our work, we use a predictive model of the expected sensor trace, that is based on the statistics of multiple task executions. This approach is similar to [21], where the statistics of sensor measurements from the successful task executions were used to construct a predictive model for online failure detection. We apply a similar method to model the region of operation (pulp or peel in our grapefruit example) using torque sensor readings. However, a key contribution of this work is the use of this classification scheme to perform boundary identification for cutting with visibility constraints, through the introduction of a control scheme for movement adaptation based on the estimated probability of being in a given medium.
Iii Problem Formulation
As discussed previously, this paper focuses on the task of precisioncutting between two mediums. Our primary interest lies in the development of an adaptive control framework, so we do not consider the use of any task specific cutting equipment or machinery, and focus on cutting using a standard kitchen knife. In addition, we allow control of the initial insertion position, and thus, we manually initialise the starting position of the knife.
The task described above can be formulated in a 2D task space. Consider two elastic mediums with different stiffness ( and ), separated by a curved boundary (see Fig. 2). Assume a strong prior over the stiffness and boundary curve (dashed line) is available, but no exact parameters are known. The objective is to steer the tip of the knife along the true boundary such that separation of mediums is maximised, while avoiding excessive deformations imposed to either of mediums by the knife. Since the exact curve of the boundary is unknown, the open loop execution of the prescribed path (based on a prior belief over the curve of the boundary) runs the risk of inserting the knife into the peel (in our grapefruit example), thereby severely restricting the knife’s maneuverability.
Iv Cutting using uncertainty feedback
We address the challenge above using a learning strategy where the desired operational region is compactly encoded in the decision boundary of a binary medium classifier. Here, the estimated likelihood of sensor readings associated with either medium guides the movement execution in the form of online trajectory correction. In summary, we 1) use the DMP framework to encode a nominal scooping trajectory, 2) learn probabilistic classification of sensor readings associated with operation in either mediums, and 3) construct a control scheme that corrects the DMP according to the estimated posterior distribution over either medium, as illustrated in Fig. 3.
Here, denotes the nominal joint trajectory, the corrected joint trajectory, and represents the nominal Cartesian trajectory. denotes the correction term applied to the Cartesian trajectory, is a gain matrix, refers to the sensed torque readings at time step , and is the probability of the knife being inserted into one of the mediums. FK and IK are Forward and Inverse Kinematics transformations for the robot arm. Note that gain matrix is time dependent, as the correction direction depends on the position along the nominal trajectory at time . We briefly discuss each element in the control framework below.
Iva Nominal trajectory modelling using DMPs
In the DMP formulation [19], any goaloriented movement primitive can be expressed as:
(1) 
where , and are the desired acceleration, velocity and position, respectively, is the goal position, is a temporal scaling factor, and are time constants, and is a nonlinear forcing function. In the above equation, the nonlinear term modulates the landscape of a global point attractor . Thus, an arbitrarily complex movement can be represented by appropriately constructing . Typically, the nonlinear function is represented using a normalized linear combination of basis functions:
(2) 
where is the number of basis functions with center , widths and weights .
Note that the forcing term does not depend on time, but does depend on phase variable that monotonically decays from 1 to 0 with a user specified rate :
(3) 
In the proposed framework, a nominal scooping trajectory is captured by kinesthetically guiding the robot arm and recording the endeffector’s Cartesian path (the time series of the endeffector’s position and orientation waypoints). The corresponding velocity and acceleration profiles of the movement, and , are obtained by twice differentiating the recorded endeffector’s path . The desired nonlinear function (expressed by rearranging (1)) is approximated by employing the Locally Weighted Regression method [22], used for optimising the weights of the basis functions. It should be noted that the nominal trajectory could be modelled using any behaviour cloning strategy, and the proposed approach is not limited to the use of DMPs.
In our approach, at each time step we add a local correction to the current point on the Cartesian path of nominal DMP (see Fig. 3). The correction term is given by
(4) 
where is timevarying positive definite gain matrix that defines the sensitivity of task variables and is the probability of the knife being in the medium at time step . Note, that the desired region of operation at each time step lies at the boundary between two mediums, where probability is equal to 0.5.
Thus, our proposed uncertainty driven control law can be formulated generally as
(5) 
where is the corrected version of the nominal trajectory .
IvB Logistic regression
In this work, we use logistic regression to model the probability of being in a given medium. In this approach, model parameters are fit by maximizing the probability of the data under a linear logistic model:
(6) 
where is the likelihood, is the number of training samples of torque readings, is the label (e.g. “Peel” or “Pulp”) of the th example of torque data, is a vector of torque readings of the th example and is a model parameter.
If the cost function is defined as the negative loglikelihood of labels , then the above expression is equivalent to minimizing:
(7) 
where is a sigmoid function and labels .
In order to discourage the optimizer from overfitting to the training data, the cost function can include an additional regularization term that penalizes extreme weight coefficients, e.g. , where denotes the regularization strength.
IvC Experimental setup
All experiments were conducted using a 7 degreeoffreedom PR2 robot arm. The PR2 arm is counterbalanced and highly compliant, and is wellsuited for kinesthetic demonstrations of flexible and fluid movements. The remaining elements of the experimental setup consisted of a chopping board clamped to the table, a halved grapefruit fixed to the chopping board, and a regular paring knife secured at the gripper (see Fig. 4). For registering the torques experienced at the joints of the arm we used standard PR2 joint effort readings (a joint torque estimate based on the joint motor current).
IvD Evaluation of nominal DMP
The learned scooping DMP was evaluated on the 10 randomly chosen segments. Before each trial, a segment was precut along the segment radii and the pose of the knife was manually adjusted, as discussed in the previous section. Successful task execution implies the complete extraction of the undamaged segment without the knife getting stuck in the peel.
The results of the trials agreed with original expectations, with only 2 successful task executions out of total 10. In 7 of the failed trials, the knife entered the peel and the execution was aborted. Moreover, an instance of tearing apart the segment during the scooping was registered. As anticipated, the main difficulty of the task was avoiding the knife’s insertion into the peel, where further knife maneuverability became limited.
V Learning the Boundary Region Using Sensed Torque
Va Dataset
The learned DMP was used to accumulate joint torque readings associated with successful (cut through the flesh) and failed task executions (cut into the peel). These traces of torque measurements were analysed and further used for training the logistic regression model to estimate the probability of the knife’s deviation from the desired region of operation, i.e. the boundary between pulp and peel, where task executions succeed. We used the learned DMP and experimental setup described in the previous section. The criteria for a successful trial remained unchanged from the preliminary evaluation of the DMP. A total of 111 scooping trials were conducted using a number of grapefruit, of which 55 trials were successful and 56 trials failed.
The nominal trajectory comprised 24 segments, at which a single snapshot of torque readings was taken. Thus the recorded data consisted of 24 timeindexed 7dimensional vectors. Fig. 5 shows the descriptive statistics of the collected data.
VB Classification
The dataset of 111 trials was randomized and split into 90 sensor traces allocated for training and validation and 21 traces held out for testing. The training and validation dataset consisted of 44 examples of “Pulp” torque traces and 46 examples of “Peel”. As discussed, each trace contained of 24 timeindexed samples of torque reading for each of the 7 joints. Thus, in total the training and validation dataset contained 1,056 and 1,104 individual examples of “Pulp” and “Peel” torque readings, respectively. The classifier’s input comprised of an 8dimensional vector (7 torque readings for each of the joints, plus the time index). The objective of the classification task was to estimate the probability of the measurement being taken inside of either medium given the current torque measurements. In our approach, a desirable property of a classifier is to be robust to the outliers and to handle the ambiguous inputs by reporting the appropriate levels of uncertainty (i.e. to avoid being overconfident). We used a logistic regression model, which we validated using the Kfold crossvalidation technique with 10 folds. Thus, each fold used 81 examples for training and 9 examples for validation. The validation and test results are given in Table I.
It should be noted that the modelling approach has an inherently noisy training process. Since all of the 24 torque samples in the sensor trace share the same label defined by the task outcome, all the intermediate phenomena are disregarded. For instance, if knife was closely following the desired boundary region throughout most of the execution but got stuck in the peel at the very end, all of the 24 torque readings would be labeled as “Peel”. However, the ability to capture the uncertainty demonstrated by Logistic Regression model alleviates this issue (see Fig. 6), and in some respects this training process forces a more conservative probabilistic model. It is important to note that despite misclassifying some intermediate samples, the trained model does not commit to any extreme beliefs over the mediums, unless the test input is strongly representative of a given class. Finally, in the case of ambiguous test inputs (i.e. where torque levels of the input trace appear uncharacteristic for a given label), the model demonstrates desirable levels of uncertainty, which is extremely important given the fact that we seek to use this for feedback control.
Vi Online Dmp Adaptation
Validation  Test  





Actual Peel  673  431  174  66  
Actual Pulp  306  750  65  199  
Sensitivity  0.61  0.73  
Specificity  0.71  0.75  

34%  26% 
We reused the nominal DMP and trained classifier from the previous section and constructed the closedloop control scheme as shown in Fig. 3. First, the nominal joint trajectory associated with learned DMP is transformed into the endeffector’s trajectory in the Cartesian space, where all the required corrections are relatively straightforward. As discussed previously, we use the probability of the knife being inserted into the peel for deriving the required motion corrections. We used a simple motion correction scheme for experimentation, in which the first half of the scooping motion (where the most dominant movement component involves pushing the knife downwards) is modulated towards the center of the grapefruit. In the second half of the movement (where the knife slides under the segment while moving towards the center of the grapefruit), the motion was modulated upwards. Thus, in each of the cases, the knife deviates from the peel region towards the pulp region, when the estimated probability of peel increases. For both cases we used a gain of 0.01m (i.e. 100 probability of peel would translate the movement 10mm away from the nominal trajectory in the prescribed direction). It should be noted that more complex schemes can be applied, e.g. modulation in the direction of the normal to the side of the knife.
Control scheme  Openloop  Closedloop 

Successful trials  55  36 
Failed trials  56  14 
Success rate  50%  72% 
For this experiment we conducted 50 trials of grapefruit scooping on randomly chosen segments in a total of 12 different grapefruits. As in the previous section, a successful trial required the complete extraction of an intact segment without the knife getting stuck inside the peel. The results are provided in the Table II. 36 out of 50 trials achieved successful task completion. In all of the 14 failed attempts, the knife entered the peel at the start of the cut and propagated deeply before the peel could be classified. In these cases, the movement corrections towards the center of the grapefruit failed, as the knife could not tear the peel with the side of the blade. Perhaps, in such cases the DMP can be reversed and reapplied with estimated corrections, as the classifier successfully reflected the event of knife being stuck in the peel.
In the successful trials, the knife visibly responded to the local increase in the resistance throughout the movement execution. It was clear that online movement adaptation improved the segment separation. Since the modulated motion acts in the direction approximately orthogonal to the boundary, it introduces a tearing effect. Similar tearing motions can be observed in human executed grapefruit scooping, where a knife’s reorientation lowers movement resistance by tearing through the fibers.
It should also be noted that the described task is strongly dependent on several factors. First, the nominal DMP plays an important role in the success of the task. Since the proposed method relies on torque readings gathered from the execution of the nominal trajectory, a poorly chosen movement can severely impair the medium classification. Second, the experiment is highly sensitive to the sharpness of the knife, as well as the position of the grapefruit relative to the initial pose of the knife. Nevertheless, these experiments highlight the promise of uncertainty driven cutting between mediums with differing stiffness properties.
Vii Conclusions and Future Work
We present an uncertainty driven feedback control law and demonstrated its performance on the task of grapefruit segmentation. This task is selected because it resembles a common surgical procedure where a hard tumour is extracted from soft tissue, and physical material properties used to guide human surgeons.
Automating tasks of this form is extremely challenging, as it requires cutting along an uncertain boundary, subject to visibility constraints.
Our experiments show that a simple movement correction scheme, where the movement of a robot arm is modulated along a single Cartesian axis in response to the probability of being in a given medium, significantly improves cutting performance. Future work involves the development of more complex movement adaptation schemes, and the extension of the approach to include higher level movement planning.
References
 [1] P. Long, W. Khalil, and P. Martinet, “Robotic deformable object cutting : From simulation to experimental validation,” in European Workshop on Deformable Object Manipulation (EWDOM)., 2014.
 [2] A. Yamaguchi and C. Atkeson, “Combining finger vision and optical tactile sensing: Reducing and handling errors while cutting vegetables,” in IEEERAS International Conference on Humanoid Robots. IEEE Computer Society, 2016, pp. 1045–1051.
 [3] X. Mu, Y. Xue, and Y. Jia, “Robotic cutting: Mechanics and control of knife motion,” in 2019 International Conference on Robotics and Automation (ICRA), May 2019, pp. 3066–3072.
 [4] V. Cristini and J. Lowengrub, Multiscale Modeling of Cancer: An Integrated Experimental and Mathematical Modeling Approach. Cambridge University Press, 2010.
 [5] A. L. McKnight, J. L. Kugel, P. J. Rossman, A. Manduca, L. C. Hartmann, and R. L. Ehman, “MR elastography of breast cancer: preliminary results,” American journal of roentgenology, vol. 178, no. 6, pp. 1411–1417, 2002.
 [6] A. Jain, H. Nguyen, M. Rath, J. Okerman, and C. C. Kemp, “The complex structure of simple devices: A survey of trajectories and forces that open doors and drawers,” in 2010 3rd IEEE RAS & EMBS International Conference on Biomedical Robotics and Biomechatronics. IEEE, 2010, pp. 184–190.
 [7] M. Kalakrishnan, L. Righetti, P. Pastor, and S. Schaal, “Learning force control policies for compliant manipulation,” in 2011 IEEE/RSJ International Conference on Intelligent Robots and Systems, Sep. 2011, pp. 4639–4644.
 [8] A. Jain and C. Kemp, “Improving robot manipulation with datadriven objectcentric models of everyday forces,” Autonomous Robots, vol. 35, no. 23, pp. 143–159, 2013.
 [9] P. Pastor, L. Righetti, M. Kalakrishnan, and S. Schaal, “Online movement adaptation based on previous sensor experiences,” in IEEE International Conference on Intelligent Robots and Systems, 2011, pp. 365–371.
 [10] J. M. Romano, K. Hsiao, G. Niemeyer, S. Chitta, and K. J. Kuchenbecker, “Humaninspired robotic grasp control with tactile sensing,” IEEE Transactions on Robotics, vol. 27, no. 6, pp. 1067–1079, Dec 2011.
 [11] A. Schneider, J. Sturm, C. Stachniss, M. Reisert, H. Burkhardt, and W. Burgard, “Object identification with tactile sensors using bagoffeatures,” in 2009 IEEE/RSJ International Conference on Intelligent Robots and Systems. IEEE, 2009, pp. 243–248.
 [12] J. A. Fishel and G. E. Loeb, “Bayesian exploration for intelligent identification of textures,” Frontiers in Neurorobotics, vol. 6, no. JUNE, 2012.
 [13] F. C. Moon and T. KalmárNagy, “Nonlinear models for complex dynamics in cutting materials,” Philosophical Transactions: Mathematical, Physical and Engineering Sciences, vol. 359, no. 1781, pp. 695–711, 2001.
 [14] G. Zeng and A. Hemami, “An adaptive control strategy for robotic cutting,” in Proceedings of International Conference on Robotics and Automation, vol. 1, April 1997, pp. 22–27 vol.1.
 [15] A. G. Atkins, X. Xu, and G. Jeronimidis, “Cutting, by ‘pressing and slicing,’ of thin floppy slices of materials illustrated by experiments on cheddar cheese and salami,” Journal of Materials Science, vol. 39, no. 8, pp. 2761–2766, Apr 2004.
 [16] I. Lenz, R. Knepper, and A. Saxena, “DeepMPC: Learning deep latent features for model predictive control,” in Robotics: Science and Systems XI, 2015.
 [17] S. Tian, F. Ebert, D. Jayaraman, M. Mudigonda, C. Finn, R. Calandra, and S. Levine, “Manipulation by feel: Touchbased control with deep predictive models,” in 2019 International Conference on Robotics and Automation (ICRA), 2019.
 [18] M. C. Gemici and A. Saxena, “Learning haptic representation for manipulating deformable food objects,” in 2014 IEEE/RSJ International Conference on Intelligent Robots and Systems, Sep. 2014, pp. 638–645.
 [19] A. J. Ijspeert, J. Nakanishi, H. Hoffmann, P. Pastor, and S. Schaal, “Dynamical movement primitives: Learning attractor models for motor behaviors,” Neural Computation, vol. 25, no. 2, pp. 328–373, Feb 2013.
 [20] P. Pastor, H. Hoffmann, T. Asfour, and S. Schaal, “Learning and generalization of motor skills by learning from demonstration,” in 2009 IEEE International Conference on Robotics and Automation, May 2009, pp. 763–768.
 [21] P. Pastor, M. Kalakrishnan, S. Chitta, E. Theodorou, and S. Schaal, “Skill learning and task outcome prediction for manipulation,” in 2011 IEEE International Conference on Robotics and Automation, May 2011, pp. 3828–3834.
 [22] S. Vijayakumar and S. Schaal, “Locally weighted projection regression: An o(n) algorithm for incremental real time learning in high dimensional space,” Proceedings of the Seventeenth International Conference on Machine Learning (ICML 2000), vol. Vol. 1, 05 2000.