Mercurial > hg > orthanc-wsi
annotate Framework/Inputs/HierarchicalTiff.cpp @ 334:f011fc199b6a default tip
back to mainline
author | Sebastien Jodogne <s.jodogne@gmail.com> |
---|---|
date | Fri, 18 Oct 2024 15:44:18 +0200 |
parents | 0683312e21ba |
children |
rev | line source |
---|---|
0 | 1 /** |
2 * Orthanc - A Lightweight, RESTful DICOM Store | |
3 * Copyright (C) 2012-2016 Sebastien Jodogne, Medical Physics | |
4 * Department, University Hospital of Liege, Belgium | |
312
0683312e21ba
updated copyright, as Orthanc Team now replaces Osimis
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
309
diff
changeset
|
5 * Copyright (C) 2017-2023 Osimis S.A., Belgium |
0683312e21ba
updated copyright, as Orthanc Team now replaces Osimis
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
309
diff
changeset
|
6 * Copyright (C) 2024-2024 Orthanc Team SRL, Belgium |
309
7020852a8fa9
updated year to 2024
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
297
diff
changeset
|
7 * Copyright (C) 2021-2024 Sebastien Jodogne, ICTEAM UCLouvain, Belgium |
0 | 8 * |
9 * This program is free software: you can redistribute it and/or | |
10 * modify it under the terms of the GNU Affero General Public License | |
11 * as published by the Free Software Foundation, either version 3 of | |
12 * the License, or (at your option) any later version. | |
13 * | |
14 * This program is distributed in the hope that it will be useful, but | |
15 * WITHOUT ANY WARRANTY; without even the implied warranty of | |
16 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU | |
17 * Affero General Public License for more details. | |
18 * | |
19 * You should have received a copy of the GNU Affero General Public License | |
20 * along with this program. If not, see <http://www.gnu.org/licenses/>. | |
21 **/ | |
22 | |
23 | |
16
7a88c614be04
preparing for precompiled headers
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
24 #include "../PrecompiledHeadersWSI.h" |
0 | 25 #include "HierarchicalTiff.h" |
26 | |
192 | 27 #include <Logging.h> |
28 #include <OrthancException.h> | |
0 | 29 |
30 #include <iostream> | |
31 #include <algorithm> | |
32 #include <cassert> | |
33 #include <string.h> | |
34 | |
35 namespace OrthancWSI | |
36 { | |
37 HierarchicalTiff::Level::Level(TIFF* tiff, | |
38 tdir_t directory, | |
39 unsigned int width, | |
40 unsigned int height) : | |
41 directory_(directory), | |
42 width_(width), | |
43 height_(height) | |
44 { | |
45 // Read the JPEG headers shared at that level, if any | |
46 uint8_t *tables = NULL; | |
47 uint32_t size; | |
48 if (TIFFGetField(tiff, TIFFTAG_JPEGTABLES, &size, &tables) && | |
49 size > 0 && | |
50 tables != NULL) | |
51 { | |
52 // Look for the EOI (end-of-image) tag == FF D9 | |
53 // https://en.wikipedia.org/wiki/JPEG_File_Interchange_Format | |
54 | |
55 bool found = false; | |
56 | |
57 for (size_t i = 0; i + 1 < size; i++) | |
58 { | |
59 if (tables[i] == 0xff && | |
60 tables[i + 1] == 0xd9) | |
61 { | |
62 headers_.assign(reinterpret_cast<const char*>(tables), i); | |
63 found = true; | |
64 } | |
65 } | |
66 | |
67 if (!found) | |
68 { | |
69 headers_.assign(reinterpret_cast<const char*>(tables), size); | |
70 } | |
71 } | |
72 } | |
73 | |
74 struct HierarchicalTiff::Comparator | |
75 { | |
76 bool operator() (const HierarchicalTiff::Level& a, | |
77 const HierarchicalTiff::Level& b) const | |
78 { | |
79 return a.width_ > b.width_; | |
80 } | |
81 }; | |
82 | |
83 | |
84 void HierarchicalTiff::CheckLevel(unsigned int level) const | |
85 { | |
86 if (level >= levels_.size()) | |
87 { | |
88 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange); | |
89 } | |
90 } | |
91 | |
92 | |
297 | 93 HierarchicalTiff::HierarchicalTiff(const std::string& path) : |
94 reader_(path), | |
95 tileWidth_(0), | |
96 tileHeight_(0) | |
0 | 97 { |
98 bool first = true; | |
99 tdir_t pos = 0; | |
100 | |
101 do | |
102 { | |
103 uint32_t w, h, tw, th; | |
104 ImageCompression compression; | |
105 Orthanc::PixelFormat pixelFormat; | |
166
f0dac1e8f736
access to photometric interpretation of source pyramids
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
157
diff
changeset
|
106 Orthanc::PhotometricInterpretation photometric; |
0 | 107 |
297 | 108 if (TIFFSetDirectory(reader_.GetTiff(), pos) && |
109 TIFFGetField(reader_.GetTiff(), TIFFTAG_IMAGEWIDTH, &w) && | |
110 TIFFGetField(reader_.GetTiff(), TIFFTAG_IMAGELENGTH, &h) && | |
111 TIFFGetField(reader_.GetTiff(), TIFFTAG_TILEWIDTH, &tw) && | |
112 TIFFGetField(reader_.GetTiff(), TIFFTAG_TILELENGTH, &th) && | |
0 | 113 w > 0 && |
114 h > 0 && | |
115 tw > 0 && | |
116 th > 0 && | |
297 | 117 reader_.GetCurrentDirectoryInformation(compression, pixelFormat, photometric)) |
0 | 118 { |
119 if (first) | |
120 { | |
121 tileWidth_ = tw; | |
122 tileHeight_ = th; | |
123 compression_ = compression; | |
124 pixelFormat_ = pixelFormat; | |
166
f0dac1e8f736
access to photometric interpretation of source pyramids
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
157
diff
changeset
|
125 photometric_ = photometric; |
0 | 126 first = false; |
127 } | |
128 else if (tw != tileWidth_ || | |
129 th != tileHeight_ || | |
130 compression_ != compression || | |
166
f0dac1e8f736
access to photometric interpretation of source pyramids
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
157
diff
changeset
|
131 pixelFormat_ != pixelFormat || |
f0dac1e8f736
access to photometric interpretation of source pyramids
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
157
diff
changeset
|
132 photometric_ != photometric) |
0 | 133 { |
297 | 134 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat, |
135 "The tile size or compression of the TIFF file varies along levels, this is not supported"); | |
0 | 136 } |
137 | |
297 | 138 levels_.push_back(Level(reader_.GetTiff(), pos, w, h)); |
0 | 139 } |
140 | |
141 pos++; | |
142 } | |
297 | 143 while (TIFFReadDirectory(reader_.GetTiff())); |
0 | 144 |
145 if (levels_.size() == 0) | |
146 { | |
297 | 147 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat, "This is not a tiled TIFF image"); |
0 | 148 } |
149 | |
150 std::sort(levels_.begin(), levels_.end(), Comparator()); | |
151 } | |
152 | |
153 | |
154 unsigned int HierarchicalTiff::GetLevelWidth(unsigned int level) const | |
155 { | |
156 CheckLevel(level); | |
157 return levels_[level].width_; | |
158 } | |
159 | |
160 | |
161 unsigned int HierarchicalTiff::GetLevelHeight(unsigned int level) const | |
162 { | |
163 CheckLevel(level); | |
164 return levels_[level].height_; | |
165 } | |
166 | |
167 | |
168 bool HierarchicalTiff::ReadRawTile(std::string& tile, | |
57
91fc9583b2de
big refactoring to support sparse tiling
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
169 ImageCompression& compression, |
0 | 170 unsigned int level, |
171 unsigned int tileX, | |
172 unsigned int tileY) | |
173 { | |
174 boost::mutex::scoped_lock lock(mutex_); | |
175 | |
176 CheckLevel(level); | |
177 | |
57
91fc9583b2de
big refactoring to support sparse tiling
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
178 compression = compression_; |
91fc9583b2de
big refactoring to support sparse tiling
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
179 |
0 | 180 // Make the TIFF context point to the level of interest |
297 | 181 if (!TIFFSetDirectory(reader_.GetTiff(), levels_[level].directory_)) |
0 | 182 { |
183 throw Orthanc::OrthancException(Orthanc::ErrorCode_CorruptedFile); | |
184 } | |
185 | |
186 // Get the index of the tile | |
297 | 187 ttile_t index = TIFFComputeTile(reader_.GetTiff(), tileX * tileWidth_, tileY * tileHeight_, 0 /*z*/, 0 /*sample*/); |
0 | 188 |
189 // Read the raw tile | |
190 toff_t *sizes; | |
297 | 191 if (!TIFFGetField(reader_.GetTiff(), TIFFTAG_TILEBYTECOUNTS, &sizes)) |
0 | 192 { |
193 throw Orthanc::OrthancException(Orthanc::ErrorCode_CorruptedFile); | |
194 } | |
195 | |
196 std::string raw; | |
197 raw.resize(sizes[index]); | |
198 | |
297 | 199 tsize_t read = TIFFReadRawTile(reader_.GetTiff(), index, &raw[0], raw.size()); |
0 | 200 if (read != static_cast<tsize_t>(sizes[index])) |
201 { | |
202 throw Orthanc::OrthancException(Orthanc::ErrorCode_CorruptedFile); | |
203 } | |
204 | |
205 const std::string& headers = levels_[level].headers_; | |
206 | |
207 // Possibly prepend the raw tile with the shared JPEG headers | |
208 if (headers.empty() || | |
209 compression_ != ImageCompression_Jpeg) | |
210 { | |
211 tile.swap(raw); // Same as "tile.assign(raw)", but optimizes memory | |
212 } | |
213 else | |
214 { | |
215 assert(compression_ == ImageCompression_Jpeg); | |
216 | |
217 // Check that the raw JPEG tile starts with the SOI (start-of-image) tag == FF D8 | |
218 if (raw.size() < 2 || | |
219 static_cast<uint8_t>(raw[0]) != 0xff || | |
220 static_cast<uint8_t>(raw[1]) != 0xd8) | |
221 { | |
222 throw Orthanc::OrthancException(Orthanc::ErrorCode_CorruptedFile); | |
223 } | |
224 | |
218
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
225 if (photometric_ == Orthanc::PhotometricInterpretation_RGB && |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
226 pixelFormat_ == Orthanc::PixelFormat_RGB24) |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
227 { |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
228 /** |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
229 * Insert an Adobe APP14 marker with the "transform" flag set |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
230 * to value 0, which indicates to the JPEG decoder that |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
231 * "3-channel images are assumed to be RGB". Section 18 of |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
232 * "Supporting the DCT Filters in PostScript Level 2 - |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
233 * Technical Note #5116": |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
234 * https://stackoverflow.com/a/9658206/881731 |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
235 * https://docs.oracle.com/javase/6/docs/api/javax/imageio/metadata/doc-files/jpeg_metadata.html |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
236 * https://www.pdfa.org/wp-content/uploads/2020/07/5116.DCT_Filter.pdf |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
237 **/ |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
238 static const uint8_t APP14[] = { |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
239 0xff, 0xee, /* JPEG Marker for Adobe segment: http://www.ozhiker.com/electronics/pjmt/jpeg_info/app_segments.html */ |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
240 0x00, 0x0e, /* Length (without the JPEG marker) == 0x0e == 14 bytes */ |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
241 0x41, 0x64, 0x6f, 0x62, 0x65, /* "Adobe" string in ASCII */ |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
242 0x00, 0x64, /* Version == Two-byte DCTEncode/DCTDecode version number == 0x64 */ |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
243 0x80, 0x00, /* Two-byte "flags0" 0x8000 bit: Encoder used Blend=1 downsampling */ |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
244 0x00, 0x00, /* Two-byte "flags1": Set to zero */ |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
245 0x00 /* One-byte color transform code == 0 <== This is the important one */ |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
246 }; |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
247 assert(sizeof(APP14) == 16); |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
248 |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
249 tile.resize(headers.size() + sizeof(APP14) + raw.size() - 2); |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
250 memcpy(&tile[0], &headers[0], headers.size()); |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
251 memcpy(&tile[0] + headers.size(), APP14, sizeof(APP14)); |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
252 memcpy(&tile[0] + headers.size() + sizeof(APP14), &raw[2], raw.size() - 2); |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
253 } |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
254 else |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
255 { |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
256 tile.resize(headers.size() + raw.size() - 2); |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
257 memcpy(&tile[0], &headers[0], headers.size()); |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
258 memcpy(&tile[0] + headers.size(), &raw[2], raw.size() - 2); |
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
259 } |
0 | 260 } |
218
c5a8b46c4cba
Fix colorspace of TIFF containing JPEG with RGB photometric interpretation
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
214
diff
changeset
|
261 |
0 | 262 return true; |
263 } | |
264 } |