2 * Copyright 2007 Google Inc.
\r
4 * Licensed under the Apache License, Version 2.0 (the "License");
\r
5 * you may not use this file except in compliance with the License.
\r
6 * You may obtain a copy of the License at
\r
8 * http://www.apache.org/licenses/LICENSE-2.0
\r
10 * Unless required by applicable law or agreed to in writing, software
\r
11 * distributed under the License is distributed on an "AS IS" BASIS,
\r
12 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
\r
13 * See the License for the specific language governing permissions and
\r
14 * limitations under the License.
\r
18 * <p>This class implements a perspective transform in two dimensions. Given four source and four destination
\r
19 * points, it will compute the transformation implied between them. The code is based directly upon section
\r
20 * 3.4.2 of George Wolberg's "Digital Image Warping"; see pages 54-56.</p>
\r
22 * @author srowen@google.com (Sean Owen)
\r
24 package com.google.zxing.qrcode.detector;
\r
26 final class PerspectiveTransform {
\r
28 private final float a11, a12, a13, a21, a22, a23, a31, a32, a33;
\r
30 private PerspectiveTransform(float a11, float a21, float a31,
\r
31 float a12, float a22, float a32,
\r
32 float a13, float a23, float a33) {
\r
44 static PerspectiveTransform quadrilateralToQuadrilateral(float x0, float y0,
\r
48 float x0p, float y0p,
\r
49 float x1p, float y1p,
\r
50 float x2p, float y2p,
\r
51 float x3p, float y3p) {
\r
53 PerspectiveTransform qToS = quadrilateralToSquare(x0, y0, x1, y1, x2, y2, x3, y3);
\r
54 PerspectiveTransform sToQ = squareToQuadrilateral(x0p, y0p, x1p, y1p, x2p, y2p, x3p, y3p);
\r
55 return sToQ.times(qToS);
\r
58 void transformPoints(float[] points) {
\r
59 int max = points.length;
\r
60 float a11 = this.a11;
\r
61 float a12 = this.a12;
\r
62 float a13 = this.a13;
\r
63 float a21 = this.a21;
\r
64 float a22 = this.a22;
\r
65 float a23 = this.a23;
\r
66 float a31 = this.a31;
\r
67 float a32 = this.a32;
\r
68 float a33 = this.a33;
\r
69 for (int i = 0; i < max; i += 2) {
\r
70 float x = points[i];
\r
71 float y = points[i + 1];
\r
72 float denominator = a13 * x + a23 * y + a33;
\r
73 points[i] = (a11 * x + a21 * y + a31) / denominator;
\r
74 points[i + 1] = (a12 * x + a22 * y + a32) / denominator;
\r
78 static PerspectiveTransform squareToQuadrilateral(float x0, float y0,
\r
81 float x3, float y3) {
\r
82 float dy2 = y3 - y2;
\r
83 float dy3 = y0 - y1 + y2 - y3;
\r
84 if (dy2 == 0.0f && dy3 == 0.0f) {
\r
85 return new PerspectiveTransform(x1 - x0, x2 - x1, x0,
\r
86 y1 - y0, y2 - y1, y0,
\r
89 float dx1 = x1 - x2;
\r
90 float dx2 = x3 - x2;
\r
91 float dx3 = x0 - x1 + x2 - x3;
\r
92 float dy1 = y1 - y2;
\r
93 float denominator = dx1 * dy2 - dx2 * dy1;
\r
94 float a13 = (dx3 * dy2 - dx2 * dy3) / denominator;
\r
95 float a23 = (dx1 * dy3 - dx3 * dy1) / denominator;
\r
96 return new PerspectiveTransform(x1 - x0 + a13 * x1, x3 - x0 + a23 * x3, x0,
\r
97 y1 - y0 + a13 * y1, y3 - y0 + a23 * y3, y0,
\r
102 private static PerspectiveTransform quadrilateralToSquare(float x0, float y0,
\r
103 float x1, float y1,
\r
104 float x2, float y2,
\r
105 float x3, float y3) {
\r
106 // Here, the adjoint serves as the inverse:
\r
107 return squareToQuadrilateral(x0, y0, x1, y1, x2, y2, x3, y3).buildAdjoint();
\r
110 PerspectiveTransform buildAdjoint() {
\r
111 // Adjoint is the transpose of the cofactor matrix:
\r
112 return new PerspectiveTransform(a22 * a33 - a23 * a32,
\r
113 a23 * a31 - a21 * a33,
\r
114 a21 * a32 - a22 * a31,
\r
115 a13 * a32 - a12 * a33,
\r
116 a11 * a33 - a13 * a31,
\r
117 a12 * a31 - a11 * a32,
\r
118 a12 * a23 - a13 * a22,
\r
119 a13 * a21 - a11 * a23,
\r
120 a11 * a22 - a12 * a21);
\r
123 PerspectiveTransform times(PerspectiveTransform other) {
\r
124 return new PerspectiveTransform(a11 * other.a11 + a21 * other.a12 + a31 * other.a13,
\r
125 a11 * other.a21 + a21 * other.a22 + a31 * other.a23,
\r
126 a11 * other.a31 + a21 * other.a32 + a31 * other.a33,
\r
127 a12 * other.a11 + a22 * other.a12 + a32 * other.a13,
\r
128 a12 * other.a21 + a22 * other.a22 + a32 * other.a23,
\r
129 a12 * other.a31 + a22 * other.a32 + a32 * other.a33,
\r
130 a13 * other.a11 + a23 * other.a12 + a33 * other.a13,
\r
131 a13 * other.a21 + a23 * other.a22 + a33 * other.a23,
\r
132 a13 * other.a31 + a23 * other.a32 + a33 * other.a33);
\r