orthanc-stone: OrthancStone/Sources/Toolbox/DicomInstanceParameters.cpp annotate

annotate OrthancStone/Sources/Toolbox/DicomInstanceParameters.cpp @ 1739:9ac2a65d4172

upgrade to year 2021

author	Sebastien Jodogne <s.jodogne@gmail.com>
date	Wed, 06 Jan 2021 17:50:00 +0100
parents	5b8b88e5bfd6
children	4bedae500652

rev	line source
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1 /**
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	2 * Stone of Orthanc
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	3 * Copyright (C) 2012-2016 Sebastien Jodogne, Medical Physics
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	4 * Department, University Hospital of Liege, Belgium
1739 9ac2a65d4172 upgrade to year 2021 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1679 diff changeset	5 * Copyright (C) 2017-2021 Osimis S.A., Belgium
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	6 *
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	7 * This program is free software: you can redistribute it and/or
1598 8563ea5d8ae4 relicensing some files, cf. osimis bm26 and chu agreement on 2020-05-20 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1596 diff changeset	8 * modify it under the terms of the GNU Lesser General Public License
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	9 * as published by the Free Software Foundation, either version 3 of
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	10 * the License, or (at your option) any later version.
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	11 *
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	12 * This program is distributed in the hope that it will be useful, but
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	13 * WITHOUT ANY WARRANTY; without even the implied warranty of
1598 8563ea5d8ae4 relicensing some files, cf. osimis bm26 and chu agreement on 2020-05-20 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1596 diff changeset	14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
8563ea5d8ae4 relicensing some files, cf. osimis bm26 and chu agreement on 2020-05-20 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1596 diff changeset	15 * Lesser General Public License for more details.
1596 4fb8fdf03314 removed annoying whitespace Sebastien Jodogne <s.jodogne@gmail.com> parents: 1571 diff changeset	16 *
1598 8563ea5d8ae4 relicensing some files, cf. osimis bm26 and chu agreement on 2020-05-20 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1596 diff changeset	17 * You should have received a copy of the GNU Lesser General Public
8563ea5d8ae4 relicensing some files, cf. osimis bm26 and chu agreement on 2020-05-20 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1596 diff changeset	18 * License along with this program. If not, see
8563ea5d8ae4 relicensing some files, cf. osimis bm26 and chu agreement on 2020-05-20 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1596 diff changeset	19 * <http://www.gnu.org/licenses/>.
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	20 **/
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	21
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	22
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	23 #include "DicomInstanceParameters.h"
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	24
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	25 #include "../Scene2D/ColorTextureSceneLayer.h"
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	26 #include "../Scene2D/FloatTextureSceneLayer.h"
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	27 #include "../Toolbox/GeometryToolbox.h"
1287 8e82fdc6200e Heavy (temporary) logging in the path that leads Benjamin Golinvaux <bgo@osimis.io> parents: 1270 diff changeset	28 #include "../Toolbox/ImageToolbox.h"
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	29
1455 30deba7bc8e2 simplifying include_directories Sebastien Jodogne <s.jodogne@gmail.com> parents: 1353 diff changeset	30 #include <Images/Image.h>
30deba7bc8e2 simplifying include_directories Sebastien Jodogne <s.jodogne@gmail.com> parents: 1353 diff changeset	31 #include <Images/ImageProcessing.h>
30deba7bc8e2 simplifying include_directories Sebastien Jodogne <s.jodogne@gmail.com> parents: 1353 diff changeset	32 #include <Logging.h>
30deba7bc8e2 simplifying include_directories Sebastien Jodogne <s.jodogne@gmail.com> parents: 1353 diff changeset	33 #include <OrthancException.h>
30deba7bc8e2 simplifying include_directories Sebastien Jodogne <s.jodogne@gmail.com> parents: 1353 diff changeset	34 #include <Toolbox.h>
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	35
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	36
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	37 namespace OrthancStone
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	38 {
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	39 static void ExtractFrameOffsets(Vector& target,
d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	40 const Orthanc::DicomMap& dicom,
d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	41 unsigned int numberOfFrames)
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	42 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	43 // http://dicom.nema.org/medical/Dicom/2016a/output/chtml/part03/sect_C.8.8.3.2.html
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	44
1634 a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	45 std::string increment;
a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	46
a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	47 if (dicom.LookupStringValue(increment, Orthanc::DICOM_TAG_FRAME_INCREMENT_POINTER, false))
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	48 {
1634 a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	49 Orthanc::Toolbox::ToUpperCase(increment);
1635 1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	50
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	51 // This is the "Grid Frame Offset Vector" tag (DICOM_TAG_GRID_FRAME_OFFSET_VECTOR)
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	52 if (increment != "3004,000C")
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	53 {
1634 a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	54 LOG(WARNING) << "Bad value for the FrameIncrementPointer tags in a multiframe image";
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	55 target.resize(0);
1634 a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	56 return;
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	57 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	58 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	59
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	60 if (!LinearAlgebra::ParseVector(target, dicom, Orthanc::DICOM_TAG_GRID_FRAME_OFFSET_VECTOR) \|\|
d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	61 target.size() != numberOfFrames)
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	62 {
1634 a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	63 LOG(INFO) << "The frame offset information is missing in a multiframe image";
1635 1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	64
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	65 // DO NOT use ".clear()" here, as the "Vector" class doesn't behave like std::vector!
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	66 target.resize(0);
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	67 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	68 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	69
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	70
1642 5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	71 DicomInstanceParameters::Data::Data(const Orthanc::DicomMap& dicom)
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	72 {
994 1f74bc3459ba fix build due to rename in Orthanc::DicomMap Sebastien Jodogne <s.jodogne@gmail.com> parents: 980 diff changeset	73 if (!dicom.LookupStringValue(studyInstanceUid_, Orthanc::DICOM_TAG_STUDY_INSTANCE_UID, false) \|\|
1f74bc3459ba fix build due to rename in Orthanc::DicomMap Sebastien Jodogne <s.jodogne@gmail.com> parents: 980 diff changeset	74 !dicom.LookupStringValue(seriesInstanceUid_, Orthanc::DICOM_TAG_SERIES_INSTANCE_UID, false) \|\|
1f74bc3459ba fix build due to rename in Orthanc::DicomMap Sebastien Jodogne <s.jodogne@gmail.com> parents: 980 diff changeset	75 !dicom.LookupStringValue(sopInstanceUid_, Orthanc::DICOM_TAG_SOP_INSTANCE_UID, false))
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	76 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	77 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat);
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	78 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	79
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	80 std::string s;
1635 1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	81 if (dicom.LookupStringValue(s, Orthanc::DICOM_TAG_SOP_CLASS_UID, false))
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	82 {
1635 1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	83 sopClassUid_ = StringToSopClassUid(s);
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	84 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	85 else
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	86 {
1635 1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	87 sopClassUid_ = SopClassUid_Other;
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	88 }
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	89
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	90 uint32_t n;
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	91 if (dicom.ParseUnsignedInteger32(n, Orthanc::DICOM_TAG_NUMBER_OF_FRAMES))
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	92 {
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	93 numberOfFrames_ = n;
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	94 }
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	95 else
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	96 {
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	97 numberOfFrames_ = 1;
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	98 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	99
1635 1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	100 if (!dicom.HasTag(Orthanc::DICOM_TAG_COLUMNS) \|\|
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	101 !dicom.GetValue(Orthanc::DICOM_TAG_COLUMNS).ParseFirstUnsignedInteger(width_))
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	102 {
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	103 width_ = 0;
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	104 }
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	105
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	106 if (!dicom.HasTag(Orthanc::DICOM_TAG_ROWS) \|\|
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	107 !dicom.GetValue(Orthanc::DICOM_TAG_ROWS).ParseFirstUnsignedInteger(height_))
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	108 {
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	109 height_ = 0;
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	110 }
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	111
1634 a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	112 if (!dicom.ParseDouble(sliceThickness_, Orthanc::DICOM_TAG_SLICE_THICKNESS))
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	113 {
1634 a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	114 sliceThickness_ = 100.0 * std::numeric_limits<double>::epsilon();
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	115 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	116
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	117 GeometryToolbox::GetPixelSpacing(pixelSpacingX_, pixelSpacingY_, dicom);
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	118
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	119 std::string position, orientation;
994 1f74bc3459ba fix build due to rename in Orthanc::DicomMap Sebastien Jodogne <s.jodogne@gmail.com> parents: 980 diff changeset	120 if (dicom.LookupStringValue(position, Orthanc::DICOM_TAG_IMAGE_POSITION_PATIENT, false) &&
1f74bc3459ba fix build due to rename in Orthanc::DicomMap Sebastien Jodogne <s.jodogne@gmail.com> parents: 980 diff changeset	121 dicom.LookupStringValue(orientation, Orthanc::DICOM_TAG_IMAGE_ORIENTATION_PATIENT, false))
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	122 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	123 geometry_ = CoordinateSystem3D(position, orientation);
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	124 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	125
1635 1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	126 // Must be AFTER setting "numberOfFrames_"
1639 5cdc5b98f14d cppcheck Sebastien Jodogne <s.jodogne@gmail.com> parents: 1637 diff changeset	127 if (numberOfFrames_ > 1)
5cdc5b98f14d cppcheck Sebastien Jodogne <s.jodogne@gmail.com> parents: 1637 diff changeset	128 {
5cdc5b98f14d cppcheck Sebastien Jodogne <s.jodogne@gmail.com> parents: 1637 diff changeset	129 ExtractFrameOffsets(frameOffsets_, dicom, numberOfFrames_);
5cdc5b98f14d cppcheck Sebastien Jodogne <s.jodogne@gmail.com> parents: 1637 diff changeset	130 }
5cdc5b98f14d cppcheck Sebastien Jodogne <s.jodogne@gmail.com> parents: 1637 diff changeset	131 else
5cdc5b98f14d cppcheck Sebastien Jodogne <s.jodogne@gmail.com> parents: 1637 diff changeset	132 {
5cdc5b98f14d cppcheck Sebastien Jodogne <s.jodogne@gmail.com> parents: 1637 diff changeset	133 frameOffsets_.resize(0);
5cdc5b98f14d cppcheck Sebastien Jodogne <s.jodogne@gmail.com> parents: 1637 diff changeset	134 }
1634 a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	135
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	136 if (sopClassUid_ == SopClassUid_RTDose)
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	137 {
1091 5a18e6a395bc Added DoseUnit tag retrieval to DicomInstanceParameters Benjamin Golinvaux <bgo@osimis.io> parents: 1027 diff changeset	138 static const Orthanc::DicomTag DICOM_TAG_DOSE_UNITS(0x3004, 0x0002);
5a18e6a395bc Added DoseUnit tag retrieval to DicomInstanceParameters Benjamin Golinvaux <bgo@osimis.io> parents: 1027 diff changeset	139
5a18e6a395bc Added DoseUnit tag retrieval to DicomInstanceParameters Benjamin Golinvaux <bgo@osimis.io> parents: 1027 diff changeset	140 if (!dicom.LookupStringValue(doseUnits_, DICOM_TAG_DOSE_UNITS, false))
5a18e6a395bc Added DoseUnit tag retrieval to DicomInstanceParameters Benjamin Golinvaux <bgo@osimis.io> parents: 1027 diff changeset	141 {
1141 7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	142 LOG(ERROR) << "Tag DoseUnits (0x3004, 0x0002) is missing in " << sopInstanceUid_;
1634 a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	143 doseUnits_.clear();
1091 5a18e6a395bc Added DoseUnit tag retrieval to DicomInstanceParameters Benjamin Golinvaux <bgo@osimis.io> parents: 1027 diff changeset	144 }
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	145 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	146
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	147 if (dicom.ParseDouble(rescaleIntercept_, Orthanc::DICOM_TAG_RESCALE_INTERCEPT) &&
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	148 dicom.ParseDouble(rescaleSlope_, Orthanc::DICOM_TAG_RESCALE_SLOPE))
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	149 {
1141 7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	150 if (sopClassUid_ == SopClassUid_RTDose)
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	151 {
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	152 LOG(INFO) << "DOSE HAS Rescale*: rescaleIntercept_ = " << rescaleIntercept_ << " rescaleSlope_ = " << rescaleSlope_;
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	153 // WE SHOULD NOT TAKE THE RESCALE VALUE INTO ACCOUNT IN THE CASE OF DOSES
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	154 hasRescale_ = false;
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	155 }
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	156 else
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	157 {
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	158 hasRescale_ = true;
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	159 }
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	160
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	161 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	162 else
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	163 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	164 hasRescale_ = false;
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	165 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	166
1141 7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	167 if (dicom.ParseDouble(doseGridScaling_, Orthanc::DICOM_TAG_DOSE_GRID_SCALING))
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	168 {
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	169 if (sopClassUid_ == SopClassUid_RTDose)
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	170 {
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	171 LOG(INFO) << "DOSE HAS DoseGridScaling: doseGridScaling_ = " << doseGridScaling_;
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	172 }
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	173 }
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	174 else
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	175 {
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	176 doseGridScaling_ = 1.0;
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	177 if (sopClassUid_ == SopClassUid_RTDose)
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	178 {
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	179 LOG(ERROR) << "Tag DoseGridScaling (0x3004, 0x000e) is missing in " << sopInstanceUid_ << " doseGridScaling_ will be set to 1.0";
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	180 }
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	181 }
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	182
1677 51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	183 bool ok = false;
51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	184
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	185 if (LinearAlgebra::ParseVector(windowingPresetCenters_, dicom, Orthanc::DICOM_TAG_WINDOW_CENTER) &&
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	186 LinearAlgebra::ParseVector(windowingPresetWidths_, dicom, Orthanc::DICOM_TAG_WINDOW_WIDTH))
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	187 {
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	188 if (windowingPresetCenters_.size() == windowingPresetWidths_.size())
1677 51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	189 {
51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	190 ok = true;
51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	191 }
51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	192 else
51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	193 {
51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	194 LOG(ERROR) << "Mismatch in the number of preset windowing widths/centers, ignoring this";
51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	195 ok = false;
51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	196 }
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	197 }
1677 51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	198
51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	199 if (!ok)
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	200 {
1677 51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	201 // Don't use "Vector::clear()", as it has not the same meaning as "std::vector::clear()"
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	202 windowingPresetCenters_.resize(0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	203 windowingPresetWidths_.resize(0);
1677 51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	204 }
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	205
980 8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	206 // This computes the "IndexInSeries" metadata from Orthanc (check
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	207 // out "Orthanc::ServerIndex::Store()")
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	208 hasIndexInSeries_ = (
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	209 dicom.ParseUnsignedInteger32(indexInSeries_, Orthanc::DICOM_TAG_INSTANCE_NUMBER) \|\|
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	210 dicom.ParseUnsignedInteger32(indexInSeries_, Orthanc::DICOM_TAG_IMAGE_INDEX));
1646 4e14735e98f8 preparing to remove FrameExtent Sebastien Jodogne <s.jodogne@gmail.com> parents: 1642 diff changeset	211
4e14735e98f8 preparing to remove FrameExtent Sebastien Jodogne <s.jodogne@gmail.com> parents: 1642 diff changeset	212 if (!dicom.LookupStringValue(
4e14735e98f8 preparing to remove FrameExtent Sebastien Jodogne <s.jodogne@gmail.com> parents: 1642 diff changeset	213 frameOfReferenceUid_, Orthanc::DICOM_TAG_FRAME_OF_REFERENCE_UID, false))
4e14735e98f8 preparing to remove FrameExtent Sebastien Jodogne <s.jodogne@gmail.com> parents: 1642 diff changeset	214 {
4e14735e98f8 preparing to remove FrameExtent Sebastien Jodogne <s.jodogne@gmail.com> parents: 1642 diff changeset	215 frameOfReferenceUid_.clear();
4e14735e98f8 preparing to remove FrameExtent Sebastien Jodogne <s.jodogne@gmail.com> parents: 1642 diff changeset	216 }
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	217 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	218
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	219
1642 5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	220 const Orthanc::DicomImageInformation& DicomInstanceParameters::GetImageInformation() const
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	221 {
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	222 assert(tags_.get() != NULL);
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	223
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	224 if (imageInformation_.get() == NULL)
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	225 {
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	226 const_cast<DicomInstanceParameters&>(*this).imageInformation_.
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	227 reset(new Orthanc::DicomImageInformation(GetTags()));
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	228
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	229 assert(imageInformation_->GetWidth() == GetWidth());
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	230 assert(imageInformation_->GetHeight() == GetHeight());
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	231 assert(imageInformation_->GetNumberOfFrames() == GetNumberOfFrames());
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	232 }
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	233
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	234 assert(imageInformation_.get() != NULL);
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	235 return *imageInformation_;
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	236 }
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	237
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	238
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	239 CoordinateSystem3D DicomInstanceParameters::GetFrameGeometry(unsigned int frame) const
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	240 {
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	241 if (frame >= data_.numberOfFrames_)
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	242 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	243 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange);
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	244 }
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	245 else if (data_.frameOffsets_.empty())
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	246 {
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	247 return data_.geometry_;
1634 a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	248 }
a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	249 else
a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	250 {
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	251 assert(data_.frameOffsets_.size() == data_.numberOfFrames_);
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	252
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	253 return CoordinateSystem3D(
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	254 data_.geometry_.GetOrigin() + data_.frameOffsets_[frame] * data_.geometry_.GetNormal(),
d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	255 data_.geometry_.GetAxisX(),
d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	256 data_.geometry_.GetAxisY());
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	257 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	258 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	259
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	260
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	261 bool DicomInstanceParameters::IsPlaneWithinSlice(unsigned int frame,
d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	262 const CoordinateSystem3D& plane) const
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	263 {
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	264 if (frame >= data_.numberOfFrames_)
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	265 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	266 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange);
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	267 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	268
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	269 CoordinateSystem3D tmp = data_.geometry_;
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	270
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	271 if (frame != 0)
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	272 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	273 tmp = GetFrameGeometry(frame);
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	274 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	275
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	276 double distance;
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	277
757 f7c236894c1a cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 746 diff changeset	278 return (CoordinateSystem3D::ComputeDistance(distance, tmp, plane) &&
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	279 distance <= data_.sliceThickness_ / 2.0);
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	280 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	281
1637 d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	282
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	283 bool DicomInstanceParameters::IsColor() const
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	284 {
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	285 Orthanc::PhotometricInterpretation photometric =
1642 5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	286 GetImageInformation().GetPhotometricInterpretation();
1637 d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	287
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	288 return (photometric != Orthanc::PhotometricInterpretation_Monochrome1 &&
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	289 photometric != Orthanc::PhotometricInterpretation_Monochrome2);
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	290 }
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	291
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	292
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	293 void DicomInstanceParameters::ApplyRescaleAndDoseScaling(Orthanc::ImageAccessor& image,
d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	294 bool useDouble) const
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	295 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	296 if (image.GetFormat() != Orthanc::PixelFormat_Float32)
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	297 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	298 throw Orthanc::OrthancException(Orthanc::ErrorCode_IncompatibleImageFormat);
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	299 }
1141 7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	300
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	301 double factor = data_.doseGridScaling_;
1141 7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	302 double offset = 0.0;
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	303
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	304 if (data_.hasRescale_)
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	305 {
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	306 factor *= data_.rescaleSlope_;
d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	307 offset = data_.rescaleIntercept_;
1141 7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	308 }
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	309
1634 a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	310 if (!LinearAlgebra::IsNear(factor, 1) \|\|
a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	311 !LinearAlgebra::IsNear(offset, 0))
1141 7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	312 {
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	313 const unsigned int width = image.GetWidth();
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	314 const unsigned int height = image.GetHeight();
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	315
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	316 for (unsigned int y = 0; y < height; y++)
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	317 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	318 float* p = reinterpret_cast<float*>(image.GetRow(y));
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	319
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	320 if (useDouble)
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	321 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	322 // Slower, accurate implementation using double
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	323 for (unsigned int x = 0; x < width; x++, p++)
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	324 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	325 double value = static_cast<double>(*p);
1141 7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	326 p = static_cast<float>(value factor + offset);
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	327 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	328 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	329 else
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	330 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	331 // Fast, approximate implementation using float
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	332 for (unsigned int x = 0; x < width; x++, p++)
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	333 {
1141 7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	334 p = (p) * static_cast<float>(factor) + static_cast<float>(offset);
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	335 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	336 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	337 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	338 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	339 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	340
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	341 double DicomInstanceParameters::GetRescaleIntercept() const
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	342 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	343 if (data_.hasRescale_)
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	344 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	345 return data_.rescaleIntercept_;
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	346 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	347 else
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	348 {
956 a7351ad54960 Made IsContextLost automatically set the flag by checking with the emscripten Benjamin Golinvaux <bgo@osimis.io> parents: 949 diff changeset	349 LOG(ERROR) << "DicomInstanceParameters::GetRescaleIntercept(): !data_.hasRescale_";
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	350 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls);
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	351 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	352 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	353
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	354
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	355 double DicomInstanceParameters::GetRescaleSlope() const
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	356 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	357 if (data_.hasRescale_)
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	358 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	359 return data_.rescaleSlope_;
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	360 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	361 else
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	362 {
956 a7351ad54960 Made IsContextLost automatically set the flag by checking with the emscripten Benjamin Golinvaux <bgo@osimis.io> parents: 949 diff changeset	363 LOG(ERROR) << "DicomInstanceParameters::GetRescaleSlope(): !data_.hasRescale_";
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	364 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls);
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	365 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	366 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	367
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	368
1637 d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	369 Orthanc::PixelFormat DicomInstanceParameters::GetExpectedPixelFormat() const
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	370 {
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	371 if (GetSopClassUid() == SopClassUid_RTDose)
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	372 {
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	373 switch (GetImageInformation().GetBitsStored())
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	374 {
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	375 case 16:
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	376 return Orthanc::PixelFormat_Grayscale16;
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	377
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	378 case 32:
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	379 return Orthanc::PixelFormat_Grayscale32;
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	380
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	381 default:
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	382 return Orthanc::PixelFormat_Grayscale16; // Rough guess
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	383 }
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	384 }
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	385 else if (IsColor())
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	386 {
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	387 return Orthanc::PixelFormat_RGB24;
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	388 }
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	389 else if (GetImageInformation().IsSigned())
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	390 {
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	391 return Orthanc::PixelFormat_SignedGrayscale16;
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	392 }
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	393 else
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	394 {
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	395 return Orthanc::PixelFormat_Grayscale16; // Rough guess
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	396 }
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	397 }
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	398
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	399
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	400 size_t DicomInstanceParameters::GetWindowingPresetsCount() const
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	401 {
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	402 assert(data_.windowingPresetCenters_.size() == data_.windowingPresetWidths_.size());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	403 return data_.windowingPresetCenters_.size();
1677 51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	404 }
51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	405
51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	406
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	407 float DicomInstanceParameters::GetWindowingPresetCenter(size_t i) const
1677 51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	408 {
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	409 if (i < GetWindowingPresetsCount())
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	410 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	411 return static_cast<float>(data_.windowingPresetCenters_[i]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	412 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	413 else
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	414 {
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	415 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	416 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	417 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	418
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	419
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	420 float DicomInstanceParameters::GetWindowingPresetWidth(size_t i) const
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	421 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	422 if (i < GetWindowingPresetsCount())
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	423 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	424 return static_cast<float>(data_.windowingPresetWidths_[i]);
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	425 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	426 else
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	427 {
1677 51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	428 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange);
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	429 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	430 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	431
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	432
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	433 static void GetWindowingBounds(float& low,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	434 float& high,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	435 double center, // in
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	436 double width) // in
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	437 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	438 low = static_cast<float>(center - width / 2.0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	439 high = static_cast<float>(center + width / 2.0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	440 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	441
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	442
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	443 void DicomInstanceParameters::GetWindowingPresetsUnion(float& center,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	444 float& width) const
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	445 {
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	446 assert(tags_.get() != NULL);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	447 size_t s = GetWindowingPresetsCount();
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	448
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	449 if (s > 0)
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	450 {
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	451 // Use the largest windowing given all the preset windowings
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	452 // that are available in the DICOM tags
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	453 float low, high;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	454 GetWindowingBounds(low, high, GetWindowingPresetCenter(0), GetWindowingPresetWidth(0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	455
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	456 for (size_t i = 1; i < s; i++)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	457 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	458 float a, b;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	459 GetWindowingBounds(a, b, GetWindowingPresetCenter(i), GetWindowingPresetWidth(i));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	460 low = std::min(low, a);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	461 high = std::max(high, b);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	462 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	463
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	464 assert(low <= high);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	465
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	466 if (LinearAlgebra::IsNear(low, high))
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	467 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	468 // Cannot infer a suitable windowing from the available tags
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	469 center = 128.0f;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	470 width = 256.0f;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	471 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	472 else
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	473 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	474 center = (low + high) / 2.0f;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	475 width = (high - low);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	476 }
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	477 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	478 else
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	479 {
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	480 float a, b;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	481 if (tags_->ParseFloat(a, Orthanc::DICOM_TAG_SMALLEST_IMAGE_PIXEL_VALUE) &&
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	482 tags_->ParseFloat(b, Orthanc::DICOM_TAG_LARGEST_IMAGE_PIXEL_VALUE) &&
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	483 a < b)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	484 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	485 center = (a + b) / 2.0f;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	486 width = (b - a);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	487 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	488 else
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	489 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	490 // Cannot infer a suitable windowing from the available tags
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	491 center = 128.0f;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	492 width = 256.0f;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	493 }
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	494 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	495 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	496
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	497
768 55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	498 Orthanc::ImageAccessor* DicomInstanceParameters::ConvertToFloat(const Orthanc::ImageAccessor& pixelData) const
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	499 {
1298 8a0a62189f46 replacing std::auto_ptr by std::unique_ptr Sebastien Jodogne <s.jodogne@gmail.com> parents: 1291 diff changeset	500 std::unique_ptr<Orthanc::Image> converted(new Orthanc::Image(Orthanc::PixelFormat_Float32,
1634 a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	501 pixelData.GetWidth(),
a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	502 pixelData.GetHeight(),
a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	503 false));
768 55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	504 Orthanc::ImageProcessing::Convert(*converted, pixelData);
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	505
1287 8e82fdc6200e Heavy (temporary) logging in the path that leads Benjamin Golinvaux <bgo@osimis.io> parents: 1270 diff changeset	506
768 55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	507 // Correct rescale slope/intercept if need be
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	508 //ApplyRescaleAndDoseScaling(*converted, (pixelData.GetFormat() == Orthanc::PixelFormat_Grayscale32));
d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	509 ApplyRescaleAndDoseScaling(*converted, false);
768 55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	510
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	511 return converted.release();
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	512 }
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	513
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	514
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	515 TextureBaseSceneLayer* DicomInstanceParameters::CreateTexture
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	516 (const Orthanc::ImageAccessor& pixelData) const
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	517 {
1353 af65bce18951 Dummy commented debug block Benjamin Golinvaux <bgo@osimis.io> parents: 1300 diff changeset	518 // {
af65bce18951 Dummy commented debug block Benjamin Golinvaux <bgo@osimis.io> parents: 1300 diff changeset	519 // const Orthanc::ImageAccessor& source = pixelData;
af65bce18951 Dummy commented debug block Benjamin Golinvaux <bgo@osimis.io> parents: 1300 diff changeset	520 // const void* sourceBuffer = source.GetConstBuffer();
af65bce18951 Dummy commented debug block Benjamin Golinvaux <bgo@osimis.io> parents: 1300 diff changeset	521 // intptr_t sourceBufferInt = reinterpret_cast<intptr_t>(sourceBuffer);
af65bce18951 Dummy commented debug block Benjamin Golinvaux <bgo@osimis.io> parents: 1300 diff changeset	522 // int sourceWidth = source.GetWidth();
af65bce18951 Dummy commented debug block Benjamin Golinvaux <bgo@osimis.io> parents: 1300 diff changeset	523 // int sourceHeight = source.GetHeight();
af65bce18951 Dummy commented debug block Benjamin Golinvaux <bgo@osimis.io> parents: 1300 diff changeset	524 // int sourcePitch = source.GetPitch();
af65bce18951 Dummy commented debug block Benjamin Golinvaux <bgo@osimis.io> parents: 1300 diff changeset	525
af65bce18951 Dummy commented debug block Benjamin Golinvaux <bgo@osimis.io> parents: 1300 diff changeset	526 // // TODO: turn error into trace below
af65bce18951 Dummy commented debug block Benjamin Golinvaux <bgo@osimis.io> parents: 1300 diff changeset	527 // LOG(ERROR) << "ConvertGrayscaleToFloat \| source:"
af65bce18951 Dummy commented debug block Benjamin Golinvaux <bgo@osimis.io> parents: 1300 diff changeset	528 // << " W = " << sourceWidth << " H = " << sourceHeight
af65bce18951 Dummy commented debug block Benjamin Golinvaux <bgo@osimis.io> parents: 1300 diff changeset	529 // << " P = " << sourcePitch << " B = " << sourceBufferInt
af65bce18951 Dummy commented debug block Benjamin Golinvaux <bgo@osimis.io> parents: 1300 diff changeset	530 // << " B % 4 == " << sourceBufferInt % 4;
af65bce18951 Dummy commented debug block Benjamin Golinvaux <bgo@osimis.io> parents: 1300 diff changeset	531 // }
af65bce18951 Dummy commented debug block Benjamin Golinvaux <bgo@osimis.io> parents: 1300 diff changeset	532
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	533 assert(sizeof(float) == 4);
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	534
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	535 Orthanc::PixelFormat sourceFormat = pixelData.GetFormat();
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	536
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	537 if (sourceFormat != GetExpectedPixelFormat())
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	538 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	539 throw Orthanc::OrthancException(Orthanc::ErrorCode_IncompatibleImageFormat);
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	540 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	541
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	542 if (sourceFormat == Orthanc::PixelFormat_RGB24)
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	543 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	544 // This is the case of a color image. No conversion has to be done.
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	545 return new ColorTextureSceneLayer(pixelData);
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	546 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	547 else
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	548 {
768 55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	549 // This is the case of a grayscale frame. Convert it to Float32.
1298 8a0a62189f46 replacing std::auto_ptr by std::unique_ptr Sebastien Jodogne <s.jodogne@gmail.com> parents: 1291 diff changeset	550 std::unique_ptr<FloatTextureSceneLayer> texture;
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	551
768 55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	552 if (pixelData.GetFormat() == Orthanc::PixelFormat_Float32)
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	553 {
768 55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	554 texture.reset(new FloatTextureSceneLayer(pixelData));
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	555 }
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	556 else
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	557 {
1298 8a0a62189f46 replacing std::auto_ptr by std::unique_ptr Sebastien Jodogne <s.jodogne@gmail.com> parents: 1291 diff changeset	558 std::unique_ptr<Orthanc::ImageAccessor> converted(ConvertToFloat(pixelData));
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	559 texture.reset(new FloatTextureSceneLayer(*converted));
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	560 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	561
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	562 if (GetWindowingPresetsCount() > 0)
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	563 {
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	564 texture->SetCustomWindowing(GetWindowingPresetCenter(0), GetWindowingPresetWidth(0));
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	565 }
1206 c93a6218f0cd Fixed photometric interpretation in texture creation Benjamin Golinvaux <bgo@osimis.io> parents: 1141 diff changeset	566
1642 5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	567 switch (GetImageInformation().GetPhotometricInterpretation())
1206 c93a6218f0cd Fixed photometric interpretation in texture creation Benjamin Golinvaux <bgo@osimis.io> parents: 1141 diff changeset	568 {
1634 a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	569 case Orthanc::PhotometricInterpretation_Monochrome1:
a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	570 texture->SetInverted(true);
a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	571 break;
a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	572
a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	573 case Orthanc::PhotometricInterpretation_Monochrome2:
a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	574 texture->SetInverted(false);
a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	575 break;
a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	576
a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	577 default:
a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	578 break;
1206 c93a6218f0cd Fixed photometric interpretation in texture creation Benjamin Golinvaux <bgo@osimis.io> parents: 1141 diff changeset	579 }
c93a6218f0cd Fixed photometric interpretation in texture creation Benjamin Golinvaux <bgo@osimis.io> parents: 1141 diff changeset	580
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	581 return texture.release();
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	582 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	583 }
768 55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	584
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	585
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	586 LookupTableTextureSceneLayer* DicomInstanceParameters::CreateLookupTableTexture
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	587 (const Orthanc::ImageAccessor& pixelData) const
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	588 {
1298 8a0a62189f46 replacing std::auto_ptr by std::unique_ptr Sebastien Jodogne <s.jodogne@gmail.com> parents: 1291 diff changeset	589 std::unique_ptr<FloatTextureSceneLayer> texture;
768 55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	590
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	591 if (pixelData.GetFormat() == Orthanc::PixelFormat_Float32)
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	592 {
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	593 return new LookupTableTextureSceneLayer(pixelData);
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	594 }
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	595 else
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	596 {
1298 8a0a62189f46 replacing std::auto_ptr by std::unique_ptr Sebastien Jodogne <s.jodogne@gmail.com> parents: 1291 diff changeset	597 std::unique_ptr<Orthanc::ImageAccessor> converted(ConvertToFloat(pixelData));
768 55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	598 return new LookupTableTextureSceneLayer(*converted);
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	599 }
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	600 }
980 8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	601
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	602
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	603 unsigned int DicomInstanceParameters::GetIndexInSeries() const
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	604 {
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	605 if (data_.hasIndexInSeries_)
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	606 {
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	607 return data_.indexInSeries_;
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	608 }
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	609 else
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	610 {
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	611 LOG(ERROR) << "DicomInstanceParameters::GetIndexInSeries(): !data_.hasIndexInSeries_";
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	612 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls);
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	613 }
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	614 }
1159 acb399643945 DicomInstanceParameters::ApplyRescale() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1141 diff changeset	615
acb399643945 DicomInstanceParameters::ApplyRescale() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1141 diff changeset	616
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	617 double DicomInstanceParameters::ApplyRescale(double value) const
1159 acb399643945 DicomInstanceParameters::ApplyRescale() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1141 diff changeset	618 {
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	619 double factor = data_.doseGridScaling_;
1159 acb399643945 DicomInstanceParameters::ApplyRescale() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1141 diff changeset	620 double offset = 0.0;
acb399643945 DicomInstanceParameters::ApplyRescale() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1141 diff changeset	621
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	622 if (data_.hasRescale_)
1159 acb399643945 DicomInstanceParameters::ApplyRescale() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1141 diff changeset	623 {
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	624 factor *= data_.rescaleSlope_;
d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	625 offset = data_.rescaleIntercept_;
1159 acb399643945 DicomInstanceParameters::ApplyRescale() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1141 diff changeset	626 }
acb399643945 DicomInstanceParameters::ApplyRescale() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1141 diff changeset	627
acb399643945 DicomInstanceParameters::ApplyRescale() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1141 diff changeset	628 return (value * factor + offset);
acb399643945 DicomInstanceParameters::ApplyRescale() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1141 diff changeset	629 }
1161 19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	630
19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	631
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	632 bool DicomInstanceParameters::ComputeRegularSpacing(double& spacing) const
1161 19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	633 {
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	634 if (data_.frameOffsets_.size() == 0) // Not a RT-DOSE
1161 19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	635 {
19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	636 return false;
19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	637 }
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	638 else if (data_.frameOffsets_.size() == 1)
1161 19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	639 {
19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	640 spacing = 1; // Edge case: RT-DOSE with one single frame
19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	641 return true;
19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	642 }
19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	643 else
19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	644 {
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	645 assert(data_.frameOffsets_.size() == GetNumberOfFrames());
d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	646
d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	647 spacing = std::abs(data_.frameOffsets_[1] - data_.frameOffsets_[0]);
1161 19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	648
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	649 for (size_t i = 1; i + 1 < data_.frameOffsets_.size(); i++)
1161 19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	650 {
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	651 double s = data_.frameOffsets_[i + 1] - data_.frameOffsets_[i];
1161 19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	652 if (!LinearAlgebra::IsNear(spacing, s, 0.001))
19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	653 {
19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	654 return false;
19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	655 }
19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	656 }
19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	657
19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	658 return true;
19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	659 }
19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	660 }
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	661 }

Mercurial > hg > orthanc-stone

annotate OrthancStone/Sources/Toolbox/DicomInstanceParameters.cpp @ 1739:9ac2a65d4172