opencv 直方图匹配的程序

void HistMatch(IplImage *histimg,double histv[])
{
int bins = 256;
int sizes[] = {bins};
CvHistogram *hist = cvCreateHist(1,sizes,CV_HIST_ARRAY);
cvCalcHist(&histimg,hist);

cvNormalizeHist(hist,1);//函数 cvNormalizeHist 通过缩放来归一化直方块，使得所有块的和等于 factor.

double val_1 = 0.0;
double val_2 = 0.0;
uchar T[256] = {0};
double S[256] = {0};
  double G[256] = {0};

for (int index = 0; index<256; ++index)
{
val_1 += cvQueryHistValue_1D(hist,index);
  val_2 += histv[index];

G[index] = val_2;
S[index] = val_1;
}

double min_val = 0.0;
int PG = 0;

for ( int i = 0; i<256; ++i)
{
min_val = 1.0;
for(int j = 0;j<256; ++j)
{
if( (G[j] - S[i]) < min_val && (G[j] - S[i]) >= 0)
{
  min_val = (G[j] - S[i]);
PG = j;
}
}
T[i] = (uchar)PG;
}
   
uchar *p = NULL;
for (int x = 0; x<histimg->height;++x)
{ 
p = (uchar*)(histimg->imageData + histimg->widthStep*x);
for (int y = 0; y<histimg->width;++y)
{
p[y] = T[p[y]];
}
}
}
图片直方图匹配的程序，哪位大虾给解释下，用的什么原理进行的匹配？
尤其是这一段
uchar *p = NULL;
for (int x = 0; x<histimg->height;++x)
{ 
p = (uchar*)(histimg->imageData + histimg->widthStep*x);
for (int y = 0; y<histimg->width;++y)
{
p[y] = T[p[y]];
}
}

中的 p[y] = T[p[y]];是什么意思？？
最好详细点~~

谢谢了~~

结贴可以走人了

看代码是要将原图像直方图以给定的直方图为标准作变换,使两图像的直方图相同和近似，从而使两幅图像具有类似的色调和反差。所要做的就是要修改某些灰度值了。
p[y] = T[p[y]]
是将要匹配的图像的(x,y)位置的像素灰度值改变为要改变的数值，T[] 在这里是看成一个查找表，有256个灰度级，
例如：
比如说原图像g(x,y)位置的灰度为102的话，目标是要改成的灰度为128的话，那么T[102]在前面几个步骤中就计算出来了是128，所以g(x,y) = T[g(x,y)] = T[102] = 128. 可以看看查找表LUT的使用。