Mercurial > hg > orthanc-stone
annotate Framework/Layers/SeriesFrameRendererFactory.cpp @ 40:7207a407bcd8
shared copyright with osimis
author | Sebastien Jodogne <s.jodogne@gmail.com> |
---|---|
date | Wed, 04 Jan 2017 16:37:42 +0100 |
parents | 517c46f527cd |
children | 28956ed68280 |
rev | line source |
---|---|
0 | 1 /** |
2 * Stone of Orthanc | |
3 * Copyright (C) 2012-2016 Sebastien Jodogne, Medical Physics | |
4 * Department, University Hospital of Liege, Belgium | |
40
7207a407bcd8
shared copyright with osimis
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
32
diff
changeset
|
5 * Copyright (C) 2017 Osimis, Belgium |
0 | 6 * |
7 * This program is free software: you can redistribute it and/or | |
8 * modify it under the terms of the GNU General Public License as | |
9 * published by the Free Software Foundation, either version 3 of the | |
10 * License, or (at your option) any later version. | |
11 * | |
12 * In addition, as a special exception, the copyright holders of this | |
13 * program give permission to link the code of its release with the | |
14 * OpenSSL project's "OpenSSL" library (or with modified versions of it | |
15 * that use the same license as the "OpenSSL" library), and distribute | |
16 * the linked executables. You must obey the GNU General Public License | |
17 * in all respects for all of the code used other than "OpenSSL". If you | |
18 * modify file(s) with this exception, you may extend this exception to | |
19 * your version of the file(s), but you are not obligated to do so. If | |
20 * you do not wish to do so, delete this exception statement from your | |
21 * version. If you delete this exception statement from all source files | |
22 * in the program, then also delete it here. | |
23 * | |
24 * This program is distributed in the hope that it will be useful, but | |
25 * WITHOUT ANY WARRANTY; without even the implied warranty of | |
26 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU | |
27 * General Public License for more details. | |
28 * | |
29 * You should have received a copy of the GNU General Public License | |
30 * along with this program. If not, see <http://www.gnu.org/licenses/>. | |
31 **/ | |
32 | |
33 | |
34 #include "SeriesFrameRendererFactory.h" | |
35 | |
36 #include "FrameRenderer.h" | |
16 | 37 #include "../../Resources/Orthanc/Core/OrthancException.h" |
32 | 38 #include "../../Resources/Orthanc/Core/Logging.h" |
39 #include "../../Resources/Orthanc/Core/Toolbox.h" | |
40 #include "../../Resources/Orthanc/Plugins/Samples/Common/OrthancPluginException.h" | |
41 #include "../../Resources/Orthanc/Plugins/Samples/Common/DicomDatasetReader.h" | |
42 | |
0 | 43 |
44 namespace OrthancStone | |
45 { | |
46 void SeriesFrameRendererFactory::ReadCurrentFrameDataset(size_t frame) | |
47 { | |
48 if (currentDataset_.get() != NULL && | |
49 (fast_ || currentFrame_ == frame)) | |
50 { | |
51 // The frame has not changed since the previous call, no need to | |
52 // update the DICOM dataset | |
53 return; | |
54 } | |
55 | |
56 currentDataset_.reset(loader_->DownloadDicom(frame)); | |
57 currentFrame_ = frame; | |
58 | |
59 if (currentDataset_.get() == NULL) | |
60 { | |
61 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError); | |
62 } | |
63 } | |
64 | |
65 | |
66 void SeriesFrameRendererFactory::GetCurrentPixelSpacing(double& spacingX, | |
67 double& spacingY) const | |
68 { | |
69 if (currentDataset_.get() == NULL) | |
70 { | |
71 // There was no previous call "ReadCurrentFrameDataset()" | |
72 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError); | |
73 } | |
74 | |
32 | 75 GeometryToolbox::GetPixelSpacing(spacingX, spacingY, *currentDataset_); |
0 | 76 } |
77 | |
78 | |
79 double SeriesFrameRendererFactory::GetCurrentSliceThickness() const | |
80 { | |
81 if (currentDataset_.get() == NULL) | |
82 { | |
83 // There was no previous call "ReadCurrentFrameDataset()" | |
84 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError); | |
85 } | |
32 | 86 |
87 try | |
0 | 88 { |
32 | 89 OrthancPlugins::DicomDatasetReader reader(*currentDataset_); |
90 | |
91 double thickness; | |
92 if (reader.GetDoubleValue(thickness, OrthancPlugins::DICOM_TAG_SLICE_THICKNESS)) | |
93 { | |
94 return thickness; | |
95 } | |
0 | 96 } |
32 | 97 catch (ORTHANC_PLUGINS_EXCEPTION_CLASS& e) |
0 | 98 { |
99 } | |
32 | 100 |
101 // Some arbitrary large slice thickness | |
102 return std::numeric_limits<double>::infinity(); | |
0 | 103 } |
104 | |
105 | |
106 SeriesFrameRendererFactory::SeriesFrameRendererFactory(ISeriesLoader* loader, // Takes ownership | |
107 bool fast) : | |
108 loader_(loader), | |
109 currentFrame_(0), | |
110 fast_(fast) | |
111 { | |
112 if (loader == NULL) | |
113 { | |
114 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError); | |
115 } | |
116 } | |
117 | |
118 | |
119 bool SeriesFrameRendererFactory::GetExtent(double& x1, | |
120 double& y1, | |
121 double& x2, | |
122 double& y2, | |
123 const SliceGeometry& viewportSlice) | |
124 { | |
125 if (currentDataset_.get() == NULL) | |
126 { | |
127 // There has been no previous call to | |
128 // "CreateLayerRenderer". Read some arbitrary DICOM frame, the | |
129 // one at the middle of the series. | |
130 unsigned int depth = loader_->GetGeometry().GetSliceCount(); | |
131 ReadCurrentFrameDataset(depth / 2); | |
132 } | |
133 | |
134 double spacingX, spacingY; | |
135 GetCurrentPixelSpacing(spacingX, spacingY); | |
136 | |
137 return FrameRenderer::ComputeFrameExtent(x1, y1, x2, y2, | |
138 viewportSlice, | |
139 loader_->GetGeometry().GetSlice(0), | |
140 loader_->GetWidth(), | |
141 loader_->GetHeight(), | |
142 spacingX, spacingY); | |
143 } | |
144 | |
145 | |
146 ILayerRenderer* SeriesFrameRendererFactory::CreateLayerRenderer(const SliceGeometry& viewportSlice) | |
147 { | |
148 size_t closest; | |
149 double distance; | |
150 | |
151 bool isOpposite; | |
152 if (!GeometryToolbox::IsParallelOrOpposite(isOpposite, loader_->GetGeometry().GetNormal(), viewportSlice.GetNormal()) || | |
153 !loader_->GetGeometry().ComputeClosestSlice(closest, distance, viewportSlice.GetOrigin())) | |
154 { | |
155 // Unable to compute the slice in the series that is the | |
156 // closest to the slice displayed by the viewport | |
157 return NULL; | |
158 } | |
159 | |
160 ReadCurrentFrameDataset(closest); | |
161 assert(currentDataset_.get() != NULL); | |
162 | |
163 double spacingX, spacingY; | |
164 GetCurrentPixelSpacing(spacingX, spacingY); | |
165 | |
166 if (distance <= GetCurrentSliceThickness() / 2.0) | |
167 { | |
168 SliceGeometry frameSlice(*currentDataset_); | |
169 return FrameRenderer::CreateRenderer(loader_->DownloadFrame(closest), | |
170 viewportSlice, | |
171 frameSlice, | |
172 *currentDataset_, | |
173 spacingX, spacingY, | |
174 true); | |
175 } | |
176 else | |
177 { | |
178 // The closest slice of the series is too far away from the | |
179 // slice displayed by the viewport | |
180 return NULL; | |
181 } | |
182 } | |
183 | |
184 | |
185 ISliceableVolume& SeriesFrameRendererFactory::GetSourceVolume() const | |
186 { | |
187 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls); | |
188 } | |
189 } |