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/ImageGeometry.h"
|
|
25 #include "../Toolbox/DicomFrameConverter.h"
|
|
26
|
|
27 #include <Core/Images/Image.h>
|
|
28 #include <Core/Images/ImageProcessing.h>
|
|
29 #include <Core/Images/PamReader.h>
|
|
30 #include <Core/Images/PamWriter.h>
|
|
31 #include <Core/Images/PngWriter.h>
|
|
32 #include <Core/OrthancException.h>
|
|
33 #include <Core/Toolbox.h>
|
|
34 #include <Plugins/Samples/Common/DicomDatasetReader.h>
|
|
35 #include <Plugins/Samples/Common/FullOrthancDataset.h>
|
|
36
|
|
37
|
|
38 namespace OrthancStone
|
|
39 {
|
|
40 static double Square(double x)
|
|
41 {
|
|
42 return x * x;
|
|
43 }
|
|
44
|
|
45
|
|
46 void RadiographyScene::Layer::UpdateTransform()
|
|
47 {
|
409
|
48 transform_ = AffineTransform2D::CreateScaling(pixelSpacingX_, pixelSpacingY_);
|
408
|
49
|
|
50 double centerX, centerY;
|
|
51 GetCenter(centerX, centerY);
|
|
52
|
409
|
53 transform_ = AffineTransform2D::Combine(
|
|
54 AffineTransform2D::CreateOffset(panX_ + centerX, panY_ + centerY),
|
|
55 AffineTransform2D::CreateRotation(angle_),
|
|
56 AffineTransform2D::CreateOffset(-centerX, -centerY),
|
408
|
57 transform_);
|
|
58
|
409
|
59 transformInverse_ = AffineTransform2D::Invert(transform_);
|
408
|
60 }
|
|
61
|
|
62
|
|
63 void RadiographyScene::Layer::AddToExtent(Extent2D& extent,
|
|
64 double x,
|
|
65 double y) const
|
|
66 {
|
409
|
67 transform_.Apply(x, y);
|
408
|
68 extent.AddPoint(x, y);
|
|
69 }
|
|
70
|
|
71
|
|
72 void RadiographyScene::Layer::GetCornerInternal(double& x,
|
|
73 double& y,
|
|
74 Corner corner,
|
|
75 unsigned int cropX,
|
|
76 unsigned int cropY,
|
|
77 unsigned int cropWidth,
|
|
78 unsigned int cropHeight) const
|
|
79 {
|
|
80 double dx = static_cast<double>(cropX);
|
|
81 double dy = static_cast<double>(cropY);
|
|
82 double dwidth = static_cast<double>(cropWidth);
|
|
83 double dheight = static_cast<double>(cropHeight);
|
|
84
|
|
85 switch (corner)
|
|
86 {
|
|
87 case Corner_TopLeft:
|
|
88 x = dx;
|
|
89 y = dy;
|
|
90 break;
|
|
91
|
|
92 case Corner_TopRight:
|
|
93 x = dx + dwidth;
|
|
94 y = dy;
|
|
95 break;
|
|
96
|
|
97 case Corner_BottomLeft:
|
|
98 x = dx;
|
|
99 y = dy + dheight;
|
|
100 break;
|
|
101
|
|
102 case Corner_BottomRight:
|
|
103 x = dx + dwidth;
|
|
104 y = dy + dheight;
|
|
105 break;
|
|
106
|
|
107 default:
|
|
108 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange);
|
|
109 }
|
|
110
|
409
|
111 transform_.Apply(x, y);
|
408
|
112 }
|
|
113
|
|
114
|
|
115 bool RadiographyScene::Layer::Contains(double x,
|
|
116 double y) const
|
|
117 {
|
409
|
118 transformInverse_.Apply(x, y);
|
408
|
119
|
|
120 unsigned int cropX, cropY, cropWidth, cropHeight;
|
|
121 GetCrop(cropX, cropY, cropWidth, cropHeight);
|
|
122
|
|
123 return (x >= cropX && x <= cropX + cropWidth &&
|
|
124 y >= cropY && y <= cropY + cropHeight);
|
|
125 }
|
|
126
|
|
127
|
|
128 void RadiographyScene::Layer::DrawBorders(CairoContext& context,
|
|
129 double zoom)
|
|
130 {
|
|
131 unsigned int cx, cy, width, height;
|
|
132 GetCrop(cx, cy, width, height);
|
|
133
|
|
134 double dx = static_cast<double>(cx);
|
|
135 double dy = static_cast<double>(cy);
|
|
136 double dwidth = static_cast<double>(width);
|
|
137 double dheight = static_cast<double>(height);
|
|
138
|
|
139 cairo_t* cr = context.GetObject();
|
|
140 cairo_set_line_width(cr, 2.0 / zoom);
|
|
141
|
|
142 double x, y;
|
|
143 x = dx;
|
|
144 y = dy;
|
409
|
145 transform_.Apply(x, y);
|
408
|
146 cairo_move_to(cr, x, y);
|
|
147
|
|
148 x = dx + dwidth;
|
|
149 y = dy;
|
409
|
150 transform_.Apply(x, y);
|
408
|
151 cairo_line_to(cr, x, y);
|
|
152
|
|
153 x = dx + dwidth;
|
|
154 y = dy + dheight;
|
409
|
155 transform_.Apply(x, y);
|
408
|
156 cairo_line_to(cr, x, y);
|
|
157
|
|
158 x = dx;
|
|
159 y = dy + dheight;
|
409
|
160 transform_.Apply(x, y);
|
408
|
161 cairo_line_to(cr, x, y);
|
|
162
|
|
163 x = dx;
|
|
164 y = dy;
|
409
|
165 transform_.Apply(x, y);
|
408
|
166 cairo_line_to(cr, x, y);
|
|
167
|
|
168 cairo_stroke(cr);
|
|
169 }
|
|
170
|
|
171
|
|
172 RadiographyScene::Layer::Layer() :
|
|
173 index_(0),
|
|
174 hasSize_(false),
|
|
175 width_(0),
|
|
176 height_(0),
|
|
177 hasCrop_(false),
|
|
178 pixelSpacingX_(1),
|
|
179 pixelSpacingY_(1),
|
|
180 panX_(0),
|
|
181 panY_(0),
|
|
182 angle_(0),
|
|
183 resizeable_(false)
|
|
184 {
|
|
185 UpdateTransform();
|
|
186 }
|
|
187
|
|
188
|
|
189 void RadiographyScene::Layer::SetCrop(unsigned int x,
|
|
190 unsigned int y,
|
|
191 unsigned int width,
|
|
192 unsigned int height)
|
|
193 {
|
|
194 if (!hasSize_)
|
|
195 {
|
|
196 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls);
|
|
197 }
|
|
198
|
|
199 if (x + width > width_ ||
|
|
200 y + height > height_)
|
|
201 {
|
|
202 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange);
|
|
203 }
|
|
204
|
|
205 hasCrop_ = true;
|
|
206 cropX_ = x;
|
|
207 cropY_ = y;
|
|
208 cropWidth_ = width;
|
|
209 cropHeight_ = height;
|
|
210
|
|
211 UpdateTransform();
|
|
212 }
|
|
213
|
|
214
|
|
215 void RadiographyScene::Layer::GetCrop(unsigned int& x,
|
|
216 unsigned int& y,
|
|
217 unsigned int& width,
|
|
218 unsigned int& height) const
|
|
219 {
|
|
220 if (hasCrop_)
|
|
221 {
|
|
222 x = cropX_;
|
|
223 y = cropY_;
|
|
224 width = cropWidth_;
|
|
225 height = cropHeight_;
|
|
226 }
|
|
227 else
|
|
228 {
|
|
229 x = 0;
|
|
230 y = 0;
|
|
231 width = width_;
|
|
232 height = height_;
|
|
233 }
|
|
234 }
|
|
235
|
|
236
|
|
237 void RadiographyScene::Layer::SetAngle(double angle)
|
|
238 {
|
|
239 angle_ = angle;
|
|
240 UpdateTransform();
|
|
241 }
|
|
242
|
|
243
|
|
244 void RadiographyScene::Layer::SetSize(unsigned int width,
|
|
245 unsigned int height)
|
|
246 {
|
|
247 if (hasSize_ &&
|
|
248 (width != width_ ||
|
|
249 height != height_))
|
|
250 {
|
|
251 throw Orthanc::OrthancException(Orthanc::ErrorCode_IncompatibleImageSize);
|
|
252 }
|
|
253
|
|
254 hasSize_ = true;
|
|
255 width_ = width;
|
|
256 height_ = height;
|
|
257
|
|
258 UpdateTransform();
|
|
259 }
|
|
260
|
|
261
|
|
262 Extent2D RadiographyScene::Layer::GetExtent() const
|
|
263 {
|
|
264 Extent2D extent;
|
|
265
|
|
266 unsigned int x, y, width, height;
|
|
267 GetCrop(x, y, width, height);
|
|
268
|
|
269 double dx = static_cast<double>(x);
|
|
270 double dy = static_cast<double>(y);
|
|
271 double dwidth = static_cast<double>(width);
|
|
272 double dheight = static_cast<double>(height);
|
|
273
|
|
274 AddToExtent(extent, dx, dy);
|
|
275 AddToExtent(extent, dx + dwidth, dy);
|
|
276 AddToExtent(extent, dx, dy + dheight);
|
|
277 AddToExtent(extent, dx + dwidth, dy + dheight);
|
|
278
|
|
279 return extent;
|
|
280 }
|
|
281
|
|
282
|
|
283 bool RadiographyScene::Layer::GetPixel(unsigned int& imageX,
|
|
284 unsigned int& imageY,
|
|
285 double sceneX,
|
|
286 double sceneY) const
|
|
287 {
|
|
288 if (width_ == 0 ||
|
|
289 height_ == 0)
|
|
290 {
|
|
291 return false;
|
|
292 }
|
|
293 else
|
|
294 {
|
409
|
295 transformInverse_.Apply(sceneX, sceneY);
|
408
|
296
|
|
297 int x = static_cast<int>(std::floor(sceneX));
|
|
298 int y = static_cast<int>(std::floor(sceneY));
|
|
299
|
|
300 if (x < 0)
|
|
301 {
|
|
302 imageX = 0;
|
|
303 }
|
|
304 else if (x >= static_cast<int>(width_))
|
|
305 {
|
|
306 imageX = width_;
|
|
307 }
|
|
308 else
|
|
309 {
|
|
310 imageX = static_cast<unsigned int>(x);
|
|
311 }
|
|
312
|
|
313 if (y < 0)
|
|
314 {
|
|
315 imageY = 0;
|
|
316 }
|
|
317 else if (y >= static_cast<int>(height_))
|
|
318 {
|
|
319 imageY = height_;
|
|
320 }
|
|
321 else
|
|
322 {
|
|
323 imageY = static_cast<unsigned int>(y);
|
|
324 }
|
|
325
|
|
326 return true;
|
|
327 }
|
|
328 }
|
|
329
|
|
330
|
|
331 void RadiographyScene::Layer::SetPan(double x,
|
|
332 double y)
|
|
333 {
|
|
334 panX_ = x;
|
|
335 panY_ = y;
|
|
336 UpdateTransform();
|
|
337 }
|
|
338
|
|
339
|
|
340 void RadiographyScene::Layer::SetPixelSpacing(double x,
|
|
341 double y)
|
|
342 {
|
|
343 pixelSpacingX_ = x;
|
|
344 pixelSpacingY_ = y;
|
|
345 UpdateTransform();
|
|
346 }
|
|
347
|
|
348
|
|
349 void RadiographyScene::Layer::GetCenter(double& centerX,
|
|
350 double& centerY) const
|
|
351 {
|
|
352 centerX = static_cast<double>(width_) / 2.0;
|
|
353 centerY = static_cast<double>(height_) / 2.0;
|
409
|
354 transform_.Apply(centerX, centerY);
|
408
|
355 }
|
|
356
|
|
357
|
|
358 void RadiographyScene::Layer::GetCorner(double& x /* out */,
|
|
359 double& y /* out */,
|
|
360 Corner corner) const
|
|
361 {
|
|
362 unsigned int cropX, cropY, cropWidth, cropHeight;
|
|
363 GetCrop(cropX, cropY, cropWidth, cropHeight);
|
|
364 GetCornerInternal(x, y, corner, cropX, cropY, cropWidth, cropHeight);
|
|
365 }
|
|
366
|
|
367
|
|
368 bool RadiographyScene::Layer::LookupCorner(Corner& corner /* out */,
|
|
369 double x,
|
|
370 double y,
|
|
371 double zoom,
|
|
372 double viewportDistance) const
|
|
373 {
|
|
374 static const Corner CORNERS[] = {
|
|
375 Corner_TopLeft,
|
|
376 Corner_TopRight,
|
|
377 Corner_BottomLeft,
|
|
378 Corner_BottomRight
|
|
379 };
|
|
380
|
|
381 unsigned int cropX, cropY, cropWidth, cropHeight;
|
|
382 GetCrop(cropX, cropY, cropWidth, cropHeight);
|
|
383
|
|
384 double threshold = Square(viewportDistance / zoom);
|
|
385
|
|
386 for (size_t i = 0; i < 4; i++)
|
|
387 {
|
|
388 double cx, cy;
|
|
389 GetCornerInternal(cx, cy, CORNERS[i], cropX, cropY, cropWidth, cropHeight);
|
|
390
|
|
391 double d = Square(cx - x) + Square(cy - y);
|
|
392
|
|
393 if (d <= threshold)
|
|
394 {
|
|
395 corner = CORNERS[i];
|
|
396 return true;
|
|
397 }
|
|
398 }
|
|
399
|
|
400 return false;
|
|
401 }
|
|
402
|
|
403
|
|
404
|
|
405 RadiographyScene::LayerAccessor::LayerAccessor(RadiographyScene& scene,
|
|
406 size_t index) :
|
|
407 scene_(scene),
|
|
408 index_(index)
|
|
409 {
|
|
410 Layers::iterator layer = scene.layers_.find(index);
|
|
411 if (layer == scene.layers_.end())
|
|
412 {
|
|
413 layer_ = NULL;
|
|
414 }
|
|
415 else
|
|
416 {
|
|
417 assert(layer->second != NULL);
|
|
418 layer_ = layer->second;
|
|
419 }
|
|
420 }
|
|
421
|
|
422
|
|
423 RadiographyScene::LayerAccessor::LayerAccessor(RadiographyScene& scene,
|
|
424 double x,
|
|
425 double y) :
|
|
426 scene_(scene),
|
|
427 index_(0) // Dummy initialization
|
|
428 {
|
|
429 if (scene.LookupLayer(index_, x, y))
|
|
430 {
|
|
431 Layers::iterator layer = scene.layers_.find(index_);
|
|
432
|
|
433 if (layer == scene.layers_.end())
|
|
434 {
|
|
435 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError);
|
|
436 }
|
|
437 else
|
|
438 {
|
|
439 assert(layer->second != NULL);
|
|
440 layer_ = layer->second;
|
|
441 }
|
|
442 }
|
|
443 else
|
|
444 {
|
|
445 layer_ = NULL;
|
|
446 }
|
|
447 }
|
|
448
|
|
449
|
|
450 RadiographyScene& RadiographyScene::LayerAccessor::GetScene() const
|
|
451 {
|
|
452 if (IsValid())
|
|
453 {
|
|
454 return scene_;
|
|
455 }
|
|
456 else
|
|
457 {
|
|
458 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls);
|
|
459 }
|
|
460 }
|
|
461
|
|
462
|
|
463 size_t RadiographyScene::LayerAccessor::GetIndex() const
|
|
464 {
|
|
465 if (IsValid())
|
|
466 {
|
|
467 return index_;
|
|
468 }
|
|
469 else
|
|
470 {
|
|
471 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls);
|
|
472 }
|
|
473 }
|
|
474
|
|
475
|
|
476 RadiographyScene::Layer& RadiographyScene::LayerAccessor::GetLayer() const
|
|
477 {
|
|
478 if (IsValid())
|
|
479 {
|
|
480 return *layer_;
|
|
481 }
|
|
482 else
|
|
483 {
|
|
484 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls);
|
|
485 }
|
|
486 }
|
|
487
|
|
488
|
|
489
|
|
490 class RadiographyScene::AlphaLayer : public Layer
|
|
491 {
|
|
492 private:
|
|
493 const RadiographyScene& scene_;
|
|
494 std::auto_ptr<Orthanc::ImageAccessor> alpha_; // Grayscale8
|
|
495 bool useWindowing_;
|
|
496 float foreground_;
|
|
497
|
|
498 public:
|
|
499 AlphaLayer(const RadiographyScene& scene) :
|
|
500 scene_(scene),
|
|
501 useWindowing_(true),
|
|
502 foreground_(0)
|
|
503 {
|
|
504 }
|
|
505
|
|
506
|
|
507 void SetForegroundValue(float foreground)
|
|
508 {
|
|
509 useWindowing_ = false;
|
|
510 foreground_ = foreground;
|
|
511 }
|
|
512
|
|
513
|
|
514 void SetAlpha(Orthanc::ImageAccessor* image)
|
|
515 {
|
|
516 std::auto_ptr<Orthanc::ImageAccessor> raii(image);
|
|
517
|
|
518 if (image == NULL)
|
|
519 {
|
|
520 throw Orthanc::OrthancException(Orthanc::ErrorCode_NullPointer);
|
|
521 }
|
|
522
|
|
523 if (image->GetFormat() != Orthanc::PixelFormat_Grayscale8)
|
|
524 {
|
|
525 throw Orthanc::OrthancException(Orthanc::ErrorCode_IncompatibleImageFormat);
|
|
526 }
|
|
527
|
|
528 SetSize(image->GetWidth(), image->GetHeight());
|
|
529 alpha_ = raii;
|
|
530 }
|
|
531
|
|
532
|
|
533 void LoadText(const Orthanc::Font& font,
|
|
534 const std::string& utf8)
|
|
535 {
|
|
536 SetAlpha(font.RenderAlpha(utf8));
|
|
537 }
|
|
538
|
|
539
|
|
540 virtual bool GetDefaultWindowing(float& center,
|
|
541 float& width) const
|
|
542 {
|
|
543 return false;
|
|
544 }
|
|
545
|
|
546
|
|
547 virtual void Render(Orthanc::ImageAccessor& buffer,
|
409
|
548 const AffineTransform2D& viewTransform,
|
408
|
549 ImageInterpolation interpolation) const
|
|
550 {
|
|
551 if (alpha_.get() == NULL)
|
|
552 {
|
|
553 return;
|
|
554 }
|
|
555
|
|
556 if (buffer.GetFormat() != Orthanc::PixelFormat_Float32)
|
|
557 {
|
|
558 throw Orthanc::OrthancException(Orthanc::ErrorCode_IncompatibleImageFormat);
|
|
559 }
|
|
560
|
|
561 unsigned int cropX, cropY, cropWidth, cropHeight;
|
|
562 GetCrop(cropX, cropY, cropWidth, cropHeight);
|
|
563
|
409
|
564 const AffineTransform2D t = AffineTransform2D::Combine(
|
|
565 viewTransform, GetTransform(),
|
|
566 AffineTransform2D::CreateOffset(cropX, cropY));
|
408
|
567
|
|
568 Orthanc::ImageAccessor cropped;
|
|
569 alpha_->GetRegion(cropped, cropX, cropY, cropWidth, cropHeight);
|
|
570
|
|
571 Orthanc::Image tmp(Orthanc::PixelFormat_Grayscale8, buffer.GetWidth(), buffer.GetHeight(), false);
|
409
|
572
|
|
573 t.Apply(tmp, cropped, interpolation, true /* clear */);
|
408
|
574
|
|
575 // Blit
|
|
576 const unsigned int width = buffer.GetWidth();
|
|
577 const unsigned int height = buffer.GetHeight();
|
|
578
|
|
579 float value = foreground_;
|
|
580
|
|
581 if (useWindowing_)
|
|
582 {
|
|
583 float center, width;
|
|
584 if (scene_.GetWindowing(center, width))
|
|
585 {
|
|
586 value = center + width / 2.0f;
|
|
587 }
|
|
588 }
|
|
589
|
|
590 for (unsigned int y = 0; y < height; y++)
|
|
591 {
|
|
592 float *q = reinterpret_cast<float*>(buffer.GetRow(y));
|
|
593 const uint8_t *p = reinterpret_cast<uint8_t*>(tmp.GetRow(y));
|
|
594
|
|
595 for (unsigned int x = 0; x < width; x++, p++, q++)
|
|
596 {
|
|
597 float a = static_cast<float>(*p) / 255.0f;
|
|
598
|
|
599 *q = (a * value + (1.0f - a) * (*q));
|
|
600 }
|
|
601 }
|
|
602 }
|
|
603
|
|
604
|
|
605 virtual bool GetRange(float& minValue,
|
|
606 float& maxValue) const
|
|
607 {
|
|
608 if (useWindowing_)
|
|
609 {
|
|
610 return false;
|
|
611 }
|
|
612 else
|
|
613 {
|
|
614 minValue = 0;
|
|
615 maxValue = 0;
|
|
616
|
|
617 if (foreground_ < 0)
|
|
618 {
|
|
619 minValue = foreground_;
|
|
620 }
|
|
621
|
|
622 if (foreground_ > 0)
|
|
623 {
|
|
624 maxValue = foreground_;
|
|
625 }
|
|
626
|
|
627 return true;
|
|
628 }
|
|
629 }
|
|
630 };
|
|
631
|
|
632
|
|
633
|
|
634 class RadiographyScene::DicomLayer : public Layer
|
|
635 {
|
|
636 private:
|
|
637 std::auto_ptr<Orthanc::ImageAccessor> source_; // Content of PixelData
|
|
638 std::auto_ptr<DicomFrameConverter> converter_;
|
|
639 std::auto_ptr<Orthanc::ImageAccessor> converted_; // Float32
|
|
640
|
|
641 static OrthancPlugins::DicomTag ConvertTag(const Orthanc::DicomTag& tag)
|
|
642 {
|
|
643 return OrthancPlugins::DicomTag(tag.GetGroup(), tag.GetElement());
|
|
644 }
|
|
645
|
|
646
|
|
647 void ApplyConverter()
|
|
648 {
|
|
649 if (source_.get() != NULL &&
|
|
650 converter_.get() != NULL)
|
|
651 {
|
|
652 converted_.reset(converter_->ConvertFrame(*source_));
|
|
653 }
|
|
654 }
|
|
655
|
|
656 public:
|
|
657 void SetDicomTags(const OrthancPlugins::FullOrthancDataset& dataset)
|
|
658 {
|
|
659 converter_.reset(new DicomFrameConverter);
|
|
660 converter_->ReadParameters(dataset);
|
|
661 ApplyConverter();
|
|
662
|
|
663 std::string tmp;
|
|
664 Vector pixelSpacing;
|
|
665
|
|
666 if (dataset.GetStringValue(tmp, ConvertTag(Orthanc::DICOM_TAG_PIXEL_SPACING)) &&
|
|
667 LinearAlgebra::ParseVector(pixelSpacing, tmp) &&
|
|
668 pixelSpacing.size() == 2)
|
|
669 {
|
|
670 SetPixelSpacing(pixelSpacing[0], pixelSpacing[1]);
|
|
671 }
|
|
672
|
|
673 //SetPan(-0.5 * GetPixelSpacingX(), -0.5 * GetPixelSpacingY());
|
|
674
|
|
675 OrthancPlugins::DicomDatasetReader reader(dataset);
|
|
676
|
|
677 unsigned int width, height;
|
|
678 if (!reader.GetUnsignedIntegerValue(width, ConvertTag(Orthanc::DICOM_TAG_COLUMNS)) ||
|
|
679 !reader.GetUnsignedIntegerValue(height, ConvertTag(Orthanc::DICOM_TAG_ROWS)))
|
|
680 {
|
|
681 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat);
|
|
682 }
|
|
683 else
|
|
684 {
|
|
685 SetSize(width, height);
|
|
686 }
|
|
687 }
|
|
688
|
|
689
|
|
690 void SetSourceImage(Orthanc::ImageAccessor* image) // Takes ownership
|
|
691 {
|
|
692 std::auto_ptr<Orthanc::ImageAccessor> raii(image);
|
|
693
|
|
694 if (image == NULL)
|
|
695 {
|
|
696 throw Orthanc::OrthancException(Orthanc::ErrorCode_NullPointer);
|
|
697 }
|
|
698
|
|
699 SetSize(image->GetWidth(), image->GetHeight());
|
|
700
|
|
701 source_ = raii;
|
|
702 ApplyConverter();
|
|
703 }
|
|
704
|
|
705
|
|
706 virtual void Render(Orthanc::ImageAccessor& buffer,
|
409
|
707 const AffineTransform2D& viewTransform,
|
408
|
708 ImageInterpolation interpolation) const
|
|
709 {
|
|
710 if (converted_.get() != NULL)
|
|
711 {
|
|
712 if (converted_->GetFormat() != Orthanc::PixelFormat_Float32)
|
|
713 {
|
|
714 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError);
|
|
715 }
|
|
716
|
|
717 unsigned int cropX, cropY, cropWidth, cropHeight;
|
|
718 GetCrop(cropX, cropY, cropWidth, cropHeight);
|
|
719
|
409
|
720 AffineTransform2D t = AffineTransform2D::Combine(
|
|
721 viewTransform, GetTransform(),
|
|
722 AffineTransform2D::CreateOffset(cropX, cropY));
|
408
|
723
|
|
724 Orthanc::ImageAccessor cropped;
|
|
725 converted_->GetRegion(cropped, cropX, cropY, cropWidth, cropHeight);
|
|
726
|
409
|
727 t.Apply(buffer, cropped, interpolation, false);
|
408
|
728 }
|
|
729 }
|
|
730
|
|
731
|
|
732 virtual bool GetDefaultWindowing(float& center,
|
|
733 float& width) const
|
|
734 {
|
|
735 if (converter_.get() != NULL &&
|
|
736 converter_->HasDefaultWindow())
|
|
737 {
|
|
738 center = static_cast<float>(converter_->GetDefaultWindowCenter());
|
|
739 width = static_cast<float>(converter_->GetDefaultWindowWidth());
|
|
740 return true;
|
|
741 }
|
|
742 else
|
|
743 {
|
|
744 return false;
|
|
745 }
|
|
746 }
|
|
747
|
|
748
|
|
749 virtual bool GetRange(float& minValue,
|
|
750 float& maxValue) const
|
|
751 {
|
|
752 if (converted_.get() != NULL)
|
|
753 {
|
|
754 if (converted_->GetFormat() != Orthanc::PixelFormat_Float32)
|
|
755 {
|
|
756 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError);
|
|
757 }
|
|
758
|
|
759 Orthanc::ImageProcessing::GetMinMaxFloatValue(minValue, maxValue, *converted_);
|
|
760 return true;
|
|
761 }
|
|
762 else
|
|
763 {
|
|
764 return false;
|
|
765 }
|
|
766 }
|
|
767 };
|
|
768
|
|
769
|
|
770 RadiographyScene::Layer& RadiographyScene::RegisterLayer(RadiographyScene::Layer* layer)
|
|
771 {
|
|
772 if (layer == NULL)
|
|
773 {
|
|
774 throw Orthanc::OrthancException(Orthanc::ErrorCode_NullPointer);
|
|
775 }
|
|
776
|
|
777 std::auto_ptr<Layer> raii(layer);
|
|
778
|
|
779 size_t index = countLayers_++;
|
|
780 raii->SetIndex(index);
|
|
781 layers_[index] = raii.release();
|
|
782
|
|
783 EmitMessage(GeometryChangedMessage(*this));
|
|
784 EmitMessage(ContentChangedMessage(*this));
|
|
785
|
|
786 return *layer;
|
|
787 }
|
|
788
|
|
789
|
|
790 RadiographyScene::RadiographyScene(MessageBroker& broker,
|
|
791 OrthancApiClient& orthanc) :
|
|
792 IObserver(broker),
|
|
793 IObservable(broker),
|
|
794 orthanc_(orthanc),
|
|
795 countLayers_(0),
|
|
796 hasWindowing_(false),
|
|
797 windowingCenter_(0), // Dummy initialization
|
|
798 windowingWidth_(0) // Dummy initialization
|
|
799 {
|
|
800 }
|
|
801
|
|
802
|
|
803 RadiographyScene::~RadiographyScene()
|
|
804 {
|
|
805 for (Layers::iterator it = layers_.begin(); it != layers_.end(); it++)
|
|
806 {
|
|
807 assert(it->second != NULL);
|
|
808 delete it->second;
|
|
809 }
|
|
810 }
|
|
811
|
|
812
|
|
813 bool RadiographyScene::GetWindowing(float& center,
|
|
814 float& width) const
|
|
815 {
|
|
816 if (hasWindowing_)
|
|
817 {
|
|
818 center = windowingCenter_;
|
|
819 width = windowingWidth_;
|
|
820 return true;
|
|
821 }
|
|
822 else
|
|
823 {
|
|
824 return false;
|
|
825 }
|
|
826 }
|
|
827
|
|
828
|
|
829 void RadiographyScene::GetWindowingWithDefault(float& center,
|
|
830 float& width) const
|
|
831 {
|
|
832 if (!GetWindowing(center, width))
|
|
833 {
|
|
834 center = 128;
|
|
835 width = 256;
|
|
836 }
|
|
837 }
|
|
838
|
|
839
|
|
840 void RadiographyScene::SetWindowing(float center,
|
|
841 float width)
|
|
842 {
|
|
843 hasWindowing_ = true;
|
|
844 windowingCenter_ = center;
|
|
845 windowingWidth_ = width;
|
|
846 }
|
|
847
|
|
848
|
|
849 RadiographyScene::Layer& RadiographyScene::LoadText(const Orthanc::Font& font,
|
|
850 const std::string& utf8)
|
|
851 {
|
|
852 std::auto_ptr<AlphaLayer> alpha(new AlphaLayer(*this));
|
|
853 alpha->LoadText(font, utf8);
|
|
854
|
|
855 return RegisterLayer(alpha.release());
|
|
856 }
|
|
857
|
|
858
|
|
859 RadiographyScene::Layer& RadiographyScene::LoadTestBlock(unsigned int width,
|
|
860 unsigned int height)
|
|
861 {
|
|
862 std::auto_ptr<Orthanc::Image> block(new Orthanc::Image(Orthanc::PixelFormat_Grayscale8, width, height, false));
|
|
863
|
|
864 for (unsigned int padding = 0;
|
|
865 (width > 2 * padding) && (height > 2 * padding);
|
|
866 padding++)
|
|
867 {
|
|
868 uint8_t color;
|
|
869 if (255 > 10 * padding)
|
|
870 {
|
|
871 color = 255 - 10 * padding;
|
|
872 }
|
|
873 else
|
|
874 {
|
|
875 color = 0;
|
|
876 }
|
|
877
|
|
878 Orthanc::ImageAccessor region;
|
|
879 block->GetRegion(region, padding, padding, width - 2 * padding, height - 2 * padding);
|
|
880 Orthanc::ImageProcessing::Set(region, color);
|
|
881 }
|
|
882
|
|
883 std::auto_ptr<AlphaLayer> alpha(new AlphaLayer(*this));
|
|
884 alpha->SetAlpha(block.release());
|
|
885
|
|
886 return RegisterLayer(alpha.release());
|
|
887 }
|
|
888
|
|
889
|
|
890 RadiographyScene::Layer& RadiographyScene::LoadDicomFrame(const std::string& instance,
|
|
891 unsigned int frame,
|
|
892 bool httpCompression)
|
|
893 {
|
|
894 Layer& layer = RegisterLayer(new DicomLayer);
|
|
895
|
|
896 {
|
|
897 IWebService::Headers headers;
|
|
898 std::string uri = "/instances/" + instance + "/tags";
|
|
899
|
|
900 orthanc_.GetBinaryAsync(
|
|
901 uri, headers,
|
|
902 new Callable<RadiographyScene, OrthancApiClient::BinaryResponseReadyMessage>
|
|
903 (*this, &RadiographyScene::OnTagsReceived), NULL,
|
|
904 new Orthanc::SingleValueObject<size_t>(layer.GetIndex()));
|
|
905 }
|
|
906
|
|
907 {
|
|
908 IWebService::Headers headers;
|
|
909 headers["Accept"] = "image/x-portable-arbitrarymap";
|
|
910
|
|
911 if (httpCompression)
|
|
912 {
|
|
913 headers["Accept-Encoding"] = "gzip";
|
|
914 }
|
|
915
|
|
916 std::string uri = ("/instances/" + instance + "/frames/" +
|
|
917 boost::lexical_cast<std::string>(frame) + "/image-uint16");
|
|
918
|
|
919 orthanc_.GetBinaryAsync(
|
|
920 uri, headers,
|
|
921 new Callable<RadiographyScene, OrthancApiClient::BinaryResponseReadyMessage>
|
|
922 (*this, &RadiographyScene::OnFrameReceived), NULL,
|
|
923 new Orthanc::SingleValueObject<size_t>(layer.GetIndex()));
|
|
924 }
|
|
925
|
|
926 return layer;
|
|
927 }
|
|
928
|
|
929
|
|
930 void RadiographyScene::OnTagsReceived(const OrthancApiClient::BinaryResponseReadyMessage& message)
|
|
931 {
|
|
932 size_t index = dynamic_cast<const Orthanc::SingleValueObject<size_t>&>
|
|
933 (message.GetPayload()).GetValue();
|
|
934
|
|
935 LOG(INFO) << "JSON received: " << message.GetUri().c_str()
|
|
936 << " (" << message.GetAnswerSize() << " bytes) for layer " << index;
|
|
937
|
|
938 Layers::iterator layer = layers_.find(index);
|
|
939 if (layer != layers_.end())
|
|
940 {
|
|
941 assert(layer->second != NULL);
|
|
942
|
|
943 OrthancPlugins::FullOrthancDataset dicom(message.GetAnswer(), message.GetAnswerSize());
|
|
944 dynamic_cast<DicomLayer*>(layer->second)->SetDicomTags(dicom);
|
|
945
|
|
946 float c, w;
|
|
947 if (!hasWindowing_ &&
|
|
948 layer->second->GetDefaultWindowing(c, w))
|
|
949 {
|
|
950 hasWindowing_ = true;
|
|
951 windowingCenter_ = c;
|
|
952 windowingWidth_ = w;
|
|
953 }
|
|
954
|
|
955 EmitMessage(GeometryChangedMessage(*this));
|
|
956 }
|
|
957 }
|
|
958
|
|
959
|
|
960 void RadiographyScene::OnFrameReceived(const OrthancApiClient::BinaryResponseReadyMessage& message)
|
|
961 {
|
|
962 size_t index = dynamic_cast<const Orthanc::SingleValueObject<size_t>&>(message.GetPayload()).GetValue();
|
|
963
|
|
964 LOG(INFO) << "DICOM frame received: " << message.GetUri().c_str()
|
|
965 << " (" << message.GetAnswerSize() << " bytes) for layer " << index;
|
|
966
|
|
967 Layers::iterator layer = layers_.find(index);
|
|
968 if (layer != layers_.end())
|
|
969 {
|
|
970 assert(layer->second != NULL);
|
|
971
|
|
972 std::string content;
|
|
973 if (message.GetAnswerSize() > 0)
|
|
974 {
|
|
975 content.assign(reinterpret_cast<const char*>(message.GetAnswer()), message.GetAnswerSize());
|
|
976 }
|
|
977
|
|
978 std::auto_ptr<Orthanc::PamReader> reader(new Orthanc::PamReader);
|
|
979 reader->ReadFromMemory(content);
|
|
980 dynamic_cast<DicomLayer*>(layer->second)->SetSourceImage(reader.release());
|
|
981
|
|
982 EmitMessage(ContentChangedMessage(*this));
|
|
983 }
|
|
984 }
|
|
985
|
|
986
|
|
987 Extent2D RadiographyScene::GetSceneExtent() const
|
|
988 {
|
|
989 Extent2D extent;
|
|
990
|
|
991 for (Layers::const_iterator it = layers_.begin();
|
|
992 it != layers_.end(); ++it)
|
|
993 {
|
|
994 assert(it->second != NULL);
|
|
995 extent.Union(it->second->GetExtent());
|
|
996 }
|
|
997
|
|
998 return extent;
|
|
999 }
|
|
1000
|
|
1001
|
|
1002 void RadiographyScene::Render(Orthanc::ImageAccessor& buffer,
|
409
|
1003 const AffineTransform2D& viewTransform,
|
408
|
1004 ImageInterpolation interpolation) const
|
|
1005 {
|
|
1006 Orthanc::ImageProcessing::Set(buffer, 0);
|
|
1007
|
|
1008 // Render layers in the background-to-foreground order
|
|
1009 for (size_t index = 0; index < countLayers_; index++)
|
|
1010 {
|
|
1011 Layers::const_iterator it = layers_.find(index);
|
|
1012 if (it != layers_.end())
|
|
1013 {
|
|
1014 assert(it->second != NULL);
|
|
1015 it->second->Render(buffer, viewTransform, interpolation);
|
|
1016 }
|
|
1017 }
|
|
1018 }
|
|
1019
|
|
1020
|
|
1021 bool RadiographyScene::LookupLayer(size_t& index /* out */,
|
|
1022 double x,
|
|
1023 double y) const
|
|
1024 {
|
|
1025 // Render layers in the foreground-to-background order
|
|
1026 for (size_t i = countLayers_; i > 0; i--)
|
|
1027 {
|
|
1028 index = i - 1;
|
|
1029 Layers::const_iterator it = layers_.find(index);
|
|
1030 if (it != layers_.end())
|
|
1031 {
|
|
1032 assert(it->second != NULL);
|
|
1033 if (it->second->Contains(x, y))
|
|
1034 {
|
|
1035 return true;
|
|
1036 }
|
|
1037 }
|
|
1038 }
|
|
1039
|
|
1040 return false;
|
|
1041 }
|
|
1042
|
|
1043
|
|
1044 void RadiographyScene::DrawBorder(CairoContext& context,
|
|
1045 unsigned int layer,
|
|
1046 double zoom)
|
|
1047 {
|
|
1048 Layers::const_iterator found = layers_.find(layer);
|
|
1049
|
|
1050 if (found != layers_.end())
|
|
1051 {
|
|
1052 context.SetSourceColor(255, 0, 0);
|
|
1053 found->second->DrawBorders(context, zoom);
|
|
1054 }
|
|
1055 }
|
|
1056
|
|
1057
|
|
1058 void RadiographyScene::GetRange(float& minValue,
|
|
1059 float& maxValue) const
|
|
1060 {
|
|
1061 bool first = true;
|
|
1062
|
|
1063 for (Layers::const_iterator it = layers_.begin();
|
|
1064 it != layers_.end(); it++)
|
|
1065 {
|
|
1066 assert(it->second != NULL);
|
|
1067
|
|
1068 float a, b;
|
|
1069 if (it->second->GetRange(a, b))
|
|
1070 {
|
|
1071 if (first)
|
|
1072 {
|
|
1073 minValue = a;
|
|
1074 maxValue = b;
|
|
1075 first = false;
|
|
1076 }
|
|
1077 else
|
|
1078 {
|
|
1079 minValue = std::min(a, minValue);
|
|
1080 maxValue = std::max(b, maxValue);
|
|
1081 }
|
|
1082 }
|
|
1083 }
|
|
1084
|
|
1085 if (first)
|
|
1086 {
|
|
1087 minValue = 0;
|
|
1088 maxValue = 0;
|
|
1089 }
|
|
1090 }
|
|
1091
|
|
1092
|
|
1093 // Export using PAM is faster than using PNG, but requires Orthanc
|
|
1094 // core >= 1.4.3
|
|
1095 void RadiographyScene::ExportDicom(const Orthanc::DicomMap& dicom,
|
|
1096 double pixelSpacingX,
|
|
1097 double pixelSpacingY,
|
|
1098 bool invert,
|
|
1099 ImageInterpolation interpolation,
|
|
1100 bool usePam)
|
|
1101 {
|
|
1102 if (pixelSpacingX <= 0 ||
|
|
1103 pixelSpacingY <= 0)
|
|
1104 {
|
|
1105 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange);
|
|
1106 }
|
|
1107
|
|
1108 LOG(INFO) << "Exporting DICOM";
|
|
1109
|
|
1110 Extent2D extent = GetSceneExtent();
|
|
1111
|
|
1112 int w = std::ceil(extent.GetWidth() / pixelSpacingX);
|
|
1113 int h = std::ceil(extent.GetHeight() / pixelSpacingY);
|
|
1114
|
|
1115 if (w < 0 || h < 0)
|
|
1116 {
|
|
1117 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError);
|
|
1118 }
|
|
1119
|
|
1120 Orthanc::Image layers(Orthanc::PixelFormat_Float32,
|
|
1121 static_cast<unsigned int>(w),
|
|
1122 static_cast<unsigned int>(h), false);
|
|
1123
|
409
|
1124 AffineTransform2D view = AffineTransform2D::Combine(
|
|
1125 AffineTransform2D::CreateScaling(1.0 / pixelSpacingX, 1.0 / pixelSpacingY),
|
|
1126 AffineTransform2D::CreateOffset(-extent.GetX1(), -extent.GetY1()));
|
408
|
1127
|
|
1128 Render(layers, view, interpolation);
|
|
1129
|
|
1130 Orthanc::Image rendered(Orthanc::PixelFormat_Grayscale16,
|
|
1131 layers.GetWidth(), layers.GetHeight(), false);
|
|
1132 Orthanc::ImageProcessing::Convert(rendered, layers);
|
|
1133
|
|
1134 std::string base64;
|
|
1135
|
|
1136 {
|
|
1137 std::string content;
|
|
1138
|
|
1139 if (usePam)
|
|
1140 {
|
|
1141 Orthanc::PamWriter writer;
|
|
1142 writer.WriteToMemory(content, rendered);
|
|
1143 }
|
|
1144 else
|
|
1145 {
|
|
1146 Orthanc::PngWriter writer;
|
|
1147 writer.WriteToMemory(content, rendered);
|
|
1148 }
|
|
1149
|
|
1150 Orthanc::Toolbox::EncodeBase64(base64, content);
|
|
1151 }
|
|
1152
|
|
1153 std::set<Orthanc::DicomTag> tags;
|
|
1154 dicom.GetTags(tags);
|
|
1155
|
|
1156 Json::Value json = Json::objectValue;
|
|
1157 json["Tags"] = Json::objectValue;
|
|
1158
|
|
1159 for (std::set<Orthanc::DicomTag>::const_iterator
|
|
1160 tag = tags.begin(); tag != tags.end(); ++tag)
|
|
1161 {
|
|
1162 const Orthanc::DicomValue& value = dicom.GetValue(*tag);
|
|
1163 if (!value.IsNull() &&
|
|
1164 !value.IsBinary())
|
|
1165 {
|
|
1166 json["Tags"][tag->Format()] = value.GetContent();
|
|
1167 }
|
|
1168 }
|
|
1169
|
|
1170 json["Tags"][Orthanc::DICOM_TAG_PHOTOMETRIC_INTERPRETATION.Format()] =
|
|
1171 (invert ? "MONOCHROME1" : "MONOCHROME2");
|
|
1172
|
|
1173 // WARNING: The order of PixelSpacing is Y/X. We use "%0.8f" to
|
|
1174 // avoid floating-point numbers to grow over 16 characters,
|
|
1175 // which would be invalid according to DICOM standard
|
|
1176 // ("dciodvfy" would complain).
|
|
1177 char buf[32];
|
|
1178 sprintf(buf, "%0.8f\\%0.8f", pixelSpacingY, pixelSpacingX);
|
|
1179
|
|
1180 json["Tags"][Orthanc::DICOM_TAG_PIXEL_SPACING.Format()] = buf;
|
|
1181
|
|
1182 float center, width;
|
|
1183 if (GetWindowing(center, width))
|
|
1184 {
|
|
1185 json["Tags"][Orthanc::DICOM_TAG_WINDOW_CENTER.Format()] =
|
|
1186 boost::lexical_cast<std::string>(boost::math::iround(center));
|
|
1187
|
|
1188 json["Tags"][Orthanc::DICOM_TAG_WINDOW_WIDTH.Format()] =
|
|
1189 boost::lexical_cast<std::string>(boost::math::iround(width));
|
|
1190 }
|
|
1191
|
|
1192 // This is Data URI scheme: https://en.wikipedia.org/wiki/Data_URI_scheme
|
|
1193 json["Content"] = ("data:" +
|
|
1194 std::string(usePam ? Orthanc::MIME_PAM : Orthanc::MIME_PNG) +
|
|
1195 ";base64," + base64);
|
|
1196
|
|
1197 orthanc_.PostJsonAsyncExpectJson(
|
|
1198 "/tools/create-dicom", json,
|
|
1199 new Callable<RadiographyScene, OrthancApiClient::JsonResponseReadyMessage>
|
|
1200 (*this, &RadiographyScene::OnDicomExported),
|
|
1201 NULL, NULL);
|
|
1202 }
|
|
1203
|
|
1204
|
|
1205 void RadiographyScene::OnDicomExported(const OrthancApiClient::JsonResponseReadyMessage& message)
|
|
1206 {
|
|
1207 LOG(INFO) << "DICOM export was successful:"
|
|
1208 << message.GetJson().toStyledString();
|
|
1209 }
|
|
1210 }
|