Mercurial > hg > orthanc-wsi
diff Framework/Inputs/DicomPyramid.cpp @ 0:4a7a53257c7d
initial commit
author | Sebastien Jodogne <s.jodogne@gmail.com> |
---|---|
date | Sat, 22 Oct 2016 21:48:33 +0200 |
parents | |
children | 7a88c614be04 |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/Framework/Inputs/DicomPyramid.cpp Sat Oct 22 21:48:33 2016 +0200 @@ -0,0 +1,237 @@ +/** + * Orthanc - A Lightweight, RESTful DICOM Store + * Copyright (C) 2012-2016 Sebastien Jodogne, Medical Physics + * Department, University Hospital of Liege, Belgium + * + * This program is free software: you can redistribute it and/or + * modify it under the terms of the GNU Affero General Public License + * as published by the Free Software Foundation, either version 3 of + * the License, or (at your option) any later version. + * + * This program is distributed in the hope that it will be useful, but + * WITHOUT ANY WARRANTY; without even the implied warranty of + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + * Affero General Public License for more details. + * + * You should have received a copy of the GNU Affero General Public License + * along with this program. If not, see <http://www.gnu.org/licenses/>. + **/ + + +#include "DicomPyramid.h" + +#include "../DicomToolbox.h" +#include "../Orthanc/Core/Logging.h" +#include "../Orthanc/Core/OrthancException.h" + +#include <algorithm> +#include <cassert> + +namespace OrthancWSI +{ + struct DicomPyramid::Comparator + { + bool operator() (DicomPyramidInstance* const& a, + DicomPyramidInstance* const& b) const + { + return a->GetTotalWidth() > b->GetTotalWidth(); + } + }; + + + void DicomPyramid::Clear() + { + for (size_t i = 0; i < levels_.size(); i++) + { + if (levels_[i] != NULL) + { + delete levels_[i]; + } + } + + for (size_t i = 0; i < instances_.size(); i++) + { + if (instances_[i] != NULL) + { + delete instances_[i]; + } + } + } + + + void DicomPyramid::RegisterInstances(const std::string& seriesId) + { + Json::Value series; + IOrthancConnection::RestApiGet(series, orthanc_, "/series/" + seriesId); + + if (series.type() != Json::objectValue) + { + throw Orthanc::OrthancException(Orthanc::ErrorCode_NetworkProtocol); + } + + const Json::Value& instances = DicomToolbox::GetSequenceTag(series, "Instances"); + instances_.reserve(instances.size()); + + for (Json::Value::ArrayIndex i = 0; i < instances.size(); i++) + { + if (instances[i].type() != Json::stringValue) + { + throw Orthanc::OrthancException(Orthanc::ErrorCode_NetworkProtocol); + } + + std::string instance = instances[i].asString(); + + try + { + instances_.push_back(new DicomPyramidInstance(orthanc_, instance)); + } + catch (Orthanc::OrthancException&) + { + LOG(ERROR) << "Skipping a DICOM instance that is not part of a whole-slide image: " << instance; + } + } + } + + + void DicomPyramid::Check(const std::string& seriesId) const + { + if (instances_.empty()) + { + LOG(ERROR) << "This series does not contain a whole-slide image: " << seriesId; + throw Orthanc::OrthancException(Orthanc::ErrorCode_UnknownResource); + } + + const DicomPyramidInstance& a = *instances_[0]; + + for (size_t i = 1; i < instances_.size(); i++) + { + const DicomPyramidInstance& b = *instances_[i]; + + if (a.GetImageCompression() != b.GetImageCompression() || + a.GetPixelFormat() != b.GetPixelFormat() || + a.GetTileWidth() != b.GetTileWidth() || + a.GetTileHeight() != b.GetTileHeight() || + a.GetTotalWidth() < b.GetTotalWidth() || + a.GetTotalHeight() < b.GetTotalHeight()) + { + throw Orthanc::OrthancException(Orthanc::ErrorCode_IncompatibleImageFormat); + } + + if (a.GetTotalWidth() == b.GetTotalWidth() && + a.GetTotalHeight() != b.GetTotalHeight()) + { + throw Orthanc::OrthancException(Orthanc::ErrorCode_IncompatibleImageFormat); + } + } + } + + + void DicomPyramid::CheckLevel(size_t level) const + { + if (level >= levels_.size()) + { + throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange); + } + } + + + DicomPyramid::DicomPyramid(IOrthancConnection& orthanc, + const std::string& seriesId) : + orthanc_(orthanc), + seriesId_(seriesId) + { + RegisterInstances(seriesId); + + // Sort the instances of the pyramid by decreasing total widths + std::sort(instances_.begin(), instances_.end(), Comparator()); + + try + { + Check(seriesId); + } + catch (Orthanc::OrthancException&) + { + Clear(); + throw; + } + + for (size_t i = 0; i < instances_.size(); i++) + { + if (i == 0 || + instances_[i - 1]->GetTotalWidth() != instances_[i]->GetTotalWidth()) + { + levels_.push_back(new DicomPyramidLevel(*instances_[i])); + } + else + { + assert(levels_.back() != NULL); + levels_.back()->AddInstance(*instances_[i]); + } + } + } + + + unsigned int DicomPyramid::GetLevelWidth(unsigned int level) const + { + CheckLevel(level); + return levels_[level]->GetTotalWidth(); + } + + + unsigned int DicomPyramid::GetLevelHeight(unsigned int level) const + { + CheckLevel(level); + return levels_[level]->GetTotalHeight(); + } + + + unsigned int DicomPyramid::GetTileWidth() const + { + assert(!levels_.empty() && levels_[0] != NULL); + return levels_[0]->GetTileWidth(); + } + + + unsigned int DicomPyramid::GetTileHeight() const + { + assert(!levels_.empty() && levels_[0] != NULL); + return levels_[0]->GetTileHeight(); + } + + + bool DicomPyramid::ReadRawTile(std::string& tile, + unsigned int level, + unsigned int tileX, + unsigned int tileY) + { + CheckLevel(level); + + ImageCompression compression; + Orthanc::PixelFormat format; + + if (levels_[level]->DownloadRawTile(compression, format, tile, orthanc_, tileX, tileY)) + { + assert(compression == GetImageCompression() && + format == GetPixelFormat()); + return true; + } + else + { + return false; + } + } + + + ImageCompression DicomPyramid::GetImageCompression() const + { + assert(!instances_.empty() && instances_[0] != NULL); + return instances_[0]->GetImageCompression(); + } + + + Orthanc::PixelFormat DicomPyramid::GetPixelFormat() const + { + assert(!instances_.empty() && instances_[0] != NULL); + return instances_[0]->GetPixelFormat(); + } +}