Mercurial > hg > orthanc-stl
annotate Sources/Plugin.cpp @ 34:bee2017f3088
fix
author | Sebastien Jodogne <s.jodogne@gmail.com> |
---|---|
date | Thu, 04 Apr 2024 20:25:23 +0200 |
parents | 2460b376d3f7 |
children | ee3bc8f7df5b |
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 |
34 | 25 #include "STLToolbox.h" |
33 | 26 #include "StructurePolygon.h" |
32 | 27 #include "VTKToolbox.h" |
28 #include "Vector3D.h" | |
29 #include "Toolbox.h" | |
30 #include "Extent2D.h" | |
31 | |
0
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
32 #include "../Resources/Orthanc/Plugins/OrthancPluginCppWrapper.h" |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
33 |
4
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
34 #include <EmbeddedResources.h> |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
35 |
28
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
36 #include <DicomFormat/DicomInstanceHasher.h> |
4
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
37 #include <DicomParsing/FromDcmtkBridge.h> |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
38 #include <DicomParsing/ParsedDicomFile.h> |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
39 #include <Images/ImageProcessing.h> |
0
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
40 #include <Logging.h> |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
41 #include <OrthancFramework.h> |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
42 #include <SerializationToolbox.h> |
0
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
43 #include <SystemToolbox.h> |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
44 |
33 | 45 #include <vtkNew.h> |
46 | |
0
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
47 #include <boost/thread/shared_mutex.hpp> |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
48 |
25
5fc3ea6acbd5
fix deprecated calls
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
24
diff
changeset
|
49 #define ORTHANC_PLUGIN_NAME "stl" |
5fc3ea6acbd5
fix deprecated calls
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
24
diff
changeset
|
50 |
6 | 51 |
0
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
52 // Forward declaration |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
53 void ReadStaticAsset(std::string& target, |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
54 const std::string& path); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
55 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
56 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
57 /** |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
58 * 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
|
59 * "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
|
60 * 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
|
61 **/ |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
62 class ResourcesCache : public boost::noncopyable |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
63 { |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
64 private: |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
65 typedef std::map<std::string, std::string*> Content; |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
66 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
67 boost::shared_mutex mutex_; |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
68 Content content_; |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
69 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
70 public: |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
71 ~ResourcesCache() |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
72 { |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
73 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
|
74 { |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
75 assert(it->second != NULL); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
76 delete it->second; |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
77 } |
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 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
80 void Answer(OrthancPluginRestOutput* output, |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
81 const std::string& path) |
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 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
|
84 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
85 { |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
86 // Check whether the cache already contains the resource |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
87 boost::shared_lock<boost::shared_mutex> lock(mutex_); |
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 Content::const_iterator found = content_.find(path); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
90 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
91 if (found != content_.end()) |
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 assert(found->second != NULL); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
94 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
|
95 return; |
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 } |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
98 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
99 // This resource has not been cached yet |
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 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
|
102 ReadStaticAsset(*item, path); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
103 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
|
104 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
105 { |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
106 // Store the resource into the cache |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
107 boost::unique_lock<boost::shared_mutex> lock(mutex_); |
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 if (content_.find(path) == content_.end()) |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
110 { |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
111 content_[path] = item.release(); |
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 } |
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 }; |
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 static ResourcesCache cache_; |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
119 static bool hasCreateDicomStl_; |
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 void ServeFile(OrthancPluginRestOutput* output, |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
122 const char* url, |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
123 const OrthancPluginHttpRequest* request) |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
124 { |
3
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
125 if (request->method != OrthancPluginHttpMethod_Get) |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
126 { |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
127 OrthancPluginSendMethodNotAllowed(OrthancPlugins::GetGlobalContext(), output, "GET"); |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
128 return; |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
129 } |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
130 |
0
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
131 std::string file = request->groups[0]; |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
132 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
133 if (file == "viewer.html") |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
134 { |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
135 std::string s; |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
136 Orthanc::EmbeddedResources::GetFileResource(s, Orthanc::EmbeddedResources::VIEWER_HTML); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
137 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
|
138 } |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
139 else if (file == "viewer.js") |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
140 { |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
141 std::string s; |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
142 Orthanc::EmbeddedResources::GetFileResource(s, Orthanc::EmbeddedResources::VIEWER_JS); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
143 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
|
144 } |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
145 else |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
146 { |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
147 cache_.Answer(output, file); |
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 } |
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 |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
152 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
153 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
154 #include <dcmtk/dcmdata/dcfilefo.h> |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
155 #include <dcmtk/dcmdata/dcsequen.h> |
29
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
156 #include <dcmtk/dcmdata/dcuid.h> |
1
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 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
159 class StructureSet : public boost::noncopyable |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
160 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
161 private: |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
162 std::vector<StructurePolygon*> polygons_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
163 std::string patientId_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
164 std::string studyInstanceUid_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
165 std::string seriesInstanceUid_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
166 std::string sopInstanceUid_; |
4
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
167 bool hasFrameOfReferenceUid_; |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
168 std::string frameOfReferenceUid_; |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
169 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
170 public: |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
171 explicit StructureSet(Orthanc::ParsedDicomFile& dicom) : |
4
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
172 hasFrameOfReferenceUid_(false) |
1
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 DcmDataset& dataset = *dicom.GetDcmtkObject().getDataset(); |
33 | 175 patientId_ = STLToolbox::GetStringValue(dataset, DCM_PatientID); |
176 studyInstanceUid_ = STLToolbox::GetStringValue(dataset, DCM_StudyInstanceUID); | |
177 seriesInstanceUid_ = STLToolbox::GetStringValue(dataset, DCM_SeriesInstanceUID); | |
178 sopInstanceUid_ = STLToolbox::GetStringValue(dataset, DCM_SOPInstanceUID); | |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
179 |
4
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
180 DcmSequenceOfItems* frame = NULL; |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
181 if (!dataset.findAndGetSequence(DCM_ReferencedFrameOfReferenceSequence, frame).good() || |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
182 frame == NULL) |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
183 { |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
184 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
185 } |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
186 |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
187 if (frame->card() == 1) |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
188 { |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
189 const char* v = NULL; |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
190 if (frame->getItem(0)->findAndGetString(DCM_FrameOfReferenceUID, v).good() && |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
191 v != NULL) |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
192 { |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
193 hasFrameOfReferenceUid_ = true; |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
194 frameOfReferenceUid_.assign(v); |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
195 } |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
196 } |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
197 |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
198 DcmSequenceOfItems* rois = NULL; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
199 if (!dataset.findAndGetSequence(DCM_ROIContourSequence, rois).good() || |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
200 rois == NULL) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
201 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
202 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
203 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
204 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
205 std::vector<DcmSequenceOfItems*> contours(rois->card()); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
206 size_t countPolygons = 0; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
207 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
208 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
|
209 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
210 DcmSequenceOfItems* contour = NULL; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
211 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
|
212 contour == NULL) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
213 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
214 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
215 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
216 else |
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 contours[i] = contour; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
219 countPolygons += contour->card(); |
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 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
222 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
223 polygons_.resize(countPolygons); |
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 size_t pos = 0; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
226 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
|
227 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
228 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
|
229 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
230 polygons_[pos] = new StructurePolygon(dicom, i, j); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
231 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
232 } |
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 assert(pos == countPolygons); |
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 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
237 ~StructureSet() |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
238 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
239 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
|
240 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
241 assert(polygons_[i] != NULL); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
242 delete polygons_[i]; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
243 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
244 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
245 |
29
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
246 const std::string& GetPatientId() const |
1
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 return patientId_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
249 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
250 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
251 const std::string& GetStudyInstanceUid() const |
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 return studyInstanceUid_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
254 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
255 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
256 const std::string& GetSeriesInstanceUid() const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
257 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
258 return seriesInstanceUid_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
259 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
260 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
261 const std::string& GetSopInstanceUid() const |
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 return sopInstanceUid_; |
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 std::string HashStudy() const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
267 { |
28
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
268 Orthanc::DicomInstanceHasher hasher(patientId_, studyInstanceUid_, seriesInstanceUid_, sopInstanceUid_); |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
269 return hasher.HashStudy(); |
1
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 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
272 size_t GetPolygonsCount() const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
273 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
274 return polygons_.size(); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
275 } |
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 const StructurePolygon& GetPolygon(size_t i) const |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
278 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
279 if (i >= polygons_.size()) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
280 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
281 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange); |
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 else |
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 assert(polygons_[i] != NULL); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
286 return *polygons_[i]; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
287 } |
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 |
4
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
290 bool HasFrameOfReferenceUid() const |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
291 { |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
292 return hasFrameOfReferenceUid_; |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
293 } |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
294 |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
295 const std::string& GetFrameOfReferenceUid() const |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
296 { |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
297 if (hasFrameOfReferenceUid_) |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
298 { |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
299 return frameOfReferenceUid_; |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
300 } |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
301 else |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
302 { |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
303 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls); |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
304 } |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
305 } |
33 | 306 |
307 // This static method is faster than constructing the full "StructureSet" object | |
308 static void ListStructuresNames(std::set<std::string>& target, | |
309 Orthanc::ParsedDicomFile& source) | |
310 { | |
311 target.clear(); | |
312 | |
313 DcmSequenceOfItems* sequence = NULL; | |
314 if (!source.GetDcmtkObject().getDataset()->findAndGetSequence(DCM_StructureSetROISequence, sequence).good() || | |
315 sequence == NULL) | |
316 { | |
317 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); | |
318 } | |
319 | |
320 for (unsigned long i = 0; i < sequence->card(); i++) | |
321 { | |
322 DcmItem* item = sequence->getItem(i); | |
323 if (item == NULL) | |
324 { | |
325 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); | |
326 } | |
327 else | |
328 { | |
329 target.insert(STLToolbox::GetStringValue(*item, DCM_ROIName)); | |
330 } | |
331 } | |
332 } | |
1
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 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
335 |
28
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
336 class StructureSetGeometry : public boost::noncopyable |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
337 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
338 private: |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
339 bool strict_; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
340 Vector3D slicesNormal_; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
341 double slicesSpacing_; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
342 double minProjectionAlongNormal_; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
343 double maxProjectionAlongNormal_; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
344 size_t slicesCount_; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
345 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
346 bool LookupProjectionIndex(size_t& index, |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
347 double z) const |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
348 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
349 if (slicesCount_ == 0) |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
350 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
351 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError); |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
352 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
353 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
354 if (z < minProjectionAlongNormal_ || |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
355 z > maxProjectionAlongNormal_) |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
356 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
357 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange); |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
358 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
359 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
360 assert(slicesSpacing_ > 0 && |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
361 minProjectionAlongNormal_ < maxProjectionAlongNormal_); |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
362 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
363 double d = (z - minProjectionAlongNormal_) / slicesSpacing_; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
364 |
33 | 365 if (STLToolbox::IsNear(d, round(d))) |
28
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
366 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
367 if (d < 0.0 || |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
368 d > static_cast<double>(slicesCount_) - 1.0) |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
369 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
370 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError); |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
371 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
372 else |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
373 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
374 index = static_cast<size_t>(round(d)); |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
375 return true; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
376 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
377 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
378 else |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
379 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
380 return false; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
381 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
382 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
383 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
384 public: |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
385 StructureSetGeometry(const StructureSet& structures, |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
386 bool strict) : |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
387 strict_(strict) |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
388 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
389 bool isValid = false; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
390 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
391 std::vector<double> projections; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
392 projections.reserve(structures.GetPolygonsCount()); |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
393 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
394 for (size_t i = 0; i < structures.GetPolygonsCount(); i++) |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
395 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
396 Vector3D normal; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
397 if (structures.GetPolygon(i).IsCoplanar(normal)) |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
398 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
399 // Initialize the normal of the whole volume, if need be |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
400 if (!isValid) |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
401 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
402 isValid = true; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
403 slicesNormal_ = normal; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
404 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
405 |
29
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
406 if (Vector3D::AreParallel(normal, slicesNormal_)) |
28
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
407 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
408 // This is a valid slice (it is parallel to the normal) |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
409 const Vector3D& point = structures.GetPolygon(i).GetPoint(0); |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
410 projections.push_back(Vector3D::DotProduct(point, slicesNormal_)); |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
411 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
412 else |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
413 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
414 // RT-STRUCT with non-parallel slices |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
415 throw Orthanc::OrthancException(Orthanc::ErrorCode_NotImplemented); |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
416 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
417 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
418 else |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
419 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
420 // Ignore slices that are not coplanar |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
421 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
422 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
423 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
424 if (projections.empty()) |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
425 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
426 throw Orthanc::OrthancException(Orthanc::ErrorCode_NotImplemented, |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
427 "Structure set without a valid geometry"); |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
428 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
429 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
430 // Only keep unique projections |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
431 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
432 std::sort(projections.begin(), projections.end()); |
33 | 433 STLToolbox::RemoveDuplicateValues(projections); |
28
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
434 assert(!projections.empty()); |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
435 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
436 if (projections.size() == 1) |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
437 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
438 // Volume with one single slice |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
439 minProjectionAlongNormal_ = projections[0]; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
440 maxProjectionAlongNormal_ = projections[0]; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
441 slicesSpacing_ = 1; // Arbitrary value |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
442 slicesCount_ = 1; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
443 return; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
444 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
445 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
446 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
447 // Compute the most probable spacing between the slices |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
448 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
449 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
450 std::vector<double> spacings; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
451 spacings.resize(projections.size() - 1); |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
452 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
453 for (size_t i = 0; i < spacings.size(); i++) |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
454 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
455 spacings[i] = projections[i + 1] - projections[i]; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
456 assert(spacings[i] > 0); |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
457 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
458 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
459 std::sort(spacings.begin(), spacings.end()); |
33 | 460 STLToolbox::RemoveDuplicateValues(spacings); |
28
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
461 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
462 if (spacings.empty()) |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
463 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
464 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError); |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
465 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
466 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
467 slicesSpacing_ = spacings[spacings.size() / 10]; // Take the 90% percentile of smallest spacings |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
468 assert(slicesSpacing_ > 0); |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
469 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
470 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
471 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
472 // Find the projection along the normal with the largest support |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
473 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
474 bool first = true; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
475 size_t bestSupport; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
476 double bestProjection; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
477 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
478 std::list<size_t> candidates; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
479 for (size_t i = 0; i < projections.size(); i++) |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
480 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
481 candidates.push_back(i); |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
482 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
483 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
484 while (!candidates.empty()) |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
485 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
486 std::list<size_t> next; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
487 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
488 size_t countSupport = 0; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
489 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
490 std::list<size_t>::const_iterator it = candidates.begin(); |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
491 size_t reference = *it; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
492 it++; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
493 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
494 while (it != candidates.end()) |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
495 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
496 double d = (projections[*it] - projections[reference]) / slicesSpacing_; |
33 | 497 if (STLToolbox::IsNear(d, round(d))) |
28
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
498 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
499 countSupport ++; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
500 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
501 else |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
502 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
503 next.push_back(*it); |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
504 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
505 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
506 it++; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
507 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
508 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
509 if (first || |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
510 countSupport > bestSupport) |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
511 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
512 first = false; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
513 bestSupport = countSupport; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
514 bestProjection = projections[reference]; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
515 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
516 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
517 if (strict && |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
518 !next.empty()) |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
519 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
520 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat, |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
521 "Structure set with multiple support, which is not allowed in Strict mode"); |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
522 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
523 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
524 candidates.swap(next); |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
525 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
526 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
527 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
528 // Compute the range of the projections |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
529 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
530 minProjectionAlongNormal_ = bestProjection; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
531 maxProjectionAlongNormal_ = bestProjection; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
532 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
533 for (size_t i = 0; i < projections.size(); i++) |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
534 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
535 double d = (projections[i] - bestProjection) / slicesSpacing_; |
33 | 536 if (STLToolbox::IsNear(d, round(d))) |
28
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
537 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
538 minProjectionAlongNormal_ = std::min(minProjectionAlongNormal_, projections[i]); |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
539 maxProjectionAlongNormal_ = std::max(maxProjectionAlongNormal_, projections[i]); |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
540 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
541 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
542 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
543 double d = (maxProjectionAlongNormal_ - minProjectionAlongNormal_) / slicesSpacing_; |
33 | 544 if (STLToolbox::IsNear(d, round(d))) |
28
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
545 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
546 slicesCount_ = static_cast<size_t>(round(d)) + 1; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
547 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
548 else |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
549 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
550 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError); |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
551 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
552 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
553 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
554 // Sanity check |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
555 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
556 size_t a, b; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
557 if (!LookupProjectionIndex(a, minProjectionAlongNormal_) || |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
558 !LookupProjectionIndex(b, maxProjectionAlongNormal_) || |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
559 a != 0 || |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
560 b + 1 != slicesCount_) |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
561 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
562 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError); |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
563 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
564 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
565 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
566 const Vector3D& GetSlicesNormal() const |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
567 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
568 return slicesNormal_; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
569 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
570 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
571 double GetSlicesSpacing() const |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
572 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
573 return slicesSpacing_; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
574 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
575 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
576 double GetMinProjectionAlongNormal() const |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
577 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
578 return minProjectionAlongNormal_; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
579 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
580 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
581 double GetMaxProjectionAlongNormal() const |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
582 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
583 return maxProjectionAlongNormal_; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
584 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
585 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
586 bool ProjectAlongNormal(double& z, |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
587 const StructurePolygon& polygon) const |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
588 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
589 Vector3D normal; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
590 if (polygon.IsCoplanar(normal) && |
29
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
591 Vector3D::AreParallel(normal, slicesNormal_)) |
28
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
592 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
593 z = Vector3D::DotProduct(polygon.GetPoint(0), slicesNormal_); |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
594 return true; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
595 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
596 else |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
597 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
598 return false; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
599 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
600 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
601 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
602 size_t GetSlicesCount() const |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
603 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
604 return slicesCount_; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
605 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
606 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
607 bool LookupSliceIndex(size_t& slice, |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
608 const StructurePolygon& polygon) const |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
609 { |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
610 double z; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
611 return (ProjectAlongNormal(z, polygon) && |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
612 LookupProjectionIndex(slice, z)); |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
613 } |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
614 }; |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
615 |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
616 |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
617 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
618 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
619 class XorFiller : public Orthanc::ImageProcessing::IPolygonFiller |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
620 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
621 private: |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
622 Orthanc::ImageAccessor& target_; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
623 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
624 public: |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
625 XorFiller(Orthanc::ImageAccessor& target) : |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
626 target_(target) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
627 { |
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 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
630 virtual void Fill(int y, |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
631 int x1, |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
632 int x2) ORTHANC_OVERRIDE |
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 assert(x1 <= x2); |
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 if (y >= 0 && |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
637 y < static_cast<int>(target_.GetHeight())) |
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 x1 = std::max(x1, 0); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
640 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
|
641 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
642 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
|
643 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
644 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
|
645 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
646 *p = (*p ^ 0xff); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
647 } |
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 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
650 }; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
651 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
652 |
29
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
653 static Orthanc::ParsedDicomFile* LoadInstance(const std::string& instanceId) |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
654 { |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
655 std::string dicom; |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
656 |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
657 if (!OrthancPlugins::RestApiGetString(dicom, "/instances/" + instanceId + "/file", false)) |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
658 { |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
659 throw Orthanc::OrthancException(Orthanc::ErrorCode_UnknownResource); |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
660 } |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
661 else |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
662 { |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
663 return new Orthanc::ParsedDicomFile(dicom); |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
664 } |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
665 } |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
666 |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
667 |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
668 static void GetReferencedVolumeAxes(Vector3D& axisX, |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
669 Vector3D& axisY, |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
670 const StructureSet& structureSet) |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
671 { |
29
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
672 // This is a rough guess for the X/Y axes |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
673 axisX = Vector3D(1, 0, 0); |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
674 axisY = Vector3D(0, 1, 0); |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
675 |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
676 // Look for one instance from the referenced volume |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
677 const std::string& sopInstanceUid = structureSet.GetPolygon(0).GetReferencedSopInstanceUid(); |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
678 |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
679 Json::Value response; |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
680 if (OrthancPlugins::RestApiPost(response, "/tools/lookup", sopInstanceUid, false)) |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
681 { |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
682 if (response.type() != Json::arrayValue) |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
683 { |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
684 throw Orthanc::OrthancException(Orthanc::ErrorCode_NetworkProtocol); |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
685 } |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
686 |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
687 bool first = true; |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
688 |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
689 for (Json::Value::ArrayIndex i = 0; i < response.size(); i++) |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
690 { |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
691 if (response[i].type() != Json::objectValue) |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
692 { |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
693 throw Orthanc::OrthancException(Orthanc::ErrorCode_NetworkProtocol); |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
694 } |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
695 |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
696 if (Orthanc::SerializationToolbox::ReadString(response[i], "Type") == "Instance") |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
697 { |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
698 if (first) |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
699 { |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
700 const std::string& instanceId = Orthanc::SerializationToolbox::ReadString(response[i], "ID"); |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
701 std::unique_ptr<Orthanc::ParsedDicomFile> reference(LoadInstance(instanceId)); |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
702 |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
703 std::string imageOrientation; |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
704 if (reference->GetTagValue(imageOrientation, Orthanc::DICOM_TAG_IMAGE_ORIENTATION_PATIENT)) |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
705 { |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
706 std::vector<std::string> items; |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
707 Orthanc::Toolbox::TokenizeString(items, imageOrientation, '\\'); |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
708 |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
709 double x1, x2, x3, y1, y2, y3; |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
710 |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
711 if (items.size() == 6 && |
33 | 712 STLToolbox::MyParseDouble(x1, items[0]) && |
713 STLToolbox::MyParseDouble(x2, items[1]) && | |
714 STLToolbox::MyParseDouble(x3, items[2]) && | |
715 STLToolbox::MyParseDouble(y1, items[3]) && | |
716 STLToolbox::MyParseDouble(y2, items[4]) && | |
717 STLToolbox::MyParseDouble(y3, items[5])) | |
29
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
718 { |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
719 axisX = Vector3D(x1, x2, x3); |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
720 axisY = Vector3D(y1, y2, y3); |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
721 } |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
722 } |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
723 } |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
724 else |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
725 { |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
726 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError, |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
727 "Multiple instances with the same SOP Instance UID"); |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
728 } |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
729 } |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
730 } |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
731 } |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
732 } |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
733 |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
734 |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
735 static bool EncodeStructureSetMesh(std::string& stl, |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
736 const StructureSet& structureSet, |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
737 const StructureSetGeometry& geometry, |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
738 const std::set<std::string>& roiNames, |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
739 unsigned int resolution, |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
740 bool smooth) |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
741 { |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
742 if (resolution < 1 || |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
743 structureSet.GetPolygonsCount() < 1) |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
744 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
745 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
746 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
747 |
33 | 748 if (!STLToolbox::IsNear(1, geometry.GetSlicesNormal().ComputeNorm())) |
1
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 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
751 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
752 |
29
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
753 Vector3D axisX, axisY; |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
754 GetReferencedVolumeAxes(axisX, axisY, structureSet); |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
755 |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
756 Vector3D axisZ = Vector3D::CrossProduct(axisX, axisY); |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
757 |
33 | 758 if (!STLToolbox::IsNear(1, axisX.ComputeNorm()) || |
759 !STLToolbox::IsNear(1, axisY.ComputeNorm()) || | |
29
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
760 !Vector3D::AreParallel(axisZ, geometry.GetSlicesNormal())) |
1
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 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError); |
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 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
765 Extent2D extent; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
766 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
|
767 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
768 structureSet.GetPolygon(i).Add(extent, axisX, axisY); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
769 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
770 |
28
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
771 const int depth = geometry.GetSlicesCount(); |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
772 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
773 vtkNew<vtkImageData> volume; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
774 volume->SetDimensions(resolution, resolution, depth); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
775 volume->AllocateScalars(VTK_UNSIGNED_CHAR, 1); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
776 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
777 assert(sizeof(unsigned char) == 1); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
778 memset(volume->GetScalarPointer(), 0, resolution * resolution * depth); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
779 |
30 | 780 const double pixelSpacingX = extent.GetWidth() / static_cast<double>(resolution); |
781 const double pixelSpacingY = extent.GetHeight() / static_cast<double>(resolution); | |
782 const double pixelSpacingZ = geometry.GetSlicesSpacing(); | |
783 | |
784 bool first = true; | |
785 | |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
786 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
|
787 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
788 const StructurePolygon& polygon = structureSet.GetPolygon(i); |
29
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
789 if (roiNames.find(polygon.GetRoiName()) != roiNames.end()) |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
790 { |
29
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
791 // This polygon corresponds to a ROI of interest |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
792 |
29
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
793 size_t z; |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
794 if (geometry.LookupSliceIndex(z, polygon)) |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
795 { |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
796 std::vector<Orthanc::ImageProcessing::ImagePoint> points; |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
797 points.reserve(polygon.GetPointsCount()); |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
798 |
29
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
799 for (size_t j = 0; j < polygon.GetPointsCount(); j++) |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
800 { |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
801 const Vector3D& point = polygon.GetPoint(j); |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
802 double x = (Vector3D::DotProduct(point, axisX) - extent.GetMinX()) / extent.GetWidth() * static_cast<double>(resolution); |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
803 double y = (Vector3D::DotProduct(point, axisY) - extent.GetMinY()) / extent.GetHeight() * static_cast<double>(resolution); |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
804 points.push_back(Orthanc::ImageProcessing::ImagePoint(static_cast<int32_t>(std::floor(x)), |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
805 static_cast<int32_t>(std::floor(y)))); |
30 | 806 |
807 if (first) | |
808 { | |
809 first = false; | |
810 volume->SetOrigin(point.GetX() - x * pixelSpacingX, | |
811 point.GetY() - y * pixelSpacingY, | |
812 point.GetZ() - z * pixelSpacingZ); | |
813 } | |
29
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
814 } |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
815 |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
816 Orthanc::ImageAccessor slice; |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
817 slice.AssignWritable(Orthanc::PixelFormat_Grayscale8, resolution, resolution, resolution /* pitch */, |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
818 reinterpret_cast<uint8_t*>(volume->GetScalarPointer()) + z * resolution * resolution); |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
819 |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
820 XorFiller filler(slice); |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
821 Orthanc::ImageProcessing::FillPolygon(filler, points); |
62abf3c523f9
big endian support, retrieve axes from referenced volume
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
28
diff
changeset
|
822 } |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
823 } |
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 |
30 | 826 volume->SetSpacing(pixelSpacingX, pixelSpacingY, pixelSpacingZ); |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
827 |
32 | 828 return VTKToolbox::EncodeVolume(stl, volume.Get(), resolution, smooth); |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
829 } |
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 void ListStructures(OrthancPluginRestOutput* output, |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
833 const char* url, |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
834 const OrthancPluginHttpRequest* request) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
835 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
836 if (request->method != OrthancPluginHttpMethod_Get) |
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 OrthancPluginSendMethodNotAllowed(OrthancPlugins::GetGlobalContext(), output, "GET"); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
839 return; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
840 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
841 |
3
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
842 const std::string instanceId(request->groups[0]); |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
843 |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
844 std::unique_ptr<Orthanc::ParsedDicomFile> dicom(LoadInstance(instanceId)); |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
845 |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
846 std::set<std::string> names; |
33 | 847 StructureSet::ListStructuresNames(names, *dicom); |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
848 |
3
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
849 Json::Value answer = Json::arrayValue; |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
850 |
3
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
851 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
|
852 { |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
853 answer.append(*it); |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
854 } |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
855 |
3
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
856 std::string s = answer.toStyledString(); |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
857 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
|
858 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
859 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
860 |
4
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
861 static void AddDefaultTagValue(Json::Value& target, |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
862 const Orthanc::DicomTag& tag, |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
863 const std::string& value) |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
864 { |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
865 if (!target.isMember(tag.Format())) |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
866 { |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
867 target[tag.Format()] = value; |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
868 } |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
869 } |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
870 |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
871 |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
872 static void AddDefaultTagValue(Json::Value& target, |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
873 const DcmTagKey& tag, |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
874 const std::string& value) |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
875 { |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
876 AddDefaultTagValue(target, Orthanc::DicomTag(tag.getGroup(), tag.getElement()), value); |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
877 } |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
878 |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
879 |
7 | 880 static void CallCreateDicom(Json::Value& answer, |
881 const std::string& stl, | |
882 const Json::Value& body, | |
883 const std::string& parentStudy, | |
884 const std::string& defaultSeriesDescription, | |
885 const std::string& defaultFrameOfReferenceUid, | |
886 const std::string& defaultTitle) | |
887 { | |
888 static const char* const KEY_TAGS = "Tags"; | |
889 | |
890 Json::Value normalized = Json::objectValue; | |
891 | |
892 if (body.isMember(KEY_TAGS)) | |
893 { | |
894 const Json::Value& tags = body[KEY_TAGS]; | |
895 | |
896 if (tags.type() != Json::objectValue) | |
897 { | |
898 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadRequest, "Tags must be provided as a JSON object"); | |
899 } | |
900 | |
901 std::vector<std::string> keys = tags.getMemberNames(); | |
902 for (size_t i = 0; i < keys.size(); i++) | |
903 { | |
904 const Orthanc::DicomTag tag = Orthanc::FromDcmtkBridge::ParseTag(keys[i]); | |
905 normalized[tag.Format()] = tags[keys[i]]; | |
906 } | |
907 } | |
908 | |
909 if (!normalized.isMember(Orthanc::DICOM_TAG_SERIES_DESCRIPTION.Format())) | |
910 { | |
911 normalized[Orthanc::DICOM_TAG_SERIES_DESCRIPTION.Format()] = defaultSeriesDescription; | |
912 } | |
913 | |
914 AddDefaultTagValue(normalized, Orthanc::DICOM_TAG_SERIES_NUMBER, "1"); | |
915 AddDefaultTagValue(normalized, Orthanc::DICOM_TAG_FRAME_OF_REFERENCE_UID, defaultFrameOfReferenceUid); | |
916 AddDefaultTagValue(normalized, Orthanc::DICOM_TAG_INSTANCE_NUMBER, "1"); | |
917 | |
918 AddDefaultTagValue(normalized, DCM_BurnedInAnnotation, "NO"); | |
919 AddDefaultTagValue(normalized, DCM_DeviceSerialNumber, ORTHANC_STL_VERSION); | |
920 AddDefaultTagValue(normalized, DCM_DocumentTitle, defaultTitle); | |
921 AddDefaultTagValue(normalized, DCM_Manufacturer, "Orthanc STL plugin"); | |
922 AddDefaultTagValue(normalized, DCM_ManufacturerModelName, "Orthanc STL plugin"); | |
923 AddDefaultTagValue(normalized, DCM_PositionReferenceIndicator, ""); | |
924 AddDefaultTagValue(normalized, DCM_SoftwareVersions, ORTHANC_STL_VERSION); | |
925 AddDefaultTagValue(normalized, DCM_ConceptNameCodeSequence, ""); | |
926 | |
927 std::string date, time; | |
928 Orthanc::SystemToolbox::GetNowDicom(date, time, true /* use UTC time (not local time) */); | |
929 AddDefaultTagValue(normalized, DCM_AcquisitionDateTime, date + time); | |
930 | |
931 const Orthanc::DicomTag MEASUREMENT_UNITS_CODE_SEQUENCE(DCM_MeasurementUnitsCodeSequence.getGroup(), | |
932 DCM_MeasurementUnitsCodeSequence.getElement()); | |
933 | |
934 if (!normalized.isMember(MEASUREMENT_UNITS_CODE_SEQUENCE.Format())) | |
935 { | |
936 Json::Value item; | |
937 item["CodeValue"] = "mm"; | |
938 item["CodingSchemeDesignator"] = "UCUM"; | |
939 item["CodeMeaning"] = defaultTitle; | |
940 | |
941 normalized[MEASUREMENT_UNITS_CODE_SEQUENCE.Format()].append(item); | |
942 } | |
943 | |
944 std::string content; | |
945 Orthanc::Toolbox::EncodeDataUriScheme(content, Orthanc::MIME_STL, stl); | |
946 | |
947 Json::Value create; | |
948 create["Content"] = content; | |
949 create["Parent"] = parentStudy; | |
950 create["Tags"] = normalized; | |
951 | |
952 if (!OrthancPlugins::RestApiPost(answer, "/tools/create-dicom", create.toStyledString(), false)) | |
953 { | |
954 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadRequest, "Cannot create DICOM from STL"); | |
955 } | |
956 } | |
957 | |
958 | |
959 void EncodeStructureSet(OrthancPluginRestOutput* output, | |
960 const char* url, | |
961 const OrthancPluginHttpRequest* request) | |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
962 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
963 static const char* const KEY_INSTANCE = "Instance"; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
964 static const char* const KEY_RESOLUTION = "Resolution"; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
965 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
|
966 static const char* const KEY_SMOOTH = "Smooth"; |
28
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
967 static const char* const KEY_STRICT = "Strict"; |
1
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 if (request->method != OrthancPluginHttpMethod_Post) |
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 OrthancPluginSendMethodNotAllowed(OrthancPlugins::GetGlobalContext(), output, "POST"); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
972 return; |
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 Json::Value body; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
976 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
|
977 { |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
978 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadRequest); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
979 } |
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 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
|
982 const bool smooth = (body.isMember(KEY_SMOOTH) ? |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
983 Orthanc::SerializationToolbox::ReadBoolean(body, KEY_SMOOTH) : |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
984 true /* smooth by default */); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
985 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
|
986 Orthanc::SerializationToolbox::ReadUnsignedInteger(body, KEY_RESOLUTION) : |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
987 256 /* default value */); |
28
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
988 const bool strict = (body.isMember(KEY_STRICT) ? |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
989 Orthanc::SerializationToolbox::ReadBoolean(body, KEY_STRICT) : |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
990 true /* strict by default */); |
1
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 std::set<std::string> roiNames; |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
993 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
|
994 |
3
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
995 std::unique_ptr<Orthanc::ParsedDicomFile> dicom(LoadInstance(instanceId)); |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
996 |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
997 StructureSet structureSet(*dicom); |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
998 |
28
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
999 StructureSetGeometry geometry(structureSet, strict); |
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
1000 |
3
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1001 std::string stl; |
28
410003c50b17
improved computation of RT-STRUCT geometry
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
25
diff
changeset
|
1002 if (!EncodeStructureSetMesh(stl, structureSet, geometry, roiNames, resolution, smooth)) |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1003 { |
8
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1004 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
|
1005 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1006 else |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1007 { |
7 | 1008 std::string seriesDescription; |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1009 |
7 | 1010 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
|
1011 { |
7 | 1012 seriesDescription += ": "; |
1013 } | |
1014 else | |
1015 { | |
1016 seriesDescription.clear(); | |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1017 } |
4
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1018 |
7 | 1019 bool first = true; |
1020 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
|
1021 { |
7 | 1022 if (first) |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1023 { |
7 | 1024 first = false; |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1025 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1026 else |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1027 { |
7 | 1028 seriesDescription += ", "; |
3
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1029 } |
2
2bdb9acb7dcf
added STL viewer button at instance level too
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1
diff
changeset
|
1030 |
7 | 1031 seriesDescription += *it; |
4
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1032 } |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1033 |
7 | 1034 std::string frameOfReferenceUid; |
4
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1035 if (structureSet.HasFrameOfReferenceUid()) |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1036 { |
7 | 1037 frameOfReferenceUid = structureSet.GetFrameOfReferenceUid(); |
4
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1038 } |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1039 else |
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1040 { |
7 | 1041 frameOfReferenceUid = Orthanc::FromDcmtkBridge::GenerateUniqueIdentifier(Orthanc::ResourceType_Instance); |
3
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1042 } |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1043 |
7 | 1044 Json::Value answer; |
1045 CallCreateDicom(answer, stl, body, structureSet.HashStudy(), seriesDescription, | |
1046 frameOfReferenceUid, "STL model generated from DICOM RT-STRUCT"); | |
4
5ee4448a8ff8
generated STL pass dciodvfy
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
3
diff
changeset
|
1047 |
7 | 1048 std::string s = answer.toStyledString(); |
1049 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
|
1050 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1051 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1052 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1053 |
3
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1054 void ExtractStl(OrthancPluginRestOutput* output, |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1055 const char* url, |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1056 const OrthancPluginHttpRequest* request) |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1057 { |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1058 if (request->method != OrthancPluginHttpMethod_Get) |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1059 { |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1060 OrthancPluginSendMethodNotAllowed(OrthancPlugins::GetGlobalContext(), output, "GET"); |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1061 return; |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1062 } |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1063 |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1064 const std::string instanceId(request->groups[0]); |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1065 |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1066 std::unique_ptr<Orthanc::ParsedDicomFile> dicom(LoadInstance(instanceId)); |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1067 DcmDataset& dataset = *dicom->GetDcmtkObject().getDataset(); |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1068 |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1069 std::string stl; |
33 | 1070 if (STLToolbox::GetStringValue(dataset, DCM_MIMETypeOfEncapsulatedDocument) != Orthanc::MIME_STL || |
1071 STLToolbox::GetStringValue(dataset, DCM_SOPClassUID) != UID_EncapsulatedSTLStorage || | |
3
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1072 !dicom->GetTagValue(stl, Orthanc::DICOM_TAG_ENCAPSULATED_DOCUMENT)) |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1073 { |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1074 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
|
1075 } |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1076 else |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1077 { |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1078 OrthancPluginAnswerBuffer(OrthancPlugins::GetGlobalContext(), output, |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1079 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
|
1080 } |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1081 } |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1082 |
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1083 |
6 | 1084 |
8
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1085 void EncodeNifti(OrthancPluginRestOutput* output, |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1086 const char* url, |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1087 const OrthancPluginHttpRequest* request) |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1088 { |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1089 static const char* const KEY_NIFTI = "Nifti"; |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1090 static const char* const KEY_RESOLUTION = "Resolution"; |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1091 static const char* const KEY_PARENT_STUDY = "ParentStudy"; |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1092 static const char* const KEY_SMOOTH = "Smooth"; |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1093 |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1094 if (request->method != OrthancPluginHttpMethod_Post) |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1095 { |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1096 OrthancPluginSendMethodNotAllowed(OrthancPlugins::GetGlobalContext(), output, "POST"); |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1097 return; |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1098 } |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1099 |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1100 Json::Value body; |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1101 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
|
1102 { |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1103 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadRequest); |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1104 } |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1105 |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1106 std::string mime, nifti; |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1107 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
|
1108 { |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1109 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
|
1110 } |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1111 |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1112 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
|
1113 const bool smooth = (body.isMember(KEY_SMOOTH) ? |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1114 Orthanc::SerializationToolbox::ReadBoolean(body, KEY_SMOOTH) : |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1115 true /* smooth by default */); |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1116 const unsigned int resolution = (body.isMember(KEY_RESOLUTION) ? |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1117 Orthanc::SerializationToolbox::ReadUnsignedInteger(body, KEY_RESOLUTION) : |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1118 256 /* default value */); |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1119 |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1120 vtkNew<vtkImageData> volume; |
32 | 1121 VTKToolbox::LoadNifti(volume.Get(), nifti); |
8
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1122 |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1123 std::string stl; |
32 | 1124 if (!VTKToolbox::EncodeVolume(stl, volume.Get(), resolution, smooth)) |
8
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1125 { |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1126 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
|
1127 } |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1128 else |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1129 { |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1130 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
|
1131 |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1132 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
|
1133 |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1134 Json::Value answer; |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1135 CallCreateDicom(answer, stl, body, parentStudy, title, frameOfReferenceUid, title); |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1136 |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1137 std::string s = answer.toStyledString(); |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1138 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
|
1139 } |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1140 } |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1141 |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1142 |
0
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1143 extern "C" |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1144 { |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1145 ORTHANC_PLUGINS_API int32_t OrthancPluginInitialize(OrthancPluginContext* context) |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1146 { |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1147 OrthancPlugins::SetGlobalContext(context); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1148 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1149 /* Check the version of the Orthanc core */ |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1150 if (OrthancPluginCheckVersion(OrthancPlugins::GetGlobalContext()) == 0) |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1151 { |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1152 char info[1024]; |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1153 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
|
1154 OrthancPlugins::GetGlobalContext()->orthancVersion, |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1155 ORTHANC_PLUGINS_MINIMAL_MAJOR_NUMBER, |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1156 ORTHANC_PLUGINS_MINIMAL_MINOR_NUMBER, |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1157 ORTHANC_PLUGINS_MINIMAL_REVISION_NUMBER); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1158 OrthancPluginLogError(OrthancPlugins::GetGlobalContext(), info); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1159 return -1; |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1160 } |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1161 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1162 #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
|
1163 Orthanc::Logging::InitializePluginContext(context); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1164 #else |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1165 Orthanc::Logging::Initialize(context); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1166 #endif |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1167 |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1168 Orthanc::InitializeFramework("", false); |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1169 |
0
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1170 hasCreateDicomStl_ = OrthancPlugins::CheckMinimalOrthancVersion(1, 12, 1); |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1171 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1172 if (!hasCreateDicomStl_) |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1173 { |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1174 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
|
1175 << ") 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
|
1176 } |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1177 |
25
5fc3ea6acbd5
fix deprecated calls
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
24
diff
changeset
|
1178 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
|
1179 |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1180 OrthancPlugins::RegisterRestCallback<ServeFile>("/stl/app/(.*)", true); |
3
0fb06c6a6c87
added route "/instances/{id}/stl"
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
2
diff
changeset
|
1181 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
|
1182 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
|
1183 |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1184 if (hasCreateDicomStl_) |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1185 { |
8
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1186 OrthancPlugins::RegisterRestCallback<EncodeStructureSet>("/stl/encode-rtstruct", true); |
d1267c6c33e1
added route /stl/encode-nifti
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
7
diff
changeset
|
1187 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
|
1188 } |
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1189 |
31
ab231760799d
added button to import STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
30
diff
changeset
|
1190 OrthancPlugins::OrthancConfiguration globalConfiguration; |
ab231760799d
added button to import STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
30
diff
changeset
|
1191 OrthancPlugins::OrthancConfiguration configuration; |
ab231760799d
added button to import STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
30
diff
changeset
|
1192 globalConfiguration.GetSection(configuration, "STL"); |
ab231760799d
added button to import STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
30
diff
changeset
|
1193 |
2
2bdb9acb7dcf
added STL viewer button at instance level too
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1
diff
changeset
|
1194 // 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
|
1195 std::string explorer; |
2bdb9acb7dcf
added STL viewer button at instance level too
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1
diff
changeset
|
1196 |
2bdb9acb7dcf
added STL viewer button at instance level too
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1
diff
changeset
|
1197 { |
2bdb9acb7dcf
added STL viewer button at instance level too
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1
diff
changeset
|
1198 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
|
1199 |
2bdb9acb7dcf
added STL viewer button at instance level too
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1
diff
changeset
|
1200 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
|
1201 dictionary["HAS_CREATE_DICOM_STL"] = (hasCreateDicomStl_ ? "true" : "false"); |
31
ab231760799d
added button to import STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
30
diff
changeset
|
1202 dictionary["SHOW_NIFTI_BUTTON"] = (configuration.GetBooleanValue("NIfTI", false) ? "true" : "false"); |
2
2bdb9acb7dcf
added STL viewer button at instance level too
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1
diff
changeset
|
1203 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
|
1204 |
25
5fc3ea6acbd5
fix deprecated calls
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
24
diff
changeset
|
1205 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
|
1206 } |
2bdb9acb7dcf
added STL viewer button at instance level too
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1
diff
changeset
|
1207 |
0
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1208 return 0; |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1209 } |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1210 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1211 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1212 ORTHANC_PLUGINS_API void OrthancPluginFinalize() |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1213 { |
1
0f03a8a0bd6f
encoding of RT-STRUCT as STL
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
1214 Orthanc::FinalizeFramework(); |
0
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1215 } |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1216 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1217 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1218 ORTHANC_PLUGINS_API const char* OrthancPluginGetName() |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1219 { |
25
5fc3ea6acbd5
fix deprecated calls
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
24
diff
changeset
|
1220 return ORTHANC_PLUGIN_NAME; |
0
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1221 } |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1222 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1223 |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1224 ORTHANC_PLUGINS_API const char* OrthancPluginGetVersion() |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1225 { |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1226 return ORTHANC_STL_VERSION; |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1227 } |
4e889a8e8be2
initial commit of the viewer plugin
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff
changeset
|
1228 } |