opencv/modules/imgproc/doc/histograms.rst

248 lines
11 KiB
ReStructuredText
Raw Normal View History

Histograms
==========
.. highlight:: cpp
.. index:: calcHist
2011-03-09 06:22:24 +08:00
.. _calcHist:
calcHist
------------
.. c:function:: void calcHist( const Mat* arrays, int narrays, const int* channels, const Mat\& mask, MatND\& hist, int dims, const int* histSize, const float** ranges, bool uniform=true, bool accumulate=false )
.. c:function:: void calcHist( const Mat* arrays, int narrays, const int* channels, const Mat\& mask, SparseMat\& hist, int dims, const int* histSize, const float** ranges, bool uniform=true, bool accumulate=false )
2011-04-25 05:02:14 +08:00
Calculates a histogram of a set of arrays.
2011-04-25 05:02:14 +08:00
:param arrays: Source arrays. They all should have the same depth, ``CV_8U`` or ``CV_32F`` , and the same size. Each of them can have an arbitrary number of channels.
2011-04-25 05:02:14 +08:00
:param narrays: Number of source arrays.
2011-04-25 05:02:14 +08:00
:param channels: List of the ``dims`` channels used to compute the histogram. The first array channels are numerated from 0 to ``arrays[0].channels()-1`` , the second array channels are counted from ``arrays[0].channels()`` to ``arrays[0].channels() + arrays[1].channels()-1`` etc.
2011-04-25 05:02:14 +08:00
:param mask: Optional mask. If the matrix is not empty, it must be an 8-bit array of the same size as ``arrays[i]`` . The non-zero mask elements mark the array elements counted in the histogram.
2011-04-25 05:02:14 +08:00
:param hist: Output histogram, which is a dense or sparse ``dims`` -dimensional array.
2011-02-26 19:05:10 +08:00
2011-04-25 05:02:14 +08:00
:param dims: Histogram dimensionality that must be positive and not greater than ``CV_MAX_DIMS`` (=32 in the current OpenCV version).
2011-04-25 05:02:14 +08:00
:param histSize: Array of histogram sizes in each dimension.
2011-04-25 05:02:14 +08:00
:param ranges: Array of the ``dims`` arrays of the histogram bin boundaries in each dimension. When the histogram is uniform ( ``uniform`` =true), then for each dimension ``i`` it is enough to specify the lower (inclusive) boundary :math:`L_0` of the 0-th histogram bin and the upper (exclusive) boundary :math:`U_{\texttt{histSize}[i]-1}` for the last histogram bin ``histSize[i]-1`` . That is, in case of a uniform histogram each of ``ranges[i]`` is an array of 2 elements. When the histogram is not uniform ( ``uniform=false`` ), then each of ``ranges[i]`` contains ``histSize[i]+1`` elements: :math:`L_0, U_0=L_1, U_1=L_2, ..., U_{\texttt{histSize[i]}-2}=L_{\texttt{histSize[i]}-1}, U_{\texttt{histSize[i]}-1}` . The array elements, that are not between :math:`L_0` and :math:`U_{\texttt{histSize[i]}-1}` , are not counted in the histogram.
2011-04-25 05:02:14 +08:00
:param uniform: Flag indicatinfg whether the histogram is uniform or not (see above).
2011-04-25 05:02:14 +08:00
:param accumulate: Accumulation flag. If it is set, the histogram is not cleared in the beginning when it is allocated. This feature enables you to compute a single histogram from several sets of arrays, or to update the histogram in time.
2011-02-26 19:05:10 +08:00
The functions ``calcHist`` calculate the histogram of one or more
2011-04-25 05:02:14 +08:00
arrays. The elements of a tuple used to increment
a histogram bin are taken from the corresponding
input arrays at the same location. The sample below shows how to compute a 2D Hue-Saturation histogram for a color image. ::
#include <cv.h>
#include <highgui.h>
2011-02-26 19:05:10 +08:00
using namespace cv;
2011-02-26 19:05:10 +08:00
int main( int argc, char** argv )
{
Mat src, hsv;
if( argc != 2 || !(src=imread(argv[1], 1)).data )
return -1;
2011-02-26 19:05:10 +08:00
cvtColor(src, hsv, CV_BGR2HSV);
2011-02-26 19:05:10 +08:00
2011-04-25 05:02:14 +08:00
// Quantize the hue to 30 levels
// and the saturation to 32 levels
int hbins = 30, sbins = 32;
int histSize[] = {hbins, sbins};
// hue varies from 0 to 179, see cvtColor
float hranges[] = { 0, 180 };
// saturation varies from 0 (black-gray-white) to
// 255 (pure spectrum color)
float sranges[] = { 0, 256 };
const float* ranges[] = { hranges, sranges };
MatND hist;
// we compute the histogram from the 0-th and 1-st channels
int channels[] = {0, 1};
2011-02-26 19:05:10 +08:00
calcHist( &hsv, 1, channels, Mat(), // do not use mask
hist, 2, histSize, ranges,
true, // the histogram is uniform
false );
double maxVal=0;
minMaxLoc(hist, 0, &maxVal, 0, 0);
2011-02-26 19:05:10 +08:00
int scale = 10;
Mat histImg = Mat::zeros(sbins*scale, hbins*10, CV_8UC3);
2011-02-26 19:05:10 +08:00
for( int h = 0; h < hbins; h++ )
for( int s = 0; s < sbins; s++ )
{
float binVal = hist.at<float>(h, s);
int intensity = cvRound(binVal*255/maxVal);
rectangle( histImg, Point(h*scale, s*scale),
Point( (h+1)*scale - 1, (s+1)*scale - 1),
Scalar::all(intensity),
CV_FILLED );
}
2011-02-26 19:05:10 +08:00
namedWindow( "Source", 1 );
imshow( "Source", src );
2011-02-26 19:05:10 +08:00
namedWindow( "H-S Histogram", 1 );
imshow( "H-S Histogram", histImg );
waitKey();
}
.. index:: calcBackProject
2011-03-09 06:22:24 +08:00
.. _calcBackProject:
calcBackProject
-------------------
.. c:function:: void calcBackProject( const Mat* arrays, int narrays, const int* channels, const MatND\& hist, Mat\& backProject, const float** ranges, double scale=1, bool uniform=true )
.. c:function:: void calcBackProject( const Mat* arrays, int narrays, const int* channels, const SparseMat\& hist, Mat\& backProject, const float** ranges, double scale=1, bool uniform=true )
Calculates the back projection of a histogram.
2011-04-25 05:02:14 +08:00
:param arrays: Source arrays. They all should have the same depth, ``CV_8U`` or ``CV_32F`` , and the same size. Each of them can have an arbitrary number of channels.
2011-02-26 19:05:10 +08:00
2011-04-25 05:02:14 +08:00
:param narrays: Number of source arrays.
2011-04-25 05:02:14 +08:00
:param channels: The list of channels that are used to compute the back projection. The number of channels must match the histogram dimensionality. The first array channels are numerated from 0 to ``arrays[0].channels()-1`` , the second array channels are counted from ``arrays[0].channels()`` to ``arrays[0].channels() + arrays[1].channels()-1`` and so on.
2011-04-25 05:02:14 +08:00
:param hist: Input histogram that can be dense or sparse.
2011-04-25 05:02:14 +08:00
:param backProject: Destination back projection aray that is a single-channel array of the same size and depth as ``arrays[0]`` .
:param ranges: Array of arrays of the histogram bin boundaries in each dimension. See :func:`calcHist` .
:param scale: Optional scale factor for the output back projection.
2011-04-25 05:02:14 +08:00
:param uniform: Flag indicating whether the histogram is uniform or not (see above).
2011-02-26 19:05:10 +08:00
2011-04-25 05:02:14 +08:00
The functions ``calcBackProject`` calculate the back project of the histogram. That is, similarly to ``calcHist`` , at each location ``(x, y)`` the function collects the values from the selected channels in the input images and finds the corresponding histogram bin. But instead of incrementing it, the function reads the bin value, scales it by ``scale`` , and stores in ``backProject(x,y)`` . In terms of statistics, the function computes probability of each element value in respect with the empirical probability distribution represented by the histogram. See how, for example, you can find and track a bright-colored object in a scene:
#.
2011-04-25 05:02:14 +08:00
Before tracking, show the object to the camera so that it covers almost the whole frame. Calculate a hue histogram. The histogram may have strong maximums, corresponding to the dominant colors in the object.
#.
2011-04-25 05:02:14 +08:00
When tracking, calculate a back projection of a hue plane of each input video frame using that pre-computed histogram. Threshold the back projection to suppress weak colors. It may also make sense to suppress pixels with non-sufficient color saturation and too dark or too bright pixels.
#.
Find connected components in the resulting picture and choose, for example, the largest component.
2011-04-25 05:02:14 +08:00
This is an approximate algorithm of the
2011-02-26 19:05:10 +08:00
:func:`CAMShift` color object tracker.
2011-04-25 05:02:14 +08:00
See Also:
2011-02-26 19:05:10 +08:00
:func:`calcHist`
2011-03-09 06:22:24 +08:00
.. index:: compareHist
2011-03-09 06:22:24 +08:00
.. _compareHist:
compareHist
2011-03-09 06:22:24 +08:00
-----------
.. c:function:: double compareHist( const MatND\& H1, const MatND\& H2, int method )
.. c:function:: double compareHist( const SparseMat\& H1, const SparseMat\& H2, int method )
2011-04-25 05:02:14 +08:00
Compares two histograms.
2011-04-25 05:02:14 +08:00
:param H1: The first compared histogram.
2011-04-25 05:02:14 +08:00
:param H2: The second compared histogram of the same size as ``H1`` .
:param method: Comparison method that could be one of the following:
2011-04-25 05:02:14 +08:00
* **CV_COMP_CORREL** Correlation
2011-04-25 05:02:14 +08:00
* **CV_COMP_CHISQR** Chi-Square
2011-04-25 05:02:14 +08:00
* **CV_COMP_INTERSECT** Intersection
2011-04-25 05:02:14 +08:00
* **CV_COMP_BHATTACHARYYA** Bhattacharyya distance
2011-02-26 19:05:10 +08:00
The functions ``compareHist`` compare two dense or two sparse histograms using the specified method:
* Correlation (method=CV\_COMP\_CORREL)
2011-02-26 19:05:10 +08:00
.. math::
2011-02-26 19:05:10 +08:00
d(H_1,H_2) = \frac{\sum_I (H_1(I) - \bar{H_1}) (H_2(I) - \bar{H_2})}{\sqrt{\sum_I(H_1(I) - \bar{H_1})^2 \sum_I(H_2(I) - \bar{H_2})^2}}
where
2011-02-26 19:05:10 +08:00
.. math::
2011-02-26 19:05:10 +08:00
\bar{H_k} = \frac{1}{N} \sum _J H_k(J)
and
2011-04-25 05:02:14 +08:00
:math:`N` is a total number of histogram bins.
* Chi-Square (method=CV\_COMP\_CHISQR)
2011-02-26 19:05:10 +08:00
.. math::
2011-02-26 19:05:10 +08:00
d(H_1,H_2) = \sum _I \frac{\left(H_1(I)-H_2(I)\right)^2}{H_1(I)+H_2(I)}
* Intersection (method=CV\_COMP\_INTERSECT)
2011-02-26 19:05:10 +08:00
.. math::
2011-02-26 19:05:10 +08:00
d(H_1,H_2) = \sum _I \min (H_1(I), H_2(I))
* Bhattacharyya distance (method=CV\_COMP\_BHATTACHARYYA)
2011-02-26 19:05:10 +08:00
.. math::
2011-02-26 19:05:10 +08:00
d(H_1,H_2) = \sqrt{1 - \frac{1}{\sqrt{\bar{H_1} \bar{H_2} N^2}} \sum_I \sqrt{H_1(I) \cdot H_2(I)}}
2011-02-26 19:05:10 +08:00
The function returns
:math:`d(H_1, H_2)` .
2011-04-25 05:02:14 +08:00
While the function works well with 1-, 2-, 3-dimensional dense histograms, it may not be suitable for high-dimensional sparse histograms. In such histograms, because of aliasing and sampling problems, the coordinates of non-zero histogram bins can slightly shift. To compare such histograms or more general sparse configurations of weighted points, consider using the
2011-02-26 19:05:10 +08:00
:func:`calcEMD` function.
2011-02-26 19:05:10 +08:00
.. index:: equalizeHist
2011-03-09 06:22:24 +08:00
.. _equalizeHist:
equalizeHist
----------------
.. c:function:: void equalizeHist( const Mat\& src, Mat\& dst )
Equalizes the histogram of a grayscale image.
2011-04-25 05:02:14 +08:00
:param src: Source 8-bit single channel image.
2011-04-25 05:02:14 +08:00
:param dst: Destination image of the same size and type as ``src`` .
2011-03-03 15:29:55 +08:00
The function equalizes the histogram of the input image using the following algorithm:
#.
2011-04-25 05:02:14 +08:00
Calculate the histogram
:math:`H` for ``src`` .
#.
2011-04-25 05:02:14 +08:00
Normalize the histogram so that the sum of histogram bins is 255.
#.
2011-04-25 05:02:14 +08:00
Compute the integral of the histogram:
2011-02-26 19:05:10 +08:00
.. math::
2011-02-26 19:05:10 +08:00
H'_i = \sum _{0 \le j < i} H(j)
#.
2011-04-25 05:02:14 +08:00
Transform the image using
2011-02-26 19:05:10 +08:00
:math:`H'` as a look-up table:
:math:`\texttt{dst}(x,y) = H'(\texttt{src}(x,y))`
2011-03-03 15:29:55 +08:00
The algorithm normalizes the brightness and increases the contrast of the image.