Mercurial > hg > orthanc-stone
annotate Framework/Toolbox/ParallelSlices.cpp @ 119:ba83e38cf3ff wasm
rendering of rt-dose
author | Sebastien Jodogne <s.jodogne@gmail.com> |
---|---|
date | Mon, 02 Oct 2017 22:01:41 +0200 |
parents | 2eca030792aa |
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:
16
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 #include "ParallelSlices.h" | |
23 | |
113
2eca030792aa
using the Orthanc Framework
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
110
diff
changeset
|
24 #include <Core/Logging.h> |
2eca030792aa
using the Orthanc Framework
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
110
diff
changeset
|
25 #include <Core/OrthancException.h> |
0 | 26 |
27 namespace OrthancStone | |
28 { | |
29 ParallelSlices::ParallelSlices() | |
30 { | |
31 GeometryToolbox::AssignVector(normal_, 0, 0, 1); | |
32 } | |
33 | |
34 | |
35 ParallelSlices::ParallelSlices(const ParallelSlices& other) | |
36 { | |
37 normal_ = other.normal_; | |
38 | |
39 slices_.resize(other.slices_.size()); | |
40 | |
41 for (size_t i = 0; i < slices_.size(); i++) | |
42 { | |
43 assert(other.slices_[i] != NULL); | |
110
53025eecbc95
renamed SliceGeometry as CoordinateSystem3D
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
47
diff
changeset
|
44 slices_[i] = new CoordinateSystem3D(*other.slices_[i]); |
0 | 45 } |
46 } | |
47 | |
48 | |
49 ParallelSlices::~ParallelSlices() | |
50 { | |
51 for (size_t i = 0; i < slices_.size(); i++) | |
52 { | |
53 if (slices_[i] != NULL) | |
54 { | |
55 delete slices_[i]; | |
56 slices_[i] = NULL; | |
57 } | |
58 } | |
59 } | |
60 | |
61 | |
110
53025eecbc95
renamed SliceGeometry as CoordinateSystem3D
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
47
diff
changeset
|
62 void ParallelSlices::AddSlice(const CoordinateSystem3D& slice) |
0 | 63 { |
64 if (slices_.empty()) | |
65 { | |
66 normal_ = slice.GetNormal(); | |
110
53025eecbc95
renamed SliceGeometry as CoordinateSystem3D
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
47
diff
changeset
|
67 slices_.push_back(new CoordinateSystem3D(slice)); |
0 | 68 } |
69 else if (GeometryToolbox::IsParallel(slice.GetNormal(), normal_)) | |
70 { | |
110
53025eecbc95
renamed SliceGeometry as CoordinateSystem3D
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
47
diff
changeset
|
71 slices_.push_back(new CoordinateSystem3D(slice)); |
0 | 72 } |
73 else | |
74 { | |
75 LOG(ERROR) << "Trying to add a slice that is not parallel to the previous ones"; | |
76 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls); | |
77 } | |
78 } | |
79 | |
80 | |
81 void ParallelSlices::AddSlice(const Vector& origin, | |
82 const Vector& axisX, | |
83 const Vector& axisY) | |
84 { | |
110
53025eecbc95
renamed SliceGeometry as CoordinateSystem3D
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
47
diff
changeset
|
85 CoordinateSystem3D slice(origin, axisX, axisY); |
0 | 86 AddSlice(slice); |
87 } | |
88 | |
89 | |
110
53025eecbc95
renamed SliceGeometry as CoordinateSystem3D
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
47
diff
changeset
|
90 const CoordinateSystem3D& ParallelSlices::GetSlice(size_t index) const |
0 | 91 { |
92 if (index >= slices_.size()) | |
93 { | |
94 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange); | |
95 } | |
96 else | |
97 { | |
98 return *slices_[index]; | |
99 } | |
100 } | |
101 | |
102 | |
103 bool ParallelSlices::ComputeClosestSlice(size_t& closestSlice, | |
104 double& closestDistance, | |
105 const Vector& origin) const | |
106 { | |
107 if (slices_.empty()) | |
108 { | |
109 return false; | |
110 } | |
111 | |
112 double reference = boost::numeric::ublas::inner_prod(origin, normal_); | |
113 | |
114 closestSlice = 0; | |
115 closestDistance = std::numeric_limits<double>::infinity(); | |
116 | |
117 for (size_t i = 0; i < slices_.size(); i++) | |
118 { | |
119 double distance = fabs(boost::numeric::ublas::inner_prod(slices_[i]->GetOrigin(), normal_) - reference); | |
120 | |
121 if (distance < closestDistance) | |
122 { | |
123 closestSlice = i; | |
124 closestDistance = distance; | |
125 } | |
126 } | |
127 | |
128 return true; | |
129 } | |
130 | |
131 | |
132 ParallelSlices* ParallelSlices::Reverse() const | |
133 { | |
134 std::auto_ptr<ParallelSlices> reversed(new ParallelSlices); | |
135 | |
136 for (size_t i = slices_.size(); i > 0; i--) | |
137 { | |
110
53025eecbc95
renamed SliceGeometry as CoordinateSystem3D
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
47
diff
changeset
|
138 const CoordinateSystem3D& slice = *slices_[i - 1]; |
0 | 139 |
140 reversed->AddSlice(slice.GetOrigin(), | |
141 -slice.GetAxisX(), | |
142 slice.GetAxisY()); | |
143 } | |
144 | |
145 return reversed.release(); | |
146 } | |
147 } |