Mercurial > hg > orthanc-wsi
annotate Framework/Inputs/DicomPyramid.cpp @ 16:7a88c614be04
preparing for precompiled headers
| author | Sebastien Jodogne <s.jodogne@gmail.com> |
|---|---|
| date | Thu, 27 Oct 2016 11:01:26 +0200 |
| parents | 4a7a53257c7d |
| 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 } |