Формулы для искажения Баррель / Подушечка

Question

Не могу понять, как получить (x ', y') оригинала (x, y) на изображении для искажения Баррель / Подушечка.

Answer 1

Раздел 2 из , этот документ объясняет преобразование.В основном:

Здесь я привел пример в Mathematica :

Answer 2

Аппроксимация модели полиномиального радиального искажения, которую вы можете найти в Fitzgibbon, 2001, равна

, где rd и ru - расстояния от центраискажения.Это также используется для фильтрации искажений изображения широкоугольной камеры для целей компьютерного зрения и обработки изображений.

Более подробное объяснение принципа и кода шейдера для реализации фильтрации без искажения (а также прямое преобразование) здесь: http://marcodiiga.github.io/radial-lens-undistortion-filtering

Я также публикую статьи, на которые вам следует обратить внимание, если вы хотите узнать математические детали для метода, который я опубликовал

• Чжан З. (1999).Гибкая калибровка камеры путем просмотра плоскости с неизвестной ориентацией
• Andrew W. Fitzgibbon (2001).Одновременная линейная оценка геометрии нескольких видов и искажения линзы

Answer 3

искажение простого ствола \ подушечка в opencv c ++

IplImage* barrel_pincusion_dist(IplImage* img, double Cx,double Cy,double kx,double ky)
{
    IplImage* mapx = cvCreateImage( cvGetSize(img), IPL_DEPTH_32F, 1 );
    IplImage* mapy = cvCreateImage( cvGetSize(img), IPL_DEPTH_32F, 1 );

    int w= img->width;
    int h= img->height;

    float* pbuf = (float*)mapx->imageData;
    for (int y = 0; y < h; y++)
    {
        for (int x = 0; x < w; x++)
        {         
            float u= Cx+(x-Cx)*(1+kx*((x-Cx)*(x-Cx)+(y-Cy)*(y-Cy)));
            *pbuf = u;
            ++pbuf;
        }
    }

    pbuf = (float*)mapy->imageData;
    for (int y = 0;y < h; y++)
    {
        for (int x = 0; x < w; x++) 
        {
            *pbuf = Cy+(y-Cy)*(1+ky*((x-Cx)*(x-Cx)+(y-Cy)*(y-Cy)));
            ++pbuf;
        }
    }

    /*float* pbuf = (float*)mapx->imageData;
    for (int y = 0; y < h; y++)
    {
        int ty= y-Cy;
        for (int x = 0; x < w; x++)
        {
            int tx= x-Cx;
            int rt= tx*tx+ty*ty;

            *pbuf = (float)(tx*(1+kx*rt)+Cx);
            ++pbuf;
        }
    }

    pbuf = (float*)mapy->imageData;
    for (int y = 0;y < h; y++)
    {
        int ty= y-Cy;
        for (int x = 0; x < w; x++) 
        {
            int tx= x-Cx;
            int rt= tx*tx+ty*ty;

            *pbuf = (float)(ty*(1+ky*rt)+Cy);
            ++pbuf;
        }
    }*/

    IplImage* temp = cvCloneImage(img);
    cvRemap( temp, img, mapx, mapy ); 
    cvReleaseImage(&temp);
    cvReleaseImage(&mapx);
    cvReleaseImage(&mapy);

    return img;
}

более сложная форма http://opencv.willowgarage.com/documentation/camera_calibration_and_3d_reconstruction.html

Answer 4

Согласно Википедии, также может быть r к степени 4. Знаки двух констант (для терминов от r до 2 и от r до 4) могут быть противоположными, что приводит к искажению руля, когда центр изображения имеет бочкообразное искажение, а край имеет искажение подушкообразной булавкой, создавая прямые линии, напоминающие усы на руле .