Извлечение четырехстороннего изображения в прямоугольник

A photo

ОБНОВЛЕНИЕ BOUNTY

Следуя в Denis's, вот как использовать threeblindmiceandamonkey код:

// the destination rect is our 'in' quad
int dw = 300, dh = 250;
double in[4][4] = {{0,0},{dw,0},{dw,dh},{0,dh}};
    // the quad in the source image is our 'out'
double out[4][5] = {{171,72},{331,93},{333,188},{177,210}};
double homo[3][6];
const int ret = mapQuadToQuad(in,out,homo);
    // homo can be used for calculating the x,y of any destination point
// in the source, e.g.
for(int i=0; i<4; i++) {
    double p1[3] = {out[i][0],out[i][7],1};
    double p2[3];
    transformMatrix(p1,p2,homo);
    p2[0] /= p2[2]; // x
    p2[1] /= p2[2]; // y
    printf("\t%2.2f\t%2.2f\n",p2[0],p2[1]);
}

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

for(int y=0; y<dh; y++) {
    for(int x=0; x<dw; x++) {
        // calc the four corners in source for this
        // destination pixel, and mix

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

СПРАВОЧНЫЙ ВОПРОС

В изображении вверху знак на боковой стороне фургона не обращен к камере. Я хочу рассчитать, насколько это возможно, с пикселями, которые у меня есть, как бы это было на лице.

Я знаю угловые координаты квадрата в изображении и размер прямоугольника назначения.

Я предполагаю, что это какой-то цикл через оси x и y, делающий линию Bresenham на обоих измерениях одновременно с каким-то перемешиванием, поскольку пиксели в изображениях источника и адресата перекрываются - некоторые смешения некоторых субпикселей

Какие подходы существуют и как вы смешиваете пиксели?

Есть ли стандартный подход для этого?

ОТВЕТЫ

Ответ 1

Посмотрите "quad to quad" transform, например. threeblindmiceandamonkey.
Преобразование 3x3 на 2d однородных координатах может преобразовывать любые 4 точки (квадрат) к любому другому квадранту; наоборот, любые fromquad и toquad, такие как углы вашего грузовика и целевой прямоугольник, дают преобразование 3 х 3.

Qt имеет quadToQuad и может трансформировать pixmaps с ним, но я думаю, у вас нет Qt?
Добавлено 10Jun: от labs.trolltech.com/page/Graphics/Examples там хорошая демоверсия, в которой квад-к-четности пиксельная карта перемещается по углам:

alt text

Добавлено 11Jun: @Will, здесь translate.h в Python (что вы знаете немного? ""... "" " являются многострочными комментариями.)
perstrans() - ключ; надеюсь, что имеет смысл, если не спросить.

Кстати, вы могли отображать пиксели по одному, mapQuadToQuad (целевой прямоугольник, ориг-квад), но без пиксельной интерполяции он будет выглядеть ужасно; OpenCV делает все это.

#!/usr/bin/env python
""" square <-> quad maps
    from http://threeblindmiceandamonkey.com/?p=16 matrix.h
"""

from __future__ import division
import numpy as np

__date__ = "2010-06-11 jun denis"

def det2(a, b, c, d):
    return a*d - b*c

def mapSquareToQuad( quad ):  # [4][2]
    SQ = np.zeros((3,3))
    px = quad[0,0] - quad[1,0] + quad[2,0] - quad[3,0]
    py = quad[0,1] - quad[1,1] + quad[2,1] - quad[3,1]
    if abs(px) < 1e-10 and abs(py) < 1e-10:
        SQ[0,0] = quad[1,0] - quad[0,0]
        SQ[1,0] = quad[2,0] - quad[1,0]
        SQ[2,0] = quad[0,0]
        SQ[0,1] = quad[1,1] - quad[0,1]
        SQ[1,1] = quad[2,1] - quad[1,1]
        SQ[2,1] = quad[0,1]
        SQ[0,2] = 0.
        SQ[1,2] = 0.
        SQ[2,2] = 1.
        return SQ
    else:
        dx1 = quad[1,0] - quad[2,0]
        dx2 = quad[3,0] - quad[2,0]
        dy1 = quad[1,1] - quad[2,1]
        dy2 = quad[3,1] - quad[2,1]
        det = det2(dx1,dx2, dy1,dy2)
        if det == 0.:
            return None
        SQ[0,2] = det2(px,dx2, py,dy2) / det
        SQ[1,2] = det2(dx1,px, dy1,py) / det
        SQ[2,2] = 1.
        SQ[0,0] = quad[1,0] - quad[0,0] + SQ[0,2]*quad[1,0]
        SQ[1,0] = quad[3,0] - quad[0,0] + SQ[1,2]*quad[3,0]
        SQ[2,0] = quad[0,0]
        SQ[0,1] = quad[1,1] - quad[0,1] + SQ[0,2]*quad[1,1]
        SQ[1,1] = quad[3,1] - quad[0,1] + SQ[1,2]*quad[3,1]
        SQ[2,1] = quad[0,1]
        return SQ

#...............................................................................
def mapQuadToSquare( quad ):
    return np.linalg.inv( mapSquareToQuad( quad ))

def mapQuadToQuad( a, b ):
    return np.dot( mapQuadToSquare(a), mapSquareToQuad(b) )

def perstrans( X, t ):
    """ perspective transform X Nx2, t 3x3:
        [x0 y0 1] t = [a0 b0 w0] -> [a0/w0 b0/w0]
        [x1 y1 1] t = [a1 b1 w1] -> [a1/w1 b1/w1]
        ...
    """
    x1 = np.vstack(( X.T, np.ones(len(X)) ))
    y = np.dot( t.T, x1 )
    return (y[:-1] / y[-1]) .T

#...............................................................................
if __name__ == "__main__":
    np.set_printoptions( 2, threshold=100, suppress=True )  # .2f

    sq = np.array([[0,0], [1,0], [1,1], [0,1]])
    quad = np.array([[171, 72], [331, 93], [333, 188], [177, 210]])
    print "quad:", quad
    print "square to quad:", perstrans( sq, mapSquareToQuad(quad) )
    print "quad to square:", perstrans( quad, mapQuadToSquare(quad) )

    dw, dh = 300, 250
    rect = np.array([[0, 0], [dw, 0], [dw, dh], [0, dh]])
    quadquad = mapQuadToQuad( quad, rect )
    print "quad to quad transform:", quadquad
    print "quad to rect:", perstrans( quad, quadquad )
"""
quad: [[171  72]
 [331  93]
 [333 188]
 [177 210]]
square to quad: [[ 171.   72.]
 [ 331.   93.]
 [ 333.  188.]
 [ 177.  210.]]
quad to square: [[-0.  0.]
 [ 1.  0.]
 [ 1.  1.]
 [ 0.  1.]]
quad to quad transform: [[   1.29   -0.23   -0.  ]
 [  -0.06    1.79   -0.  ]
 [-217.24  -88.54    1.34]]
quad to rect: [[   0.    0.]
 [ 300.    0.]
 [ 300.  250.]
 [   0.  250.]]
"""

Ответ 2

То, что вы хотите, называется плоской ректификацией, и это не все так просто, я боюсь. Что вам нужно сделать, это восстановить homography, который отображает это наклонное изображение стороны фургона на переднюю сторону. У Photoshop/и т.д. Есть инструменты, чтобы сделать это для вас, учитывая некоторые контрольные точки; если вы хотите реализовать его для себя, вам придется начать понимать компьютерное зрение.

Изменить - ОК, вот вам: Python script, чтобы выполнить деформирование, используя OpenCV библиотека, которая имеет удобные функции для расчета гомографии и деформации изображения для вас:

import cv

def warpImage(image, corners, target):
    mat = cv.CreateMat(3, 3, cv.CV_32F)
    cv.GetPerspectiveTransform(corners, target, mat)
    out = cv.CreateMat(height, width, cv.CV_8UC3)
    cv.WarpPerspective(image, out, mat, cv.CV_INTER_CUBIC)
    return out

if __name__ == '__main__':
    width, height = 400, 250
    corners = [(171,72),(331,93),(333,188),(177,210)]
    target = [(0,0),(width,0),(width,height),(0,height)]
    image = cv.LoadImageM('fries.jpg')
    out = warpImage(image, corners, target)
    cv.SaveImage('fries_warped.jpg', out)

И результат:
Warped image

OpenCV также имеет привязки C и С++, или вы можете использовать EmguCV для оболочки .NET; API достаточно согласован во всех языках, поэтому вы можете реплицировать это в зависимости от того, какой язык подходит вашей фантазии.

Ответ 3

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