Next Article in Journal
Schizophrenia Detection Using Machine Learning Approach from Social Media Content
Previous Article in Journal
A Systematic Survey on the Role of Cloud, Fog, and Edge Computing Combination in Smart Agriculture
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Object Positioning Algorithm Based on Multidimensional Scaling and Optimization for Synthetic Gesture Data Generation

by
Borja Saez-Mingorance
1,2,
Antonio Escobar-Molero
1,
Javier Mendez-Gomez
1,2,
Encarnacion Castillo-Morales
2 and
Diego P. Morales-Santos
2,3,*
1
Infineon Technologies AG, Am Campeon 1-15, 85579 Neubiberg, Germany
2
Department of Electronic and Computer Technology, University of Granada, Avenida de Fuente Nueva s/n, 18071 Granada, Spain
3
RedNodeLabs UG, 80469 Munich, Germany
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(17), 5923; https://doi.org/10.3390/s21175923
Submission received: 30 July 2021 / Revised: 26 August 2021 / Accepted: 31 August 2021 / Published: 3 September 2021
(This article belongs to the Section Navigation and Positioning)

Abstract

:
This work studies the feasibility of a novel two-step algorithm for infrastructure and object positioning, using pairwise distances. The proposal is based on the optimization algorithms, Scaling-by-Majorizing-a-Complicated-Function and the Limited-Memory-Broyden-Fletcher-Goldfarb-Shannon. A qualitative evaluation of these algorithms is performed for 3D positioning. As the final stage, smoothing filtering techniques are applied to estimate the trajectory, from the previously obtained positions. This approach can also be used as a synthetic gesture data generator framework. This framework is independent from the hardware and can be used to simulate the estimation of trajectories from noisy distances gathered with a large range of sensors by modifying the noise properties of the initial distances. The framework is validated, using a system of ultrasound transceivers. The results show this framework to be an efficient and simple positioning and filtering approach, accurately reconstructing the real path followed by the mobile object while maintaining low latency. Furthermore, these capabilities can be exploited by using the proposed algorithms for synthetic data generation, as demonstrated in this work, where synthetic ultrasound gesture data are generated.

1. Introduction

Ultrasound technology is widely used for object positioning. Applications such as robot navigation [1], indoor navigation [2], human-device interface systems [3], body-tracking [4] or medical-probes tracking [5] are just a few examples of the many potential applications based on ultrasonic waves.
Positioning systems usually have a set of fixed anchor nodes that defines the infrastructure for the location. To locate the target, there are mainly two different approaches:
  • Locating an active object, one able to emit and/or receive ultrasonic signals [6,7].
  • Locating a passive object, one which just reflects the incoming ultrasonic wave emitted by the anchors [8].
Within the active-object alternative, one option is to use the anchors as receivers and the mobile node as the signal emitter. Based on the Time of Flight (ToF) or on the Received Signal Strength Indicator (RSSI), the anchors calculate their distance to the object [9]. Another alternative is using the Angle of Arrival (AoA) [10], where the position is obtained from the direction of arrival of the signal to the receiver. This work focuses on mechanisms based on ToF measurements, which are generally more robust and accurate, by relying on the predictable velocity of the ultrasonic wave in the air. If a ToF-based mechanism is used, all the anchors require an additional synchronization mechanism to have a common clock reference. Similarly, the roles can be inverted and the anchors can synchronously transmit beacons, with the mobile node as the receiver, which computes the distances locally. To achieve the time synchronization between the nodes, a combination of different technologies may be used in the same system, such as ultra wide band (UWB) and ultrasounds [11]. Another popular approach to avoid the requirement of having a tight synchronization mechanism between the anchors is to use two-way ranging mechanisms [12] in which either the mobile node or the anchors reply with another signal after a fixed amount of time and the one-to-one distances are computed based on the individual round-trip times.
There are works that use the active-object approach, based on ultrasound technology, for positioning and tracking: Chen H. et al. [9] proposed a system, where the positioning is based in a fixed receiver array performing the localization of a transmitter array attached to the hand of the user. Chen J. et al. [13] described using ultrasonic signal and radio signal together to develop a transmitting 3D pen, and the algorithm to position the pen based on a set of receiving nodes covering the writing plane.
In the passive-object alternative, just the echo, or reflected wave, is detected back by the anchors. This is typically feasible for very short-range applications, such as gesture recognition [14], in which the surface to locate is interference-free and has a reflective surface large enough to be easily recognized. The same distance measuring techniques used in the active alternative can be used with the passive alternative [15], taking into account the characteristics of the passive approach.
Ultimately, the accuracy and robustness of the system rely on the dependability of the distance measurements. It is critical to recognize the incoming signal (either reflected or actively transmitted by another node) over the ultrasonic background noise. Several methods are proposed in the literature based on different criteria (time, frequency, phase) [16], with the most popular being the cross correlation of the received and expected signal. It requires low computational power, introduces low delay and offers higher robustness against noise when detecting an echo [13,16].
This technique enables the emission of different signals (i.e., in the case where the anchors play the emitter role) to differentiate between incoming pulses, such as Direct Sequence Code-Division Multiple Access (DS-CDMA) [17,18].
Once the distances to the anchors are obtained, the location of the object can be determined, using a positioning algorithm, one based on the trilateration concept [19,20]. Knowing the position of three anchors A 1 ( x 1 , y 1 , z 1 ) , A 2 ( x 2 , y 2 , z 2 ) , A 3 ( x 3 , y 3 , z 3 ) and the pairwise distances ( d 1 , d 2 and d 3 ), the coordinates of the object, O ( x , y , z ) can be calculated solving the following system of equations:
( x x 1 ) 2 + ( y y 1 ) 2 + ( z z 1 ) 2 = d 1 2 ( x x 2 ) 2 + ( y y 2 ) 2 + ( z z 2 ) 2 = d 2 2 ( x x 3 ) 2 + ( y y 3 ) 2 + ( z z 3 ) 2 = d 3 2
This algebraic solution corresponds to the cross points of the three spheres with center A 1 , A 2 and A 3 , and radii d 1 , d 2 and d 3 , respectively.
Three anchor nodes, and the distance from all three anchors to the object, are needed as a minimum requirement to obtain the 3D location of the object. If there are fewer than three distances to the anchor nodes (i.e., there is no direct acoustic channel between the object and one anchor), it is not possible to determine the location.
When there are more than three anchors involved in the location, we have an overdetermined system, and the method is called multilateration. Its advantage is a potentially increased robustness against inaccurate or missing distances. With N anchors, it is required to solve a system with N equations, making necessary the use of recursive algorithms to obtain an optimal solution [21]:
( x x 1 ) 2 + ( y y 1 ) 2 + ( z z 1 ) 2 = d 1 2 ( x x 2 ) 2 + ( y y 2 ) 2 + ( z z 2 ) 2 = d 2 2 ( x x N ) 2 + ( y y N ) 2 + ( z z N ) 2 = d N 2
To compute the coordinates of the object, and even to previously or simultaneously position the anchor infrastructure, fast and robust algorithms are required. They should be able to easily adapt to a varying number of noisy distances, and therefore are not totally reliable. Furthermore, if trajectories are to be obtained, a further processing step is useful to smooth out the path of the object and improve the accuracy of the estimated track.
In this paper, several approaches to achieve the object location and tracking are proposed using multidimensional scaling (MDS) and optimization algorithms. A qualitative evaluation of these algorithms is performed in this work. In addition, the integration of the algorithms in a synthetic data generation framework is discussed. This use case shows how the dataset creation task, i.e., ultrasound gesture dataset, could benefit from these algorithms due to the high flexibility to configure the desired output with different noise levels and gesture options. At the same time, since the desired data are configured by the user, this framework would generate simultaneously data and labels. By applying this framework, the possibility of incurring human error is reduced, as is the required time to generate synthetic labeled datasets.
MDS localization techniques have been previously researched, mostly for technologies such as Wireless Sensor Networks (WSN), Radio or 5G [22]. However, to the best of our knowledge, these techniques have not been evaluated in emerging techniques, such as ultrasound for airborne applications. Because of this, the aim of this work is the usage of this algorithm for ultrasound data for target localization.
This work is structured as follows: Section 2 presents the objectives to cover in this work. Section 3 explains the proposed new algorithms to perform both the infrastructure and target positioning. Section 3.4 explains the filter techniques studied in this work for smoothing out the trajectory, and Section 4 describes the simulation performed. Section 5 summarizes the results obtained, focusing on different parameters of each algorithm, and methods to improve the results via filtering or changing the infrastructure layout. Finally, Section 6 presents the conclusions of this work.

2. Envisioned System

The goal of the present work is to analyze the feasibility and performance of a synthetic data generation framework based on the researched algorithms, due to its capabilities to accurately generate numerical samples. The required input for the data generation is an initial selection of the followed path (equation or time series of the desired movement). At the same time, this framework would enable the user to generate a more varied dataset since the noise level can be controlled as well as different modifications of the initial data (including rotating, scaling and translating the samples) in the 3D axis, which can later be converted to different formats to fit the specific application, i.e., images or voxels.
This framework could ease data gathering tasks, as real sensors are not required for this process, and it can generate numerous relevant samples that emulate different scenarios/technologies based on the configuration selected by the user, such as the anchor distribution and noise levels.
This system would be beneficial for tasks, such as gesture recognition based on multiple technologies, which numerous authors are researching. Most of the studies in this field are focused on radar [15,23], Wi-Fi [24] and ultrasound sensors [9,13,25]. In this paper, the framework will be evaluated for the generation of ultrasound data for gesture recognition. This technology is selected due to the emerging techniques with ultrasound sensors, which could be implemented directly on simple microcontroller-based devices, like that proposed in [25].
The system to be emulated with the proposed framework is assumed to perform the following tasks (Figure 1):
  • Distance estimation. The devices use ultrasound transceiver(s) to locally compute their distances to an object, e.g., the user’s hand, typically using ToF-based measurements. The pairwise distances between the anchors are also computed (with a lower frequency) to self-locate the anchor infrastructure.
  • Positioning algorithms. Using the pairwise distances between the anchors obtained in the previous point, the position of each anchor is computed. Then, using these positions and the distances between the user’s hand and all the anchors, the current position of the object is computed.
  • Tracking algorithms. The position of the object is periodically updated, effectively obtaining an estimation of its trajectory. This trajectory is filtered to improve its accuracy.
  • Recognition. The estimated trajectory is used as input for a gesture recognition stage, e.g., implemented with a neural network.
The current work focuses on the second and third steps, in which we transform from a temporal series of distances to the 3D trajectory of the object and the 3D positions of the anchors. It is important to say that, even when in this paper, the localization algorithms are tested with synthetic data, the proposed algorithms could also be deployed in a real scenario for target positioning.
To evaluate the applicability of the proposed algorithms, the following criteria are used:
  • The computational requirements of the positioning and tracking algorithms must be low enough to be executed in real time on low-power devices. Furthermore, they must be flexible enough to adapt to time-varying and noisy conditions, with a potentially variable number of anchors in range.
  • Analyze the accuracy of both the estimated object’s trajectory and the anchor’s position. The precision of the measured data will directly affect the results when using a classification algorithm to study the data. Because of this, it is important to ensure the high performance of the localization algorithms as well as the proposed filtering techniques. To evaluate this, noise—as typically encountered in ultrasounds systems in this case—is added to the raw distances. The positioning and tracking algorithms must provide optimal estimations and a robust behavior in the presence of noise, missing distances and outliers.

3. Infrastructure and Object Positioning Using Pairwise Distances

As presented in Section 1, classical trilateration techniques calculate the position of an object based on the measurement of its distance from different reference points, or anchors, which define the location infrastructure. The algorithms used for infrastructure and object positioning using pairwise distances will be described in this section, presenting the novel two-step approach proposed in this work.

3.1. Infrastructure Positioning Using Multidimensional Scaling

For many applications, in which the infrastructure may be portable and flexible, or a quick and seamless deployment is desired, the positions of the anchors may not be known beforehand.
A self-positioning infrastructure, in which relative coordinates of the anchors are obtained from their pairwise distances, can be achieved by using metric multidimensional scaling (mMDS) techniques. All the pairwise Euclidean distances among the anchors shape the dissimilarity matrix, which is then used to calculate the relative coordinates of the anchors by minimizing a stress function based on iterative metric-preserving techniques [26]. The scaling by majorizing a complicated function (SMACOF) algorithm is proposed for a computationally-efficient resolution of the problem [27]. Distances, or equivalently dissimilarities, are considered noisy, and some pairs may be missing. To account for different degrees of confidence in the dissimilarities, they can be weighted differently, e.g., from zero (distance is considered missing and ignored during the stress computation) to one. Due to the noisy nature of real-world distance measurements, an analytical exact solution is usually not available, and iterative techniques, such as SMACOF, are more suitable.
Different approaches based on MDS are proposed in the literature to improve the computation of the coordinates [28], such as matrix completion, in which missing distances are estimated (e.g., with Dijkstra’s shortest path algorithm) instead of being given zero weight. Another is out-of-sample MDS [29,30], in which the position of a subset of anchors (landmarks) can be fixed and only the remaining positions are computed. Furthermore, mixing different steps of non-metric [31] and metric MDS computations can be beneficial, particularly if the dissimilarities are not directly proportional to the Euclidean distances (e.g., they are based on RSSI measurements [32]). The initialization of the SMACOF algorithm may also impact the accuracy of the solution, and it is usual to run it with multiple random initializations and keep the solution with lower stress [33].
The result of the MDS algorithm is a cloud of points, one for every anchor. There is a set of transformations (translation, rotation and reflection) that can be applied to these points, without modifying the stress, which results in an infinite amount of equally valid solutions. The last step is to use physical constraints or general knowledge about how and where the anchors are deployed to apply these transformations and fix the coordinate system to go from relative to absolute positions [34].

3.2. Object Positioning Using Multidimensional Scaling

Once the coordinates of the anchors are computed, the second step is the addition of the pairwise distances of the moving objects to the dissimilarity matrix, ideally using the out-of-sample variation of the MDS algorithm, in order to keep the anchor positions fixed [29]. Nevertheless, MDS is computationally expensive, and while it might be the optimal solution for low-frequency infrastructure positioning, it could be too slow for high-frequency positioning of moving objects, particularly in edge computing and low power environments.
As an alternative, the position of the fixed anchors and the moving objects can be simultaneously computed. The main disadvantage of this approach is that normally the distance measurements between the fixed anchors are more reliable than the measurements between the anchors and the moving objects, e.g., distances between the anchors can be heavily averaged for noise reduction. Introducing noisier dissimilarities in the matrix affects the accuracy of the overall positioning, including that of the infrastructure, resulting in worse results than with the two-step (first MDS without mobile objects, then out-of-sample MDS) approach. Furthermore, the fixed infrastructure does not need to be re-positioned as fast as the objects, so it is useful to decouple both computations.

3.3. Object Positioning Using Optimization Algorithms

A more efficient approach for object positioning is to compute only the coordinates of the moving objects at a faster rate, using a classical optimization algorithm, based on the anchor coordinates previously obtained with MDS. In our case, we choose the Limited-Memory Broyden Fletcher Goldfarb Shanno (LM-BFGS) algorithm [35,36] with the mean squared error as the objective function to be minimized. The processing and memory requirements of LM-BFGS are low-enough to be run in real-time in low-power edge devices with a typical number of distances to the anchors (fewer than a dozen) [37].
The accuracy of this approach is comparable to the out-of-sample MDS one since the positions of the anchors are considered fixed during the optimization iterations, but it is faster for independently obtaining the coordinates of individual objects. Furthermore, the error is typically lower than when using the one-step MDS approach, in which the coordinates of the anchors and mobile objects are computed at the same time since the pairwise distances involving the mobile objects are normally noisier.

3.4. Trajectory Optimization Using Smoothing Techniques

The coordinates of the moving object create a (typically noisy) trajectory that benefits from proper filtering in order to provide a more accurate estimation for the final application, which could allow real-time localization or hand-gesture recognition. Different techniques are widely used for low-pass filtering and outlier detection. In our approach, we compare simple moving-average and moving-median filters [38] with a fixed window length determined heuristically. They provide optimal results, while keeping both the computational cost and low complexity.
To obtain the trajectory of a moving object, we compute its position periodically, building a discrete time-series of successive equally spaced points in time. Since additional information about the expected path is usually known, such as the maximum velocity of the object, we can exploit this to further filter the trajectory, smoothing it out, restoring missing points by interpolation, and decreasing or removing the effect of outliers in the position time-series. We explore two filtering alternatives:
  • Moving-average filter. It is a low-pass filter that provides effective noise reduction, particularly in applications where the focus is on time-response (instead of frequency-response) analysis. It smooths the signal, but the predicted trajectory may fail to respond to quick movements.
  • Moving-median filter. The median filter is a non-linear filter that replaces the values in data with the moving median of the filtered and neighboring points. It is very robust against outliers and in suppressing spiky noise, but as with the moving-average filter filtering, it may lead to an underestimation of the path, particularly in sharp corners.

4. Simulation Setup

After presenting the different alternatives for object and infrastructure positioning based on pairwise distance measurements, we estimate the performance of the proposed algorithms in terms of accuracy and execution speed. We compare the two different approaches depicted in Figure 2:
  • One-step approach (Figure 2a). The SMACOF MDS algorithm is used to simultaneously obtain the positions of the anchors and the moving objects. It is expected to be slower and less accurate if noisy dissimilarities (such as those between the moving objects) are introduced in the computation, but all the anchors and object positions are computed simultaneously.
  • Two-step approach (Figure 2b). The SMACOF MDS algorithm is used once to obtain the positions of the anchors. Then, the LM-BFGS optimization algorithm is used to compute only the coordinates of the moving object, and repeated periodically to update its position. This approach is faster, but relies on an accurate initial estimation of the anchor positions.
By using these two algorithms, it is possible to design the proposed framework for synthetic data generation. It executes the SMACOF MDS algorithm periodically to ensure that the position of the anchors is correct while locating simultaneously the target. Between these anchors check, the LM-BFGS algorithm is used due to its low latency and high accuracy when the position of the anchors is known.This approach is faster and results in a smaller error, as we will discuss in Section 5.2 and Section 5.4. The SMACOF MDS and LM-BFGS optimization computation steps can be done in a central processing node to which all the anchors report, or it can be done locally in the anchors or the mobile object, if they have access to all the distances. The particular communication scheme to disseminate the distances and the positions is out of the scope of this work. Finally, if we want to estimate a path and not only the single positions, a smoothing filter is used to compute the trajectory of the object.
Consequently, this framework can be used to generate synthetic trajectories for an arbitrary number of anchor configurations and gestures to fit multiple scenarios and applications, as shown in the synthetic data creation block in Figure 2. At the same time, data augmentation for a single gesture and anchors setup is possible by varying the random initialization seeds of the noise for the SMACOF MDS and the LM-BFGS optimization algorithms, as shown in the data estimation block in Figure 2. Consequently, this framework can efficiently generate numerous samples of the desired data to contemplate all the possible results of measurements with real devices. Furthermore, different noise models and strengths can be injected to the raw distances, emulating different disturbances and inaccuracies that the ultrasound distance gathering system can experience in a real deployment.

4.1. System Modeling

In this work, the framework emulates a system of nine ultrasound anchors sending a sinusoidal pulse (reference pulse). The echoes are then sampled, and the ToF is obtained with classical cross-correlation techniques, using a reference signal. Every pairwise distance is computed at 20 Hz, i.e., a new target position is computed every 50 ms. A two-step approach to compute the anchor and target positions as depicted in Figure 2.
The infrastructure of the anchors for the proposed framework is shaped as a 2D array of nine anchors (ultrasound transceivers) located in the same surface (in the XY plane with z = 0 ), which represents a plausible configuration for future applications. Specifically, the anchors are located as seen in Figure 3. The positioning is limited to the space in front of said surface ( z > 0 ) since the sign of the z coordinate cannot be defined when all the anchors are in the same plane. Furthermore, the ultrasound transceivers sensors used as an experimental support for the simulations have a detection range limited to 180 degrees in the Z-axis.
The anchor array is able to transmit and receive ultrasonic signals and locate passive objects based on ToF measurements. It has two operating modes:
  • To calculate the pairwise distances between the anchors, they actively exchange ultrasonic signals (two-way ranging).
  • To calculate the pairwise distances between the anchors and the mobile object, they actively transmit and then sense the reflected echo. Anchors can be synchronized, in which case, only one of the transceivers needs to transmit and they all can receive the echo and timestamp it based on a common clock. Otherwise, they can all transmit and sense only the echo coming from their own transmission; in such a case, time synchronization is not required.

4.2. Noise Modeling

The noise in the distances obtained with an ultrasound-based measurement system depends on the accuracy of the ToF samples. There are different factors that impact the performance, such as the bandwidth of the transmitted pulse and the sampling rate of the acquisition stage.
Based on our experimental measurements using the system of Figure 1, the noise, N, in the computed distances, d ( t ) , can be modeled as unbiased (zero average) additive white Gaussian noise (AWGN), with a given standard deviation, σ , and a probability density function, p d f ( N ) , given by the following:
d ( t ) = d r e a l ( t ) + N ( t )
p d f ( N ) = 1 σ 2 π exp 1 2 x σ 2
In a representative x-y-z point, P = ( 0 cm , 0 cm , 50 cm ) , the measured equivalent noise in the Euclidean distance, after acquiring 50,000 samples, can be fitted with a σ = 3.2 mm, as seen in Figure 4. This provides a good estimation of the scale of the expected noise in a real system, and it is used as a reference to model the noise in the simulations.
As a summary, in this section, the two-steps algorithm proposed and evaluated in this work is described, as well as the steps performed to particularize the framework for validate the use of ultrasonic system as technology for that algorithm.

5. Performance Results

In this section, the performance of the proposed algorithms for the object position and anchors localization is analyzed. At the same time, the utility of the proposed system for data generation tasks is discussed.

5.1. Infrastructure-Positioning Accuracy

First, we characterize the performance of SMACOF MDS locating the infrastructure by building the dissimilarity matrix, D, with the metric distances measured between the anchors, adding different realistic levels of AWGN noise, according to the experimental results obtained in Section 4.2.
Since the positions computed by SMACOF MDS are equally valid if they are translated, rotated or mirrored, we need additional constraints to fix the coordinate system. We use an approach to agree on a common absolute reference that requires one translation and three rotations:
  • We translate the points so that the first anchor, m 0 , fixes the origin of the coordinates system, P m 0 = ( 0 , 0 , 0 ) .
  • We rotate the points around the X axis at an angle such that the second anchor, m 1 , is in z = 0 , P m 1 = ( a , b , 0 ) .
  • We rotate the points around the X axis at an angle such that the second anchor, m 1 , is in y = 0 and c > 0 , P m 1 = ( c , 0 , 0 ) .
  • We rotate the points around the X axis at an angle such that the third anchor, m 2 , is in z = 0 and e > 0 , P m 2 = ( d , e , 0 ) .
  • In our configuration, all the anchors are in the same surface ( z = 0 ), so no additional condition is required. In general, it is required to fix the positive direction of the X-axis. If the anchors constitute a three-dimensional shape, another anchor (e.g., m 3 ), located in a surface different to the previous three ( m 0 to m 2 ), is used to define the positive X direction.
Following these steps, we obtain a consistent reference system every time, removing any ambiguity in the positions. The results locating the anchor infrastructure depending on the strength, σ , of the AWGN noise in the distances are shown in Figure 5. The error on the X and Y axes between the real and calculated positions can be observed. Table 1 provides the quantitative comparison for this error in the X, Y and Z axes, and also the total displacement. We obtain an optimal accuracy, with an expected error comparable to the deviation introduced by the noise strength.

5.2. Object-Positioning Accuracy

Once the anchor infrastructure is positioned, either with the MDS-based self-locating mechanisms explained in Section 5.1, or because the positions of the anchors are previously known, we can either apply LM-BFGS optimization to locate the mobile objects or reapply SMACOF MDS, adding the distance of the mobile object to the dissimilarity matrix. For the simulations, we use the anchor infrastructure of Figure 3, and repeat the positioning of 250,000 objects randomly distributed in the surface given by 0 < z < 1000 mm; 500 mm < x < 500 mm and y = 0 . As shown in Figure 6 and Table 2, since the anchor infrastructure is symmetrical with respect to the z = 0 surface, the averaged error in the X and Y axes is similar, while it is generally much larger for the Z axis. The chosen anchor infrastructure is flat, with no variability in the X axis. Such a configuration works well for locating in the X and Y axes, but it struggles in the Z axis. The error also increases as we move away from the anchors, being minimal in the region close to them. LM-BFGS optimization, besides being faster, has fewer errors since it only tries to optimize the position of the mobile object, without trying to minimize the error by also relocating the fixed anchors. The results further validate our simple two-step approach of using SMACOF MDS just for positioning the infrastructure and then employing LM-BFGS optimization to update the position of the mobile object.

5.3. Trajectory-Positioning Accuracy

As a relevant application of proposed system for synthetic data generation, we simulated the data acquisition process based on ultrasound transceivers for gesture recognition. This framework was evaluated by creating three different gestures based on some initial time series of the 3D coordinates of the trajectory.
The chosen gestures are presented in Figure 7. They simulate the imperfections of real gestures. The first one (Figure 7a, c i r c l e X Y ) is a circular shape in a surface with low X axis variability ( 100 < z < 120 ), and mostly parallel to the anchor infrastructure, to simulate a gesture in the optimal recognizing conditions. The second one (Figure 7b, l o o p X Y ) is a gesture with sharp corners, in the same surface, to evaluate the performance of the smoothing filters in gestures with sharp edges. The third one (Figure 7c, c i r c l e X Z ) is another circular shape, this time located in a surface mostly perpendicular to the anchor infrastructure ( 100 < y < 120 ), to check the worst-case performance (strong z variability). The first gesture has 400 samples, while the second and third have 250 samples each. They are all sampled at 20 Hz.
First of all, the frame tests different filter orders, M, to heuristically choose a value in which the error is properly minimized. In Figure 8, the error dependency with the filter order is shown, for both filters described (moving-median and moving-average) and for both positioning algorithms (SMACOF MDS and LM-BFGS optimization). The error without filtering is added too, as a reference. Based on those results, M = 11 is selected as a good trade-off between added delay and error suppression. With the considered AWGN noise model, the moving-average filtering performs consistently better than the moving-median filter, which may change if spiky noise and outliers are introduced in the model.
Without added noise ( σ = 0 mm), we can see in Table 3 that the smoothing filters actually deteriorate the performance since the edges are underestimated. In a realistic noisy scenario, the smoothing filter greatly reduces the error, as shown in Table 4, Table 5 and Table 6, where the average error (in each axis and in total) is compared for each algorithm, and each filtering technique previously described. As seen in Figure 9, as expected, the noisiest gesture in the surface of interest is the c i r c l e X Z , due to it being performed in a plane normal to the anchor surface. The smoothing filtering also struggles in the sharp edges of the l o o p X Y gesture, while it estimates the path with very good accuracy for the inherently smooth circular gestures.
Nevertheless, the achieved results show how the proposed algorithms are able to accurately generate the desired gesture/trajectory data with different noise levels. This would enable the framework to simulate a large range of possible scenarios and sensors. The noise level can be defined by the user to adapt the specific application and scenario that the framework is emulating.

5.4. Execution Time

We use the implementation of the SMACOF MDS and LM-BFGS optimization algorithms included in the Python library for machine learning scikit-learn 0.23, and run them in a typical Windows laptop (Intel Core [email protected]., 8GB RAM). Since absolute values depend on the processing capabilities of the particular machine, only relative time differences are evaluated. In this case, the evaluation is performed over the gesture represented in Figure 7a. As seen in Table 7, the SMACOF MDS algorithm is about ten times slower, which might be relevant for real-time applications, particularly if the position inference needs to be done directly in an embedded edge device, which typically has constrained resources. These constrains could lead to a larger difference between these algorithms, due to the computing power required to execute the SMACOF MDS algorithm. In the context of continuously positioning a slowly (relative to the positioning sampling rate) moving object, initializing the algorithm with the previous position greatly reduces (it takes half the time) the execution time of SMACOF MDS, while for LM-BFGS optimization, it barely changes.
Although SMACOF MDS is slower, it needs to be stressed that its output includes the position of all the anchors, while LM-BFGS optimization only computes the position of a single moving object. Therefore, SMACOF MDS is relevant for self-calibrating anchor infrastructures or by simultaneously locating multiple moving objects. Once the infrastructure is properly positioned during the initial set-up phase, LM-BFGS optimization can be used for quick updates of the object position.

6. Conclusions

This work presents a novel two-step technique to perform general infrastructure and moving-object positioning based on measured pairwise distances. In the first step, MDS is used to obtain the coordinates of the anchors, repeated with a low frequency, e.g., to correct minor and infrequent potential displacements of the anchors. We use the SMACOF variant of the mMDS family of algorithms. With the coordinates of the computed anchors computed, a fast optimization algorithm is used to obtain the unknown coordinates of the objects. This step is repeated with a high frequency. The LM-BFGS optimization algorithm is used for this step. Its performance is thoroughly analyzed with simulations, particularized to the use case of a system with ultrasound transceivers. The distribution and shape of the anchor infrastructure, the size of the region in which the positioning takes place and the strength of the noise are realistically modeled after such a system.
This two-step approach described in the work would be optimal in scenarios where the position of the anchors does not change frequently through time. Therefore, the one-step approach described in Section 4 in which all the positions are computed at the same time, is reserved for special situations, e.g., when there are no anchors (all the objects are considered mobile) or when we want to simultaneously obtain the position of several (more than a dozen) mobile objects. For the rest of the scenarios, our approach performs the localization with low computational time, making it suitable for use in real-time systems and even in constrained edge devices.
Efficient and simple filtering techniques significantly reduce the error and improve the reconstruction of the real path followed by the mobile object. This feature can be exploited when using the proposed algorithms for synthetic data generation. The current dataset creation step for applications, such as AI models, are time consuming, due to the complexity of the recording and labeling tasks, which could be reduced by using the proposed system as a synthetic data generation framework. This framework is independent from hardware and it could simulate trajectories/movement from a large range of sensors. The parameters of this framework (noise, gesture and anchors number and position) are defined by the user through the initial configurations.
The use of ultrasonic signals for target positioning is widely researched, but to the best of our knowledge, our two-step approach inspired by wireless sensor network’s positioning algorithms has not been used or described. The proposed technique enables using an arbitrary number of ultrasound transceivers, and removes the constraint of knowing the position of the anchors beforehand, while providing an optimal AWGN rejection. This could drive the adoption of ultrasound technology in the positioning field and foster the research of novel applications and electronic components based on non-audible acoustic waves.

Author Contributions

Conceptualization, B.S.-M.; methodology, B.S.-M. and A.E.-M.; validation, B.S.-M., A.E.-M. and J.M.-G.; formal analysis, E.C.-M. and D.P.M.-S.; investigation, B.S.-M.; writing—original draft preparation, B.S.-M., A.E.-M. and J.M.-G.; writing—review and editing, E.C.-M. and D.P.M.-S.; supervision, E.C.-M. and D.P.M.-S.; project administration, E.C.-M. and D.P.M.-S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Part of this work has been developed in the context of the project “SEMULIN” (German project number 19A20012D) and financed by the German Federal Ministry for Economic Affairs and Energy (BMWi), and the project P20-00265 funded by the “Junta de Andalucia” of Spain.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Leonard, J.J.; Durrant-Whyte, H.F. Mobile robot localization by tracking geometric beacons. IEEE Trans. Robot. Autom. 1991, 7, 376–382. [Google Scholar] [CrossRef]
  2. Smith, A.; Balakrishnan, H.; Goraczko, M.; Priyantha, N. Tracking moving devices with the cricket location system. In Proceedings of the 2nd International Conference on Mobile Systems, Applications, and Services, Boston, MA, USA, 6 June 2004; pp. 190–202. [Google Scholar]
  3. Hand, C. A survey of 3D interaction techniques. In Computer Graphics Forum; Wiley Online Library: Hoboken, NJ, USA, 1997; Volume 16, pp. 269–281. [Google Scholar]
  4. Wren, C.R.; Azarbayejani, A.; Darrell, T.; Pentland, A.P. Pfinder: Real-time tracking of the human body. IEEE Trans. Pattern Anal. Mach. Intell. 1997, 19, 780–785. [Google Scholar] [CrossRef] [Green Version]
  5. Gobbi, D.G.; Comeau, R.M.; Peters, T.M. Ultrasound probe tracking for real-time ultrasound/MRI overlay and visualization of brain shift. In International Conference on Medical Image Computing and Computer-Assisted Intervention; Springer: Berlin/Heidelberg, Germany, 1999; pp. 920–927. [Google Scholar]
  6. Qi, J.; Liu, G.P. A robust high-accuracy ultrasound indoor positioning system based on a wireless sensor network. Sensors 2017, 17, 2554. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Medina, C.; Segura, J.C.; De la Torre, A. Ultrasound indoor positioning system based on a low-power wireless sensor network providing sub-centimeter accuracy. Sensors 2013, 13, 3501–3526. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Povolny, A.; Kikura, H.; Ihara, T. Ultrasound pulse-echo coupled with a tracking technique for simultaneous measurement of multiple bubbles. Sensors 2018, 18, 1327. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Chen, H.; Ballal, T.; Muqaibel, A.H.; Zhang, X.; Al-Naffouri, T.Y. Air Writing via Receiver Array-Based Ultrasonic Source Localization. IEEE Trans. Instrum. Meas. 2020, 69, 8088–8101. [Google Scholar] [CrossRef]
  10. Cobos, M.; Antonacci, F.; Alexandridis, A.; Mouchtaris, A.; Lee, B. A survey of sound source localization methods in wireless acoustic sensor networks. Wirel. Commun. Mob. Comput. 2017, 2017, 3956282. [Google Scholar] [CrossRef]
  11. Putz, V.; Mayer, J.; Fenzl, H.; Schmidt, R.; Pichler-Scheder, M.; Kastl, C. Cyber–Physical Mobile Arm Gesture Recognition using Ultrasound and Motion Data. In Proceedings of the 3rd IEEE International Conference on Industrial Cyber—Physical Systems, Tampere, Finland, 10–12 June 2020. [Google Scholar]
  12. Laadung, T.; Ulp, S.; Alam, M.M.; Le Moullec, Y. Active-Passive Two-Way Ranging Using UWB. In Proceedings of the 2020 14th International Conference on Signal Processing and Communication Systems (ICSPCS), Adelaide, Australia, 14–16 December 2020; pp. 1–5. [Google Scholar] [CrossRef]
  13. Chen, J.; Yu, F.; Yu, J.; Lin, L. A Three-Dimensional Pen-Like Ultrasonic Positioning System Based on Quasi-Spherical PVDF Ultrasonic Transmitter. IEEE Sensors J. 2020, 21, 1756–1763. [Google Scholar] [CrossRef]
  14. Yiallourides, C.; Parada, P.P. Low Power Ultrasonic Gesture Recognition for Mobile Handsets. In Proceedings of the ICASSP 2019-2019 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Brighton, UK, 12–17 May 2019; pp. 2697–2701. [Google Scholar] [CrossRef]
  15. Arsalan, M.; Santra, A. Character recognition in air-writing based on network of radars for human-machine interface. IEEE Sensors J. 2019, 19, 8855–8864. [Google Scholar] [CrossRef]
  16. Jackson, J.C.; Summan, R.; Dobie, G.I.; Whiteley, S.M.; Pierce, S.G.; Hayward, G. Time-of-flight measurement techniques for airborne ultrasonic ranging. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2013, 60, 343–355. [Google Scholar] [CrossRef] [PubMed]
  17. Lazik, P.; Rowe, A. Indoor pseudo-ranging of mobile devices using ultrasonic chirps. In Proceedings of the 10th ACM Conference on Embedded Network Sensor Systems, Toronto, ON, Canada, 6 –9 November 2012; pp. 99–112. [Google Scholar]
  18. Lindo, A.; Garcia, E.; Urena, J.; del Carmen, M.; Hernandez, A. Multiband waveform design for an ultrasonic indoor positioning system. IEEE Sensors J. 2015, 15, 7190–7199. [Google Scholar] [CrossRef]
  19. Martın, J.M.; Jiménez, A.R.; Seco, F.; Calderón, L.; Pons, J.L.; Ceres, R. Estimating the 3D-Position from Time Delay Data of US-Waves: Experimental Analysis and a New Processing Algorithm; Elsevier: Amsterdam, The Netherlands, 2002; Volume 101, pp. 311–321. [Google Scholar]
  20. Ruiz, D.; Ureña, J.; García, J.C.; Pérez, C.; Villadangos, J.M.; García, E. Efficient Trilateration Algorithm Using Time Differences of Arrival; Elsevier: Amsterdam, The Netherlands, 2013; Volume 193, pp. 220–232. [Google Scholar]
  21. Norrdine, A. An algebraic solution to the multilateration problem. In Proceedings of the 15th International Conference on Indoor Positioning and Indoor Navigation, Sydney, Australia, 13–15 November 2012; Volume 1315. [Google Scholar]
  22. Saeed, N.; Nam, H.; Al-Naffouri, T.Y.; Alouini, M.S. A state-of-the-art survey on multidimensional scaling-based localization techniques. IEEE Commun. Surv. Tutorials 2019, 21, 3565–3583. [Google Scholar] [CrossRef] [Green Version]
  23. Leem, S.K.; Khan, F.; Cho, S.H. Detecting mid-air gestures for digit writing with radio sensors and a CNN. IEEE Trans. Instrum. Meas. 2019, 69, 1066–1081. [Google Scholar] [CrossRef]
  24. Fang, Y.; Xu, Y.; Li, H.; He, X.; Kang, L. Writing in the air: Recognize Letters Using Deep Learning Through WiFi Signals. In Proceedings of the 2020 6th International Conference on Big Data Computing and Communications (BIGCOM), Deqing, China, 24–25 July 2020; pp. 8–14. [Google Scholar] [CrossRef]
  25. Saez, B.; Mendez, J.; Molina, M.; Castillo, E.; Pegalajar, M.; Morales, D.P. Gesture Recognition with Ultrasounds and Edge Computing. IEEE Access 2021, 9, 38999–39008. [Google Scholar] [CrossRef]
  26. Abdi, H. Metric multidimensional scaling (MDS): Analyzing distance matrices. Encycl. Meas. Stat. 2007, 1–13. [Google Scholar]
  27. De Leeuw, J.; Mair, P. Multidimensional scaling using majorization: SMACOF in R. J. Stat. Softw. 2009, 31, 1–30. [Google Scholar] [CrossRef] [Green Version]
  28. Saeed, N.; Nam, H.; Haq, M.I.U.; Muhammad Saqib, D.B. A Survey on Multidimensional Scaling; ACM: New York, NY, USA, 2018; pp. 1–25. [Google Scholar]
  29. Trosset, M.W.; Priebe, C.E. The Out-of-Sample Problem for Classical Multidimensional Scaling; Elsevier: Amsterdam, The Netherlands, 2008; Volume 52, pp. 4635–4642. [Google Scholar]
  30. De Silva, V.; Tenenbaum, J.B. Sparse Multidimensional Scaling Using Landmark Points; Technical Report; Stanford University: Stanford, CA, USA, 2004. [Google Scholar]
  31. Xiao, L.; Li, R.; Luo, J. Sensor localization based on nonmetric multidimensional scaling. STRESS 2006, 2, 3225–3230. [Google Scholar]
  32. Wang, S. Wireless Network Indoor Positioning Method Using Nonmetric Multidimensional Scaling and RSSI in the Internet of Things Environment. Hindawi 2020, 2020, 8830891. [Google Scholar] [CrossRef]
  33. Patanè, M.; Rossi, B.; Fragneto, P.; Fusiello, A. Wireless Sensor Networks localization with outliers and structured missing data. In Proceedings of the 2017 IEEE 28th Annual International Symposium on Personal, Indoor, and Mobile Radio Communications (PIMRC), Montreal, QC, Canada, 8–13 October 2017; pp. 1–7. [Google Scholar]
  34. Li, X. Collaborative localization with received-signal strength in wireless sensor networks. IEEE Trans. Veh. Technol. 2007, 56, 3807–3817. [Google Scholar]
  35. Liu, D.C.; Nocedal, J. On the Limited Memory BFGS Method for Large Scale Optimization; Springer: Berlin/Heidelberg, Germany, 1989; Volume 45, pp. 503–528. [Google Scholar]
  36. Morales, J.L. A Numerical Study of Limited Memory BFGS Methods; Elsevier: Amsterdam, The Netherlands, 2002; Volume 15, pp. 481–487. [Google Scholar]
  37. Mokhtari, A.; Ribeiro, A. Global convergence of online limited memory BFGS. J. Mach. Learn. Res. 2015, 16, 3151–3181. [Google Scholar]
  38. Singhal, D.; Garimella, R.M. Simple Median based information fusion in wireless sensor network. In Proceedings of the 2012 International Conference on Computer Communication and Informatics, Chennai, India, 3–5 August 2012; pp. 1–7. [Google Scholar]
Figure 1. Ultrasound positioning system used as a reference for the simulations.
Figure 1. Ultrasound positioning system used as a reference for the simulations.
Sensors 21 05923 g001
Figure 2. Proposed simulation framework, with (a) One-step approach and (b) Two-step approach, for the generation of synthetic gesture data.
Figure 2. Proposed simulation framework, with (a) One-step approach and (b) Two-step approach, for the generation of synthetic gesture data.
Sensors 21 05923 g002
Figure 3. Position of the ultrasound anchors.
Figure 3. Position of the ultrasound anchors.
Sensors 21 05923 g003
Figure 4. Histogram of the measured noise in experimentally-computed distances and normal distribution fit ( σ = 3.2 mm).
Figure 4. Histogram of the measured noise in experimentally-computed distances and normal distribution fit ( σ = 3.2 mm).
Sensors 21 05923 g004
Figure 5. Simulated infrastructure-positioning results for different noise levels using SMACOF MDS (MDS).
Figure 5. Simulated infrastructure-positioning results for different noise levels using SMACOF MDS (MDS).
Sensors 21 05923 g005
Figure 6. Average object-positioning error (Euclidean distance between the real and the computed point) in different regions and for different noise levels, with SMACOF MDS (MDS) and LM-BFGS optimization (OPT).
Figure 6. Average object-positioning error (Euclidean distance between the real and the computed point) in different regions and for different noise levels, with SMACOF MDS (MDS) and LM-BFGS optimization (OPT).
Sensors 21 05923 g006
Figure 7. Gestures used in the simulations.
Figure 7. Gestures used in the simulations.
Sensors 21 05923 g007
Figure 8. Comparison of the average trajectory error for the cilclexy gesture, with a fixed noise level (σ = 10 mm), with moving-median (Med.), moving-average (Avg.) smoothing filter and without (Org.) smoothing, for SMACOF MDS (MDS) and LM-BFGS optimization (OPT) algorithms.
Figure 8. Comparison of the average trajectory error for the cilclexy gesture, with a fixed noise level (σ = 10 mm), with moving-median (Med.), moving-average (Avg.) smoothing filter and without (Org.) smoothing, for SMACOF MDS (MDS) and LM-BFGS optimization (OPT) algorithms.
Sensors 21 05923 g008
Figure 9. Comparison between the real and estimated trajectories for SMACOF MDS (MDS) and LM-BFGS optimization (OPT) algorithms for the three different gestures, with a fixed noise level (σ = 10 mm), without smoothing (sub-figures ac), with moving-average (Avg.) smoothing (subfigures df) and with moving-median (Med.) smoothing (sub-figures gi).
Figure 9. Comparison between the real and estimated trajectories for SMACOF MDS (MDS) and LM-BFGS optimization (OPT) algorithms for the three different gestures, with a fixed noise level (σ = 10 mm), without smoothing (sub-figures ac), with moving-average (Avg.) smoothing (subfigures df) and with moving-median (Med.) smoothing (sub-figures gi).
Sensors 21 05923 g009
Table 1. Simulated infrastructure-positioning average error (of the nine anchors) for different noise levels with SMACOF MDS, per-axis and Euclidean distance.
Table 1. Simulated infrastructure-positioning average error (of the nine anchors) for different noise levels with SMACOF MDS, per-axis and Euclidean distance.
Average Positioning Error (mm)
Algorithm σ = 0 mm σ = 2.5 mm σ = 5 mm σ = 10 mm
X Ax.SMACOF MDS00.020.080.11
Y Ax.SMACOF MDS00.030.080.10
Z Ax.SMACOF MDS00.010.010.03
Euc. Dis.SMACOF MDS00.030.080.14
Table 2. Per-axis object-positioning error and Euclidean distance between the real and the computed point for different noise levels, with SMACOF MDS and LM-BFGS optimization, averaged for the 0 < z < 1000 mm; 500 mm < x < 500 mm and y = 0 surface.
Table 2. Per-axis object-positioning error and Euclidean distance between the real and the computed point for different noise levels, with SMACOF MDS and LM-BFGS optimization, averaged for the 0 < z < 1000 mm; 500 mm < x < 500 mm and y = 0 surface.
Average Positioning Error (mm)
Algorithm σ = 0 mm σ = 2.5 mm σ = 5 mm σ = 10 mm
X Ax.SMACOF MDS0.0316.3528.3948.21
LM-BFGS0.017.6415.2128.69
Y Ax.SMACOF MDS0.0315.6726.9845.54
LM-BFGS0.017.5314.8128.88
Z Ax.SMACOF MDS0.1736.9456.9179.65
LM-BFGS0.0423.2542.0567.50
Euc. Dis.SMACOF MDS0.1947.4475.85114.79
LM-BFGS0.0527.6850.9785.91
Table 3. Average trajectory error for the three different gestures without noise ( σ = 0 mm), with no smoothing (Original) and with 11th order ( M = 11 ) moving-average (Average) and moving-median (Median) smoothing filters, for SMACOF MDS and LM-BFGS optimization algorithms.
Table 3. Average trajectory error for the three different gestures without noise ( σ = 0 mm), with no smoothing (Original) and with 11th order ( M = 11 ) moving-average (Average) and moving-median (Median) smoothing filters, for SMACOF MDS and LM-BFGS optimization algorithms.
Average Positioning Error (mm)
c i r c l e X Y l o o p X Y c i r c l e X Z
X Ax.SMACOF MDSOriginal0.030.020.04
Average0.692.711.67
Median0.191.020.14
LM-BFGSOriginal0.010.010.01
Average0.692.721.67
Median0.171.010.11
Y Ax.SMACOF MDSOriginal0.030.030.04
Average1.41.790.29
Median0.070.670.19
LM-BFGSOriginal0.010.010.01
Average1.41.790.29
Median0.050.660.18
Z Ax.SMACOF MDSOriginal0.080.090.02
Average0.340.362.17
Median0.220.280.38
LM-BFGSOriginal0.020.020.01
Average0.340.342.17
Median0.190.260.37
Euc. Dis.SMACOF MDSOriginal0.090.10.07
Average1.763.553.07
Median0.431.490.62
LM-BFGSOriginal0.030.030.02
Average1.763.543.07
Median0.391.460.6
Table 4. Average trajectory error for the c i r c l e X Y gesture, with no smoothing (Original) and with 11th order ( M = 11 ) moving-average (Average) and moving-median (Median) smoothing filters, for SMACOF MDS and LM-BFGS optimization algorithms.
Table 4. Average trajectory error for the c i r c l e X Y gesture, with no smoothing (Original) and with 11th order ( M = 11 ) moving-average (Average) and moving-median (Median) smoothing filters, for SMACOF MDS and LM-BFGS optimization algorithms.
Average Positioning Error (mm)
σ = 2.5 mm σ = 5 mm σ = 10 mm
X Ax.SMACOF MDSOriginal14.5623.6840.88
Average5.058.3514.30
Median6.269.8415.84
LM-BFGSOriginal6.6813.2425.73
Average2.194.027.11
Median3.365.458.08
Y Ax.SMACOF MDSOriginal13.4221.7235.29
Average4.849.8119.26
Median6.239.5017.76
LM-BFGSOriginal6.4312.8223.41
Average2.544.228.04
Median3.175.4410.36
Z Ax.SMACOF MDSOriginal33.3348.2463.19
Average10.4615.5222.56
Median12.6519.1231.46
LM-BFGSOriginal18.6138.2263.7
Average7.5415.7322.16
Median7.2113.7230.55
Euc. Dis.SMACOF MDSOriginal41.9362.7291.61
Average13.7222.0236.5
Median17.0125.7943.27
LM-BFGSOriginal22.2845.3478.02
Average9.0218.0227.12
Median9.6817.3935.56
Table 5. Average trajectory error for the l o o p X Y gesture, with no smoothing (Original) and with 11th order ( M = 11 ) moving-average (Average) and moving-median (Median) smoothing filters, for SMACOF MDS and LM-BFGS optimization algorithms.
Table 5. Average trajectory error for the l o o p X Y gesture, with no smoothing (Original) and with 11th order ( M = 11 ) moving-average (Average) and moving-median (Median) smoothing filters, for SMACOF MDS and LM-BFGS optimization algorithms.
Average Positioning Error (mm)
σ = 2.5 mm σ = 5 mm σ = 10 mm
X Ax.SMACOF MDSOriginal11.4419.3333.62
Average6.288.2318.59
Median8.9512.0119.29
LM-BFGSOriginal5.3010.1821.5
Average3.264.5511.02
Median5.107.8813.77
Y Ax.SMACOF MDSOriginal11.1118.3234.4
Average4.987.2918.63
Median6.648.5921.52
LM-BFGSOriginal5.7510.6422.07
Average2.453.958.75
Median3.486.1211.34
Z Ax.SMACOF MDSOriginal29.6242.0462.45
Average11.7313.7537.2
Median15.1416.3742.33
LM-BFGSOriginal18.932.8553.58
Average6.7113.0717.92
Median8.8013.1424.65
Euc. Dis.SMACOF MDSOriginal36.0953.8585.89
Average15.4819.3348.43
Median20.7724.955.18
LM-BFGSOriginal21.9238.8266.6
Average8.9615.6525.77
Median12.1918.9933.85
Table 6. Average trajectory error for the c i r c l e X Z gesture, with no smoothing (Original) and with 11th order ( M = 11 ) moving-average (Average) and moving-median (Median) smoothing filters, for SMACOF MDS and LM-BFGS optimization algorithms.
Table 6. Average trajectory error for the c i r c l e X Z gesture, with no smoothing (Original) and with 11th order ( M = 11 ) moving-average (Average) and moving-median (Median) smoothing filters, for SMACOF MDS and LM-BFGS optimization algorithms.
Average Positioning Error (mm)
σ = 2.5 mm σ = 5 mm σ = 10 mm
X Ax.SMACOF MDSOriginal32.4665.271.5
Average8.9819.1232.04
Median12.2525.0435.4
LM-BFGSOriginal11.7425.6147.97
Average3.977.8016.2
Median5.319.4517.85
Y Ax.SMACOF MDSOriginal34.4663.7688.26
Average12.1619.2533.57
Median13.7623.634.62
LM-BFGSOriginal11.0123.1449.88
Average3.948.7615.84
Median4.8411.5522.7
Z Ax.SMACOF MDSOriginal12.8526.4647.18
Average5.5914.1433.14
Median7.3612.9632.24
LM-BFGSOriginal4.719.9319.13
Average2.614.079.21
Median3.915.809.80
Euc. Dis.SMACOF MDSOriginal53.97105.00158.170
Average18.5333.7262.72
Median23.2641.2766.34
LM-BFGSOriginal18.6240.1380.59
Average7.0914.3327.79
Median9.618.4734.19
Table 7. Average computation time of individual positions after 1000 executions of the gesture of Figure 7a using SMACOF MDS and LM-BFGS optimization algorithms.
Table 7. Average computation time of individual positions after 1000 executions of the gesture of Figure 7a using SMACOF MDS and LM-BFGS optimization algorithms.
Time (ms)
SMACOF MDSRandom initialization162
Previous-point initialization77
LM-BFGSRandom initialization16
Previous-point initialization15
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Saez-Mingorance, B.; Escobar-Molero, A.; Mendez-Gomez, J.; Castillo-Morales, E.; Morales-Santos, D.P. Object Positioning Algorithm Based on Multidimensional Scaling and Optimization for Synthetic Gesture Data Generation. Sensors 2021, 21, 5923. https://doi.org/10.3390/s21175923

AMA Style

Saez-Mingorance B, Escobar-Molero A, Mendez-Gomez J, Castillo-Morales E, Morales-Santos DP. Object Positioning Algorithm Based on Multidimensional Scaling and Optimization for Synthetic Gesture Data Generation. Sensors. 2021; 21(17):5923. https://doi.org/10.3390/s21175923

Chicago/Turabian Style

Saez-Mingorance, Borja, Antonio Escobar-Molero, Javier Mendez-Gomez, Encarnacion Castillo-Morales, and Diego P. Morales-Santos. 2021. "Object Positioning Algorithm Based on Multidimensional Scaling and Optimization for Synthetic Gesture Data Generation" Sensors 21, no. 17: 5923. https://doi.org/10.3390/s21175923

APA Style

Saez-Mingorance, B., Escobar-Molero, A., Mendez-Gomez, J., Castillo-Morales, E., & Morales-Santos, D. P. (2021). Object Positioning Algorithm Based on Multidimensional Scaling and Optimization for Synthetic Gesture Data Generation. Sensors, 21(17), 5923. https://doi.org/10.3390/s21175923

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop