Mercurial > hg > orthanc-stone
annotate Framework/Radiography/RadiographyDicomLayer.cpp @ 1227:a1c0c9c9f9af broker
deprecating SliceImageQuality enum
author | Sebastien Jodogne <s.jodogne@gmail.com> |
---|---|
date | Sun, 08 Dec 2019 12:06:44 +0100 |
parents | b519c1c878f1 |
children | 7ec8fea061b9 |
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 | |
5 * Copyright (C) 2017-2018 Osimis S.A., Belgium | |
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 | |
1066
b537002f83a9
removing broker from deprecated classes
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
739
diff
changeset
|
51 RadiographyDicomLayer::RadiographyDicomLayer(const RadiographyScene& scene) : |
b537002f83a9
removing broker from deprecated classes
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
739
diff
changeset
|
52 RadiographyLayer(scene) |
655
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 } |
1e26bb5f2a02
Fixed truncating conversion warnings + fixed deletion of incomplete type
Benjamin Golinvaux <bgo@osimis.io>
parents:
649
diff
changeset
|
56 |
430 | 57 void RadiographyDicomLayer::SetDicomTags(const OrthancPlugins::FullOrthancDataset& dataset) |
58 { | |
714
d2c0e347ddc2
deprecating DicomFrameConverter
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
655
diff
changeset
|
59 converter_.reset(new Deprecated::DicomFrameConverter); |
430 | 60 converter_->ReadParameters(dataset); |
61 ApplyConverter(); | |
62 | |
63 std::string tmp; | |
64 Vector pixelSpacing; | |
65 | |
66 if (dataset.GetStringValue(tmp, ConvertTag(Orthanc::DICOM_TAG_PIXEL_SPACING)) && | |
67 LinearAlgebra::ParseVector(pixelSpacing, tmp) && | |
68 pixelSpacing.size() == 2) | |
69 { | |
70 SetPixelSpacing(pixelSpacing[0], pixelSpacing[1]); | |
71 } | |
72 | |
73 OrthancPlugins::DicomDatasetReader reader(dataset); | |
74 | |
75 unsigned int width, height; | |
76 if (!reader.GetUnsignedIntegerValue(width, ConvertTag(Orthanc::DICOM_TAG_COLUMNS)) || | |
77 !reader.GetUnsignedIntegerValue(height, ConvertTag(Orthanc::DICOM_TAG_ROWS))) | |
78 { | |
79 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); | |
80 } | |
81 else | |
82 { | |
83 SetSize(width, height); | |
84 } | |
432
4eb96c6b4e96
improved handling of MONOCHROME1, background and invertion
am@osimis.io
parents:
430
diff
changeset
|
85 |
4eb96c6b4e96
improved handling of MONOCHROME1, background and invertion
am@osimis.io
parents:
430
diff
changeset
|
86 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
|
87 { |
4eb96c6b4e96
improved handling of MONOCHROME1, background and invertion
am@osimis.io
parents:
430
diff
changeset
|
88 if (tmp == "MONOCHROME1") |
4eb96c6b4e96
improved handling of MONOCHROME1, background and invertion
am@osimis.io
parents:
430
diff
changeset
|
89 { |
739
be9c1530d40a
deprecating enum SliceImageQuality
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
732
diff
changeset
|
90 SetPreferredPhotomotricDisplayMode(RadiographyPhotometricDisplayMode_Monochrome1); |
432
4eb96c6b4e96
improved handling of MONOCHROME1, background and invertion
am@osimis.io
parents:
430
diff
changeset
|
91 } |
4eb96c6b4e96
improved handling of MONOCHROME1, background and invertion
am@osimis.io
parents:
430
diff
changeset
|
92 else if (tmp == "MONOCHROME2") |
4eb96c6b4e96
improved handling of MONOCHROME1, background and invertion
am@osimis.io
parents:
430
diff
changeset
|
93 { |
739
be9c1530d40a
deprecating enum SliceImageQuality
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
732
diff
changeset
|
94 SetPreferredPhotomotricDisplayMode(RadiographyPhotometricDisplayMode_Monochrome2); |
432
4eb96c6b4e96
improved handling of MONOCHROME1, background and invertion
am@osimis.io
parents:
430
diff
changeset
|
95 } |
4eb96c6b4e96
improved handling of MONOCHROME1, background and invertion
am@osimis.io
parents:
430
diff
changeset
|
96 } |
430 | 97 } |
98 | |
99 void RadiographyDicomLayer::SetSourceImage(Orthanc::ImageAccessor* image) // Takes ownership | |
100 { | |
101 std::auto_ptr<Orthanc::ImageAccessor> raii(image); | |
102 | |
103 if (image == NULL) | |
104 { | |
105 throw Orthanc::OrthancException(Orthanc::ErrorCode_NullPointer); | |
106 } | |
107 | |
108 SetSize(image->GetWidth(), image->GetHeight()); | |
109 | |
110 source_ = raii; | |
111 ApplyConverter(); | |
503
77e0eb83ff63
layers are now Observable and emitting LayerEdited messages
amazy
parents:
432
diff
changeset
|
112 |
623
42dadae61fa9
renamed IObservable::EmitMessage() as BroadcastMessage()
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
503
diff
changeset
|
113 BroadcastMessage(RadiographyLayer::LayerEditedMessage(*this)); |
430 | 114 } |
115 | |
1118
35e798b16b65
RadiographyDicomLayer: allow replacing low res image by high res image
Alain Mazy <alain@mazy.be>
parents:
739
diff
changeset
|
116 void RadiographyDicomLayer::SetSourceImage(Orthanc::ImageAccessor* image, double newPixelSpacingX, double newPixelSpacingY) // Takes ownership |
35e798b16b65
RadiographyDicomLayer: allow replacing low res image by high res image
Alain Mazy <alain@mazy.be>
parents:
739
diff
changeset
|
117 { |
35e798b16b65
RadiographyDicomLayer: allow replacing low res image by high res image
Alain Mazy <alain@mazy.be>
parents:
739
diff
changeset
|
118 std::auto_ptr<Orthanc::ImageAccessor> raii(image); |
35e798b16b65
RadiographyDicomLayer: allow replacing low res image by high res image
Alain Mazy <alain@mazy.be>
parents:
739
diff
changeset
|
119 |
35e798b16b65
RadiographyDicomLayer: allow replacing low res image by high res image
Alain Mazy <alain@mazy.be>
parents:
739
diff
changeset
|
120 if (image == NULL) |
35e798b16b65
RadiographyDicomLayer: allow replacing low res image by high res image
Alain Mazy <alain@mazy.be>
parents:
739
diff
changeset
|
121 { |
35e798b16b65
RadiographyDicomLayer: allow replacing low res image by high res image
Alain Mazy <alain@mazy.be>
parents:
739
diff
changeset
|
122 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
|
123 } |
35e798b16b65
RadiographyDicomLayer: allow replacing low res image by high res image
Alain Mazy <alain@mazy.be>
parents:
739
diff
changeset
|
124 |
35e798b16b65
RadiographyDicomLayer: allow replacing low res image by high res image
Alain Mazy <alain@mazy.be>
parents:
739
diff
changeset
|
125 SetSize(image->GetWidth(), image->GetHeight()); |
35e798b16b65
RadiographyDicomLayer: allow replacing low res image by high res image
Alain Mazy <alain@mazy.be>
parents:
739
diff
changeset
|
126 |
35e798b16b65
RadiographyDicomLayer: allow replacing low res image by high res image
Alain Mazy <alain@mazy.be>
parents:
739
diff
changeset
|
127 source_ = raii; |
35e798b16b65
RadiographyDicomLayer: allow replacing low res image by high res image
Alain Mazy <alain@mazy.be>
parents:
739
diff
changeset
|
128 ApplyConverter(); |
35e798b16b65
RadiographyDicomLayer: allow replacing low res image by high res image
Alain Mazy <alain@mazy.be>
parents:
739
diff
changeset
|
129 |
35e798b16b65
RadiographyDicomLayer: allow replacing low res image by high res image
Alain Mazy <alain@mazy.be>
parents:
739
diff
changeset
|
130 SetPixelSpacing(newPixelSpacingX, newPixelSpacingY, false); |
35e798b16b65
RadiographyDicomLayer: allow replacing low res image by high res image
Alain Mazy <alain@mazy.be>
parents:
739
diff
changeset
|
131 |
35e798b16b65
RadiographyDicomLayer: allow replacing low res image by high res image
Alain Mazy <alain@mazy.be>
parents:
739
diff
changeset
|
132 BroadcastMessage(RadiographyLayer::LayerEditedMessage(*this)); |
35e798b16b65
RadiographyDicomLayer: allow replacing low res image by high res image
Alain Mazy <alain@mazy.be>
parents:
739
diff
changeset
|
133 } |
35e798b16b65
RadiographyDicomLayer: allow replacing low res image by high res image
Alain Mazy <alain@mazy.be>
parents:
739
diff
changeset
|
134 |
649
f1bfe3d1759f
Moved RadiographyDicomLayer::SetDicomFrameConverter to cpp file to prevent
Benjamin Golinvaux <bgo@osimis.io>
parents:
623
diff
changeset
|
135 |
714
d2c0e347ddc2
deprecating DicomFrameConverter
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
655
diff
changeset
|
136 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
|
137 { |
f1bfe3d1759f
Moved RadiographyDicomLayer::SetDicomFrameConverter to cpp file to prevent
Benjamin Golinvaux <bgo@osimis.io>
parents:
623
diff
changeset
|
138 converter_.reset(converter); |
f1bfe3d1759f
Moved RadiographyDicomLayer::SetDicomFrameConverter to cpp file to prevent
Benjamin Golinvaux <bgo@osimis.io>
parents:
623
diff
changeset
|
139 } |
f1bfe3d1759f
Moved RadiographyDicomLayer::SetDicomFrameConverter to cpp file to prevent
Benjamin Golinvaux <bgo@osimis.io>
parents:
623
diff
changeset
|
140 |
430 | 141 void RadiographyDicomLayer::Render(Orthanc::ImageAccessor& buffer, |
142 const AffineTransform2D& viewTransform, | |
1196
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
143 ImageInterpolation interpolation, |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
144 float windowCenter, |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
145 float windowWidth, |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
146 bool applyWindowing) const |
430 | 147 { |
148 if (converted_.get() != NULL) | |
149 { | |
150 if (converted_->GetFormat() != Orthanc::PixelFormat_Float32) | |
151 { | |
152 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError); | |
153 } | |
154 | |
155 unsigned int cropX, cropY, cropWidth, cropHeight; | |
156 GetCrop(cropX, cropY, cropWidth, cropHeight); | |
157 | |
158 AffineTransform2D t = AffineTransform2D::Combine( | |
159 viewTransform, GetTransform(), | |
160 AffineTransform2D::CreateOffset(cropX, cropY)); | |
161 | |
162 Orthanc::ImageAccessor cropped; | |
163 converted_->GetRegion(cropped, cropX, cropY, cropWidth, cropHeight); | |
164 | |
1201 | 165 unsigned int x1, y1, x2, y2; |
166 if (!OrthancStone::GetProjectiveTransformExtent(x1, y1, x2, y2, | |
167 t.GetHomogeneousMatrix(), | |
168 cropped.GetWidth(), | |
169 cropped.GetHeight(), | |
170 buffer.GetWidth(), | |
171 buffer.GetHeight())) | |
172 { | |
173 return; // layer is outside the buffer | |
174 } | |
175 | |
430 | 176 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
|
177 |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
178 if (applyWindowing) |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
179 { |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
180 // 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
|
181 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
|
182 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
|
183 |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
184 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
|
185 { |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
186 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
|
187 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
|
188 { |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
189 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
|
190 |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
191 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
|
192 { |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
193 if (*p >= w1) |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
194 { |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
195 *p = 65535.0; |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
196 } |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
197 else if (*p <= w0) |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
198 { |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
199 *p = 0; |
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 else |
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 // 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
|
204 *p = scaling * (*p - w0); // (*) |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
205 } |
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 } |
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 } |
a5f2a6b04a31
RadiographyScene: windowing is now only applied to the Dicom layer
Alain Mazy <alain@mazy.be>
parents:
1118
diff
changeset
|
210 |
430 | 211 } |
212 } | |
213 | |
214 | |
215 bool RadiographyDicomLayer::GetDefaultWindowing(float& center, | |
216 float& width) const | |
217 { | |
218 if (converter_.get() != NULL && | |
219 converter_->HasDefaultWindow()) | |
220 { | |
221 center = static_cast<float>(converter_->GetDefaultWindowCenter()); | |
222 width = static_cast<float>(converter_->GetDefaultWindowWidth()); | |
223 return true; | |
224 } | |
225 else | |
226 { | |
227 return false; | |
228 } | |
229 } | |
230 | |
231 | |
232 bool RadiographyDicomLayer::GetRange(float& minValue, | |
233 float& maxValue) const | |
234 { | |
235 if (converted_.get() != NULL) | |
236 { | |
237 if (converted_->GetFormat() != Orthanc::PixelFormat_Float32) | |
238 { | |
239 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError); | |
240 } | |
241 | |
242 Orthanc::ImageProcessing::GetMinMaxFloatValue(minValue, maxValue, *converted_); | |
243 return true; | |
244 } | |
245 else | |
246 { | |
247 return false; | |
248 } | |
249 } | |
250 | |
251 } |