Mercurial > hg > orthanc-wsi
annotate Framework/Inputs/DicomPyramidLevel.cpp @ 254:20a730889ae2
upgrade to 2023
| author | Sebastien Jodogne <s.jodogne@gmail.com> |
|---|---|
| date | Fri, 07 Jul 2023 17:38:07 +0200 |
| parents | 49f647ed1b4c |
| children | 7020852a8fa9 |
| 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 | |
| 254 | 5 * Copyright (C) 2017-2023 Osimis S.A., Belgium |
| 6 * Copyright (C) 2021-2023 Sebastien Jodogne, ICTEAM UCLouvain, Belgium | |
| 0 | 7 * |
| 8 * This program is free software: you can redistribute it and/or | |
| 9 * modify it under the terms of the GNU Affero General Public License | |
| 10 * as published by the Free Software Foundation, either version 3 of | |
| 11 * the License, or (at your option) any later version. | |
| 12 * | |
| 13 * This program is distributed in the hope that it will be useful, but | |
| 14 * WITHOUT ANY WARRANTY; without even the implied warranty of | |
| 15 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU | |
| 16 * Affero General Public License for more details. | |
| 17 * | |
| 18 * You should have received a copy of the GNU Affero General Public License | |
| 19 * along with this program. If not, see <http://www.gnu.org/licenses/>. | |
| 20 **/ | |
| 21 | |
| 22 | |
|
16
7a88c614be04
preparing for precompiled headers
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
0
diff
changeset
|
23 #include "../PrecompiledHeadersWSI.h" |
| 0 | 24 #include "DicomPyramidLevel.h" |
| 25 | |
| 192 | 26 #include <Logging.h> |
| 27 #include <OrthancException.h> | |
| 0 | 28 |
| 29 #include <boost/lexical_cast.hpp> | |
| 30 | |
| 31 namespace OrthancWSI | |
| 32 { | |
|
55
b6432a00b103
use of simpler data structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
33 DicomPyramidLevel::TileContent& DicomPyramidLevel::GetTileContent(unsigned int tileX, |
|
b6432a00b103
use of simpler data structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
34 unsigned int tileY) |
|
b6432a00b103
use of simpler data structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
35 { |
|
b6432a00b103
use of simpler data structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
36 if (tileX >= countTilesX_ || |
|
b6432a00b103
use of simpler data structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
37 tileY >= countTilesY_) |
|
b6432a00b103
use of simpler data structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
38 { |
|
b6432a00b103
use of simpler data structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
39 LOG(ERROR) << "Tile location (" << tileX << "," << tileY << ") is outside the image"; |
|
b6432a00b103
use of simpler data structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
40 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); |
|
b6432a00b103
use of simpler data structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
41 } |
|
b6432a00b103
use of simpler data structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
42 |
|
b6432a00b103
use of simpler data structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
43 return tiles_[tileY * countTilesX_ + tileX]; |
|
b6432a00b103
use of simpler data structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
44 } |
|
b6432a00b103
use of simpler data structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
45 |
|
56
83cd735c885d
speedup the loading of DICOM sources
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
55
diff
changeset
|
46 void DicomPyramidLevel::RegisterFrame(DicomPyramidInstance& instance, |
| 0 | 47 unsigned int frame) |
| 48 { | |
|
55
b6432a00b103
use of simpler data structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
49 unsigned int tileX = instance.GetFrameLocationX(frame); |
|
b6432a00b103
use of simpler data structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
50 unsigned int tileY = instance.GetFrameLocationY(frame); |
|
b6432a00b103
use of simpler data structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
51 TileContent& tile = GetTileContent(tileX, tileY); |
| 0 | 52 |
|
55
b6432a00b103
use of simpler data structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
53 if (tile.instance_ != NULL) |
| 0 | 54 { |
|
55
b6432a00b103
use of simpler data structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
55 LOG(ERROR) << "Tile with location (" << tileX << "," |
|
b6432a00b103
use of simpler data structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
56 << tileY << ") is indexed twice in level of size " |
| 0 | 57 << totalWidth_ << "x" << totalHeight_; |
| 58 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadFileFormat); | |
| 59 } | |
| 60 | |
|
55
b6432a00b103
use of simpler data structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
61 tile.instance_ = &instance; |
|
b6432a00b103
use of simpler data structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
62 tile.frame_ = frame; |
| 0 | 63 } |
| 64 | |
| 65 | |
| 66 bool DicomPyramidLevel::LookupTile(TileContent& tile, | |
| 67 unsigned int tileX, | |
| 68 unsigned int tileY) const | |
| 69 { | |
|
55
b6432a00b103
use of simpler data structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
70 const TileContent& tmp = const_cast<DicomPyramidLevel&>(*this).GetTileContent(tileX, tileY); |
|
b6432a00b103
use of simpler data structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
71 |
|
b6432a00b103
use of simpler data structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
72 if (tmp.instance_ == NULL) |
| 0 | 73 { |
| 74 return false; | |
| 75 } | |
| 76 else | |
| 77 { | |
|
55
b6432a00b103
use of simpler data structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
78 tile = tmp; |
| 0 | 79 return true; |
| 80 } | |
| 81 } | |
| 82 | |
| 83 | |
|
56
83cd735c885d
speedup the loading of DICOM sources
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
55
diff
changeset
|
84 DicomPyramidLevel::DicomPyramidLevel(DicomPyramidInstance& instance) : |
| 0 | 85 totalWidth_(instance.GetTotalWidth()), |
| 86 totalHeight_(instance.GetTotalHeight()), | |
| 87 tileWidth_(instance.GetTileWidth()), | |
| 88 tileHeight_(instance.GetTileHeight()) | |
| 89 { | |
| 90 if (totalWidth_ == 0 || | |
| 91 totalHeight_ == 0) | |
| 92 { | |
| 93 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange); | |
| 94 } | |
|
55
b6432a00b103
use of simpler data structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
95 |
|
b6432a00b103
use of simpler data structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
96 countTilesX_ = CeilingDivision(totalWidth_, tileWidth_); |
|
b6432a00b103
use of simpler data structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
97 countTilesY_ = CeilingDivision(totalHeight_, tileHeight_); |
|
b6432a00b103
use of simpler data structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
98 tiles_.resize(countTilesX_ * countTilesY_); |
| 0 | 99 |
| 100 AddInstance(instance); | |
| 101 } | |
| 102 | |
| 103 | |
|
56
83cd735c885d
speedup the loading of DICOM sources
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
55
diff
changeset
|
104 void DicomPyramidLevel::AddInstance(DicomPyramidInstance& instance) |
| 0 | 105 { |
| 106 if (instance.GetTotalWidth() != totalWidth_ || | |
| 107 instance.GetTotalHeight() != totalHeight_ || | |
| 108 instance.GetTileWidth() != tileWidth_ || | |
| 109 instance.GetTileHeight() != tileHeight_) | |
| 110 { | |
| 111 throw Orthanc::OrthancException(Orthanc::ErrorCode_IncompatibleImageSize); | |
| 112 } | |
| 113 | |
| 114 for (size_t frame = 0; frame < instance.GetFrameCount(); frame++) | |
| 115 { | |
| 116 RegisterFrame(instance, frame); | |
| 117 } | |
| 118 } | |
| 119 | |
| 120 | |
|
57
91fc9583b2de
big refactoring to support sparse tiling
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
56
diff
changeset
|
121 bool DicomPyramidLevel::DownloadRawTile(std::string& raw /* out */, |
| 0 | 122 Orthanc::PixelFormat& format /* out */, |
|
57
91fc9583b2de
big refactoring to support sparse tiling
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
56
diff
changeset
|
123 ImageCompression& compression /* out */, |
|
196
b0bd22077cd8
sharing code with orthanc-stone
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
192
diff
changeset
|
124 OrthancStone::IOrthancConnection& orthanc, |
| 0 | 125 unsigned int tileX, |
| 126 unsigned int tileY) const | |
| 127 { | |
| 128 TileContent tile; | |
| 129 if (LookupTile(tile, tileX, tileY)) | |
| 130 { | |
|
55
b6432a00b103
use of simpler data structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
131 assert(tile.instance_ != NULL); |
|
56
83cd735c885d
speedup the loading of DICOM sources
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
55
diff
changeset
|
132 DicomPyramidInstance& instance = *tile.instance_; |
| 0 | 133 |
| 134 std::string uri = ("/instances/" + instance.GetInstanceId() + | |
|
55
b6432a00b103
use of simpler data structures
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
16
diff
changeset
|
135 "/frames/" + boost::lexical_cast<std::string>(tile.frame_) + "/raw"); |
| 0 | 136 |
| 137 orthanc.RestApiGet(raw, uri); | |
| 138 | |
|
56
83cd735c885d
speedup the loading of DICOM sources
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
55
diff
changeset
|
139 compression = instance.GetImageCompression(orthanc); |
| 0 | 140 format = instance.GetPixelFormat(); |
| 141 | |
| 142 return true; | |
| 143 } | |
| 144 else | |
| 145 { | |
| 146 return false; | |
| 147 } | |
| 148 } | |
| 149 } |