M-tree

Short description: Tree data structure

In computer science, M-trees are tree data structures that are similar to R-trees and B-trees. It is constructed using a metric and relies on the triangle inequality for efficient range and k-nearest neighbor (k-NN) queries. While M-trees can perform well in many conditions, the tree can also have large overlap and there is no clear strategy on how to best avoid overlap. In addition, it can only be used for distance functions that satisfy the triangle inequality, while many advanced dissimilarity functions used in information retrieval do not satisfy this.

Overview

2D M-Tree visualized using ELKI. Every blue sphere (leaf) is contained in a red sphere (directory nodes). Leaves overlap, but not too much; directory nodes overlap much more here.

As in any tree-based data structure, the M-tree is composed of nodes and leaves. In each node there is a data object that identifies it uniquely and a pointer to a sub-tree where its children reside. Every leaf has several data objects. For each node there is a radius [math]\displaystyle{ r }[/math] that defines a Ball in the desired metric space. Thus, every node [math]\displaystyle{ n }[/math] and leaf [math]\displaystyle{ l }[/math] residing in a particular node [math]\displaystyle{ N }[/math] is at most distance [math]\displaystyle{ r }[/math] from [math]\displaystyle{ N }[/math], and every node [math]\displaystyle{ n }[/math] and leaf [math]\displaystyle{ l }[/math] with node parent [math]\displaystyle{ N }[/math] keep the distance from it.

M-tree construction

Components

An M-tree has these components and sub-components:

1. Non-leaf nodes
   1. A set of routing objects N_RO.
   2. Pointer to Node's parent object O_p.
2. Leaf nodes
   1. A set of objects N_O.
   2. Pointer to Node's parent object O_p.
3. Routing Object
   1. (Feature value of) routing object O_r.
   2. Covering radius r(O_r).
   3. Pointer to covering tree T(O_r).
   4. Distance of O_r from its parent object d(O_r,P(O_r))
4. Object
   1. (Feature value of the) object O_j.
   2. Object identifier oid(O_j).
   3. Distance of O_j from its parent object d(O_j,P(O_j))

Insert

The main idea is first to find a leaf node N where the new object O belongs. If N is not full then just attach it to N. If N is full then invoke a method to split N. The algorithm is as follows:

Algorithm Insert
  Input: Node N of M-Tree MT, Entry [math]\displaystyle{ O_{n} }[/math]
  Output: A new instance of MT containing all entries in original MT plus [math]\displaystyle{ O_{n} }[/math]

  [math]\displaystyle{ N_{e} \gets N }[/math]'s routing objects or objects
  if N is not a leaf then
  {
       /* Look for entries that the new object fits into */
       let [math]\displaystyle{ N_{in} }[/math] be routing objects from [math]\displaystyle{ N_{e} }[/math]'s set of routing objects [math]\displaystyle{ N_{RO} }[/math] such that [math]\displaystyle{ d(O_{r}, O_{n}) \le r(O_{r}) }[/math]
       if [math]\displaystyle{ N_{in} }[/math] is not empty then
       {
          /* If there are one or more entry, then look for an entry such that is closer to the new object */
          [math]\displaystyle{ O_{r}^{*} \gets \min_{O_{r}\in N_{in}} d(O_{r}, O_{n}) }[/math]
       }
       else
       {
          /* If there are no such entry, then look for an object with minimal distance from */ 
          /* its covering radius's edge to the new object */
          [math]\displaystyle{ O_{r}^{*} \gets \min_{O_{r}\in N_{in}} d(O_{r}, O_{n}) - r(O_{r}) }[/math]
          /* Upgrade the new radii of the entry */
          [math]\displaystyle{ r(O_{r}^{*}) \gets d(O_{r}^{*}, O_{n}) }[/math]
       }
       /* Continue inserting in the next level */
       return insert([math]\displaystyle{ T(O_{r}^{*}), O_{n} }[/math]);
  else
  {
       /* If the node has capacity then just insert the new object */
       if N is not full then
       { store([math]\displaystyle{ N, O_{n} }[/math]) }
       /* The node is at full capacity, then it is needed to do a new split in this level */
       else
       { split([math]\displaystyle{ N, O_{n} }[/math]) }
  }

• "←" denotes assignment. For instance, "largest ← item" means that the value of largest changes to the value of item.
• "return" terminates the algorithm and outputs the following value.

Split

If the split method arrives to the root of the tree, then it choose two routing objects from N, and creates two new nodes containing all the objects in original N, and store them into the new root. If split methods arrives to a node N that is not the root of the tree, the method choose two new routing objects from N, re-arrange every routing object in N in two new nodes [math]\displaystyle{ N_{1} }[/math] and [math]\displaystyle{ N_{2} }[/math], and store these new nodes in the parent node [math]\displaystyle{ N_{p} }[/math] of original N. The split must be repeated if [math]\displaystyle{ N_{p} }[/math] has not enough capacity to store [math]\displaystyle{ N_{2} }[/math]. The algorithm is as follow:

Algorithm Split
  Input: Node N of M-Tree MT, Entry [math]\displaystyle{ O_{n} }[/math]
  Output: A new instance of MT containing a new partition.

  /* The new routing objects are now all those in the node plus the new routing object */
  let be NN entries of [math]\displaystyle{ N \cup O }[/math]
  if N is not the root then
  {
     /*Get the parent node and the parent routing object*/
     let [math]\displaystyle{ O_{p} }[/math] be the parent routing object of N
     let [math]\displaystyle{ N_{p} }[/math] be the parent node of N
  }
  /* This node will contain part of the objects of the node to be split */
  Create a new node N'
  /* Promote two routing objects from the node to be split, to be new routing objects */
  Create new objects [math]\displaystyle{ O_{p1} }[/math] and [math]\displaystyle{ O_{p2} }[/math].
  Promote([math]\displaystyle{ N, O_{p1}, O_{p2} }[/math])
  /* Choose which objects from the node being split will act as new routing objects */
  Partition([math]\displaystyle{ N, O_{p1}, O_{p2}, N_{1}, N_{2} }[/math])
  /* Store entries in each new routing object */
  Store [math]\displaystyle{ N_{1} }[/math]'s entries in N and [math]\displaystyle{ N_{2} }[/math]'s entries in N'
  if N is the current root then
  {
      /* Create a new node and set it as new root and store the new routing objects */
      Create a new root node [math]\displaystyle{ N_{p} }[/math]
      Store [math]\displaystyle{ O_{p1} }[/math] and [math]\displaystyle{ O_{p2} }[/math] in [math]\displaystyle{ N_{p} }[/math]
  }
  else
  {
      /* Now use the parent routing object to store one of the new objects */
      Replace entry [math]\displaystyle{ O_{p} }[/math] with entry [math]\displaystyle{ O_{p1} }[/math] in [math]\displaystyle{ N_{p} }[/math]
      if [math]\displaystyle{ N_{p} }[/math] is no full then
      {
          /* The second routing object is stored in the parent only if it has free capacity */
          Store [math]\displaystyle{ O_{p2} }[/math] in [math]\displaystyle{ N_{p} }[/math]
      }
      else
      {
           /*If there is no free capacity then split the level up*/
           split([math]\displaystyle{ N_{p}, O_{p2} }[/math])
      }
  }

• "←" denotes assignment. For instance, "largest ← item" means that the value of largest changes to the value of item.
• "return" terminates the algorithm and outputs the following value.

M-tree queries

Range query

A range query is where a minimum similarity/maximum distance value is specified. For a given query object [math]\displaystyle{ Q \in D }[/math] and a maximum search distance [math]\displaystyle{ r(Q) }[/math], the range query range(Q, r(Q)) selects all the indexed objects [math]\displaystyle{ O_j }[/math] such that [math]\displaystyle{ d(O_j, Q) \le r(Q) }[/math].

Algorithm RangeSearch starts from the root node and recursively traverses all the paths which cannot be excluded from leading to qualifying objects.

Algorithm RangeSearch
Input: Node N of M-Tree MT,  Q: query object, [math]\displaystyle{ r(Q) }[/math]: search radius

Output: all the DB objects such that [math]\displaystyle{ d(Oj, Q) \le r(Q) }[/math]

{ 
  let [math]\displaystyle{ O_{p} }[/math] be the parent object of node N;

  if N is not a leaf then { 
    for each entry([math]\displaystyle{ O_{r} }[/math]) in N do {
          if [math]\displaystyle{ | d(O_{p}, Q) - d(O_{r}, O_{p}) | \le r(Q) + r(O_{r}) }[/math] then { 
            Compute [math]\displaystyle{ d(O_{r}, Q) }[/math];
            if [math]\displaystyle{ d(O_{r}, Q) \le r(Q) +r(O_{r}) }[/math] then
              RangeSearch(*ptr([math]\displaystyle{ T(O_{r} }[/math])),Q,[math]\displaystyle{ r(Q) }[/math]); 
          }
    }
  }
  else { 
    for each entry([math]\displaystyle{ O_{j} }[/math]) in N do {
          if [math]\displaystyle{ | d(O_{p}, Q) - d(O_{j}, O_{p}) | \le r(Q) }[/math] then { 
            Compute [math]\displaystyle{ d(O_{j}, Q) }[/math];
            if [math]\displaystyle{ d(O_{j}, Q) }[/math] ≤ [math]\displaystyle{ r(Q) }[/math] then
              add [math]\displaystyle{ oid(O_{j}) }[/math] to the result;
          }
    }
  }
}

• "←" denotes assignment. For instance, "largest ← item" means that the value of largest changes to the value of item.
• "return" terminates the algorithm and outputs the following value.

• [math]\displaystyle{ oid(O_{j}) }[/math] is the identifier of the object which resides on a separate data file.
• [math]\displaystyle{ T(O_{r}) }[/math] is a sub-tree – the covering tree of [math]\displaystyle{ O_{r} }[/math]

k-NN queries

k-nearest neighbor (k-NN) query takes the cardinality of the input set as an input parameter. For a given query object Q ∈ D and an integer k ≥ 1, the k-NN query NN(Q, k) selects the k indexed objects which have the shortest distance from Q, according to the distance function d.

See also

• Segment tree
• Interval tree - A degenerate R-tree for one dimension (usually time)
• Bounding volume hierarchy
• Spatial index
• GiST
• Cover tree