fix / workaround

MagneticResonanceImaging · May 3, 2021 · 70b91c8 · 70b91c8 · tknopp · May 6, 2021
1 parent 35b0aba
commit 70b91c8
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git 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
@@ -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