Mercurial > hg > orthanc-neuro
view Sources/Framework/CSAHeader.cpp @ 11:12e36b1b84f3
testing with dcmtk
author | Sebastien Jodogne <s.jodogne@gmail.com> |
---|---|
date | Sat, 30 Apr 2022 11:26:21 +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 "CSAHeader.h" #include "BufferReader.h" #include <OrthancException.h> #include <cassert> namespace Neuro { void CSAHeader::Clear() { for (Content::iterator it = content_.begin(); it != content_.end(); ++it) { assert(it->second != NULL); delete it->second; } } void CSAHeader::Load(const std::string& tag) { // https://nipy.org/nibabel/dicom/siemens_csa.html Clear(); BufferReader reader(tag); if (reader.ReadUInt32() != 0x30315653) // This is the "SV10" header { throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); } reader.ReadUInt32(); // Unused, often equals to 0x01020304 const uint32_t n_tags = reader.ReadUInt32(); if (n_tags == 0 || n_tags > 128) { // This should in the range 1..128 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); } if (reader.ReadUInt32() != 77) { throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); } for (uint32_t i = 0; i < n_tags; i++) { const std::string name = reader.ReadNullTerminatedString(); if (name.size() >= 63) { throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); } reader.Skip(64 - name.size() - 1); const uint32_t vm = reader.ReadUInt32(); const std::string vr = reader.ReadNullTerminatedString(); if (vr.size() >= 4) { throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); } reader.Skip(4 - vr.size() - 1); reader.ReadUInt32(); // "syngodt" = syngo.via data type const uint32_t nitems = reader.ReadUInt32(); const uint32_t sync = reader.ReadUInt32(); if (sync != 77 && sync != 205) { throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); } std::unique_ptr<CSATag> tag2(new CSATag(vr)); for (uint32_t j = 0; j < nitems; j++) { reader.ReadUInt32(); const uint32_t item_len = reader.ReadUInt32(); reader.ReadUInt32(); reader.ReadUInt32(); if (vm == 0 || j < vm) { tag2->AddValue(reader.ReadBlock(item_len)); } else { reader.Skip(item_len); } // Set the stream position to the next 4 byte boundary if (reader.GetPosition() % 4 != 0) { reader.Skip(4 - reader.GetPosition() % 4); } } if (content_.find(name) != content_.end()) { throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat, "Tag is repeated in CSA header: " + name); } else { content_[name] = tag2.release(); } } } const CSATag& CSAHeader::GetTag(const std::string& name) const { Content::const_iterator found = content_.find(name); if (found == content_.end()) { throw Orthanc::OrthancException(Orthanc::ErrorCode_InexistentItem); } else { assert(found->second != NULL); return *found->second; } } void CSAHeader::ListTags(std::list<std::string>& tags) const { for (Content::const_iterator it = content_.begin(); it != content_.end(); ++it) { assert(it->second != NULL); tags.push_back(it->first); } } bool CSAHeader::ParseUnsignedInteger32(uint32_t& target, const std::string& tagName) const { Content::const_iterator found = content_.find(tagName); if (found == content_.end()) { return false; } else if (found->second->GetSize() != 1) { return false; } else { return found->second->ParseUnsignedInteger32(target, 0); } } CSATag& CSAHeader::AddTag(const std::string& name, const std::string& vr) { if (content_.find(name) != content_.end()) { throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange, "Tag already exists: " + name); } else { CSATag* tag = new CSATag(vr); content_[name] = tag; return *tag; } } void CSAHeader::AddValue(const std::string& tagName, const std::string& value) { Content::iterator found = content_.find(tagName); if (found == content_.end()) { throw Orthanc::OrthancException(Orthanc::ErrorCode_InexistentItem); } else { assert(found->second != NULL); found->second->AddValue(value); } } }