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

annotate OrthancStone/Sources/Toolbox/DicomInstanceParameters.cpp @ 2140:04211f6fd304 default tip

todo

author	Sebastien Jodogne <s.jodogne@gmail.com>
date	Mon, 24 Jun 2024 21:01:33 +0200
parents	16c01cc201e7
children

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
2124 16c01cc201e7 updated copyright, as Osimis is not active on Orthanc anymore Sebastien Jodogne <s.jodogne@gmail.com> parents: 2114 diff changeset	5 * Copyright (C) 2017-2023 Osimis S.A., Belgium
2114 c23eef785569 update year to 2024 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2077 diff changeset	6 * Copyright (C) 2021-2024 Sebastien Jodogne, ICTEAM UCLouvain, Belgium
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	7 *
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	8 * This program is free software: you can redistribute it and/or
1598 8563ea5d8ae4 relicensing some files, cf. osimis bm26 and chu agreement on 2020-05-20 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1596 diff changeset	9 * modify it under the terms of the GNU Lesser General Public License
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	10 * as published by the Free Software Foundation, either version 3 of
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	11 * the License, or (at your option) any later version.
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	12 *
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	13 * This program is distributed in the hope that it will be useful, but
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	14 * WITHOUT ANY WARRANTY; without even the implied warranty of
1598 8563ea5d8ae4 relicensing some files, cf. osimis bm26 and chu agreement on 2020-05-20 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1596 diff changeset	15 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
8563ea5d8ae4 relicensing some files, cf. osimis bm26 and chu agreement on 2020-05-20 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1596 diff changeset	16 * Lesser General Public License for more details.
1596 4fb8fdf03314 removed annoying whitespace Sebastien Jodogne <s.jodogne@gmail.com> parents: 1571 diff changeset	17 *
1598 8563ea5d8ae4 relicensing some files, cf. osimis bm26 and chu agreement on 2020-05-20 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1596 diff changeset	18 * You should have received a copy of the GNU Lesser General Public
8563ea5d8ae4 relicensing some files, cf. osimis bm26 and chu agreement on 2020-05-20 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1596 diff changeset	19 * License along with this program. If not, see
8563ea5d8ae4 relicensing some files, cf. osimis bm26 and chu agreement on 2020-05-20 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1596 diff changeset	20 * <http://www.gnu.org/licenses/>.
746 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
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	24 #include "DicomInstanceParameters.h"
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	25
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	26 #include "../Scene2D/ColorTextureSceneLayer.h"
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	27 #include "../Scene2D/FloatTextureSceneLayer.h"
1750 5f74ebe2516b fixed includes Benjamin Golinvaux <bgo@osimis.io> parents: 1746 diff changeset	28 #include "GeometryToolbox.h"
5f74ebe2516b fixed includes Benjamin Golinvaux <bgo@osimis.io> parents: 1746 diff changeset	29 #include "ImageToolbox.h"
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	30
1455 30deba7bc8e2 simplifying include_directories Sebastien Jodogne <s.jodogne@gmail.com> parents: 1353 diff changeset	31 #include <Images/Image.h>
30deba7bc8e2 simplifying include_directories Sebastien Jodogne <s.jodogne@gmail.com> parents: 1353 diff changeset	32 #include <Images/ImageProcessing.h>
30deba7bc8e2 simplifying include_directories Sebastien Jodogne <s.jodogne@gmail.com> parents: 1353 diff changeset	33 #include <Logging.h>
30deba7bc8e2 simplifying include_directories Sebastien Jodogne <s.jodogne@gmail.com> parents: 1353 diff changeset	34 #include <OrthancException.h>
30deba7bc8e2 simplifying include_directories Sebastien Jodogne <s.jodogne@gmail.com> parents: 1353 diff changeset	35 #include <Toolbox.h>
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	36
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	37
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	38 namespace OrthancStone
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	39 {
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	40 static void ExtractFrameOffsets(Vector& target,
d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	41 const Orthanc::DicomMap& dicom,
d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	42 unsigned int numberOfFrames)
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	43 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	44 // 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	45
1634 a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	46 std::string increment;
a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	47
a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	48 if (dicom.LookupStringValue(increment, Orthanc::DICOM_TAG_FRAME_INCREMENT_POINTER, false))
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	49 {
1634 a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	50 Orthanc::Toolbox::ToUpperCase(increment);
1635 1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	51
1746 4bedae500652 Automatically compute SliceThickness from Benjamin Golinvaux <bgo@osimis.io> parents: 1739 diff changeset	52 // We only support volumes where the FrameIncrementPointer (0028,0009) (required) contains
4bedae500652 Automatically compute SliceThickness from Benjamin Golinvaux <bgo@osimis.io> parents: 1739 diff changeset	53 // the "Grid Frame Offset Vector" tag (DICOM_TAG_GRID_FRAME_OFFSET_VECTOR)
1635 1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	54 if (increment != "3004,000C")
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	55 {
1634 a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	56 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	57 target.resize(0);
1634 a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	58 return;
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	59 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	60 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	61
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	62 if (!LinearAlgebra::ParseVector(target, dicom, Orthanc::DICOM_TAG_GRID_FRAME_OFFSET_VECTOR) \|\|
d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	63 target.size() != numberOfFrames)
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	64 {
1746 4bedae500652 Automatically compute SliceThickness from Benjamin Golinvaux <bgo@osimis.io> parents: 1739 diff changeset	65 LOG(ERROR) << "The frame offset information (GridFrameOffsetVector (3004,000C)) is missing in a multiframe image";
1635 1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	66
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	67 // 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	68 target.resize(0);
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	69 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	70 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	71
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	72
1642 5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	73 DicomInstanceParameters::Data::Data(const Orthanc::DicomMap& dicom)
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	74 {
994 1f74bc3459ba fix build due to rename in Orthanc::DicomMap Sebastien Jodogne <s.jodogne@gmail.com> parents: 980 diff changeset	75 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	76 !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	77 !dicom.LookupStringValue(sopInstanceUid_, Orthanc::DICOM_TAG_SOP_INSTANCE_UID, false))
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	78 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	79 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat);
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	80 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	81
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	82 std::string s;
1635 1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	83 if (dicom.LookupStringValue(s, Orthanc::DICOM_TAG_SOP_CLASS_UID, false))
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	84 {
1635 1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	85 sopClassUid_ = StringToSopClassUid(s);
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	86 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	87 else
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	88 {
1635 1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	89 sopClassUid_ = SopClassUid_Other;
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	90 }
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	91
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	92 uint32_t n;
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	93 if (dicom.ParseUnsignedInteger32(n, Orthanc::DICOM_TAG_NUMBER_OF_FRAMES))
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	94 {
1846 2789b4b0e0a8 DicomInstanceParameters::HasNumberOfFrames() and DicomInstanceParameters::GetInstanceNumber() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1822 diff changeset	95 hasNumberOfFrames_ = true;
1635 1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	96 numberOfFrames_ = n;
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	97 }
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	98 else
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	99 {
1846 2789b4b0e0a8 DicomInstanceParameters::HasNumberOfFrames() and DicomInstanceParameters::GetInstanceNumber() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1822 diff changeset	100 hasNumberOfFrames_ = false;
1635 1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	101 numberOfFrames_ = 1;
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	102 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	103
1635 1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	104 if (!dicom.HasTag(Orthanc::DICOM_TAG_COLUMNS) \|\|
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	105 !dicom.GetValue(Orthanc::DICOM_TAG_COLUMNS).ParseFirstUnsignedInteger(width_))
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	106 {
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	107 width_ = 0;
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
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	110 if (!dicom.HasTag(Orthanc::DICOM_TAG_ROWS) \|\|
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	111 !dicom.GetValue(Orthanc::DICOM_TAG_ROWS).ParseFirstUnsignedInteger(height_))
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	112 {
1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	113 height_ = 0;
1772 6c246f862b00 unit test VolumeRendering.Basic Sebastien Jodogne <s.jodogne@gmail.com> parents: 1750 diff changeset	114 }
1635 1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	115
1746 4bedae500652 Automatically compute SliceThickness from Benjamin Golinvaux <bgo@osimis.io> parents: 1739 diff changeset	116
2076 990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	117 if (dicom.ParseDouble(sliceThickness_, Orthanc::DICOM_TAG_SLICE_THICKNESS))
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	118 {
2076 990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	119 hasSliceThickness_ = true;
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	120 }
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	121 else
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	122 {
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	123 hasSliceThickness_ = false;
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	124
1772 6c246f862b00 unit test VolumeRendering.Basic Sebastien Jodogne <s.jodogne@gmail.com> parents: 1750 diff changeset	125 if (numberOfFrames_ > 1)
6c246f862b00 unit test VolumeRendering.Basic Sebastien Jodogne <s.jodogne@gmail.com> parents: 1750 diff changeset	126 {
6c246f862b00 unit test VolumeRendering.Basic Sebastien Jodogne <s.jodogne@gmail.com> parents: 1750 diff changeset	127 LOG(INFO) << "The (non-madatory) slice thickness information is missing in a multiframe image";
6c246f862b00 unit test VolumeRendering.Basic Sebastien Jodogne <s.jodogne@gmail.com> parents: 1750 diff changeset	128 }
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	129 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	130
1821 36430d73e36c introducing measure units in AnnotationsSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 1817 diff changeset	131 hasPixelSpacing_ = GeometryToolbox::GetPixelSpacing(pixelSpacingX_, pixelSpacingY_, dicom);
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	132
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	133 std::string position, orientation;
994 1f74bc3459ba fix build due to rename in Orthanc::DicomMap Sebastien Jodogne <s.jodogne@gmail.com> parents: 980 diff changeset	134 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	135 dicom.LookupStringValue(orientation, Orthanc::DICOM_TAG_IMAGE_ORIENTATION_PATIENT, false))
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	136 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	137 geometry_ = CoordinateSystem3D(position, orientation);
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	138 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	139
1635 1a714e21ea7c start refactoring DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1634 diff changeset	140 // Must be AFTER setting "numberOfFrames_"
1639 5cdc5b98f14d cppcheck Sebastien Jodogne <s.jodogne@gmail.com> parents: 1637 diff changeset	141 if (numberOfFrames_ > 1)
5cdc5b98f14d cppcheck Sebastien Jodogne <s.jodogne@gmail.com> parents: 1637 diff changeset	142 {
5cdc5b98f14d cppcheck Sebastien Jodogne <s.jodogne@gmail.com> parents: 1637 diff changeset	143 ExtractFrameOffsets(frameOffsets_, dicom, numberOfFrames_);
5cdc5b98f14d cppcheck Sebastien Jodogne <s.jodogne@gmail.com> parents: 1637 diff changeset	144 }
1634 a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	145
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	146 if (sopClassUid_ == SopClassUid_RTDose)
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	147 {
1091 5a18e6a395bc Added DoseUnit tag retrieval to DicomInstanceParameters Benjamin Golinvaux <bgo@osimis.io> parents: 1027 diff changeset	148 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	149
5a18e6a395bc Added DoseUnit tag retrieval to DicomInstanceParameters Benjamin Golinvaux <bgo@osimis.io> parents: 1027 diff changeset	150 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	151 {
1141 7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	152 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	153 doseUnits_.clear();
1091 5a18e6a395bc Added DoseUnit tag retrieval to DicomInstanceParameters Benjamin Golinvaux <bgo@osimis.io> parents: 1027 diff changeset	154 }
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	155 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	156
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	157 if (dicom.ParseDouble(rescaleIntercept_, Orthanc::DICOM_TAG_RESCALE_INTERCEPT) &&
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	158 dicom.ParseDouble(rescaleSlope_, Orthanc::DICOM_TAG_RESCALE_SLOPE))
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	159 {
1141 7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	160 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	161 {
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	162 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	163 // 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	164 hasRescale_ = false;
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	165 }
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	166 else
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	167 {
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	168 hasRescale_ = true;
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	169 }
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	170
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	171 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	172 else
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	173 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	174 hasRescale_ = false;
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	175 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	176
1141 7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	177 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	178 {
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	179 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	180 {
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	181 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	182 }
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	183 }
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	184 else
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	185 {
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	186 doseGridScaling_ = 1.0;
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	187 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	188 {
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	189 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	190 }
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	191 }
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	192
1677 51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	193 bool ok = false;
51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	194
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	195 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	196 LinearAlgebra::ParseVector(windowingPresetWidths_, dicom, Orthanc::DICOM_TAG_WINDOW_WIDTH))
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	197 {
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	198 if (windowingPresetCenters_.size() == windowingPresetWidths_.size())
1677 51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	199 {
51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	200 ok = true;
51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	201 }
51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	202 else
51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	203 {
51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	204 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	205 ok = false;
51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	206 }
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	207 }
1677 51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	208
51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	209 if (!ok)
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	210 {
1677 51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	211 // 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	212 windowingPresetCenters_.resize(0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	213 windowingPresetWidths_.resize(0);
1677 51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	214 }
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	215
980 8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	216 // This computes the "IndexInSeries" metadata from Orthanc (check
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	217 // out "Orthanc::ServerIndex::Store()")
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	218 hasIndexInSeries_ = (
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	219 dicom.ParseUnsignedInteger32(indexInSeries_, Orthanc::DICOM_TAG_INSTANCE_NUMBER) \|\|
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	220 dicom.ParseUnsignedInteger32(indexInSeries_, Orthanc::DICOM_TAG_IMAGE_INDEX));
1646 4e14735e98f8 preparing to remove FrameExtent Sebastien Jodogne <s.jodogne@gmail.com> parents: 1642 diff changeset	221
4e14735e98f8 preparing to remove FrameExtent Sebastien Jodogne <s.jodogne@gmail.com> parents: 1642 diff changeset	222 if (!dicom.LookupStringValue(
4e14735e98f8 preparing to remove FrameExtent Sebastien Jodogne <s.jodogne@gmail.com> parents: 1642 diff changeset	223 frameOfReferenceUid_, Orthanc::DICOM_TAG_FRAME_OF_REFERENCE_UID, false))
4e14735e98f8 preparing to remove FrameExtent Sebastien Jodogne <s.jodogne@gmail.com> parents: 1642 diff changeset	224 {
4e14735e98f8 preparing to remove FrameExtent Sebastien Jodogne <s.jodogne@gmail.com> parents: 1642 diff changeset	225 frameOfReferenceUid_.clear();
4e14735e98f8 preparing to remove FrameExtent Sebastien Jodogne <s.jodogne@gmail.com> parents: 1642 diff changeset	226 }
1846 2789b4b0e0a8 DicomInstanceParameters::HasNumberOfFrames() and DicomInstanceParameters::GetInstanceNumber() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1822 diff changeset	227
2789b4b0e0a8 DicomInstanceParameters::HasNumberOfFrames() and DicomInstanceParameters::GetInstanceNumber() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1822 diff changeset	228 if (!dicom.HasTag(Orthanc::DICOM_TAG_INSTANCE_NUMBER) \|\|
2789b4b0e0a8 DicomInstanceParameters::HasNumberOfFrames() and DicomInstanceParameters::GetInstanceNumber() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1822 diff changeset	229 !dicom.GetValue(Orthanc::DICOM_TAG_INSTANCE_NUMBER).ParseInteger32(instanceNumber_))
2789b4b0e0a8 DicomInstanceParameters::HasNumberOfFrames() and DicomInstanceParameters::GetInstanceNumber() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1822 diff changeset	230 {
2789b4b0e0a8 DicomInstanceParameters::HasNumberOfFrames() and DicomInstanceParameters::GetInstanceNumber() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1822 diff changeset	231 instanceNumber_ = 0;
2789b4b0e0a8 DicomInstanceParameters::HasNumberOfFrames() and DicomInstanceParameters::GetInstanceNumber() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1822 diff changeset	232 }
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	233 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	234
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	235
2076 990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	236 double DicomInstanceParameters::GetSliceThickness() const
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	237 {
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	238 if (data_.hasSliceThickness_)
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	239 {
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	240 return data_.sliceThickness_;
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	241 }
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	242 else
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	243 {
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	244 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls);
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	245 }
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	246 }
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	247
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	248
1642 5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	249 const Orthanc::DicomImageInformation& DicomInstanceParameters::GetImageInformation() const
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	250 {
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	251 assert(tags_.get() != NULL);
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	252
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	253 if (imageInformation_.get() == NULL)
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	254 {
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	255 const_cast<DicomInstanceParameters&>(*this).imageInformation_.
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	256 reset(new Orthanc::DicomImageInformation(GetTags()));
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	257
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	258 assert(imageInformation_->GetWidth() == GetWidth());
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	259 assert(imageInformation_->GetHeight() == GetHeight());
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	260 assert(imageInformation_->GetNumberOfFrames() == GetNumberOfFrames());
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	261 }
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	262
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	263 assert(imageInformation_.get() != NULL);
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	264 return *imageInformation_;
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	265 }
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	266
5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	267
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	268 CoordinateSystem3D DicomInstanceParameters::GetFrameGeometry(unsigned int frame) const
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	269 {
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	270 if (frame >= data_.numberOfFrames_)
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	271 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	272 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange);
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	273 }
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	274 else if (data_.frameOffsets_.empty())
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	275 {
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	276 return data_.geometry_;
1634 a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	277 }
a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	278 else
a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	279 {
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	280 assert(data_.frameOffsets_.size() == data_.numberOfFrames_);
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	281
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	282 return CoordinateSystem3D(
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	283 data_.geometry_.GetOrigin() + data_.frameOffsets_[frame] * data_.geometry_.GetNormal(),
d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	284 data_.geometry_.GetAxisX(),
d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	285 data_.geometry_.GetAxisY());
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	286 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	287 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	288
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	289
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	290 bool DicomInstanceParameters::IsPlaneWithinSlice(unsigned int frame,
d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	291 const CoordinateSystem3D& plane) const
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	292 {
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	293 if (frame >= data_.numberOfFrames_)
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	294 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	295 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange);
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	296 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	297
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	298 CoordinateSystem3D tmp = data_.geometry_;
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	299
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	300 if (frame != 0)
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	301 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	302 tmp = GetFrameGeometry(frame);
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	303 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	304
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	305 double distance;
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	306
757 f7c236894c1a cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 746 diff changeset	307 return (CoordinateSystem3D::ComputeDistance(distance, tmp, plane) &&
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	308 distance <= data_.sliceThickness_ / 2.0);
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	309 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	310
1637 d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	311
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	312 bool DicomInstanceParameters::IsColor() const
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	313 {
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	314 Orthanc::PhotometricInterpretation photometric =
1642 5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	315 GetImageInformation().GetPhotometricInterpretation();
1637 d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	316
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	317 return (photometric != Orthanc::PhotometricInterpretation_Monochrome1 &&
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	318 photometric != Orthanc::PhotometricInterpretation_Monochrome2);
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	319 }
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	320
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	321
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	322 void DicomInstanceParameters::ApplyRescaleAndDoseScaling(Orthanc::ImageAccessor& image,
d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	323 bool useDouble) const
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	324 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	325 if (image.GetFormat() != Orthanc::PixelFormat_Float32)
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	326 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	327 throw Orthanc::OrthancException(Orthanc::ErrorCode_IncompatibleImageFormat);
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	328 }
1141 7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	329
1924 f4050908c6bc display of overlays Sebastien Jodogne <s.jodogne@gmail.com> parents: 1919 diff changeset	330 double scaling = data_.doseGridScaling_;
1141 7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	331 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	332
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	333 if (data_.hasRescale_)
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	334 {
1924 f4050908c6bc display of overlays Sebastien Jodogne <s.jodogne@gmail.com> parents: 1919 diff changeset	335 scaling *= data_.rescaleSlope_;
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	336 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	337 }
7681f3943748 Changed handling of DoseGridScaling: before this commit, rescaleSlope was set to Benjamin Golinvaux <bgo@osimis.io> parents: 1091 diff changeset	338
1924 f4050908c6bc display of overlays Sebastien Jodogne <s.jodogne@gmail.com> parents: 1919 diff changeset	339 Orthanc::ImageProcessing::ShiftScale2(image, offset, scaling, false);
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	340 }
1924 f4050908c6bc display of overlays Sebastien Jodogne <s.jodogne@gmail.com> parents: 1919 diff changeset	341
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	342
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	343 double DicomInstanceParameters::GetRescaleIntercept() const
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	344 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	345 if (data_.hasRescale_)
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	346 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	347 return data_.rescaleIntercept_;
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	348 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	349 else
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	350 {
956 a7351ad54960 Made IsContextLost automatically set the flag by checking with the emscripten Benjamin Golinvaux <bgo@osimis.io> parents: 949 diff changeset	351 LOG(ERROR) << "DicomInstanceParameters::GetRescaleIntercept(): !data_.hasRescale_";
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	352 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls);
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
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	356
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	357 double DicomInstanceParameters::GetRescaleSlope() const
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	358 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	359 if (data_.hasRescale_)
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	360 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	361 return data_.rescaleSlope_;
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	362 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	363 else
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	364 {
956 a7351ad54960 Made IsContextLost automatically set the flag by checking with the emscripten Benjamin Golinvaux <bgo@osimis.io> parents: 949 diff changeset	365 LOG(ERROR) << "DicomInstanceParameters::GetRescaleSlope(): !data_.hasRescale_";
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	366 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls);
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	367 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	368 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	369
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	370
1637 d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	371 Orthanc::PixelFormat DicomInstanceParameters::GetExpectedPixelFormat() const
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 if (GetSopClassUid() == SopClassUid_RTDose)
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 switch (GetImageInformation().GetBitsStored())
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	376 {
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	377 case 16:
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	378 return Orthanc::PixelFormat_Grayscale16;
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	379
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	380 case 32:
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	381 return Orthanc::PixelFormat_Grayscale32;
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	382
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	383 default:
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	384 return Orthanc::PixelFormat_Grayscale16; // Rough guess
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	385 }
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 else if (IsColor())
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 return Orthanc::PixelFormat_RGB24;
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 else if (GetImageInformation().IsSigned())
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 return Orthanc::PixelFormat_SignedGrayscale16;
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 else
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 return Orthanc::PixelFormat_Grayscale16; // Rough guess
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 }
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	400
d569effcd433 cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1636 diff changeset	401
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	402 size_t DicomInstanceParameters::GetWindowingPresetsCount() const
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	403 {
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	404 assert(data_.windowingPresetCenters_.size() == data_.windowingPresetWidths_.size());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	405 return data_.windowingPresetCenters_.size();
1677 51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	406 }
51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	407
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 float DicomInstanceParameters::GetWindowingPresetCenter(size_t i) const
1677 51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	410 {
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	411 if (i < GetWindowingPresetsCount())
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 return static_cast<float>(data_.windowingPresetCenters_[i]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	414 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	415 else
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	416 {
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	417 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange);
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
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 float DicomInstanceParameters::GetWindowingPresetWidth(size_t i) const
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 if (i < GetWindowingPresetsCount())
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	425 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	426 return static_cast<float>(data_.windowingPresetWidths_[i]);
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	427 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	428 else
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	429 {
1677 51bab5188a13 start multiple preset windowings Sebastien Jodogne <s.jodogne@gmail.com> parents: 1646 diff changeset	430 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange);
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	431 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	432 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	433
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	434
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	435 static void GetWindowingBounds(float& low,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	436 float& high,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	437 double center, // in
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	438 double width) // in
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	439 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	440 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	441 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	442 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	443
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	444
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	445 void DicomInstanceParameters::GetWindowingPresetsUnion(float& center,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	446 float& width) const
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	447 {
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	448 assert(tags_.get() != NULL);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	449 size_t s = GetWindowingPresetsCount();
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	450
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	451 if (s > 0)
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	452 {
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	453 // 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	454 // that are available in the DICOM tags
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	455 float low, high;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	456 GetWindowingBounds(low, high, GetWindowingPresetCenter(0), GetWindowingPresetWidth(0));
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 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	459 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	460 float a, b;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	461 GetWindowingBounds(a, b, GetWindowingPresetCenter(i), GetWindowingPresetWidth(i));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	462 low = std::min(low, a);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	463 high = std::max(high, b);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	464 }
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 assert(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 if (LinearAlgebra::IsNear(low, high))
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	469 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	470 // 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	471 center = 128.0f;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	472 width = 256.0f;
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 else
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	475 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	476 center = (low + high) / 2.0f;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	477 width = (high - low);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	478 }
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	479 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	480 else
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	481 {
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	482 float a, b;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	483 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	484 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	485 a < b)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	486 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	487 center = (a + b) / 2.0f;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	488 width = (b - a);
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 else
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	491 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	492 // 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	493 center = 128.0f;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	494 width = 256.0f;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	495 }
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	496 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	497 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	498
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	499
768 55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	500 Orthanc::ImageAccessor* DicomInstanceParameters::ConvertToFloat(const Orthanc::ImageAccessor& pixelData) const
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	501 {
1298 8a0a62189f46 replacing std::auto_ptr by std::unique_ptr Sebastien Jodogne <s.jodogne@gmail.com> parents: 1291 diff changeset	502 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	503 pixelData.GetWidth(),
a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	504 pixelData.GetHeight(),
a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	505 false));
768 55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	506 Orthanc::ImageProcessing::Convert(*converted, pixelData);
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	507
1287 8e82fdc6200e Heavy (temporary) logging in the path that leads Benjamin Golinvaux <bgo@osimis.io> parents: 1270 diff changeset	508
768 55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	509 // Correct rescale slope/intercept if need be
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	510 //ApplyRescaleAndDoseScaling(*converted, (pixelData.GetFormat() == Orthanc::PixelFormat_Grayscale32));
d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	511 ApplyRescaleAndDoseScaling(*converted, false);
768 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 return converted.release();
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
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	516
1919 ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	517 TextureBaseSceneLayer* DicomInstanceParameters::CreateTexture(
ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	518 const Orthanc::ImageAccessor& pixelData) const
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	519 {
1353 af65bce18951 Dummy commented debug block Benjamin Golinvaux <bgo@osimis.io> parents: 1300 diff changeset	520 // {
af65bce18951 Dummy commented debug block Benjamin Golinvaux <bgo@osimis.io> parents: 1300 diff changeset	521 // const Orthanc::ImageAccessor& source = pixelData;
af65bce18951 Dummy commented debug block Benjamin Golinvaux <bgo@osimis.io> parents: 1300 diff changeset	522 // const void* sourceBuffer = source.GetConstBuffer();
af65bce18951 Dummy commented debug block Benjamin Golinvaux <bgo@osimis.io> parents: 1300 diff changeset	523 // intptr_t sourceBufferInt = reinterpret_cast<intptr_t>(sourceBuffer);
af65bce18951 Dummy commented debug block Benjamin Golinvaux <bgo@osimis.io> parents: 1300 diff changeset	524 // int sourceWidth = source.GetWidth();
af65bce18951 Dummy commented debug block Benjamin Golinvaux <bgo@osimis.io> parents: 1300 diff changeset	525 // int sourceHeight = source.GetHeight();
af65bce18951 Dummy commented debug block Benjamin Golinvaux <bgo@osimis.io> parents: 1300 diff changeset	526 // int sourcePitch = source.GetPitch();
af65bce18951 Dummy commented debug block Benjamin Golinvaux <bgo@osimis.io> parents: 1300 diff changeset	527
af65bce18951 Dummy commented debug block Benjamin Golinvaux <bgo@osimis.io> parents: 1300 diff changeset	528 // // TODO: turn error into trace below
af65bce18951 Dummy commented debug block Benjamin Golinvaux <bgo@osimis.io> parents: 1300 diff changeset	529 // LOG(ERROR) << "ConvertGrayscaleToFloat \| source:"
af65bce18951 Dummy commented debug block Benjamin Golinvaux <bgo@osimis.io> parents: 1300 diff changeset	530 // << " W = " << sourceWidth << " H = " << sourceHeight
af65bce18951 Dummy commented debug block Benjamin Golinvaux <bgo@osimis.io> parents: 1300 diff changeset	531 // << " P = " << sourcePitch << " B = " << sourceBufferInt
af65bce18951 Dummy commented debug block Benjamin Golinvaux <bgo@osimis.io> parents: 1300 diff changeset	532 // << " B % 4 == " << sourceBufferInt % 4;
af65bce18951 Dummy commented debug block Benjamin Golinvaux <bgo@osimis.io> parents: 1300 diff changeset	533 // }
af65bce18951 Dummy commented debug block Benjamin Golinvaux <bgo@osimis.io> parents: 1300 diff changeset	534
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	535 assert(sizeof(float) == 4);
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	536
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	537 Orthanc::PixelFormat sourceFormat = pixelData.GetFormat();
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	538
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	539 if (sourceFormat != GetExpectedPixelFormat())
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	540 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	541 throw Orthanc::OrthancException(Orthanc::ErrorCode_IncompatibleImageFormat);
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	542 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	543
1911 898774b4e02d DicomInstanceParameters::CreateTexture() now sets pixel spacing Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	544 std::unique_ptr<TextureBaseSceneLayer> texture;
898774b4e02d DicomInstanceParameters::CreateTexture() now sets pixel spacing Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	545
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	546 if (sourceFormat == Orthanc::PixelFormat_RGB24)
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	547 {
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	548 // This is the case of a color image. No conversion has to be done.
1911 898774b4e02d DicomInstanceParameters::CreateTexture() now sets pixel spacing Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	549 texture.reset(new ColorTextureSceneLayer(pixelData));
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	550 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	551 else
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	552 {
768 55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	553 // This is the case of a grayscale frame. Convert it to Float32.
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	554 if (pixelData.GetFormat() == Orthanc::PixelFormat_Float32)
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	555 {
768 55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	556 texture.reset(new FloatTextureSceneLayer(pixelData));
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	557 }
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	558 else
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	559 {
1298 8a0a62189f46 replacing std::auto_ptr by std::unique_ptr Sebastien Jodogne <s.jodogne@gmail.com> parents: 1291 diff changeset	560 std::unique_ptr<Orthanc::ImageAccessor> converted(ConvertToFloat(pixelData));
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	561 texture.reset(new FloatTextureSceneLayer(*converted));
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	562 }
d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	563
1911 898774b4e02d DicomInstanceParameters::CreateTexture() now sets pixel spacing Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	564 FloatTextureSceneLayer& floatTexture = dynamic_cast<FloatTextureSceneLayer&>(*texture);
898774b4e02d DicomInstanceParameters::CreateTexture() now sets pixel spacing Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	565
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: 1677 diff changeset	566 if (GetWindowingPresetsCount() > 0)
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	567 {
1911 898774b4e02d DicomInstanceParameters::CreateTexture() now sets pixel spacing Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	568 floatTexture.SetCustomWindowing(GetWindowingPresetCenter(0), GetWindowingPresetWidth(0));
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	569 }
1206 c93a6218f0cd Fixed photometric interpretation in texture creation Benjamin Golinvaux <bgo@osimis.io> parents: 1141 diff changeset	570
1642 5cc589bfb385 lazy computation of DicomInstanceParameters::GetImageInformation() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1639 diff changeset	571 switch (GetImageInformation().GetPhotometricInterpretation())
1206 c93a6218f0cd Fixed photometric interpretation in texture creation Benjamin Golinvaux <bgo@osimis.io> parents: 1141 diff changeset	572 {
1634 a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	573 case Orthanc::PhotometricInterpretation_Monochrome1:
1911 898774b4e02d DicomInstanceParameters::CreateTexture() now sets pixel spacing Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	574 floatTexture.SetInverted(true);
1634 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 case Orthanc::PhotometricInterpretation_Monochrome2:
1911 898774b4e02d DicomInstanceParameters::CreateTexture() now sets pixel spacing Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	578 floatTexture.SetInverted(false);
1634 a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	579 break;
a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	580
a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	581 default:
a4418a489e86 improving robustness of DicomInstanceParameters Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	582 break;
1206 c93a6218f0cd Fixed photometric interpretation in texture creation Benjamin Golinvaux <bgo@osimis.io> parents: 1141 diff changeset	583 }
1911 898774b4e02d DicomInstanceParameters::CreateTexture() now sets pixel spacing Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	584 }
1206 c93a6218f0cd Fixed photometric interpretation in texture creation Benjamin Golinvaux <bgo@osimis.io> parents: 1141 diff changeset	585
1911 898774b4e02d DicomInstanceParameters::CreateTexture() now sets pixel spacing Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	586 if (HasPixelSpacing())
898774b4e02d DicomInstanceParameters::CreateTexture() now sets pixel spacing Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	587 {
898774b4e02d DicomInstanceParameters::CreateTexture() now sets pixel spacing Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	588 texture->SetPixelSpacing(GetPixelSpacingX(), GetPixelSpacingY());
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	589 }
1911 898774b4e02d DicomInstanceParameters::CreateTexture() now sets pixel spacing Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	590
898774b4e02d DicomInstanceParameters::CreateTexture() now sets pixel spacing Sebastien Jodogne <s.jodogne@gmail.com> parents: 1871 diff changeset	591 return texture.release();
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	592 }
768 55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	593
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	594
1919 ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	595 LookupTableTextureSceneLayer* DicomInstanceParameters::CreateLookupTableTexture(
ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	596 const Orthanc::ImageAccessor& pixelData) const
768 55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	597 {
1919 ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	598 std::unique_ptr<LookupTableTextureSceneLayer> texture;
ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	599
768 55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	600 if (pixelData.GetFormat() == Orthanc::PixelFormat_Float32)
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	601 {
1919 ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	602 texture.reset(new LookupTableTextureSceneLayer(pixelData));
768 55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	603 }
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	604 else
55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	605 {
1298 8a0a62189f46 replacing std::auto_ptr by std::unique_ptr Sebastien Jodogne <s.jodogne@gmail.com> parents: 1291 diff changeset	606 std::unique_ptr<Orthanc::ImageAccessor> converted(ConvertToFloat(pixelData));
1919 ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	607 texture.reset(new LookupTableTextureSceneLayer(*converted));
768 55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	608 }
1919 ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	609
ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	610 if (HasPixelSpacing())
ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	611 {
ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	612 texture->SetPixelSpacing(GetPixelSpacingX(), GetPixelSpacingY());
ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	613 }
ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	614
ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	615 return texture.release();
768 55411e7da2f7 LookupTableTextureSceneLayer Sebastien Jodogne <s.jodogne@gmail.com> parents: 762 diff changeset	616 }
980 8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	617
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	618
1919 ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	619 LookupTableTextureSceneLayer* DicomInstanceParameters::CreateOverlayTexture(
ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	620 int originX,
ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	621 int originY,
ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	622 const Orthanc::ImageAccessor& overlay) const
ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	623 {
ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	624 if (overlay.GetFormat() != Orthanc::PixelFormat_Grayscale8)
ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	625 {
ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	626 throw Orthanc::OrthancException(Orthanc::ErrorCode_IncompatibleImageFormat);
ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	627 }
ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	628
ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	629 std::unique_ptr<LookupTableTextureSceneLayer> texture(CreateLookupTableTexture(overlay));
ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	630
ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	631 texture->SetOrigin(static_cast<double>(originX - 1) * texture->GetPixelSpacingX(),
ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	632 static_cast<double>(originY - 1) * texture->GetPixelSpacingY());
ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	633
1924 f4050908c6bc display of overlays Sebastien Jodogne <s.jodogne@gmail.com> parents: 1919 diff changeset	634 std::vector<uint8_t> lut(4 * 256);
f4050908c6bc display of overlays Sebastien Jodogne <s.jodogne@gmail.com> parents: 1919 diff changeset	635 for (size_t i = 0; i < 256; i++)
f4050908c6bc display of overlays Sebastien Jodogne <s.jodogne@gmail.com> parents: 1919 diff changeset	636 {
f4050908c6bc display of overlays Sebastien Jodogne <s.jodogne@gmail.com> parents: 1919 diff changeset	637 if (i < 127)
f4050908c6bc display of overlays Sebastien Jodogne <s.jodogne@gmail.com> parents: 1919 diff changeset	638 {
f4050908c6bc display of overlays Sebastien Jodogne <s.jodogne@gmail.com> parents: 1919 diff changeset	639 // Black pixels are converted to transparent pixels
f4050908c6bc display of overlays Sebastien Jodogne <s.jodogne@gmail.com> parents: 1919 diff changeset	640 lut[4 * i] = 0;
f4050908c6bc display of overlays Sebastien Jodogne <s.jodogne@gmail.com> parents: 1919 diff changeset	641 lut[4 * i + 1] = 0;
f4050908c6bc display of overlays Sebastien Jodogne <s.jodogne@gmail.com> parents: 1919 diff changeset	642 lut[4 * i + 2] = 0;
f4050908c6bc display of overlays Sebastien Jodogne <s.jodogne@gmail.com> parents: 1919 diff changeset	643 lut[4 * i + 3] = 0; // alpha
f4050908c6bc display of overlays Sebastien Jodogne <s.jodogne@gmail.com> parents: 1919 diff changeset	644 }
f4050908c6bc display of overlays Sebastien Jodogne <s.jodogne@gmail.com> parents: 1919 diff changeset	645 else
f4050908c6bc display of overlays Sebastien Jodogne <s.jodogne@gmail.com> parents: 1919 diff changeset	646 {
f4050908c6bc display of overlays Sebastien Jodogne <s.jodogne@gmail.com> parents: 1919 diff changeset	647 // White pixels are converted to opaque white
f4050908c6bc display of overlays Sebastien Jodogne <s.jodogne@gmail.com> parents: 1919 diff changeset	648 lut[4 * i] = 255;
f4050908c6bc display of overlays Sebastien Jodogne <s.jodogne@gmail.com> parents: 1919 diff changeset	649 lut[4 * i + 1] = 255;
f4050908c6bc display of overlays Sebastien Jodogne <s.jodogne@gmail.com> parents: 1919 diff changeset	650 lut[4 * i + 2] = 255;
f4050908c6bc display of overlays Sebastien Jodogne <s.jodogne@gmail.com> parents: 1919 diff changeset	651 lut[4 * i + 3] = 255; // alpha
f4050908c6bc display of overlays Sebastien Jodogne <s.jodogne@gmail.com> parents: 1919 diff changeset	652 }
f4050908c6bc display of overlays Sebastien Jodogne <s.jodogne@gmail.com> parents: 1919 diff changeset	653 }
f4050908c6bc display of overlays Sebastien Jodogne <s.jodogne@gmail.com> parents: 1919 diff changeset	654
f4050908c6bc display of overlays Sebastien Jodogne <s.jodogne@gmail.com> parents: 1919 diff changeset	655 texture->SetLookupTable(lut);
f4050908c6bc display of overlays Sebastien Jodogne <s.jodogne@gmail.com> parents: 1919 diff changeset	656
1919 ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	657 return texture.release();
ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	658 }
ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	659
ed4831e08961 added DicomInstanceParameters::CreateOverlayTexture() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1911 diff changeset	660
980 8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	661 unsigned int DicomInstanceParameters::GetIndexInSeries() const
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	662 {
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	663 if (data_.hasIndexInSeries_)
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	664 {
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	665 return data_.indexInSeries_;
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	666 }
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	667 else
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	668 {
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	669 LOG(ERROR) << "DicomInstanceParameters::GetIndexInSeries(): !data_.hasIndexInSeries_";
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	670 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls);
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	671 }
8e497a4e3d96 DicomInstanceParameters::GetIndexInSeries() Sebastien Jodogne <s.jodogne@gmail.com> parents: 956 diff changeset	672 }
1159 acb399643945 DicomInstanceParameters::ApplyRescale() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1141 diff changeset	673
acb399643945 DicomInstanceParameters::ApplyRescale() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1141 diff changeset	674
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	675 double DicomInstanceParameters::ApplyRescale(double value) const
1159 acb399643945 DicomInstanceParameters::ApplyRescale() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1141 diff changeset	676 {
1924 f4050908c6bc display of overlays Sebastien Jodogne <s.jodogne@gmail.com> parents: 1919 diff changeset	677 double scaling = data_.doseGridScaling_;
1159 acb399643945 DicomInstanceParameters::ApplyRescale() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1141 diff changeset	678 double offset = 0.0;
acb399643945 DicomInstanceParameters::ApplyRescale() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1141 diff changeset	679
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	680 if (data_.hasRescale_)
1159 acb399643945 DicomInstanceParameters::ApplyRescale() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1141 diff changeset	681 {
1924 f4050908c6bc display of overlays Sebastien Jodogne <s.jodogne@gmail.com> parents: 1919 diff changeset	682 scaling *= data_.rescaleSlope_;
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	683 offset = data_.rescaleIntercept_;
1159 acb399643945 DicomInstanceParameters::ApplyRescale() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1141 diff changeset	684 }
acb399643945 DicomInstanceParameters::ApplyRescale() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1141 diff changeset	685
1924 f4050908c6bc display of overlays Sebastien Jodogne <s.jodogne@gmail.com> parents: 1919 diff changeset	686 return (value * scaling + offset);
1159 acb399643945 DicomInstanceParameters::ApplyRescale() Sebastien Jodogne <s.jodogne@gmail.com> parents: 1141 diff changeset	687 }
1161 19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	688
19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	689
2076 990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	690 bool DicomInstanceParameters::ComputeFrameOffsetsSpacing(double& spacing) const
1161 19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	691 {
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	692 if (data_.frameOffsets_.size() == 0) // Not a RT-DOSE
1161 19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	693 {
19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	694 return false;
19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	695 }
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	696 else if (data_.frameOffsets_.size() == 1)
1161 19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	697 {
19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	698 spacing = 1; // Edge case: RT-DOSE with one single frame
19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	699 return true;
19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	700 }
19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	701 else
19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	702 {
2076 990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	703 static double THRESHOLD = 0.001;
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	704
2076 990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	705 if (data_.frameOffsets_.size() != GetNumberOfFrames())
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	706 {
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	707 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError);
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	708 }
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	709
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	710 double a = data_.frameOffsets_[1] - data_.frameOffsets_[0];
1161 19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	711
1636 d1e0b08b809d cont Sebastien Jodogne <s.jodogne@gmail.com> parents: 1635 diff changeset	712 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	713 {
2076 990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	714 double b = data_.frameOffsets_[i + 1] - data_.frameOffsets_[i];
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	715 if (!LinearAlgebra::IsNear(a, b, THRESHOLD))
1161 19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	716 {
2076 990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	717 LOG(ERROR) << "Unable to extract slice thickness from GridFrameOffsetVector (3004,000C) (reason: varying spacing)";
1161 19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	718 return false;
19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	719 }
19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	720 }
2076 990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	721
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	722 spacing = std::abs(a);
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	723
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	724 if (HasSliceThickness() &&
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	725 !LinearAlgebra::IsNear(spacing, GetSliceThickness(), THRESHOLD))
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	726 {
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	727 LOG(WARNING) << "SliceThickness and GridFrameOffsetVector (3004,000C) do not match";
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	728 }
1161 19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	729
19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	730 return true;
19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	731 }
19b1c8caade4 fix sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1159 diff changeset	732 }
1822 0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	733
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	734
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	735 void DicomInstanceParameters::SetPixelSpacing(double pixelSpacingX,
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	736 double pixelSpacingY)
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	737 {
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	738 data_.hasPixelSpacing_ = true;
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	739 data_.pixelSpacingX_ = pixelSpacingX;
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	740 data_.pixelSpacingY_ = pixelSpacingY;
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	741 }
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	742
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	743
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	744 static const Json::Value* LookupDicomWebSingleValue(const Json::Value& dicomweb,
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	745 const std::string& tag,
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	746 const std::string& vr)
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	747 {
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	748 static const char* const VALUE = "Value";
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	749 static const char* const VR = "vr";
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	750
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	751 if (dicomweb.type() == Json::objectValue &&
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	752 dicomweb.isMember(tag) &&
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	753 dicomweb[tag].type() == Json::objectValue &&
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	754 dicomweb[tag].isMember(VALUE) &&
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	755 dicomweb[tag].isMember(VR) &&
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	756 dicomweb[tag][VR].type() == Json::stringValue &&
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	757 dicomweb[tag][VR].asString() == vr &&
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	758 dicomweb[tag][VALUE].type() == Json::arrayValue &&
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	759 dicomweb[tag][VALUE].size() == 1u)
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	760 {
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	761 return &dicomweb[tag][VALUE][0];
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	762 }
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	763 else
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	764 {
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	765 return NULL;
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	766 }
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	767 }
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	768
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	769
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	770 void DicomInstanceParameters::EnrichUsingDicomWeb(const Json::Value& dicomweb)
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	771 {
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	772 /**
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	773 * Use DICOM tag "SequenceOfUltrasoundRegions" (0018,6011) in
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	774 * order to derive the pixel spacing on ultrasound (US) images
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	775 **/
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	776
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	777 if (!data_.hasPixelSpacing_)
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	778 {
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	779 const Json::Value* region = LookupDicomWebSingleValue(dicomweb, "00186011", "SQ");
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	780 if (region != NULL)
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	781 {
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	782 const Json::Value* physicalUnitsXDirection = LookupDicomWebSingleValue(*region, "00186024", "US");
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	783 const Json::Value* physicalUnitsYDirection = LookupDicomWebSingleValue(*region, "00186026", "US");
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	784 const Json::Value* physicalDeltaX = LookupDicomWebSingleValue(*region, "0018602C", "FD");
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	785 const Json::Value* physicalDeltaY = LookupDicomWebSingleValue(*region, "0018602E", "FD");
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	786
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	787 if (physicalUnitsXDirection != NULL &&
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	788 physicalUnitsYDirection != NULL &&
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	789 physicalDeltaX != NULL &&
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	790 physicalDeltaY != NULL &&
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	791 physicalUnitsXDirection->type() == Json::intValue &&
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	792 physicalUnitsYDirection->type() == Json::intValue &&
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	793 physicalUnitsXDirection->asInt() == 0x0003 && // Centimeters
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	794 physicalUnitsYDirection->asInt() == 0x0003 && // Centimeters
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	795 physicalDeltaX->isNumeric() &&
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	796 physicalDeltaY->isNumeric())
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	797 {
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	798 // Scene coordinates are expressed in millimeters => multiplication by 10
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	799 SetPixelSpacing(10.0 * physicalDeltaX->asDouble(),
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	800 10.0 * physicalDeltaY->asDouble());
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	801 }
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	802 }
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	803 }
0489fe25ce48 support of pixel spacing in ultrasound images from tag SequenceOfUltrasoundRegions Sebastien Jodogne <s.jodogne@gmail.com> parents: 1821 diff changeset	804 }
2076 990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	805
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	806
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	807 CoordinateSystem3D DicomInstanceParameters::GetMultiFrameGeometry() const
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	808 {
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	809 if (data_.frameOffsets_.empty())
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	810 {
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	811 return data_.geometry_;
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	812 }
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	813 else
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	814 {
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	815 assert(data_.frameOffsets_.size() == data_.numberOfFrames_);
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	816
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	817 double lowest = data_.frameOffsets_[0];
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	818 for (size_t i = 1; i < data_.frameOffsets_.size(); i++)
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	819 {
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	820 lowest = std::min(lowest, data_.frameOffsets_[i]);
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	821 }
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	822
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	823 return CoordinateSystem3D(
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	824 data_.geometry_.GetOrigin() + lowest * data_.geometry_.GetNormal(),
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	825 data_.geometry_.GetAxisX(),
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	826 data_.geometry_.GetAxisY());
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	827 }
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	828 }
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	829
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	830
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	831 bool DicomInstanceParameters::IsReversedFrameOffsets() const
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	832 {
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	833 if (data_.frameOffsets_.size() < 2)
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	834 {
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	835 return false;
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	836 }
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	837 else
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	838 {
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	839 return (data_.frameOffsets_[0] > data_.frameOffsets_[1]);
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	840 }
990f396484b1 fix rendering of RT-DOSE with negative GridFrameOffsetVector Sebastien Jodogne <s.jodogne@gmail.com> parents: 2075 diff changeset	841 }
746 d716bfb3e07c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	842 }

Mercurial > hg > orthanc-stone

annotate OrthancStone/Sources/Toolbox/DicomInstanceParameters.cpp @ 2140:04211f6fd304 default tip