Mercurial > hg > orthanc-stl
annotate Sources/Plugin.cpp @ 27:bd269cfb0da7
update year to 2024
author | Sebastien Jodogne <s.jodogne@gmail.com> |
---|---|
date | Wed, 03 Apr 2024 10:17:16 +0200 |
parents | 5fc3ea6acbd5 |
children | 410003c50b17 |
rev | line source |
---|---|
0
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1 /** |
24 | 2 * SPDX-FileCopyrightText: 2023-2024 Sebastien Jodogne, UCLouvain, Belgium |
0
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
3 * SPDX-License-Identifier: GPL-3.0-or-later |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
4 */ |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
5 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
6 /** |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
7 * STL plugin for Orthanc |
24 | 8 * Copyright (C) 2023-2024 Sebastien Jodogne, UCLouvain, Belgium |
0
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
9 * |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
10 * This program is free software: you can redistribute it and/or |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
11 * modify it under the terms of the GNU General Public License as |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
12 * published by the Free Software Foundation, either version 3 of the |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
13 * License, or (at your option) any later version. |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
14 * |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
15 * This program is distributed in the hope that it will be useful, but |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
16 * WITHOUT ANY WARRANTY; without even the implied warranty of |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
17 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
18 * General Public License for more details. |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
19 * |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
20 * You should have received a copy of the GNU General Public License |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
21 * along with this program. If not, see <http://www.gnu.org/licenses/>. |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
22 **/ |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
23 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
24 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
25 #include "../Resources/Orthanc/Plugins/OrthancPluginCppWrapper.h" |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
26 |
4
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
27 #include <EmbeddedResources.h> |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
28 |
6 | 29 #include <Compression/GzipCompressor.h> |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
30 #include <ChunkedBuffer.h> |
4
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
31 #include <DicomParsing/FromDcmtkBridge.h> |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
32 #include <DicomParsing/ParsedDicomFile.h> |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
33 #include <Images/ImageProcessing.h> |
0
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
34 #include <Logging.h> |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
35 #include <OrthancFramework.h> |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
36 #include <SerializationToolbox.h> |
0
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
37 #include <SystemToolbox.h> |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
38 |
4
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
39 #include <vtkImageConstantPad.h> |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
40 #include <vtkImageData.h> |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
41 #include <vtkImageResize.h> |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
42 #include <vtkMarchingCubes.h> |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
43 #include <vtkNew.h> |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
44 #include <vtkPolyData.h> |
4
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
45 #include <vtkPolyDataNormals.h> |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
46 #include <vtkSmoothPolyDataFilter.h> |
4
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
47 #include <vtkTriangle.h> |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
48 |
0
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
49 #include <boost/thread/shared_mutex.hpp> |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
50 |
6 | 51 #include <nifti1_io.h> |
52 | |
25
5fc3ea6acbd5
fix deprecated calls
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
24
diff
changeset
|
53 #define ORTHANC_PLUGIN_NAME "stl" |
5fc3ea6acbd5
fix deprecated calls
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
24
diff
changeset
|
54 |
6 | 55 |
0
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
56 // Forward declaration |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
57 void ReadStaticAsset(std::string& target, |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
58 const std::string& path); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
59 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
60 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
61 /** |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
62 * As the Three.js static assets are gzipped by the |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
63 * "EmbedStaticAssets.py" script, we use a cache to maintain the |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
64 * uncompressed assets in order to avoid multiple gzip decodings. |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
65 **/ |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
66 class ResourcesCache : public boost::noncopyable |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
67 { |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
68 private: |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
69 typedef std::map<std::string, std::string*> Content; |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
70 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
71 boost::shared_mutex mutex_; |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
72 Content content_; |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
73 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
74 public: |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
75 ~ResourcesCache() |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
76 { |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
77 for (Content::iterator it = content_.begin(); it != content_.end(); ++it) |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
78 { |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
79 assert(it->second != NULL); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
80 delete it->second; |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
81 } |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
82 } |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
83 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
84 void Answer(OrthancPluginRestOutput* output, |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
85 const std::string& path) |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
86 { |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
87 const std::string mime = Orthanc::EnumerationToString(Orthanc::SystemToolbox::AutodetectMimeType(path)); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
88 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
89 { |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
90 // Check whether the cache already contains the resource |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
91 boost::shared_lock<boost::shared_mutex> lock(mutex_); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
92 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
93 Content::const_iterator found = content_.find(path); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
94 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
95 if (found != content_.end()) |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
96 { |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
97 assert(found->second != NULL); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
98 OrthancPluginAnswerBuffer(OrthancPlugins::GetGlobalContext(), output, found->second->c_str(), found->second->size(), mime.c_str()); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
99 return; |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
100 } |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
101 } |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
102 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
103 // This resource has not been cached yet |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
104 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
105 std::unique_ptr<std::string> item(new std::string); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
106 ReadStaticAsset(*item, path); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
107 OrthancPluginAnswerBuffer(OrthancPlugins::GetGlobalContext(), output, item->c_str(), item->size(), mime.c_str()); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
108 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
109 { |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
110 // Store the resource into the cache |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
111 boost::unique_lock<boost::shared_mutex> lock(mutex_); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
112 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
113 if (content_.find(path) == content_.end()) |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
114 { |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
115 content_[path] = item.release(); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
116 } |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
117 } |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
118 } |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
119 }; |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
120 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
121 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
122 static ResourcesCache cache_; |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
123 static bool hasCreateDicomStl_; |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
124 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
125 void ServeFile(OrthancPluginRestOutput* output, |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
126 const char* url, |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
127 const OrthancPluginHttpRequest* request) |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
128 { |
3
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
129 if (request->method != OrthancPluginHttpMethod_Get) |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
130 { |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
131 OrthancPluginSendMethodNotAllowed(OrthancPlugins::GetGlobalContext(), output, "GET"); |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
132 return; |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
133 } |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
134 |
0
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
135 std::string file = request->groups[0]; |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
136 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
137 if (file == "viewer.html") |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
138 { |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
139 std::string s; |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
140 Orthanc::EmbeddedResources::GetFileResource(s, Orthanc::EmbeddedResources::VIEWER_HTML); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
141 OrthancPluginAnswerBuffer(OrthancPlugins::GetGlobalContext(), output, s.c_str(), s.size(), Orthanc::EnumerationToString(Orthanc::MimeType_Html)); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
142 } |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
143 else if (file == "viewer.js") |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
144 { |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
145 std::string s; |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
146 Orthanc::EmbeddedResources::GetFileResource(s, Orthanc::EmbeddedResources::VIEWER_JS); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
147 OrthancPluginAnswerBuffer(OrthancPlugins::GetGlobalContext(), output, s.c_str(), s.size(), Orthanc::EnumerationToString(Orthanc::MimeType_JavaScript)); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
148 } |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
149 else |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
150 { |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
151 cache_.Answer(output, file); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
152 } |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
153 } |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
154 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
155 |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
156 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
157 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
158 #include <dcmtk/dcmdata/dcdeftag.h> |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
159 #include <dcmtk/dcmdata/dcfilefo.h> |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
160 #include <dcmtk/dcmdata/dcitem.h> |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
161 #include <dcmtk/dcmdata/dcsequen.h> |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
162 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
163 class Extent2D : public boost::noncopyable |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
164 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
165 private: |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
166 bool isEmpty_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
167 double x1_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
168 double y1_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
169 double x2_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
170 double y2_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
171 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
172 void CheckNotEmpty() const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
173 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
174 if (isEmpty_) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
175 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
176 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
177 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
178 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
179 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
180 public: |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
181 Extent2D() : |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
182 isEmpty_(true), |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
183 x1_(0), |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
184 y1_(0), |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
185 x2_(0), |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
186 y2_(0) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
187 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
188 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
189 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
190 bool IsEmpty() const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
191 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
192 return isEmpty_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
193 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
194 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
195 double GetMinX() const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
196 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
197 CheckNotEmpty(); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
198 return x1_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
199 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
200 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
201 double GetMaxX() const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
202 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
203 CheckNotEmpty(); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
204 return x2_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
205 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
206 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
207 double GetMinY() const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
208 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
209 CheckNotEmpty(); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
210 return y1_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
211 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
212 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
213 double GetMaxY() const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
214 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
215 CheckNotEmpty(); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
216 return y2_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
217 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
218 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
219 double GetWidth() const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
220 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
221 CheckNotEmpty(); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
222 return x2_ - x1_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
223 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
224 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
225 double GetHeight() const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
226 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
227 CheckNotEmpty(); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
228 return y2_ - y1_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
229 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
230 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
231 void Add(double x, |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
232 double y) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
233 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
234 if (isEmpty_) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
235 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
236 x1_ = x2_ = x; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
237 y1_ = y2_ = y; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
238 isEmpty_ = false; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
239 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
240 else |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
241 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
242 x1_ = std::min(x1_, x); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
243 x2_ = std::max(x2_, x); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
244 y1_ = std::min(y1_, y); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
245 y2_ = std::max(y2_, y); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
246 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
247 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
248 }; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
249 |
3
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
250 static std::string GetStringValue(DcmItem& item, |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
251 const DcmTagKey& key) |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
252 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
253 const char* s = NULL; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
254 if (!item.findAndGetString(key, s).good() || |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
255 s == NULL) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
256 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
257 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
258 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
259 else |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
260 { |
3
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
261 return Orthanc::Toolbox::StripSpaces(s); |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
262 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
263 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
264 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
265 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
266 static void ListStructuresNames(std::set<std::string>& target, |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
267 Orthanc::ParsedDicomFile& source) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
268 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
269 target.clear(); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
270 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
271 DcmSequenceOfItems* sequence = NULL; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
272 if (!source.GetDcmtkObject().getDataset()->findAndGetSequence(DCM_StructureSetROISequence, sequence).good() || |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
273 sequence == NULL) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
274 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
275 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
276 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
277 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
278 for (unsigned long i = 0; i < sequence->card(); i++) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
279 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
280 DcmItem* item = sequence->getItem(i); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
281 if (item == NULL) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
282 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
283 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
284 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
285 else |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
286 { |
3
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
287 target.insert(GetStringValue(*item, DCM_ROIName)); |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
288 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
289 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
290 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
291 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
292 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
293 static bool IsNear(double a, |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
294 double b) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
295 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
296 return std::abs(a - b) < 10.0 * std::numeric_limits<double>::epsilon(); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
297 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
298 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
299 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
300 class Vector3D |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
301 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
302 private: |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
303 double x_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
304 double y_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
305 double z_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
306 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
307 public: |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
308 Vector3D() : |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
309 x_(0), |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
310 y_(0), |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
311 z_(0) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
312 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
313 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
314 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
315 Vector3D(double x, |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
316 double y, |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
317 double z) : |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
318 x_(x), |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
319 y_(y), |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
320 z_(z) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
321 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
322 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
323 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
324 Vector3D(const Vector3D& from, |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
325 const Vector3D& to) : |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
326 x_(to.x_ - from.x_), |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
327 y_(to.y_ - from.y_), |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
328 z_(to.z_ - from.z_) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
329 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
330 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
331 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
332 double GetX() const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
333 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
334 return x_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
335 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
336 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
337 double GetY() const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
338 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
339 return y_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
340 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
341 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
342 double GetZ() const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
343 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
344 return z_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
345 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
346 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
347 double ComputeNorm() const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
348 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
349 return sqrt(x_ * x_ + y_ * y_ + z_ * z_); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
350 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
351 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
352 void Normalize() |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
353 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
354 double norm = ComputeNorm(); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
355 if (!IsNear(norm, 0)) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
356 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
357 x_ /= norm; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
358 y_ /= norm; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
359 z_ /= norm; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
360 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
361 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
362 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
363 static Vector3D CrossProduct(const Vector3D& u, |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
364 const Vector3D& v) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
365 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
366 return Vector3D(u.GetY() * v.GetZ() - u.GetZ() * v.GetY(), |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
367 u.GetZ() * v.GetX() - u.GetX() * v.GetZ(), |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
368 u.GetX() * v.GetY() - u.GetY() * v.GetX()); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
369 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
370 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
371 static double DotProduct(const Vector3D& a, |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
372 const Vector3D& b) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
373 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
374 return a.GetX() * b.GetX() + a.GetY() * b.GetY() + a.GetZ() * b.GetZ(); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
375 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
376 }; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
377 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
378 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
379 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
380 static bool MyParseDouble(double& value, |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
381 const std::string& s) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
382 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
383 #if 1 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
384 char* end = NULL; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
385 value = strtod(s.c_str(), &end); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
386 return (end == s.c_str() + s.size()); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
387 #else |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
388 return Orthanc::SerializationToolbox::ParseDouble(value, s); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
389 #endif |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
390 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
391 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
392 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
393 class StructurePolygon : public boost::noncopyable |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
394 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
395 private: |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
396 std::string roiName_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
397 std::string referencedSopInstanceUid_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
398 uint8_t red_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
399 uint8_t green_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
400 uint8_t blue_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
401 std::vector<Vector3D> points_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
402 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
403 public: |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
404 StructurePolygon(Orthanc::ParsedDicomFile& dicom, |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
405 unsigned long roiIndex, |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
406 unsigned long contourIndex) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
407 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
408 DcmDataset& dataset = *dicom.GetDcmtkObject().getDataset(); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
409 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
410 DcmItem* structure = NULL; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
411 DcmItem* roi = NULL; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
412 DcmItem* contour = NULL; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
413 DcmSequenceOfItems* referenced = NULL; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
414 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
415 if (!dataset.findAndGetSequenceItem(DCM_StructureSetROISequence, structure, roiIndex).good() || |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
416 structure == NULL || |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
417 !dataset.findAndGetSequenceItem(DCM_ROIContourSequence, roi, roiIndex).good() || |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
418 roi == NULL || |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
419 !roi->findAndGetSequenceItem(DCM_ContourSequence, contour, contourIndex).good() || |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
420 contour == NULL || |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
421 !contour->findAndGetSequence(DCM_ContourImageSequence, referenced).good() || |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
422 referenced == NULL || |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
423 referenced->card() != 1) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
424 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
425 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
426 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
427 |
3
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
428 roiName_ = GetStringValue(*structure, DCM_ROIName); |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
429 referencedSopInstanceUid_ = GetStringValue(*referenced->getItem(0), DCM_ReferencedSOPInstanceUID); |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
430 |
3
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
431 if (GetStringValue(*contour, DCM_ContourGeometricType) != "CLOSED_PLANAR") |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
432 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
433 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
434 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
435 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
436 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
437 std::vector<std::string> tokens; |
3
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
438 Orthanc::Toolbox::TokenizeString(tokens, GetStringValue(*roi, DCM_ROIDisplayColor), '\\'); |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
439 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
440 uint32_t r, g, b; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
441 if (tokens.size() != 3 || |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
442 !Orthanc::SerializationToolbox::ParseFirstUnsignedInteger32(r, tokens[0]) || |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
443 !Orthanc::SerializationToolbox::ParseFirstUnsignedInteger32(g, tokens[1]) || |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
444 !Orthanc::SerializationToolbox::ParseFirstUnsignedInteger32(b, tokens[2]) || |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
445 r > 255 || |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
446 g > 255 || |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
447 b > 255) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
448 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
449 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
450 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
451 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
452 red_ = r; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
453 green_ = g; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
454 blue_ = b; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
455 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
456 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
457 { |
3
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
458 std::vector<std::string> tokens; |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
459 Orthanc::Toolbox::TokenizeString(tokens, GetStringValue(*contour, DCM_ContourData), '\\'); |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
460 |
3
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
461 const std::string s = GetStringValue(*contour, DCM_NumberOfContourPoints); |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
462 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
463 uint32_t countPoints; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
464 if (!Orthanc::SerializationToolbox::ParseUnsignedInteger32(countPoints, s) || |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
465 tokens.size() != 3 * countPoints) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
466 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
467 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
468 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
469 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
470 points_.reserve(countPoints); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
471 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
472 for (size_t i = 0; i < tokens.size(); i += 3) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
473 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
474 double x, y, z; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
475 if (!MyParseDouble(x, tokens[i]) || |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
476 !MyParseDouble(y, tokens[i + 1]) || |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
477 !MyParseDouble(z, tokens[i + 2])) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
478 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
479 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
480 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
481 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
482 points_.push_back(Vector3D(x, y, z)); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
483 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
484 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
485 assert(points_.size() == countPoints); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
486 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
487 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
488 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
489 const std::string& GetRoiName() const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
490 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
491 return roiName_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
492 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
493 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
494 const std::string& GetReferencedSopInstanceUid() const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
495 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
496 return referencedSopInstanceUid_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
497 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
498 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
499 size_t GetPointsCount() const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
500 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
501 return points_.size(); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
502 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
503 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
504 const Vector3D& GetPoint(size_t i) const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
505 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
506 if (i >= points_.size()) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
507 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
508 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
509 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
510 else |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
511 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
512 return points_[i]; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
513 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
514 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
515 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
516 bool IsCoplanar(Vector3D& normal) const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
517 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
518 if (points_.size() < 3) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
519 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
520 return false; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
521 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
522 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
523 bool hasNormal = false; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
524 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
525 for (size_t i = 0; i < points_.size(); i++) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
526 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
527 normal = Vector3D::CrossProduct(Vector3D(points_[1], points_[0]), |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
528 Vector3D(points_[2], points_[0])); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
529 if (!IsNear(normal.ComputeNorm(), 0)) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
530 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
531 normal.Normalize(); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
532 hasNormal = true; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
533 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
534 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
535 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
536 if (!hasNormal) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
537 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
538 return false; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
539 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
540 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
541 double a = Vector3D::DotProduct(points_[0], normal); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
542 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
543 for (size_t i = 1; i < points_.size(); i++) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
544 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
545 double b = Vector3D::DotProduct(points_[i], normal); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
546 if (!IsNear(a, b)) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
547 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
548 return false; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
549 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
550 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
551 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
552 return true; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
553 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
554 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
555 void Add(Extent2D& extent, |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
556 const Vector3D& axisX, |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
557 const Vector3D& axisY) const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
558 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
559 assert(IsNear(1, axisX.ComputeNorm())); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
560 assert(IsNear(1, axisY.ComputeNorm())); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
561 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
562 for (size_t i = 0; i < points_.size(); i++) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
563 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
564 extent.Add(Vector3D::DotProduct(axisX, points_[i]), |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
565 Vector3D::DotProduct(axisY, points_[i])); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
566 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
567 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
568 }; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
569 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
570 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
571 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
572 struct IsNearPredicate |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
573 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
574 bool operator() (const double& a, |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
575 const double& b) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
576 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
577 return IsNear(a, b); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
578 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
579 }; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
580 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
581 static void RemoveDuplicateValues(std::vector<double>& v) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
582 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
583 IsNearPredicate predicate; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
584 std::vector<double>::iterator last = std::unique(v.begin(), v.end(), predicate); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
585 v.erase(last, v.end()); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
586 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
587 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
588 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
589 class StructureSet : public boost::noncopyable |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
590 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
591 private: |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
592 std::vector<StructurePolygon*> polygons_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
593 bool hasGeometry_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
594 Vector3D slicesNormal_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
595 double slicesSpacing_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
596 double minProjectionAlongNormal_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
597 double maxProjectionAlongNormal_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
598 std::string patientId_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
599 std::string studyInstanceUid_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
600 std::string seriesInstanceUid_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
601 std::string sopInstanceUid_; |
4
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
602 bool hasFrameOfReferenceUid_; |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
603 std::string frameOfReferenceUid_; |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
604 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
605 void ComputeGeometry() |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
606 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
607 std::list<double> positionsList; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
608 hasGeometry_ = false; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
609 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
610 for (size_t i = 0; i < polygons_.size(); i++) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
611 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
612 assert(polygons_[i] != NULL); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
613 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
614 Vector3D n; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
615 if (polygons_[i]->IsCoplanar(n)) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
616 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
617 const Vector3D& point = polygons_[i]->GetPoint(0); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
618 double z = Vector3D::DotProduct(point, n); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
619 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
620 if (!hasGeometry_) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
621 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
622 hasGeometry_ = true; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
623 slicesNormal_ = n; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
624 minProjectionAlongNormal_ = z; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
625 maxProjectionAlongNormal_ = z; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
626 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
627 else if (!IsNear(std::abs(Vector3D::DotProduct(n, slicesNormal_)), 1)) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
628 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
629 hasGeometry_ = false; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
630 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
631 // RT-STRUCT with non-parallel slices |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
632 throw Orthanc::OrthancException(Orthanc::ErrorCode_NotImplemented); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
633 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
634 else |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
635 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
636 minProjectionAlongNormal_ = std::min(minProjectionAlongNormal_, z); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
637 maxProjectionAlongNormal_ = std::max(maxProjectionAlongNormal_, z); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
638 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
639 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
640 positionsList.push_back(Vector3D::DotProduct(n, point)); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
641 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
642 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
643 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
644 if (hasGeometry_) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
645 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
646 std::vector<double> positions(positionsList.begin(), positionsList.end()); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
647 assert(!positions.empty()); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
648 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
649 std::sort(positions.begin(), positions.end()); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
650 RemoveDuplicateValues(positions); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
651 assert(!positions.empty()); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
652 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
653 if (positions.size() == 1) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
654 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
655 hasGeometry_ = false; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
656 return; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
657 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
658 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
659 std::vector<double> offsets; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
660 offsets.resize(positions.size() - 1); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
661 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
662 for (size_t i = 0; i < offsets.size(); i++) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
663 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
664 offsets[i] = positions[i + 1] - positions[i]; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
665 assert(offsets[i] > 0); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
666 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
667 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
668 std::sort(offsets.begin(), offsets.end()); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
669 RemoveDuplicateValues(offsets); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
670 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
671 slicesSpacing_ = offsets[0]; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
672 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
673 for (size_t i = 1; i < offsets.size(); i++) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
674 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
675 double d = offsets[i] / slicesSpacing_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
676 if (!IsNear(d, round(d))) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
677 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
678 // Irregular spacing between the slices |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
679 hasGeometry_ = false; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
680 break; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
681 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
682 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
683 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
684 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
685 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
686 void CheckHasGeometry() const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
687 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
688 if (!hasGeometry_) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
689 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
690 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
691 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
692 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
693 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
694 public: |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
695 explicit StructureSet(Orthanc::ParsedDicomFile& dicom) : |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
696 hasGeometry_(false), |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
697 slicesSpacing_(0), |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
698 minProjectionAlongNormal_(0), |
4
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
699 maxProjectionAlongNormal_(0), |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
700 hasFrameOfReferenceUid_(false) |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
701 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
702 DcmDataset& dataset = *dicom.GetDcmtkObject().getDataset(); |
3
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
703 patientId_ = GetStringValue(dataset, DCM_PatientID); |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
704 studyInstanceUid_ = GetStringValue(dataset, DCM_StudyInstanceUID); |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
705 seriesInstanceUid_ = GetStringValue(dataset, DCM_SeriesInstanceUID); |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
706 sopInstanceUid_ = GetStringValue(dataset, DCM_SOPInstanceUID); |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
707 |
4
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
708 DcmSequenceOfItems* frame = NULL; |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
709 if (!dataset.findAndGetSequence(DCM_ReferencedFrameOfReferenceSequence, frame).good() || |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
710 frame == NULL) |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
711 { |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
712 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
713 } |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
714 |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
715 if (frame->card() == 1) |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
716 { |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
717 const char* v = NULL; |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
718 if (frame->getItem(0)->findAndGetString(DCM_FrameOfReferenceUID, v).good() && |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
719 v != NULL) |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
720 { |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
721 hasFrameOfReferenceUid_ = true; |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
722 frameOfReferenceUid_.assign(v); |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
723 } |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
724 } |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
725 |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
726 DcmSequenceOfItems* rois = NULL; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
727 if (!dataset.findAndGetSequence(DCM_ROIContourSequence, rois).good() || |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
728 rois == NULL) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
729 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
730 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
731 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
732 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
733 std::vector<DcmSequenceOfItems*> contours(rois->card()); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
734 size_t countPolygons = 0; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
735 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
736 for (unsigned long i = 0; i < rois->card(); i++) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
737 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
738 DcmSequenceOfItems* contour = NULL; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
739 if (!rois->getItem(i)->findAndGetSequence(DCM_ContourSequence, contour).good() || |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
740 contour == NULL) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
741 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
742 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
743 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
744 else |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
745 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
746 contours[i] = contour; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
747 countPolygons += contour->card(); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
748 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
749 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
750 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
751 polygons_.resize(countPolygons); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
752 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
753 size_t pos = 0; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
754 for (unsigned long i = 0; i < contours.size(); i++) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
755 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
756 for (unsigned long j = 0; j < contours[i]->card(); j++, pos++) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
757 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
758 polygons_[pos] = new StructurePolygon(dicom, i, j); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
759 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
760 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
761 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
762 assert(pos == countPolygons); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
763 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
764 ComputeGeometry(); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
765 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
766 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
767 ~StructureSet() |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
768 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
769 for (size_t i = 0; i < polygons_.size(); i++) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
770 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
771 assert(polygons_[i] != NULL); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
772 delete polygons_[i]; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
773 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
774 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
775 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
776 const std::string& GetPatient() const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
777 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
778 return patientId_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
779 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
780 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
781 const std::string& GetStudyInstanceUid() const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
782 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
783 return studyInstanceUid_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
784 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
785 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
786 const std::string& GetSeriesInstanceUid() const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
787 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
788 return seriesInstanceUid_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
789 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
790 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
791 const std::string& GetSopInstanceUid() const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
792 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
793 return sopInstanceUid_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
794 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
795 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
796 std::string HashStudy() const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
797 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
798 std::string s; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
799 Orthanc::Toolbox::ComputeSHA1(s, patientId_ + "|" + studyInstanceUid_); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
800 return s; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
801 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
802 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
803 size_t GetPolygonsCount() const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
804 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
805 return polygons_.size(); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
806 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
807 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
808 const StructurePolygon& GetPolygon(size_t i) const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
809 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
810 if (i >= polygons_.size()) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
811 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
812 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
813 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
814 else |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
815 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
816 assert(polygons_[i] != NULL); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
817 return *polygons_[i]; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
818 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
819 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
820 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
821 bool HasGeometry() const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
822 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
823 return hasGeometry_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
824 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
825 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
826 const Vector3D& GetSlicesNormal() const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
827 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
828 CheckHasGeometry(); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
829 return slicesNormal_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
830 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
831 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
832 double GetSlicesSpacing() const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
833 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
834 CheckHasGeometry(); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
835 return slicesSpacing_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
836 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
837 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
838 double GetMinProjectionAlongNormal() const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
839 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
840 CheckHasGeometry(); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
841 return minProjectionAlongNormal_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
842 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
843 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
844 double GetMaxProjectionAlongNormal() const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
845 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
846 CheckHasGeometry(); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
847 return maxProjectionAlongNormal_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
848 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
849 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
850 double ProjectAlongNormal(const StructurePolygon& polygon) const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
851 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
852 CheckHasGeometry(); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
853 return Vector3D::DotProduct(slicesNormal_, polygon.GetPoint(0)); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
854 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
855 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
856 size_t GetSlicesCount() const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
857 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
858 CheckHasGeometry(); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
859 double c = (maxProjectionAlongNormal_ - minProjectionAlongNormal_) / slicesSpacing_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
860 assert(c >= 0); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
861 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
862 if (!IsNear(c, round(c))) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
863 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
864 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
865 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
866 else |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
867 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
868 return static_cast<size_t>(round(c)) + 1; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
869 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
870 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
871 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
872 bool LookupSliceIndex(size_t& slice, |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
873 const StructurePolygon& polygon) const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
874 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
875 CheckHasGeometry(); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
876 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
877 double z = ProjectAlongNormal(polygon); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
878 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
879 if (z < minProjectionAlongNormal_ || |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
880 z > maxProjectionAlongNormal_) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
881 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
882 return false; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
883 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
884 else |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
885 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
886 double c = (z - minProjectionAlongNormal_) / slicesSpacing_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
887 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
888 if (IsNear(c, round(c))) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
889 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
890 slice = static_cast<size_t>(round(c)); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
891 return true; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
892 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
893 else |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
894 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
895 return false; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
896 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
897 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
898 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
899 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
900 bool LookupReferencedSopInstanceUid(std::string& sopInstanceUid) const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
901 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
902 if (HasGeometry()) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
903 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
904 for (size_t i = 0; i < polygons_.size(); i++) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
905 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
906 assert(polygons_[i] != NULL); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
907 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
908 Vector3D n; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
909 if (polygons_[i]->IsCoplanar(n) && |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
910 Vector3D::DotProduct(n, slicesNormal_)) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
911 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
912 sopInstanceUid = polygons_[i]->GetReferencedSopInstanceUid(); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
913 return true; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
914 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
915 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
916 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
917 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
918 return false; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
919 } |
4
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
920 |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
921 bool HasFrameOfReferenceUid() const |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
922 { |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
923 return hasFrameOfReferenceUid_; |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
924 } |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
925 |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
926 const std::string& GetFrameOfReferenceUid() const |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
927 { |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
928 if (hasFrameOfReferenceUid_) |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
929 { |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
930 return frameOfReferenceUid_; |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
931 } |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
932 else |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
933 { |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
934 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls); |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
935 } |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
936 } |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
937 }; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
938 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
939 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
940 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
941 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
942 class XorFiller : public Orthanc::ImageProcessing::IPolygonFiller |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
943 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
944 private: |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
945 Orthanc::ImageAccessor& target_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
946 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
947 public: |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
948 XorFiller(Orthanc::ImageAccessor& target) : |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
949 target_(target) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
950 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
951 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
952 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
953 virtual void Fill(int y, |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
954 int x1, |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
955 int x2) ORTHANC_OVERRIDE |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
956 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
957 assert(x1 <= x2); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
958 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
959 if (y >= 0 && |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
960 y < static_cast<int>(target_.GetHeight())) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
961 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
962 x1 = std::max(x1, 0); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
963 x2 = std::min(x2, static_cast<int>(target_.GetWidth()) - 1); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
964 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
965 uint8_t* p = reinterpret_cast<uint8_t*>(target_.GetRow(y)) + x1; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
966 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
967 for (int i = x1; i <= x2; i++, p++) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
968 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
969 *p = (*p ^ 0xff); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
970 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
971 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
972 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
973 }; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
974 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
975 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
976 static void EncodeSTL(std::string& target /* out */, |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
977 vtkPolyData& mesh /* in */) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
978 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
979 // TODO - Conversion to little endian on big endian |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
980 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
981 Orthanc::ChunkedBuffer buffer; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
982 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
983 uint8_t header[80]; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
984 memset(header, 0, sizeof(header)); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
985 buffer.AddChunk(header, sizeof(header)); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
986 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
987 uint32_t n = mesh.GetNumberOfCells(); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
988 buffer.AddChunk(&n, sizeof(n)); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
989 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
990 for (vtkIdType i = 0; i < mesh.GetNumberOfCells(); i++) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
991 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
992 vtkCell* cell = mesh.GetCell(i); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
993 vtkTriangle* triangle = dynamic_cast<vtkTriangle*>(cell); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
994 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
995 double p0[3]; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
996 double p1[3]; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
997 double p2[3]; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
998 triangle->GetPoints()->GetPoint(0, p0); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
999 triangle->GetPoints()->GetPoint(1, p1); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1000 triangle->GetPoints()->GetPoint(2, p2); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1001 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1002 double normal[3]; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1003 vtkTriangle::ComputeNormal(p0, p1, p2, normal); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1004 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1005 float d[4 * 3] = { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1006 static_cast<float>(normal[0]), static_cast<float>(normal[1]), static_cast<float>(normal[2]), |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1007 static_cast<float>(p0[0]), static_cast<float>(p0[1]), static_cast<float>(p0[2]), |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1008 static_cast<float>(p1[0]), static_cast<float>(p1[1]), static_cast<float>(p1[2]), |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1009 static_cast<float>(p2[0]), static_cast<float>(p2[1]), static_cast<float>(p2[2]) }; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1010 buffer.AddChunk(d, sizeof(d)); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1011 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1012 uint16_t a = 0; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1013 buffer.AddChunk(&a, sizeof(a)); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1014 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1015 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1016 buffer.Flatten(target); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1017 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1018 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1019 |
8
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1020 bool EncodeVolume(std::string& stl, |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1021 vtkImageData* volume, |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1022 unsigned int resolution, |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1023 bool smooth) |
7 | 1024 { |
1025 if (volume == NULL) | |
1026 { | |
1027 throw Orthanc::OrthancException(Orthanc::ErrorCode_NullPointer); | |
1028 } | |
1029 | |
1030 vtkNew<vtkImageResize> resize; | |
1031 resize->SetOutputDimensions(resolution, resolution, resolution); | |
1032 resize->SetInputData(volume); | |
1033 resize->Update(); | |
1034 | |
1035 vtkNew<vtkImageConstantPad> padding; | |
1036 padding->SetConstant(0); | |
1037 padding->SetOutputNumberOfScalarComponents(1); | |
1038 padding->SetOutputWholeExtent(-1, resolution, -1, resolution, -1, resolution); | |
1039 padding->SetInputData(resize->GetOutput()); | |
1040 padding->Update(); | |
1041 | |
1042 double range[2]; | |
1043 padding->GetOutput()->GetScalarRange(range); | |
1044 | |
1045 const double isoValue = (range[0] + range[1]) / 2.0; | |
1046 | |
1047 vtkNew<vtkMarchingCubes> surface; | |
1048 surface->SetInputData(padding->GetOutput()); | |
1049 surface->ComputeNormalsOn(); | |
1050 surface->SetValue(0, isoValue); | |
1051 surface->Update(); | |
1052 | |
1053 if (smooth) | |
1054 { | |
1055 vtkNew<vtkSmoothPolyDataFilter> smoothFilter; | |
1056 // Apply volume smoothing | |
1057 // https://examples.vtk.org/site/Cxx/PolyData/SmoothPolyDataFilter/ | |
1058 smoothFilter->SetInputConnection(surface->GetOutputPort()); | |
1059 smoothFilter->SetNumberOfIterations(15); | |
1060 smoothFilter->SetRelaxationFactor(0.1); | |
1061 smoothFilter->FeatureEdgeSmoothingOff(); | |
1062 smoothFilter->BoundarySmoothingOn(); | |
1063 smoothFilter->Update(); | |
1064 | |
1065 vtkNew<vtkPolyDataNormals> normalGenerator; | |
1066 normalGenerator->SetInputConnection(smoothFilter->GetOutputPort()); | |
1067 normalGenerator->ComputePointNormalsOn(); | |
1068 normalGenerator->ComputeCellNormalsOn(); | |
1069 normalGenerator->Update(); | |
1070 | |
1071 EncodeSTL(stl, *normalGenerator->GetOutput()); | |
1072 } | |
1073 else | |
1074 { | |
1075 EncodeSTL(stl, *surface->GetOutput()); | |
1076 } | |
1077 | |
1078 return true; | |
1079 } | |
1080 | |
1081 | |
1082 bool EncodeStructureSetMesh(std::string& stl, | |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1083 const StructureSet& structureSet, |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1084 const std::set<std::string>& roiNames, |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1085 unsigned int resolution, |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1086 bool smooth) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1087 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1088 if (!structureSet.HasGeometry()) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1089 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1090 return false; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1091 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1092 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1093 if (resolution < 1) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1094 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1095 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1096 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1097 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1098 if (!IsNear(1, structureSet.GetSlicesNormal().ComputeNorm())) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1099 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1100 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1101 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1102 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1103 // TODO - Axes could be retrieved from the referenced CT volume |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1104 Vector3D axisX(1, 0, 0); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1105 Vector3D axisY = Vector3D::CrossProduct(structureSet.GetSlicesNormal(), axisX); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1106 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1107 if (!IsNear(1, axisX.ComputeNorm()) || |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1108 !IsNear(1, axisY.ComputeNorm())) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1109 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1110 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1111 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1112 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1113 Extent2D extent; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1114 for (size_t i = 0; i < structureSet.GetPolygonsCount(); i++) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1115 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1116 structureSet.GetPolygon(i).Add(extent, axisX, axisY); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1117 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1118 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1119 const int depth = structureSet.GetSlicesCount(); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1120 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1121 vtkNew<vtkImageData> volume; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1122 volume->SetDimensions(resolution, resolution, depth); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1123 volume->AllocateScalars(VTK_UNSIGNED_CHAR, 1); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1124 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1125 assert(sizeof(unsigned char) == 1); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1126 memset(volume->GetScalarPointer(), 0, resolution * resolution * depth); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1127 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1128 for (size_t i = 0; i < structureSet.GetPolygonsCount(); i++) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1129 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1130 const StructurePolygon& polygon = structureSet.GetPolygon(i); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1131 if (roiNames.find(polygon.GetRoiName()) == roiNames.end()) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1132 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1133 // This polygon doesn't correspond to a ROI of interest |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1134 continue; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1135 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1136 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1137 size_t j; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1138 if (!structureSet.LookupSliceIndex(j, polygon)) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1139 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1140 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1141 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1142 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1143 std::vector<Orthanc::ImageProcessing::ImagePoint> points; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1144 points.reserve(polygon.GetPointsCount()); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1145 for (size_t j = 0; j < polygon.GetPointsCount(); j++) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1146 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1147 const Vector3D& point = polygon.GetPoint(j); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1148 double x = (Vector3D::DotProduct(point, axisX) - extent.GetMinX()) / extent.GetWidth() * static_cast<double>(resolution); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1149 double y = (Vector3D::DotProduct(point, axisY) - extent.GetMinY()) / extent.GetHeight() * static_cast<double>(resolution); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1150 points.push_back(Orthanc::ImageProcessing::ImagePoint(static_cast<int32_t>(std::floor(x)), |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1151 static_cast<int32_t>(std::floor(y)))); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1152 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1153 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1154 Orthanc::ImageAccessor slice; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1155 slice.AssignWritable(Orthanc::PixelFormat_Grayscale8, resolution, resolution, resolution /* pitch */, |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1156 reinterpret_cast<uint8_t*>(volume->GetScalarPointer()) + j * resolution * resolution); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1157 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1158 XorFiller filler(slice); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1159 Orthanc::ImageProcessing::FillPolygon(filler, points); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1160 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1161 |
7 | 1162 volume->SetSpacing( |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1163 extent.GetWidth() / static_cast<double>(resolution), |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1164 extent.GetHeight() / static_cast<double>(resolution), |
7 | 1165 (structureSet.GetMaxProjectionAlongNormal() - structureSet.GetMinProjectionAlongNormal()) / static_cast<double>(depth)); |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1166 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1167 // TODO |
7 | 1168 // volume->SetOrigin() |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1169 |
8
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1170 return EncodeVolume(stl, volume.Get(), resolution, smooth); |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1171 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1172 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1173 |
3
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1174 static Orthanc::ParsedDicomFile* LoadInstance(const std::string& instanceId) |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1175 { |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1176 std::string dicom; |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1177 |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1178 if (!OrthancPlugins::RestApiGetString(dicom, "/instances/" + instanceId + "/file", false)) |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1179 { |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1180 throw Orthanc::OrthancException(Orthanc::ErrorCode_UnknownResource); |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1181 } |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1182 else |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1183 { |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1184 return new Orthanc::ParsedDicomFile(dicom); |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1185 } |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1186 } |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1187 |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1188 |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1189 void ListStructures(OrthancPluginRestOutput* output, |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1190 const char* url, |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1191 const OrthancPluginHttpRequest* request) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1192 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1193 if (request->method != OrthancPluginHttpMethod_Get) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1194 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1195 OrthancPluginSendMethodNotAllowed(OrthancPlugins::GetGlobalContext(), output, "GET"); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1196 return; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1197 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1198 |
3
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1199 const std::string instanceId(request->groups[0]); |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1200 |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1201 std::unique_ptr<Orthanc::ParsedDicomFile> dicom(LoadInstance(instanceId)); |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1202 |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1203 std::set<std::string> names; |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1204 ListStructuresNames(names, *dicom); |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1205 |
3
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1206 Json::Value answer = Json::arrayValue; |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1207 |
3
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1208 for (std::set<std::string>::const_iterator it = names.begin(); it != names.end(); ++it) |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1209 { |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1210 answer.append(*it); |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1211 } |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1212 |
3
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1213 std::string s = answer.toStyledString(); |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1214 OrthancPluginAnswerBuffer(OrthancPlugins::GetGlobalContext(), output, s.c_str(), s.size(), Orthanc::MIME_JSON); |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1215 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1216 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1217 |
4
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1218 static void AddDefaultTagValue(Json::Value& target, |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1219 const Orthanc::DicomTag& tag, |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1220 const std::string& value) |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1221 { |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1222 if (!target.isMember(tag.Format())) |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1223 { |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1224 target[tag.Format()] = value; |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1225 } |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1226 } |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1227 |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1228 |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1229 static void AddDefaultTagValue(Json::Value& target, |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1230 const DcmTagKey& tag, |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1231 const std::string& value) |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1232 { |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1233 AddDefaultTagValue(target, Orthanc::DicomTag(tag.getGroup(), tag.getElement()), value); |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1234 } |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1235 |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1236 |
7 | 1237 static void CallCreateDicom(Json::Value& answer, |
1238 const std::string& stl, | |
1239 const Json::Value& body, | |
1240 const std::string& parentStudy, | |
1241 const std::string& defaultSeriesDescription, | |
1242 const std::string& defaultFrameOfReferenceUid, | |
1243 const std::string& defaultTitle) | |
1244 { | |
1245 static const char* const KEY_TAGS = "Tags"; | |
1246 | |
1247 Json::Value normalized = Json::objectValue; | |
1248 | |
1249 if (body.isMember(KEY_TAGS)) | |
1250 { | |
1251 const Json::Value& tags = body[KEY_TAGS]; | |
1252 | |
1253 if (tags.type() != Json::objectValue) | |
1254 { | |
1255 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadRequest, "Tags must be provided as a JSON object"); | |
1256 } | |
1257 | |
1258 std::vector<std::string> keys = tags.getMemberNames(); | |
1259 for (size_t i = 0; i < keys.size(); i++) | |
1260 { | |
1261 const Orthanc::DicomTag tag = Orthanc::FromDcmtkBridge::ParseTag(keys[i]); | |
1262 normalized[tag.Format()] = tags[keys[i]]; | |
1263 } | |
1264 } | |
1265 | |
1266 if (!normalized.isMember(Orthanc::DICOM_TAG_SERIES_DESCRIPTION.Format())) | |
1267 { | |
1268 normalized[Orthanc::DICOM_TAG_SERIES_DESCRIPTION.Format()] = defaultSeriesDescription; | |
1269 } | |
1270 | |
1271 AddDefaultTagValue(normalized, Orthanc::DICOM_TAG_SERIES_NUMBER, "1"); | |
1272 AddDefaultTagValue(normalized, Orthanc::DICOM_TAG_FRAME_OF_REFERENCE_UID, defaultFrameOfReferenceUid); | |
1273 AddDefaultTagValue(normalized, Orthanc::DICOM_TAG_INSTANCE_NUMBER, "1"); | |
1274 | |
1275 AddDefaultTagValue(normalized, DCM_BurnedInAnnotation, "NO"); | |
1276 AddDefaultTagValue(normalized, DCM_DeviceSerialNumber, ORTHANC_STL_VERSION); | |
1277 AddDefaultTagValue(normalized, DCM_DocumentTitle, defaultTitle); | |
1278 AddDefaultTagValue(normalized, DCM_Manufacturer, "Orthanc STL plugin"); | |
1279 AddDefaultTagValue(normalized, DCM_ManufacturerModelName, "Orthanc STL plugin"); | |
1280 AddDefaultTagValue(normalized, DCM_PositionReferenceIndicator, ""); | |
1281 AddDefaultTagValue(normalized, DCM_SoftwareVersions, ORTHANC_STL_VERSION); | |
1282 AddDefaultTagValue(normalized, DCM_ConceptNameCodeSequence, ""); | |
1283 | |
1284 std::string date, time; | |
1285 Orthanc::SystemToolbox::GetNowDicom(date, time, true /* use UTC time (not local time) */); | |
1286 AddDefaultTagValue(normalized, DCM_AcquisitionDateTime, date + time); | |
1287 | |
1288 const Orthanc::DicomTag MEASUREMENT_UNITS_CODE_SEQUENCE(DCM_MeasurementUnitsCodeSequence.getGroup(), | |
1289 DCM_MeasurementUnitsCodeSequence.getElement()); | |
1290 | |
1291 if (!normalized.isMember(MEASUREMENT_UNITS_CODE_SEQUENCE.Format())) | |
1292 { | |
1293 Json::Value item; | |
1294 item["CodeValue"] = "mm"; | |
1295 item["CodingSchemeDesignator"] = "UCUM"; | |
1296 item["CodeMeaning"] = defaultTitle; | |
1297 | |
1298 normalized[MEASUREMENT_UNITS_CODE_SEQUENCE.Format()].append(item); | |
1299 } | |
1300 | |
1301 std::string content; | |
1302 Orthanc::Toolbox::EncodeDataUriScheme(content, Orthanc::MIME_STL, stl); | |
1303 | |
1304 Json::Value create; | |
1305 create["Content"] = content; | |
1306 create["Parent"] = parentStudy; | |
1307 create["Tags"] = normalized; | |
1308 | |
1309 if (!OrthancPlugins::RestApiPost(answer, "/tools/create-dicom", create.toStyledString(), false)) | |
1310 { | |
1311 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadRequest, "Cannot create DICOM from STL"); | |
1312 } | |
1313 } | |
1314 | |
1315 | |
1316 void EncodeStructureSet(OrthancPluginRestOutput* output, | |
1317 const char* url, | |
1318 const OrthancPluginHttpRequest* request) | |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1319 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1320 static const char* const KEY_INSTANCE = "Instance"; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1321 static const char* const KEY_RESOLUTION = "Resolution"; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1322 static const char* const KEY_ROI_NAMES = "RoiNames"; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1323 static const char* const KEY_SMOOTH = "Smooth"; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1324 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1325 if (request->method != OrthancPluginHttpMethod_Post) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1326 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1327 OrthancPluginSendMethodNotAllowed(OrthancPlugins::GetGlobalContext(), output, "POST"); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1328 return; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1329 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1330 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1331 Json::Value body; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1332 if (!Orthanc::Toolbox::ReadJson(body, request->body, request->bodySize)) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1333 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1334 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadRequest); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1335 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1336 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1337 const std::string instanceId = Orthanc::SerializationToolbox::ReadString(body, KEY_INSTANCE); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1338 const bool smooth = (body.isMember(KEY_SMOOTH) ? |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1339 Orthanc::SerializationToolbox::ReadBoolean(body, KEY_SMOOTH) : |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1340 true /* smooth by default */); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1341 const unsigned int resolution = (body.isMember(KEY_RESOLUTION) ? |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1342 Orthanc::SerializationToolbox::ReadUnsignedInteger(body, KEY_RESOLUTION) : |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1343 256 /* default value */); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1344 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1345 std::set<std::string> roiNames; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1346 Orthanc::SerializationToolbox::ReadSetOfStrings(roiNames, body, KEY_ROI_NAMES); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1347 |
3
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1348 std::unique_ptr<Orthanc::ParsedDicomFile> dicom(LoadInstance(instanceId)); |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1349 |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1350 StructureSet structureSet(*dicom); |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1351 |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1352 std::string stl; |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1353 if (!EncodeStructureSetMesh(stl, structureSet, roiNames, resolution, smooth)) |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1354 { |
8
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1355 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat, "Cannot encode STL from RT-STRUCT"); |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1356 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1357 else |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1358 { |
7 | 1359 std::string seriesDescription; |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1360 |
7 | 1361 if (dicom->GetTagValue(seriesDescription, Orthanc::DICOM_TAG_SERIES_DESCRIPTION)) |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1362 { |
7 | 1363 seriesDescription += ": "; |
1364 } | |
1365 else | |
1366 { | |
1367 seriesDescription.clear(); | |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1368 } |
4
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1369 |
7 | 1370 bool first = true; |
1371 for (std::set<std::string>::const_iterator it = roiNames.begin(); it != roiNames.end(); ++it) | |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1372 { |
7 | 1373 if (first) |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1374 { |
7 | 1375 first = false; |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1376 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1377 else |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1378 { |
7 | 1379 seriesDescription += ", "; |
3
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1380 } |
2
2bdb9acb7dcf
added STL viewer button at instance level too
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1
diff
changeset
|
1381 |
7 | 1382 seriesDescription += *it; |
4
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1383 } |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1384 |
7 | 1385 std::string frameOfReferenceUid; |
4
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1386 if (structureSet.HasFrameOfReferenceUid()) |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1387 { |
7 | 1388 frameOfReferenceUid = structureSet.GetFrameOfReferenceUid(); |
4
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1389 } |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1390 else |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1391 { |
7 | 1392 frameOfReferenceUid = Orthanc::FromDcmtkBridge::GenerateUniqueIdentifier(Orthanc::ResourceType_Instance); |
3
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1393 } |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1394 |
7 | 1395 Json::Value answer; |
1396 CallCreateDicom(answer, stl, body, structureSet.HashStudy(), seriesDescription, | |
1397 frameOfReferenceUid, "STL model generated from DICOM RT-STRUCT"); | |
4
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1398 |
7 | 1399 std::string s = answer.toStyledString(); |
1400 OrthancPluginAnswerBuffer(OrthancPlugins::GetGlobalContext(), output, s.c_str(), s.size(), Orthanc::MIME_JSON); | |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1401 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1402 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1403 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1404 |
3
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1405 void ExtractStl(OrthancPluginRestOutput* output, |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1406 const char* url, |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1407 const OrthancPluginHttpRequest* request) |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1408 { |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1409 if (request->method != OrthancPluginHttpMethod_Get) |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1410 { |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1411 OrthancPluginSendMethodNotAllowed(OrthancPlugins::GetGlobalContext(), output, "GET"); |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1412 return; |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1413 } |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1414 |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1415 const std::string instanceId(request->groups[0]); |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1416 |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1417 std::unique_ptr<Orthanc::ParsedDicomFile> dicom(LoadInstance(instanceId)); |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1418 DcmDataset& dataset = *dicom->GetDcmtkObject().getDataset(); |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1419 |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1420 std::string stl; |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1421 if (GetStringValue(dataset, DCM_MIMETypeOfEncapsulatedDocument) != Orthanc::MIME_STL || |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1422 GetStringValue(dataset, DCM_SOPClassUID) != "1.2.840.10008.5.1.4.1.1.104.3" || |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1423 !dicom->GetTagValue(stl, Orthanc::DICOM_TAG_ENCAPSULATED_DOCUMENT)) |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1424 { |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1425 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadRequest, "DICOM instance not encapsulating a STL model: " + instanceId); |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1426 } |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1427 else |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1428 { |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1429 OrthancPluginAnswerBuffer(OrthancPlugins::GetGlobalContext(), output, |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1430 stl.empty() ? NULL : stl.c_str(), stl.size(), Orthanc::MIME_STL); |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1431 } |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1432 } |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1433 |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1434 |
6 | 1435 |
1436 class NiftiHeader : public boost::noncopyable | |
1437 { | |
1438 private: | |
1439 nifti_image* image_; | |
1440 | |
1441 public: | |
1442 NiftiHeader(const std::string& nifti) | |
1443 { | |
1444 nifti_1_header header; | |
1445 if (nifti.size() < sizeof(header)) | |
1446 { | |
1447 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); | |
1448 } | |
1449 | |
1450 memcpy(&header, nifti.c_str(), sizeof(header)); | |
1451 if (!nifti_hdr_looks_good(&header)) | |
1452 { | |
1453 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); | |
1454 } | |
1455 | |
1456 image_ = nifti_convert_nhdr2nim(header, "dummy_filename"); | |
1457 if (image_ == NULL) | |
1458 { | |
1459 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); | |
1460 } | |
1461 } | |
1462 | |
1463 ~NiftiHeader() | |
1464 { | |
1465 nifti_image_free(image_); | |
1466 } | |
1467 | |
1468 const nifti_image& GetInfo() const | |
1469 { | |
1470 assert(image_ != NULL); | |
1471 return *image_; | |
1472 } | |
1473 }; | |
1474 | |
1475 | |
1476 static void LoadNifti(vtkImageData* volume, | |
1477 std::string& nifti) | |
1478 { | |
1479 if (volume == NULL) | |
1480 { | |
1481 throw Orthanc::OrthancException(Orthanc::ErrorCode_NullPointer); | |
1482 } | |
1483 | |
1484 const uint8_t* p = reinterpret_cast<const uint8_t*>(nifti.c_str()); | |
1485 | |
1486 if (nifti.size() >= 2 && | |
1487 p[0] == 0x1f && | |
1488 p[1] == 0x8b) | |
1489 { | |
1490 Orthanc::GzipCompressor compressor; | |
1491 std::string uncompressed; | |
1492 Orthanc::IBufferCompressor::Uncompress(uncompressed, compressor, nifti); | |
1493 nifti.swap(uncompressed); | |
1494 } | |
1495 | |
1496 NiftiHeader header(nifti); | |
1497 | |
1498 if (header.GetInfo().ndim != 3) | |
1499 { | |
1500 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat, | |
1501 "Only 3D NIfTI volumes are allowed"); | |
1502 } | |
1503 | |
8
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1504 size_t itemSize; |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1505 int vtkType; |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1506 |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1507 switch (header.GetInfo().datatype) |
6 | 1508 { |
8
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1509 case DT_UNSIGNED_CHAR: |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1510 itemSize = 1; |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1511 vtkType = VTK_UNSIGNED_CHAR; |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1512 break; |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1513 |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1514 case DT_FLOAT: |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1515 itemSize = sizeof(float); |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1516 vtkType = VTK_FLOAT; |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1517 break; |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1518 |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1519 case DT_DOUBLE: |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1520 itemSize = sizeof(double); |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1521 vtkType = VTK_DOUBLE; |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1522 break; |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1523 |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1524 default: |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1525 throw Orthanc::OrthancException(Orthanc::ErrorCode_NotImplemented); |
6 | 1526 } |
1527 | |
1528 assert(static_cast<int>(header.GetInfo().nvox) == header.GetInfo().nx * header.GetInfo().ny * header.GetInfo().nz); | |
1529 | |
1530 const size_t pixelDataOffset = sizeof(nifti_1_header) + 4 /* extension */; | |
1531 | |
8
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1532 if (nifti.size() != pixelDataOffset + header.GetInfo().nvox * itemSize) |
6 | 1533 { |
1534 throw Orthanc::OrthancException(Orthanc::ErrorCode_CorruptedFile); | |
1535 } | |
1536 | |
1537 volume->SetDimensions(header.GetInfo().nx, header.GetInfo().ny, header.GetInfo().nz); | |
8
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1538 volume->AllocateScalars(vtkType, 1); |
6 | 1539 volume->SetSpacing(header.GetInfo().dx, header.GetInfo().dy, header.GetInfo().dz); |
8
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1540 memcpy(volume->GetScalarPointer(), &nifti[pixelDataOffset], header.GetInfo().nvox * itemSize); |
6 | 1541 } |
1542 | |
8
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1543 |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1544 void EncodeNifti(OrthancPluginRestOutput* output, |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1545 const char* url, |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1546 const OrthancPluginHttpRequest* request) |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1547 { |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1548 static const char* const KEY_NIFTI = "Nifti"; |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1549 static const char* const KEY_RESOLUTION = "Resolution"; |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1550 static const char* const KEY_PARENT_STUDY = "ParentStudy"; |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1551 static const char* const KEY_SMOOTH = "Smooth"; |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1552 |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1553 if (request->method != OrthancPluginHttpMethod_Post) |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1554 { |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1555 OrthancPluginSendMethodNotAllowed(OrthancPlugins::GetGlobalContext(), output, "POST"); |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1556 return; |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1557 } |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1558 |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1559 Json::Value body; |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1560 if (!Orthanc::Toolbox::ReadJson(body, request->body, request->bodySize)) |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1561 { |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1562 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadRequest); |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1563 } |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1564 |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1565 std::string mime, nifti; |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1566 if (!Orthanc::Toolbox::DecodeDataUriScheme(mime, nifti, Orthanc::SerializationToolbox::ReadString(body, KEY_NIFTI))) |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1567 { |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1568 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadRequest, "Missing the \"Nifti\" argument containing the NIfTI file"); |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1569 } |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1570 |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1571 const std::string parentStudy = Orthanc::SerializationToolbox::ReadString(body, KEY_PARENT_STUDY); |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1572 const bool smooth = (body.isMember(KEY_SMOOTH) ? |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1573 Orthanc::SerializationToolbox::ReadBoolean(body, KEY_SMOOTH) : |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1574 true /* smooth by default */); |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1575 const unsigned int resolution = (body.isMember(KEY_RESOLUTION) ? |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1576 Orthanc::SerializationToolbox::ReadUnsignedInteger(body, KEY_RESOLUTION) : |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1577 256 /* default value */); |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1578 |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1579 vtkNew<vtkImageData> volume; |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1580 LoadNifti(volume.Get(), nifti); |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1581 |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1582 std::string stl; |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1583 if (!EncodeVolume(stl, volume.Get(), resolution, smooth)) |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1584 { |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1585 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat, "Cannot encode STL from NIfTI"); |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1586 } |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1587 else |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1588 { |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1589 const std::string title = "STL model generated from NIfTI"; |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1590 |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1591 const std::string frameOfReferenceUid = Orthanc::FromDcmtkBridge::GenerateUniqueIdentifier(Orthanc::ResourceType_Instance); |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1592 |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1593 Json::Value answer; |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1594 CallCreateDicom(answer, stl, body, parentStudy, title, frameOfReferenceUid, title); |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1595 |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1596 std::string s = answer.toStyledString(); |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1597 OrthancPluginAnswerBuffer(OrthancPlugins::GetGlobalContext(), output, s.c_str(), s.size(), Orthanc::MIME_JSON); |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1598 } |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1599 } |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1600 |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1601 |
0
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1602 extern "C" |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1603 { |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1604 ORTHANC_PLUGINS_API int32_t OrthancPluginInitialize(OrthancPluginContext* context) |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1605 { |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1606 OrthancPlugins::SetGlobalContext(context); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1607 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1608 /* Check the version of the Orthanc core */ |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1609 if (OrthancPluginCheckVersion(OrthancPlugins::GetGlobalContext()) == 0) |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1610 { |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1611 char info[1024]; |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1612 sprintf(info, "Your version of Orthanc (%s) must be above %d.%d.%d to run this plugin", |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1613 OrthancPlugins::GetGlobalContext()->orthancVersion, |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1614 ORTHANC_PLUGINS_MINIMAL_MAJOR_NUMBER, |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1615 ORTHANC_PLUGINS_MINIMAL_MINOR_NUMBER, |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1616 ORTHANC_PLUGINS_MINIMAL_REVISION_NUMBER); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1617 OrthancPluginLogError(OrthancPlugins::GetGlobalContext(), info); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1618 return -1; |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1619 } |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1620 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1621 #if ORTHANC_FRAMEWORK_VERSION_IS_ABOVE(1, 7, 2) |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1622 Orthanc::Logging::InitializePluginContext(context); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1623 #else |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1624 Orthanc::Logging::Initialize(context); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1625 #endif |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1626 |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1627 Orthanc::InitializeFramework("", false); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1628 |
0
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1629 hasCreateDicomStl_ = OrthancPlugins::CheckMinimalOrthancVersion(1, 12, 1); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1630 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1631 if (!hasCreateDicomStl_) |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1632 { |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1633 LOG(WARNING) << "Your version of Orthanc (" << std::string(context->orthancVersion) |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1634 << ") is insufficient to create DICOM STL, it should be above 1.12.1"; |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1635 } |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1636 |
25
5fc3ea6acbd5
fix deprecated calls
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
24
diff
changeset
|
1637 OrthancPlugins::SetDescription(ORTHANC_PLUGIN_NAME, "STL plugin for Orthanc."); |
0
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1638 |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1639 OrthancPlugins::RegisterRestCallback<ServeFile>("/stl/app/(.*)", true); |
3
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1640 OrthancPlugins::RegisterRestCallback<ExtractStl>("/instances/([0-9a-f-]+)/stl", true); |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1641 OrthancPlugins::RegisterRestCallback<ListStructures>("/stl/rt-struct/([0-9a-f-]+)", true); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1642 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1643 if (hasCreateDicomStl_) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1644 { |
8
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1645 OrthancPlugins::RegisterRestCallback<EncodeStructureSet>("/stl/encode-rtstruct", true); |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1646 OrthancPlugins::RegisterRestCallback<EncodeNifti>("/stl/encode-nifti", true); |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1647 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1648 |
2
2bdb9acb7dcf
added STL viewer button at instance level too
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1
diff
changeset
|
1649 // Extend the default Orthanc Explorer with custom JavaScript for STL |
2bdb9acb7dcf
added STL viewer button at instance level too
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1
diff
changeset
|
1650 std::string explorer; |
2bdb9acb7dcf
added STL viewer button at instance level too
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1
diff
changeset
|
1651 |
2bdb9acb7dcf
added STL viewer button at instance level too
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1
diff
changeset
|
1652 { |
2bdb9acb7dcf
added STL viewer button at instance level too
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1
diff
changeset
|
1653 Orthanc::EmbeddedResources::GetFileResource(explorer, Orthanc::EmbeddedResources::ORTHANC_EXPLORER); |
2bdb9acb7dcf
added STL viewer button at instance level too
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1
diff
changeset
|
1654 |
2bdb9acb7dcf
added STL viewer button at instance level too
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1
diff
changeset
|
1655 std::map<std::string, std::string> dictionary; |
2bdb9acb7dcf
added STL viewer button at instance level too
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1
diff
changeset
|
1656 dictionary["HAS_CREATE_DICOM_STL"] = (hasCreateDicomStl_ ? "true" : "false"); |
2bdb9acb7dcf
added STL viewer button at instance level too
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1
diff
changeset
|
1657 explorer = Orthanc::Toolbox::SubstituteVariables(explorer, dictionary); |
2bdb9acb7dcf
added STL viewer button at instance level too
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1
diff
changeset
|
1658 |
25
5fc3ea6acbd5
fix deprecated calls
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
24
diff
changeset
|
1659 OrthancPlugins::ExtendOrthancExplorer(ORTHANC_PLUGIN_NAME, explorer); |
2
2bdb9acb7dcf
added STL viewer button at instance level too
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1
diff
changeset
|
1660 } |
2bdb9acb7dcf
added STL viewer button at instance level too
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1
diff
changeset
|
1661 |
0
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1662 return 0; |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1663 } |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1664 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1665 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1666 ORTHANC_PLUGINS_API void OrthancPluginFinalize() |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1667 { |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1668 Orthanc::FinalizeFramework(); |
0
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1669 } |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1670 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1671 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1672 ORTHANC_PLUGINS_API const char* OrthancPluginGetName() |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1673 { |
25
5fc3ea6acbd5
fix deprecated calls
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
24
diff
changeset
|
1674 return ORTHANC_PLUGIN_NAME; |
0
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1675 } |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1676 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1677 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1678 ORTHANC_PLUGINS_API const char* OrthancPluginGetVersion() |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1679 { |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1680 return ORTHANC_STL_VERSION; |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1681 } |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1682 } |