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