Mercurial > hg > orthanc-webviewer
annotate Orthanc/Core/Images/ImageAccessor.cpp @ 178:6547e6941b15
note
author | Sebastien Jodogne <s.jodogne@gmail.com> |
---|---|
date | Thu, 04 Jan 2018 09:49:10 +0100 |
parents | dd1ad819ca33 |
children |
rev | line source |
---|---|
78 | 1 /** |
2 * Orthanc - A Lightweight, RESTful DICOM Store | |
122 | 3 * Copyright (C) 2012-2016 Sebastien Jodogne, Medical Physics |
78 | 4 * Department, University Hospital of Liege, Belgium |
175 | 5 * Copyright (C) 2017-2018 Osimis S.A., Belgium |
78 | 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 "../PrecompiledHeaders.h" | |
35 #include "ImageAccessor.h" | |
36 | |
37 #include "../Logging.h" | |
38 #include "../OrthancException.h" | |
39 #include "../ChunkedBuffer.h" | |
40 | |
41 #include <stdint.h> | |
42 #include <cassert> | |
43 #include <boost/lexical_cast.hpp> | |
44 | |
45 | |
46 | |
47 namespace Orthanc | |
48 { | |
49 template <typename PixelType> | |
50 static void ToMatlabStringInternal(ChunkedBuffer& target, | |
51 const ImageAccessor& source) | |
52 { | |
53 target.AddChunk("double([ "); | |
54 | |
55 for (unsigned int y = 0; y < source.GetHeight(); y++) | |
56 { | |
57 const PixelType* p = reinterpret_cast<const PixelType*>(source.GetConstRow(y)); | |
58 | |
59 std::string s; | |
60 if (y > 0) | |
61 { | |
62 s = "; "; | |
63 } | |
64 | |
65 s.reserve(source.GetWidth() * 8); | |
66 | |
67 for (unsigned int x = 0; x < source.GetWidth(); x++, p++) | |
68 { | |
129 | 69 s += boost::lexical_cast<std::string>(static_cast<double>(*p)) + " "; |
78 | 70 } |
71 | |
72 target.AddChunk(s); | |
73 } | |
74 | |
75 target.AddChunk("])"); | |
76 } | |
77 | |
78 | |
79 static void RGB24ToMatlabString(ChunkedBuffer& target, | |
80 const ImageAccessor& source) | |
81 { | |
82 assert(source.GetFormat() == PixelFormat_RGB24); | |
83 | |
84 target.AddChunk("double(permute(reshape([ "); | |
85 | |
86 for (unsigned int y = 0; y < source.GetHeight(); y++) | |
87 { | |
88 const uint8_t* p = reinterpret_cast<const uint8_t*>(source.GetConstRow(y)); | |
89 | |
90 std::string s; | |
91 s.reserve(source.GetWidth() * 3 * 8); | |
92 | |
93 for (unsigned int x = 0; x < 3 * source.GetWidth(); x++, p++) | |
94 { | |
95 s += boost::lexical_cast<std::string>(static_cast<int>(*p)) + " "; | |
96 } | |
97 | |
98 target.AddChunk(s); | |
99 } | |
100 | |
101 target.AddChunk("], [ 3 " + boost::lexical_cast<std::string>(source.GetHeight()) + | |
102 " " + boost::lexical_cast<std::string>(source.GetWidth()) + " ]), [ 3 2 1 ]))"); | |
103 } | |
104 | |
105 | |
106 void* ImageAccessor::GetBuffer() const | |
107 { | |
108 if (readOnly_) | |
109 { | |
110 #if ORTHANC_ENABLE_LOGGING == 1 | |
111 LOG(ERROR) << "Trying to write on a read-only image"; | |
112 #endif | |
113 | |
114 throw OrthancException(ErrorCode_ReadOnly); | |
115 } | |
116 | |
117 return buffer_; | |
118 } | |
119 | |
120 | |
121 const void* ImageAccessor::GetConstRow(unsigned int y) const | |
122 { | |
123 if (buffer_ != NULL) | |
124 { | |
128
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
125 return buffer_ + y * pitch_; |
78 | 126 } |
127 else | |
128 { | |
129 return NULL; | |
130 } | |
131 } | |
132 | |
133 | |
134 void* ImageAccessor::GetRow(unsigned int y) const | |
135 { | |
136 if (readOnly_) | |
137 { | |
138 #if ORTHANC_ENABLE_LOGGING == 1 | |
139 LOG(ERROR) << "Trying to write on a read-only image"; | |
140 #endif | |
141 | |
142 throw OrthancException(ErrorCode_ReadOnly); | |
143 } | |
144 | |
145 if (buffer_ != NULL) | |
146 { | |
128
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
147 return buffer_ + y * pitch_; |
78 | 148 } |
149 else | |
150 { | |
151 return NULL; | |
152 } | |
153 } | |
154 | |
155 | |
156 void ImageAccessor::AssignEmpty(PixelFormat format) | |
157 { | |
158 readOnly_ = false; | |
159 format_ = format; | |
160 width_ = 0; | |
161 height_ = 0; | |
162 pitch_ = 0; | |
163 buffer_ = NULL; | |
164 } | |
165 | |
166 | |
167 void ImageAccessor::AssignReadOnly(PixelFormat format, | |
168 unsigned int width, | |
169 unsigned int height, | |
170 unsigned int pitch, | |
171 const void *buffer) | |
172 { | |
173 readOnly_ = true; | |
174 format_ = format; | |
175 width_ = width; | |
176 height_ = height; | |
177 pitch_ = pitch; | |
128
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
178 buffer_ = reinterpret_cast<uint8_t*>(const_cast<void*>(buffer)); |
78 | 179 |
128
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
180 if (GetBytesPerPixel() * width_ > pitch_) |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
181 { |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
182 throw OrthancException(ErrorCode_ParameterOutOfRange); |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
183 } |
78 | 184 } |
185 | |
186 | |
187 void ImageAccessor::AssignWritable(PixelFormat format, | |
188 unsigned int width, | |
189 unsigned int height, | |
190 unsigned int pitch, | |
191 void *buffer) | |
192 { | |
193 readOnly_ = false; | |
194 format_ = format; | |
195 width_ = width; | |
196 height_ = height; | |
197 pitch_ = pitch; | |
128
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
198 buffer_ = reinterpret_cast<uint8_t*>(buffer); |
78 | 199 |
128
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
200 if (GetBytesPerPixel() * width_ > pitch_) |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
201 { |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
202 throw OrthancException(ErrorCode_ParameterOutOfRange); |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
203 } |
78 | 204 } |
205 | |
206 | |
207 void ImageAccessor::ToMatlabString(std::string& target) const | |
208 { | |
209 ChunkedBuffer buffer; | |
210 | |
211 switch (GetFormat()) | |
212 { | |
213 case PixelFormat_Grayscale8: | |
214 ToMatlabStringInternal<uint8_t>(buffer, *this); | |
215 break; | |
216 | |
217 case PixelFormat_Grayscale16: | |
218 ToMatlabStringInternal<uint16_t>(buffer, *this); | |
219 break; | |
220 | |
172 | 221 case PixelFormat_Grayscale32: |
222 ToMatlabStringInternal<uint32_t>(buffer, *this); | |
223 break; | |
224 | |
78 | 225 case PixelFormat_SignedGrayscale16: |
226 ToMatlabStringInternal<int16_t>(buffer, *this); | |
227 break; | |
228 | |
129 | 229 case PixelFormat_Float32: |
230 ToMatlabStringInternal<float>(buffer, *this); | |
231 break; | |
232 | |
78 | 233 case PixelFormat_RGB24: |
234 RGB24ToMatlabString(buffer, *this); | |
235 break; | |
236 | |
237 default: | |
238 throw OrthancException(ErrorCode_NotImplemented); | |
239 } | |
240 | |
241 buffer.Flatten(target); | |
242 } | |
243 | |
128
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
244 |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
245 |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
246 ImageAccessor ImageAccessor::GetRegion(unsigned int x, |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
247 unsigned int y, |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
248 unsigned int width, |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
249 unsigned int height) const |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
250 { |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
251 if (x + width > width_ || |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
252 y + height > height_) |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
253 { |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
254 throw OrthancException(ErrorCode_ParameterOutOfRange); |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
255 } |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
256 |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
257 ImageAccessor result; |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
258 |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
259 if (width == 0 || |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
260 height == 0) |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
261 { |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
262 result.AssignWritable(format_, 0, 0, 0, NULL); |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
263 } |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
264 else |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
265 { |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
266 uint8_t* p = (buffer_ + |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
267 y * pitch_ + |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
268 x * GetBytesPerPixel()); |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
269 |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
270 if (readOnly_) |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
271 { |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
272 result.AssignReadOnly(format_, width, height, pitch_, p); |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
273 } |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
274 else |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
275 { |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
276 result.AssignWritable(format_, width, height, pitch_, p); |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
277 } |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
278 } |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
279 |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
280 return result; |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
281 } |
e8cfda4c8a2f
Sync + support of JPEG2000 lossless images with YBR_RCT photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
122
diff
changeset
|
282 |
145 | 283 |
284 void ImageAccessor::SetFormat(PixelFormat format) | |
285 { | |
286 if (readOnly_) | |
287 { | |
288 #if ORTHANC_ENABLE_LOGGING == 1 | |
289 LOG(ERROR) << "Trying to modify the format of a read-only image"; | |
290 #endif | |
291 throw OrthancException(ErrorCode_ReadOnly); | |
292 } | |
293 | |
294 if (::Orthanc::GetBytesPerPixel(format) != ::Orthanc::GetBytesPerPixel(format_)) | |
295 { | |
296 throw OrthancException(ErrorCode_IncompatibleImageFormat); | |
297 } | |
298 | |
299 format_ = format; | |
300 } | |
78 | 301 } |