Mercurial > hg > orthanc-stone
annotate Framework/Toolbox/AffineTransform2D.cpp @ 574:911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
author | Sebastien Jodogne <s.jodogne@gmail.com> |
---|---|
date | Thu, 18 Apr 2019 18:06:29 +0200 |
parents | b70e9be013e4 |
children | 919226caca82 |
rev | line source |
---|---|
409 | 1 /** |
2 * Stone of Orthanc | |
3 * Copyright (C) 2012-2016 Sebastien Jodogne, Medical Physics | |
4 * Department, University Hospital of Liege, Belgium | |
439 | 5 * Copyright (C) 2017-2019 Osimis S.A., Belgium |
409 | 6 * |
7 * This program is free software: you can redistribute it and/or | |
8 * modify it under the terms of the GNU Affero General Public License | |
9 * as published by the Free Software Foundation, either version 3 of | |
10 * the License, or (at your option) any later version. | |
11 * | |
12 * This program is distributed in the hope that it will be useful, but | |
13 * WITHOUT ANY WARRANTY; without even the implied warranty of | |
14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU | |
15 * Affero General Public License for more details. | |
16 * | |
17 * You should have received a copy of the GNU Affero General Public License | |
18 * along with this program. If not, see <http://www.gnu.org/licenses/>. | |
19 **/ | |
20 | |
21 | |
22 #include "AffineTransform2D.h" | |
23 | |
24 #include "ImageGeometry.h" | |
25 | |
26 #include <Core/Logging.h> | |
27 #include <Core/OrthancException.h> | |
28 | |
29 namespace OrthancStone | |
30 { | |
31 AffineTransform2D::AffineTransform2D() : | |
32 matrix_(LinearAlgebra::IdentityMatrix(3)) | |
33 { | |
34 } | |
35 | |
36 | |
37 AffineTransform2D::AffineTransform2D(const Matrix& m) | |
38 { | |
39 if (m.size1() != 3 || | |
40 m.size2() != 3) | |
41 { | |
42 throw Orthanc::OrthancException(Orthanc::ErrorCode_IncompatibleImageSize); | |
43 } | |
44 | |
45 if (!LinearAlgebra::IsCloseToZero(m(2, 0)) || | |
46 !LinearAlgebra::IsCloseToZero(m(2, 1)) || | |
47 LinearAlgebra::IsCloseToZero(m(2, 2))) | |
48 { | |
49 LOG(ERROR) << "Cannot setup an AffineTransform2D with perspective effects"; | |
50 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange); | |
51 } | |
52 | |
53 matrix_ = m / m(2, 2); | |
54 } | |
55 | |
56 | |
57 void AffineTransform2D::Apply(double& x /* inout */, | |
58 double& y /* inout */) const | |
59 { | |
60 Vector p; | |
61 LinearAlgebra::AssignVector(p, x, y, 1); | |
62 | |
63 Vector q = LinearAlgebra::Product(matrix_, p); | |
64 | |
65 if (!LinearAlgebra::IsNear(q[2], 1.0)) | |
66 { | |
67 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError); | |
68 } | |
69 else | |
70 { | |
71 x = q[0]; | |
72 y = q[1]; | |
73 } | |
74 } | |
75 | |
76 | |
77 void AffineTransform2D::Apply(Orthanc::ImageAccessor& target, | |
78 const Orthanc::ImageAccessor& source, | |
79 ImageInterpolation interpolation, | |
80 bool clear) const | |
81 { | |
82 assert(LinearAlgebra::IsNear(matrix_(2, 0), 0) && | |
83 LinearAlgebra::IsNear(matrix_(2, 1), 0) && | |
84 LinearAlgebra::IsNear(matrix_(2, 2), 1)); | |
85 | |
86 ApplyAffineTransform(target, source, | |
87 matrix_(0, 0), matrix_(0, 1), matrix_(0, 2), | |
88 matrix_(1, 0), matrix_(1, 1), matrix_(1, 2), | |
89 interpolation, clear); | |
90 } | |
91 | |
92 | |
574
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
93 void AffineTransform2D::ConvertToOpenGLMatrix(float target[16], |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
94 unsigned int canvasWidth, |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
95 unsigned int canvasHeight) const |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
96 { |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
97 const AffineTransform2D t = AffineTransform2D::Combine( |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
98 CreateOpenGLClipspace(canvasWidth, canvasHeight), *this); |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
99 |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
100 const Matrix source = t.GetHomogeneousMatrix(); |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
101 |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
102 if (source.size1() != 3 || |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
103 source.size2() != 3) |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
104 { |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
105 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError); |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
106 } |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
107 |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
108 // "z" must be in the [-1,1] range, otherwise the texture does not show up |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
109 float z = 0; |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
110 |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
111 // Embed the 3x3 affine transform of the 2D plane into a 4x4 |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
112 // matrix (3D) for OpenGL. The matrix must be transposed. |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
113 |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
114 target[0] = static_cast<float>(source(0, 0)); |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
115 target[1] = static_cast<float>(source(1, 0)); |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
116 target[2] = 0; |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
117 target[3] = static_cast<float>(source(2, 0)); |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
118 target[4] = static_cast<float>(source(0, 1)); |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
119 target[5] = static_cast<float>(source(1, 1)); |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
120 target[6] = 0; |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
121 target[7] = static_cast<float>(source(2, 1)); |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
122 target[8] = 0; |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
123 target[9] = 0; |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
124 target[10] = -1; |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
125 target[11] = 0; |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
126 target[12] = static_cast<float>(source(0, 2)); |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
127 target[13] = static_cast<float>(source(1, 2)); |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
128 target[14] = -z; |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
129 target[15] = static_cast<float>(source(2, 2)); |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
130 } |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
131 |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
132 |
409 | 133 AffineTransform2D AffineTransform2D::Invert(const AffineTransform2D& a) |
134 { | |
135 AffineTransform2D t; | |
136 LinearAlgebra::InvertMatrix(t.matrix_, a.matrix_); | |
137 return t; | |
138 } | |
139 | |
140 | |
141 AffineTransform2D AffineTransform2D::Combine(const AffineTransform2D& a, | |
142 const AffineTransform2D& b) | |
143 { | |
144 return AffineTransform2D(LinearAlgebra::Product(a.GetHomogeneousMatrix(), | |
145 b.GetHomogeneousMatrix())); | |
146 } | |
147 | |
148 | |
149 AffineTransform2D AffineTransform2D::Combine(const AffineTransform2D& a, | |
150 const AffineTransform2D& b, | |
151 const AffineTransform2D& c) | |
152 { | |
153 return AffineTransform2D(LinearAlgebra::Product(a.GetHomogeneousMatrix(), | |
154 b.GetHomogeneousMatrix(), | |
155 c.GetHomogeneousMatrix())); | |
156 } | |
157 | |
158 | |
159 AffineTransform2D AffineTransform2D::Combine(const AffineTransform2D& a, | |
160 const AffineTransform2D& b, | |
161 const AffineTransform2D& c, | |
162 const AffineTransform2D& d) | |
163 { | |
164 return AffineTransform2D(LinearAlgebra::Product(a.GetHomogeneousMatrix(), | |
165 b.GetHomogeneousMatrix(), | |
166 c.GetHomogeneousMatrix(), | |
167 d.GetHomogeneousMatrix())); | |
168 } | |
169 | |
170 | |
171 AffineTransform2D AffineTransform2D::CreateOffset(double dx, | |
172 double dy) | |
173 { | |
174 AffineTransform2D t; | |
175 t.matrix_(0, 2) = dx; | |
176 t.matrix_(1, 2) = dy; | |
177 | |
178 return t; | |
179 } | |
180 | |
181 | |
182 AffineTransform2D AffineTransform2D::CreateScaling(double sx, | |
183 double sy) | |
184 { | |
185 AffineTransform2D t; | |
186 t.matrix_(0, 0) = sx; | |
187 t.matrix_(1, 1) = sy; | |
188 | |
189 return t; | |
190 } | |
191 | |
192 | |
193 AffineTransform2D AffineTransform2D::CreateRotation(double angle) | |
194 { | |
195 double cosine = cos(angle); | |
196 double sine = sin(angle); | |
197 | |
198 AffineTransform2D t; | |
199 t.matrix_(0, 0) = cosine; | |
200 t.matrix_(0, 1) = -sine; | |
201 t.matrix_(1, 0) = sine; | |
202 t.matrix_(1, 1) = cosine; | |
203 | |
204 return t; | |
205 } | |
574
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
206 |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
207 |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
208 AffineTransform2D AffineTransform2D::CreateOpenGLClipspace(unsigned int canvasWidth, |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
209 unsigned int canvasHeight) |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
210 { |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
211 AffineTransform2D t; |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
212 t.matrix_(0, 0) = 2.0 / static_cast<double>(canvasWidth); |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
213 t.matrix_(0, 2) = -1.0; |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
214 t.matrix_(1, 1) = -2.0 / static_cast<double>(canvasHeight); |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
215 t.matrix_(1, 2) = 1.0; |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
216 |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
217 return t; |
911297a277c4
AffineTransform2D::ConvertToOpenGLMatrix()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
218 } |
409 | 219 } |