Преобразование цвета в оттенки серого - PullRequest
0 голосов
/ 09 октября 2019

Я использую C ++ openCV программу для первых принципов. Разработка алгоритма для обнаружения объектов изображения HDL (Verilog). Наконец-то мне удалось довести версию HDL до такой степени, что ее можно обнаружить. Для проверки двух, оба должны иметь идентичный вывод. Я обнаружил, что есть тонкие отличия, которые, как мне кажется, вносят вклад в преобразование цвета изображения в оттенки серого в смещение зеленого цвета openCV. Сглаженное изображение в целом более яркое в методе openCV C ++. Глядя на метод rgb2gray, кажется, что openCV использовал смещение, т.е. (R X + G Y + B * Z) / 3, в то время как в HDL я использовал (R + G + B) / 3 в качестве Iтребует, чтобы он завершил фильтры Гаусса, Собеля и Канни. Человеческая визуализация является вторичной, а умножение на не-int нежелательно.

Существует ли стандартное линейное преобразование в оттенках серого для преобразования или средство для отмены существующего метода? ...

int main()
{
            int thold = 15;

            clock_t start;
            double duration;
            const int sobelX[3][3] = { {-1, 0, 1}, {-2, 0, 2}, {-1, 0, 1} };  //Where origionally floats in python
            const int sobelY[3][3] = { {-1, -2, -1}, {0, 0, 0}, {1, 2, 1} }; //Where origionally floats in python
            const int kernel[5][5] = { {1,6,12,6,1},
                                                                                                {6,42,79,42,6},
                                                                                                                        {12,79,148,79,12},
                                                                                                                        {6,42,79,42,6},
                                                                                                                        {1,6,12,6,1} };// 1/732
            // Above normalised kernal for smoothing,  see origional python script for method 
            start = std::clock();
            int height, width, intPixel, tSx, tSy, tS, dirE, dirEE, maxDir, curPoint, contDirection, cannyImgPix, nd, tl, tm, tr, mr, br, bm, bl, ml = 0;
            int contNum = 128;
            int contPixCount = 0;
            int curContNum = 0;
            int contPlace = 0;
            int oldContPlace = 0;
            int g = 0;
            bool maxPoint;
            struct pixel {
                        int number;
                        int h;
                        int w;

            };
            std::vector<pixel> contourList;
            //double floatPixel = 0.0;
            int kernalCumulator = 0;
            const int mp = 3;
        //  Scalar color(0, 0, 255);
            //          duration = ((clock()) - start) / (double)CLOCKS_PER_SEC;
            //          start = clock();
            //          cout << "Start image in" << duration << '\n';
            //          Mat dst;
            Mat rawImg = imread("C:\\Users\\&&&\\Documents\\pycode\\paddedGS.png",0);
            printf("%d",rawImg.type());

//          Mat rawImg = imread("C:\\Users\\&&&\\Documents\\openCV_Master\\openCVexample\\openCVexample\\brace200.jpg ", 0);
            height = rawImg.rows;
            width = rawImg.cols;
            cout << "Height of image " << height << '\n';
            cout << "Width of image " << width << '\n';
            Mat filteredImg = Mat::zeros(height, width, CV_8U);
            printf("%d", filteredImg.type());
            Mat sobelImg = Mat::zeros(height, width, CV_8U);
            Mat directionImg = Mat::zeros(height, width, CV_8U);
            Mat cannyImg = Mat::zeros(height, width, CV_8U);
            Mat contourImg = Mat::zeros(height, width, CV_16U);

//          rawImg.convertTo(rawImg, CV_8UC1);

            duration = ((clock()) - start) / (double)CLOCKS_PER_SEC;
            start = clock();
            cout << "Start image in" << duration << '\n';
            // Loop to threshold already grayscaled image           
            /*
            for (int h = 0; h < (height); h++)
            {
                        for (int w = 0; w < (width); w++)
                        {
                                    g = (int)rawImg.at<uchar>(h, w,0);
                                    cout << g << "g";
                                    g+= (int)rawImg.at<uchar>(h, w, 1);
                                    cout << g << "g";
                                    g+= (int)rawImg.at<uchar>(h, w, 2);
                                    cout << g << "g";
                                    g = g/3;
                                    rawGImg.at<uchar>(h,w) = g;
                        }
            }

            */
            //          imshow("thresholded Image", rawImg);
            //          waitKey();
            // Loop to smooth using Gausian 5 x 5 kernal

//          imshow("raw Image", rawImg);


            for (int h = 3; h < (height - 3); h++)
            {
                        for (int w = 3; w < (width - 3); w++)
                        {
                                    if (rawImg.at<uchar>(h, w) >=6 )//Thresholding included
                                    {
                                                for (int xk = 0; xk < 5; xk++)
                                                {
                                                            for (int yk = 0; yk < 5; yk++)
                                                            {
                                                                        intPixel = rawImg.at<uchar>((h + (xk - mp)), (w + (yk - mp)));
                                                                        kernalCumulator += intPixel*(kernel[xk][yk]);//Mutiplier required as rounding is making number go above 255,  better solution?
                                                            }
                                                }
                                    }
                                    else
                                                kernalCumulator = 0;

                                    kernalCumulator = kernalCumulator / 732;
                                    if (kernalCumulator < 0 || kernalCumulator > 255)
                                    {
        //                                      cout << "kernal Value: " << kernalCumulator;
            //                                  cout << " intPixel:" << intPixel << '\n';
                                    }
                                    filteredImg.at<uchar>(h, w) = (uchar)kernalCumulator;
                                    kernalCumulator = 0;
                        }
            }

Ответы [ 2 ]

0 голосов
/ 09 октября 2019

Согласно совету Paul92 выше

cv::Mat linearRgbToGray(const cv::Mat &color) {
            cv::Mat gray(color.size(), CV_8UC1);
            for (int i = 0; i < color.rows; i++)
                        for (int j = 0; j < color.cols; j++)
                                    gray.at<uchar>(i, j) = ((color.at<cv::Vec3b>(i, j)[0] + color.at<cv::Vec3b>(i, j)[1] + color.at<cv::Vec3b>(i, j)[2]) / 3);
            return gray;
}

Приведенный выше код сработал и преодолел предельные ошибки, с которыми я столкнулся ранее. Спасибо, Роб.

0 голосов
/ 09 октября 2019

Наше зрение не воспринимает яркость линейно, поэтому для обычных приложений имеет смысл использовать какое-то преобразование, которое пытается имитировать человеческое восприятие.

Для вашего приложения у вас есть 2 варианта: либо использоватьаналогичное преобразование в HDL (которое может быть непростым или нежелательным) или создание настраиваемого rgb в оттенки серого для OpenCV, в котором используется то же преобразование, что и вы.

Короткий фрагмент (более похожий на псевдокодчтобы выяснить детали) для этого было бы что-то вроде:

cv::Mat linearRgbToGray(const cv::Mat &color) {
    cv::Mat gray(color.size(), CV_8UC1);
    for (int i = 0; i < color.rows; i++)
        for (int j = 0; j < color.cols; j++)
           gray.at(i, j) = (color.at(i, j)[0] + color.at(i, j)[1] + color.at(i, j)[2]) / 3;
}
Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...