Ranklet

Short description: Statistical method for pattern recognition.

In statistics, a ranklet is an orientation-selective non-parametric feature which is based on the computation of Mann–Whitney–Wilcoxon (MWW) rank-sum test statistics. Ranklets achieve similar response to Haar wavelets as they share the same pattern of orientation-selectivity, multi-scale nature and a suitable notion of completeness. There were invented by Fabrizio Smeralhi in 2002.

Rank-based (non-parametric) features have become popular in the field of image processing for their robustness in detecting outliers and invariance to monotonic transformations such as brightness, contrast changes and gamma correction.

The MWW is a combination of Wilcoxon rank-sum test and Mann–Whitney U-test. It is a non-parametric alternative to the t-test used to test the hypothesis for the comparison of two independent distributions. It assesses whether two samples of observations, usually referred as Treatment T and Control C, come from the same distribution but do not have to be normally distributed.

The Wilcoxon rank-sum statistics W_s is determined as:

[math]\displaystyle{ W_s=\sum^N_{i=1}\pi_i V_i \text{ where }\pi_i = \text{rank of element }i\text{ and }V_i = \begin{cases} 0 & \text{ for }\pi_i\in C \\[3pt] 1 & \text{ for }\pi_i\in T \end{cases} }[/math]

Subsequently, let MW be the Mann–Whitney statistics defined by:

[math]\displaystyle{ MW=W_s-\frac{m(m+1)}{2} }[/math]

where m is the number of Treatment values.

A ranklet R is defined as the normalization of MW in the range [−1, +1]:

[math]\displaystyle{ R = \frac{MW}{mn/2}-1 }[/math]

where a positive value means that the Treatment region is brighter than the Control region, and a negative value otherwise.

Example

Suppose [math]\displaystyle{ T=\lbrace 5,9,1,10,15\rbrace }[/math] and [math]\displaystyle{ C=\lbrace 20,4,7,13,19,11\rbrace }[/math] then

Intensity	1	4	5	7	9	10	11	13	15	19	20
Sample	T	C	T	C	T	T	C	C	T	C	C
Rank	1	2	3	4	5	6	7	8	9	10	11

• [math]\displaystyle{ W_s=\Big\lbrace 1+3+5+6+9\Big\rbrace =24 }[/math]
• [math]\displaystyle{ MW=24-[5 \times (5+1)/2]=9 }[/math]
• [math]\displaystyle{ R=[9/[5 \times 6/2]]-1=-0.4 }[/math]

Hence, in the above example the Control region was a little bit brighter than the Treatment region.

Method

Since Ranklets are non-linear filters, they can only be applied in the spatial domain. Filtering with Ranklets involves dividing an image window W into Treatment and Control regions as shown in the image below:

Subsequently, Wilcoxon rank-sum test statistics are computed in order to determine the intensity variations among conveniently chosen regions (according to the required orientation) of the samples in W. The intensity values of both regions are then replaced by the respective ranking scores. These ranking scores determine a pairwise comparison between the T and C regions. This means that a ranklet essentially counts the number of TxC pairs which are brighter in the T set. Hence a positive value means that the Treatment values are brighter than the Control values, and vice versa.

External links

• Matlab RankletFilter.m -> source file to decompose an image into Intensity Ranklets