Mercurial > hg > orthanc-stone
annotate Framework/Radiography/RadiographyDicomLayer.cpp @ 1298:8a0a62189f46
replacing std::auto_ptr by std::unique_ptr
author | Sebastien Jodogne <s.jodogne@gmail.com> |
---|---|
date | Mon, 02 Mar 2020 16:31:30 +0100 |
parents | 1c7ae79c426d |
children | 257f2c9a02ac a72c2c9af49a |
rev | line source |
---|---|
430 | 1 /** |
2 * Stone of Orthanc | |
3 * Copyright (C) 2012-2016 Sebastien Jodogne, Medical Physics | |
4 * Department, University Hospital of Liege, Belgium | |
1280
1c7ae79c426d
fix copyright years
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1272
diff
changeset
|
5 * Copyright (C) 2017-2020 Osimis S.A., Belgium |
430 | 6 * |
7 * This program is free software: you can redistribute it and/or | |
8 * modify it under the terms of the GNU Affero General Public License | |
9 * as published by the Free Software Foundation, either version 3 of | |
10 * the License, or (at your option) any later version. | |
11 * | |
12 * This program is distributed in the hope that it will be useful, but | |
13 * WITHOUT ANY WARRANTY; without even the implied warranty of | |
14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU | |
15 * Affero General Public License for more details. | |
16 * | |
17 * You should have received a copy of the GNU Affero General Public License | |
18 * along with this program. If not, see <http://www.gnu.org/licenses/>. | |
19 **/ | |
20 | |
21 | |
22 #include "RadiographyDicomLayer.h" | |
23 | |
24 #include "RadiographyScene.h" | |
732
c35e98d22764
move Deprecated classes to a separate folder
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
714
diff
changeset
|
25 #include "../Deprecated/Toolbox/DicomFrameConverter.h" |
430 | 26 |
27 #include <Core/OrthancException.h> | |
28 #include <Core/Images/Image.h> | |
29 #include <Core/Images/ImageProcessing.h> | |
30 #include <Plugins/Samples/Common/DicomDatasetReader.h> | |
1196
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
31 #include "../Toolbox/ImageGeometry.h" |
430 | 32 |
33 static OrthancPlugins::DicomTag ConvertTag(const Orthanc::DicomTag& tag) | |
34 { | |
35 return OrthancPlugins::DicomTag(tag.GetGroup(), tag.GetElement()); | |
36 } | |
37 | |
38 namespace OrthancStone | |
39 { | |
40 | |
41 void RadiographyDicomLayer::ApplyConverter() | |
42 { | |
43 if (source_.get() != NULL && | |
44 converter_.get() != NULL) | |
45 { | |
46 converted_.reset(converter_->ConvertFrame(*source_)); | |
47 } | |
48 } | |
49 | |
50 | |
655
1e26bb5f2a02
Fixed truncating conversion warnings + fixed deletion of incomplete type
Benjamin Golinvaux <bgo@osimis.io>
parents:
649
diff
changeset
|
51 RadiographyDicomLayer::RadiographyDicomLayer(MessageBroker& broker, const RadiographyScene& scene) : RadiographyLayer(broker, scene) |
1e26bb5f2a02
Fixed truncating conversion warnings + fixed deletion of incomplete type
Benjamin Golinvaux <bgo@osimis.io>
parents:
649
diff
changeset
|
52 { |
1e26bb5f2a02
Fixed truncating conversion warnings + fixed deletion of incomplete type
Benjamin Golinvaux <bgo@osimis.io>
parents:
649
diff
changeset
|
53 |
1e26bb5f2a02
Fixed truncating conversion warnings + fixed deletion of incomplete type
Benjamin Golinvaux <bgo@osimis.io>
parents:
649
diff
changeset
|
54 } |
1e26bb5f2a02
Fixed truncating conversion warnings + fixed deletion of incomplete type
Benjamin Golinvaux <bgo@osimis.io>
parents:
649
diff
changeset
|
55 |
430 | 56 void RadiographyDicomLayer::SetDicomTags(const OrthancPlugins::FullOrthancDataset& dataset) |
57 { | |
714
d2c0e347ddc2
deprecating DicomFrameConverter
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
655
diff
changeset
|
58 converter_.reset(new Deprecated::DicomFrameConverter); |
430 | 59 converter_->ReadParameters(dataset); |
60 ApplyConverter(); | |
61 | |
62 std::string tmp; | |
63 Vector pixelSpacing; | |
64 | |
65 if (dataset.GetStringValue(tmp, ConvertTag(Orthanc::DICOM_TAG_PIXEL_SPACING)) && | |
66 LinearAlgebra::ParseVector(pixelSpacing, tmp) && | |
67 pixelSpacing.size() == 2) | |
68 { | |
69 SetPixelSpacing(pixelSpacing[0], pixelSpacing[1]); | |
70 } | |
71 | |
72 OrthancPlugins::DicomDatasetReader reader(dataset); | |
73 | |
74 unsigned int width, height; | |
75 if (!reader.GetUnsignedIntegerValue(width, ConvertTag(Orthanc::DICOM_TAG_COLUMNS)) || | |
76 !reader.GetUnsignedIntegerValue(height, ConvertTag(Orthanc::DICOM_TAG_ROWS))) | |
77 { | |
78 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); | |
79 } | |
80 else | |
81 { | |
82 SetSize(width, height); | |
83 } | |
432
4eb96c6b4e96
improved handling of MONOCHROME1, background and invertion
am@osimis.io
parents:
430
diff
changeset
|
84 |
4eb96c6b4e96
improved handling of MONOCHROME1, background and invertion
am@osimis.io
parents:
430
diff
changeset
|
85 if (dataset.GetStringValue(tmp, ConvertTag(Orthanc::DICOM_TAG_PHOTOMETRIC_INTERPRETATION))) |
4eb96c6b4e96
improved handling of MONOCHROME1, background and invertion
am@osimis.io
parents:
430
diff
changeset
|
86 { |
4eb96c6b4e96
improved handling of MONOCHROME1, background and invertion
am@osimis.io
parents:
430
diff
changeset
|
87 if (tmp == "MONOCHROME1") |
4eb96c6b4e96
improved handling of MONOCHROME1, background and invertion
am@osimis.io
parents:
430
diff
changeset
|
88 { |
739
be9c1530d40a
deprecating enum SliceImageQuality
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
732
diff
changeset
|
89 SetPreferredPhotomotricDisplayMode(RadiographyPhotometricDisplayMode_Monochrome1); |
432
4eb96c6b4e96
improved handling of MONOCHROME1, background and invertion
am@osimis.io
parents:
430
diff
changeset
|
90 } |
4eb96c6b4e96
improved handling of MONOCHROME1, background and invertion
am@osimis.io
parents:
430
diff
changeset
|
91 else if (tmp == "MONOCHROME2") |
4eb96c6b4e96
improved handling of MONOCHROME1, background and invertion
am@osimis.io
parents:
430
diff
changeset
|
92 { |
739
be9c1530d40a
deprecating enum SliceImageQuality
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
732
diff
changeset
|
93 SetPreferredPhotomotricDisplayMode(RadiographyPhotometricDisplayMode_Monochrome2); |
432
4eb96c6b4e96
improved handling of MONOCHROME1, background and invertion
am@osimis.io
parents:
430
diff
changeset
|
94 } |
4eb96c6b4e96
improved handling of MONOCHROME1, background and invertion
am@osimis.io
parents:
430
diff
changeset
|
95 } |
430 | 96 } |
97 | |
98 void RadiographyDicomLayer::SetSourceImage(Orthanc::ImageAccessor* image) // Takes ownership | |
99 { | |
1298
8a0a62189f46
replacing std::auto_ptr by std::unique_ptr
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1280
diff
changeset
|
100 std::unique_ptr<Orthanc::ImageAccessor> raii(image); |
430 | 101 |
102 if (image == NULL) | |
103 { | |
104 throw Orthanc::OrthancException(Orthanc::ErrorCode_NullPointer); | |
105 } | |
106 | |
107 SetSize(image->GetWidth(), image->GetHeight()); | |
108 | |
1298
8a0a62189f46
replacing std::auto_ptr by std::unique_ptr
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1280
diff
changeset
|
109 #if __cplusplus < 201103L |
8a0a62189f46
replacing std::auto_ptr by std::unique_ptr
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1280
diff
changeset
|
110 source_.reset(raii.release()); |
8a0a62189f46
replacing std::auto_ptr by std::unique_ptr
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1280
diff
changeset
|
111 #else |
8a0a62189f46
replacing std::auto_ptr by std::unique_ptr
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1280
diff
changeset
|
112 source_ = std::move(raii); |
8a0a62189f46
replacing std::auto_ptr by std::unique_ptr
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1280
diff
changeset
|
113 #endif |
8a0a62189f46
replacing std::auto_ptr by std::unique_ptr
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1280
diff
changeset
|
114 |
430 | 115 ApplyConverter(); |
503
77e0eb83ff63
layers are now Observable and emitting LayerEdited messages
amazy
parents:
432
diff
changeset
|
116 |
623
42dadae61fa9
renamed IObservable::EmitMessage() as BroadcastMessage()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
503
diff
changeset
|
117 BroadcastMessage(RadiographyLayer::LayerEditedMessage(*this)); |
430 | 118 } |
119 | |
1272
a989c7d46b9a
options to avoid multiple LayerEditedMessage
Alain Mazy <alain@mazy.be>
parents:
1201
diff
changeset
|
120 void RadiographyDicomLayer::SetSourceImage(Orthanc::ImageAccessor* image, double newPixelSpacingX, double newPixelSpacingY, bool emitLayerEditedEvent) // Takes ownership |
1118
35e798b16b65
RadiographyDicomLayer: allow replacing low res image by high res image
Alain Mazy <alain@mazy.be>
parents:
739
diff
changeset
|
121 { |
1298
8a0a62189f46
replacing std::auto_ptr by std::unique_ptr
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1280
diff
changeset
|
122 std::unique_ptr<Orthanc::ImageAccessor> raii(image); |
1118
35e798b16b65
RadiographyDicomLayer: allow replacing low res image by high res image
Alain Mazy <alain@mazy.be>
parents:
739
diff
changeset
|
123 |
35e798b16b65
RadiographyDicomLayer: allow replacing low res image by high res image
Alain Mazy <alain@mazy.be>
parents:
739
diff
changeset
|
124 if (image == NULL) |
35e798b16b65
RadiographyDicomLayer: allow replacing low res image by high res image
Alain Mazy <alain@mazy.be>
parents:
739
diff
changeset
|
125 { |
35e798b16b65
RadiographyDicomLayer: allow replacing low res image by high res image
Alain Mazy <alain@mazy.be>
parents:
739
diff
changeset
|
126 throw Orthanc::OrthancException(Orthanc::ErrorCode_NullPointer); |
35e798b16b65
RadiographyDicomLayer: allow replacing low res image by high res image
Alain Mazy <alain@mazy.be>
parents:
739
diff
changeset
|
127 } |
35e798b16b65
RadiographyDicomLayer: allow replacing low res image by high res image
Alain Mazy <alain@mazy.be>
parents:
739
diff
changeset
|
128 |
1272
a989c7d46b9a
options to avoid multiple LayerEditedMessage
Alain Mazy <alain@mazy.be>
parents:
1201
diff
changeset
|
129 SetSize(image->GetWidth(), image->GetHeight(), false); |
1118
35e798b16b65
RadiographyDicomLayer: allow replacing low res image by high res image
Alain Mazy <alain@mazy.be>
parents:
739
diff
changeset
|
130 |
1298
8a0a62189f46
replacing std::auto_ptr by std::unique_ptr
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1280
diff
changeset
|
131 #if __cplusplus < 201103L |
8a0a62189f46
replacing std::auto_ptr by std::unique_ptr
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1280
diff
changeset
|
132 source_.reset(raii.release()); |
8a0a62189f46
replacing std::auto_ptr by std::unique_ptr
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1280
diff
changeset
|
133 #else |
8a0a62189f46
replacing std::auto_ptr by std::unique_ptr
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1280
diff
changeset
|
134 source_ = std::move(raii); |
8a0a62189f46
replacing std::auto_ptr by std::unique_ptr
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1280
diff
changeset
|
135 #endif |
8a0a62189f46
replacing std::auto_ptr by std::unique_ptr
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1280
diff
changeset
|
136 |
1118
35e798b16b65
RadiographyDicomLayer: allow replacing low res image by high res image
Alain Mazy <alain@mazy.be>
parents:
739
diff
changeset
|
137 ApplyConverter(); |
35e798b16b65
RadiographyDicomLayer: allow replacing low res image by high res image
Alain Mazy <alain@mazy.be>
parents:
739
diff
changeset
|
138 |
35e798b16b65
RadiographyDicomLayer: allow replacing low res image by high res image
Alain Mazy <alain@mazy.be>
parents:
739
diff
changeset
|
139 SetPixelSpacing(newPixelSpacingX, newPixelSpacingY, false); |
35e798b16b65
RadiographyDicomLayer: allow replacing low res image by high res image
Alain Mazy <alain@mazy.be>
parents:
739
diff
changeset
|
140 |
1272
a989c7d46b9a
options to avoid multiple LayerEditedMessage
Alain Mazy <alain@mazy.be>
parents:
1201
diff
changeset
|
141 if (emitLayerEditedEvent) |
a989c7d46b9a
options to avoid multiple LayerEditedMessage
Alain Mazy <alain@mazy.be>
parents:
1201
diff
changeset
|
142 { |
a989c7d46b9a
options to avoid multiple LayerEditedMessage
Alain Mazy <alain@mazy.be>
parents:
1201
diff
changeset
|
143 BroadcastMessage(RadiographyLayer::LayerEditedMessage(*this)); |
a989c7d46b9a
options to avoid multiple LayerEditedMessage
Alain Mazy <alain@mazy.be>
parents:
1201
diff
changeset
|
144 } |
1118
35e798b16b65
RadiographyDicomLayer: allow replacing low res image by high res image
Alain Mazy <alain@mazy.be>
parents:
739
diff
changeset
|
145 } |
35e798b16b65
RadiographyDicomLayer: allow replacing low res image by high res image
Alain Mazy <alain@mazy.be>
parents:
739
diff
changeset
|
146 |
649
f1bfe3d1759f
Moved RadiographyDicomLayer::SetDicomFrameConverter to cpp file to prevent
Benjamin Golinvaux <bgo@osimis.io>
parents:
623
diff
changeset
|
147 |
714
d2c0e347ddc2
deprecating DicomFrameConverter
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
655
diff
changeset
|
148 void RadiographyDicomLayer::SetDicomFrameConverter(Deprecated::DicomFrameConverter* converter) |
649
f1bfe3d1759f
Moved RadiographyDicomLayer::SetDicomFrameConverter to cpp file to prevent
Benjamin Golinvaux <bgo@osimis.io>
parents:
623
diff
changeset
|
149 { |
f1bfe3d1759f
Moved RadiographyDicomLayer::SetDicomFrameConverter to cpp file to prevent
Benjamin Golinvaux <bgo@osimis.io>
parents:
623
diff
changeset
|
150 converter_.reset(converter); |
f1bfe3d1759f
Moved RadiographyDicomLayer::SetDicomFrameConverter to cpp file to prevent
Benjamin Golinvaux <bgo@osimis.io>
parents:
623
diff
changeset
|
151 } |
f1bfe3d1759f
Moved RadiographyDicomLayer::SetDicomFrameConverter to cpp file to prevent
Benjamin Golinvaux <bgo@osimis.io>
parents:
623
diff
changeset
|
152 |
430 | 153 void RadiographyDicomLayer::Render(Orthanc::ImageAccessor& buffer, |
154 const AffineTransform2D& viewTransform, | |
1196
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
155 ImageInterpolation interpolation, |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
156 float windowCenter, |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
157 float windowWidth, |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
158 bool applyWindowing) const |
430 | 159 { |
160 if (converted_.get() != NULL) | |
161 { | |
162 if (converted_->GetFormat() != Orthanc::PixelFormat_Float32) | |
163 { | |
164 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError); | |
165 } | |
166 | |
167 unsigned int cropX, cropY, cropWidth, cropHeight; | |
168 GetCrop(cropX, cropY, cropWidth, cropHeight); | |
169 | |
170 AffineTransform2D t = AffineTransform2D::Combine( | |
171 viewTransform, GetTransform(), | |
172 AffineTransform2D::CreateOffset(cropX, cropY)); | |
173 | |
174 Orthanc::ImageAccessor cropped; | |
175 converted_->GetRegion(cropped, cropX, cropY, cropWidth, cropHeight); | |
176 | |
1201 | 177 unsigned int x1, y1, x2, y2; |
178 if (!OrthancStone::GetProjectiveTransformExtent(x1, y1, x2, y2, | |
179 t.GetHomogeneousMatrix(), | |
180 cropped.GetWidth(), | |
181 cropped.GetHeight(), | |
182 buffer.GetWidth(), | |
183 buffer.GetHeight())) | |
184 { | |
185 return; // layer is outside the buffer | |
186 } | |
187 | |
430 | 188 t.Apply(buffer, cropped, interpolation, false); |
1196
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
189 |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
190 if (applyWindowing) |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
191 { |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
192 // apply windowing but stay in the range [0.0, 65535.0] |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
193 float w0 = windowCenter - windowWidth / 2.0f; |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
194 float w1 = windowCenter + windowWidth / 2.0f; |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
195 |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
196 if (windowWidth >= 0.001f) // Avoid division by zero at (*) |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
197 { |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
198 float scaling = 1.0f / (w1 - w0) * 65535.0f; |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
199 for (unsigned int y = y1; y <= y2; y++) |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
200 { |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
201 float* p = reinterpret_cast<float*>(buffer.GetRow(y)) + x1; |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
202 |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
203 for (unsigned int x = x1; x <= x2; x++, p++) |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
204 { |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
205 if (*p >= w1) |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
206 { |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
207 *p = 65535.0; |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
208 } |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
209 else if (*p <= w0) |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
210 { |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
211 *p = 0; |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
212 } |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
213 else |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
214 { |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
215 // https://en.wikipedia.org/wiki/Linear_interpolation |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
216 *p = scaling * (*p - w0); // (*) |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
217 } |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
218 } |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
219 } |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
220 } |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
221 } |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
222 |
430 | 223 } |
224 } | |
225 | |
226 | |
227 bool RadiographyDicomLayer::GetDefaultWindowing(float& center, | |
228 float& width) const | |
229 { | |
230 if (converter_.get() != NULL && | |
231 converter_->HasDefaultWindow()) | |
232 { | |
233 center = static_cast<float>(converter_->GetDefaultWindowCenter()); | |
234 width = static_cast<float>(converter_->GetDefaultWindowWidth()); | |
235 return true; | |
236 } | |
237 else | |
238 { | |
239 return false; | |
240 } | |
241 } | |
242 | |
243 | |
244 bool RadiographyDicomLayer::GetRange(float& minValue, | |
245 float& maxValue) const | |
246 { | |
247 if (converted_.get() != NULL) | |
248 { | |
249 if (converted_->GetFormat() != Orthanc::PixelFormat_Float32) | |
250 { | |
251 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError); | |
252 } | |
253 | |
254 Orthanc::ImageProcessing::GetMinMaxFloatValue(minValue, maxValue, *converted_); | |
255 return true; | |
256 } | |
257 else | |
258 { | |
259 return false; | |
260 } | |
261 } | |
262 | |
263 } |