Mercurial > hg > orthanc-wsi
annotate Framework/Inputs/DicomPyramid.cpp @ 57:91fc9583b2de
big refactoring to support sparse tiling
| author | Sebastien Jodogne <s.jodogne@gmail.com> |
|---|---|
| date | Thu, 24 Nov 2016 17:48:24 +0100 |
| parents | 83cd735c885d |
| children | 7a3853d51c45 |
| 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, | |
|
57
91fc9583b2de
big refactoring to support sparse tiling
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
56
diff
changeset
|
203 ImageCompression& compression, |
| 0 | 204 unsigned int level, |
| 205 unsigned int tileX, | |
| 206 unsigned int tileY) | |
| 207 { | |
| 208 CheckLevel(level); | |
| 209 | |
| 210 Orthanc::PixelFormat format; | |
| 211 | |
|
57
91fc9583b2de
big refactoring to support sparse tiling
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
56
diff
changeset
|
212 if (levels_[level]->DownloadRawTile(tile, format, compression, orthanc_, tileX, tileY)) |
| 0 | 213 { |
|
57
91fc9583b2de
big refactoring to support sparse tiling
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
56
diff
changeset
|
214 if (format != GetPixelFormat()) |
|
91fc9583b2de
big refactoring to support sparse tiling
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
56
diff
changeset
|
215 { |
|
91fc9583b2de
big refactoring to support sparse tiling
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
56
diff
changeset
|
216 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); |
|
91fc9583b2de
big refactoring to support sparse tiling
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
56
diff
changeset
|
217 } |
|
91fc9583b2de
big refactoring to support sparse tiling
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
56
diff
changeset
|
218 |
| 0 | 219 return true; |
| 220 } | |
| 221 else | |
| 222 { | |
| 223 return false; | |
| 224 } | |
| 225 } | |
| 226 | |
| 227 | |
| 228 Orthanc::PixelFormat DicomPyramid::GetPixelFormat() const | |
| 229 { | |
| 230 assert(!instances_.empty() && instances_[0] != NULL); | |
| 231 return instances_[0]->GetPixelFormat(); | |
| 232 } | |
| 233 } |