[AVR] Suche Formel(n) zur Kissenentzerrung

[Juergen (Gast)]

Hallo.

Ich suche für mein aktuelles Projekt Formeln um eine Kissenentzerrung zu 
machen.
Es handelt sich um einen Laser und 2 Galvanometer im XY-Betrieb.
Leider ist das 'Koordinatensystem', aufgrund des geringen Abstands zur 
Projektionsfläche, stark Kissen- bzw. Tonnenförmig verzerrt.

Ich wollte es erst mit einer Koordinatenliste machen, mit deren Hilfe 
dann die Werte interpoliert werden, aber da komme ich schnell auf 
mehrere tausend Werte, die ich alle einyeln ausmessen muss. Da habe ich 
nicht so den Antrieb zu. ;)

Kann mir da vielleicht jemand weiter helfen?

[Joe F. (easylife)]

http://marcodiiga.github.io/radial-lens-undistortion-filtering

[Juergen (Gast)]

Herzlichen Dank. Werde es sofort versuchen.

[Juergen (Gast)]

Kurze Rückmeldung:
Das klappt einwandfrei mit dem im Link vorhandenm Beispielcode.

Für meine Zwecke reichen die beiden auskommentierten Zeilen. Die 4 
Zeilen mit x2 und x3 habe ioch auskommmentiert und die beiden mit x2 
reingenommen.
Falls der Link mal nicht mehr geht, hier der Code, so wie ich ihn jetzt 
drin habe:

// Normalize the u,v coordinates in the range [-1;+1]
volatile float x0 = ((2.0 * x_32) - SOURCE_XSIZE) / SOURCE_XSIZE;
volatile float y0 = ((2.0 * y_32) - SOURCE_YSIZE) / SOURCE_YSIZE;

// Calculate l2 norm
float r = x0 * x0 + y0 * y0;

// Calculate the deflated or inflated new coordinate (reverse transform)
//         float x1 = x0 / (1.0 - (PARSER_ALPHAX * r));
//         float y1 = y0 / (1.0 - (PARSER_ALPHAY * r));
//         float x2 = x0 / (1.0 - (PARSER_ALPHAX * ((x1 * x1) + (y1 * y1))));
//         float y2 = y0 / (1.0 - (PARSER_ALPHAY * ((x1 * x1) + (y1 * y1))));

// Forward transform
float x2 = x0 * (1.0 - PARSER_ALPHAX * r);
float y2 = y0 * (1.0 - PARSER_ALPHAY * r);

// De-normalize to the original range
float k = (x2 + 1.0) * DEST_XSIZE / 2.0;
float l = (y2 + 1.0) * DEST_YSIZE / 2.0;

Kurze Erklärung hierzu:
SOURCE_nSIZE und DEST_nSIZE sind wg. unterschiedlichen Größen des 
Koordinatensystems. Einmal wird in mm eingelesen und als Ausgabe habe 
ich einen 16-Bit DAC...
x_32 und y_32 sind die Werte aus der Datei als 32-Bit Variable 
umgerechnet. Zur Info: 250,0000mm (im maximum) entsprechen, der 
Einfachheit halber, 2500000.
k und l sind das Ergebnis.

Hoffe mich klar ausgedrüückt zu haben. ;)