Mercurial > hg > orthanc-stone
annotate Framework/Toolbox/GeometryToolbox.h @ 122:e3433dabfb8d wasm
refactoring DicomStructureSet
author | Sebastien Jodogne <s.jodogne@gmail.com> |
---|---|
date | Fri, 06 Oct 2017 17:25:08 +0200 |
parents | a4d0b6c82b29 |
children | e2fe9352f240 |
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 | |
40
7207a407bcd8
shared copyright with osimis
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
32
diff
changeset
|
5 * Copyright (C) 2017 Osimis, 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 | |
66 | 24 // Patch for ublas in Boost 1.64.0 |
65
885932a893de
OrthancFrameLayerSource
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
47
diff
changeset
|
25 // https://github.com/dealii/dealii/issues/4302 |
66 | 26 #include <boost/version.hpp> |
65
885932a893de
OrthancFrameLayerSource
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
47
diff
changeset
|
27 #if BOOST_VERSION >= 106300 // or 64, need to check |
885932a893de
OrthancFrameLayerSource
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
47
diff
changeset
|
28 # include <boost/serialization/array_wrapper.hpp> |
885932a893de
OrthancFrameLayerSource
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
47
diff
changeset
|
29 #endif |
885932a893de
OrthancFrameLayerSource
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
47
diff
changeset
|
30 |
0 | 31 #include <boost/numeric/ublas/vector.hpp> |
32 | |
118
a4d0b6c82b29
using Orthanc::DicomMap instead of OrthancPlugins::DicomDatasetReader
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
113
diff
changeset
|
33 #include <Core/DicomFormat/DicomMap.h> |
0 | 34 |
35 namespace OrthancStone | |
36 { | |
37 typedef boost::numeric::ublas::vector<double> Vector; | |
38 | |
39 namespace GeometryToolbox | |
40 { | |
41 void Print(const Vector& v); | |
42 | |
43 bool ParseVector(Vector& target, | |
32 | 44 const std::string& s); |
45 | |
46 bool ParseVector(Vector& target, | |
118
a4d0b6c82b29
using Orthanc::DicomMap instead of OrthancPlugins::DicomDatasetReader
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
113
diff
changeset
|
47 const Orthanc::DicomMap& dataset, |
a4d0b6c82b29
using Orthanc::DicomMap instead of OrthancPlugins::DicomDatasetReader
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
113
diff
changeset
|
48 const Orthanc::DicomTag& tag); |
0 | 49 |
50 void AssignVector(Vector& v, | |
51 double v1, | |
52 double v2); | |
53 | |
54 void AssignVector(Vector& v, | |
55 double v1, | |
56 double v2, | |
57 double v3); | |
58 | |
59 inline bool IsNear(double x, | |
60 double y, | |
61 double threshold) | |
62 { | |
63 return fabs(x - y) < threshold; | |
64 } | |
65 | |
66 bool IsNear(double x, | |
67 double y); | |
68 | |
69 inline bool IsCloseToZero(double x) | |
70 { | |
71 return IsNear(x, 0.0); | |
72 } | |
73 | |
74 void NormalizeVector(Vector& u); | |
75 | |
76 void CrossProduct(Vector& result, | |
77 const Vector& u, | |
78 const Vector& v); | |
79 | |
80 void ProjectPointOntoPlane(Vector& result, | |
81 const Vector& point, | |
82 const Vector& planeNormal, | |
83 const Vector& planeOrigin); | |
84 | |
85 bool IsParallel(const Vector& u, | |
86 const Vector& v); | |
87 | |
88 bool IsParallelOrOpposite(bool& isOpposite, | |
89 const Vector& u, | |
90 const Vector& v); | |
91 | |
92 bool IntersectTwoPlanes(Vector& p, | |
93 Vector& direction, | |
94 const Vector& origin1, | |
95 const Vector& normal1, | |
96 const Vector& origin2, | |
97 const Vector& normal2); | |
98 | |
99 bool ClipLineToRectangle(double& x1, // Coordinates of the clipped line (out) | |
100 double& y1, | |
101 double& x2, | |
102 double& y2, | |
103 const double ax, // Two points defining the line (in) | |
104 const double ay, | |
105 const double bx, | |
106 const double by, | |
107 const double& xmin, // Coordinates of the rectangle (in) | |
108 const double& ymin, | |
109 const double& xmax, | |
110 const double& ymax); | |
32 | 111 |
112 void GetPixelSpacing(double& spacingX, | |
113 double& spacingY, | |
118
a4d0b6c82b29
using Orthanc::DicomMap instead of OrthancPlugins::DicomDatasetReader
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
113
diff
changeset
|
114 const Orthanc::DicomMap& dicom); |
122
e3433dabfb8d
refactoring DicomStructureSet
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
118
diff
changeset
|
115 |
e3433dabfb8d
refactoring DicomStructureSet
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
118
diff
changeset
|
116 inline double ProjectAlongNormal(const Vector& point, |
e3433dabfb8d
refactoring DicomStructureSet
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
118
diff
changeset
|
117 const Vector& normal) |
e3433dabfb8d
refactoring DicomStructureSet
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
118
diff
changeset
|
118 { |
e3433dabfb8d
refactoring DicomStructureSet
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
118
diff
changeset
|
119 return boost::numeric::ublas::inner_prod(point, normal); |
e3433dabfb8d
refactoring DicomStructureSet
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
118
diff
changeset
|
120 } |
0 | 121 }; |
122 } |