Mercurial > hg > orthanc-stone

annotate Framework/Radiography/RadiographyMaskLayer.cpp @ 478:017a66a2999b am-touch-events

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
cleanup + avoid conversions
author am@osimis.io
date Wed, 13 Feb 2019 14:12:41 +0100
parents baf8c8d68bbc
children e3d316ba34ba
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
475
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
1 /**
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
2 * Stone of Orthanc
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
3 * Copyright (C) 2012-2016 Sebastien Jodogne, Medical Physics
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
4 * Department, University Hospital of Liege, Belgium
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
5 * Copyright (C) 2017-2018 Osimis S.A., Belgium
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
6 *
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
7 * This program is free software: you can redistribute it and/or
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
8 * modify it under the terms of the GNU Affero General Public License
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
9 * as published by the Free Software Foundation, either version 3 of
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
10 * the License, or (at your option) any later version.
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
11 *
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
12 * This program is distributed in the hope that it will be useful, but
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
13 * WITHOUT ANY WARRANTY; without even the implied warranty of
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
15 * Affero General Public License for more details.
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
16 *
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
17 * You should have received a copy of the GNU Affero General Public License
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
18 * along with this program. If not, see <http://www.gnu.org/licenses/>.
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
19 **/
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
20
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
21
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
22 #include "RadiographyMaskLayer.h"
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
23 #include "RadiographyDicomLayer.h"
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
24
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
25 #include "RadiographyScene.h"
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
26 #include "Core/Images/Image.h"
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
27 #include "Core/Images/ImageProcessing.h"
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
28 #include <Core/OrthancException.h>
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
29
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
30 namespace OrthancStone
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
31 {
476
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
32 const unsigned char IN_MASK_VALUE = 0x00;
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
33 const unsigned char OUT_MASK_VALUE = 0xFF;
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
34
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
35 const AffineTransform2D& RadiographyMaskLayer::GetTransform() const
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
36 {
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
37 return dicomLayer_.GetTransform();
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
38 }
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
39
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
40 const AffineTransform2D& RadiographyMaskLayer::GetTransformInverse() const
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
41 {
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
42 return dicomLayer_.GetTransformInverse();
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
43 }
475
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
44
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
45 void ComputeMaskExtent(unsigned int& left, unsigned int& right, unsigned int& top, unsigned int& bottom, const std::vector<MaskPoint>& corners)
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
46 {
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
47 left = std::numeric_limits<unsigned int>::max();
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
48 right = std::numeric_limits<unsigned int>::min();
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
49 top = std::numeric_limits<unsigned int>::max();
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
50 bottom = std::numeric_limits<unsigned int>::min();
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
51
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
52 for (size_t i = 0; i < corners.size(); i++)
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
53 {
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
54 const MaskPoint& p = corners[i];
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
55 left = std::min(p.x, left);
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
56 right = std::max(p.x, right);
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
57 bottom = std::max(p.y, bottom);
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
58 top = std::min(p.y, top);
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
59 }
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
60 }
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
61
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
62 void RadiographyMaskLayer::SetCorners(const std::vector<MaskPoint>& corners)
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
63 {
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
64 corners_ = corners;
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
65 invalidated_ = true;
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
66 }
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
67
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
68 void RadiographyMaskLayer::Render(Orthanc::ImageAccessor& buffer,
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
69 const AffineTransform2D& viewTransform,
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
70 ImageInterpolation interpolation) const
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
71 {
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
72 if (dicomLayer_.GetWidth() == 0) // nothing to do if the DICOM layer is not displayed (or not loaded)
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
73 return;
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
74
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
75 if (invalidated_)
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
76 {
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
77 mask_.reset(new Orthanc::Image(Orthanc::PixelFormat_Grayscale8, dicomLayer_.GetWidth(), dicomLayer_.GetHeight(), false));
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
78
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
79 DrawMask();
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
80
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
81 invalidated_ = false;
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
82 }
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
83
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
84 {// rendering
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
85 if (buffer.GetFormat() != Orthanc::PixelFormat_Float32)
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
86 {
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
87 throw Orthanc::OrthancException(Orthanc::ErrorCode_IncompatibleImageFormat);
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
88 }
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
89
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
90 unsigned int cropX, cropY, cropWidth, cropHeight;
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
91 dicomLayer_.GetCrop(cropX, cropY, cropWidth, cropHeight);
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
92
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
93 const AffineTransform2D t = AffineTransform2D::Combine(
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
94 viewTransform, dicomLayer_.GetTransform(),
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
95 AffineTransform2D::CreateOffset(cropX, cropY));
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
96
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
97 Orthanc::ImageAccessor cropped;
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
98 mask_->GetRegion(cropped, cropX, cropY, cropWidth, cropHeight);
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
99
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
100 Orthanc::Image tmp(Orthanc::PixelFormat_Grayscale8, buffer.GetWidth(), buffer.GetHeight(), false);
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
101
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
102 t.Apply(tmp, cropped, interpolation, true /* clear */);
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
103
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
104 // Blit
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
105 const unsigned int width = buffer.GetWidth();
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
106 const unsigned int height = buffer.GetHeight();
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
107
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
108 for (unsigned int y = 0; y < height; y++)
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
109 {
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
110 float *q = reinterpret_cast<float*>(buffer.GetRow(y));
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
111 const uint8_t *p = reinterpret_cast<uint8_t*>(tmp.GetRow(y));
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
112
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
113 for (unsigned int x = 0; x < width; x++, p++, q++)
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
114 {
476
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
115 if (*p == OUT_MASK_VALUE)
475
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
116 *q = foreground_;
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
117 // else keep the underlying pixel value
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
118 }
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
119 }
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
120
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
121 }
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
122 }
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
123
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
124 // from https://www.geeksforgeeks.org/how-to-check-if-a-given-point-lies-inside-a-polygon/
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
125 // Given three colinear points p, q, r, the function checks if
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
126 // point q lies on line segment 'pr'
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
127 bool onSegment(const MaskPoint& p, const MaskPoint& q, const MaskPoint& r)
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
128 {
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
129 if (q.x <= std::max(p.x, r.x) && q.x >= std::min(p.x, r.x) &&
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
130 q.y <= std::max(p.y, r.y) && q.y >= std::min(p.y, r.y))
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
131 return true;
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
132 return false;
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
133 }
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
134
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
135 // To find orientation of ordered triplet (p, q, r).
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
136 // The function returns following values
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
137 // 0 --> p, q and r are colinear
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
138 // 1 --> Clockwise
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
139 // 2 --> Counterclockwise
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
140 int orientation(const MaskPoint& p, const MaskPoint& q, const MaskPoint& r)
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
141 {
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
142 int val = (q.y - p.y) * (r.x - q.x) -
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
143 (q.x - p.x) * (r.y - q.y);
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
144
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
145 if (val == 0) return 0; // colinear
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
146 return (val > 0)? 1: 2; // clock or counterclock wise
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
147 }
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
148
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
149 // The function that returns true if line segment 'p1q1'
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
150 // and 'p2q2' intersect.
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
151 bool doIntersect(const MaskPoint& p1, const MaskPoint& q1, const MaskPoint& p2, const MaskPoint& q2)
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
152 {
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
153 // Find the four orientations needed for general and
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
154 // special cases
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
155 int o1 = orientation(p1, q1, p2);
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
156 int o2 = orientation(p1, q1, q2);
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
157 int o3 = orientation(p2, q2, p1);
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
158 int o4 = orientation(p2, q2, q1);
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
159
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
160 // General case
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
161 if (o1 != o2 && o3 != o4)
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
162 return true;
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
163
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
164 // Special Cases
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
165 // p1, q1 and p2 are colinear and p2 lies on segment p1q1
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
166 if (o1 == 0 && onSegment(p1, p2, q1)) return true;
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
167
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
168 // p1, q1 and p2 are colinear and q2 lies on segment p1q1
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
169 if (o2 == 0 && onSegment(p1, q2, q1)) return true;
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
170
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
171 // p2, q2 and p1 are colinear and p1 lies on segment p2q2
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
172 if (o3 == 0 && onSegment(p2, p1, q2)) return true;
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
173
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
174 // p2, q2 and q1 are colinear and q1 lies on segment p2q2
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
175 if (o4 == 0 && onSegment(p2, q1, q2)) return true;
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
176
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
177 return false; // Doesn't fall in any of the above cases
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
178 }
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
179
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
180 // Define Infinite (Using INT_MAX caused overflow problems)
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
181 #define MASK_INF 1000000
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
182
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
183 // Returns true if the point p lies inside the polygon[] with n vertices
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
184 bool isInside(const std::vector<MaskPoint>& polygon, const MaskPoint& p)
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
185 {
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
186 // There must be at least 3 vertices in polygon[]
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
187 if (polygon.size() < 3) return false;
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
188
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
189 // Create a point for line segment from p to infinite
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
190 MaskPoint extreme = {MASK_INF, p.y};
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
191
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
192 // Count intersections of the above line with sides of polygon
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
193 int count = 0, i = 0;
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
194 do
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
195 {
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
196 int next = (i+1) % polygon.size();
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
197
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
198 // Check if the line segment from 'p' to 'extreme' intersects
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
199 // with the line segment from 'polygon[i]' to 'polygon[next]'
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
200 if (doIntersect(polygon[i], polygon[next], p, extreme))
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
201 {
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
202 // If the point 'p' is colinear with line segment 'i-next',
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
203 // then check if it lies on segment. If it lies, return true,
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
204 // otherwise false
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
205 if (orientation(polygon[i], p, polygon[next]) == 0)
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
206 return onSegment(polygon[i], p, polygon[next]);
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
207
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
208 count++;
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
209 }
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
210 i = next;
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
211 } while (i != 0);
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
212
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
213 // Return true if count is odd, false otherwise
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
214 return count&1; // Same as (count%2 == 1)
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
215 }
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
216
476
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
217 void RadiographyMaskLayer::DrawLine(const MaskPoint& start, const MaskPoint& end) const
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
218 {
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
219 int dx = (int)(end.x) - (int)(start.x);
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
220 int dy = (int)(end.y) - (int)(start.y);
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
221
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
222 if (std::abs(dx) > std::abs(dy))
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
223 { // the line is closer to horizontal
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
224
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
225 int incx = dx / std::abs(dx);
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
226 double incy = (double)(dy)/(double)(dx);
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
227 double y = (double)(start.y);
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
228 for (int x = (int)(start.x); x != (int)(end.x); x += incx, y += incy)
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
229 {
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
230 unsigned char* p = reinterpret_cast<unsigned char*>(mask_->GetRow((int)(y + 0.5))) + x;
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
231 *p = IN_MASK_VALUE;
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
232 }
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
233 }
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
234 else
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
235 { // the line is closer to vertical
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
236 int incy = dy / std::abs(dy);
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
237 double incx = (double)(dx)/(double)(dy);
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
238 double x = (double)(start.x);
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
239 for (int y = (int)(start.y); y != (int)(end.y); y += incy, x += incx)
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
240 {
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
241 unsigned char* p = reinterpret_cast<unsigned char*>(mask_->GetRow(y)) + (int)(x + 0.5);
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
242 *p = IN_MASK_VALUE;
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
243 }
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
244 }
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
245 }
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
246
475
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
247
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
248 void RadiographyMaskLayer::DrawMask() const
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
249 {
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
250 unsigned int left;
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
251 unsigned int right;
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
252 unsigned int top;
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
253 unsigned int bottom;
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
254
477
baf8c8d68bbc cleanup
am@osimis.io
parents: 476
diff changeset
255 ComputeMaskExtent(left, right, top, bottom, corners_);
475
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
256
476
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
257 // first fill the complete image
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
258 Orthanc::ImageProcessing::Set(*mask_, OUT_MASK_VALUE);
475
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
259
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
260 {
476
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
261 // from http://alienryderflex.com/polygon_fill/
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
262 std::auto_ptr<int> raiiNodeX(new int(corners_.size()));
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
263
478
017a66a2999b cleanup + avoid conversions
am@osimis.io
parents: 477
diff changeset
264 // convert all control points to double only once
017a66a2999b cleanup + avoid conversions
am@osimis.io
parents: 477
diff changeset
265 std::vector<double> cpx;
017a66a2999b cleanup + avoid conversions
am@osimis.io
parents: 477
diff changeset
266 std::vector<double> cpy;
017a66a2999b cleanup + avoid conversions
am@osimis.io
parents: 477
diff changeset
267 int cpSize = corners_.size();
017a66a2999b cleanup + avoid conversions
am@osimis.io
parents: 477
diff changeset
268 for (size_t i = 0; i < corners_.size(); i++)
017a66a2999b cleanup + avoid conversions
am@osimis.io
parents: 477
diff changeset
269 {
017a66a2999b cleanup + avoid conversions
am@osimis.io
parents: 477
diff changeset
270 cpx.push_back((double)corners_[i].x);
017a66a2999b cleanup + avoid conversions
am@osimis.io
parents: 477
diff changeset
271 cpy.push_back((double)corners_[i].y);
017a66a2999b cleanup + avoid conversions
am@osimis.io
parents: 477
diff changeset
272 }
017a66a2999b cleanup + avoid conversions
am@osimis.io
parents: 477
diff changeset
273
476
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
274 std::vector<int> nodeX;
478
017a66a2999b cleanup + avoid conversions
am@osimis.io
parents: 477
diff changeset
275 nodeX.resize(cpSize);
476
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
276 int nodes, pixelX, pixelY, i, j, swap ;
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
277
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
278 // Loop through the rows of the image.
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
279 for (pixelY = (int)top; pixelY < (int)bottom; pixelY++)
475
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
280 {
478
017a66a2999b cleanup + avoid conversions
am@osimis.io
parents: 477
diff changeset
281 double y = (double)pixelY;
476
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
282 // Build a list of nodes.
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
283 nodes = 0;
478
017a66a2999b cleanup + avoid conversions
am@osimis.io
parents: 477
diff changeset
284 j = cpSize - 1;
476
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
285
478
017a66a2999b cleanup + avoid conversions
am@osimis.io
parents: 477
diff changeset
286 for (i = 0; i < cpSize; i++)
476
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
287 {
478
017a66a2999b cleanup + avoid conversions
am@osimis.io
parents: 477
diff changeset
288 if ((cpy[i] <= y && cpy[j] >= y)
017a66a2999b cleanup + avoid conversions
am@osimis.io
parents: 477
diff changeset
289 || (cpy[j] <= y && cpy[i] >= y))
476
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
290 {
478
017a66a2999b cleanup + avoid conversions
am@osimis.io
parents: 477
diff changeset
291 nodeX[nodes++]= (int)(cpx[i] + (y - cpy[i])/(cpy[j] - cpy[i]) *(cpx[j] - cpx[i]));
476
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
292 }
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
293 j=i;
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
294 }
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
295
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
296 // Sort the nodes, via a simple “Bubble” sort.
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
297 i=0;
477
baf8c8d68bbc cleanup
am@osimis.io
parents: 476
diff changeset
298 while (i < nodes-1)
475
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
299 {
477
baf8c8d68bbc cleanup
am@osimis.io
parents: 476
diff changeset
300 if (nodeX[i] > nodeX[i+1])
476
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
301 {
477
baf8c8d68bbc cleanup
am@osimis.io
parents: 476
diff changeset
302 swap = nodeX[i];
baf8c8d68bbc cleanup
am@osimis.io
parents: 476
diff changeset
303 nodeX[i] = nodeX[i+1];
baf8c8d68bbc cleanup
am@osimis.io
parents: 476
diff changeset
304 nodeX[i+1] = swap;
baf8c8d68bbc cleanup
am@osimis.io
parents: 476
diff changeset
305 if (i)
baf8c8d68bbc cleanup
am@osimis.io
parents: 476
diff changeset
306 i--;
476
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
307 }
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
308 else
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
309 {
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
310 i++;
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
311 }
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
312 }
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
313
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
314 unsigned char* row = reinterpret_cast<unsigned char*>(mask_->GetRow(pixelY));
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
315 // Fill the pixels between node pairs.
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
316 for (i=0; i<nodes; i+=2)
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
317 {
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
318 if (nodeX[i ]>=(int)right)
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
319 break;
477
baf8c8d68bbc cleanup
am@osimis.io
parents: 476
diff changeset
320 if (nodeX[i+1]>= (int)left)
476
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
321 {
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
322 if (nodeX[i ]< (int)left )
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
323 nodeX[i ]=(int)left ;
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
324 if (nodeX[i+1]> (int)right)
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
325 nodeX[i+1]=(int)right;
477
baf8c8d68bbc cleanup
am@osimis.io
parents: 476
diff changeset
326 for (pixelX = nodeX[i]; pixelX <= nodeX[i+1]; pixelX++)
476
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
327 {
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
328 *(row + pixelX) = IN_MASK_VALUE;
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
329 }
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
330 }
475
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
331 }
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
332 }
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
333 }
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
334
476
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
335 // // draw lines
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
336 // for (size_t i = 1; i < corners_.size(); i++)
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
337 // {
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
338 // DrawLine(corners_[i-1], corners_[i]);
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
339 // }
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
340 // DrawLine(corners_[corners_.size()-1], corners_[0]);
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
341
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
342 // // fill between lines
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
343 // MaskPoint p(left, top);
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
344 // for (p.y = top; p.y <= bottom; p.y++)
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
345 // {
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
346 // unsigned char* q = reinterpret_cast<unsigned char*>(mask_->GetRow(p.y)) + left;
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
347 // unsigned char previousPixelValue1 = OUT_MASK_VALUE;
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
348 // unsigned char previousPixelValue2 = OUT_MASK_VALUE;
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
349 // for (p.x = left; p.x <= right; p.x++, q++)
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
350 // {
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
351 // if (*p == OUT_MASK_VALUE && previousPixelValue1 == IN_MASK_VALUE && previousPixelValue2 == OUT_MASK_VALUE) // we just passed over a single one pixel line => start filling
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
352 // {
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
353
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
354 // *q = IN_MASK_VALUE;
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
355 // }
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
356 // }
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
357 // }
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
358
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
359
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
360 // MaskPoint p(left, top);
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
361 // for (p.y = top; p.y <= bottom; p.y++)
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
362 // {
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
363 // unsigned char* q = reinterpret_cast<unsigned char*>(mask_->GetRow(p.y)) + left;
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
364 // for (p.x = left; p.x <= right; p.x++, q++)
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
365 // {
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
366 // if (isInside(corners_, p))
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
367 // {
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
368 // *q = IN_MASK_VALUE;
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
369 // }
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
370 // }
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
371 // }
a95090305dd4 Introduced ControlPoint instead of Corner in the trackers and layers + drawing mask from the ControlPoints
am@osimis.io
parents: 475
diff changeset
372
475
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
373 }
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
374
3c28542229a3 added a mask layer in the RadiographyWidget (to be cleaned)
am@osimis.io
parents:
diff changeset
375 }