Mercurial > hg > orthanc-neuro
view Sources/Framework/IDicomFrameDecoder.cpp @ 12:8747857b539e
Support Philips private tags for rescale slope/intercept
author | Sebastien Jodogne <s.jodogne@gmail.com> |
---|---|
date | Wed, 04 May 2022 17:52:56 +0200 |
parents | a2d79b456440 |
children | 2c8f8168c1ab |
line wrap: on
line source
/** * Neuroimaging plugin for Orthanc * Copyright (C) 2021-2022 Sebastien Jodogne, UCLouvain, Belgium * * This program is free software: you can redistribute it and/or * modify it under the terms of the GNU 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 * General Public License for more details. * * You should have received a copy of the GNU General Public License * along with this program. If not, see <http://www.gnu.org/licenses/>. **/ #include "IDicomFrameDecoder.h" #include <OrthancException.h> namespace Neuro { void IDicomFrameDecoder::Apply(NiftiWriter& writer, IDicomFrameDecoder& decoder, const std::vector<Slice>& slices) { for (size_t i = 1; i < slices.size(); i++) { if (slices[0].GetWidth() != slices[i].GetWidth() || slices[0].GetHeight() != slices[i].GetHeight()) { throw Orthanc::OrthancException(Orthanc::ErrorCode_NotImplemented, "The slices have varying dimensions"); } } size_t currentInstanceIndex; unsigned int currentFrameNumber; std::unique_ptr<IDecodedFrame> currentFrame; bool first = true; Orthanc::PixelFormat format; for (size_t i = 0; i < slices.size(); i++) { if (currentFrame.get() == NULL || currentInstanceIndex != slices[i].GetInstanceIndexInCollection() || currentFrameNumber != slices[i].GetFrameNumber()) { currentFrame.reset(decoder.DecodeFrame(slices[i])); if (currentFrame.get() == NULL) { throw Orthanc::OrthancException(Orthanc::ErrorCode_NullPointer); } currentInstanceIndex = slices[i].GetInstanceIndexInCollection(); currentFrameNumber = slices[i].GetFrameNumber(); } Orthanc::ImageAccessor region; currentFrame->GetRegion(region, slices[i].GetX(), slices[i].GetY(), slices[i].GetWidth(), slices[i].GetHeight()); if (region.GetWidth() != slices[i].GetWidth() || region.GetHeight() != slices[i].GetHeight()) { throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError); } if (first) { first = false; format = region.GetFormat(); } if (region.GetFormat() != format) { throw Orthanc::OrthancException(Orthanc::ErrorCode_IncompatibleImageFormat, "The slices have varying pixel formats"); } writer.AddSlice(region); } } }