Проблема с другим результатом преобразования DCT на изображении при преобразовании C ++ в C # - PullRequest
Новая
       29

Проблема с другим результатом преобразования DCT на изображении при преобразовании C ++ в C #

0 голосов
/ 18 сентября 2018

Я попытался написать приложение, которое будет выполнять преобразование DCT-2 для фрагмента изображения, а затем преобразовывать обратно с помощью обратного DCT-2.Я нашел код в c ++ / opencv, который я пытался преобразовать в C #, но в какой-то момент у меня был другой результат.Вот код, который я пытался преобразовать: 

for (unsigned v = 0; v < BLOCK_SIZE; ++v)
{   
    for (unsigned u = 0; u < BLOCK_SIZE; ++u)
    {   
        const double cu = (u == 0) ? 1.0 / sqrt(2) : 1.0;
        const double cv = (v == 0) ? 1.0 / sqrt(2) : 1.0;
        double dctCoeff = 0;

        for (unsigned y = 0; y < BLOCK_SIZE; ++y)
        {   
            for (unsigned x = 0; x < BLOCK_SIZE; ++x)
            {   
                double uCosFactor = cos((double)(2 * x + 1) * M_PI * (double)u / (2 * (double) BLOCK_SIZE));
                double vCosFactor = cos((double)(2 * y + 1) * M_PI * (double)v / (2 * (double) BLOCK_SIZE));
                double pixel = (double)(lenaNoseGrey.at<unsigned char>(cv::Point(x,y)));
                dctCoeff += pixel * uCosFactor * vCosFactor;
            }
        }
        dctCoeff *= (2 / (double) BLOCK_SIZE) * cu * cv;
        lenaNoseDct.at<double>(cv::Point(u,v)) = dctCoeff;
    }
}

А вот мой: 

for (int v = 0; v < BLOCK_SIZE; ++v)
{   
    for (int u = 0; u < BLOCK_SIZE; ++u)
    {   

        double cu = (u == 0) ? 1.0 / Math.Sqrt(2) : 1.0;
        double cv = (v == 0) ? 1.0 / Math.Sqrt(2) : 1.0;
        double dctCoeff = 0;
        double dctCoeffAlpha = 0;

        for (int y1 = 0; y1 < BLOCK_SIZE; ++y1)
        {   
            for (int x1 = 0; x1 < BLOCK_SIZE; ++x1)
            {   
                double uCosFactor = Math.Cos((2 * x1 + 1) * Math.PI * u / (2 * (double)BLOCK_SIZE));
                double vCosFactor = Math.Cos((2 * y1 + 1) * Math.PI * v / (2 * (double)BLOCK_SIZE));
                double pixel = (double)bitmapaWy1.GetPixel((x1 + 284), (y1 + 313)).R;
                double pixelalpha = (double)bitmapaWy1.GetPixel((x1 + 284), (y1 + 313)).A;
                dctCoeff += pixel * uCosFactor * vCosFactor;
                //dctCoeffAlpha += pixelalpha * uCosFactor * vCosFactor;
                dctCoeffAlpha = pixelalpha;

            }
        }
        dctCoeffAlpha *= (2 / (double)BLOCK_SIZE) * cu * cv;
        dctCoeff *= (2 / (double)BLOCK_SIZE) * cu * cv;
        macierz[u, v] = dctCoeff;

    }
}

У меня другой результат в моей матрице, но когда я конвертирую матрицу из кода C ++, описанного вышес моим обратным кодом, это работает хорошо.Можете ли вы найти, что я сделал не так?Единственное отличие, которое я могу заметить, - это метод getpixel, но он выполняется для точно такого же изображения в градациях серого и фрагмента изображения, который был выполнен в коде C ++.

Когда я снова обращаю свою матрицу к изображению, я могу видеть изображение, но в нем много случайных пикселей, которые слишком белые или слишком черные.

Ответы [ 2 ]

0 голосов
/ 19 сентября 2018

Проблема решена, проблема в том, что слишком много кода, фактический код, который я опубликовал, работает, проблема была в том, что у меня был код, который менял dctcoeff на байт, и это было внутри цикла, поэтому значение dct не сбрасывалось в 0;

int BLOCK_SIZE = 16;
        double[,] macierz = new double[16, 16];

        for (int v = 0; v < BLOCK_SIZE; ++v)
        {
            for (int u = 0; u < BLOCK_SIZE; ++u)
            {

                 double cu = (u == 0) ? 1.0 / Math.Sqrt(2) : 1.0;
                double cv = (v == 0) ? 1.0 / Math.Sqrt(2) : 1.0;
                double dctCoeff = 0;


                for (int y1 = 0; y1 < BLOCK_SIZE; ++y1)
                {
                    for (int x1 = 0; x1 < BLOCK_SIZE; ++x1)
                    {
                        double uCosFactor = Math.Cos((2 * x1 + 1) * Math.PI * u / (2 * (double)BLOCK_SIZE));
                        double vCosFactor = Math.Cos((2 * y1 + 1) * Math.PI * v / (2 * (double)BLOCK_SIZE));
                        double pixel = bitmapaWy1.GetPixel((x1 + 284), (y1 + 313)).R;

                        dctCoeff += pixel * uCosFactor * vCosFactor;
                        //dctCoeffAlpha += pixelalpha * uCosFactor * vCosFactor;


                    }
                }

                dctCoeff *= (2 / (double)BLOCK_SIZE) * cu * cv;
                macierz[u, v] = dctCoeff;

            }
        }

Этот код работает как брелок.Я надеюсь, что какая-то хорошая душа будет использовать его, потому что я потратил около 20 часов, чтобы заставить его работать.

0 голосов
/ 19 сентября 2018

Это явно проблема точности.

Скорее всего, проблема в следующей строке: 

 macierz[u, v] = dctCoeff;

У меня также есть сомнения в: 

            double pixel = (double)bitmapaWy1.GetPixel((x1 + 284), (y1 + 313)).R;
            double pixelalpha = (double)bitmapaWy1.GetPixel((x1 + 284), (y1 + 313)).A;
            dctCoeff += pixel * uCosFactor * vCosFactor;
            //dctCoeffAlpha += pixelalpha * uCosFactor * vCosFactor;
            dctCoeffAlpha = pixelalpha;
...