Mercurial > hg > orthanc-wsi

annotate Framework/Inputs/OpenSlidePyramid.cpp @ 317:f611fb47d0e8

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
reuse base class members in OpenSlidePyramid
author Sebastien Jodogne <s.jodogne@gmail.com>
date Wed, 11 Sep 2024 13:43:39 +0200
parents 072968f00d26
children 8ad12abde290
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
1 /**
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
2 * Orthanc - A Lightweight, RESTful DICOM Store
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
3 * Copyright (C) 2012-2016 Sebastien Jodogne, Medical Physics
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
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: 279
diff changeset
7 * Copyright (C) 2021-2024 Sebastien Jodogne, ICTEAM UCLouvain, Belgium
0
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
8 *
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
9 * This program is free software: you can redistribute it and/or
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
10 * modify it under the terms of the GNU Affero General Public License
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
11 * as published by the Free Software Foundation, either version 3 of
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
12 * the License, or (at your option) any later version.
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
13 *
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
14 * This program is distributed in the hope that it will be useful, but
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
15 * WITHOUT ANY WARRANTY; without even the implied warranty of
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
16 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
17 * Affero General Public License for more details.
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
18 *
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
19 * You should have received a copy of the GNU Affero General Public License
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
20 * along with this program. If not, see <http://www.gnu.org/licenses/>.
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
21 **/
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
22
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
23
16
7a88c614be04 preparing for precompiled headers
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 0
diff changeset
24 #include "../PrecompiledHeadersWSI.h"
0
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
25 #include "OpenSlidePyramid.h"
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
26
199
a1c265cb2174 replacing deprecated std::auto_ptr by std::unique_ptr
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 192
diff changeset
27 #include <Compatibility.h> // For std::unique_ptr
192
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 175
diff changeset
28 #include <Images/ImageProcessing.h>
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 175
diff changeset
29 #include <OrthancException.h>
279
77afef2cf64b automated extraction of the imaged volume if using OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 254
diff changeset
30 #include <SerializationToolbox.h>
192
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 175
diff changeset
31 #include <Logging.h>
140
a0f9a3df1110 resort to Orthanc framework
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 115
diff changeset
32
315
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
33 #include <boost/math/special_functions/round.hpp>
140
a0f9a3df1110 resort to Orthanc framework
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 115
diff changeset
34 #include <memory>
0
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
35
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
36 namespace OrthancWSI
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
37 {
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
38 void OpenSlidePyramid::ReadRegion(Orthanc::ImageAccessor& target,
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
39 unsigned int level,
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
40 unsigned int x,
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
41 unsigned int y)
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
42 {
199
a1c265cb2174 replacing deprecated std::auto_ptr by std::unique_ptr
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 192
diff changeset
43 std::unique_ptr<Orthanc::ImageAccessor> source(image_.ReadRegion(level, x, y, target.GetWidth(), target.GetHeight()));
315
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
44
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
45 if (target.GetWidth() != source->GetWidth() ||
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
46 target.GetHeight() != source->GetHeight())
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
47 {
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
48 throw Orthanc::OrthancException(Orthanc::ErrorCode_IncompatibleImageSize);
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
49 }
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
50
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
51 const unsigned int width = source->GetWidth();
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
52 const unsigned int height = source->GetHeight();
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
53
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
54 if (target.GetFormat() == Orthanc::PixelFormat_RGB24 &&
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
55 source->GetFormat() == Orthanc::PixelFormat_BGRA32)
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
56 {
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
57 // Implements alpha blending: https://en.wikipedia.org/wiki/Alpha_compositing#Alpha_blending
317
f611fb47d0e8 reuse base class members in OpenSlidePyramid
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 315
diff changeset
58
f611fb47d0e8 reuse base class members in OpenSlidePyramid
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 315
diff changeset
59 uint8_t backgroundRed, backgroundGreen, backgroundBlue;
f611fb47d0e8 reuse base class members in OpenSlidePyramid
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 315
diff changeset
60 GetBackgroundColor(backgroundRed, backgroundGreen, backgroundBlue);
f611fb47d0e8 reuse base class members in OpenSlidePyramid
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 315
diff changeset
61
315
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
62 for (unsigned int y = 0; y < height; y++)
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
63 {
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
64 const uint8_t* p = reinterpret_cast<const uint8_t*>(source->GetConstRow(y));
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
65 uint8_t* q = reinterpret_cast<uint8_t*>(target.GetRow(y));
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
66 for (unsigned int x = 0; x < width; x++)
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
67 {
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
68 /**
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
69 Alpha blending using integer arithmetics only (16 bits avoids overflows)
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
70
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
71 p = (1 - alpha) * background + alpha * value
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
72 <=> p = (1 - p[3] / 255) * background + p[3] / 255 * value
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
73 <=> p = ((255 - p[3]) * background + p[3] * value) / 255
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
74
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
75 **/
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
76
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
77 uint16_t alpha = p[3];
317
f611fb47d0e8 reuse base class members in OpenSlidePyramid
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 315
diff changeset
78 q[0] = static_cast<uint8_t>(((255 - alpha) * backgroundRed + alpha * p[2]) / 255);
f611fb47d0e8 reuse base class members in OpenSlidePyramid
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 315
diff changeset
79 q[1] = static_cast<uint8_t>(((255 - alpha) * backgroundGreen + alpha * p[1]) / 255);
f611fb47d0e8 reuse base class members in OpenSlidePyramid
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 315
diff changeset
80 q[2] = static_cast<uint8_t>(((255 - alpha) * backgroundBlue + alpha * p[0]) / 255);
315
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
81
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
82 p += 4;
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
83 q += 3;
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
84 }
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
85 }
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
86 }
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
87 else
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
88 {
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
89 Orthanc::ImageProcessing::Convert(target, *source);
072968f00d26 support of transparency in OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 312
diff changeset
90 }
0
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
91 }
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
92
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
93
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
94 OpenSlidePyramid::OpenSlidePyramid(const std::string& path,
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
95 unsigned int tileWidth,
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
96 unsigned int tileHeight) :
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
97 image_(path),
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
98 tileWidth_(tileWidth),
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
99 tileHeight_(tileHeight)
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
100 {
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
101 }
279
77afef2cf64b automated extraction of the imaged volume if using OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 254
diff changeset
102
77afef2cf64b automated extraction of the imaged volume if using OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 254
diff changeset
103
77afef2cf64b automated extraction of the imaged volume if using OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 254
diff changeset
104 bool OpenSlidePyramid::LookupImagedVolumeSize(float& width,
77afef2cf64b automated extraction of the imaged volume if using OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 254
diff changeset
105 float& height) const
77afef2cf64b automated extraction of the imaged volume if using OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 254
diff changeset
106 {
77afef2cf64b automated extraction of the imaged volume if using OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 254
diff changeset
107 std::string s;
77afef2cf64b automated extraction of the imaged volume if using OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 254
diff changeset
108 double mppx;
77afef2cf64b automated extraction of the imaged volume if using OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 254
diff changeset
109 double mppy;
77afef2cf64b automated extraction of the imaged volume if using OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 254
diff changeset
110
77afef2cf64b automated extraction of the imaged volume if using OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 254
diff changeset
111 if (image_.LookupProperty(s, "openslide.mpp-x") &&
77afef2cf64b automated extraction of the imaged volume if using OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 254
diff changeset
112 Orthanc::SerializationToolbox::ParseDouble(mppx, s) &&
77afef2cf64b automated extraction of the imaged volume if using OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 254
diff changeset
113 image_.LookupProperty(s, "openslide.mpp-y") &&
77afef2cf64b automated extraction of the imaged volume if using OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 254
diff changeset
114 Orthanc::SerializationToolbox::ParseDouble(mppy, s))
77afef2cf64b automated extraction of the imaged volume if using OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 254
diff changeset
115 {
77afef2cf64b automated extraction of the imaged volume if using OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 254
diff changeset
116 // In the 2 lines below, remember to switch X/Y when going from physical to pixel coordinates!
77afef2cf64b automated extraction of the imaged volume if using OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 254
diff changeset
117 width = mppy / 1000.0 * static_cast<double>(image_.GetLevelHeight(0));
77afef2cf64b automated extraction of the imaged volume if using OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 254
diff changeset
118 height = mppx / 1000.0 * static_cast<double>(image_.GetLevelWidth(0));
77afef2cf64b automated extraction of the imaged volume if using OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 254
diff changeset
119 return true;
77afef2cf64b automated extraction of the imaged volume if using OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 254
diff changeset
120 }
77afef2cf64b automated extraction of the imaged volume if using OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 254
diff changeset
121 else
77afef2cf64b automated extraction of the imaged volume if using OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 254
diff changeset
122 {
77afef2cf64b automated extraction of the imaged volume if using OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 254
diff changeset
123 return false;
77afef2cf64b automated extraction of the imaged volume if using OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 254
diff changeset
124 }
77afef2cf64b automated extraction of the imaged volume if using OpenSlide
Sebastien Jodogne <s.jodogne@gmail.com>
parents: 254
diff changeset
125 }
0
4a7a53257c7d initial commit
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
diff changeset
126 }