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