Abstract
Multidimensional projection techniques provide graphical representations computed based on instance similarities to enable the analysis of abstract and possibly large data sets. However, when the data set size grows these graphical representations can hardly avoid overlap among markers. To overcome this issue, while some techniques attempt to remove overlap after multidimensional projection, some projection techniques were developed considering non-overlapping constraints. In this work, we present an analysis of four overlap removal techniques and two projection techniques considering non-overlapping. The evaluation was performed according to five metrics that consider structural and similarity relations, and based on the results we provide a guide to use a technique according to the data set and analysis goals.
Keywords
Introduction
In the last years, much effort has been deposited towards improving data exploration due to the amount and complexity of data generated daily. One of the approaches used to aid in the decision-making process is the integration of different areas of knowledge and a visual approach. This process is called visual analytics. 1 In this context, visualization uses visual metaphors to provide insights about data sets. Such visual metaphors can preserve and enhance the characteristics and relationship among instances or attributes.
For multidimensional data, one class of technique widely explored by the visualization community is multidimensional projections. A multidimensional projection2–5 technique aims at mapping data from high-dimensional space
Although multidimensional projections provide a way to explore a data set by means of similarity, an inherent problem of the dimensionality reduction is the overlap among markers when a visual metaphor is created, that is, when the circles or images used to represent the instances of a data set are obstructed by other markers. This problem impairs the user exploration and detailed analysis of instances similarity, requiring several interactions with the graphical representation. From the literature, several techniques were proposed to remove markers overlap from bi-dimensional spaces8–10 or to place markers directly without overlap11–13—each technique works in a distinct way in order to generate a graphical representation without overlap. Some works presented an initial evaluation of these techniques,8,14 but they did not detail the quality of dimensionality reduction preservation when comparing the original with the modified layout, as well as if the quality of clusters are affected by layout arrangements from the initial projection, that is, the projection in which no overlap removal algorithm has yet been applied. Gansner and Hu 9 performed a comparison, but only according to few techniques focusing on graph structures, in which results cannot be transferred to multidimensional projections, and Nieto et al. 8 did not provide results that consider the scalability of the methods discussed here. Moreover, these works did not compare the overlap removal techniques with point placement technique that provide a layout arrangement without markers overlap, such as IncBoard 11 and HexBoard. 12 Therefore, in this work, we performed a detailed analysis of the quality of these techniques to choose the most accurate and precise technique, that is, the technique capable of maintaining the relationship between original multidimensional space, the original layout from a projection technique, and cluster formation. In addition, we also perform comparisons of time spent on execution. Thus, the main contributions of this article are the detailed analysis of overlap removal techniques and a guide for experts that intend to use or propose overlap removal techniques.
The remainder of this article is organized as follows. In section “Related works,” we provide a brief discussion on related works. In section “Overlap removal techniques,” we provide a description about the analyzed techniques. In section “Evaluation of overlap removal techniques,” we show the evaluation of the discussed techniques using different data sets. In section “Discussion,” we provide a discussion on the results. Finally, in section “Conclusion,” we conclude our work.
Related works
Although it does not present a new evaluation technique, our work is related to the works well discussed in the literature in which the objective is to present techniques that can be used to evaluate multidimensional projection according to class separation. However, in this work, we present such evaluation in the context of overlap removal techniques.
Tatu et al. 15 presented a study where they investigate the relationship between clusters interpretation in two-dimensional (2D) scatter plots and measures automatically extracted out of them. The measures used were Class Consistency Measure (CCM), 16 which is based on the distance of points to the centroid; Histogram Density Measure (one-dimensional (1D) and 2D), 17 which is a quality measure for scatter plots and it considers the class distribution of points in 2D scatter plots when they are projected on the axes; and Class Density Measure (CDM), 17 which evaluates scatter plots according to separation properties. Other work, presented by Seldmair et al., 18 performs a systematic evaluation of separation measures with respect to human judgment. In addition, they also provide a taxonomy of factors of data sets that influence the visual separation of clusters. Sedlmair et al. 18 used the Distance Consistency, 16 Distribution Consistency, and 2D-Histogram Density 17 measures. Recently, Wang et al. 19 proposed a technique where the dimensionality reduction process is performed according to visual separation measures that simulate the human perception of class separation. The approach, perception-driven dimensionality reduction, aims at generating truthful 2D linear projections with maximal class separation.
Unlike the previous works, in this article, we present an evaluation of the behavior of overlap removal techniques using data sets with a sort of characteristics, that is, we evaluate multidimensional projections according to separation and similarity measures before and after the overlap removal process.
Overlap removal techniques
Overlap removal techniques use different approaches to remove overlap among markers—while some techniques create well-defined structures in the plane, other techniques only consider a graph structure among the instances. By minimizing overlap among instances, these techniques have to maintain the relations imposed by multidimensional projection as much as possible.
In this section, two techniques that perform dimensionality reduction and point placement without markers overlap are presented, as well as four techniques intended to remove markers overlap from bi-dimensional spaces, that is, the techniques that attempt to remove the overlap of circles or scatter plot images in the plane.
Incremental board
The Incremental Board/IncBoard 11 technique is a projection technique that projects instance without overlap by creating a mosaic of rectangles so that each instance is placed in an individual cell of the mosaic. The IncBoard technique considers the relative similarity of elements in an n-dimensional space to project them in the 2D space and it does not allow data overlap because it does not allow more than one element in the same cell. Using the chess board structure, IncBoard allows an effective progressive visualization of dynamic data sets, resulting in a space that maintains a coherent layout. 11
HexBoard
Despite of the advantages offered by IncBoard technique, it does not convey the similarity of pairs between neighboring elements effectively. Since the IncBoard technique places elements inside a grid of fixed size cells, the neighborhood of a given element is always placed at same distance, even if not all of elements in its neighborhood are equally similar. The HexBoard 12 technique improves the ability to convey similarity information of each pair of neighbors in the projection, besides supporting identification of clusters and subsets of elements. This improvement is given by the fact that the HexBoard technique uses a mosaic of hexagons instead of the chess board structure.
ProjSnippet
The ProjSnippet
8
is a technique to visualize Web search results. The technique has an overlap removal approach which is described by the sum of two components, one that quantifies the overlapping snippets
where the parameter
The
RWordle
The Wordle technique is used to define positions of tags in a tag cloud visualization. Each tag is placed in a random point and when there is overlap, the new tag is moved in spiral fashion to its initial position, in order to find a valid position (with no overlap).
In the RWordle 13 technique, the main idea is to provide an orthogonal ordering for the Wordle technique. In this case, two geometry-directed ordering can be used:
Linear Ordering (RWordle-L), where items are classified along a scan line operating over an axis defined by an angle
Concentric Ordering (RWordle-C), where items are classified along their distances to the geometric center of the scene.
Variable placement with separation constraints
The variable placement with separation constraint (VPSC) 10 technique attempts to remove overlap among instances using two steps: first, a set of non-overlap constraints is generated; second, the objective is to find an optimal solution based on this set.
The non-overlap constraints are generated with a line-sweep algorithm that can detect overlap among rectangles. Based on the non-overlap constraints set, the VPSC technique solves the VPSC problem: Let there be n weighted variables
Note that the VPSC considers the problem for one axis at a time, and
PRISM
The PRISM 9 technique attempts to remove nodes overlap by defining a proximity graph that maintains proximity relations. At first, an overlap factor is computed for each edge of the proximity graph as follows
where
Graph nodes that do not overlap take the value
or, generating a layout in which each edge of the proximity graph has the length close to
Evaluation of overlap removal techniques
This section presents a detailed evaluation of the techniques presented in the section “Overlap Removal Techniques.” The metrics used in the evaluation are based on the similarity between the original and the new layout without overlap, the similarity preservation of multidimensional space and the bi-dimensional space without overlap, and clustering of instances belonging to a same class. These metrics are described in the following.
Evaluation metrics
Layout similarity
Gansner and Hu
9
proposed a technique to assess the similarity preservation between two layouts. The Layout Similarity metric attempts to quantify how much of the neighborhood structures are affected by the overlap removal algorithm. The technique calculates the length of the edges of the proximity graph created using a Delaunay Triangulation. Given the initial layout
where
Stress
The Stress 2 metric attempts to measure how much information was lost during the dimensionality reduction process. Thus, the lower the Stress value, the better the preservation of the structures.
An approach to compute the Stress (S) is shown in equation (6)
where d is the dissimilarities in the n-dimensional space and
Silhouette coefficient
The Silhouette Coefficient 20 metric is used to interpret clusters consistency and to provide information on how well each instance participates in its cluster. In this case, the clusters were generated by the projection technique, which defines the clusters according to the instances’ classes.
Given an instance i, its silhouette can be computed as follows:
Compute the mean of the distances
Compute the smallest distance
The silhouette of i is
The value of Silhouette Coefficient is given by the silhouette mean of all instances varying from −1 to 1, where better cohesion and cluster separation are indicated by values closer to 1.
Neighborhood Preservation
The Neighborhood Preservation (NP)
21
metric is used to assess the projection techniques according to the NP. The k-nearest neighbors of each instance in both multidimensional and projected space are used to compute the preservation ratio, that is, how many neighbors are preserved in the visual space. For example, suppose
The NP value for a projection is the mean of the NP applied to each instance in the data set. Notice that, as the value of k grows, the NP technique shows better results since more instances are maintained as nearest neighbors.
Neighborhood Hit
The Neighborhood Hit (NH)
5
metric is used to assess the class perception in the projection. In this way, a projection technique shows good results if it is able to consistently separate the classes in the data set. For each instance, its k-nearest neighbors are found to compute the ratio of neighbors belonging to the same class of the evaluated instance. For example, suppose
The NH value for a projection is the mean of the NH applied for each instance in the whole projection. Unlike the NP technique, as the value of k grows, the lower is the value of NH technique since more instances will have neighbors with different classes.
Results
The experiments were performed upon eight data sets, described in the following:
Diabetes—contains 169 instances representing patients described by eight features. The data set has 27 kB;
Iris—contains 150 instances representing types of iris plants divided in three classes of 50 instances. The data set has 4 kB;
Mammals—contains 183 instances automatically generated to represent different mammal features. Each mammal is described by 72 features. The data set has 4189 kB;
medical12—contains 72 instances of magnetic resonance images (MRI) described by 28 features. It is a reduced data set from original one with 540 instances. The data set has 21 kB;
medical6—contains 270 instances of MRI images described by 28 features. All instances are classified according to 12 classes. The data set has 75 kB;
CBR-ILP-IR—contains 574 instances of documents described by 43 terms. Each instance represents a paper of three fields: Case-Based-Reasoning (CBR), Intuitive Logic Programming (ILP), and Information Retrieval (IR). The data set has 154 kB;
Corel-1000—contains 1000 images from natural scenes and artificial objects separated in 10 classes. The data set has 1183 kB;
Segmentation—contains 2100 instances of features of 3 × 3 regions of a set of seven outdoor manually segmented images described by 19 features. The data set has 404 kB.
Notice that all data sets are structured and the greater the number of instances, the greater the probability of overlapping markers when projecting on the 2D space. Thus, the number of instances is what matters for evaluating the overlap removal techniques. After the dimensionality reduction, all data sets will have only two coordinates and the number of instances will dictate their size. The tests were performed on a computer with the following configuration: Intel Core i5-4440 CPU, 8 GB RAM, Windows 7 64 bits. The initial projections of the data sets and algorithms’ results can be visualized in Figure 1. By initial projection, we mean the data set projected on 2D space and no overlap removal technique has been applied to it.

Layouts generated by application of overlap removal techniques.
While the Corel-1000 and segmentation data sets were projected using the least square projection (LSP)
5
technique, the remaining data sets were projected using the Interactive Document Map (IDMAP)
4
technique. We used LSP only for those two data sets due to the complexity of the projection technique. The IDMAP is a force-based technique that tries to make the relations among instances in the multidimensional space and in the projected space to be constant for each pair of points (
In this work, we used a scan line for the RWordle-L technique with inclination angle of
The following sections represent the results according to the selected metrics and time evaluation.
Layout similarity
The Layout Similarity metric was applied among the initial projection and the layout generated by overlap removal techniques. Figure 2 shows the results for each technique and data set.

Heatmap conveying the Layout Similarity values for each technique and data set. The closer the color is to blue, the better the layout preservation.
It is possible to realize that the ProjSnippet technique was able to reach the best result for all of the data sets. This is explained by the fact that besides being able to maintain the relations imposed by the projected technique, the ProjSnippet technique could not remove all overlap in the initial layout. Such issue helps the ProjSnippet technique to generate layouts more similar to the initial projections. In the second place is the PRISM technique. Since it was able to remove all overlap in the initial layout, its result according to Layout Similarity was higher than the ProjSnippet, that is, the generated layouts were less similar to the original layouts if compared with ProjSnippet layouts. The VPSC technique removes overlaps according to the axes x and y, which does not provide much flexibility in the overlap removal process and leads to a stretched layout. Finally, since the RWordle-C and RWordle-L do not consider any structural relation, both techniques presented higher values of Layout Similarity. Notice that, for small data sets the RWordle-C and RWordle-L techniques performed well.
Silhouette coefficient
As already defined, the Silhouette Coefficient is a measure to evaluate the separation and cohesion of instances belonging to a same class—the better layouts are that whose values are close to

Heatmap conveying the Silhouette Coefficient values for each technique and data set. The closer the color to red, the better the class understanding.
Unlike what happened with the results of Layout Similarity, the ProjSnippet technique presented similar results to the PRISM technique. Notice that these techniques are benefited from the fact that the Silhouette Coefficient metric presents better results when blank space is added between clusters of different classes. The IncBoard and HexBoard technique, in general, presented lower result mainly because they place instances of different classes in a close neighborhood. However, notice that for the medical6 data set, the IncBoard and HexBoard techniques were able to project instances separated by classes, leading to better results.
Stress
The Stress is commonly used to measure how much the instances similarity was preserved from the original multidimensional space in the projected space. Figure 4 shows the Stress values for all new layouts computed with the overlap removal techniques.

Heatmap conveying the Stress values for each technique and data set. The closer the color is to blue, the better the similarity preservation among instances.
In general, all overlap removal techniques got satisfactory results for the Stress technique. On the other hand, since the IncBoard and HexBoard technique compact instances which were distant to each other in the multidimensional space, the Stress values for these techniques are higher.
NP
The NP metric is another way to identify how much the instance relationships were preserved in the dimensionality reduction. The metric was computed for all projections considering up to 30 neighbors. These results are presented in Figures 5 and 6. Note that, for NP and NH graphics, the initial projection is labeled as the name of the data set and the line that represents its analysis is colored by black.

Performance of each technique on Neighborhood Preservation metric for the data sets (a) diabetes, (b) iris, (c) mammals, and (d) medical12. The neighborhood preservation rate was calculated varying the number of neighbors from 1 to 30. Notice that the black line is the original projection.

Performance of each technique on Neighborhood Preservation metric for the data sets (a) medical6, (b) CBR-ILP-IR, (c) Corel-1000, and (d) segmentation. The neighborhood preservation rate was calculated varying the number of neighbors from 1 to 30. Notice that the black line is the original projection.
In general, the techniques presented similar results for all data sets. As we commented on the evaluation of Layout Similarity, since the ProjSnippet technique was not capable of removing all overlap among instances for the majority of data sets, its result according to NP follows the initial projection result—logically, the ability of preserving structural and similarity relations contributed to the good result as well. The RWordle-C and RWordle-L technique reached the worst results for the data sets in which classes are too closer to each other, that is, the medical6 and segmentation data sets. For small data sets, we can see that the IncBoard and HexBoard technique affect much more the structures imposed by a multidimensional projection.
NH
NH can be used to measure the quality of groups of instances belonging to a same class. Thus, we can measure if the removal overlap technique messed the groups computed in the original projection. The NH results are presented in Figures 7 and 8.

Performance of each technique on Neighborhood Hit metric for the data sets (a) diabetes, (b) iris, (c) mammals, and (d) medical12. The metric values were calculated varying the number of neighbors from 1 to 30. Notice that the black line is the original projection.

Performance of each technique on Neighborhood Hit metric for the data sets (a) medical6, (b) CBR-ILP-IR, (c) Corel-1000, and (d) segmentation. The metric values were calculated varying the number of neighbors from 1 to 30. Notice that the black line is the original projection.
As happened in the NP metric, the techniques presented similar results for all data sets according to the NH metric. Again, since the RWordle-C and RWordle-L technique do not consider any structural relation among instances, such techniques presented difficult and even bigger results for the medical6 and segmentation data sets. Besides that, by compacting layouts, the IncBoard and HexBoard impair their results according to the NH metric since instance of different classes are placed closer to each other. This effect can be verified for the mammals data set (see Figure 7(c)).
It is worth noticing that the ProjSnippet technique follows the results of initial projection. Figure 9 shows the distribution for each technique according to the NH and NP metrics. The PRISM and ProjSnippet techniques presented similar aspects; however, the ProjSnippet technique gave more similar results to the initial projection. Notice that according to NH metric, the distribution given by IncBoard and HexBoard techniques are smaller since the metric is only affected by boundaries of classes, which are treated by IncBoard and HexBoard techniques.

Distributions of the values of (a) Neighborhood Preservation and (b) Neighborhood Hit metrics, varying the number of neighbors from 1 to 30. Each colored point represents a Neighborhood Hit/Neighborhood Preservation rate value referring to the number of neighbors.
Running time
The running time of each technique can be verified in Figure 10. The time spent by each technique is indicated in milliseconds using logarithmic scale

Heatmap conveying the time execution in logarithmic scale for each technique and data set. The close the color is to yellow, the better the technique according to time execution.
We can see that only VPSC and ProjSnippet were able to maintain an acceptable time execution for all data sets. The number of instances influenced the time required to complete the process for each algorithm; however, notice that the RWordle-L and RWordle-C techniques were also impaired by the amount of overlap in the initial projection. Since the number of edges in the proximity graph is proportional to the overlap rate, the PRISM technique was also affected by this issue.
Experiment with users
We performed an experiment with 26 users to classify the overlap removal techniques according to the similarity of the layout generated by them and the original layout. For each data set, we presented the initial layout and the rearranged layout for each overlap removal technique. The users chose a score between 0 and 10 to evaluate the generated projection according to its similarity.
By analyzing the boxplots of Figure 11, we can see that the scores given to the techniques are similar, but with slightly better results to the PRISM and RWordle-C techniques. These results give us the insight that the user sometimes cannot rely only on the visualization aspects to define how good or not a projection is after overlap removal, that is, they cannot realize how much the generated layout approaches the initial layout.

Boxplots of the scores given by the users to the overlap removal techniques.
Discussion
Based on experiments, it is possible to emphasize some characteristics for each technique evaluated in this article:
The IncBoard and HexBoard techniques, in general, got lower performance if compared with the remaining techniques. However, these techniques can be a good option when instances of different classes share boundaries since these techniques can be used to emphasize the separation between such instances. Besides that, the layout generated is more compact;
The RWordle-C and RWordle-L techniques are less sophisticated, therefore, their performances in critical layout characteristics, such as NP, are impaired if they are applied in layouts with high overlap rate or where groups of instances are near each other. However, they are efficient in the cases where the instances are relatively separate;
The ProjSnippet technique is recommended for layouts with low overlap rate and in data sets with a limited number of instances. The technique has the characteristic of preserving the structures, clusters, and the initial neighborhood effectively, besides being possible to control the contribution of the overlap removal and NP. It is interesting to note that the assessment results for the ProjSnippet technique follow the results given by the initial projection;
The VPSC technique presents layouts that preserve the general structure of the projection, besides having a good performance considering the processing time. It can remove the overlap with effectivity, but it presents medium performance if compared with the other techniques;
The PRISM technique presents good results according to the metrics, being able to preserve the relations imposed by the initial projection and removing any overlap. However, this technique has the difficulty of maintaining the layout compact, and the time execution can be impaired by the number of instances and overlap rate increase.
Considering these characteristics, Figure 12 shows a simple guide to use an algorithm based on the requirement of structures preservation and size of a data set—notice that a data set with moderate size and high overlap rate can be evaluated as a big data set for some techniques, such as PRISM, RWordle-L, and RWordle-C. Moreover, it is not recommended to apply IncBoard and HexBoard techniques on data sets with small or moderate sizes since the trade-off between compact shape and structures preservation is not worth it. Although the classification of data set sizes according to small, medium, and big can be empirical, we can define small data sets as those with a maximum of 200 instances, medium data sets as those with size between 201 and 500, and big data sets as those with size above 501. However, note that the overlap rate can influence the choice of a technique, that is, the higher the overlap rate, the bigger a data set must be considered.

Visual guide to apply an overlap removal algorithm. The structures preservation means how much structures are needed to be preserved, that is, if one wants to preserve as many structures as possible, the option “High” must be chosen. If the final projection can have a balanced structures preservation, the option “Medium” must be chosen, and one can choose “Low” if no or low structures must be preserved. On the data set size, the choice must be made based on the complexity of the algorithms. If the data set size can be executed through
Conclusion
Multidimensional projections are a powerful tool for analyzing instances similarity of a multidimensional data set. Such techniques provide an overview about the data and allow identification of trends, patterns, and anomalies. However, the overlap problem faced by the visual representation impairs effective analysis since it is not possible to visualize the data instances that are obstructed. So that, overlap removal techniques can be applied to decrease the overlap problem and to improve the exploratory process.
In this work, we compared four approaches to remove overlap in the visual metaphor used to show multidimensional projections, as well as two projection techniques that already consider overlap among markers. These techniques were evaluated according to class, neighborhood, and structure preservation. Through the assessment of these techniques, it was possible to highlight their characteristics and to provide a simple guide to choose an overlap removal or a multidimensional projection technique that has non-overlap constraints.
As future works, it is necessary to investigate these algorithms by means of clusters, in other words, the overlap removal process can be applied on small portions of data. As it was possible to realize based on results, the algorithms tend to present better layouts when applied on small data sets. If the overlap removal process is applied selectively, the issues faced by some algorithms can be decreased, such as display screen increment and processing delay.
