opencv实现matlab的median()函数

总结一下之前在做项目时需要用到把matlab的一些代码转换成c/C++代码，下面是实现matlab中的median()函数的代码，供大家借鉴。
matlab中median（）函数的实现。
int cmp(const void a, const void b)//排序函数，
{
return (float )a>(float)b ? 1 : -1;
}
///
int median(vector<IplImage> LR_I, IplImageresult)
{
int LR_value = LR_I.size();
float last_gray = 0;
float p[20] = { 0 };
for (int i = 0;i< LR_I[0]->width; i++)
{
for (int j = 0;j< LR_I[0]->height; j++)
{
for (int k = 0; k < LR_value; k++)
{
float gray = cvGetReal2D(LR_I[k], i, j);
p[k] = gray;
}
qsort(p, LR_value, sizeof(p[0]), cmp);
if (LR_value%2==0)
{
int m = LR_value / 2;
int n = (LR_value / 2)-1;
last_gray = (p[m] + p[n]) / 2.0;
}
else
{
int m = LR_value / 2;
last_gray = p[m];
}
cvSetReal2D(result, i, j, last_gray);
}
}
return 0;
}

这里需要注意的是排序函数qsort(),在对float或是double型进行排序时，要用三目运算符，
写返回值，return *(int )a >(int *)b?1:-1; 因为cmp排序函数返回值为int型，如果两个数很相近，则会认为返回值为0，无法正确排序。
//整型数组排序
int num[100];
int cmp ( const void *a , const void *b )
{
return *(int *)a - *(int *)b;
}
qsort(num,100,sizeof(num[0]),cmp);
///double型数组排序
Double num[100];
int cmp ( const void a , const void b )
{
return (double)a >(double)b? 1:-1;
}
qsort(num,100,sizeof(num[0]),cmp);