Mercurial > hg > orthanc-stone
annotate Framework/Radiography/RadiographyScene.cpp @ 428:751fb354149e am-vsol-upgrade
ability to change the scene of the RadiographyWidget
author | am@osimis.io |
---|---|
date | Wed, 28 Nov 2018 10:44:28 +0100 |
parents | 660fe6f6bf4a |
children | b85f635f1eb5 |
rev | line source |
---|---|
408 | 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 "RadiographyScene.h" | |
23 | |
24 #include "../Toolbox/DicomFrameConverter.h" | |
25 | |
26 #include <Core/Images/Image.h> | |
27 #include <Core/Images/ImageProcessing.h> | |
28 #include <Core/Images/PamReader.h> | |
29 #include <Core/Images/PamWriter.h> | |
30 #include <Core/Images/PngWriter.h> | |
31 #include <Core/OrthancException.h> | |
32 #include <Core/Toolbox.h> | |
33 #include <Plugins/Samples/Common/DicomDatasetReader.h> | |
34 #include <Plugins/Samples/Common/FullOrthancDataset.h> | |
35 | |
412 | 36 #include <boost/math/special_functions/round.hpp> |
37 | |
408 | 38 |
39 namespace OrthancStone | |
40 { | |
41 RadiographyScene::LayerAccessor::LayerAccessor(RadiographyScene& scene, | |
42 size_t index) : | |
43 scene_(scene), | |
44 index_(index) | |
45 { | |
46 Layers::iterator layer = scene.layers_.find(index); | |
47 if (layer == scene.layers_.end()) | |
48 { | |
49 layer_ = NULL; | |
50 } | |
51 else | |
52 { | |
53 assert(layer->second != NULL); | |
54 layer_ = layer->second; | |
55 } | |
56 } | |
57 | |
426 | 58 |
408 | 59 RadiographyScene::LayerAccessor::LayerAccessor(RadiographyScene& scene, |
60 double x, | |
61 double y) : | |
62 scene_(scene), | |
63 index_(0) // Dummy initialization | |
64 { | |
65 if (scene.LookupLayer(index_, x, y)) | |
66 { | |
67 Layers::iterator layer = scene.layers_.find(index_); | |
426 | 68 |
408 | 69 if (layer == scene.layers_.end()) |
70 { | |
71 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError); | |
72 } | |
73 else | |
74 { | |
75 assert(layer->second != NULL); | |
76 layer_ = layer->second; | |
77 } | |
78 } | |
79 else | |
80 { | |
81 layer_ = NULL; | |
82 } | |
83 } | |
84 | |
85 | |
86 RadiographyScene& RadiographyScene::LayerAccessor::GetScene() const | |
87 { | |
88 if (IsValid()) | |
89 { | |
90 return scene_; | |
91 } | |
92 else | |
93 { | |
94 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls); | |
95 } | |
96 } | |
97 | |
98 | |
99 size_t RadiographyScene::LayerAccessor::GetIndex() const | |
100 { | |
101 if (IsValid()) | |
102 { | |
103 return index_; | |
104 } | |
105 else | |
106 { | |
107 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls); | |
108 } | |
109 } | |
110 | |
111 | |
410
6decc0ba9da5
rename RadiographyScene::Layer as RadiographyLayer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
409
diff
changeset
|
112 RadiographyLayer& RadiographyScene::LayerAccessor::GetLayer() const |
408 | 113 { |
114 if (IsValid()) | |
115 { | |
116 return *layer_; | |
117 } | |
118 else | |
119 { | |
120 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls); | |
121 } | |
426 | 122 } |
408 | 123 |
124 | |
125 | |
410
6decc0ba9da5
rename RadiographyScene::Layer as RadiographyLayer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
409
diff
changeset
|
126 class RadiographyScene::AlphaLayer : public RadiographyLayer |
408 | 127 { |
128 private: | |
129 const RadiographyScene& scene_; | |
130 std::auto_ptr<Orthanc::ImageAccessor> alpha_; // Grayscale8 | |
131 bool useWindowing_; | |
132 float foreground_; | |
133 | |
134 public: | |
135 AlphaLayer(const RadiographyScene& scene) : | |
136 scene_(scene), | |
137 useWindowing_(true), | |
138 foreground_(0) | |
139 { | |
140 } | |
141 | |
142 | |
143 void SetForegroundValue(float foreground) | |
144 { | |
145 useWindowing_ = false; | |
146 foreground_ = foreground; | |
147 } | |
426 | 148 |
149 | |
408 | 150 void SetAlpha(Orthanc::ImageAccessor* image) |
151 { | |
152 std::auto_ptr<Orthanc::ImageAccessor> raii(image); | |
426 | 153 |
408 | 154 if (image == NULL) |
155 { | |
156 throw Orthanc::OrthancException(Orthanc::ErrorCode_NullPointer); | |
157 } | |
158 | |
159 if (image->GetFormat() != Orthanc::PixelFormat_Grayscale8) | |
160 { | |
161 throw Orthanc::OrthancException(Orthanc::ErrorCode_IncompatibleImageFormat); | |
162 } | |
163 | |
164 SetSize(image->GetWidth(), image->GetHeight()); | |
165 alpha_ = raii; | |
166 } | |
167 | |
168 | |
169 void LoadText(const Orthanc::Font& font, | |
170 const std::string& utf8) | |
171 { | |
172 SetAlpha(font.RenderAlpha(utf8)); | |
426 | 173 } |
408 | 174 |
175 | |
176 virtual bool GetDefaultWindowing(float& center, | |
177 float& width) const | |
178 { | |
179 return false; | |
180 } | |
426 | 181 |
408 | 182 |
183 virtual void Render(Orthanc::ImageAccessor& buffer, | |
409 | 184 const AffineTransform2D& viewTransform, |
408 | 185 ImageInterpolation interpolation) const |
186 { | |
187 if (alpha_.get() == NULL) | |
188 { | |
189 return; | |
190 } | |
426 | 191 |
408 | 192 if (buffer.GetFormat() != Orthanc::PixelFormat_Float32) |
193 { | |
194 throw Orthanc::OrthancException(Orthanc::ErrorCode_IncompatibleImageFormat); | |
195 } | |
196 | |
197 unsigned int cropX, cropY, cropWidth, cropHeight; | |
198 GetCrop(cropX, cropY, cropWidth, cropHeight); | |
199 | |
409 | 200 const AffineTransform2D t = AffineTransform2D::Combine( |
426 | 201 viewTransform, GetTransform(), |
202 AffineTransform2D::CreateOffset(cropX, cropY)); | |
408 | 203 |
204 Orthanc::ImageAccessor cropped; | |
205 alpha_->GetRegion(cropped, cropX, cropY, cropWidth, cropHeight); | |
426 | 206 |
408 | 207 Orthanc::Image tmp(Orthanc::PixelFormat_Grayscale8, buffer.GetWidth(), buffer.GetHeight(), false); |
409 | 208 |
209 t.Apply(tmp, cropped, interpolation, true /* clear */); | |
408 | 210 |
211 // Blit | |
212 const unsigned int width = buffer.GetWidth(); | |
213 const unsigned int height = buffer.GetHeight(); | |
214 | |
215 float value = foreground_; | |
426 | 216 |
408 | 217 if (useWindowing_) |
218 { | |
219 float center, width; | |
220 if (scene_.GetWindowing(center, width)) | |
221 { | |
222 value = center + width / 2.0f; | |
223 } | |
224 } | |
426 | 225 |
408 | 226 for (unsigned int y = 0; y < height; y++) |
227 { | |
228 float *q = reinterpret_cast<float*>(buffer.GetRow(y)); | |
229 const uint8_t *p = reinterpret_cast<uint8_t*>(tmp.GetRow(y)); | |
230 | |
231 for (unsigned int x = 0; x < width; x++, p++, q++) | |
232 { | |
233 float a = static_cast<float>(*p) / 255.0f; | |
426 | 234 |
408 | 235 *q = (a * value + (1.0f - a) * (*q)); |
236 } | |
426 | 237 } |
408 | 238 } |
239 | |
426 | 240 |
408 | 241 virtual bool GetRange(float& minValue, |
242 float& maxValue) const | |
243 { | |
244 if (useWindowing_) | |
245 { | |
246 return false; | |
247 } | |
248 else | |
249 { | |
250 minValue = 0; | |
251 maxValue = 0; | |
252 | |
253 if (foreground_ < 0) | |
254 { | |
255 minValue = foreground_; | |
256 } | |
257 | |
258 if (foreground_ > 0) | |
259 { | |
260 maxValue = foreground_; | |
261 } | |
262 | |
263 return true; | |
264 } | |
265 } | |
266 }; | |
426 | 267 |
268 | |
408 | 269 |
410
6decc0ba9da5
rename RadiographyScene::Layer as RadiographyLayer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
409
diff
changeset
|
270 class RadiographyScene::DicomLayer : public RadiographyLayer |
408 | 271 { |
272 private: | |
273 std::auto_ptr<Orthanc::ImageAccessor> source_; // Content of PixelData | |
274 std::auto_ptr<DicomFrameConverter> converter_; | |
275 std::auto_ptr<Orthanc::ImageAccessor> converted_; // Float32 | |
276 | |
277 static OrthancPlugins::DicomTag ConvertTag(const Orthanc::DicomTag& tag) | |
278 { | |
279 return OrthancPlugins::DicomTag(tag.GetGroup(), tag.GetElement()); | |
280 } | |
426 | 281 |
408 | 282 |
283 void ApplyConverter() | |
284 { | |
285 if (source_.get() != NULL && | |
286 converter_.get() != NULL) | |
287 { | |
288 converted_.reset(converter_->ConvertFrame(*source_)); | |
289 } | |
290 } | |
426 | 291 |
408 | 292 public: |
293 void SetDicomTags(const OrthancPlugins::FullOrthancDataset& dataset) | |
294 { | |
295 converter_.reset(new DicomFrameConverter); | |
296 converter_->ReadParameters(dataset); | |
297 ApplyConverter(); | |
298 | |
299 std::string tmp; | |
300 Vector pixelSpacing; | |
426 | 301 |
408 | 302 if (dataset.GetStringValue(tmp, ConvertTag(Orthanc::DICOM_TAG_PIXEL_SPACING)) && |
303 LinearAlgebra::ParseVector(pixelSpacing, tmp) && | |
304 pixelSpacing.size() == 2) | |
305 { | |
306 SetPixelSpacing(pixelSpacing[0], pixelSpacing[1]); | |
307 } | |
308 | |
309 //SetPan(-0.5 * GetPixelSpacingX(), -0.5 * GetPixelSpacingY()); | |
310 | |
311 OrthancPlugins::DicomDatasetReader reader(dataset); | |
312 | |
313 unsigned int width, height; | |
314 if (!reader.GetUnsignedIntegerValue(width, ConvertTag(Orthanc::DICOM_TAG_COLUMNS)) || | |
315 !reader.GetUnsignedIntegerValue(height, ConvertTag(Orthanc::DICOM_TAG_ROWS))) | |
316 { | |
317 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); | |
318 } | |
319 else | |
320 { | |
321 SetSize(width, height); | |
322 } | |
323 } | |
324 | |
426 | 325 |
408 | 326 void SetSourceImage(Orthanc::ImageAccessor* image) // Takes ownership |
327 { | |
328 std::auto_ptr<Orthanc::ImageAccessor> raii(image); | |
426 | 329 |
408 | 330 if (image == NULL) |
331 { | |
332 throw Orthanc::OrthancException(Orthanc::ErrorCode_NullPointer); | |
333 } | |
334 | |
335 SetSize(image->GetWidth(), image->GetHeight()); | |
426 | 336 |
408 | 337 source_ = raii; |
338 ApplyConverter(); | |
339 } | |
340 | |
426 | 341 |
408 | 342 virtual void Render(Orthanc::ImageAccessor& buffer, |
409 | 343 const AffineTransform2D& viewTransform, |
408 | 344 ImageInterpolation interpolation) const |
345 { | |
346 if (converted_.get() != NULL) | |
347 { | |
348 if (converted_->GetFormat() != Orthanc::PixelFormat_Float32) | |
349 { | |
350 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError); | |
351 } | |
352 | |
353 unsigned int cropX, cropY, cropWidth, cropHeight; | |
354 GetCrop(cropX, cropY, cropWidth, cropHeight); | |
355 | |
409 | 356 AffineTransform2D t = AffineTransform2D::Combine( |
426 | 357 viewTransform, GetTransform(), |
358 AffineTransform2D::CreateOffset(cropX, cropY)); | |
408 | 359 |
360 Orthanc::ImageAccessor cropped; | |
361 converted_->GetRegion(cropped, cropX, cropY, cropWidth, cropHeight); | |
362 | |
409 | 363 t.Apply(buffer, cropped, interpolation, false); |
408 | 364 } |
365 } | |
366 | |
367 | |
368 virtual bool GetDefaultWindowing(float& center, | |
369 float& width) const | |
370 { | |
371 if (converter_.get() != NULL && | |
372 converter_->HasDefaultWindow()) | |
373 { | |
374 center = static_cast<float>(converter_->GetDefaultWindowCenter()); | |
375 width = static_cast<float>(converter_->GetDefaultWindowWidth()); | |
376 return true; | |
377 } | |
378 else | |
379 { | |
380 return false; | |
381 } | |
382 } | |
383 | |
384 | |
385 virtual bool GetRange(float& minValue, | |
386 float& maxValue) const | |
387 { | |
388 if (converted_.get() != NULL) | |
389 { | |
390 if (converted_->GetFormat() != Orthanc::PixelFormat_Float32) | |
391 { | |
392 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError); | |
393 } | |
394 | |
395 Orthanc::ImageProcessing::GetMinMaxFloatValue(minValue, maxValue, *converted_); | |
396 return true; | |
397 } | |
398 else | |
399 { | |
400 return false; | |
401 } | |
402 } | |
403 }; | |
404 | |
405 | |
410
6decc0ba9da5
rename RadiographyScene::Layer as RadiographyLayer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
409
diff
changeset
|
406 RadiographyLayer& RadiographyScene::RegisterLayer(RadiographyLayer* layer) |
408 | 407 { |
408 if (layer == NULL) | |
409 { | |
410 throw Orthanc::OrthancException(Orthanc::ErrorCode_NullPointer); | |
411 } | |
412 | |
410
6decc0ba9da5
rename RadiographyScene::Layer as RadiographyLayer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
409
diff
changeset
|
413 std::auto_ptr<RadiographyLayer> raii(layer); |
426 | 414 |
428
751fb354149e
ability to change the scene of the RadiographyWidget
am@osimis.io
parents:
426
diff
changeset
|
415 LOG(INFO) << "Registering layer: " << countLayers_; |
751fb354149e
ability to change the scene of the RadiographyWidget
am@osimis.io
parents:
426
diff
changeset
|
416 |
408 | 417 size_t index = countLayers_++; |
418 raii->SetIndex(index); | |
419 layers_[index] = raii.release(); | |
420 | |
421 EmitMessage(GeometryChangedMessage(*this)); | |
422 EmitMessage(ContentChangedMessage(*this)); | |
423 | |
424 return *layer; | |
425 } | |
426 | 426 |
408 | 427 |
417
aee3d7941c9b
preparing to load images using DICOMweb
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
412
diff
changeset
|
428 RadiographyScene::RadiographyScene(MessageBroker& broker) : |
408 | 429 IObserver(broker), |
430 IObservable(broker), | |
431 countLayers_(0), | |
432 hasWindowing_(false), | |
433 windowingCenter_(0), // Dummy initialization | |
434 windowingWidth_(0) // Dummy initialization | |
435 { | |
436 } | |
437 | |
438 | |
439 RadiographyScene::~RadiographyScene() | |
440 { | |
441 for (Layers::iterator it = layers_.begin(); it != layers_.end(); it++) | |
442 { | |
443 assert(it->second != NULL); | |
444 delete it->second; | |
445 } | |
446 } | |
447 | |
425 | 448 void RadiographyScene::RemoveLayer(size_t layerIndex) |
449 { | |
428
751fb354149e
ability to change the scene of the RadiographyWidget
am@osimis.io
parents:
426
diff
changeset
|
450 LOG(INFO) << "Removing layer: " << layerIndex; |
751fb354149e
ability to change the scene of the RadiographyWidget
am@osimis.io
parents:
426
diff
changeset
|
451 |
425 | 452 if (layerIndex > countLayers_) |
453 { | |
454 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange); | |
455 } | |
456 delete layers_[layerIndex]; | |
457 layers_.erase(layerIndex); | |
458 countLayers_--; | |
428
751fb354149e
ability to change the scene of the RadiographyWidget
am@osimis.io
parents:
426
diff
changeset
|
459 LOG(INFO) << "Removing layer, there are now : " << countLayers_ << " layers"; |
425 | 460 } |
461 | |
462 RadiographyLayer& RadiographyScene::GetLayer(size_t layerIndex) | |
463 { | |
464 if (layerIndex > countLayers_) | |
465 { | |
466 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange); | |
467 } | |
468 return *(layers_[layerIndex]); | |
469 } | |
408 | 470 |
471 bool RadiographyScene::GetWindowing(float& center, | |
472 float& width) const | |
473 { | |
474 if (hasWindowing_) | |
475 { | |
476 center = windowingCenter_; | |
477 width = windowingWidth_; | |
478 return true; | |
479 } | |
480 else | |
481 { | |
482 return false; | |
483 } | |
484 } | |
485 | |
486 | |
487 void RadiographyScene::GetWindowingWithDefault(float& center, | |
488 float& width) const | |
489 { | |
490 if (!GetWindowing(center, width)) | |
491 { | |
492 center = 128; | |
493 width = 256; | |
494 } | |
495 } | |
496 | |
497 | |
498 void RadiographyScene::SetWindowing(float center, | |
499 float width) | |
500 { | |
501 hasWindowing_ = true; | |
502 windowingCenter_ = center; | |
503 windowingWidth_ = width; | |
504 } | |
505 | |
506 | |
410
6decc0ba9da5
rename RadiographyScene::Layer as RadiographyLayer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
409
diff
changeset
|
507 RadiographyLayer& RadiographyScene::LoadText(const Orthanc::Font& font, |
6decc0ba9da5
rename RadiographyScene::Layer as RadiographyLayer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
409
diff
changeset
|
508 const std::string& utf8) |
408 | 509 { |
510 std::auto_ptr<AlphaLayer> alpha(new AlphaLayer(*this)); | |
511 alpha->LoadText(font, utf8); | |
512 | |
513 return RegisterLayer(alpha.release()); | |
514 } | |
515 | |
426 | 516 |
410
6decc0ba9da5
rename RadiographyScene::Layer as RadiographyLayer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
409
diff
changeset
|
517 RadiographyLayer& RadiographyScene::LoadTestBlock(unsigned int width, |
6decc0ba9da5
rename RadiographyScene::Layer as RadiographyLayer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
409
diff
changeset
|
518 unsigned int height) |
408 | 519 { |
520 std::auto_ptr<Orthanc::Image> block(new Orthanc::Image(Orthanc::PixelFormat_Grayscale8, width, height, false)); | |
521 | |
522 for (unsigned int padding = 0; | |
523 (width > 2 * padding) && (height > 2 * padding); | |
524 padding++) | |
525 { | |
526 uint8_t color; | |
527 if (255 > 10 * padding) | |
528 { | |
529 color = 255 - 10 * padding; | |
530 } | |
531 else | |
532 { | |
533 color = 0; | |
534 } | |
535 | |
536 Orthanc::ImageAccessor region; | |
537 block->GetRegion(region, padding, padding, width - 2 * padding, height - 2 * padding); | |
538 Orthanc::ImageProcessing::Set(region, color); | |
539 } | |
540 | |
541 std::auto_ptr<AlphaLayer> alpha(new AlphaLayer(*this)); | |
542 alpha->SetAlpha(block.release()); | |
543 | |
544 return RegisterLayer(alpha.release()); | |
545 } | |
546 | |
426 | 547 |
417
aee3d7941c9b
preparing to load images using DICOMweb
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
412
diff
changeset
|
548 RadiographyLayer& RadiographyScene::LoadDicomFrame(OrthancApiClient& orthanc, |
aee3d7941c9b
preparing to load images using DICOMweb
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
412
diff
changeset
|
549 const std::string& instance, |
410
6decc0ba9da5
rename RadiographyScene::Layer as RadiographyLayer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
409
diff
changeset
|
550 unsigned int frame, |
6decc0ba9da5
rename RadiographyScene::Layer as RadiographyLayer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
409
diff
changeset
|
551 bool httpCompression) |
408 | 552 { |
410
6decc0ba9da5
rename RadiographyScene::Layer as RadiographyLayer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
409
diff
changeset
|
553 RadiographyLayer& layer = RegisterLayer(new DicomLayer); |
408 | 554 |
555 { | |
417
aee3d7941c9b
preparing to load images using DICOMweb
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
412
diff
changeset
|
556 IWebService::HttpHeaders headers; |
408 | 557 std::string uri = "/instances/" + instance + "/tags"; |
426 | 558 |
417
aee3d7941c9b
preparing to load images using DICOMweb
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
412
diff
changeset
|
559 orthanc.GetBinaryAsync( |
426 | 560 uri, headers, |
561 new Callable<RadiographyScene, OrthancApiClient::BinaryResponseReadyMessage> | |
562 (*this, &RadiographyScene::OnTagsReceived), NULL, | |
563 new Orthanc::SingleValueObject<size_t>(layer.GetIndex())); | |
408 | 564 } |
565 | |
566 { | |
417
aee3d7941c9b
preparing to load images using DICOMweb
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
412
diff
changeset
|
567 IWebService::HttpHeaders headers; |
408 | 568 headers["Accept"] = "image/x-portable-arbitrarymap"; |
569 | |
570 if (httpCompression) | |
571 { | |
572 headers["Accept-Encoding"] = "gzip"; | |
573 } | |
426 | 574 |
408 | 575 std::string uri = ("/instances/" + instance + "/frames/" + |
576 boost::lexical_cast<std::string>(frame) + "/image-uint16"); | |
426 | 577 |
417
aee3d7941c9b
preparing to load images using DICOMweb
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
412
diff
changeset
|
578 orthanc.GetBinaryAsync( |
426 | 579 uri, headers, |
580 new Callable<RadiographyScene, OrthancApiClient::BinaryResponseReadyMessage> | |
581 (*this, &RadiographyScene::OnFrameReceived), NULL, | |
582 new Orthanc::SingleValueObject<size_t>(layer.GetIndex())); | |
408 | 583 } |
584 | |
585 return layer; | |
586 } | |
587 | |
417
aee3d7941c9b
preparing to load images using DICOMweb
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
412
diff
changeset
|
588 |
aee3d7941c9b
preparing to load images using DICOMweb
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
412
diff
changeset
|
589 RadiographyLayer& RadiographyScene::LoadDicomWebFrame(IWebService& web) |
aee3d7941c9b
preparing to load images using DICOMweb
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
412
diff
changeset
|
590 { |
aee3d7941c9b
preparing to load images using DICOMweb
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
412
diff
changeset
|
591 RadiographyLayer& layer = RegisterLayer(new DicomLayer); |
aee3d7941c9b
preparing to load images using DICOMweb
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
412
diff
changeset
|
592 |
426 | 593 |
417
aee3d7941c9b
preparing to load images using DICOMweb
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
412
diff
changeset
|
594 return layer; |
aee3d7941c9b
preparing to load images using DICOMweb
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
412
diff
changeset
|
595 } |
aee3d7941c9b
preparing to load images using DICOMweb
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
412
diff
changeset
|
596 |
aee3d7941c9b
preparing to load images using DICOMweb
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
412
diff
changeset
|
597 |
426 | 598 |
408 | 599 void RadiographyScene::OnTagsReceived(const OrthancApiClient::BinaryResponseReadyMessage& message) |
600 { | |
601 size_t index = dynamic_cast<const Orthanc::SingleValueObject<size_t>&> | |
426 | 602 (message.GetPayload()).GetValue(); |
408 | 603 |
604 LOG(INFO) << "JSON received: " << message.GetUri().c_str() | |
605 << " (" << message.GetAnswerSize() << " bytes) for layer " << index; | |
426 | 606 |
408 | 607 Layers::iterator layer = layers_.find(index); |
608 if (layer != layers_.end()) | |
609 { | |
610 assert(layer->second != NULL); | |
426 | 611 |
408 | 612 OrthancPlugins::FullOrthancDataset dicom(message.GetAnswer(), message.GetAnswerSize()); |
613 dynamic_cast<DicomLayer*>(layer->second)->SetDicomTags(dicom); | |
614 | |
615 float c, w; | |
616 if (!hasWindowing_ && | |
617 layer->second->GetDefaultWindowing(c, w)) | |
618 { | |
619 hasWindowing_ = true; | |
620 windowingCenter_ = c; | |
621 windowingWidth_ = w; | |
622 } | |
623 | |
624 EmitMessage(GeometryChangedMessage(*this)); | |
625 } | |
626 } | |
426 | 627 |
408 | 628 |
629 void RadiographyScene::OnFrameReceived(const OrthancApiClient::BinaryResponseReadyMessage& message) | |
630 { | |
631 size_t index = dynamic_cast<const Orthanc::SingleValueObject<size_t>&>(message.GetPayload()).GetValue(); | |
426 | 632 |
408 | 633 LOG(INFO) << "DICOM frame received: " << message.GetUri().c_str() |
634 << " (" << message.GetAnswerSize() << " bytes) for layer " << index; | |
426 | 635 |
408 | 636 Layers::iterator layer = layers_.find(index); |
637 if (layer != layers_.end()) | |
638 { | |
639 assert(layer->second != NULL); | |
640 | |
641 std::string content; | |
642 if (message.GetAnswerSize() > 0) | |
643 { | |
644 content.assign(reinterpret_cast<const char*>(message.GetAnswer()), message.GetAnswerSize()); | |
645 } | |
426 | 646 |
408 | 647 std::auto_ptr<Orthanc::PamReader> reader(new Orthanc::PamReader); |
648 reader->ReadFromMemory(content); | |
649 dynamic_cast<DicomLayer*>(layer->second)->SetSourceImage(reader.release()); | |
650 | |
651 EmitMessage(ContentChangedMessage(*this)); | |
652 } | |
653 } | |
654 | |
655 | |
656 Extent2D RadiographyScene::GetSceneExtent() const | |
657 { | |
658 Extent2D extent; | |
659 | |
660 for (Layers::const_iterator it = layers_.begin(); | |
661 it != layers_.end(); ++it) | |
662 { | |
663 assert(it->second != NULL); | |
664 extent.Union(it->second->GetExtent()); | |
665 } | |
666 | |
667 return extent; | |
668 } | |
426 | 669 |
408 | 670 |
671 void RadiographyScene::Render(Orthanc::ImageAccessor& buffer, | |
409 | 672 const AffineTransform2D& viewTransform, |
408 | 673 ImageInterpolation interpolation) const |
674 { | |
675 Orthanc::ImageProcessing::Set(buffer, 0); | |
676 | |
677 // Render layers in the background-to-foreground order | |
678 for (size_t index = 0; index < countLayers_; index++) | |
679 { | |
680 Layers::const_iterator it = layers_.find(index); | |
681 if (it != layers_.end()) | |
682 { | |
683 assert(it->second != NULL); | |
684 it->second->Render(buffer, viewTransform, interpolation); | |
685 } | |
686 } | |
687 } | |
688 | |
689 | |
690 bool RadiographyScene::LookupLayer(size_t& index /* out */, | |
691 double x, | |
692 double y) const | |
693 { | |
694 // Render layers in the foreground-to-background order | |
695 for (size_t i = countLayers_; i > 0; i--) | |
696 { | |
697 index = i - 1; | |
698 Layers::const_iterator it = layers_.find(index); | |
699 if (it != layers_.end()) | |
700 { | |
701 assert(it->second != NULL); | |
702 if (it->second->Contains(x, y)) | |
703 { | |
704 return true; | |
705 } | |
706 } | |
707 } | |
708 | |
709 return false; | |
710 } | |
711 | |
426 | 712 |
408 | 713 void RadiographyScene::DrawBorder(CairoContext& context, |
714 unsigned int layer, | |
715 double zoom) | |
716 { | |
717 Layers::const_iterator found = layers_.find(layer); | |
426 | 718 |
408 | 719 if (found != layers_.end()) |
720 { | |
721 context.SetSourceColor(255, 0, 0); | |
722 found->second->DrawBorders(context, zoom); | |
723 } | |
724 } | |
725 | |
726 | |
727 void RadiographyScene::GetRange(float& minValue, | |
728 float& maxValue) const | |
729 { | |
730 bool first = true; | |
426 | 731 |
408 | 732 for (Layers::const_iterator it = layers_.begin(); |
733 it != layers_.end(); it++) | |
734 { | |
735 assert(it->second != NULL); | |
736 | |
737 float a, b; | |
738 if (it->second->GetRange(a, b)) | |
739 { | |
740 if (first) | |
741 { | |
742 minValue = a; | |
743 maxValue = b; | |
744 first = false; | |
745 } | |
746 else | |
747 { | |
748 minValue = std::min(a, minValue); | |
749 maxValue = std::max(b, maxValue); | |
750 } | |
751 } | |
752 } | |
753 | |
754 if (first) | |
755 { | |
756 minValue = 0; | |
757 maxValue = 0; | |
758 } | |
759 } | |
760 | |
761 | |
417
aee3d7941c9b
preparing to load images using DICOMweb
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
412
diff
changeset
|
762 void RadiographyScene::ExportDicom(OrthancApiClient& orthanc, |
aee3d7941c9b
preparing to load images using DICOMweb
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
412
diff
changeset
|
763 const Orthanc::DicomMap& dicom, |
426 | 764 const std::string& parentOrthancId, |
408 | 765 double pixelSpacingX, |
766 double pixelSpacingY, | |
767 bool invert, | |
768 ImageInterpolation interpolation, | |
769 bool usePam) | |
770 { | |
426 | 771 Json::Value createDicomRequestContent; |
772 | |
773 Export(createDicomRequestContent, dicom, pixelSpacingX, pixelSpacingY, invert, interpolation, usePam); | |
774 | |
775 if (!parentOrthancId.empty()) | |
776 { | |
777 createDicomRequestContent["Parent"] = parentOrthancId; | |
778 } | |
779 | |
780 orthanc.PostJsonAsyncExpectJson( | |
781 "/tools/create-dicom", createDicomRequestContent, | |
782 new Callable<RadiographyScene, OrthancApiClient::JsonResponseReadyMessage> | |
783 (*this, &RadiographyScene::OnDicomExported), | |
784 NULL, NULL); | |
785 } | |
786 | |
787 // Export using PAM is faster than using PNG, but requires Orthanc | |
788 // core >= 1.4.3 | |
789 void RadiographyScene::Export(Json::Value& createDicomRequestContent, | |
790 const Orthanc::DicomMap& dicom, | |
791 double pixelSpacingX, | |
792 double pixelSpacingY, | |
793 bool invert, | |
794 ImageInterpolation interpolation, | |
795 bool usePam) | |
796 { | |
408 | 797 if (pixelSpacingX <= 0 || |
798 pixelSpacingY <= 0) | |
799 { | |
800 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange); | |
801 } | |
426 | 802 |
408 | 803 LOG(INFO) << "Exporting DICOM"; |
804 | |
805 Extent2D extent = GetSceneExtent(); | |
806 | |
807 int w = std::ceil(extent.GetWidth() / pixelSpacingX); | |
808 int h = std::ceil(extent.GetHeight() / pixelSpacingY); | |
809 | |
810 if (w < 0 || h < 0) | |
811 { | |
812 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError); | |
813 } | |
814 | |
815 Orthanc::Image layers(Orthanc::PixelFormat_Float32, | |
816 static_cast<unsigned int>(w), | |
817 static_cast<unsigned int>(h), false); | |
818 | |
409 | 819 AffineTransform2D view = AffineTransform2D::Combine( |
426 | 820 AffineTransform2D::CreateScaling(1.0 / pixelSpacingX, 1.0 / pixelSpacingY), |
821 AffineTransform2D::CreateOffset(-extent.GetX1(), -extent.GetY1())); | |
822 | |
408 | 823 Render(layers, view, interpolation); |
824 | |
825 Orthanc::Image rendered(Orthanc::PixelFormat_Grayscale16, | |
826 layers.GetWidth(), layers.GetHeight(), false); | |
827 Orthanc::ImageProcessing::Convert(rendered, layers); | |
828 | |
829 std::string base64; | |
830 | |
831 { | |
832 std::string content; | |
833 | |
834 if (usePam) | |
835 { | |
836 Orthanc::PamWriter writer; | |
837 writer.WriteToMemory(content, rendered); | |
838 } | |
839 else | |
840 { | |
841 Orthanc::PngWriter writer; | |
842 writer.WriteToMemory(content, rendered); | |
843 } | |
844 | |
845 Orthanc::Toolbox::EncodeBase64(base64, content); | |
846 } | |
847 | |
848 std::set<Orthanc::DicomTag> tags; | |
849 dicom.GetTags(tags); | |
850 | |
426 | 851 createDicomRequestContent["Tags"] = Json::objectValue; |
852 | |
408 | 853 for (std::set<Orthanc::DicomTag>::const_iterator |
426 | 854 tag = tags.begin(); tag != tags.end(); ++tag) |
408 | 855 { |
856 const Orthanc::DicomValue& value = dicom.GetValue(*tag); | |
857 if (!value.IsNull() && | |
858 !value.IsBinary()) | |
859 { | |
426 | 860 createDicomRequestContent["Tags"][tag->Format()] = value.GetContent(); |
408 | 861 } |
862 } | |
863 | |
426 | 864 createDicomRequestContent["Tags"][Orthanc::DICOM_TAG_PHOTOMETRIC_INTERPRETATION.Format()] = |
865 (invert ? "MONOCHROME1" : "MONOCHROME2"); | |
408 | 866 |
867 // WARNING: The order of PixelSpacing is Y/X. We use "%0.8f" to | |
868 // avoid floating-point numbers to grow over 16 characters, | |
869 // which would be invalid according to DICOM standard | |
870 // ("dciodvfy" would complain). | |
871 char buf[32]; | |
872 sprintf(buf, "%0.8f\\%0.8f", pixelSpacingY, pixelSpacingX); | |
426 | 873 |
874 createDicomRequestContent["Tags"][Orthanc::DICOM_TAG_PIXEL_SPACING.Format()] = buf; | |
408 | 875 |
876 float center, width; | |
877 if (GetWindowing(center, width)) | |
878 { | |
426 | 879 createDicomRequestContent["Tags"][Orthanc::DICOM_TAG_WINDOW_CENTER.Format()] = |
880 boost::lexical_cast<std::string>(boost::math::iround(center)); | |
408 | 881 |
426 | 882 createDicomRequestContent["Tags"][Orthanc::DICOM_TAG_WINDOW_WIDTH.Format()] = |
883 boost::lexical_cast<std::string>(boost::math::iround(width)); | |
408 | 884 } |
885 | |
426 | 886 |
408 | 887 // This is Data URI scheme: https://en.wikipedia.org/wiki/Data_URI_scheme |
426 | 888 createDicomRequestContent["Content"] = ("data:" + |
889 std::string(usePam ? Orthanc::MIME_PAM : Orthanc::MIME_PNG) + | |
890 ";base64," + base64); | |
408 | 891 } |
892 | |
893 | |
894 void RadiographyScene::OnDicomExported(const OrthancApiClient::JsonResponseReadyMessage& message) | |
895 { | |
417
aee3d7941c9b
preparing to load images using DICOMweb
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
412
diff
changeset
|
896 LOG(INFO) << "DICOM export was successful: " |
408 | 897 << message.GetJson().toStyledString(); |
898 } | |
417
aee3d7941c9b
preparing to load images using DICOMweb
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
412
diff
changeset
|
899 |
aee3d7941c9b
preparing to load images using DICOMweb
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
412
diff
changeset
|
900 |
aee3d7941c9b
preparing to load images using DICOMweb
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
412
diff
changeset
|
901 void RadiographyScene::OnDicomWebReceived(const IWebService::HttpRequestSuccessMessage& message) |
aee3d7941c9b
preparing to load images using DICOMweb
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
412
diff
changeset
|
902 { |
aee3d7941c9b
preparing to load images using DICOMweb
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
412
diff
changeset
|
903 LOG(INFO) << "DICOMweb WADO-RS received: " << message.GetAnswerSize() << " bytes"; |
418 | 904 |
905 const IWebService::HttpHeaders& h = message.GetAnswerHttpHeaders(); | |
906 for (IWebService::HttpHeaders::const_iterator | |
426 | 907 it = h.begin(); it != h.end(); ++it) |
418 | 908 { |
909 printf("[%s] = [%s]\n", it->first.c_str(), it->second.c_str()); | |
910 } | |
417
aee3d7941c9b
preparing to load images using DICOMweb
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
412
diff
changeset
|
911 } |
aee3d7941c9b
preparing to load images using DICOMweb
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
412
diff
changeset
|
912 |
408 | 913 } |