orthanc-wsi: Framework/Inputs/DicomPyramid.cpp comparison

comparison 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

comparison

equal deleted inserted replaced

--1:000000000000
+:4a7a53257c7d
+/**
+* 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();
+}
+}

Mercurial > hg > orthanc-wsi

comparison Framework/Inputs/DicomPyramid.cpp @ 0:4a7a53257c7d