From 70b91c8262d4c34b1e090e200da47853249bcbbe Mon Sep 17 00:00:00 2001
From: Tobias Knopp <tobias@knoppweb.de>
Date: Mon, 3 May 2021 22:52:02 +0200
Subject: [PATCH] fix / workaround

---
 src/Datatypes/RawAcqData.jl | 2 +-
 src/Operators/EncodingOp.jl | 4 ++--
 2 files changed, 3 insertions(+), 3 deletions(-)

diff --git a/src/Datatypes/RawAcqData.jl b/src/Datatypes/RawAcqData.jl
index 856f9c7b..51afd0bd 100644
--- a/src/Datatypes/RawAcqData.jl
+++ b/src/Datatypes/RawAcqData.jl
@@ -157,7 +157,7 @@ function subsampleIndices(f::RawAcquisitionData; slice::Int=1, contrast::Int=1,
   encSt2 = encSteps2(f)
   numEncSamp, numProf, numSl = f.params["encodedSize"]
   for i=1:length(f.profiles)
-    # only consider data for the specified slice, constrast and repetition
+    # only consider data for the specified slice, contrast and repetition
     if f.profiles[i].head.idx.slice+1 != slice || f.profiles[i].head.idx.contrast+1 != contrast || f.profiles[i].head.idx.repetition != 0
       continue
     end
diff --git a/src/Operators/EncodingOp.jl b/src/Operators/EncodingOp.jl
index 0206ac6c..7864e373 100644
--- a/src/Operators/EncodingOp.jl
+++ b/src/Operators/EncodingOp.jl
@@ -230,7 +230,7 @@ function fourierEncodingOp2d(shape::NTuple{2,Int64}, tr::Trajectory, opName::Str
   end
 
   # subsampling
-  if !isempty(subsampleIdx) && (subsampleIdx != collect(1:size(tr,2))) 
+  if !isempty(subsampleIdx) && (subsampleIdx != collect(1:size(tr,2))) && isCartesian(tr)
     S = SamplingOp(subsampleIdx,(tr.numSamplingPerProfile,tr.numProfiles))
     return S ∘ ftOp
   else
@@ -263,7 +263,7 @@ function fourierEncodingOp3d(shape::NTuple{3,Int64}, tr::Trajectory, opName::Str
   end
 
   # subsampling
-  if !isempty(subsampleIdx) && length(subsampleIdx)!=size(tr,2)
+  if !isempty(subsampleIdx) && (subsampleIdx != collect(1:size(tr,2))) && isCartesian(tr)
     S = SamplingOp(subsampleIdx,shape)
     return S ∘ ftOp
   else