Mercurial > hg > orthanc-wsi
annotate Framework/Inputs/DicomPyramid.cpp @ 56:83cd735c885d
speedup the loading of DICOM sources
author | Sebastien Jodogne <s.jodogne@gmail.com> |
---|---|
date | Thu, 24 Nov 2016 15:41:21 +0100 |
parents | 7a88c614be04 |
children | 91fc9583b2de |
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 | |
56
83cd735c885d
speedup the loading of DICOM sources
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
111 if (a.GetPixelFormat() != b.GetPixelFormat() || |
0 | 112 a.GetTileWidth() != b.GetTileWidth() || |
113 a.GetTileHeight() != b.GetTileHeight() || | |
114 a.GetTotalWidth() < b.GetTotalWidth() || | |
115 a.GetTotalHeight() < b.GetTotalHeight()) | |
116 { | |
117 throw Orthanc::OrthancException(Orthanc::ErrorCode_IncompatibleImageFormat); | |
118 } | |
119 | |
120 if (a.GetTotalWidth() == b.GetTotalWidth() && | |
121 a.GetTotalHeight() != b.GetTotalHeight()) | |
122 { | |
123 throw Orthanc::OrthancException(Orthanc::ErrorCode_IncompatibleImageFormat); | |
124 } | |
125 } | |
126 } | |
127 | |
128 | |
129 void DicomPyramid::CheckLevel(size_t level) const | |
130 { | |
131 if (level >= levels_.size()) | |
132 { | |
133 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange); | |
134 } | |
135 } | |
136 | |
137 | |
138 DicomPyramid::DicomPyramid(IOrthancConnection& orthanc, | |
139 const std::string& seriesId) : | |
140 orthanc_(orthanc), | |
141 seriesId_(seriesId) | |
142 { | |
143 RegisterInstances(seriesId); | |
144 | |
145 // Sort the instances of the pyramid by decreasing total widths | |
146 std::sort(instances_.begin(), instances_.end(), Comparator()); | |
147 | |
148 try | |
149 { | |
150 Check(seriesId); | |
151 } | |
152 catch (Orthanc::OrthancException&) | |
153 { | |
154 Clear(); | |
155 throw; | |
156 } | |
157 | |
158 for (size_t i = 0; i < instances_.size(); i++) | |
159 { | |
160 if (i == 0 || | |
161 instances_[i - 1]->GetTotalWidth() != instances_[i]->GetTotalWidth()) | |
162 { | |
163 levels_.push_back(new DicomPyramidLevel(*instances_[i])); | |
164 } | |
165 else | |
166 { | |
167 assert(levels_.back() != NULL); | |
168 levels_.back()->AddInstance(*instances_[i]); | |
169 } | |
170 } | |
171 } | |
172 | |
173 | |
174 unsigned int DicomPyramid::GetLevelWidth(unsigned int level) const | |
175 { | |
176 CheckLevel(level); | |
177 return levels_[level]->GetTotalWidth(); | |
178 } | |
179 | |
180 | |
181 unsigned int DicomPyramid::GetLevelHeight(unsigned int level) const | |
182 { | |
183 CheckLevel(level); | |
184 return levels_[level]->GetTotalHeight(); | |
185 } | |
186 | |
187 | |
188 unsigned int DicomPyramid::GetTileWidth() const | |
189 { | |
190 assert(!levels_.empty() && levels_[0] != NULL); | |
191 return levels_[0]->GetTileWidth(); | |
192 } | |
193 | |
194 | |
195 unsigned int DicomPyramid::GetTileHeight() const | |
196 { | |
197 assert(!levels_.empty() && levels_[0] != NULL); | |
198 return levels_[0]->GetTileHeight(); | |
199 } | |
200 | |
201 | |
202 bool DicomPyramid::ReadRawTile(std::string& tile, | |
203 unsigned int level, | |
204 unsigned int tileX, | |
205 unsigned int tileY) | |
206 { | |
207 CheckLevel(level); | |
208 | |
209 ImageCompression compression; | |
210 Orthanc::PixelFormat format; | |
211 | |
212 if (levels_[level]->DownloadRawTile(compression, format, tile, orthanc_, tileX, tileY)) | |
213 { | |
214 assert(compression == GetImageCompression() && | |
215 format == GetPixelFormat()); | |
216 return true; | |
217 } | |
218 else | |
219 { | |
220 return false; | |
221 } | |
222 } | |
223 | |
224 | |
225 ImageCompression DicomPyramid::GetImageCompression() const | |
226 { | |
227 assert(!instances_.empty() && instances_[0] != NULL); | |
56
83cd735c885d
speedup the loading of DICOM sources
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
228 return instances_[0]->GetImageCompression(orthanc_); |
0 | 229 } |
230 | |
231 | |
232 Orthanc::PixelFormat DicomPyramid::GetPixelFormat() const | |
233 { | |
234 assert(!instances_.empty() && instances_[0] != NULL); | |
235 return instances_[0]->GetPixelFormat(); | |
236 } | |
237 } |