Mercurial > hg > orthanc-stone
annotate Framework/Toolbox/SlicesSorter.cpp @ 700:059e1fd05fd6 refactor-viewport-controller
Introduced the ViewportController that sits between the application and the
Scene2D to handle the trackers and measuring tools. This is a work in progress.
The Scene2D is no longer an observable. Message sending is managed by the
ViewportController.
Move some refs to shared and weak to prevent lifetime issues.
| author | Benjamin Golinvaux <bgo@osimis.io> |
|---|---|
| date | Sun, 19 May 2019 16:26:17 +0200 |
| parents | e9339f2b5de7 |
| children | f7c236894c1a |
| rev | line source |
|---|---|
| 73 | 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 |
| 73 | 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 "SlicesSorter.h" | |
| 23 | |
|
159
0a73d76333db
populating LinearAlgebra
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
135
diff
changeset
|
24 #include "GeometryToolbox.h" |
|
0a73d76333db
populating LinearAlgebra
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
135
diff
changeset
|
25 |
|
212
5412adf19980
resort to OrthancFramework
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
201
diff
changeset
|
26 #include <Core/OrthancException.h> |
| 73 | 27 |
| 28 namespace OrthancStone | |
| 29 { | |
| 30 class SlicesSorter::SliceWithDepth : public boost::noncopyable | |
| 31 { | |
| 32 private: | |
|
647
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
33 CoordinateSystem3D geometry_; |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
34 double depth_; |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
35 |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
36 std::auto_ptr<Orthanc::IDynamicObject> payload_; |
| 73 | 37 |
| 38 public: | |
|
647
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
39 SliceWithDepth(const CoordinateSystem3D& geometry, |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
40 Orthanc::IDynamicObject* payload) : |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
41 geometry_(geometry), |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
42 depth_(0), |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
43 payload_(payload) |
| 73 | 44 { |
| 45 } | |
| 46 | |
| 47 void SetNormal(const Vector& normal) | |
| 48 { | |
|
647
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
49 depth_ = boost::numeric::ublas::inner_prod(geometry_.GetOrigin(), normal); |
| 73 | 50 } |
| 51 | |
| 52 double GetDepth() const | |
| 53 { | |
| 54 return depth_; | |
| 55 } | |
| 56 | |
|
647
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
57 const CoordinateSystem3D& GetGeometry() const |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
58 { |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
59 return geometry_; |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
60 } |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
61 |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
62 bool HasPayload() const |
| 73 | 63 { |
|
647
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
64 return (payload_.get() != NULL); |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
65 } |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
66 |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
67 const Orthanc::IDynamicObject& GetPayload() const |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
68 { |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
69 if (HasPayload()) |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
70 { |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
71 return *payload_; |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
72 } |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
73 else |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
74 { |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
75 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls); |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
76 } |
| 73 | 77 } |
| 78 }; | |
| 79 | |
| 80 | |
| 81 struct SlicesSorter::Comparator | |
| 82 { | |
| 83 bool operator() (const SliceWithDepth* const& a, | |
| 84 const SliceWithDepth* const& b) const | |
| 85 { | |
| 86 return a->GetDepth() < b->GetDepth(); | |
| 87 } | |
| 88 }; | |
| 89 | |
| 90 | |
| 91 SlicesSorter::~SlicesSorter() | |
| 92 { | |
| 93 for (size_t i = 0; i < slices_.size(); i++) | |
| 94 { | |
| 95 assert(slices_[i] != NULL); | |
| 96 delete slices_[i]; | |
| 97 } | |
| 98 } | |
| 99 | |
| 100 | |
|
647
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
101 void SlicesSorter::AddSlice(const CoordinateSystem3D& slice, |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
102 Orthanc::IDynamicObject* payload) |
| 73 | 103 { |
|
647
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
104 slices_.push_back(new SliceWithDepth(slice, payload)); |
| 73 | 105 } |
| 106 | |
| 107 | |
|
647
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
108 const SlicesSorter::SliceWithDepth& SlicesSorter::GetSlice(size_t i) const |
| 73 | 109 { |
| 110 if (i >= slices_.size()) | |
| 111 { | |
| 112 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange); | |
| 113 } | |
|
647
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
114 else |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
115 { |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
116 assert(slices_[i] != NULL); |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
117 return *slices_[i]; |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
118 } |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
119 } |
| 73 | 120 |
|
647
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
121 |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
122 const CoordinateSystem3D& SlicesSorter::GetSliceGeometry(size_t i) const |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
123 { |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
124 return GetSlice(i).GetGeometry(); |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
125 } |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
126 |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
127 |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
128 bool SlicesSorter::HasSlicePayload(size_t i) const |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
129 { |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
130 return GetSlice(i).HasPayload(); |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
131 } |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
132 |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
133 |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
134 const Orthanc::IDynamicObject& SlicesSorter::GetSlicePayload(size_t i) const |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
135 { |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
136 return GetSlice(i).GetPayload(); |
| 73 | 137 } |
| 138 | |
| 139 | |
| 140 void SlicesSorter::SetNormal(const Vector& normal) | |
| 141 { | |
| 142 for (size_t i = 0; i < slices_.size(); i++) | |
| 143 { | |
| 144 slices_[i]->SetNormal(normal); | |
| 145 } | |
| 146 | |
| 147 hasNormal_ = true; | |
| 148 } | |
| 149 | |
| 150 | |
| 648 | 151 void SlicesSorter::SortInternal() |
| 73 | 152 { |
| 153 if (!hasNormal_) | |
| 154 { | |
| 155 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls); | |
| 156 } | |
| 157 | |
| 158 Comparator comparator; | |
| 159 std::sort(slices_.begin(), slices_.end(), comparator); | |
| 160 } | |
| 161 | |
| 162 | |
| 163 void SlicesSorter::FilterNormal(const Vector& normal) | |
| 164 { | |
| 165 size_t pos = 0; | |
| 166 | |
| 167 for (size_t i = 0; i < slices_.size(); i++) | |
| 168 { | |
|
647
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
169 if (GeometryToolbox::IsParallel(normal, slices_[i]->GetGeometry().GetNormal())) |
| 73 | 170 { |
| 171 // This slice is compatible with the selected normal | |
| 172 slices_[pos] = slices_[i]; | |
| 173 pos += 1; | |
| 174 } | |
| 175 else | |
| 176 { | |
| 177 delete slices_[i]; | |
| 178 slices_[i] = NULL; | |
| 179 } | |
| 180 } | |
| 181 | |
| 182 slices_.resize(pos); | |
| 183 } | |
| 184 | |
| 185 | |
| 186 bool SlicesSorter::SelectNormal(Vector& normal) const | |
| 187 { | |
| 188 std::vector<Vector> normalCandidates; | |
| 189 std::vector<unsigned int> normalCount; | |
| 190 | |
| 191 bool found = false; | |
| 192 | |
|
647
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
193 for (size_t i = 0; !found && i < GetSlicesCount(); i++) |
| 73 | 194 { |
| 195 const Vector& normal = GetSlice(i).GetGeometry().GetNormal(); | |
| 196 | |
| 197 bool add = true; | |
| 198 for (size_t j = 0; add && j < normalCandidates.size(); j++) // (*) | |
| 199 { | |
| 200 if (GeometryToolbox::IsParallel(normal, normalCandidates[j])) | |
| 201 { | |
| 202 normalCount[j] += 1; | |
| 203 add = false; | |
| 204 } | |
| 205 } | |
| 206 | |
| 207 if (add) | |
| 208 { | |
| 209 if (normalCount.size() > 2) | |
| 210 { | |
| 211 // To get linear-time complexity in (*). This heuristics | |
| 212 // allows the series to have one single frame that is | |
| 213 // not parallel to the others (such a frame could be a | |
| 214 // generated preview) | |
| 215 found = false; | |
| 216 } | |
| 217 else | |
| 218 { | |
| 219 normalCandidates.push_back(normal); | |
| 220 normalCount.push_back(1); | |
| 221 } | |
| 222 } | |
| 223 } | |
| 224 | |
| 225 for (size_t i = 0; !found && i < normalCandidates.size(); i++) | |
| 226 { | |
| 227 unsigned int count = normalCount[i]; | |
|
647
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
228 if (count == GetSlicesCount() || |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
229 count + 1 == GetSlicesCount()) |
| 73 | 230 { |
| 231 normal = normalCandidates[i]; | |
| 232 found = true; | |
| 233 } | |
| 234 } | |
| 235 | |
| 236 return found; | |
| 237 } | |
| 77 | 238 |
| 239 | |
| 648 | 240 bool SlicesSorter::Sort() |
| 241 { | |
| 242 if (GetSlicesCount() > 0) | |
| 243 { | |
| 244 Vector normal; | |
| 245 if (SelectNormal(normal)) | |
| 246 { | |
| 247 FilterNormal(normal); | |
| 248 SetNormal(normal); | |
| 249 SortInternal(); | |
| 250 return true; | |
| 251 } | |
| 252 } | |
| 253 | |
| 254 return false; | |
| 255 } | |
| 256 | |
| 257 | |
|
647
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
258 bool SlicesSorter::LookupClosestSlice(size_t& index, |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
259 double& distance, |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
260 const CoordinateSystem3D& slice) const |
| 77 | 261 { |
| 262 // TODO Turn this linear-time lookup into a log-time lookup, | |
| 263 // keeping track of whether the slices are sorted along the normal | |
| 264 | |
|
647
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
265 bool found = false; |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
266 |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
267 distance = std::numeric_limits<double>::infinity(); |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
268 |
| 77 | 269 for (size_t i = 0; i < slices_.size(); i++) |
| 270 { | |
|
647
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
271 assert(slices_[i] != NULL); |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
272 |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
273 double tmp; |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
274 if (CoordinateSystem3D::GetDistance(tmp, slices_[i]->GetGeometry(), slice)) |
| 77 | 275 { |
|
647
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
276 if (!found || |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
277 tmp < distance) |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
278 { |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
279 index = i; |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
280 distance = tmp; |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
281 found = true; |
|
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
282 } |
| 77 | 283 } |
| 284 } | |
| 285 | |
|
647
6af3099ed8da
uncoupling OrthancStone::SlicesSorter from OrthancStone::Slice
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
439
diff
changeset
|
286 return found; |
| 77 | 287 } |
|
667
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
288 |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
289 |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
290 double SlicesSorter::ComputeSpacingBetweenSlices() const |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
291 { |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
292 if (GetSlicesCount() <= 1) |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
293 { |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
294 // This is a volume that is empty or that contains one single |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
295 // slice: Choose a dummy z-dimension for voxels |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
296 return 1.0; |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
297 } |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
298 |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
299 const OrthancStone::CoordinateSystem3D& reference = GetSliceGeometry(0); |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
300 |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
301 double referencePosition = reference.ProjectAlongNormal(reference.GetOrigin()); |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
302 |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
303 double p = reference.ProjectAlongNormal(GetSliceGeometry(1).GetOrigin()); |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
304 double spacingZ = p - referencePosition; |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
305 |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
306 if (spacingZ <= 0) |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
307 { |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
308 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls, |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
309 "Please call the Sort() method before"); |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
310 } |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
311 |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
312 for (size_t i = 1; i < GetSlicesCount(); i++) |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
313 { |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
314 OrthancStone::Vector p = reference.GetOrigin() + spacingZ * static_cast<double>(i) * reference.GetNormal(); |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
315 double d = boost::numeric::ublas::norm_2(p - GetSliceGeometry(i).GetOrigin()); |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
316 |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
317 if (!OrthancStone::LinearAlgebra::IsNear(d, 0, 0.001 /* tolerance expressed in mm */)) |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
318 { |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
319 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadGeometry, |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
320 "The origins of the slices of a volume image are not regularly spaced"); |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
321 } |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
322 } |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
323 |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
324 return spacingZ; |
|
e9339f2b5de7
refactoring of VolumeImage
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
648
diff
changeset
|
325 } |
| 73 | 326 } |