Mercurial > hg > orthanc-stone
annotate OrthancStone/Sources/Volumes/DicomVolumeImageReslicer.cpp @ 1771:f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
| author | Sebastien Jodogne <s.jodogne@gmail.com> |
|---|---|
| date | Wed, 12 May 2021 15:09:32 +0200 |
| parents | 9ac2a65d4172 |
| children | fca942f4b4a7 |
| rev | line source |
|---|---|
| 815 | 1 /** |
| 2 * Stone of Orthanc | |
| 3 * Copyright (C) 2012-2016 Sebastien Jodogne, Medical Physics | |
| 4 * Department, University Hospital of Liege, Belgium | |
|
1739
9ac2a65d4172
upgrade to year 2021
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1598
diff
changeset
|
5 * Copyright (C) 2017-2021 Osimis S.A., Belgium |
| 815 | 6 * |
| 7 * This program is free software: you can redistribute it and/or | |
|
1598
8563ea5d8ae4
relicensing some files, cf. osimis bm26 and chu agreement on 2020-05-20
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1596
diff
changeset
|
8 * modify it under the terms of the GNU Lesser General Public License |
| 815 | 9 * as published by the Free Software Foundation, either version 3 of |
| 10 * the License, or (at your option) any later version. | |
| 11 * | |
| 12 * This program is distributed in the hope that it will be useful, but | |
| 13 * WITHOUT ANY WARRANTY; without even the implied warranty of | |
|
1598
8563ea5d8ae4
relicensing some files, cf. osimis bm26 and chu agreement on 2020-05-20
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1596
diff
changeset
|
14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU |
|
8563ea5d8ae4
relicensing some files, cf. osimis bm26 and chu agreement on 2020-05-20
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1596
diff
changeset
|
15 * Lesser General Public License for more details. |
|
1596
4fb8fdf03314
removed annoying whitespace
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1571
diff
changeset
|
16 * |
|
1598
8563ea5d8ae4
relicensing some files, cf. osimis bm26 and chu agreement on 2020-05-20
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1596
diff
changeset
|
17 * You should have received a copy of the GNU Lesser General Public |
|
8563ea5d8ae4
relicensing some files, cf. osimis bm26 and chu agreement on 2020-05-20
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1596
diff
changeset
|
18 * License along with this program. If not, see |
|
8563ea5d8ae4
relicensing some files, cf. osimis bm26 and chu agreement on 2020-05-20
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1596
diff
changeset
|
19 * <http://www.gnu.org/licenses/>. |
| 815 | 20 **/ |
| 21 | |
| 22 | |
| 23 #include "DicomVolumeImageReslicer.h" | |
| 24 | |
|
1455
30deba7bc8e2
simplifying include_directories
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1298
diff
changeset
|
25 #include <OrthancException.h> |
| 815 | 26 |
| 27 namespace OrthancStone | |
| 28 { | |
| 29 class DicomVolumeImageReslicer::Slice : public IVolumeSlicer::IExtractedSlice | |
| 30 { | |
| 31 private: | |
| 32 DicomVolumeImageReslicer& that_; | |
| 33 CoordinateSystem3D cuttingPlane_; | |
| 34 | |
| 35 public: | |
| 36 Slice(DicomVolumeImageReslicer& that, | |
| 37 const CoordinateSystem3D& cuttingPlane) : | |
| 38 that_(that), | |
| 39 cuttingPlane_(cuttingPlane) | |
| 40 { | |
| 41 } | |
| 42 | |
| 1571 | 43 virtual bool IsValid() ORTHANC_OVERRIDE |
| 815 | 44 { |
| 45 return true; | |
| 46 } | |
| 47 | |
| 1571 | 48 virtual uint64_t GetRevision() ORTHANC_OVERRIDE |
| 815 | 49 { |
| 50 return that_.volume_->GetRevision(); | |
| 51 } | |
| 52 | |
| 53 virtual ISceneLayer* CreateSceneLayer(const ILayerStyleConfigurator* configurator, | |
| 1571 | 54 const CoordinateSystem3D& cuttingPlane) ORTHANC_OVERRIDE |
| 815 | 55 { |
| 56 VolumeReslicer& reslicer = that_.reslicer_; | |
| 57 | |
| 58 if (configurator == NULL) | |
| 59 { | |
| 60 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError, | |
| 61 "Must provide a layer style configurator"); | |
| 62 } | |
|
1771
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
63 |
| 815 | 64 reslicer.SetOutputFormat(that_.volume_->GetPixelData().GetFormat()); |
| 65 reslicer.Apply(that_.volume_->GetPixelData(), | |
| 66 that_.volume_->GetGeometry(), | |
| 67 cuttingPlane); | |
| 68 | |
| 69 if (reslicer.IsSuccess()) | |
| 70 { | |
|
1771
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
71 std::unique_ptr<TextureBaseSceneLayer> texture |
| 815 | 72 (configurator->CreateTextureFromDicom(reslicer.GetOutputSlice(), |
| 73 that_.volume_->GetDicomParameters())); | |
|
1771
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
74 if (texture.get() == NULL) |
| 815 | 75 { |
| 76 return NULL; | |
| 77 } | |
| 78 | |
|
1771
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
79 const double s = reslicer.GetPixelSpacing(); |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
80 |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
81 #if 1 |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
82 const double x1 = reslicer.GetOutputExtent().GetX1(); |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
83 //const double x2 = reslicer.GetOutputExtent().GetX2(); |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
84 const double y1 = reslicer.GetOutputExtent().GetY1(); |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
85 const double y2 = reslicer.GetOutputExtent().GetY2(); |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
86 |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
87 const Vector p1 = cuttingPlane.MapSliceToWorldCoordinates(x1, y1); |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
88 const Vector p2 = cuttingPlane.MapSliceToWorldCoordinates(x1, y2); |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
89 |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
90 if (1) |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
91 { |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
92 texture->SetCuttingPlaneTransform(cuttingPlane, p1, |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
93 s * cuttingPlane.GetAxisX(), |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
94 s * cuttingPlane.GetAxisY()); |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
95 } |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
96 else |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
97 { |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
98 /** |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
99 * TODO - ONE WAS TO SOMETIMES FLIP the Y axis. Is it also |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
100 * possible for the X axis? |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
101 **/ |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
102 |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
103 texture->SetCuttingPlaneTransform(cuttingPlane, p2, |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
104 s * cuttingPlane.GetAxisX(), |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
105 -s * cuttingPlane.GetAxisY()); |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
106 } |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
107 |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
108 #else |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
109 texture->SetPixelSpacing(s, s); |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
110 texture->SetOrigin(reslicer.GetOutputExtent().GetX1() + 0.5 * s, |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
111 reslicer.GetOutputExtent().GetY1() + 0.5 * s); |
|
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
112 //texture->SetFlipY(true); |
| 815 | 113 |
| 114 // TODO - Angle!! | |
|
1771
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
115 #endif |
| 815 | 116 |
|
1771
f302bbddf94d
sync, trying to fix DicomVolumeImageReslicer
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
1739
diff
changeset
|
117 return texture.release(); |
| 815 | 118 } |
| 119 else | |
| 120 { | |
| 121 return NULL; | |
| 122 } | |
| 123 } | |
| 124 }; | |
| 125 | |
| 126 | |
| 127 DicomVolumeImageReslicer::DicomVolumeImageReslicer(const boost::shared_ptr<DicomVolumeImage>& volume) : | |
| 128 volume_(volume) | |
| 129 { | |
| 130 if (volume.get() == NULL) | |
| 131 { | |
| 132 throw Orthanc::OrthancException(Orthanc::ErrorCode_NullPointer); | |
| 133 } | |
| 134 } | |
| 135 | |
| 136 | |
| 137 IVolumeSlicer::IExtractedSlice* DicomVolumeImageReslicer::ExtractSlice(const CoordinateSystem3D& cuttingPlane) | |
| 138 { | |
| 139 if (volume_->HasGeometry()) | |
| 140 { | |
| 141 return new Slice(*this, cuttingPlane); | |
| 142 } | |
| 143 else | |
| 144 { | |
| 145 return new IVolumeSlicer::InvalidSlice; | |
| 146 } | |
| 147 } | |
| 148 } |