# Image Classification Using PSO-SVM and an RGB-D Sensor

Authors: Carlos López-Franco; Luis Villavicencio; Nancy Arana-Daniel; Alma Y. Alanis
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

Image classification is a process that depends on the descriptor used to represent an object. To create such descriptors we use object models with rich information of the distribution of points. The object model stage is improved with an optimization process by spreading the point that conforms the mesh. In this paper, particle swarm optimization (PSO) is used to improve the model generation, while for the classification problem a support vector machine (SVM) is used. In order to measure the performance of the proposed method a group of objects from a public RGB-D object data set has been used. Experimental results show that our approach improves the distribution on the feature space of the model, which allows to reduce the number of support vectors obtained in the training process.


---

## Body

## 1. Introduction

Over the past years, there has been an increasing interest in object recognition. Object recognition can be divided in two major tasks: object localization and image classification.Object localization detects instances of a given category in the image.Image classification can be defined as the task of defining labels to an image, depending on the presence of an object.In this paper we propose an image classification system based on invariant moment descriptor that includes depth information. The 3D data allows producing small and robust descriptors that will improve the image classification. These descriptors are constructed using object models with rich information of the distribution of points. The model generation stage requires that best points are selected; therefore this stage can be defined as an optimization problem.Mathematical optimization is the selection of the best element with regard to some criteria. In the simplest case, it is consisted of maximizing or minimizing a fitness function [1]. Metaheuristic designates a computational method that optimizes a problem by iteratively trying to improve a candidate solution [2]. Many metaheuristics implement nature-inspired stochastic optimization. One of these algorithms is particle swarm optimization (PSO) developed by [3] and inspired by social behavior of bird flocking or fish schooling. It has been applied in many fields such as tremor analysis for biomedical engineering, trajectory planning [4], electric power [5], and image processing [6]. Optimization algorithms are often used in computer vision tasks such as image classification, which finds a relation between an input image and a set of previously known models [7].The sensor used in this work is a Kinect [8] which is an RGB-D sensor providing synchronized color and depth images. This sensor is widely used by the computer vision community due to its capabilities.In this work we analyze the inclusion of 3D information (provided by an RGB-D sensor) and the use of the PSO algorithm to create robust object models. Using this approach we can construct small and robust descriptors that improve image classification.The rest of the paper is organized as follows. The next section will present the proposed image classification system. In Section4 we present the mesh optimization process, with a brief introduction to the PSO algorithm. In Section 5 the invariant moment descriptor and classification are explained. In Section 6 we show the results of the proposed approach. Finally, in Section 8 we give the conclusions.


## 2. Related Work

In recent years image classification has become a very active research field. The goal of the classification task is to develop an algorithm which will assign a class-label to each one of the samples in the training data set. There are two main approaches for classification in the literature, namely, (a) supervised approach and (b) unsupervised approach, where the former uses a set of samples to train the classification, and the latter performs the classification by exploring the data [9]. Many popular image classification methods [10, 11] are based on local descriptor and support vector machines as basic techniques.One of the most widely used detectors is the Harris Corner detector [12]. However the Harris detector is not invariant to affine transformations and scale [13]. In [14], an automatic scale detection was proposed. Later in [15], the authors present an algorithm for interest point detection which is invariant to important scale changes. Then in [16], the SIFT detector is proposed; the approach approximates the Laplacian of Gaussians (LoG) by difference of Gaussians (DoG) filter to identify the potential keypoints. Later in [17] the SURF detector was presented, and SURF uses an approach similar to SIFT; however the approach uses integral images instead of DoGs; this allows a fast computation of approximate LoG using a box filter.Many different feature descriptors have been proposed in the literature: steerable filters [18], moment invariants [19], phase-based local features [20], Gaussian derivatives [21], and descriptors that represent the distribution of smaller-scale features [10]. The SIFT detector produces robust descriptors; however the speed of the detector is not very high. In [17], the SURF detector was presented; this approach is faster than SIFT and it has a similar performance on real images [22].Histogram of oriented gradients (HoG) was proposed in [23]. The HoG local patters are based on gradient histogram. The HoG and SIFT detectors share the same concept of implementation in the sense that local histograms of oriented patches defined across the image are used to produce the descriptor. In [23] the authors used HoG descriptors as feature vectors for a linear SVM. In [24] the authors propose the use of HoG for object detection and show that HoG descriptors have robust and effective features for such application. In [25] the authors build an efficient detector using AdaBoost to train a chain of progressively more complex region rejection rules based on Haar descriptors for pedestrian detection.Our Approach. In the last decades, SVMs have been proven to be an effective classifier approach for small sample sets or high dimensional problems [26, 27]. Such result is very important since both problems are too difficult to be solved by classical paradigms.One of the problems encountered during the design of an image classification technique is data overfitting, which arises when the training samples are small in comparison with the number of features. To overcome this problem we design a descriptor with an optimized 3D point distribution.The contribution of this work is the development of a small and robust descriptor based on invariant moments and 3D information in order to improve the classification process. The 3D information is incorporated from depth data obtained from an RGB-D sensor. The object model is optimized through the use of a PSO algorithm; this optimization allows improving the classification.


## 3. Image Classification System

The proposed classification system uses 3D information and is based on local features, invariant moments, contour generation, and a reduction of depth information using a mesh grid. It consists of five main steps: feature extraction, contour creation, mesh reduction, mesh optimization, and invariant moment descriptor formation; see Figure1.Figure 1
Image classification system.


### 3.1. Feature Extraction

The first step is the extraction of SURF keypoints from the input image (Figure5(b)). The speeded up robust features algorithm (SURF) [28] finds features using integral images and Haar-like filters. SURF features provide robustness and speed and are able to be detected despite scale, translation, and rotation of objects or changes in illumination. The following steps rely on this stability to group points and create object models, and they are described in the following subsections.SURF and SIFT detectors employ slightly different ways of detecting features. SIFT detector builds image pyramids and filters each layer with a Gaussian of increasing sigma values and then it takes the difference. SURF detector uses a box filter approximation. A comparison of the SIFT and SURF detectors is presented in [22]. In such work, the authors conclude thatSURF is as good as SIFT on most tests, except for scaling, large blur, and viewpoint. With respect to this, we have to mention that the proposed descriptor only uses SURF detector in the first step of the classifier to detect interest points. Then the proposed descriptor is constructed using additional 3D information and invariant moments. In the experiments results show that the proposed descriptor is small and robust. In [22], the authors also mention thaton real image data sets there is little to separate the different SIFTs and SURF except for efficiency. In this work we choose SURF detector because we are interested in implementing a real time classifier algorithm that could be used by a mobile robot. In addition, as we will see in the results’ section, the proposed descriptor is improved thanks to the spreading step in conjunction with the SURF detector, and this provides good results that are comparable and overcome similar approaches.


### 3.2. Contour Creation

The RGB-D sensor provides a huge amount of data; thus the information must be reduced into a contour-mesh model to decrease computational cost. To construct the contour-mesh we use keypoints. First, they are scaled and translated; then, we compute the magnitude of the keypoints with respect to their centroid(x-,y-) as
(1)m(x,y)=(x-x-)2+(y-y-)2
and its orientation with
(2)θ(x,y)=tan-1(y-y-x-x-).
After that, we divide the 360 degrees into orientation bins and use a sliding window to take one point for every bin according to their magnitude. If any of the bins remains empty, then its value is linearly interpolated using the previous and next known points in the contour. Examples for this stage of the process are shown in Figures 2 and 3.Figure 2
Totem object, with SURF keypoints.Figure 3
Sample contour for the totem object.In the next sections we will explain the last two steps of the proposed approach, namely, the mesh reduction step and the construction of the descriptor.


### 3.3. Mesh Reduction

In this step, the 3D data that belongs to the object contour is segmented, and then a cloud of 3D points is obtained. Due to the large number of points in the cloud, a reduction of the data will be required in order to keep a low computational cost. Therefore, in this step a mesh that covers the 3D point cloud is constructed, which allows to reduce the number of 3D points in the cloud, without losing important information.First, we proceed to extract depth information contained within the boundaries of the contour and reduce the points that will be taken to compute the moments. The reduction of points aims to generate a smaller set with rich information by adjusting a mesh grid over the object. The initial position of the points is obtained sectioning the bounding box, into equally separated cells, generating(x,y) coordinates. This creates a set of bidimensional points, that is,
(3)Grid={(x1,y1),(x2,y2),…,(xn,yn)},
where n is a predefined constant. Besides,
(4)xi∈[xmin⁡,xmax⁡],yi∈[ymin⁡,ymax⁡],x1=xmin⁡,y1=ymin⁡,xn=xmax⁡,yn=ymax⁡,xi+1=xi+xinc,yi+1=yi+yinc,xinc=xmax⁡-xmin⁡n,yinc=ymax⁡-ymin⁡n
minimum and maximum values are defined by the bounding box. The remaining points are tested using the point in polygon (PiP) algorithm, if the points are inside the polygon then they are considered valid.Then, for each invalid point, we take at random two valid points and move the outlier to a position between them,; this can be seen as a biased migration. Later, we attach thez coordinate (depth) to the (x,y) points. Examples for this stage of the process are shown in Figures 4 and 5(d).Figure 4
Mesh creation example.An image through the different stages of the procedure. (a) Original input image. (b) Feature extraction. (c) Contour created for the object. (d) Depth information attached in a mesh grid.
(a)
(b)
(c)
(d)It is important to mention that a point will be considered valid if and only if the coordinates(x,y) of the point are inside the contour, and the depth is not zero.These steps generate an object model for which we can extract information that we can use for classification; however the simple migration of points produces a model in which points are not equally distributed. This problem is solved by applying evolutionary computation.


## 3.1. Feature Extraction

The first step is the extraction of SURF keypoints from the input image (Figure5(b)). The speeded up robust features algorithm (SURF) [28] finds features using integral images and Haar-like filters. SURF features provide robustness and speed and are able to be detected despite scale, translation, and rotation of objects or changes in illumination. The following steps rely on this stability to group points and create object models, and they are described in the following subsections.SURF and SIFT detectors employ slightly different ways of detecting features. SIFT detector builds image pyramids and filters each layer with a Gaussian of increasing sigma values and then it takes the difference. SURF detector uses a box filter approximation. A comparison of the SIFT and SURF detectors is presented in [22]. In such work, the authors conclude thatSURF is as good as SIFT on most tests, except for scaling, large blur, and viewpoint. With respect to this, we have to mention that the proposed descriptor only uses SURF detector in the first step of the classifier to detect interest points. Then the proposed descriptor is constructed using additional 3D information and invariant moments. In the experiments results show that the proposed descriptor is small and robust. In [22], the authors also mention thaton real image data sets there is little to separate the different SIFTs and SURF except for efficiency. In this work we choose SURF detector because we are interested in implementing a real time classifier algorithm that could be used by a mobile robot. In addition, as we will see in the results’ section, the proposed descriptor is improved thanks to the spreading step in conjunction with the SURF detector, and this provides good results that are comparable and overcome similar approaches.


## 3.2. Contour Creation

The RGB-D sensor provides a huge amount of data; thus the information must be reduced into a contour-mesh model to decrease computational cost. To construct the contour-mesh we use keypoints. First, they are scaled and translated; then, we compute the magnitude of the keypoints with respect to their centroid(x-,y-) as
(1)m(x,y)=(x-x-)2+(y-y-)2
and its orientation with
(2)θ(x,y)=tan-1(y-y-x-x-).
After that, we divide the 360 degrees into orientation bins and use a sliding window to take one point for every bin according to their magnitude. If any of the bins remains empty, then its value is linearly interpolated using the previous and next known points in the contour. Examples for this stage of the process are shown in Figures 2 and 3.Figure 2
Totem object, with SURF keypoints.Figure 3
Sample contour for the totem object.In the next sections we will explain the last two steps of the proposed approach, namely, the mesh reduction step and the construction of the descriptor.


## 3.3. Mesh Reduction

In this step, the 3D data that belongs to the object contour is segmented, and then a cloud of 3D points is obtained. Due to the large number of points in the cloud, a reduction of the data will be required in order to keep a low computational cost. Therefore, in this step a mesh that covers the 3D point cloud is constructed, which allows to reduce the number of 3D points in the cloud, without losing important information.First, we proceed to extract depth information contained within the boundaries of the contour and reduce the points that will be taken to compute the moments. The reduction of points aims to generate a smaller set with rich information by adjusting a mesh grid over the object. The initial position of the points is obtained sectioning the bounding box, into equally separated cells, generating(x,y) coordinates. This creates a set of bidimensional points, that is,
(3)Grid={(x1,y1),(x2,y2),…,(xn,yn)},
where n is a predefined constant. Besides,
(4)xi∈[xmin⁡,xmax⁡],yi∈[ymin⁡,ymax⁡],x1=xmin⁡,y1=ymin⁡,xn=xmax⁡,yn=ymax⁡,xi+1=xi+xinc,yi+1=yi+yinc,xinc=xmax⁡-xmin⁡n,yinc=ymax⁡-ymin⁡n
minimum and maximum values are defined by the bounding box. The remaining points are tested using the point in polygon (PiP) algorithm, if the points are inside the polygon then they are considered valid.Then, for each invalid point, we take at random two valid points and move the outlier to a position between them,; this can be seen as a biased migration. Later, we attach thez coordinate (depth) to the (x,y) points. Examples for this stage of the process are shown in Figures 4 and 5(d).Figure 4
Mesh creation example.An image through the different stages of the procedure. (a) Original input image. (b) Feature extraction. (c) Contour created for the object. (d) Depth information attached in a mesh grid.
(a)
(b)
(c)
(d)It is important to mention that a point will be considered valid if and only if the coordinates(x,y) of the point are inside the contour, and the depth is not zero.These steps generate an object model for which we can extract information that we can use for classification; however the simple migration of points produces a model in which points are not equally distributed. This problem is solved by applying evolutionary computation.


## 4. Mesh Optimization

The mesh reduction step produces an object model; however the points are not equally distributed and thus an optimization step is required. For the optimization step, we have chosen an evolutionary computation (EC) technique, due to all the constraints of the mesh optimization problem.Evolutionary algorithms (EA) are stochastic search methods inspired by the behavior of swarms’ populations or natural biological behavior. In general there are five popular algorithms: genetic algorithms [29], memetic algorithms [30], particle swarm optimization [3], ant-colony optimization [31], and shuffled frog leaping algorithm [32]. In [33], the authors present a comparison study of these five EA; they conclude that the PSO algorithm performs better in general, with respect to the quality of the solution and the success rate. For these reasons, we decided to choose the PSO algorithm among the others.To solve the mesh optimization problem, the PSO algorithm (Algorithm1) was adapted in order to spread the set of (x,y) points that conform the mesh and to obtain better models with rich information about the distribution of the points. The purpose is to maximize the distance between points while maintaining them inside the boundaries conformed by the object contour.Algorithm 1:PSO algorithm.
(1)repeat
(2)  Calculate fitness value of each particle using the fitness function
(3)  Update local best if the current fitness value is better
(4)  Determine global best: take the best fitness particle and compare it to current best
(5)  For each particle
(6)  Calculate particle velocity according to (1)
(7)  Update particle position according to (2)
(8)until Stop criteria is metIn our approach, each PSO particle represents a point in the mesh. The problem has multiple boundaries; every point of the contour is one. Thus, we have to check if particles lie inside the polygon and clamp them after every update. Instead of gathering particles to a global best position they take positions separated uniformly from each other. This is obtained through the fitness function and a modification in the updating rules.


### 4.1. Description of the Method

The objective of this method is the construction of a mesh with the best distribution of points. The first step in the mesh construction is the determination of the object contour (Section3.2). This step uses features detected with SURF and a sampling technique to find the object contour. The second step is the mesh reduction step (Section 3.3); in this step the 3D data is used to enhance the object model. However, if we use all the 3D data that belong to the obtained contour, then the computational cost, required to process it, will be high. Therefore, the objective of this step is the construction of a mesh that covers the 3D point cloud, but with fewer points. The third step is the mesh optimization; this process is required since the mesh constructed in the previous step does not distribute efficiently. For this purpose, the PSO algorithm was adapted to spread the points over the object. The adapted PSO minimizes the distance of each particle and its neighbors with respect to the mean distance of all the particles with respect to its neighbors. The PSO particles are initialized with the coordinates of the points. The number of particles is therefore equal to the number of points. Finally, the object model is recovered from the best local particle value (Pi).


### 4.2. Particle Swarm Optimization Algorithm

PSO is a stochastic search method inspired from the behavior of swarm animals, like bird flocking and fish schooling. In PSO, particles, or solution candidates, move over the search space towards the zone with the best conditions using a cognitive component of the relative particle and a social component generated by the swarm (best local and global positions). This lets PSO to evolve social behavior and relative movement into global optimum solutions [34, 35].In the iterative process, the positionXi and velocity Vi of particles are updated according to the cognition component Pi and the social component G with
(5)Vi(t+1)=ωVi(t)+c1φ1(Pi-Xi)+c2φ2(G-Xi),Xi(t+1)=Xi(t)+Vi(t+1),
where c1, c2 are positive constants, φ1, φ2 are two random variables with uniform distribution between 0 and 1, and ω is the inertia weight which balances the effect of the previous velocity vector on the new one. The cognitive component, Pi, is updated by each particle when a better position is obtained. The social component, G, is updated when a new best position within the whole swarm is found. After initializing the swarm, in each iteration, the PSO basic steps are performed until the stop criterion is reached [36].For more details on PSO, the interested reader is referred to view [2, 3, 36, 37].


### 4.3. PSO Fitness Function

Instead of a single fitness function we undertake three steps to get a fitness value. First, we measure the distance of each particle to its nearest neighbor; this measure gives information about how separated is every particle. We only take the distance between the currenti-particle and a neighbor j-particle with(6)Di=(xi-xj)2+(yi-yj)2,
to keep a low computational load.Then we calculate the mean of distances(7)D¯=1N∑k=1N‍Dk,
where N is equal to the swarm size. This can be seen as the global value to maximize, since we want the points to be uniformly separated and yet lie inside the boundaries; that is, they have the same distance to each other and cover the space uniformly.Next, we compute the difference between the local distance of the particle to its neighbor and the mean global distance. Thus, the fitness ofith particle is defined as
(8)fi=(Di-D¯)2.This is the fitness value, by minimizing (8) every particle tries to minimize its own distance with respect to the mean distance of the swarm and by doing this we would obtain an approximation to a geometric uniform distribution of the swarm on the space; that is, we will obtain a swarm in which every particle is separated from its neighbors by a distance of D¯. Thus, the distance of the particle to its neighbor is the closest value to the mean estimated distance. This value is stored in by the cognitive component, Pi.Since our function is multiconstrained, we add a penalty term to the fitness function. Particles that go out of the boundaries are penalized according to the distance that they have to the centroid. We add a constantα which determines the influence of the penalization in the fitness value computation [36]; the fitness function with the penalization term is defined as
(9)fi=(Di-D¯)2+αT(X)|X-X-|,
with
(10)T(X)={1Particleisoutsidecontour0Particleisinsidecontour,
where X is a point (x,y) in the mesh represented by a particle, X- is the centroid (x-,y-), and T(X) is a point-in-polygon function added so only points that got outside the contour are penalized.Additional to the penalization term, we also use a preserve feasibility approach as explained by [36], since all our initialization particles are feasible we want to preserve the final result like that; thus Pi is only updated when the particle lies inside the contour.


### 4.4. PSO Update Formulas

The position update includes a variableβ that multiplies velocity:
(11)Xi(t+1)=Xi(t)+βVi(t+1).β is a term that defines how close are particles from boundaries; we want particles to be sufficiently spread, but we also want some particles close to the boundaries. Clamping and updates may cause particles to go out or be repositioned. If the particle is close to boundaries, β will be small and the update effect will be less. This value is thresholded so after certain value it becomes one and the effect of velocity applies normally; this way, particles far from boundaries are not affected. In conclusion, β can be seen as a function of the distance of the particle to the boundaries with
(12)Db=(xi-xboundary)2+(yi-yboundary)2
the distance of a point to the object boundary, and
(13)β={DbDb<Threshold1Otherwise.The threshold value depends on the range of the data set; it establishes how far to the border the data can be without being affected. In our case feature points are scaled to[-10,10], and β values in the range [0.2,1.0] were tested, resulting in 0.7 as the value that best performed.In the velocity update, the global coefficient has the effect on particles to move towards the best position of the swarm. Since in this application we do not need that effect, the term is replaced by one that makes particles get closer or farther from their neighbor as needed. This behavior is accomplished by using the sign of the distance between particles, and a constant valueρ to determine how fast particles are going to move towards or away from others; that is,
(14)Vi(t+1)=ωVi(t)+c1φ1(Pi-xi)+ρsign⁡(Xi-Xj),
where Xi,Xj are two particles (xi,yi), (xj,yj), with Xj being the closest particle to Xi. Another option is the inclusion of a scale factor of the current fitness value, since it gives a value of how much the particles are separated.


### 4.5. Details of the Modified PSO

The fitness function of the algorithm is defined to minimize the distance between the particlei, with respect to the mean distance of all the particles with respect to their corresponding neighbors; therefore at the end of the PSO iterations we obtain the best spread particles found by the algorithm. Although many different fitness functions can be defined, for this particular application the chosen fitness function shows a good performance, as it was demonstrated on the experimental results.With respect to the stop criteria of the algorithm, we cannot force PSO to stop when only the best valueG is close to zero, since it does not guarantee that the distances between each particle with respect to their neighbors are close to the mean distance. Instead, the algorithm stops when the particle best value Pi of each particle is close to zero, or in practice smaller than a certain threshold. In addition we can stop the algorithm if a certain threshold of iterations has been met.


### 4.6. Point Spreading Algorithm (PSA)

The inclusion of PSO in the mesh reduction step of the image classification system aims to spread points and create models that describe the entire object surface better. After the initial generation of points in a grid over the object and the migration of points that lie outside the object, the PSO variation is applied to the mesh of points. Particles are initialized with the coordinate values of points. The number of particles is therefore determined by the number of points that we want over the object. A maximum number of iterations is set and the algorithm is executed, until the stop criterion is reached. An object model is recovered from the local best position(Pi) that each particle generated. The main steps of the algorithm are portrayed in Algorithm 2.Algorithm 2:Mesh construction algorithm.
(1) Obtain object contour
(2) Generate starting point coordinates from bounding box
(3) Migrate non-valid points
(4) PSO spreadingIn Figures6 and 7 we can see an example of this procedure. The first picture shows the original image. The second image shows the initial coordinates of the mesh. The third image shows migration of points; some points are migrated to very close places and thus they are covering the object poorly. The fourth image shows the points after the PSO algorithm is applied; we can see that points are moved and the object surface is covered in a better way.Different stages of the mesh grid procedure for the hair dryer object. (a) Original input image. (b) Feature extraction. (c) Contour created for the object. (d) Depth information attached in a mesh grid.
(a)
(b)
(c)
(d)Different stages of the mesh grid procedure for the earth globe object. (a) Original input image. (b) Feature extraction. (c) Contour created for the object. (d) Depth information attached in a mesh grid.
(a)
(b)
(c)
(d)


## 4.1. Description of the Method

The objective of this method is the construction of a mesh with the best distribution of points. The first step in the mesh construction is the determination of the object contour (Section3.2). This step uses features detected with SURF and a sampling technique to find the object contour. The second step is the mesh reduction step (Section 3.3); in this step the 3D data is used to enhance the object model. However, if we use all the 3D data that belong to the obtained contour, then the computational cost, required to process it, will be high. Therefore, the objective of this step is the construction of a mesh that covers the 3D point cloud, but with fewer points. The third step is the mesh optimization; this process is required since the mesh constructed in the previous step does not distribute efficiently. For this purpose, the PSO algorithm was adapted to spread the points over the object. The adapted PSO minimizes the distance of each particle and its neighbors with respect to the mean distance of all the particles with respect to its neighbors. The PSO particles are initialized with the coordinates of the points. The number of particles is therefore equal to the number of points. Finally, the object model is recovered from the best local particle value (Pi).


## 4.2. Particle Swarm Optimization Algorithm

PSO is a stochastic search method inspired from the behavior of swarm animals, like bird flocking and fish schooling. In PSO, particles, or solution candidates, move over the search space towards the zone with the best conditions using a cognitive component of the relative particle and a social component generated by the swarm (best local and global positions). This lets PSO to evolve social behavior and relative movement into global optimum solutions [34, 35].In the iterative process, the positionXi and velocity Vi of particles are updated according to the cognition component Pi and the social component G with
(5)Vi(t+1)=ωVi(t)+c1φ1(Pi-Xi)+c2φ2(G-Xi),Xi(t+1)=Xi(t)+Vi(t+1),
where c1, c2 are positive constants, φ1, φ2 are two random variables with uniform distribution between 0 and 1, and ω is the inertia weight which balances the effect of the previous velocity vector on the new one. The cognitive component, Pi, is updated by each particle when a better position is obtained. The social component, G, is updated when a new best position within the whole swarm is found. After initializing the swarm, in each iteration, the PSO basic steps are performed until the stop criterion is reached [36].For more details on PSO, the interested reader is referred to view [2, 3, 36, 37].


## 4.3. PSO Fitness Function

Instead of a single fitness function we undertake three steps to get a fitness value. First, we measure the distance of each particle to its nearest neighbor; this measure gives information about how separated is every particle. We only take the distance between the currenti-particle and a neighbor j-particle with(6)Di=(xi-xj)2+(yi-yj)2,
to keep a low computational load.Then we calculate the mean of distances(7)D¯=1N∑k=1N‍Dk,
where N is equal to the swarm size. This can be seen as the global value to maximize, since we want the points to be uniformly separated and yet lie inside the boundaries; that is, they have the same distance to each other and cover the space uniformly.Next, we compute the difference between the local distance of the particle to its neighbor and the mean global distance. Thus, the fitness ofith particle is defined as
(8)fi=(Di-D¯)2.This is the fitness value, by minimizing (8) every particle tries to minimize its own distance with respect to the mean distance of the swarm and by doing this we would obtain an approximation to a geometric uniform distribution of the swarm on the space; that is, we will obtain a swarm in which every particle is separated from its neighbors by a distance of D¯. Thus, the distance of the particle to its neighbor is the closest value to the mean estimated distance. This value is stored in by the cognitive component, Pi.Since our function is multiconstrained, we add a penalty term to the fitness function. Particles that go out of the boundaries are penalized according to the distance that they have to the centroid. We add a constantα which determines the influence of the penalization in the fitness value computation [36]; the fitness function with the penalization term is defined as
(9)fi=(Di-D¯)2+αT(X)|X-X-|,
with
(10)T(X)={1Particleisoutsidecontour0Particleisinsidecontour,
where X is a point (x,y) in the mesh represented by a particle, X- is the centroid (x-,y-), and T(X) is a point-in-polygon function added so only points that got outside the contour are penalized.Additional to the penalization term, we also use a preserve feasibility approach as explained by [36], since all our initialization particles are feasible we want to preserve the final result like that; thus Pi is only updated when the particle lies inside the contour.


## 4.4. PSO Update Formulas

The position update includes a variableβ that multiplies velocity:
(11)Xi(t+1)=Xi(t)+βVi(t+1).β is a term that defines how close are particles from boundaries; we want particles to be sufficiently spread, but we also want some particles close to the boundaries. Clamping and updates may cause particles to go out or be repositioned. If the particle is close to boundaries, β will be small and the update effect will be less. This value is thresholded so after certain value it becomes one and the effect of velocity applies normally; this way, particles far from boundaries are not affected. In conclusion, β can be seen as a function of the distance of the particle to the boundaries with
(12)Db=(xi-xboundary)2+(yi-yboundary)2
the distance of a point to the object boundary, and
(13)β={DbDb<Threshold1Otherwise.The threshold value depends on the range of the data set; it establishes how far to the border the data can be without being affected. In our case feature points are scaled to[-10,10], and β values in the range [0.2,1.0] were tested, resulting in 0.7 as the value that best performed.In the velocity update, the global coefficient has the effect on particles to move towards the best position of the swarm. Since in this application we do not need that effect, the term is replaced by one that makes particles get closer or farther from their neighbor as needed. This behavior is accomplished by using the sign of the distance between particles, and a constant valueρ to determine how fast particles are going to move towards or away from others; that is,
(14)Vi(t+1)=ωVi(t)+c1φ1(Pi-xi)+ρsign⁡(Xi-Xj),
where Xi,Xj are two particles (xi,yi), (xj,yj), with Xj being the closest particle to Xi. Another option is the inclusion of a scale factor of the current fitness value, since it gives a value of how much the particles are separated.


## 4.5. Details of the Modified PSO

The fitness function of the algorithm is defined to minimize the distance between the particlei, with respect to the mean distance of all the particles with respect to their corresponding neighbors; therefore at the end of the PSO iterations we obtain the best spread particles found by the algorithm. Although many different fitness functions can be defined, for this particular application the chosen fitness function shows a good performance, as it was demonstrated on the experimental results.With respect to the stop criteria of the algorithm, we cannot force PSO to stop when only the best valueG is close to zero, since it does not guarantee that the distances between each particle with respect to their neighbors are close to the mean distance. Instead, the algorithm stops when the particle best value Pi of each particle is close to zero, or in practice smaller than a certain threshold. In addition we can stop the algorithm if a certain threshold of iterations has been met.


## 4.6. Point Spreading Algorithm (PSA)

The inclusion of PSO in the mesh reduction step of the image classification system aims to spread points and create models that describe the entire object surface better. After the initial generation of points in a grid over the object and the migration of points that lie outside the object, the PSO variation is applied to the mesh of points. Particles are initialized with the coordinate values of points. The number of particles is therefore determined by the number of points that we want over the object. A maximum number of iterations is set and the algorithm is executed, until the stop criterion is reached. An object model is recovered from the local best position(Pi) that each particle generated. The main steps of the algorithm are portrayed in Algorithm 2.Algorithm 2:Mesh construction algorithm.
(1) Obtain object contour
(2) Generate starting point coordinates from bounding box
(3) Migrate non-valid points
(4) PSO spreadingIn Figures6 and 7 we can see an example of this procedure. The first picture shows the original image. The second image shows the initial coordinates of the mesh. The third image shows migration of points; some points are migrated to very close places and thus they are covering the object poorly. The fourth image shows the points after the PSO algorithm is applied; we can see that points are moved and the object surface is covered in a better way.Different stages of the mesh grid procedure for the hair dryer object. (a) Original input image. (b) Feature extraction. (c) Contour created for the object. (d) Depth information attached in a mesh grid.
(a)
(b)
(c)
(d)Different stages of the mesh grid procedure for the earth globe object. (a) Original input image. (b) Feature extraction. (c) Contour created for the object. (d) Depth information attached in a mesh grid.
(a)
(b)
(c)
(d)


## 5. Invariant Moment Descriptor and Classification

Moments provide useful and compact information of a data set, such as its spread or dispersion. A pattern may be represented by a density distribution function, moments can be obtained for a set of points representing an object, and they can then be used to discriminate between objects [38]. The first order moments can be used to locate the centroid of the points’ distribution. If we compute the moments considering a translation to the centroid, we generate central moments which can be made scale invariant [39, 40]. The general equation for three-dimensional central moment (for short, 3D moment) is defined as
(15)μpqr=∑x‍∑y‍∑z‍(x-x-)p(y-y-)q(z-z-)rf(x,y,z),
where f(x,y,z) is a distribution function of the variables and (x-,y-,z-) is the centroid. The scaling is performed using
(16)ηpqr=μpqrμ000γ,γ=[p+q+r3]+1.In particular, [38] defines seven values, computed by normalizing central moments through order three that are invariant to object scale, position, and orientation. Tests were performed using different combinations of these seven values and then using the moments of order one to four. Similar results were obtained using the three first values defined by Hu and the moments of order one to four. Therefore, invariant moments can be computed from the reduced 3D set. In the proposed method we use the first four moments, over the three dimensions, then the values for p, q, and r are defined over the interval from 1 to 4, with these moments we define a descriptor of 12 elements to represent each object model. Such vector has the form
(17)d=[η100,η200,η300,η400,η010,η020,η030,η040,kkkη001,η002,η003,η004].Finally, the descriptor is given as input to a SVM and we get the result on whether the image contains the target object. The SVM [41] is an algorithm to solve classification and regression problems. It defines a subset on the train data composed by those samples that are closer to the decision area. It is based in the maximization of the margin of separation between classes. The SVM algorithm is able to perform nonlinear classification thanks to the use of kernel functions.


## 6. Results

In the following, we validate the proposed approach and compare it with histograms of oriented gradients (HOG) [42], scale invariant feature transform (SIFT) detector [10], and a detection system using cascades of HAAR-like features [43].First, a data set composed by 5 objects was defined (cups, hair dryers, irons, cereal boxes, and soda cans) with around 50 images for each object; see Figure8. The images included changes in the scene conditions such as illumination, object orientation, and position, besides partially occluded objects.Figure 8
Objects conforming the data set.In addition to the house made data set, the RGB-D object data set from the University of Washington [44] was used to perform tests. We used a set consisting of eight different objects (bowl, cap, flashlight, coffee mug, cereal box, soda can, camera, and pitcher) similar to [45] (Figure 9). This data set was selected since it contains similar objects to our set and they were acquired through the same range sensor. For each of the classes, the set includes a group of around four objects that belong to the same class with three different views and 200 images for each view. The bounding box, cropped image, and object mask are provided.Figure 9
Objects from the RGB-D data set of the University of Washington [44].


### 6.1. Algorithm Parameters

We worked with images containing a single object with a discriminative background. As explained before, we start by extracting key points from the image. The contours were formed by 72 bins of 5 degrees each. The mesh was composed by a9×9 grid, generating 81 points. The algorithmic control parameters of PSO, coefficient of cognition (c1), and inertia weight (ω) (14) were set to ω=0.729844 and c1=1.49618 as suggested in [46]. Different values for the parameters α (9) and ρ (14) were tested; the best results were obtained with α=1.230 and ρ=3.751, though this might be application dependent.In our case, forρ, values greater than 5 resulted too aggressive bringing neighbor particles too close rapidly, and values smaller than 2 had little effect on the movement of particles. Values for α within the range [0.5,1.5] caused the expected penalty result while other values caused fitness values to boost or to be decreased abnormally. Finally, as explained before β (11) was set to 0.7 defining that particles closer to the boundaries will be less affected by the velocity update.


### 6.2. Training

The cross-validation method was used to validate the training process [47]. Using a 4-fold approach in the five-object set and a 5-fold approach in the six-object set, the data was randomly separated in different subsets with equal number of elements; three of the subsets were used to train different SVMs and the rest were used for validation. The classifiers were trained and tested 5 times rotating the subsets used for training; the estimate of accuracy is the overall number of correct classifications divided by the number of instances in the dataset; finally, the classifier performance is measured with the average of the accuracy throughout the rotation process.To define the SVM model, different SVMs with linear,xi′·xj, and polynomial, (αxi′·xj+b)d, kernel functions were tested varying the cost parameter C using values in the range (102,108), and the degree of the polynomial d using values in the range (3,10); cross-validation was used to define which value is the best. A supervised learning approach was considered in the whole training and validation processes; the classifiers were trained using a labeled dataset of images, in each image the information about which object is present and the bounding box of the object was also provided. Therefore, the class to which each descriptor belongs is known. To classify each object 100 tests were performed.


### 6.3. Classification

The first experiment consisted of binary classification of one object being discriminated from another, 1 versus 1 classification; see Table1. The results are summarized in a result table. Each row indicates the percentage of (a) correct recognitions (CR) and (b) class one errors (C1E); objects of class one are classified as objects of class two and (C2E) class two errors; objects of class two are classified as objects of class one.Table 1
Binary classification using 2 objects.
1 versus 1 classification percentages
Test
PSA-moment descriptor
Name
CR
C1E
C2E
Dryer and box
100%
0%
0%
Can and irons
94%
3%
3%
Cup and can
86%
10%
4%
Iron and dryer
88%
4%
8%
Box and cup
90%
0%
10%The second test consisted of binary classification using the 5 objects (an object being discriminated from the rest); see Table2. In this case, result table contains the percentages of (a) correct recognitions (CR), (b) false positives (FP), a different object is classified as the target object, and (c) false negatives (FN), the target object is classified as a different object.Table 2
Binary classification using 5 objects.
1 versus all classification percentages
Test
PSA-moment descriptor
Name
CR
FP
FN
Dryer
96%
2%
2%
Iron
88%
0%
12%
Cup
92%
2%
6%
Can
90%
10%
0%
Box
96%
2%
2%Finally, a test was made in multiclass classification where objects were classified all at once. Each object was set as a different class and the classification delivered the class to which the objects belong. This test generated a 79% of correct classification. An 8% of incorrect classification was present in class 1 (cups) and 9% in class 2 (irons); all objects of classes 3 and 4 (irons and cans) were correctly classified, and 4% of incorrect classification was in class 5 (box).After testing the descriptors with the small data set, classification tests were performed using a group of objects from the RGB-D object data set of the University of Washington. All objects were discriminated from the rest in binary classification. Tests were performed using first the moment descriptor obtained when points in the net are simply migrated, then the same tests were done using the point spreading approach. Later, these results were compared with histograms of oriented gradients (HOG) [42], scale invariant feature transform (SIFT) detector [10], and a detection system using cascades of HAAR-like features [43].The classification has been made for each descriptor (moment, moment-PSA, HOG, and SIFT) using SVMs. From cross-validation the best results obtained were a linear kernel function with a cost parameterC=103. Results of the tests for the four objects are summarized in Table 3: object one is the bowl, object two is the pitcher, object three is the cap, and object four is the soda can.Table 3
Tests using objects from the Washington University data set.
Classification percentages
Object
Object 1
Object 2
Object 3
Object 4
Result
CR
FN
FP
CR
FN
FP
CR
FN
FP
CR
FN
FP
Moment descriptor
80
8
12
96
0
4
92
4
4
79
7
14
Moment-PSA descriptor
85
5
10
98
0
2
97
1
2
82
8
10
HOG descriptor
99
1
0
98
2
0
90
8
2
85
10
5
SIFT descriptor
74
22
4
85
8
7
76
10
14
70
27
3
HAAR cascade
90
9
1
89
7
4
75
10
15
72
27
1We can see from these results that spreading the points over the object surface improves the descriptor computed and then the classification of the objects. With this approach we have similar results to the well-known HOG descriptors which are used in many of the state-of-the-art recognition systems. These two descriptors achieved the highest correct recognition percentages in all the tests.The accuracy of the classification system has been assessed in terms of receiving operating characteristics curves (ROC) which relate the positive and false acceptance rates according to an acceptance thresholdδ varying in the rate [0,1]. Figure 10 shows the ROC curves for the classification of the bowl and Figure 11 the ROC curves for the classification of the coffee mug.ROC plots for the classification of the bowl. (a) HOG ROC, (b) moments ROC, (c) moments-PSA ROC, (d) HAAR ROC.
(a)
(b)
(c)
(d)ROC plots for the classification of the coffee mug. (a) HOG ROC, (b) moments ROC, (c) moments-PSA ROC, (d) HAAR ROC.
(a)
(b)
(c)
(d)With the ROC curves, we can also note that the moments-PSA classifiers show improvement over the one with the grid that is not optimized, and also it shows similar performance to the classification through HOG descriptors.The efficiency of the descriptors using a SVM classifier was also evaluated in terms of the number of support vectors obtained by the SVM algorithm. This was done in the cases of moment, moment-PSA, HOG, and SIFT descriptors. Table4 shows the average and minimum number of support vectors for 100 tests on the classification of the four objects used in Table 3.Table 4
Average and minimum number of support vectors.
Support vectors
Object
Object 1
Object 2
Object 3
Object 4
Support vectors
avg.
min.
avg.
min.
avg.
min.
avg.
min.
Moment descriptor
20
15
25
17
101
73
27
20
Moment-PSA descriptor
8
5
24
12
73
52
12
10
HOG descriptor
300
23
400
60
370
43
418
70
SIFT descriptor
310
22
404
65
401
84
450
72With these results we can see that although HOG and our approach have similar classification results, the number of support vectors needed in the training process is smaller for our approach, and when the distribution of the points is optimized it also results in patterns that are better distributed in the feature space so we need even less support vectors.In this paper we propose the design of a new small and robust feature descriptor based on SURF features, 3D data obtained by RGB-D sensor, an optimization based on a modified PSO, and moment invariants. The experimental results presented in this section show that even that the new descriptors are small the addition of 3D data is advantageous and provides robust features. The importance of the mesh optimization step lies on the reduction of points and more distributed points over the object. Table4 shows that the proposed approach requires less support vectors; therefore the classification is faster. From the ROC curves we can also note that the proposed approach has a performance compared with the HOG approach, but as we mention, with fewer support vectors.In order to determine the time required to construct the descriptor, we have tested each step 100 times for each object to estimate its average processing time, a resume of this results is presented in Table5. In this table, we can note that the most time consuming process is the feature extraction process and that the computation of the invariants moments of the optimized mesh is insignificant. With respect to the classification time, we use the same process used to estimate the descriptor construction time, and we obtained an average of 0.007565 s. The implementation of the algorithm has been made in Matlab, and therefore it can be improved by a C or C++ implementation.Table 5
Descriptor construction time. Where the steps are feature extraction (FE), contour creation (CC), mesh reduction (MR), mesh optimization (MO), and invariant moments (IM).
Step
FE
CC
MR
MO
IM
Time
0.0929 s
0.0419 s
0.0202 s
0.0566 s
7.85
×
10
-
4


## 6.1. Algorithm Parameters

We worked with images containing a single object with a discriminative background. As explained before, we start by extracting key points from the image. The contours were formed by 72 bins of 5 degrees each. The mesh was composed by a9×9 grid, generating 81 points. The algorithmic control parameters of PSO, coefficient of cognition (c1), and inertia weight (ω) (14) were set to ω=0.729844 and c1=1.49618 as suggested in [46]. Different values for the parameters α (9) and ρ (14) were tested; the best results were obtained with α=1.230 and ρ=3.751, though this might be application dependent.In our case, forρ, values greater than 5 resulted too aggressive bringing neighbor particles too close rapidly, and values smaller than 2 had little effect on the movement of particles. Values for α within the range [0.5,1.5] caused the expected penalty result while other values caused fitness values to boost or to be decreased abnormally. Finally, as explained before β (11) was set to 0.7 defining that particles closer to the boundaries will be less affected by the velocity update.


## 6.2. Training

The cross-validation method was used to validate the training process [47]. Using a 4-fold approach in the five-object set and a 5-fold approach in the six-object set, the data was randomly separated in different subsets with equal number of elements; three of the subsets were used to train different SVMs and the rest were used for validation. The classifiers were trained and tested 5 times rotating the subsets used for training; the estimate of accuracy is the overall number of correct classifications divided by the number of instances in the dataset; finally, the classifier performance is measured with the average of the accuracy throughout the rotation process.To define the SVM model, different SVMs with linear,xi′·xj, and polynomial, (αxi′·xj+b)d, kernel functions were tested varying the cost parameter C using values in the range (102,108), and the degree of the polynomial d using values in the range (3,10); cross-validation was used to define which value is the best. A supervised learning approach was considered in the whole training and validation processes; the classifiers were trained using a labeled dataset of images, in each image the information about which object is present and the bounding box of the object was also provided. Therefore, the class to which each descriptor belongs is known. To classify each object 100 tests were performed.


## 6.3. Classification

The first experiment consisted of binary classification of one object being discriminated from another, 1 versus 1 classification; see Table1. The results are summarized in a result table. Each row indicates the percentage of (a) correct recognitions (CR) and (b) class one errors (C1E); objects of class one are classified as objects of class two and (C2E) class two errors; objects of class two are classified as objects of class one.Table 1
Binary classification using 2 objects.
1 versus 1 classification percentages
Test
PSA-moment descriptor
Name
CR
C1E
C2E
Dryer and box
100%
0%
0%
Can and irons
94%
3%
3%
Cup and can
86%
10%
4%
Iron and dryer
88%
4%
8%
Box and cup
90%
0%
10%The second test consisted of binary classification using the 5 objects (an object being discriminated from the rest); see Table2. In this case, result table contains the percentages of (a) correct recognitions (CR), (b) false positives (FP), a different object is classified as the target object, and (c) false negatives (FN), the target object is classified as a different object.Table 2
Binary classification using 5 objects.
1 versus all classification percentages
Test
PSA-moment descriptor
Name
CR
FP
FN
Dryer
96%
2%
2%
Iron
88%
0%
12%
Cup
92%
2%
6%
Can
90%
10%
0%
Box
96%
2%
2%Finally, a test was made in multiclass classification where objects were classified all at once. Each object was set as a different class and the classification delivered the class to which the objects belong. This test generated a 79% of correct classification. An 8% of incorrect classification was present in class 1 (cups) and 9% in class 2 (irons); all objects of classes 3 and 4 (irons and cans) were correctly classified, and 4% of incorrect classification was in class 5 (box).After testing the descriptors with the small data set, classification tests were performed using a group of objects from the RGB-D object data set of the University of Washington. All objects were discriminated from the rest in binary classification. Tests were performed using first the moment descriptor obtained when points in the net are simply migrated, then the same tests were done using the point spreading approach. Later, these results were compared with histograms of oriented gradients (HOG) [42], scale invariant feature transform (SIFT) detector [10], and a detection system using cascades of HAAR-like features [43].The classification has been made for each descriptor (moment, moment-PSA, HOG, and SIFT) using SVMs. From cross-validation the best results obtained were a linear kernel function with a cost parameterC=103. Results of the tests for the four objects are summarized in Table 3: object one is the bowl, object two is the pitcher, object three is the cap, and object four is the soda can.Table 3
Tests using objects from the Washington University data set.
Classification percentages
Object
Object 1
Object 2
Object 3
Object 4
Result
CR
FN
FP
CR
FN
FP
CR
FN
FP
CR
FN
FP
Moment descriptor
80
8
12
96
0
4
92
4
4
79
7
14
Moment-PSA descriptor
85
5
10
98
0
2
97
1
2
82
8
10
HOG descriptor
99
1
0
98
2
0
90
8
2
85
10
5
SIFT descriptor
74
22
4
85
8
7
76
10
14
70
27
3
HAAR cascade
90
9
1
89
7
4
75
10
15
72
27
1We can see from these results that spreading the points over the object surface improves the descriptor computed and then the classification of the objects. With this approach we have similar results to the well-known HOG descriptors which are used in many of the state-of-the-art recognition systems. These two descriptors achieved the highest correct recognition percentages in all the tests.The accuracy of the classification system has been assessed in terms of receiving operating characteristics curves (ROC) which relate the positive and false acceptance rates according to an acceptance thresholdδ varying in the rate [0,1]. Figure 10 shows the ROC curves for the classification of the bowl and Figure 11 the ROC curves for the classification of the coffee mug.ROC plots for the classification of the bowl. (a) HOG ROC, (b) moments ROC, (c) moments-PSA ROC, (d) HAAR ROC.
(a)
(b)
(c)
(d)ROC plots for the classification of the coffee mug. (a) HOG ROC, (b) moments ROC, (c) moments-PSA ROC, (d) HAAR ROC.
(a)
(b)
(c)
(d)With the ROC curves, we can also note that the moments-PSA classifiers show improvement over the one with the grid that is not optimized, and also it shows similar performance to the classification through HOG descriptors.The efficiency of the descriptors using a SVM classifier was also evaluated in terms of the number of support vectors obtained by the SVM algorithm. This was done in the cases of moment, moment-PSA, HOG, and SIFT descriptors. Table4 shows the average and minimum number of support vectors for 100 tests on the classification of the four objects used in Table 3.Table 4
Average and minimum number of support vectors.
Support vectors
Object
Object 1
Object 2
Object 3
Object 4
Support vectors
avg.
min.
avg.
min.
avg.
min.
avg.
min.
Moment descriptor
20
15
25
17
101
73
27
20
Moment-PSA descriptor
8
5
24
12
73
52
12
10
HOG descriptor
300
23
400
60
370
43
418
70
SIFT descriptor
310
22
404
65
401
84
450
72With these results we can see that although HOG and our approach have similar classification results, the number of support vectors needed in the training process is smaller for our approach, and when the distribution of the points is optimized it also results in patterns that are better distributed in the feature space so we need even less support vectors.In this paper we propose the design of a new small and robust feature descriptor based on SURF features, 3D data obtained by RGB-D sensor, an optimization based on a modified PSO, and moment invariants. The experimental results presented in this section show that even that the new descriptors are small the addition of 3D data is advantageous and provides robust features. The importance of the mesh optimization step lies on the reduction of points and more distributed points over the object. Table4 shows that the proposed approach requires less support vectors; therefore the classification is faster. From the ROC curves we can also note that the proposed approach has a performance compared with the HOG approach, but as we mention, with fewer support vectors.In order to determine the time required to construct the descriptor, we have tested each step 100 times for each object to estimate its average processing time, a resume of this results is presented in Table5. In this table, we can note that the most time consuming process is the feature extraction process and that the computation of the invariants moments of the optimized mesh is insignificant. With respect to the classification time, we use the same process used to estimate the descriptor construction time, and we obtained an average of 0.007565 s. The implementation of the algorithm has been made in Matlab, and therefore it can be improved by a C or C++ implementation.Table 5
Descriptor construction time. Where the steps are feature extraction (FE), contour creation (CC), mesh reduction (MR), mesh optimization (MO), and invariant moments (IM).
Step
FE
CC
MR
MO
IM
Time
0.0929 s
0.0419 s
0.0202 s
0.0566 s
7.85
×
10
-
4


## 7. Discussion

### 7.1. Difficulties of Each Step of the System

In the previous section the different steps of the proposed approach had been presented. The approach produces a descriptor based on moment invariants which has been computed using 3D data provided by an RGB-D sensor. A new method has been introduced for the extraction of object models from the 3D data.The first step of the procedure is the feature extraction; in this approach the SURF detector was used, and the difficulties of obtaining the SURF features are the traditional of any image feature detector. The second step is the contour creation: this step takes as input the features detected in the previous step and extracts the contour of the object. In this step we must take into account that some interest points do not belong to the object; they are close to the contour but they are outside of it. This point can be seen as noise, and they are eliminated if their distance to the centroid is the double of the mean distance.After the construction of the object contour, the depth information is added. If the area occupied by the object is big, then there would be a lot of data; thus the next steps will require a lot of computational process. For this reason, the third step of the proposed approach tries to optimize the mesh in order to reduce the computational time. In this step a modified version of PSO has been used to optimize the points on the mesh. One of the challenges encountered on this step was the stop criteria; we cannot force PSO to stop when the best value is zero, and therefore the algorithm stops when the best value of each particlei  (Pi) is smaller than a certain threshold. With respect to the SVM, the only problem was the estimation of the cost parameter C, but this was solved by using cross-validation.


### 7.2. Computation Complexity of the System

Taking into account the steps illustrated in Figure1, next, a complexity analysis of the classification system per step is presented.With respect to the feature extraction step, as we mentioned before, it is made by using a SURF detector. The theoretical complexity of SURF was determined and validated through experimentation in [48]; this is O(MN+K), where M and N are the image height and width, respectively, and K stands for the number of detected points. It was proven that SURF is a very attractive algorithm concerning its performance when usual size and high rates of frames by seconds are used; this is our case; the images’ size is around 640×480 pixels. This step produces a dense 3D point cloud that is greatly reduced in the next steps.The contour creation step is a linear process with respect to the number of points of the cloud produced by the SURF algorithm; nevertheless when the contour is created this number is reduced and only the points which belongs to the silhouette of the object are taken to be processed by the mesh reduction step, which is also linear with respect to this reduced number of points.The mesh optimization step is performed by using the PSO algorithm. The number of computations to complete a PSO run is the computations of the cost function and the position and velocity updates, which are directly proportional to the number of particles and iterations of the algorithm. The computational complexity of evaluating the cost function depends basically on the Euclidean distance of the particles, which requires2N additions, 2N multiplications, and N square roots, where N represents the swarm size. Whereas the position update requires 2N multiplications and 2N additions, and the velocity update requires 7N multiplications and 4N additions. Therefore in total the PSO requires 6N additions, 11N multiplications, and N square roots. Again this step reduces even more the number of points which belong to the contour of the object.The invariant moment descriptor generation seems to be the more expensive step which is executed in real time by the system because it is a process with a computation complexity ofO(N3), with N the number of points of each contour, but this number is small (because it was reduced by the above three steps of the system) and the descriptors are very fast computed as it is shown in Table 5.The last step of the system is the classification of the descriptors obtained; it is made using a SVM. A SVM solves a quadratic programming problem (QP) and therefore its computation complexity depends on the QP solver used. The computation complexity of the traditional algorithm of the SVM classifier isO(Nsv3+LNsv2+dlLNsv) [49], where dl is the dimension of the input space, L is the number of training points, the number of the samples in the support vector set is Nsv. Nevertheless there are several efficient implementations, such as incremental learning algorithm [50] which reports a computation complexity of O(Nsv3+dlNsv2), or sequential minimal optimization (SMO) [51] which for a real-world test sets are 1200 times faster for linear SVMs and 15 times faster for nonlinear SVMs against another quadratic programming solver techniques. SMO was the training method used in our system. It is important to note that, even when the system includes a training step with the complexity that was already determined for the SVM algorithm, this step is made offline and, in the real time process, the classification step consists of evaluating the classification function obtained with the training of the SVM and this is a linear process with respect to the number of support vectors Nsv, which is usually related with the number of training data N as follows Nsv≪N.


## 7.1. Difficulties of Each Step of the System

In the previous section the different steps of the proposed approach had been presented. The approach produces a descriptor based on moment invariants which has been computed using 3D data provided by an RGB-D sensor. A new method has been introduced for the extraction of object models from the 3D data.The first step of the procedure is the feature extraction; in this approach the SURF detector was used, and the difficulties of obtaining the SURF features are the traditional of any image feature detector. The second step is the contour creation: this step takes as input the features detected in the previous step and extracts the contour of the object. In this step we must take into account that some interest points do not belong to the object; they are close to the contour but they are outside of it. This point can be seen as noise, and they are eliminated if their distance to the centroid is the double of the mean distance.After the construction of the object contour, the depth information is added. If the area occupied by the object is big, then there would be a lot of data; thus the next steps will require a lot of computational process. For this reason, the third step of the proposed approach tries to optimize the mesh in order to reduce the computational time. In this step a modified version of PSO has been used to optimize the points on the mesh. One of the challenges encountered on this step was the stop criteria; we cannot force PSO to stop when the best value is zero, and therefore the algorithm stops when the best value of each particlei  (Pi) is smaller than a certain threshold. With respect to the SVM, the only problem was the estimation of the cost parameter C, but this was solved by using cross-validation.


## 7.2. Computation Complexity of the System

Taking into account the steps illustrated in Figure1, next, a complexity analysis of the classification system per step is presented.With respect to the feature extraction step, as we mentioned before, it is made by using a SURF detector. The theoretical complexity of SURF was determined and validated through experimentation in [48]; this is O(MN+K), where M and N are the image height and width, respectively, and K stands for the number of detected points. It was proven that SURF is a very attractive algorithm concerning its performance when usual size and high rates of frames by seconds are used; this is our case; the images’ size is around 640×480 pixels. This step produces a dense 3D point cloud that is greatly reduced in the next steps.The contour creation step is a linear process with respect to the number of points of the cloud produced by the SURF algorithm; nevertheless when the contour is created this number is reduced and only the points which belongs to the silhouette of the object are taken to be processed by the mesh reduction step, which is also linear with respect to this reduced number of points.The mesh optimization step is performed by using the PSO algorithm. The number of computations to complete a PSO run is the computations of the cost function and the position and velocity updates, which are directly proportional to the number of particles and iterations of the algorithm. The computational complexity of evaluating the cost function depends basically on the Euclidean distance of the particles, which requires2N additions, 2N multiplications, and N square roots, where N represents the swarm size. Whereas the position update requires 2N multiplications and 2N additions, and the velocity update requires 7N multiplications and 4N additions. Therefore in total the PSO requires 6N additions, 11N multiplications, and N square roots. Again this step reduces even more the number of points which belong to the contour of the object.The invariant moment descriptor generation seems to be the more expensive step which is executed in real time by the system because it is a process with a computation complexity ofO(N3), with N the number of points of each contour, but this number is small (because it was reduced by the above three steps of the system) and the descriptors are very fast computed as it is shown in Table 5.The last step of the system is the classification of the descriptors obtained; it is made using a SVM. A SVM solves a quadratic programming problem (QP) and therefore its computation complexity depends on the QP solver used. The computation complexity of the traditional algorithm of the SVM classifier isO(Nsv3+LNsv2+dlLNsv) [49], where dl is the dimension of the input space, L is the number of training points, the number of the samples in the support vector set is Nsv. Nevertheless there are several efficient implementations, such as incremental learning algorithm [50] which reports a computation complexity of O(Nsv3+dlNsv2), or sequential minimal optimization (SMO) [51] which for a real-world test sets are 1200 times faster for linear SVMs and 15 times faster for nonlinear SVMs against another quadratic programming solver techniques. SMO was the training method used in our system. It is important to note that, even when the system includes a training step with the complexity that was already determined for the SVM algorithm, this step is made offline and, in the real time process, the classification step consists of evaluating the classification function obtained with the training of the SVM and this is a linear process with respect to the number of support vectors Nsv, which is usually related with the number of training data N as follows Nsv≪N.


## 8. Conclusions

In this paper we have presented an image classification technique based on an invariant moment descriptor that includes depth information. The inclusion of 3D data enables invariant moments to produce small and robust descriptors improving image classification. To create such descriptors, we used object models with rich information of the distribution of points. The application of the optimization algorithm PSO to the model generation stage improved the computed descriptor and the object recognition. From experimental results, it is clear that these descriptors have achieved high correct recognition percentages. Furthermore, the number of support vectors obtained in the training process is smaller for our approach, due to the fact that the points are optimized and thus the patterns are better distributed in the feature space.



---

*Source: 695910-2014-07-10.xml*