orthanc-stone: OrthancStone/Sources/Loaders/OrthancMultiframeVolumeLoader.cpp annotate

annotate OrthancStone/Sources/Loaders/OrthancMultiframeVolumeLoader.cpp @ 2077:07964689cb0b

upgrade to year 2023

author	Sebastien Jodogne <s.jodogne@gmail.com>
date	Wed, 12 Jul 2023 21:20:22 +0200
parents	990f396484b1
children	c23eef785569

rev	line source
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	1 /**
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	2 * Stone of Orthanc
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	3 * Copyright (C) 2012-2016 Sebastien Jodogne, Medical Physics
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	4 * Department, University Hospital of Liege, Belgium
2077 07964689cb0b upgrade to year 2023 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2076 diff changeset	5 * Copyright (C) 2017-2023 Osimis S.A., Belgium
07964689cb0b upgrade to year 2023 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2076 diff changeset	6 * Copyright (C) 2021-2023 Sebastien Jodogne, ICTEAM UCLouvain, Belgium
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	7 *
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	8 * This program is free software: you can redistribute it and/or
1598 8563ea5d8ae4 relicensing some files, cf. osimis bm26 and chu agreement on 2020-05-20 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1596 diff changeset	9 * modify it under the terms of the GNU Lesser General Public License
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	10 * as published by the Free Software Foundation, either version 3 of
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	11 * the License, or (at your option) any later version.
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	12 *
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	13 * This program is distributed in the hope that it will be useful, but
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	14 * WITHOUT ANY WARRANTY; without even the implied warranty of
1598 8563ea5d8ae4 relicensing some files, cf. osimis bm26 and chu agreement on 2020-05-20 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1596 diff changeset	15 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
8563ea5d8ae4 relicensing some files, cf. osimis bm26 and chu agreement on 2020-05-20 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1596 diff changeset	16 * Lesser General Public License for more details.
1596 4fb8fdf03314 removed annoying whitespace Sebastien Jodogne <s.jodogne@gmail.com> parents: 1571 diff changeset	17 *
1598 8563ea5d8ae4 relicensing some files, cf. osimis bm26 and chu agreement on 2020-05-20 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1596 diff changeset	18 * You should have received a copy of the GNU Lesser General Public
8563ea5d8ae4 relicensing some files, cf. osimis bm26 and chu agreement on 2020-05-20 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1596 diff changeset	19 * License along with this program. If not, see
8563ea5d8ae4 relicensing some files, cf. osimis bm26 and chu agreement on 2020-05-20 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1596 diff changeset	20 * <http://www.gnu.org/licenses/>.
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	21 **/
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	22
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	23
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	24 #include "OrthancMultiframeVolumeLoader.h"
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	25
1455 30deba7bc8e2 simplifying include_directories Sebastien Jodogne <s.jodogne@gmail.com> parents: 1381 diff changeset	26 #include <Endianness.h>
30deba7bc8e2 simplifying include_directories Sebastien Jodogne <s.jodogne@gmail.com> parents: 1381 diff changeset	27 #include <Toolbox.h>
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	28
1751 946eb7200b82 FastParseVector usage + timing instrumentation guarded by STONE_TIME_BLOCKING_OPS Benjamin Golinvaux <bgo@osimis.io> parents: 1739 diff changeset	29 #if STONE_TIME_BLOCKING_OPS
946eb7200b82 FastParseVector usage + timing instrumentation guarded by STONE_TIME_BLOCKING_OPS Benjamin Golinvaux <bgo@osimis.io> parents: 1739 diff changeset	30 # include <boost/date_time/posix_time/posix_time.hpp>
946eb7200b82 FastParseVector usage + timing instrumentation guarded by STONE_TIME_BLOCKING_OPS Benjamin Golinvaux <bgo@osimis.io> parents: 1739 diff changeset	31 #endif
946eb7200b82 FastParseVector usage + timing instrumentation guarded by STONE_TIME_BLOCKING_OPS Benjamin Golinvaux <bgo@osimis.io> parents: 1739 diff changeset	32
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	33 namespace OrthancStone
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	34 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	35 class OrthancMultiframeVolumeLoader::LoadRTDoseGeometry : public LoaderStateMachine::State
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	36 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	37 private:
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	38 std::unique_ptr<Orthanc::DicomMap> dicom_;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	39
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	40 public:
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	41 LoadRTDoseGeometry(OrthancMultiframeVolumeLoader& that,
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	42 Orthanc::DicomMap* dicom) :
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	43 State(that),
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	44 dicom_(dicom)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	45 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	46 if (dicom == NULL)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	47 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	48 throw Orthanc::OrthancException(Orthanc::ErrorCode_NullPointer);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	49 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	50
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	51 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	52
1640 52b8b96cb55f cleaning namespaces Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	53 virtual void Handle(const OrthancRestApiCommand::SuccessMessage& message) ORTHANC_OVERRIDE
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	54 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	55 // Complete the DICOM tags with just-received "Grid Frame Offset Vector"
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	56 std::string s = Orthanc::Toolbox::StripSpaces(message.GetAnswer());
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	57 dicom_->SetValue(Orthanc::DICOM_TAG_GRID_FRAME_OFFSET_VECTOR, s, false);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	58
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	59 GetLoader<OrthancMultiframeVolumeLoader>().SetGeometry(*dicom_);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	60 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	61 };
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	62
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	63
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	64 static std::string GetSopClassUid(const Orthanc::DicomMap& dicom)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	65 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	66 std::string s;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	67 if (!dicom.LookupStringValue(s, Orthanc::DICOM_TAG_SOP_CLASS_UID, false))
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	68 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	69 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat,
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	70 "DICOM file without SOP class UID");
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	71 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	72 else
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	73 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	74 return s;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	75 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	76 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	77
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	78
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	79 class OrthancMultiframeVolumeLoader::LoadGeometry : public State
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	80 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	81 public:
1571 85e117739eca cppcheck Sebastien Jodogne <s.jodogne@gmail.com> parents: 1512 diff changeset	82 explicit LoadGeometry(OrthancMultiframeVolumeLoader& that) :
85e117739eca cppcheck Sebastien Jodogne <s.jodogne@gmail.com> parents: 1512 diff changeset	83 State(that)
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	84 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	85 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	86
1640 52b8b96cb55f cleaning namespaces Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	87 virtual void Handle(const OrthancRestApiCommand::SuccessMessage& message)
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	88 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	89 OrthancMultiframeVolumeLoader& loader = GetLoader<OrthancMultiframeVolumeLoader>();
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	90
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	91 Json::Value body;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	92 message.ParseJsonBody(body);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	93
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	94 if (body.type() != Json::objectValue)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	95 {
1758 4ee11b8773e2 Better error description in case of network json payload errors Benjamin Golinvaux <bgo@osimis.io> parents: 1757 diff changeset	96 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadJson, "JSON body should be an object value");
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	97 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	98
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	99 std::unique_ptr<Orthanc::DicomMap> dicom(new Orthanc::DicomMap);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	100 dicom->FromDicomAsJson(body);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	101
1640 52b8b96cb55f cleaning namespaces Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	102 if (StringToSopClassUid(GetSopClassUid(*dicom)) == SopClassUid_RTDose)
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	103 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	104 // Download the "Grid Frame Offset Vector" DICOM tag, that is
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	105 // mandatory for RT-DOSE, but is too long to be returned by default
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	106
1640 52b8b96cb55f cleaning namespaces Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	107 std::unique_ptr<OrthancRestApiCommand> command(new OrthancRestApiCommand);
1831 a03260db8835 Added SetCallerName to IOracleCommand to ease debugging Benjamin Golinvaux <bgo@osimis.io> parents: 1758 diff changeset	108 command->SetCallerName("OrthancMultiframeVolumeLoader::LoadGeometry");
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	109 command->SetUri("/instances/" + loader.GetInstanceId() + "/content/" +
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	110 Orthanc::DICOM_TAG_GRID_FRAME_OFFSET_VECTOR.Format());
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	111 command->AcquirePayload(new LoadRTDoseGeometry(loader, dicom.release()));
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	112
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	113 Schedule(command.release());
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	114 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	115 else
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	116 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	117 loader.SetGeometry(*dicom);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	118 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	119 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	120 };
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	121
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	122 class OrthancMultiframeVolumeLoader::LoadTransferSyntax : public State
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	123 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	124 public:
1571 85e117739eca cppcheck Sebastien Jodogne <s.jodogne@gmail.com> parents: 1512 diff changeset	125 explicit LoadTransferSyntax(OrthancMultiframeVolumeLoader& that) :
85e117739eca cppcheck Sebastien Jodogne <s.jodogne@gmail.com> parents: 1512 diff changeset	126 State(that)
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	127 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	128 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	129
1640 52b8b96cb55f cleaning namespaces Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	130 virtual void Handle(const OrthancRestApiCommand::SuccessMessage& message)
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	131 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	132 GetLoader<OrthancMultiframeVolumeLoader>().SetTransferSyntax(message.GetAnswer());
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	133 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	134 };
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	135
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	136 class OrthancMultiframeVolumeLoader::LoadUncompressedPixelData : public State
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	137 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	138 public:
1571 85e117739eca cppcheck Sebastien Jodogne <s.jodogne@gmail.com> parents: 1512 diff changeset	139 explicit LoadUncompressedPixelData(OrthancMultiframeVolumeLoader& that) :
85e117739eca cppcheck Sebastien Jodogne <s.jodogne@gmail.com> parents: 1512 diff changeset	140 State(that)
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	141 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	142 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	143
1640 52b8b96cb55f cleaning namespaces Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	144 virtual void Handle(const OrthancRestApiCommand::SuccessMessage& message)
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	145 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	146 GetLoader<OrthancMultiframeVolumeLoader>().SetUncompressedPixelData(message.GetAnswer());
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	147 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	148 };
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	149
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	150 const std::string& OrthancMultiframeVolumeLoader::GetInstanceId() const
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	151 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	152 if (IsActive())
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	153 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	154 return instanceId_;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	155 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	156 else
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	157 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	158 LOG(ERROR) << "OrthancMultiframeVolumeLoader::GetInstanceId(): (!IsActive())";
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	159 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	160 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	161 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	162
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	163 void OrthancMultiframeVolumeLoader::ScheduleFrameDownloads()
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	164 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	165 if (transferSyntaxUid_.empty() \|\|
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	166 !volume_->HasGeometry())
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	167 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	168 return;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	169 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	170 /*
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	171 1.2.840.10008.1.2 Implicit VR Endian: Default Transfer Syntax for DICOM
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	172 1.2.840.10008.1.2.1 Explicit VR Little Endian
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	173 1.2.840.10008.1.2.2 Explicit VR Big Endian
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	174
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	175 See https://www.dicomlibrary.com/dicom/transfer-syntax/
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	176 */
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	177 if (transferSyntaxUid_ == "1.2.840.10008.1.2" \|\|
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	178 transferSyntaxUid_ == "1.2.840.10008.1.2.1" \|\|
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	179 transferSyntaxUid_ == "1.2.840.10008.1.2.2")
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	180 {
1640 52b8b96cb55f cleaning namespaces Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	181 std::unique_ptr<OrthancRestApiCommand> command(new OrthancRestApiCommand);
1831 a03260db8835 Added SetCallerName to IOracleCommand to ease debugging Benjamin Golinvaux <bgo@osimis.io> parents: 1758 diff changeset	182 command->SetCallerName("OrthancMultiframeVolumeLoader::ScheduleFrameDownloads");
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	183 command->SetHttpHeader("Accept-Encoding", "gzip");
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	184 command->SetUri("/instances/" + instanceId_ + "/content/" +
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	185 Orthanc::DICOM_TAG_PIXEL_DATA.Format() + "/0");
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	186 command->AcquirePayload(new LoadUncompressedPixelData(*this));
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	187 Schedule(command.release());
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	188 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	189 else
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	190 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	191 throw Orthanc::OrthancException(
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	192 Orthanc::ErrorCode_NotImplemented,
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	193 "No support for multiframe instances with transfer syntax: " + transferSyntaxUid_);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	194 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	195 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	196
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	197 void OrthancMultiframeVolumeLoader::SetTransferSyntax(const std::string& transferSyntax)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	198 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	199 transferSyntaxUid_ = Orthanc::Toolbox::StripSpaces(transferSyntax);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	200 ScheduleFrameDownloads();
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	201 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	202
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	203 void OrthancMultiframeVolumeLoader::SetGeometry(const Orthanc::DicomMap& dicom)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	204 {
1640 52b8b96cb55f cleaning namespaces Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	205 DicomInstanceParameters parameters(dicom);
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	206 volume_->SetDicomParameters(parameters);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	207
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	208 Orthanc::PixelFormat format;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	209 if (!parameters.GetImageInformation().ExtractPixelFormat(format, true))
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	210 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	211 throw Orthanc::OrthancException(Orthanc::ErrorCode_NotImplemented);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	212 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	213
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	214 double spacingZ;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	215 switch (parameters.GetSopClassUid())
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	216 {
1640 52b8b96cb55f cleaning namespaces Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	217 case SopClassUid_RTDose:
2076 990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	218 if (!parameters.ComputeFrameOffsetsSpacing(spacingZ))
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	219 {
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	220 throw Orthanc::OrthancException(Orthanc::ErrorCode_NotImplemented,
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	221 "Cannot load RT-DOSE with incorrect GridFrameOffsetVector (3004,000C)");
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	222 }
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	223 break;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	224
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	225 default:
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	226 throw Orthanc::OrthancException(
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	227 Orthanc::ErrorCode_NotImplemented,
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	228 "No support for multiframe instances with SOP class UID: " + GetSopClassUid(dicom));
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	229 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	230
2076 990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	231 isReversedFrameOffsets_ = parameters.IsReversedFrameOffsets();
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	232
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	233 const unsigned int width = parameters.GetImageInformation().GetWidth();
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	234 const unsigned int height = parameters.GetImageInformation().GetHeight();
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	235 const unsigned int depth = parameters.GetImageInformation().GetNumberOfFrames();
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	236
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	237 {
1640 52b8b96cb55f cleaning namespaces Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	238 VolumeImageGeometry geometry;
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	239 geometry.SetSizeInVoxels(width, height, depth);
2076 990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	240 geometry.SetAxialGeometry(parameters.GetMultiFrameGeometry());
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	241 geometry.SetVoxelDimensions(parameters.GetPixelSpacingX(),
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	242 parameters.GetPixelSpacingY(), spacingZ);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	243 volume_->Initialize(geometry, format, true /* Do compute range */);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	244 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	245
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	246 volume_->GetPixelData().Clear();
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	247
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	248 ScheduleFrameDownloads();
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	249
1640 52b8b96cb55f cleaning namespaces Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	250 BroadcastMessage(DicomVolumeImage::GeometryReadyMessage(*volume_));
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	251 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	252
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	253
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	254 ORTHANC_FORCE_INLINE
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	255 static void CopyPixel(uint32_t& target, const void* source)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	256 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	257 // TODO - check alignement?
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	258 target = le32toh(reinterpret_cast<const uint32_t>(source));
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	259 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	260
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	261 ORTHANC_FORCE_INLINE
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	262 static void CopyPixel(uint16_t& target, const void* source)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	263 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	264 // TODO - check alignement?
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	265 target = le16toh(reinterpret_cast<const uint16_t>(source));
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	266 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	267
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	268 ORTHANC_FORCE_INLINE
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	269 static void CopyPixel(int16_t& target, const void* source)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	270 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	271 // byte swapping is the same for unsigned and signed integers
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	272 // (the sign bit is always stored with the MSByte)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	273 uint16_t* targetUp = reinterpret_cast<uint16_t*>(&target);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	274 CopyPixel(*targetUp, source);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	275 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	276
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	277 template <typename T>
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	278 void OrthancMultiframeVolumeLoader::CopyPixelDataAndComputeDistribution(
1757 28979b77ce90 Suppressed one pass on the dose distribution computation bgo@SHARKNADO.localdomain parents: 1751 diff changeset	279 const std::string& pixelData, std::map<T, PixelCount>& distribution)
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	280 {
1751 946eb7200b82 FastParseVector usage + timing instrumentation guarded by STONE_TIME_BLOCKING_OPS Benjamin Golinvaux <bgo@osimis.io> parents: 1739 diff changeset	281 #if STONE_TIME_BLOCKING_OPS
946eb7200b82 FastParseVector usage + timing instrumentation guarded by STONE_TIME_BLOCKING_OPS Benjamin Golinvaux <bgo@osimis.io> parents: 1739 diff changeset	282 boost::posix_time::ptime timerStart = boost::posix_time::microsec_clock::universal_time();
946eb7200b82 FastParseVector usage + timing instrumentation guarded by STONE_TIME_BLOCKING_OPS Benjamin Golinvaux <bgo@osimis.io> parents: 1739 diff changeset	283 #endif
946eb7200b82 FastParseVector usage + timing instrumentation guarded by STONE_TIME_BLOCKING_OPS Benjamin Golinvaux <bgo@osimis.io> parents: 1739 diff changeset	284
1640 52b8b96cb55f cleaning namespaces Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	285 ImageBuffer3D& target = volume_->GetPixelData();
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	286
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	287 const unsigned int bpp = target.GetBytesPerPixel();
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	288 const unsigned int width = target.GetWidth();
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	289 const unsigned int height = target.GetHeight();
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	290 const unsigned int depth = target.GetDepth();
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	291
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	292 if (pixelData.size() != bpp * width * height * depth)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	293 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	294 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat,
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	295 "The pixel data has not the proper size");
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	296 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	297
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	298 if (pixelData.empty())
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	299 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	300 return;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	301 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	302
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	303 // first pass to initialize map
1757 28979b77ce90 Suppressed one pass on the dose distribution computation bgo@SHARKNADO.localdomain parents: 1751 diff changeset	304 #if 0
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	305 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	306 const uint8_t* source = reinterpret_cast<const uint8_t*>(pixelData.c_str());
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	307
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	308 for (unsigned int z = 0; z < depth; z++)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	309 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	310 for (unsigned int y = 0; y < height; y++)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	311 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	312 for (unsigned int x = 0; x < width; x++)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	313 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	314 T value;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	315 CopyPixel(value, source);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	316 distribution[value] = 0;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	317 source += bpp;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	318 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	319 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	320 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	321 }
1757 28979b77ce90 Suppressed one pass on the dose distribution computation bgo@SHARKNADO.localdomain parents: 1751 diff changeset	322 #endif
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	323
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	324 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	325 const uint8_t* source = reinterpret_cast<const uint8_t*>(pixelData.c_str());
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	326
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	327 for (unsigned int z = 0; z < depth; z++)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	328 {
2076 990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	329 unsigned targetZ;
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	330 if (isReversedFrameOffsets_)
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	331 {
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	332 targetZ = depth - 1 - z;
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	333 }
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	334 else
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	335 {
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	336 targetZ = z;
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	337 }
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	338
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	339 ImageBuffer3D::SliceWriter writer(target, VolumeProjection_Axial, targetZ);
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	340
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	341 assert(writer.GetAccessor().GetWidth() == width &&
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	342 writer.GetAccessor().GetHeight() == height);
1362 104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	343 #if 0
104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	344 for (unsigned int y = 0; y < height; y++)
104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	345 {
104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	346 assert(sizeof(T) == Orthanc::GetBytesPerPixel(target.GetFormat()));
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	347
1362 104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	348 T* target = reinterpret_cast<T*>(writer.GetAccessor().GetRow(y));
104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	349
104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	350 for (unsigned int x = 0; x < width; x++)
104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	351 {
104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	352 CopyPixel(*target, source);
104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	353
104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	354 distribution[*target] += 1;
104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	355
104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	356 target++;
104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	357 source += bpp;
104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	358 }
104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	359 }
104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	360 #else
104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	361 // optimized version (fixed) as of 2020-04-15
1349 b1e6bef86955 Tentative loop optimization in CopyPixelDataAndComputeDistribution Benjamin Golinvaux <bgo@osimis.io> parents: 1337 diff changeset	362 unsigned int pitch = writer.GetAccessor().GetPitch();
1362 104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	363 T* targetAddrLine = reinterpret_cast<T*>(writer.GetAccessor().GetRow(0));
1349 b1e6bef86955 Tentative loop optimization in CopyPixelDataAndComputeDistribution Benjamin Golinvaux <bgo@osimis.io> parents: 1337 diff changeset	364 assert(sizeof(T) == Orthanc::GetBytesPerPixel(target.GetFormat()));
b1e6bef86955 Tentative loop optimization in CopyPixelDataAndComputeDistribution Benjamin Golinvaux <bgo@osimis.io> parents: 1337 diff changeset	365
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	366 for (unsigned int y = 0; y < height; y++)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	367 {
1362 104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	368 T* targetAddrPix = targetAddrLine;
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	369 for (unsigned int x = 0; x < width; x++)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	370 {
1362 104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	371 CopyPixel(*targetAddrPix, source);
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	372
1757 28979b77ce90 Suppressed one pass on the dose distribution computation bgo@SHARKNADO.localdomain parents: 1751 diff changeset	373 distribution[*targetAddrPix].count_ += 1;
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	374
1362 104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	375 targetAddrPix++;
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	376 source += bpp;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	377 }
1362 104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	378 uint8_t* targetAddrLineBytes = reinterpret_cast<uint8_t*>(targetAddrLine) + pitch;
104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	379 targetAddrLine = reinterpret_cast<T*>(targetAddrLineBytes);
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	380 }
1362 104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	381 #endif
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	382 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	383 }
1751 946eb7200b82 FastParseVector usage + timing instrumentation guarded by STONE_TIME_BLOCKING_OPS Benjamin Golinvaux <bgo@osimis.io> parents: 1739 diff changeset	384 #if STONE_TIME_BLOCKING_OPS
946eb7200b82 FastParseVector usage + timing instrumentation guarded by STONE_TIME_BLOCKING_OPS Benjamin Golinvaux <bgo@osimis.io> parents: 1739 diff changeset	385 boost::posix_time::ptime timerEnd = boost::posix_time::microsec_clock::universal_time();
946eb7200b82 FastParseVector usage + timing instrumentation guarded by STONE_TIME_BLOCKING_OPS Benjamin Golinvaux <bgo@osimis.io> parents: 1739 diff changeset	386 boost::posix_time::time_duration duration = timerEnd - timerStart;
946eb7200b82 FastParseVector usage + timing instrumentation guarded by STONE_TIME_BLOCKING_OPS Benjamin Golinvaux <bgo@osimis.io> parents: 1739 diff changeset	387 int64_t durationMs = duration.total_milliseconds();
946eb7200b82 FastParseVector usage + timing instrumentation guarded by STONE_TIME_BLOCKING_OPS Benjamin Golinvaux <bgo@osimis.io> parents: 1739 diff changeset	388 LOG(WARNING) << "OrthancMultiframeVolumeLoader::CopyPixelDataAndComputeDistribution took " << durationMs << " ms";
946eb7200b82 FastParseVector usage + timing instrumentation guarded by STONE_TIME_BLOCKING_OPS Benjamin Golinvaux <bgo@osimis.io> parents: 1739 diff changeset	389 #endif
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	390 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	391
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	392 template <typename T>
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	393 void OrthancMultiframeVolumeLoader::ComputeMinMaxWithOutlierRejection(
1757 28979b77ce90 Suppressed one pass on the dose distribution computation bgo@SHARKNADO.localdomain parents: 1751 diff changeset	394 const std::map<T, PixelCount>& distribution)
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	395 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	396 if (distribution.size() == 0)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	397 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	398 LOG(ERROR) << "ComputeMinMaxWithOutlierRejection -- Volume image empty.";
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	399 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	400 else
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	401 {
1640 52b8b96cb55f cleaning namespaces Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	402 ImageBuffer3D& target = volume_->GetPixelData();
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	403
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	404 const uint64_t width = target.GetWidth();
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	405 const uint64_t height = target.GetHeight();
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	406 const uint64_t depth = target.GetDepth();
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	407 const uint64_t voxelCount = width * height * depth;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	408
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	409 // now that we have distribution[pixelValue] == numberOfPixelsWithValue
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	410 // compute number of values and check (assertion) that it is equal to
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	411 // width * height * depth
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	412 {
1757 28979b77ce90 Suppressed one pass on the dose distribution computation bgo@SHARKNADO.localdomain parents: 1751 diff changeset	413 typename std::map<T, PixelCount>::const_iterator it = distribution.begin();
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	414 uint64_t totalCount = 0;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	415 distributionRawMin_ = static_cast<float>(it->first);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	416
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	417 while (it != distribution.end())
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	418 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	419 T pixelValue = it->first;
1757 28979b77ce90 Suppressed one pass on the dose distribution computation bgo@SHARKNADO.localdomain parents: 1751 diff changeset	420 totalCount += it->second.count_;
1571 85e117739eca cppcheck Sebastien Jodogne <s.jodogne@gmail.com> parents: 1512 diff changeset	421 ++it;
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	422 if (it == distribution.end())
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	423 distributionRawMax_ = static_cast<float>(pixelValue);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	424 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	425 LOG(INFO) << "Volume image. First distribution value = "
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	426 << static_cast<float>(distributionRawMin_)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	427 << " \| Last distribution value = "
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	428 << static_cast<float>(distributionRawMax_);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	429
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	430 if (totalCount != voxelCount)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	431 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	432 LOG(ERROR) << "Internal error in dose distribution computation. TC ("
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	433 << totalCount << ") != VoxC (" << voxelCount;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	434 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	435 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	436 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	437
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	438 // compute the number of voxels to reject at each end of the distribution
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	439 uint64_t endRejectionCount = static_cast<uint64_t>(
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	440 outliersHalfRejectionRate_ * voxelCount);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	441
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	442 if (endRejectionCount > voxelCount)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	443 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	444 LOG(ERROR) << "Internal error in dose distribution computation."
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	445 << " endRejectionCount = " << endRejectionCount
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	446 << " \| voxelCount = " << voxelCount;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	447 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	448 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	449
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	450 // this will contain the actual distribution minimum after outlier
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	451 // rejection
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	452 T resultMin = 0;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	453
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	454 // then start from start and remove pixel values up to
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	455 // endRejectionCount voxels rejected
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	456 {
1757 28979b77ce90 Suppressed one pass on the dose distribution computation bgo@SHARKNADO.localdomain parents: 1751 diff changeset	457 typename std::map<T, PixelCount>::const_iterator it = distribution.begin();
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	458
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	459 uint64_t currentCount = 0;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	460
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	461 while (it != distribution.end())
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	462 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	463 T pixelValue = it->first;
1757 28979b77ce90 Suppressed one pass on the dose distribution computation bgo@SHARKNADO.localdomain parents: 1751 diff changeset	464 uint64_t count = it->second.count_;
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	465
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	466 // if this pixelValue crosses the rejection threshold, let's set it
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	467 // and exit the loop
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	468 if ((currentCount <= endRejectionCount) &&
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	469 (currentCount + count > endRejectionCount))
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	470 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	471 resultMin = pixelValue;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	472 break;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	473 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	474 else
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	475 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	476 currentCount += count;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	477 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	478 // and continue walking along the distribution
1571 85e117739eca cppcheck Sebastien Jodogne <s.jodogne@gmail.com> parents: 1512 diff changeset	479 ++it;
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	480 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	481 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	482
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	483 // this will contain the actual distribution maximum after outlier
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	484 // rejection
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	485 T resultMax = 0;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	486 // now start from END and remove pixel values up to
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	487 // endRejectionCount voxels rejected
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	488 {
1757 28979b77ce90 Suppressed one pass on the dose distribution computation bgo@SHARKNADO.localdomain parents: 1751 diff changeset	489 typename std::map<T, PixelCount>::const_reverse_iterator it = distribution.rbegin();
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	490
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	491 uint64_t currentCount = 0;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	492
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	493 while (it != distribution.rend())
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	494 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	495 T pixelValue = it->first;
1757 28979b77ce90 Suppressed one pass on the dose distribution computation bgo@SHARKNADO.localdomain parents: 1751 diff changeset	496 uint64_t count = it->second.count_;
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	497
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	498 if ((currentCount <= endRejectionCount) &&
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	499 (currentCount + count > endRejectionCount))
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	500 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	501 resultMax = pixelValue;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	502 break;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	503 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	504 else
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	505 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	506 currentCount += count;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	507 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	508 // and continue walking along the distribution
1571 85e117739eca cppcheck Sebastien Jodogne <s.jodogne@gmail.com> parents: 1512 diff changeset	509 ++it;
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	510 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	511 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	512 if (resultMin > resultMax)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	513 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	514 LOG(ERROR) << "Internal error in dose distribution computation! " <<
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	515 "resultMin (" << resultMin << ") > resultMax (" << resultMax << ")";
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	516 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	517 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	518 computedDistributionMin_ = static_cast<float>(resultMin);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	519 computedDistributionMax_ = static_cast<float>(resultMax);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	520 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	521 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	522
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	523 template <typename T>
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	524 void OrthancMultiframeVolumeLoader::CopyPixelDataAndComputeMinMax(
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	525 const std::string& pixelData)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	526 {
1757 28979b77ce90 Suppressed one pass on the dose distribution computation bgo@SHARKNADO.localdomain parents: 1751 diff changeset	527 std::map<T, PixelCount> distribution;
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	528 CopyPixelDataAndComputeDistribution(pixelData, distribution);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	529 ComputeMinMaxWithOutlierRejection(distribution);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	530 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	531
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	532 void OrthancMultiframeVolumeLoader::SetUncompressedPixelData(const std::string& pixelData)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	533 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	534 switch (volume_->GetPixelData().GetFormat())
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	535 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	536 case Orthanc::PixelFormat_Grayscale32:
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	537 CopyPixelDataAndComputeMinMax<uint32_t>(pixelData);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	538 break;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	539 case Orthanc::PixelFormat_Grayscale16:
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	540 CopyPixelDataAndComputeMinMax<uint16_t>(pixelData);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	541 break;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	542 case Orthanc::PixelFormat_SignedGrayscale16:
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	543 CopyPixelDataAndComputeMinMax<int16_t>(pixelData);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	544 break;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	545 default:
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	546 throw Orthanc::OrthancException(Orthanc::ErrorCode_NotImplemented);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	547 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	548
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	549 volume_->IncrementRevision();
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	550
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	551 pixelDataLoaded_ = true;
1640 52b8b96cb55f cleaning namespaces Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	552 BroadcastMessage(DicomVolumeImage::ContentUpdatedMessage(*volume_));
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	553 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	554
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	555 bool OrthancMultiframeVolumeLoader::HasGeometry() const
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	556 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	557 return volume_->HasGeometry();
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	558 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	559
1640 52b8b96cb55f cleaning namespaces Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	560 const VolumeImageGeometry& OrthancMultiframeVolumeLoader::GetImageGeometry() const
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	561 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	562 return volume_->GetGeometry();
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	563 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	564
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	565 OrthancMultiframeVolumeLoader::OrthancMultiframeVolumeLoader(
1640 52b8b96cb55f cleaning namespaces Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	566 ILoadersContext& loadersContext,
52b8b96cb55f cleaning namespaces Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	567 boost::shared_ptr<DicomVolumeImage> volume,
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	568 float outliersHalfRejectionRate)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	569 : LoaderStateMachine(loadersContext)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	570 , volume_(volume)
2076 990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	571 , isReversedFrameOffsets_(false)
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	572 , pixelDataLoaded_(false)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	573 , outliersHalfRejectionRate_(outliersHalfRejectionRate)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	574 , distributionRawMin_(0)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	575 , distributionRawMax_(0)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	576 , computedDistributionMin_(0)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	577 , computedDistributionMax_(0)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	578 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	579 if (volume.get() == NULL)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	580 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	581 throw Orthanc::OrthancException(Orthanc::ErrorCode_NullPointer);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	582 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	583 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	584
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	585
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	586 boost::shared_ptr<OrthancMultiframeVolumeLoader>
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	587 OrthancMultiframeVolumeLoader::Create(
1640 52b8b96cb55f cleaning namespaces Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	588 ILoadersContext& loadersContext,
52b8b96cb55f cleaning namespaces Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	589 boost::shared_ptr<DicomVolumeImage> volume,
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	590 float outliersHalfRejectionRate /= 0.0005/)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	591 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	592 boost::shared_ptr<OrthancMultiframeVolumeLoader> obj(
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	593 new OrthancMultiframeVolumeLoader(
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	594 loadersContext,
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	595 volume,
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	596 outliersHalfRejectionRate));
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	597 obj->LoaderStateMachine::PostConstructor();
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	598 return obj;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	599 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	600
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	601 OrthancMultiframeVolumeLoader::~OrthancMultiframeVolumeLoader()
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	602 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	603 LOG(TRACE) << "OrthancMultiframeVolumeLoader::~OrthancMultiframeVolumeLoader()";
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	604 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	605
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	606 void OrthancMultiframeVolumeLoader::GetDistributionMinMax
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	607 (float& minValue, float& maxValue) const
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	608 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	609 if (distributionRawMin_ == 0 && distributionRawMax_ == 0)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	610 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	611 LOG(WARNING) << "GetDistributionMinMaxWithOutliersRejection called before computation!";
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	612 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	613 minValue = distributionRawMin_;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	614 maxValue = distributionRawMax_;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	615 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	616
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	617 void OrthancMultiframeVolumeLoader::GetDistributionMinMaxWithOutliersRejection
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	618 (float& minValue, float& maxValue) const
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	619 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	620 if (computedDistributionMin_ == 0 && computedDistributionMax_ == 0)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	621 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	622 LOG(WARNING) << "GetDistributionMinMaxWithOutliersRejection called before computation!";
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	623 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	624 minValue = computedDistributionMin_;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	625 maxValue = computedDistributionMax_;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	626 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	627
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	628 void OrthancMultiframeVolumeLoader::LoadInstance(const std::string& instanceId)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	629 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	630 Start();
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	631
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	632 instanceId_ = instanceId;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	633
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	634 {
1640 52b8b96cb55f cleaning namespaces Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	635 std::unique_ptr<OrthancRestApiCommand> command(new OrthancRestApiCommand);
1831 a03260db8835 Added SetCallerName to IOracleCommand to ease debugging Benjamin Golinvaux <bgo@osimis.io> parents: 1758 diff changeset	636 command->SetCallerName("OrthancMultiframeVolumeLoader::LoadInstance");
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	637 command->SetHttpHeader("Accept-Encoding", "gzip");
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	638 command->SetUri("/instances/" + instanceId + "/tags");
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	639 command->AcquirePayload(new LoadGeometry(*this));
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	640 Schedule(command.release());
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	641 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	642
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	643 {
1640 52b8b96cb55f cleaning namespaces Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	644 std::unique_ptr<OrthancRestApiCommand> command(new OrthancRestApiCommand);
1831 a03260db8835 Added SetCallerName to IOracleCommand to ease debugging Benjamin Golinvaux <bgo@osimis.io> parents: 1758 diff changeset	645 command->SetCallerName("OrthancMultiframeVolumeLoader::LoadInstance");
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	646 command->SetUri("/instances/" + instanceId + "/metadata/TransferSyntax");
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	647 command->AcquirePayload(new LoadTransferSyntax(*this));
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	648 Schedule(command.release());
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	649 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	650 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	651 }

Mercurial > hg > orthanc-stone

annotate OrthancStone/Sources/Loaders/OrthancMultiframeVolumeLoader.cpp @ 2077:07964689cb0b