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