Mercurial > hg > orthanc-stone
annotate Framework/Toolbox/GeometryToolbox.h @ 177:83200c4d07ca wasm
fix interpolation
| author | Sebastien Jodogne <s.jodogne@gmail.com> |
|---|---|
| date | Fri, 09 Mar 2018 12:32:01 +0100 |
| parents | 15d92d93738b |
| children | 964118e7e6de |
| rev | line source |
|---|---|
| 0 | 1 /** |
| 2 * Stone of Orthanc | |
| 3 * Copyright (C) 2012-2016 Sebastien Jodogne, Medical Physics | |
| 4 * Department, University Hospital of Liege, Belgium | |
|
135
e2fe9352f240
upgrade to year 2018
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
5 * Copyright (C) 2017-2018 Osimis S.A., Belgium |
| 0 | 6 * |
| 7 * This program is free software: you can redistribute it and/or | |
| 47 | 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. | |
| 0 | 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 | |
| 47 | 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 | |
| 0 | 18 * along with this program. If not, see <http://www.gnu.org/licenses/>. |
| 19 **/ | |
| 20 | |
| 21 | |
| 22 #pragma once | |
| 23 | |
|
158
a053ca7fa5c6
LinearAlgebra toolbox
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
157
diff
changeset
|
24 #include "LinearAlgebra.h" |
| 0 | 25 |
| 26 namespace OrthancStone | |
| 27 { | |
| 28 namespace GeometryToolbox | |
| 29 { | |
| 30 void ProjectPointOntoPlane(Vector& result, | |
| 31 const Vector& point, | |
| 32 const Vector& planeNormal, | |
| 33 const Vector& planeOrigin); | |
| 34 | |
| 35 bool IsParallel(const Vector& u, | |
| 36 const Vector& v); | |
| 37 | |
| 38 bool IsParallelOrOpposite(bool& isOpposite, | |
| 39 const Vector& u, | |
| 40 const Vector& v); | |
| 41 | |
| 42 bool IntersectTwoPlanes(Vector& p, | |
| 43 Vector& direction, | |
| 44 const Vector& origin1, | |
| 45 const Vector& normal1, | |
| 46 const Vector& origin2, | |
| 47 const Vector& normal2); | |
| 48 | |
| 49 bool ClipLineToRectangle(double& x1, // Coordinates of the clipped line (out) | |
| 50 double& y1, | |
| 51 double& x2, | |
| 52 double& y2, | |
| 53 const double ax, // Two points defining the line (in) | |
| 54 const double ay, | |
| 55 const double bx, | |
| 56 const double by, | |
| 57 const double& xmin, // Coordinates of the rectangle (in) | |
| 58 const double& ymin, | |
| 59 const double& xmax, | |
| 60 const double& ymax); | |
| 32 | 61 |
| 62 void GetPixelSpacing(double& spacingX, | |
| 63 double& spacingY, | |
|
118
a4d0b6c82b29
using Orthanc::DicomMap instead of OrthancPlugins::DicomDatasetReader
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
113
diff
changeset
|
64 const Orthanc::DicomMap& dicom); |
|
122
e3433dabfb8d
refactoring DicomStructureSet
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
118
diff
changeset
|
65 |
|
e3433dabfb8d
refactoring DicomStructureSet
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
118
diff
changeset
|
66 inline double ProjectAlongNormal(const Vector& point, |
|
e3433dabfb8d
refactoring DicomStructureSet
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
118
diff
changeset
|
67 const Vector& normal) |
|
e3433dabfb8d
refactoring DicomStructureSet
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
118
diff
changeset
|
68 { |
|
e3433dabfb8d
refactoring DicomStructureSet
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
118
diff
changeset
|
69 return boost::numeric::ublas::inner_prod(point, normal); |
|
e3433dabfb8d
refactoring DicomStructureSet
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
118
diff
changeset
|
70 } |
|
140
2115530d3703
OrientedBoundingBox
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
135
diff
changeset
|
71 |
|
2115530d3703
OrientedBoundingBox
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
135
diff
changeset
|
72 Matrix CreateRotationMatrixAlongX(double a); |
|
2115530d3703
OrientedBoundingBox
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
135
diff
changeset
|
73 |
|
2115530d3703
OrientedBoundingBox
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
135
diff
changeset
|
74 Matrix CreateRotationMatrixAlongY(double a); |
|
2115530d3703
OrientedBoundingBox
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
135
diff
changeset
|
75 |
|
144
9b83f30fc1c0
bilinear and trilinear interpolation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
140
diff
changeset
|
76 Matrix CreateRotationMatrixAlongZ(double a); |
|
9b83f30fc1c0
bilinear and trilinear interpolation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
140
diff
changeset
|
77 |
|
165
8d50e6be565d
LinearAlgebra toolbox
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
158
diff
changeset
|
78 Matrix CreateTranslationMatrix(double dx, |
|
8d50e6be565d
LinearAlgebra toolbox
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
158
diff
changeset
|
79 double dy, |
|
8d50e6be565d
LinearAlgebra toolbox
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
158
diff
changeset
|
80 double dz); |
|
8d50e6be565d
LinearAlgebra toolbox
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
158
diff
changeset
|
81 |
|
8d50e6be565d
LinearAlgebra toolbox
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
158
diff
changeset
|
82 Matrix CreateScalingMatrix(double sx, |
|
8d50e6be565d
LinearAlgebra toolbox
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
158
diff
changeset
|
83 double sy, |
|
8d50e6be565d
LinearAlgebra toolbox
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
158
diff
changeset
|
84 double sz); |
|
8d50e6be565d
LinearAlgebra toolbox
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
158
diff
changeset
|
85 |
|
151
c5044bbfc303
CoordinateSystem3D::IntersectSegment()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
150
diff
changeset
|
86 bool IntersectPlaneAndSegment(Vector& p, |
|
c5044bbfc303
CoordinateSystem3D::IntersectSegment()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
150
diff
changeset
|
87 const Vector& normal, |
|
c5044bbfc303
CoordinateSystem3D::IntersectSegment()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
150
diff
changeset
|
88 double d, |
|
c5044bbfc303
CoordinateSystem3D::IntersectSegment()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
150
diff
changeset
|
89 const Vector& edgeFrom, |
|
c5044bbfc303
CoordinateSystem3D::IntersectSegment()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
150
diff
changeset
|
90 const Vector& edgeTo); |
|
c5044bbfc303
CoordinateSystem3D::IntersectSegment()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
150
diff
changeset
|
91 |
|
157
2309e8d86efe
IntersectPlaneAndLine
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
156
diff
changeset
|
92 bool IntersectPlaneAndLine(Vector& p, |
|
2309e8d86efe
IntersectPlaneAndLine
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
156
diff
changeset
|
93 const Vector& normal, |
|
2309e8d86efe
IntersectPlaneAndLine
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
156
diff
changeset
|
94 double d, |
|
2309e8d86efe
IntersectPlaneAndLine
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
156
diff
changeset
|
95 const Vector& origin, |
|
2309e8d86efe
IntersectPlaneAndLine
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
156
diff
changeset
|
96 const Vector& direction); |
|
2309e8d86efe
IntersectPlaneAndLine
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
156
diff
changeset
|
97 |
| 177 | 98 inline float ComputeBilinearInterpolationUnitSquare(float x, |
| 99 float y, | |
| 100 float f00, // source(0, 0) | |
| 101 float f01, // source(1, 0) | |
| 102 float f10, // source(0, 1) | |
| 103 float f11); // source(1, 1) | |
|
144
9b83f30fc1c0
bilinear and trilinear interpolation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
140
diff
changeset
|
104 |
| 177 | 105 inline float ComputeTrilinearInterpolationUnitSquare(float x, |
| 106 float y, | |
| 107 float z, | |
| 108 float f000, // source(0, 0, 0) | |
| 109 float f001, // source(1, 0, 0) | |
| 110 float f010, // source(0, 1, 0) | |
| 111 float f011, // source(1, 1, 0) | |
| 112 float f100, // source(0, 0, 1) | |
| 113 float f101, // source(1, 0, 1) | |
| 114 float f110, // source(0, 1, 1) | |
| 115 float f111); // source(1, 1, 1) | |
| 0 | 116 }; |
| 117 } | |
|
144
9b83f30fc1c0
bilinear and trilinear interpolation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
140
diff
changeset
|
118 |
|
9b83f30fc1c0
bilinear and trilinear interpolation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
140
diff
changeset
|
119 |
| 177 | 120 float OrthancStone::GeometryToolbox::ComputeBilinearInterpolationUnitSquare(float x, |
| 121 float y, | |
| 122 float f00, | |
| 123 float f01, | |
| 124 float f10, | |
| 125 float f11) | |
|
144
9b83f30fc1c0
bilinear and trilinear interpolation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
140
diff
changeset
|
126 { |
| 177 | 127 // This function only works within the unit square |
| 128 assert(x >= 0 && y >= 0 && x <= 1 && y <= 1); | |
|
144
9b83f30fc1c0
bilinear and trilinear interpolation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
140
diff
changeset
|
129 |
|
9b83f30fc1c0
bilinear and trilinear interpolation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
140
diff
changeset
|
130 // https://en.wikipedia.org/wiki/Bilinear_interpolation#Unit_square |
|
175
15d92d93738b
fix interpolation with negative coordinates
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
172
diff
changeset
|
131 return (f00 * (1.0f - x) * (1.0f - y) + |
|
15d92d93738b
fix interpolation with negative coordinates
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
172
diff
changeset
|
132 f01 * x * (1.0f - y) + |
|
15d92d93738b
fix interpolation with negative coordinates
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
172
diff
changeset
|
133 f10 * (1.0f - x) * y + |
|
15d92d93738b
fix interpolation with negative coordinates
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
172
diff
changeset
|
134 f11 * x * y); |
|
144
9b83f30fc1c0
bilinear and trilinear interpolation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
140
diff
changeset
|
135 } |
|
9b83f30fc1c0
bilinear and trilinear interpolation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
140
diff
changeset
|
136 |
|
9b83f30fc1c0
bilinear and trilinear interpolation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
140
diff
changeset
|
137 |
| 177 | 138 float OrthancStone::GeometryToolbox::ComputeTrilinearInterpolationUnitSquare(float x, |
| 139 float y, | |
| 140 float z, | |
| 141 float f000, | |
| 142 float f001, | |
| 143 float f010, | |
| 144 float f011, | |
| 145 float f100, | |
| 146 float f101, | |
| 147 float f110, | |
| 148 float f111) | |
|
144
9b83f30fc1c0
bilinear and trilinear interpolation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
140
diff
changeset
|
149 { |
|
9b83f30fc1c0
bilinear and trilinear interpolation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
140
diff
changeset
|
150 // "In practice, a trilinear interpolation is identical to two |
|
9b83f30fc1c0
bilinear and trilinear interpolation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
140
diff
changeset
|
151 // bilinear interpolation combined with a linear interpolation" |
|
9b83f30fc1c0
bilinear and trilinear interpolation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
140
diff
changeset
|
152 // https://en.wikipedia.org/wiki/Trilinear_interpolation#Method |
| 177 | 153 float a = ComputeBilinearInterpolationUnitSquare(x, y, f000, f001, f010, f011); |
| 154 float b = ComputeBilinearInterpolationUnitSquare(x, y, f100, f101, f110, f111); | |
|
144
9b83f30fc1c0
bilinear and trilinear interpolation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
140
diff
changeset
|
155 |
| 177 | 156 return (1.0f - z) * a + z * b; |
|
144
9b83f30fc1c0
bilinear and trilinear interpolation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
140
diff
changeset
|
157 } |