# HG changeset patch
# User Sebastien Jodogne <s.jodogne@gmail.com>
# Date 1558359841 -7200
# Node ID e63c8b9b7b02c969508c60e6a0c051e7774d4516
# Parent  f334b098b243284833b999e393765ce37aa47d04
setting texture geometry

diff -r f334b098b243 -r e63c8b9b7b02 Samples/Sdl/Loader.cpp
--- a/Samples/Sdl/Loader.cpp	Mon May 20 14:56:03 2019 +0200
+++ b/Samples/Sdl/Loader.cpp	Mon May 20 15:44:01 2019 +0200
@@ -1562,7 +1562,7 @@
             new Refactoring::GetOrthancWebViewerJpegCommand);
           tmp->SetHttpHeader("Accept-Encoding", "gzip");
           tmp->SetInstance(instance);
-          tmp->SetQuality((quality == 0 ? 50 : 95));
+          tmp->SetQuality((quality == 0 ? 50 : 90));
           tmp->SetExpectedPixelFormat(slice.GetExpectedPixelFormat());
           command.reset(tmp.release());
         }
@@ -1794,6 +1794,7 @@
   class DicomVolumeMPRSlicer : public IVolumeSlicer
   {
   private:
+    bool                            linearInterpolation_;
     OrthancStone::Scene2D&          scene_;
     int                             layerDepth_;
     IDicomVolumeSource&             source_;
@@ -1806,12 +1807,23 @@
     DicomVolumeMPRSlicer(OrthancStone::Scene2D& scene,
                          int layerDepth,
                          IDicomVolumeSource& source) :
+      linearInterpolation_(false),
       scene_(scene),
       layerDepth_(layerDepth),
       source_(source),
       first_(true)
     {
     }
+
+    void SetLinearInterpolation(bool enabled)
+    {
+      linearInterpolation_ = enabled;
+    }
+
+    bool IsLinearInterpolation() const
+    {
+      return linearInterpolation_;
+    }
     
     virtual void SetViewportPlane(const OrthancStone::CoordinateSystem3D& plane)
     {
@@ -1822,9 +1834,11 @@
         return;
       }
 
+      const OrthancStone::VolumeImageGeometry& geometry = source_.GetVolume().GetImage().GetGeometry();
+
       OrthancStone::VolumeProjection projection;
       unsigned int sliceIndex;
-      if (!source_.GetVolume().GetImage().GetGeometry().DetectSlice(projection, sliceIndex, plane))
+      if (!geometry.DetectSlice(projection, sliceIndex, plane))
       {
         // The cutting plane is neither axial, nor coronal, nor
         // sagittal. Could use "VolumeReslicer" here.
@@ -1874,12 +1888,26 @@
           texture.reset(parameters.CreateTexture(reader.GetAccessor()));
         }
 
-        // TODO - 
-        // void SetOrigin(double x, double y);
-        // void SetPixelSpacing(double sx, double sy);
-        // void SetAngle(double angle);
-        // void SetLinearInterpolation(bool isLinearInterpolation);
+        const OrthancStone::CoordinateSystem3D& system = geometry.GetProjectionGeometry(projection);
+
+        double x0, y0, x1, y1;
+        system.ProjectPoint(x0, y0, system.GetOrigin());
+        system.ProjectPoint(x0, y0, system.GetOrigin() + system.GetAxisX());
+        texture->SetOrigin(x0, y0);
 
+        double dx = x1 - x0;
+        double dy = y1 - y0;
+        if (!OrthancStone::LinearAlgebra::IsCloseToZero(dx) ||
+            !OrthancStone::LinearAlgebra::IsCloseToZero(dy))
+        {
+          texture->SetAngle(atan2(dy, dx));
+        }
+        
+        OrthancStone::Vector tmp;
+        geometry.GetVoxelDimensions(projection);
+        texture->SetPixelSpacing(tmp[0], tmp[1]);
+
+        texture->SetLinearInterpolation(linearInterpolation_);
     
         scene_.SetLayer(layerDepth_, texture.release());    
       }