--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/Framework/Volumes/DicomVolumeImageMPRSlicer.cpp	Tue May 28 21:16:39 2019 +0200
@@ -0,0 +1,120 @@
+/**
+ * Stone of Orthanc
+ * Copyright (C) 2012-2016 Sebastien Jodogne, Medical Physics
+ * Department, University Hospital of Liege, Belgium
+ * Copyright (C) 2017-2019 Osimis S.A., Belgium
+ *
+ * This program is free software: you can redistribute it and/or
+ * modify it under the terms of the GNU Affero General Public License
+ * as published by the Free Software Foundation, either version 3 of
+ * the License, or (at your option) any later version.
+ *
+ * This program is distributed in the hope that it will be useful, but
+ * WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
+ * Affero General Public License for more details.
+ * 
+ * You should have received a copy of the GNU Affero General Public License
+ * along with this program. If not, see <http://www.gnu.org/licenses/>.
+ **/
+
+
+#include "DicomVolumeImageMPRSlicer.h"
+
+#include <Core/OrthancException.h>
+
+namespace OrthancStone
+{
+  void DicomVolumeImageMPRSlicer::Slice::CheckValid() const
+  {
+    if (!valid_)
+    {
+      throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls);
+    }
+  }
+
+
+  DicomVolumeImageMPRSlicer::Slice::Slice(const DicomVolumeImage& volume,
+                                          const CoordinateSystem3D& cuttingPlane) :
+    volume_(volume)
+  {
+    valid_ = (volume_.HasDicomParameters() &&
+              volume_.GetGeometry().DetectSlice(projection_, sliceIndex_, cuttingPlane));
+  }
+
+
+  VolumeProjection DicomVolumeImageMPRSlicer::Slice::GetProjection() const
+  {
+    CheckValid();
+    return projection_;
+  }
+
+
+  unsigned int DicomVolumeImageMPRSlicer::Slice::GetSliceIndex() const
+  {
+    CheckValid();
+    return sliceIndex_;
+  }
+
+  uint64_t DicomVolumeImageMPRSlicer::Slice::GetRevision()
+  {
+    CheckValid();
+    return GetRevisionInternal(projection_, sliceIndex_);
+  }
+
+
+  ISceneLayer* DicomVolumeImageMPRSlicer::Slice::CreateSceneLayer(const ILayerStyleConfigurator* configurator,
+                                                                  const CoordinateSystem3D& cuttingPlane)
+  {
+    CheckValid();
+
+    if (configurator == NULL)
+    {
+      throw Orthanc::OrthancException(Orthanc::ErrorCode_NullPointer,
+                                      "A style configurator is mandatory for textures");
+    }
+
+    std::auto_ptr<TextureBaseSceneLayer> texture;
+        
+    {
+      const DicomInstanceParameters& parameters = volume_.GetDicomParameters();
+      ImageBuffer3D::SliceReader reader(volume_.GetPixelData(), projection_, sliceIndex_);
+      texture.reset(dynamic_cast<TextureBaseSceneLayer*>
+                    (configurator->CreateTextureFromDicom(reader.GetAccessor(), parameters)));
+    }
+
+    const CoordinateSystem3D& system = volume_.GetGeometry().GetProjectionGeometry(projection_);
+      
+    double x0, y0, x1, y1;
+    cuttingPlane.ProjectPoint(x0, y0, system.GetOrigin());
+    cuttingPlane.ProjectPoint(x1, y1, system.GetOrigin() + system.GetAxisX());
+    texture->SetOrigin(x0, y0);
+
+    double dx = x1 - x0;
+    double dy = y1 - y0;
+    if (!LinearAlgebra::IsCloseToZero(dx) ||
+        !LinearAlgebra::IsCloseToZero(dy))
+    {
+      texture->SetAngle(atan2(dy, dx));
+    }
+        
+    Vector tmp = volume_.GetGeometry().GetVoxelDimensions(projection_);
+    texture->SetPixelSpacing(tmp[0], tmp[1]);
+
+    return texture.release();
+  }
+
+
+  IVolumeSlicer::IExtractedSlice* 
+  DicomVolumeImageMPRSlicer::ExtractSlice(const CoordinateSystem3D& cuttingPlane)
+  {
+    if (volume_->HasGeometry())
+    {
+      return new Slice(*volume_, cuttingPlane);
+    }
+    else
+    {
+      return new IVolumeSlicer::InvalidSlice;
+    }
+  }
+}