Mercurial > hg > orthanc-wsi
annotate Framework/Inputs/DicomPyramid.cpp @ 55:b6432a00b103
use of simpler data structures
author | Sebastien Jodogne <s.jodogne@gmail.com> |
---|---|
date | Thu, 24 Nov 2016 14:44:11 +0100 |
parents | 7a88c614be04 |
children | 83cd735c885d |
rev | line source |
---|---|
0 | 1 /** |
2 * Orthanc - A Lightweight, RESTful DICOM Store | |
3 * Copyright (C) 2012-2016 Sebastien Jodogne, Medical Physics | |
4 * Department, University Hospital of Liege, Belgium | |
5 * | |
6 * This program is free software: you can redistribute it and/or | |
7 * modify it under the terms of the GNU Affero General Public License | |
8 * as published by the Free Software Foundation, either version 3 of | |
9 * the License, or (at your option) any later version. | |
10 * | |
11 * This program is distributed in the hope that it will be useful, but | |
12 * WITHOUT ANY WARRANTY; without even the implied warranty of | |
13 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU | |
14 * Affero General Public License for more details. | |
15 * | |
16 * You should have received a copy of the GNU Affero General Public License | |
17 * along with this program. If not, see <http://www.gnu.org/licenses/>. | |
18 **/ | |
19 | |
20 | |
16
7a88c614be04
preparing for precompiled headers
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
21 #include "../PrecompiledHeadersWSI.h" |
0 | 22 #include "DicomPyramid.h" |
23 | |
24 #include "../DicomToolbox.h" | |
25 #include "../Orthanc/Core/Logging.h" | |
26 #include "../Orthanc/Core/OrthancException.h" | |
27 | |
28 #include <algorithm> | |
29 #include <cassert> | |
30 | |
31 namespace OrthancWSI | |
32 { | |
33 struct DicomPyramid::Comparator | |
34 { | |
35 bool operator() (DicomPyramidInstance* const& a, | |
36 DicomPyramidInstance* const& b) const | |
37 { | |
38 return a->GetTotalWidth() > b->GetTotalWidth(); | |
39 } | |
40 }; | |
41 | |
42 | |
43 void DicomPyramid::Clear() | |
44 { | |
45 for (size_t i = 0; i < levels_.size(); i++) | |
46 { | |
47 if (levels_[i] != NULL) | |
48 { | |
49 delete levels_[i]; | |
50 } | |
51 } | |
52 | |
53 for (size_t i = 0; i < instances_.size(); i++) | |
54 { | |
55 if (instances_[i] != NULL) | |
56 { | |
57 delete instances_[i]; | |
58 } | |
59 } | |
60 } | |
61 | |
62 | |
63 void DicomPyramid::RegisterInstances(const std::string& seriesId) | |
64 { | |
65 Json::Value series; | |
66 IOrthancConnection::RestApiGet(series, orthanc_, "/series/" + seriesId); | |
67 | |
68 if (series.type() != Json::objectValue) | |
69 { | |
70 throw Orthanc::OrthancException(Orthanc::ErrorCode_NetworkProtocol); | |
71 } | |
72 | |
73 const Json::Value& instances = DicomToolbox::GetSequenceTag(series, "Instances"); | |
74 instances_.reserve(instances.size()); | |
75 | |
76 for (Json::Value::ArrayIndex i = 0; i < instances.size(); i++) | |
77 { | |
78 if (instances[i].type() != Json::stringValue) | |
79 { | |
80 throw Orthanc::OrthancException(Orthanc::ErrorCode_NetworkProtocol); | |
81 } | |
82 | |
83 std::string instance = instances[i].asString(); | |
84 | |
85 try | |
86 { | |
87 instances_.push_back(new DicomPyramidInstance(orthanc_, instance)); | |
88 } | |
89 catch (Orthanc::OrthancException&) | |
90 { | |
91 LOG(ERROR) << "Skipping a DICOM instance that is not part of a whole-slide image: " << instance; | |
92 } | |
93 } | |
94 } | |
95 | |
96 | |
97 void DicomPyramid::Check(const std::string& seriesId) const | |
98 { | |
99 if (instances_.empty()) | |
100 { | |
101 LOG(ERROR) << "This series does not contain a whole-slide image: " << seriesId; | |
102 throw Orthanc::OrthancException(Orthanc::ErrorCode_UnknownResource); | |
103 } | |
104 | |
105 const DicomPyramidInstance& a = *instances_[0]; | |
106 | |
107 for (size_t i = 1; i < instances_.size(); i++) | |
108 { | |
109 const DicomPyramidInstance& b = *instances_[i]; | |
110 | |
111 if (a.GetImageCompression() != b.GetImageCompression() || | |
112 a.GetPixelFormat() != b.GetPixelFormat() || | |
113 a.GetTileWidth() != b.GetTileWidth() || | |
114 a.GetTileHeight() != b.GetTileHeight() || | |
115 a.GetTotalWidth() < b.GetTotalWidth() || | |
116 a.GetTotalHeight() < b.GetTotalHeight()) | |
117 { | |
118 throw Orthanc::OrthancException(Orthanc::ErrorCode_IncompatibleImageFormat); | |
119 } | |
120 | |
121 if (a.GetTotalWidth() == b.GetTotalWidth() && | |
122 a.GetTotalHeight() != b.GetTotalHeight()) | |
123 { | |
124 throw Orthanc::OrthancException(Orthanc::ErrorCode_IncompatibleImageFormat); | |
125 } | |
126 } | |
127 } | |
128 | |
129 | |
130 void DicomPyramid::CheckLevel(size_t level) const | |
131 { | |
132 if (level >= levels_.size()) | |
133 { | |
134 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange); | |
135 } | |
136 } | |
137 | |
138 | |
139 DicomPyramid::DicomPyramid(IOrthancConnection& orthanc, | |
140 const std::string& seriesId) : | |
141 orthanc_(orthanc), | |
142 seriesId_(seriesId) | |
143 { | |
144 RegisterInstances(seriesId); | |
145 | |
146 // Sort the instances of the pyramid by decreasing total widths | |
147 std::sort(instances_.begin(), instances_.end(), Comparator()); | |
148 | |
149 try | |
150 { | |
151 Check(seriesId); | |
152 } | |
153 catch (Orthanc::OrthancException&) | |
154 { | |
155 Clear(); | |
156 throw; | |
157 } | |
158 | |
159 for (size_t i = 0; i < instances_.size(); i++) | |
160 { | |
161 if (i == 0 || | |
162 instances_[i - 1]->GetTotalWidth() != instances_[i]->GetTotalWidth()) | |
163 { | |
164 levels_.push_back(new DicomPyramidLevel(*instances_[i])); | |
165 } | |
166 else | |
167 { | |
168 assert(levels_.back() != NULL); | |
169 levels_.back()->AddInstance(*instances_[i]); | |
170 } | |
171 } | |
172 } | |
173 | |
174 | |
175 unsigned int DicomPyramid::GetLevelWidth(unsigned int level) const | |
176 { | |
177 CheckLevel(level); | |
178 return levels_[level]->GetTotalWidth(); | |
179 } | |
180 | |
181 | |
182 unsigned int DicomPyramid::GetLevelHeight(unsigned int level) const | |
183 { | |
184 CheckLevel(level); | |
185 return levels_[level]->GetTotalHeight(); | |
186 } | |
187 | |
188 | |
189 unsigned int DicomPyramid::GetTileWidth() const | |
190 { | |
191 assert(!levels_.empty() && levels_[0] != NULL); | |
192 return levels_[0]->GetTileWidth(); | |
193 } | |
194 | |
195 | |
196 unsigned int DicomPyramid::GetTileHeight() const | |
197 { | |
198 assert(!levels_.empty() && levels_[0] != NULL); | |
199 return levels_[0]->GetTileHeight(); | |
200 } | |
201 | |
202 | |
203 bool DicomPyramid::ReadRawTile(std::string& tile, | |
204 unsigned int level, | |
205 unsigned int tileX, | |
206 unsigned int tileY) | |
207 { | |
208 CheckLevel(level); | |
209 | |
210 ImageCompression compression; | |
211 Orthanc::PixelFormat format; | |
212 | |
213 if (levels_[level]->DownloadRawTile(compression, format, tile, orthanc_, tileX, tileY)) | |
214 { | |
215 assert(compression == GetImageCompression() && | |
216 format == GetPixelFormat()); | |
217 return true; | |
218 } | |
219 else | |
220 { | |
221 return false; | |
222 } | |
223 } | |
224 | |
225 | |
226 ImageCompression DicomPyramid::GetImageCompression() const | |
227 { | |
228 assert(!instances_.empty() && instances_[0] != NULL); | |
229 return instances_[0]->GetImageCompression(); | |
230 } | |
231 | |
232 | |
233 Orthanc::PixelFormat DicomPyramid::GetPixelFormat() const | |
234 { | |
235 assert(!instances_.empty() && instances_[0] != NULL); | |
236 return instances_[0]->GetPixelFormat(); | |
237 } | |
238 } |