Mercurial > hg > orthanc-stone

comparison Framework/Volumes/DicomStructureSetSlicer2.cpp @ 998:38b6bb0bdd72

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
added a new set of classes that correctly handle non-convex polygons (not used yet because of limitations in coordinates computing): DicomStructure2, DicomStructureSet2, DicomStructurePolygon2, DicomStructureSetSlicer2. Too many shortcuts have been taken when computing the actual position.
author Benjamin Golinvaux <bgo@osimis.io>
date Fri, 20 Sep 2019 11:58:00 +0200
parents
children 29f5f2031310
comparison
equal deleted inserted replaced
995:9893fa8cd7a6 998:38b6bb0bdd72
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-2019 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 #include "DicomStructureSetSlicer2.h"
22
23 #include "../Toolbox/GeometryToolbox.h"
24 #include "../Volumes/IVolumeSlicer.h"
25 #include "../Scene2D/PolylineSceneLayer.h"
26
27 namespace OrthancStone
28 {
29 DicomStructureSetSlicer2::DicomStructureSetSlicer2(boost::shared_ptr<DicomStructureSet2> structureSet)
30 : structureSet_(structureSet)
31 {}
32
33 IVolumeSlicer::IExtractedSlice* DicomStructureSetSlicer2::ExtractSlice(const CoordinateSystem3D& cuttingPlane)
34 {
35 // revision is always the same, hence 0
36 return new DicomStructureSetSlice2(structureSet_, 0, cuttingPlane);
37 }
38
39 DicomStructureSetSlice2::DicomStructureSetSlice2(
40 boost::weak_ptr<DicomStructureSet2> structureSet,
41 uint64_t revision,
42 const CoordinateSystem3D& cuttingPlane)
43 : structureSet_(structureSet.lock())
44 , isValid_(false)
45 {
46 bool opposite = false;
47
48 if (structureSet_->GetStructuresCount() == 0)
49 {
50 isValid_ = false;
51 }
52 else
53 {
54 // some structures seen in real life have no polygons. We must be
55 // careful
56 bool found = false;
57 size_t curStructure = 0;
58 while (!found && curStructure < structureSet_->GetStructuresCount())
59 {
60 if (structureSet_->GetStructure(curStructure).IsValid())
61 {
62 found = true;
63 const Vector normal = structureSet_->GetStructure(0).GetNormal();
64 isValid_ = (
65 GeometryToolbox::IsParallelOrOpposite(opposite, normal, cuttingPlane.GetNormal()) ||
66 GeometryToolbox::IsParallelOrOpposite(opposite, normal, cuttingPlane.GetAxisX()) ||
67 GeometryToolbox::IsParallelOrOpposite(opposite, normal, cuttingPlane.GetAxisY()));
68 }
69 }
70 }
71 }
72
73 ISceneLayer* DicomStructureSetSlice2::CreateSceneLayer(
74 const ILayerStyleConfigurator* configurator,
75 const CoordinateSystem3D& cuttingPlane)
76 {
77 assert(isValid_);
78
79 std::auto_ptr<PolylineSceneLayer> layer(new PolylineSceneLayer);
80 layer->SetThickness(2); // thickness of the on-screen line
81
82 for (size_t i = 0; i < structureSet_->GetStructuresCount(); i++)
83 {
84 const DicomStructure2& structure = structureSet_->GetStructure(i);
85 if (structure.IsValid())
86 {
87 const Color& color = structure.GetColor();
88
89 std::vector< std::pair<Point2D, Point2D> > segments;
90
91 if (structure.Project(segments, cuttingPlane))
92 {
93 for (size_t j = 0; j < segments.size(); j++)
94 {
95 PolylineSceneLayer::Chain chain;
96 chain.resize(2);
97
98 chain[0] = ScenePoint2D(segments[j].first.x, segments[j].first.y);
99 chain[1] = ScenePoint2D(segments[j].second.x, segments[j].second.y);
100
101 layer->AddChain(chain, false /* NOT closed */, color);
102 }
103 }
104 }
105 }
106 return layer.release();
107 }
108 }
109