Maximally stable extremal regions

In computer vision, maximally stable extremal regions (MSER) are used as a method of blob detection in images. This technique was proposed by Matas et al. to find correspondences between image elements from two images with different viewpoints. This method of extracting a comprehensive number of corresponding image elements contributes to the wide-baseline matching, and it has led to better stereo matching and object recognition algorithms.

Terms and definitions

Image [math]\displaystyle{ I }[/math] is a mapping [math]\displaystyle{ I : D \subset \mathbb{Z}^2 \to S }[/math]. Extremal regions are well defined on images if:

1. [math]\displaystyle{ S }[/math] is totally ordered (total, antisymmetric and transitive binary relations [math]\displaystyle{ \le }[/math] exist).
2. An adjacency relation [math]\displaystyle{ A \subset D \times D }[/math] is defined. We will denote that two points are adjacent as [math]\displaystyle{ pAq }[/math].

Region [math]\displaystyle{ Q }[/math] is a contiguous (aka connected) subset of [math]\displaystyle{ D }[/math]. (For each [math]\displaystyle{ p,q \in Q }[/math] there is a sequence [math]\displaystyle{ p, a_1, a_2, .., a_n, q }[/math] such as [math]\displaystyle{ pAa_1, a_1Aa_2, \dots, a_{n-1}Aa_{n}, a_nAq }[/math].) Note that under this definition the region can contain "holes" (for example, a ring-shaped region is connected, but its internal circle is not the part of [math]\displaystyle{ Q }[/math]).

(Outer) region boundary [math]\displaystyle{ \partial Q = \{ q \in D \setminus Q: \exists p \in Q : qAp \} }[/math], which means the boundary [math]\displaystyle{ \partial Q }[/math] of [math]\displaystyle{ Q }[/math] is the set of pixels adjacent to at least one pixel of [math]\displaystyle{ Q }[/math] but not belonging to [math]\displaystyle{ Q }[/math]. Again, in case of regions with "holes", the region boundary is not obliged to be connected subset of [math]\displaystyle{ D }[/math] (a ring has inner bound and outer bound which do not intersect).

Extremal region [math]\displaystyle{ Q \subset D }[/math] is a region such that either for all [math]\displaystyle{ p \in Q, q \in \partial Q : I(p) \gt I(q) }[/math] (maximum intensity region) or for all [math]\displaystyle{ p \in Q, q \in \partial Q : I(p) \lt I(q) }[/math] (minimum intensity region). As far as [math]\displaystyle{ S }[/math] is totally ordered, we can reformulate these conditions as [math]\displaystyle{ \min(I(p)) \gt \max(I(q)) }[/math] for maximum intensity region and [math]\displaystyle{ \max(I(p)) \lt \min(I(q)) }[/math] for minimum intensity region, respectively. In this form we can use a notion of a threshold intensity value which separates the region and its boundary.

Maximally stable extremal region Let [math]\displaystyle{ Q_i }[/math] an extremal region such as all points on it have an intensity smaller than [math]\displaystyle{ i \in S }[/math]. Note [math]\displaystyle{ Q_i \subset Q_{i+\Delta} }[/math] for all positive [math]\displaystyle{ \Delta \in S }[/math]. Extremal region [math]\displaystyle{ Q_{i*} }[/math] is maximally stable if and only if [math]\displaystyle{ | Q_{i+\Delta} \setminus Q_{i-\Delta} | / |Q_i| }[/math] has a local minimum at [math]\displaystyle{ i* }[/math]. (Here [math]\displaystyle{ | \cdot | }[/math] denotes cardinality). [math]\displaystyle{ \Delta \in S }[/math] is here a parameter of the method.

The equation checks for regions that remain stable over a certain number of thresholds. If a region [math]\displaystyle{ Q_{i+\Delta} }[/math] is not significantly larger than a region [math]\displaystyle{ Q_{i-\Delta} }[/math], region [math]\displaystyle{ Q_i }[/math] is taken as a maximally stable region.

The concept more simply can be explained by thresholding. All the pixels below a given threshold are 'black' and all those above or equal are 'white'. Given a source image, if a sequence of thresholded result images [math]\displaystyle{ I_t }[/math] is generated where each image [math]\displaystyle{ t }[/math] corresponds to an increasing threshold t, first a white image would be seen, then 'black' spots corresponding to local intensity minima will appear then grow larger. A maximally stable extremal region is found when size of one of these black areas is the same (or near the same) than in previous image.

These 'black' spots will eventually merge, until the whole image is black. The set of all connected components in the sequence is the set of all extremal regions. In that sense, the concept of MSER is linked to the one of component tree of the image. The component tree indeed provide an easy way for implementing MSER.

Extremal regions

Extremal regions in this context have two important properties, that the set is closed under...

1. continuous transformation of image coordinates. This means it is affine invariant and it doesn't matter if the image is warped or skewed.
2. monotonic transformation of image intensities. The approach is of course sensitive to natural lighting effects as change of day light or moving shadows.

Advantages of MSER

Because the regions are defined exclusively by the intensity function in the region and the outer border, this leads to many key characteristics of the regions which make them useful. Over a large range of thresholds, the local binarization is stable in certain regions, and have the properties listed below.

• Invariance to affine transformation of image intensities
• Covariance to adjacency preserving (continuous) transformation [math]\displaystyle{ T : D \to D }[/math] on the image domain
• Stability: only regions whose support is nearly the same over a range of thresholds is selected.
• Multi-scale detection without any smoothing involved, both fine and large structure is detected.
  Note, however, that detection of MSERs in a scale pyramid improves repeatability, and number of correspondences across scale changes.
• The set of all extremal regions can be enumerated in worst-case [math]\displaystyle{ O(n) }[/math], where [math]\displaystyle{ n }[/math] is the number of pixels in the image.

Comparison to other region detectors

In Mikolajczyk et al., six region detectors are studied (Harris-affine, Hessian-affine, MSER, edge-based regions, intensity extrema, and salient regions). A summary of MSER performance in comparison to the other five follows.

• Region density – in comparison to the others MSER offers the most variety detecting about 2600 regions for a textured blur scene and 230 for a light changed scene, and variety is generally considered to be good. Also MSER had a repeatability of 92% for this test.
• Region size – MSER tended to detect many small regions, versus large regions which are more likely to be occluded or to not cover a planar part of the scene. Though large regions may be slightly easier to match.
• Viewpoint change – MSER outperforms the five other region detectors in both the original images and those with repeated texture motifs.
• Scale change – Following Hessian-affine detector, MSER comes in second under a scale change and in-plane rotation.
• Blur – MSER proved to be the most sensitive to this type of change in image, the only area in which this type of detection is lacking.
  Note however that this evaluation did not make use of multi-resolution detection, which has been shown to improve repeatability under blur.
• Light change – MSER showed the highest repeatability score for this type of scene, with all the other having good robustness as well.

MSER consistently resulted in the highest score through many tests, proving it to be a reliable region detector.

Implementation

The original algorithm of Matas et al. is [math]\displaystyle{ O(n\,\log(\log(n))) }[/math] in the number [math]\displaystyle{ n\, }[/math] of pixels. It proceeds by first sorting the pixels by intensity. This would take [math]\displaystyle{ O(n)\, }[/math] time, using BINSORT. After sorting, pixels are marked in the image, and the list of growing and merging connected components and their areas is maintained using the union-find algorithm. This would take [math]\displaystyle{ O(n\,\log(\log(n))) }[/math] time. In practice these steps are very fast. During this process, the area of each connected component as a function of intensity is stored producing a data structure. A merge of two components is viewed as termination of existence of the smaller component and an insertion of all pixels of the smaller component into the larger one. In the extremal regions, the 'maximally stable' ones are those corresponding to thresholds where the relative area change as a function of relative change of threshold is at a local minimum, i.e. the MSER are the parts of the image where local binarization is stable over a large range of thresholds.

The component tree is the set of all connected components of the thresholds of the image, ordered by inclusion. Efficient (quasi-linear whatever the range of the weights) algorithms for computing it do exist. Thus this structure offers an easy way for implementing MSER.

More recently, Nister and Stewenius have proposed a truly (if the weight are small integers) worst-case [math]\displaystyle{ O(n)\, }[/math] method in, which is also much faster in practice. This algorithm is similar to the one of Ph. Salembier et al.

Robust wide-baseline algorithm

The purpose of this algorithm is to match MSERs to establish correspondence points between images. First MSER regions are computed on the intensity image (MSER+) and on the inverted image (MSER-). Measurement regions are selected at multiple scales: the size of the actual region, 1.5x, 2x, and 3x scaled convex hull of the region. Matching is accomplished in a robust manner, so it is better to increase the distinctiveness of large regions without being severely affected by clutter or non-planarity of the region's pre-image. A measurement taken from an almost planar patch of the scene with stable invariant description are called a 'good measurement'. Unstable ones or those on non-planar surfaces or discontinuities are called 'corrupted measurements'. The robust similarity is computed: For each [math]\displaystyle{ M_A^i }[/math] on region [math]\displaystyle{ A, k }[/math] regions [math]\displaystyle{ B_1,\dots, B_k }[/math] from the other image with the corresponding i-th measurement [math]\displaystyle{ M_{B_1}^i ,\dots, M_{B_k}^i }[/math] nearest to [math]\displaystyle{ M_A^i }[/math] are found and a vote is cast suggesting correspondence of A and each of [math]\displaystyle{ B_1, \dots , B_k }[/math]. Votes are summed over all measurements, and using probability analysis, 'good measurements' can be picked out as the 'corrupt measurements' will likely spread their votes randomly. By applying RANSAC to the centers of gravity of the regions, a rough epipolar geometry can be computed. An affine transformation between pairs of potentially corresponding regions is computed, and correspondences define it up to a rotation, which is then determined by epipolar lines. The regions are then filtered, and the ones with correlation of their transformed images above a threshold are chosen. RANSAC is applied again with a more narrow threshold, and the final epipolar geometry is estimated by the eight-point algorithm.

This algorithm can be tested here (Epipolar or homography geometry constrained matches): WBS Image Matcher

Use in text detection

The MSER algorithm has been used in text detection by Chen by combining MSER with Canny edges. Canny edges are used to help cope with the weakness of MSER to blur. MSER is first applied to the image in question to determine the character regions. To enhance the MSER regions any pixels outside the boundaries formed by Canny edges are removed. The separation of the later provided by the edges greatly increase the usability of MSER in the extraction of blurred text.

An alternative use of MSER in text detection is the work by Shi using a graph model. This method again applies MSER to the image to generate preliminary regions. These are then used to construct a graph model based on the position distance and color distance between each MSER, which is treated as a node. Next the nodes are separated into foreground and background using cost functions. One cost function is to relate the distance from the node to the foreground and background. The other penalizes nodes for being significantly different from its neighbor. When these are minimized the graph is then cut to separate the text nodes from the non-text nodes.

To enable text detection in a general scene, Neumann uses the MSER algorithm in a variety of projections. In addition to the greyscale intensity projection, he uses the red, blue, and green color channels to detect text regions that are color distinct but not necessarily distinct in greyscale intensity. This method allows for detection of more text than solely using the MSER+ and MSER- functions discussed above.

Extensions and adaptations

• The MSER algorithm has been adapted to colour images, by replacing thresholding of the intensity function with agglomerative clustering, based on colour gradients.
• The MSER algorithm can be used to detect regions based on color as opposed to intensity. This is done by Chavez by creating an intensity function for red, green, and blue in the HSV color space. The MSER algorithm is then run five times; over the three color pseudo-intensities and then over the grey scale intensities using the standard MSER+ and MSER- functions.
• The MSER algorithm can be used to track colour objects, by performing MSER detection on the Mahalanobis distance to a colour distribution.
• By detecting MSERs in multiple resolutions, robustness to blur, and scale change can be improved.

Other applications

• Shape Descriptors for Maximally Stable Extremal Regions
• Efficient Maximally Stable Extremal Region (MSER) Tracking
• N-tree Disjoint-Set Forests for Maximally Stable Extremal Regions
• Video Google and Object Level Grouping for Video Shots
• Real-Time Extraction of Maximally Stable Extremal Regions on an FPGA
• Maximally Stable Colour Regions for Recognition and Matching

See also

• Blob detection
• Feature detection (computer vision)

External links

• VLFeat, an open source computer vision library in C (with a MEX interface to MATLAB), including an implementation of MSER
• OpenCV, an open source computer vision library in C/C++, including an implementation of Linear Time MSER
• Detector Repeatabilty Study, Kristian Mikolajczyk Binaries (Win/Linux to compute MSER/HarrisAffine... . Binary used in his repeatability study.
• Linear Time MSER Implementation, Charles Dubout, C++ implementation of MSER as a blob detector