Alpha-trimmed 过滤器问题答案

【问题标题】：Alpha-trimmed filter troublesAlpha-trimmed 过滤器问题
【发布时间】：2016-10-18 18:32:30
【问题描述】：

我正在尝试在 openCV 库中制作一个 alphatrimmed 过滤器。我的代码无法正常工作，并且生成的图像在过滤后看起来不像图像。过滤器应按以下方式工作。

在我的示例中选择一些（数组）像素，它是 9 像素的“3x3”窗口。
以递增的方式对它们进行排序。
为 alpha-2 切割我们的“阵列”两侧。
计算剩余像素的算术平均值并将它们插入到适当的位置。

int alphatrimmed(Mat img, int alpha)
{
Mat img9 = img.clone();
const int start = alpha/2 ;
const int end = 9 - (alpha/2);

//going through whole image
for (int i = 1; i < img.rows - 1; i++)
{
for (int j = 1; j < img.cols - 1; j++)
{
    uchar element[9];
    Vec3b element3[9];
    int k = 0;
    int a = 0;
    //selecting elements for window 3x3
        for (int m = i -1; m < i + 2; m++)
        {
            for (int n = j - 1; n < j + 2; n++)
            {
                element3[a] = img.at<Vec3b>(m*img.cols + n);
                a++;
                    for (int c = 0; c < img.channels(); c++)
                    {
                        element[k] += img.at<Vec3b>(m*img.cols + n)[c];
                }
                    k++;
            }
        }
    //comparing and sorting elements in window (uchar element [9])
    for (int b = 0; b < end; b++)
    {
        int min = b;
        for (int d = b + 1; d < 9; d++)
        {
            if (element[d] < element[min])
            {
                min = d;
                const   uchar temp = element[b];
                element[b] = element[min];
                element[min] = temp;
                const   Vec3b temporary = element3[b];
                element3[b] = element3[min];
                element3[min] = temporary;
            }
        }

    }

//  index in resultant image( after alpha-trimmed filter)
int result = (i - 1) * (img.cols - 2) + j - 1;
    for (int l = start ; l < end; l++)
        img9.at<Vec3b>(result) += element3[l];
    img9.at<Vec3b>(result) /= (9 - alpha);
}
}
 namedWindow("AlphaTrimmed Filter", WINDOW_AUTOSIZE);
 imshow("AlphaTrimmed Filter", img9);
return 0;
}

【问题讨论】：

留给opencv秃鹰。
所以你有element[9]和element3，后者硬编码为Vec3b，然后当你知道时使用c < img.channels()它必须是3 ?那是相当令人费解的。排序也很复杂。我明白你为什么“需要”它——你试图对两个数组进行排序——但这只是早期糟糕设计的结果。

标签： c++ opencv

【解决方案1】：

没有实际数据，这有点猜测，但uchar 无法容纳 3 个通道的总和。它以 256 为模（至少在 OpenCV 支持的任何平台上）。

正确的解决方案是std::sort，并为您的Vec3b 提供适当的比较器：

void L1(Vec3b a, Vec3b b) { return a[0]+a[1]+a[2] < b[0]+b[1]+b[2]; }

【讨论】：