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 } |