Mercurial > hg > orthanc-stone
annotate Framework/Toolbox/GeometryToolbox.h @ 140:2115530d3703 wasm
OrientedBoundingBox
| author | Sebastien Jodogne <s.jodogne@gmail.com> |
|---|---|
| date | Thu, 18 Jan 2018 14:42:33 +0100 |
| parents | e2fe9352f240 |
| children | 9b83f30fc1c0 |
| 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 | |
| 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 |
|
140
2115530d3703
OrientedBoundingBox
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
135
diff
changeset
|
31 #include <boost/numeric/ublas/matrix.hpp> |
| 0 | 32 #include <boost/numeric/ublas/vector.hpp> |
| 33 | |
|
118
a4d0b6c82b29
using Orthanc::DicomMap instead of OrthancPlugins::DicomDatasetReader
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
113
diff
changeset
|
34 #include <Core/DicomFormat/DicomMap.h> |
| 0 | 35 |
| 36 namespace OrthancStone | |
| 37 { | |
|
140
2115530d3703
OrientedBoundingBox
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
135
diff
changeset
|
38 typedef boost::numeric::ublas::matrix<double> Matrix; |
| 0 | 39 typedef boost::numeric::ublas::vector<double> Vector; |
| 40 | |
| 41 namespace GeometryToolbox | |
| 42 { | |
| 43 void Print(const Vector& v); | |
| 44 | |
| 45 bool ParseVector(Vector& target, | |
| 32 | 46 const std::string& s); |
| 47 | |
| 48 bool ParseVector(Vector& target, | |
|
118
a4d0b6c82b29
using Orthanc::DicomMap instead of OrthancPlugins::DicomDatasetReader
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
113
diff
changeset
|
49 const Orthanc::DicomMap& dataset, |
|
a4d0b6c82b29
using Orthanc::DicomMap instead of OrthancPlugins::DicomDatasetReader
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
113
diff
changeset
|
50 const Orthanc::DicomTag& tag); |
| 0 | 51 |
| 52 void AssignVector(Vector& v, | |
| 53 double v1, | |
| 54 double v2); | |
| 55 | |
| 56 void AssignVector(Vector& v, | |
| 57 double v1, | |
| 58 double v2, | |
| 59 double v3); | |
| 60 | |
| 61 inline bool IsNear(double x, | |
| 62 double y, | |
| 63 double threshold) | |
| 64 { | |
| 65 return fabs(x - y) < threshold; | |
| 66 } | |
| 67 | |
| 68 bool IsNear(double x, | |
| 69 double y); | |
| 70 | |
| 71 inline bool IsCloseToZero(double x) | |
| 72 { | |
| 73 return IsNear(x, 0.0); | |
| 74 } | |
| 75 | |
| 76 void NormalizeVector(Vector& u); | |
| 77 | |
| 78 void CrossProduct(Vector& result, | |
| 79 const Vector& u, | |
| 80 const Vector& v); | |
| 81 | |
| 82 void ProjectPointOntoPlane(Vector& result, | |
| 83 const Vector& point, | |
| 84 const Vector& planeNormal, | |
| 85 const Vector& planeOrigin); | |
| 86 | |
| 87 bool IsParallel(const Vector& u, | |
| 88 const Vector& v); | |
| 89 | |
| 90 bool IsParallelOrOpposite(bool& isOpposite, | |
| 91 const Vector& u, | |
| 92 const Vector& v); | |
| 93 | |
| 94 bool IntersectTwoPlanes(Vector& p, | |
| 95 Vector& direction, | |
| 96 const Vector& origin1, | |
| 97 const Vector& normal1, | |
| 98 const Vector& origin2, | |
| 99 const Vector& normal2); | |
| 100 | |
| 101 bool ClipLineToRectangle(double& x1, // Coordinates of the clipped line (out) | |
| 102 double& y1, | |
| 103 double& x2, | |
| 104 double& y2, | |
| 105 const double ax, // Two points defining the line (in) | |
| 106 const double ay, | |
| 107 const double bx, | |
| 108 const double by, | |
| 109 const double& xmin, // Coordinates of the rectangle (in) | |
| 110 const double& ymin, | |
| 111 const double& xmax, | |
| 112 const double& ymax); | |
| 32 | 113 |
| 114 void GetPixelSpacing(double& spacingX, | |
| 115 double& spacingY, | |
|
118
a4d0b6c82b29
using Orthanc::DicomMap instead of OrthancPlugins::DicomDatasetReader
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
113
diff
changeset
|
116 const Orthanc::DicomMap& dicom); |
|
122
e3433dabfb8d
refactoring DicomStructureSet
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
118
diff
changeset
|
117 |
|
e3433dabfb8d
refactoring DicomStructureSet
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
118
diff
changeset
|
118 inline double ProjectAlongNormal(const Vector& point, |
|
e3433dabfb8d
refactoring DicomStructureSet
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
118
diff
changeset
|
119 const Vector& normal) |
|
e3433dabfb8d
refactoring DicomStructureSet
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
118
diff
changeset
|
120 { |
|
e3433dabfb8d
refactoring DicomStructureSet
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
118
diff
changeset
|
121 return boost::numeric::ublas::inner_prod(point, normal); |
|
e3433dabfb8d
refactoring DicomStructureSet
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
118
diff
changeset
|
122 } |
|
140
2115530d3703
OrientedBoundingBox
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
135
diff
changeset
|
123 |
|
2115530d3703
OrientedBoundingBox
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
135
diff
changeset
|
124 Matrix CreateRotationMatrixAlongX(double a); |
|
2115530d3703
OrientedBoundingBox
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
135
diff
changeset
|
125 |
|
2115530d3703
OrientedBoundingBox
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
135
diff
changeset
|
126 Matrix CreateRotationMatrixAlongY(double a); |
|
2115530d3703
OrientedBoundingBox
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
135
diff
changeset
|
127 |
|
2115530d3703
OrientedBoundingBox
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
135
diff
changeset
|
128 Matrix CreateRotationMatrixAlongZ(double a); |
| 0 | 129 }; |
| 130 } |