iCLAP: Shape Recognition by Combining Proprioception and Touch Sensing

For humans, both the proprioception and touch sensing are highly utilized when performing haptic perception. However, most approaches in robotics use only either proprioceptive data or touch data in haptic object recognition. In this paper, we present a novel method named Iterative Closest Labeled Point (iCLAP) to link the kinesthetic cues and tactile patterns fundamentally and also introduce its extensions to recognize object shapes. In the training phase, the iCLAP first clusters the features of tactile readings into a codebook and assigns these features with distinct label numbers. A 4D point cloud of the object is then formed by taking the label numbers of the tactile features as an additional dimension to the 3D sensor positions; hence, the two sensing modalities are merged to achieve a synthesized perception of the touched object. Furthermore, we developed and validated hybrid fusion strategies, product based and weighted sum based, to combine decisions obtained from iCLAP and single sensing modalities. Extensive experimentation demonstrates a dramatic improvement of object recognition using the proposed methods and it shows great potential to enhance robot perception ability.

sensing are highly utilized when performing haptic perception. However, most approaches in robotics use only either proprioceptive data or touch data in haptic object recognition. In this paper, we present a novel method named Iterative Closest Labeled Point (iCLAP) to link the kinesthetic cues and tactile patterns fundamentally and also introduce its extensions to recognize object shapes. In the training phase, the iCLAP first clusters the features of tactile readings into a codebook and assigns these features with distinct label numbers. A 4D point cloud of the object is then formed by taking the label numbers of the tactile features as an additional dimension to the 3D sensor positions; hence, the two sensing modalities are merged to achieve a synthesized perception of the touched object. Furthermore, we developed and validated hybrid fusion strategies, product based and weighted sum based, to combine decisions obtained from iCLAP and single sensing modalities. Extensive experimentation demonstrates a dramatic improvement of object recognition using the proposed methods and it shows great potential to enhance robot perception ability.
Index Terms-Object recognition, robot systems, data fusion, tactile sensing.

I. INTRODUCTION
T HE sense of touch plays an important role in robot perception and many tactile sensors have been developed in the last few decades [1]- [3]. In addition to obvious applications of collision detection and avoidance, tactile sensors can be applied in multiple tasks such as object recognition [4]- [6], dexterous manipulation [7], [8] and localization [9], [10]. The haptic object recognition can be considered at two scales, i.e., local and global shapes [11], [12]. The local object shape, e.g., shapes that can fit into fingertips, can be recovered or recognized by a single touch, analogous to human cutaneous sense. The global shapes, e.g., contours that extend beyond the fingertip scale, usually require the contribution of both cutaneous and kinesthetic inputs. In such case, mechanoreceptors in joints are also utilized to acquire the movement of the fingers/end-effectors in space, with the assistance of local tactile features, to recognize the object identity. The kinesthetic inputs here are similar to human proprioception that refers to the awareness of positions and movements of body parts.
In this paper, an algorithm named Iterative Closest Labeled Point (iCLAP) and its extensions are presented to incorporate tactile and kinesthetic cues for haptic shape recognition. With only tactile readings, a dictionary of tactile features can be first formed. By searching for the nearest word in the dictionary, each tactile feature is assigned a label number. A four-dimensional data point is then obtained by concatenating the 3D position of the object-tactile sensor interaction point and the word label number. In this manner, 4D point clouds of objects can be obtained from training and test data. The partial 4D point cloud obtained from test set iteratively matches with all the reference point clouds in the training set and the identity of the best-fit reference model is assigned to the test object. Furthermore, weighted sum based and product based fusion strategies have also been developed for haptic shape recognition.
The contributions of this paper are: • A novel iCLAP algorithm that incorporates tactile readings and kinesthetic cues for object shape recognition is proposed; • Extensions of iCLAP algorithm based on different fusion strategies are created to enhance the recognition performance; • Extensive experiments are conducted to demonstrate the difference in the recognition performances between different approaches. This paper extends our previous work [13] by including extension methods of the iCLAP algorithm and introducing more thorough experiments. The remainder of this paper is organized as follows. The literature in the haptic object shape processing is reviewed in Section II. The proposed tactilekinesthetic shape recognition system is described in Section III. A series of the fusion approaches, both product-based and weighted sum based, are then introduced in Section IV. In Section V, the experimental setup and the data acquisition are described. The experiment results are then provided and analyzed in Section VI. Finally, in Section VII the conclusions are drawn and possible applications and future research directions are presented.

II. LITERATURE REVIEW
Thanks to the development of tactile sensor technologies [14]- [16], haptic object shape processing has received increasing attention in the recent years [17]. Some research arXiv:1806.07507v1 [cs.RO] 19 Jun 2018 produced point clouds by collecting contact points to constrain the geometry of the object [18], [19] whereas some others relied on tactile appearances by extracting features from tactile readings [5], [20], [21]. In contrast, for humans, the sense of touch consists of both kinesthetic (position and movement) and cutaneous (tactile) sensing [11]. Therefore, the fusion of spatial information of sensor movements and features of tactile appearance could be beneficial for object recognition tasks.
The methods based on contact points often employ techniques from computer graphics to fit the obtained cloud of contact points to a geometric model and outline the object contour. These methods were widely used by early researchers due to the low resolution of tactile sensors and prevalence of single-contact force sensors [19], [22], [23]. The resultant points from tactile readings can be fit to either super-quadric surfaces [19] or a a polyhedral model [18] in order to reconstruct unknown shapes. Different from the point cloud based approaches, a non-linear model-based inversion is proposed in [24] to recover surface curvatures by using a cylindrical tactile sensor. In more recent works [25]- [27], the curvatures at curve intersection points are analyzed and thus a patch is described through polynomial fitting. Tactile array sensors have also been utilized to obtain the spatial distribution of the object in space. In [10] an object representation is constructed based on mosaics of tactile measurements, in which the objects are a set of raised letter shapes. Kalman filters are applied in [28] to generate 3D representations of objects from contact point clouds collected by tactile array sensors, and the objects are then classified with the Iterative Closest Point (ICP) algorithm. Through utilizing these methods, arbitrary contact shapes can be retrieved, however, it can be time consuming when exploring a large object surface as excessive contacts are required to recognize the global object shape.
Another approach is to recognize the contact shapes by extracting shape features from pressure distributions within tactile images. The image descriptors from computer vision have been applied to represent the local contact patterns, e.g., image moments [5], [29], [30], SIFT based features [5], [21] and raw readings [20], [31]. However, there is only a limited number of approaches available for recovering the global object shape by analyzing pressure distributions in tactile images collected at different contact locations. One popular method is to generate a codebook of tactile features and use it to classify objects [20], [21], [32], [33]; a particular paradigm is the bag-of-words model. In this framework, only local contact features are taken to generate a fixed length feature occurrence vector to represent the object whereas the three-dimensional distribution information is not incorporated.
For humans, haptic perception makes use of both kinesthetic and cutaneous sensing [11], which can also be beneficial to robots. In [34], [35], a series of local "tactile probe images" is assembled and concatenated together to obtain a "global tactile image" using 2D correlation techniques with the assistance of kinesthetic data. However, in this work tactile features are not extracted whereas raw tactile readings are utilized instead, which would bring high computational cost when investigating large object surfaces. In [36], the tensor product operation is applied to code proprioceptive and cutaneous information, using Self-Organizing Maps (SOMs) and Neural Networks. In [37], the tactile and kinesthetic data are integrated by decision fusion and description fusion methods. In the former, classification is done with two sensing modalities independently and recognition results are combined into one decision afterwards. In the latter, the descriptors of kinesthetic data (finger configurations/positions) and tactile features for a palpation sequence are concatenated into one vector for classification. In other words, the information of the positions where specific tactile features are collected is lost. In both methods, the tactile and kinesthetic information is not fundamentally linked. In a similar manner, in [38] the tactile and kinesthetic modalities are fused in a decision fusion fashion. Both tactile features and joint configurations are clustered by SOMs and classified by ANNs separately and the classification results are merged to achieve a final decision. In a more recent work [39], the actuator positions of robot fingers and tactile sensor values form the feature space to classify object classes using random forests but there are no exploratory motions involved, with data acquired during a single and unplanned grasp.
In our paper, a novel iCLAP algorithm is proposed to incorporate tactile feature information into the location description of the sensor-object contact in a four dimensional space. In this way, the kinesthetic cues and tactile patterns are fundamentally linked. The experiments of classifying 20 real objects show that the classification performance is improved by up to 14.76% by using iCLAP compared to methods based on single sensing sources and a high recognition rate up to 95.52% can be achieved when a hybrid fusion is applied.

III. ITERATIVE CLOSEST LABELED POINT
In the haptic exploration, at each data collection point, both the tactile reading and the 3D sensor location can be recorded simultaneously. An illustration of data extracted from a pair of scissors is depicted in Fig. 1. Our proposed Iterative Closest Labeled Point (iCLAP) algorithm is based on two sensing sources, i.e., appearance features from obtained tactile images and spatial distributions of objects in space.
To create a label for each tactile feature, a dictionary of tactile features is first formed from training tactile readings and each tactile feature is then assigned a label by "indexing" the dictionary. The dictionary is formed by clustering the features extracted from tactile images using k-means clustering, where k is the dictionary size, namely, the number of clusters. The features of training/test objects are then assigned to their nearest clusters in Euclidean distance and are therefore labeled with the cluster numbers µ ∈ {1, 2,. . . , k}.
With the feature labels created from tactile readings and the sensor locations in 3D space, a single 4D point p can be represented by a tuple p = (p x , p y , p z , p µ ), where p x , p y , p z are the x, y, z coordinates of the tactile sensor in 3D space and p µ is the word label assigned to this location respectively. In this manner, the object can be represented in a four-dimensional space. To calculate the mutual distance between 4D sparse data point cloud P in the test set and the model (reference) point clouds Q in the training set, the Iterative Closest Point (ICP) algorithm [40] is extended to To find the closest point in the model point cloud to each transformed test data point p i , a k-d tree of the model point cloud is constructed [41], [42]. An error metric E iCLAP is defined to evaluate the mean square root distance of the associated point pairs and it is minimized with the optimal rigid rotation matrix R and translation vector t. In physical terms, the error metric E iCLAP can be visualized as the total potential energy of springs attached between matched point pairs. E iCLAP of point clouds P and Q to be matched is The centroids of P and Q are c p = 1 n n i=1 p i and c q = 1 m m i=1 q i , where n and m are the number of test and model data points respectively. In general case, n is equal to the number of the matched point pairs N. The point deviations from the centroids of P and Q can be obtained as To minimize the error metric, t is set as t = c q − Rc p . Therefore, the error metric is simplified as where R is orthogonal for orthogonal transformation, therefore, In an expanded form, where M ab = N i=1 p ia q ib , and a, b ∈ {x, y, z, µ}. To minimize E iCLAP the trace tr(RH) has to be maximized. Let the columns of H and the rows of R be h j and r j respectively, where j ∈ {1, 2, 3, 4}. The trace of RH can be expanded as where the inequality is just a reformulation of the Cauchy-Schwarz inequality. Since the rotation matrix R is orthogonal, its row vectors all have unit length. This implies where the square root is taken in the operator sense and tr( which is the maximum according to Eq. 5. It means that E iCLAP is minimized with the resulting optimal rotation matrix R = V U T and translation vector t = c q − Rc p . The iCLAP is iterated until any of the termination conditions is reached: error metric E iCLAP > preset tolerance; number of iterations > preset maximum number of iterations n max ; the relative change in the error metrics of two consecutive iterations falls below a predefined threshold. The obtained distances between the test and reference point clouds are then normalized by L2 norm. A reference point cloud with the minimum E iCLAP can be found and its identity is assigned to the test object by comparing E iCLAP .

IV. FUSION METHODS FOR HAPTIC OBJECT RECOGNITION
In general, the fusion of multiple modalities can provide complementary knowledge and improve the decision-making performance; it can be performed at either the feature level or the decision level [43], [44]. In feature fusion, the features extracted from different modalities are combined into high dimensional feature vectors prior to the classification, which are fed into a single classifier. In essence, our iCLAP algorithm is performed at the feature fusion level. However, due to the distinct representations of information sources, i.e., sensor locations and tactile labels in our case, normalization is hard to be performed and it needs to find the best normalization parameters by trial and error. Instead, decision fusion combines the decisions made based on individual modalities and makes a final decision in the semantic space where individual decisions usually have the same representation. Therefore, the decision level strategy is also adopted. Two different decision fusion methods are developed, i.e., weighted sum-based and productbased. To take advantage of both feature and decision fusion methods, hybrid fusion approaches are developed to combine the decisions of iCLAP and methods based on single sensing modalities. In total, nine synthesis methods, including one at feature level (iCLAP), two at decision level and six at hybrid level, are developed. There are two recognition pipelines, i.e., tactile based and kinesthetics based.
1) Tactile based object recognition: As described in Section III, a dictionary is formed by clustering tactile features extracted from the training tactile images using k-means clustering. The descriptors are then assigned to their nearest clusters in Euclidean distance. In this way, both training and test objects can be represented by histograms of word occurrences h class and h test respectively. Therefore, the distance E BoW between the test and reference objects can be computed using the histogram intersection metric of their word occurrence histograms h test and h class [21]: where k is the dictionary size.
2) Kinesthetics based object recognition: In each exploratory procedure, the locations of contact points form a point cloud. To calculate the mutual distance between 3D sparse data point cloud P 3 in the test set and the reference point clouds Q 3 in the training set, the classic Iterative Closest Point (ICP) [42] algorithm is employed. Similar to the proposed iCLAP algorithm, a point-to-point error metric E ICP between P 3 and Q 3 is defined: where R p and t p denote the rotation matrix and translation vector in 3D space, respectively, and N p is the number of matched point pairs. It is minimized iteratively to achieve an optimal transformation from the source points in P 3 to corresponding target points in Q 3 . The obtained distances between the test point cloud and the reference models in the training set are then normalized by L2 norm.

A. Decision fusion methods
Two decision fusion methods are proposed to synthesize both the tactile-based and kinesthetic-based recognition results. One is to calculate the distance between the test object and reference objects by a weighted sum of the two distances obtained from the two pipelines I B E + , which is named as the weighted sum fusion. In this method, the distances obtained from tactile sensing and sensor movements are combined in a linear fashion. It combines the scores/decisions in two modalities and ranks the combined results. Here I and B denote ICP and BoW respectively. Let w p1 be the weight assigned to the kinesthetic sensing source, I B E + can be formed as: The other is to acquire the distance by the product of E ICP and E BoW , which is named as product fusion. This method is based on the probability analysis of the likelihoods of the object classes based on the two modalities. The distance I B E × can be formed as: The identity of the reference object with the nearest distance to the test object is therefore assigned to the test object.

B. Hybrid fusion methods
To exploit the advantages of both feature and decision fusion strategies, a hybrid fusion strategy can be applied. Accordingly, the recognition results of iCLAP algorithm (feature fusion) are further integrated with decisions made by methods using single sensing modalities by a decision fusion to obtain the final decisions. In total, six hybrid fusion methods are developed, combining different decision fusion manners and various number of sensing modalities used. It can be divided into two groups with regards to the decision fusion manners. The first group is to compute the distances in a weighted sum manner of the two distances obtained from the kinesthetic only based pipeline ( I+ I E + ), tactile only based pipeline ( I+ B E + ) or both of them ( All E + ) and iCLAP algorithm. Here I+ stands for iCLAP algorithm. Let w p2 be the weight assigned to the kinesthetic sensing source, I+ I E + can be formed as: Let w w3 be the weight assigned to the tactile sensing source, I+ B E + can be formed as: Let w p4 and w w4 be the weights assigned to the kinesthetic and tactile sensing sources respectively, All E + can be formed as: The second group is to evaluate the distance between the test object and the reference objects in a product manner of the distances obtained from the kinesthetic only based pipeline ( I+ I E × ), tactile only based pipeline ( I+ B E × ) or both of them ( All E × ) and iCLAP algorithm. I+ I E × , I+ B E × , All E × can be formed respectively as: Fig. 2: A Phantom Omni device is used as a robotic arm to explore the object and a tactile sensor is attached to its stylus.

V. DATA COLLECTION
As illustrated in Fig. ??, the experimental setup consists of a tactile sensor and a positioning device. A Weiss tactile sensor WTS 0614-34A 1 was used that has an array of 14×6 sensing elements. It is of 51mm × 24mm and has a spatial resolution of 3.4mm × 3.4mm for each element. The sensing array is covered by elastic rubber foam to conduct the exerted pressures. A rate of 5 fps was used as we used in our initial studies [21], [32]. A Phantom Omni device 2 with six degrees of freedom was used for positioning the tactile sensor and its stylus served as a robotic manipulator. The stylus endeffector position can be obtained that has a nominal position resolution of around 0.055mm. The tactile sensor is attached to the stylus and its surface center is aligned with the end point of the stylus, therefore the end-effector position can be taken as tactile sensor position in the world coordinates.
During the data collection, each object was explored five times. Each exploration procedure was initialized without object-sensor interaction. The stylus was controlled with a speed of around 5 mm/s to explore the object while keeping sensor plane normal to the object surface; in this manner, the object surface was covered while a number of tactile observations and movement data of the tactile sensor were collected. Following [20], [21], an uninformed exploration strategy was employed. In total, 8492 tactile images with corresponding contact locations for 20 objects were collected, as shown in Fig. ??. It can be found that some objects are of similar appearances or spatial distributions. For example, pliers 1 and 2 are of similar size and have a similar frame, whereas they have different local appearances, i.e., the shape of jaws. On the other hand, some objects have similar local appearances but have different spatial distributions, for instance, fixed wrenches 1 and 2.

VI. RESULTS AND ANALYSIS
To evaluate the performance of the proposed iCLAP algorithm and extended methods, they are utilized to classify the 20 objects in the experiments. A leave-one-out cross validation is taken and averages of cross validation results are used. The general objective is to achieve a high recognition rate while minimizing the amount of samples needed. Following [13], the dictionary size k is set to 50 through the experiments. In [13], we also compared different features and Tactile-SIFT features [21] were found to perform best that are therefore also used in this paper.

A. iCLAP vs methods using single sensing modalities
The classification results by applying BoW (tactile only), ICP (kinesthetics only) and iCLAP with different number of object-sensor contacts, from 1 to 20, are shown in Fig. 4. As the number of contacts increases, all the performances of three approaches are enhanced. When the tactile sensor contacts the test object for less than 3 touches, the tactile sensing can achieve a better performance than the kinesthetic cues as tactile images are more likely to capture key appearance features. In addition, iCLAP outperforms the ICP by up 14.76%, while performing similarly to BoW. As the number of contacts increases, the performance of iCLAP improves dramatically and it performs much better than those with only one modality. It means that our proposed iCLAP algorithm exploits the benefits of both tactile and kinesthetic sensing channels and achieves a better perception of interacted objects. When the number of contacts is greater than 12, the recognition rates of all the three methods grow slightly and iCLAP still outperforms the other two with single sensing modalities. With 20 touches, an average recognition rate of 80.28% can be achieved by iCLAP.

B. iCLAP vs decision fusion approaches
Both product and weighted sum based decision fusion methods combine recognition decisions of BoW and ICP. The weight w p1 assigned in the weighted sum-based method (Eq. 9) has been investigated to find the optimal combination by brutal force search. The recognition rates with weights w p1 from 0.1 to 0.9 at an interval of 0.1 are shown in Fig. 5. It is found that a good recognition rate of around 90% can be achieved with 15 touches if w p1 is in the range of from 0.5 to 0.8. And it is observed that the best recognition performance can be achieved with w p1 =0.7, i.e., As shown in Fig. 6, the weighted sum-based decision fusion method surpasses iCLAP consistently; the product-based decision fusion approach falls behind iCLAP when limited touches are obtained (<5) whereas outperforms it when more touches are acquired. Moreover, the weighted sum-based method leads at first (<15) yet is caught up by product-based approach when more data are gathered. A good recognition performance, around 90%, can be achieved by applying either decision fusion approach. It can be estimated that a hybrid fusion strategy, combining decisions of iCLAP and separate sensing sources, can further enhance the classification accuracy.

C. Hybrid fusion of iCLAP and classic ICP
The classification results of iCLAP are first fused with the decisions of classic ICP using only kinesthetic cues to achieve    a hybrid conclusion of the object identity. The weight w p2 assigned in the weighted sum-based fusion method (Eq. 11) has been studied using brutal force search, set from 0.1 to 0.9 at an interval of 0.1. The changes of the recognition rate with weight w p2 are shown in Fig. ??. It is found that a recognition rate of around 80% can be achieved with 15 touches if w p2 is in the range of from 0.1 to 0.4. And it is observed that the best recognition rate can be achieved with w p2 =0.1, i.e., The recognition rates are shown in Fig. ??, against the number of contacts. It can be noticed that iCLAP consistently outperforms both the hybrid fusion methods, i.e., product and weighted sum-based approaches. The probable reason is that the inaccurate matching of spatial distributions deteriorates the overall hybrid fusion performance. Nevertheless, the trend of performance enhancement follows that of decision fusion methods, i.e., the weighted sum-based hybrid method outperforms the product-based hybrid approach when limited (<5) touches are available whereas the gap is narrowed when more inputs are supplemented.

D. Hybrid fusion of iCLAP and BoW classification pipeline
The decisions obtained from iCLAP are also combined with those achieved by BoW framework using only local pressure patterns, in a hybrid fusion manner. The weight w w3 assigned in the weighted sum-based method (Eq. 12) has been investigated by brutal force search, set from 0.1 to 0.9 at an interval of 0.1. The changes of the recognition rate with regards to the variance of weight w w3 are illustrated in Fig. 9. It is found that a good recognition rate of around 90% can be achieved with 15 touches if w w3 is in the range of from 0.1 to 0.5. And it is observed that the best recognition performance can be achieved when w w3 =0.2, i.e., In Fig. 10, the recognition rates of iCLAP and hybrid methods, both product-based and weighted sum-based, are compared  against the number of touches. It can be observed that both hybrid approaches achieve better recognition performance than original iCLAP algorithm consistently, which means the inclusion of the decisions made by BoW framework brings benefit to enhancing the recognition performance of iCLAP algorithm. Similar to the decision fusion methods, the weighted sumbased hybrid pattern out-stands to a large extent when limited (<5) touches are available yet is caught up by product-based hybrid approach when more data are collected. In addition, it can also be found that a satisfactory recognition rate, appropriately 90%, could be reached with only 10 sensorobject contacts with the hybrid fusion strategy.

E. Hybrid fusion of iCLAP, ICP and BoW
The recognition results based on iCLAP are then merged with both decisions acquired by BoW and ICP. As the same, product-based and weighted sum-based hybrid fusion strategies are employed. By employing brutal force search, the weight combinations of w p4 and w w4 assigned in the weighted sum-based method (Eq. 13) has been studied, set from 0.1 to 0.8 at an interval of 0.1 with the sum of the weights of three parts as 1. The recognition rates with various weights w p4 and w w4 are shown in Fig. 11. It can be seen that the series with w w4 =0.2 perform better than ones with the other w w4 values. The combination of w p4 =0.2 and w w4 =0.2 was found to have the best recognition performance, i.e., As shown in Fig. 12, the performances of both fusion strategies have been improved compared to the hybrid fusion approaches with only ICP classification pipeline. The probable reason is that the complement of tactile sensing pipeline makes up the inaccurate matching caused by kinesthetic sensing channel. However, the situation alters when compared to the hybrid fusion approaches with decisions of only BoW framework. The recognition rates of weighted sum-based hybrid fusion method are enhanced whereas the performance of productbased hybrid fusion method is deteriorated. The possible reason is that the inclusion of inaccurate kinesthetic based recognition undermines the tactile sensing based classification performance, especially when employing the product-based strategy. A decent recognition rate can be reached with only 10 touches; a comprehensive highest rate of 95.52% (averaged by cross validation results) can be achieved with 20 sensor-object contacts when a weighted sum based hybrid fusion strategy is employed to integrate decisions of iCLAP, BoW and ICP.

F. Discussions
Multi-dimensional ICP. The obvious unit mismatch might have a strong effect on the recognition performance using an integer index directly as the 4th dimension of points in the ICP algorithm, e.g., features of clusters 1 and 2 may be more distinct than features of clusters 1 and 7. It is an open issue how to normalize the 4D points due to the different nature of geometrical point coordinates and tactile feature labels. It is proposed that the cluster labels could be ranked by the mutual similarity of different clusters; both the coordinate values and the feature labels could be normalized into the range of from 0 to 1 and weighted during the closest point search. As an alternative way, each data point can be represented as the concatenation of its three positional coordinates with the feature descriptor obtained in this location to form a point in the multi-dimensional space. This method introduces more information of local pressure patterns in the multi-dimensional description of the objects and is potential to improve the algorithm robustness and the recognition performance but this method will introduce more computational burden.
Optimal weight assignment in fusion methods. Based on the experiment results using different fusion methods, it can be found that the weighted sum fusion can achieve a better recognition performance than the product based fusion. However, on the other hand, the complexity of finding an optimal weight assignment is the major drawback of the weighted sum fusion. The brutal force approach is straightforward and easy to be implemented but it brings about additional computational cost and the selection of weights needs to be done by trial and error. The issue of search for appropriate weights for different modalities is still an open question as discussed in [43] and beyond the scope of this paper.
Haptic exploration. In the autonomous haptic exploration with the multi-fingered hands with skins, there will have multiple contacts on the object when the object is held inhand. With the collected tactile patterns and contact locations, the algorithm proposed in this paper can be implemented for the object recognition task. To explore the unknown object, the planning methods in the literature [5], [45] can be applied to design the exploration strategies. In the haptic exploration, a confidence level of object recognition can be defined, and multiple grasps are implemented to explore the object surfaces iteratively until the defined confidence level is achieved.

VII. CONCLUSION AND FUTURE DIRECTIONS
In this paper we propose a novel algorithm named Iterative Closest Labeled Point to integrate tactile patterns with kinesthetic cues applied in the object recognition task. A bagof-words framework is first applied to exploit the collected training tactile readings and a dictionary of tactile features is formed by clustering tactile features. By indexing the dictionary, each tactile feature is assigned a word label. The numerical label is then appended to the 3D location of tactile sensor where it is obtained, to compose a 4D descriptive point cloud of the objects. During the object recognition, the test object point cloud is transformed iteratively to find the optimal reference model in the training set and the distance between them serves as the similarity criterion. The experimental results of classifying 20 objects show that iCLAP can improve the recognition performance (up to 14.76%) compared to the methods using only one sensing channel.
We also extend iCLAP to a series of approaches by employing hybrid fusion strategies. The recognition performance has been further improved when the decisions of iCLAP and BoW framework (or also with classic ICP) are integrated using either product based or weighted sum based fusion strategy. Besides, it can be observed that the weighted sum based fusion method outstands when limited number of contacts is available whereas both strategies perform quite well when considerable readings are collected (greater than 10 touches). A satisfactory recognition rate of 95.52% can be achieved when 20 touches are used and a weighted sum based hybrid fusion strategy is employed to integrate decisions of iCLAP, BoW and classic ICP.
The proposed iCLAP and its extended approaches can be applied to several other fields such as computer vision related applications. In the view of scene classification [46], as the landscape observations are correlated with the locations where they are collected, the proposed iCLAP combining the two sensing modalities is expected to enhance the classification performance. It can also be applied to medical applications such as Minimally Invasive Surgeries (MIS). If the local tactile patterns or visual appearances are merged with the spatial distributions of these observations by iCLAP, it is estimated to better recognize the interacted workspace within the body.
There are several directions to extend our work. As only the word label is utilized in the iCLAP to represent the tactile data, there is information loss to certain extent. Therefore, it is planned to include more clues of tactile patterns in the future designed algorithm. In addition, it will also be studied to recognize objects with multiple tactile sensing pads. It is proposed to implement our proposed algorithms on an instrumented robotic hand with multiple tactile array sensors on the fingers and the palm. And it will be explored how to minimize the number of touches to recognize objects in such cases.