orthanc-stone: Framework/Loaders/OrthancMultiframeVolumeLoader.cpp annotate

annotate Framework/Loaders/OrthancMultiframeVolumeLoader.cpp @ 1487:75ac66d5f4b2 bgo-timing-tests

Added timing test for dose loading and histogram computation

author	Benjamin Golinvaux <bgo@osimis.io>
date	Tue, 09 Jun 2020 09:50:03 +0200
parents	f4a06ad1580b
children

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
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	5 * Copyright (C) 2017-2020 Osimis S.A., Belgium
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	6 *
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	7 * This program is free software: you can redistribute it and/or
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	8 * modify it under the terms of the GNU Affero General Public License
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	9 * 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	10 * 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	11 *
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	12 * 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	13 * WITHOUT ANY WARRANTY; without even the implied warranty of
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	15 * Affero General Public License for more details.
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	16 *
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	17 * You should have received a copy of the GNU Affero General Public License
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	18 * along with this program. If not, see <http://www.gnu.org/licenses/>.
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	19 **/
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	20
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 #include "OrthancMultiframeVolumeLoader.h"
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 <Core/Endianness.h>
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	25 #include <Core/Toolbox.h>
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	26
1487 75ac66d5f4b2 Added timing test for dose loading and histogram computation Benjamin Golinvaux <bgo@osimis.io> parents: 1381 diff changeset	27 #include <boost/date_time/posix_time/posix_time.hpp>
75ac66d5f4b2 Added timing test for dose loading and histogram computation Benjamin Golinvaux <bgo@osimis.io> parents: 1381 diff changeset	28
75ac66d5f4b2 Added timing test for dose loading and histogram computation Benjamin Golinvaux <bgo@osimis.io> parents: 1381 diff changeset	29
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	30 namespace OrthancStone
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	31 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	32 class OrthancMultiframeVolumeLoader::LoadRTDoseGeometry : public LoaderStateMachine::State
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	33 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	34 private:
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	35 std::unique_ptr<Orthanc::DicomMap> dicom_;
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 public:
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	38 LoadRTDoseGeometry(OrthancMultiframeVolumeLoader& that,
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	39 Orthanc::DicomMap* dicom) :
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	40 State(that),
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	41 dicom_(dicom)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	42 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	43 if (dicom == NULL)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	44 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	45 throw Orthanc::OrthancException(Orthanc::ErrorCode_NullPointer);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	46 }
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 }
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 virtual void Handle(const OrthancStone::OrthancRestApiCommand::SuccessMessage& message)
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 // 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	53 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	54 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	55
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	56 GetLoader<OrthancMultiframeVolumeLoader>().SetGeometry(*dicom_);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	57 }
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
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 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	62 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	63 std::string s;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	64 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	65 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	66 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat,
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	67 "DICOM file without SOP class UID");
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 else
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	70 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	71 return s;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	72 }
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
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 class OrthancMultiframeVolumeLoader::LoadGeometry : public State
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 public:
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	79 LoadGeometry(OrthancMultiframeVolumeLoader& that) :
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	80 State(that)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	81 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	82 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	83
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	84 virtual void Handle(const OrthancStone::OrthancRestApiCommand::SuccessMessage& message)
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 OrthancMultiframeVolumeLoader& loader = GetLoader<OrthancMultiframeVolumeLoader>();
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	87
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	88 Json::Value body;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	89 message.ParseJsonBody(body);
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 if (body.type() != Json::objectValue)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	92 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	93 throw Orthanc::OrthancException(Orthanc::ErrorCode_NetworkProtocol);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	94 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	95
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	96 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	97 dicom->FromDicomAsJson(body);
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 if (OrthancStone::StringToSopClassUid(GetSopClassUid(*dicom)) == OrthancStone::SopClassUid_RTDose)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	100 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	101 // 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	102 // 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	103
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	104 std::unique_ptr<OrthancStone::OrthancRestApiCommand> command(new OrthancStone::OrthancRestApiCommand);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	105 command->SetUri("/instances/" + loader.GetInstanceId() + "/content/" +
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	106 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	107 command->AcquirePayload(new LoadRTDoseGeometry(loader, dicom.release()));
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	108
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	109 Schedule(command.release());
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	110 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	111 else
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 loader.SetGeometry(*dicom);
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 }
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
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	118 class OrthancMultiframeVolumeLoader::LoadTransferSyntax : public State
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 public:
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	121 LoadTransferSyntax(OrthancMultiframeVolumeLoader& that) :
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	122 State(that)
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 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	125
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	126 virtual void Handle(const OrthancStone::OrthancRestApiCommand::SuccessMessage& message)
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 GetLoader<OrthancMultiframeVolumeLoader>().SetTransferSyntax(message.GetAnswer());
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	129 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	130 };
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 class OrthancMultiframeVolumeLoader::LoadUncompressedPixelData : public State
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 public:
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	135 LoadUncompressedPixelData(OrthancMultiframeVolumeLoader& that) :
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	136 State(that)
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 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	139
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	140 virtual void Handle(const OrthancStone::OrthancRestApiCommand::SuccessMessage& message)
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 GetLoader<OrthancMultiframeVolumeLoader>().SetUncompressedPixelData(message.GetAnswer());
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	143 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	144 };
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 const std::string& OrthancMultiframeVolumeLoader::GetInstanceId() const
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 if (IsActive())
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 return instanceId_;
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 else
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 LOG(ERROR) << "OrthancMultiframeVolumeLoader::GetInstanceId(): (!IsActive())";
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	155 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	156 }
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
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	159 void OrthancMultiframeVolumeLoader::ScheduleFrameDownloads()
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 if (transferSyntaxUid_.empty() \|\|
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	162 !volume_->HasGeometry())
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	163 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	164 return;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	165 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	166 /*
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	167 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	168 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	169 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	170
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	171 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	172 */
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	173 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	174 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	175 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	176 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	177 std::unique_ptr<OrthancStone::OrthancRestApiCommand> command(new OrthancStone::OrthancRestApiCommand);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	178 command->SetHttpHeader("Accept-Encoding", "gzip");
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	179 command->SetUri("/instances/" + instanceId_ + "/content/" +
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	180 Orthanc::DICOM_TAG_PIXEL_DATA.Format() + "/0");
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	181 command->AcquirePayload(new LoadUncompressedPixelData(*this));
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	182 Schedule(command.release());
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	183 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	184 else
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	185 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	186 throw Orthanc::OrthancException(
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	187 Orthanc::ErrorCode_NotImplemented,
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	188 "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	189 }
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
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	192 void OrthancMultiframeVolumeLoader::SetTransferSyntax(const std::string& transferSyntax)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	193 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	194 transferSyntaxUid_ = Orthanc::Toolbox::StripSpaces(transferSyntax);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	195 ScheduleFrameDownloads();
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
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	198 void OrthancMultiframeVolumeLoader::SetGeometry(const Orthanc::DicomMap& dicom)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	199 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	200 OrthancStone::DicomInstanceParameters parameters(dicom);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	201 volume_->SetDicomParameters(parameters);
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 Orthanc::PixelFormat format;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	204 if (!parameters.GetImageInformation().ExtractPixelFormat(format, true))
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	205 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	206 throw Orthanc::OrthancException(Orthanc::ErrorCode_NotImplemented);
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
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	209 double spacingZ;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	210 switch (parameters.GetSopClassUid())
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	211 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	212 case OrthancStone::SopClassUid_RTDose:
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	213 spacingZ = parameters.GetThickness();
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	214 break;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	215
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	216 default:
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	217 throw Orthanc::OrthancException(
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	218 Orthanc::ErrorCode_NotImplemented,
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	219 "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	220 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	221
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	222 const unsigned int width = parameters.GetImageInformation().GetWidth();
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	223 const unsigned int height = parameters.GetImageInformation().GetHeight();
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	224 const unsigned int depth = parameters.GetImageInformation().GetNumberOfFrames();
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	225
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	226 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	227 OrthancStone::VolumeImageGeometry geometry;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	228 geometry.SetSizeInVoxels(width, height, depth);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	229 geometry.SetAxialGeometry(parameters.GetGeometry());
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	230 geometry.SetVoxelDimensions(parameters.GetPixelSpacingX(),
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	231 parameters.GetPixelSpacingY(), spacingZ);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	232 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	233 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	234
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	235 volume_->GetPixelData().Clear();
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 ScheduleFrameDownloads();
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	238
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	239
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	240
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	241 BroadcastMessage(OrthancStone::DicomVolumeImage::GeometryReadyMessage(*volume_));
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	242 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	243
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 ORTHANC_FORCE_INLINE
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	246 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	247 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	248 // TODO - check alignement?
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	249 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	250 }
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 ORTHANC_FORCE_INLINE
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	253 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	254 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	255 // TODO - check alignement?
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	256 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	257 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	258
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	259 ORTHANC_FORCE_INLINE
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	260 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	261 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	262 // 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	263 // (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	264 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	265 CopyPixel(*targetUp, 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
1487 75ac66d5f4b2 Added timing test for dose loading and histogram computation Benjamin Golinvaux <bgo@osimis.io> parents: 1381 diff changeset	268
75ac66d5f4b2 Added timing test for dose loading and histogram computation Benjamin Golinvaux <bgo@osimis.io> parents: 1381 diff changeset	269 #define TIME_CHECKPOINT(name) \
75ac66d5f4b2 Added timing test for dose loading and histogram computation Benjamin Golinvaux <bgo@osimis.io> parents: 1381 diff changeset	270 { \
75ac66d5f4b2 Added timing test for dose loading and histogram computation Benjamin Golinvaux <bgo@osimis.io> parents: 1381 diff changeset	271 pt::ptime nowTime = pt::second_clock::local_time(); \
75ac66d5f4b2 Added timing test for dose loading and histogram computation Benjamin Golinvaux <bgo@osimis.io> parents: 1381 diff changeset	272 pt::time_duration diff = nowTime - initialTime; \
75ac66d5f4b2 Added timing test for dose loading and histogram computation Benjamin Golinvaux <bgo@osimis.io> parents: 1381 diff changeset	273 uint64_t us = diff.total_microseconds(); \
75ac66d5f4b2 Added timing test for dose loading and histogram computation Benjamin Golinvaux <bgo@osimis.io> parents: 1381 diff changeset	274 timingUSecMap_[name] = us; \
75ac66d5f4b2 Added timing test for dose loading and histogram computation Benjamin Golinvaux <bgo@osimis.io> parents: 1381 diff changeset	275 }
75ac66d5f4b2 Added timing test for dose loading and histogram computation Benjamin Golinvaux <bgo@osimis.io> parents: 1381 diff changeset	276
1337 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(
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	279 const std::string& pixelData, std::map<T,uint64_t>& distribution)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	280 {
1487 75ac66d5f4b2 Added timing test for dose loading and histogram computation Benjamin Golinvaux <bgo@osimis.io> parents: 1381 diff changeset	281 namespace pt = boost::posix_time;
75ac66d5f4b2 Added timing test for dose loading and histogram computation Benjamin Golinvaux <bgo@osimis.io> parents: 1381 diff changeset	282
75ac66d5f4b2 Added timing test for dose loading and histogram computation Benjamin Golinvaux <bgo@osimis.io> parents: 1381 diff changeset	283 pt::ptime initialTime = pt::second_clock::local_time();
75ac66d5f4b2 Added timing test for dose loading and histogram computation Benjamin Golinvaux <bgo@osimis.io> parents: 1381 diff changeset	284
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	285 OrthancStone::ImageBuffer3D& target = volume_->GetPixelData();
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
1487 75ac66d5f4b2 Added timing test for dose loading and histogram computation Benjamin Golinvaux <bgo@osimis.io> parents: 1381 diff changeset	303 TIME_CHECKPOINT("start")
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	304 // first pass to initialize map
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 }
1487 75ac66d5f4b2 Added timing test for dose loading and histogram computation Benjamin Golinvaux <bgo@osimis.io> parents: 1381 diff changeset	322 TIME_CHECKPOINT("map_initialized")
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 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	325
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	326 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	327 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	328 OrthancStone::ImageBuffer3D::SliceWriter writer(target, OrthancStone::VolumeProjection_Axial, z);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	329
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	330 assert(writer.GetAccessor().GetWidth() == width &&
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	331 writer.GetAccessor().GetHeight() == height);
1362 104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	332 #if 0
104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	333 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	334 {
104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	335 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	336
1362 104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	337 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	338
104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	339 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	340 {
104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	341 CopyPixel(*target, source);
104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	342
104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	343 distribution[*target] += 1;
104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	344
104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	345 target++;
104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	346 source += bpp;
104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	347 }
104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	348 }
104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	349 #else
104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	350 // optimized version (fixed) as of 2020-04-15
1349 b1e6bef86955 Tentative loop optimization in CopyPixelDataAndComputeDistribution Benjamin Golinvaux <bgo@osimis.io> parents: 1337 diff changeset	351 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	352 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	353 assert(sizeof(T) == Orthanc::GetBytesPerPixel(target.GetFormat()));
b1e6bef86955 Tentative loop optimization in CopyPixelDataAndComputeDistribution Benjamin Golinvaux <bgo@osimis.io> parents: 1337 diff changeset	354
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	355 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	356 {
1362 104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	357 T* targetAddrPix = targetAddrLine;
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	358 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	359 {
1362 104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	360 CopyPixel(*targetAddrPix, source);
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	361
1362 104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	362 distribution[*targetAddrPix] += 1;
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	363
1362 104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	364 targetAddrPix++;
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	365 source += bpp;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	366 }
1362 104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	367 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	368 targetAddrLine = reinterpret_cast<T*>(targetAddrLineBytes);
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	369 }
1362 104e0b0f2316 Fixed the broken function, and left the unoptimized Benjamin Golinvaux <bgo@osimis.io> parents: 1349 diff changeset	370 #endif
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	371 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	372 }
1487 75ac66d5f4b2 Added timing test for dose loading and histogram computation Benjamin Golinvaux <bgo@osimis.io> parents: 1381 diff changeset	373 TIME_CHECKPOINT("pixels_and_distribution_ok")
1337 b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	374 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	375
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	376 template <typename T>
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	377 void OrthancMultiframeVolumeLoader::ComputeMinMaxWithOutlierRejection(
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	378 const std::map<T, uint64_t>& distribution)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	379 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	380 if (distribution.size() == 0)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	381 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	382 LOG(ERROR) << "ComputeMinMaxWithOutlierRejection -- Volume image empty.";
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	383 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	384 else
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	385 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	386 OrthancStone::ImageBuffer3D& target = volume_->GetPixelData();
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	387
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	388 const uint64_t width = target.GetWidth();
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	389 const uint64_t height = target.GetHeight();
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	390 const uint64_t depth = target.GetDepth();
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	391 const uint64_t voxelCount = width * height * depth;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	392
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	393 // now that we have distribution[pixelValue] == numberOfPixelsWithValue
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	394 // 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	395 // width * height * depth
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	396 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	397 typename std::map<T, uint64_t>::const_iterator it = distribution.begin();
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	398 uint64_t totalCount = 0;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	399 distributionRawMin_ = static_cast<float>(it->first);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	400
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	401 while (it != distribution.end())
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	402 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	403 T pixelValue = it->first;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	404 uint64_t count = it->second;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	405 totalCount += count;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	406 it++;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	407 if (it == distribution.end())
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	408 distributionRawMax_ = static_cast<float>(pixelValue);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	409 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	410 LOG(INFO) << "Volume image. First distribution value = "
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	411 << static_cast<float>(distributionRawMin_)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	412 << " \| Last distribution value = "
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	413 << static_cast<float>(distributionRawMax_);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	414
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	415 if (totalCount != voxelCount)
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 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	418 << totalCount << ") != VoxC (" << voxelCount;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	419 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	420 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	421 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	422
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	423 // 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	424 uint64_t endRejectionCount = static_cast<uint64_t>(
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	425 outliersHalfRejectionRate_ * voxelCount);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	426
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	427 if (endRejectionCount > voxelCount)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	428 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	429 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	430 << " endRejectionCount = " << endRejectionCount
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	431 << " \| voxelCount = " << voxelCount;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	432 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	433 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	434
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	435 // 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	436 // rejection
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	437 T resultMin = 0;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	438
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	439 // 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	440 // endRejectionCount voxels rejected
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 typename std::map<T, uint64_t>::const_iterator it = distribution.begin();
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 uint64_t currentCount = 0;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	445
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	446 while (it != distribution.end())
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	447 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	448 T pixelValue = it->first;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	449 uint64_t count = it->second;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	450
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	451 // 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	452 // and exit the loop
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	453 if ((currentCount <= endRejectionCount) &&
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	454 (currentCount + count > endRejectionCount))
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	455 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	456 resultMin = pixelValue;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	457 break;
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 else
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 currentCount += count;
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 // and continue walking along the distribution
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	464 it++;
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 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	467
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	468 // 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	469 // rejection
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	470 T resultMax = 0;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	471 // 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	472 // endRejectionCount voxels rejected
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 typename std::map<T, uint64_t>::const_reverse_iterator it = distribution.rbegin();
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 uint64_t currentCount = 0;
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 while (it != distribution.rend())
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	479 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	480 T pixelValue = it->first;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	481 uint64_t count = it->second;
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 if ((currentCount <= endRejectionCount) &&
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	484 (currentCount + count > endRejectionCount))
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	485 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	486 resultMax = pixelValue;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	487 break;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	488 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	489 else
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 currentCount += count;
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 // and continue walking along the distribution
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	494 it++;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	495 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	496 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	497 if (resultMin > resultMax)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	498 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	499 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	500 "resultMin (" << resultMin << ") > resultMax (" << resultMax << ")";
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	501 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	502 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	503 computedDistributionMin_ = static_cast<float>(resultMin);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	504 computedDistributionMax_ = static_cast<float>(resultMax);
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 }
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 template <typename T>
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	509 void OrthancMultiframeVolumeLoader::CopyPixelDataAndComputeMinMax(
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	510 const std::string& pixelData)
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 std::map<T, uint64_t> distribution;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	513 CopyPixelDataAndComputeDistribution(pixelData, distribution);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	514 ComputeMinMaxWithOutlierRejection(distribution);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	515 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	516
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	517 void OrthancMultiframeVolumeLoader::SetUncompressedPixelData(const std::string& pixelData)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	518 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	519 switch (volume_->GetPixelData().GetFormat())
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 case Orthanc::PixelFormat_Grayscale32:
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	522 CopyPixelDataAndComputeMinMax<uint32_t>(pixelData);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	523 break;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	524 case Orthanc::PixelFormat_Grayscale16:
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	525 CopyPixelDataAndComputeMinMax<uint16_t>(pixelData);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	526 break;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	527 case Orthanc::PixelFormat_SignedGrayscale16:
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	528 CopyPixelDataAndComputeMinMax<int16_t>(pixelData);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	529 break;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	530 default:
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	531 throw Orthanc::OrthancException(Orthanc::ErrorCode_NotImplemented);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	532 }
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 volume_->IncrementRevision();
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 pixelDataLoaded_ = true;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	537 BroadcastMessage(OrthancStone::DicomVolumeImage::ContentUpdatedMessage(*volume_));
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	538 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	539
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	540 bool OrthancMultiframeVolumeLoader::HasGeometry() const
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	541 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	542 return volume_->HasGeometry();
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	543 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	544
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	545 const OrthancStone::VolumeImageGeometry& OrthancMultiframeVolumeLoader::GetImageGeometry() const
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	546 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	547 return volume_->GetGeometry();
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
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	550 OrthancMultiframeVolumeLoader::OrthancMultiframeVolumeLoader(
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	551 OrthancStone::ILoadersContext& loadersContext,
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	552 boost::shared_ptr<OrthancStone::DicomVolumeImage> volume,
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	553 float outliersHalfRejectionRate)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	554 : LoaderStateMachine(loadersContext)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	555 , volume_(volume)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	556 , pixelDataLoaded_(false)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	557 , outliersHalfRejectionRate_(outliersHalfRejectionRate)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	558 , distributionRawMin_(0)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	559 , distributionRawMax_(0)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	560 , computedDistributionMin_(0)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	561 , computedDistributionMax_(0)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	562 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	563 if (volume.get() == NULL)
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 throw Orthanc::OrthancException(Orthanc::ErrorCode_NullPointer);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	566 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	567 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	568
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	569
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	570 boost::shared_ptr<OrthancMultiframeVolumeLoader>
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	571 OrthancMultiframeVolumeLoader::Create(
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	572 OrthancStone::ILoadersContext& loadersContext,
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	573 boost::shared_ptr<OrthancStone::DicomVolumeImage> volume,
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	574 float outliersHalfRejectionRate /= 0.0005/)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	575 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	576 boost::shared_ptr<OrthancMultiframeVolumeLoader> obj(
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	577 new OrthancMultiframeVolumeLoader(
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	578 loadersContext,
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	579 volume,
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	580 outliersHalfRejectionRate));
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	581 obj->LoaderStateMachine::PostConstructor();
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	582 return obj;
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 OrthancMultiframeVolumeLoader::~OrthancMultiframeVolumeLoader()
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	586 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	587 LOG(TRACE) << "OrthancMultiframeVolumeLoader::~OrthancMultiframeVolumeLoader()";
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	588 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	589
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	590 void OrthancMultiframeVolumeLoader::GetDistributionMinMax
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	591 (float& minValue, float& maxValue) const
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	592 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	593 if (distributionRawMin_ == 0 && distributionRawMax_ == 0)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	594 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	595 LOG(WARNING) << "GetDistributionMinMaxWithOutliersRejection called before computation!";
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	596 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	597 minValue = distributionRawMin_;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	598 maxValue = distributionRawMax_;
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 void OrthancMultiframeVolumeLoader::GetDistributionMinMaxWithOutliersRejection
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	602 (float& minValue, float& maxValue) const
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	603 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	604 if (computedDistributionMin_ == 0 && computedDistributionMax_ == 0)
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 LOG(WARNING) << "GetDistributionMinMaxWithOutliersRejection called before computation!";
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	607 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	608 minValue = computedDistributionMin_;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	609 maxValue = computedDistributionMax_;
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
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	612 void OrthancMultiframeVolumeLoader::LoadInstance(const std::string& instanceId)
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	613 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	614 Start();
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 instanceId_ = instanceId;
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	617
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	618 {
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	619 std::unique_ptr<OrthancStone::OrthancRestApiCommand> command(new OrthancStone::OrthancRestApiCommand);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	620 command->SetHttpHeader("Accept-Encoding", "gzip");
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	621 command->SetUri("/instances/" + instanceId + "/tags");
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	622 command->AcquirePayload(new LoadGeometry(*this));
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	623 Schedule(command.release());
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	624 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	625
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 std::unique_ptr<OrthancStone::OrthancRestApiCommand> command(new OrthancStone::OrthancRestApiCommand);
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	628 command->SetUri("/instances/" + instanceId + "/metadata/TransferSyntax");
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	629 command->AcquirePayload(new LoadTransferSyntax(*this));
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	630 Schedule(command.release());
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 }
b1396be5aa27 Moved the fixed loaders back from the dead Benjamin Golinvaux <bgo@osimis.io> parents: diff changeset	633 }

Mercurial > hg > orthanc-stone

annotate Framework/Loaders/OrthancMultiframeVolumeLoader.cpp @ 1487:75ac66d5f4b2 bgo-timing-tests