Add check for valid points after downsampling

spinicist · May 12, 2023 · 7b571bd · 7b571bd
1 parent d3fc6fe
commit 7b571bd
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 1 deletion.
diff --git a/src/cmd/downsamp.cpp b/src/cmd/downsamp.cpp
@@ -25,4 +25,4 @@ int main_downsamp(args::Subparser &parser)
   writer.writeTensor(HD5::Keys::Noncartesian, ks2.dimensions(), ks2.data());
 
   return EXIT_SUCCESS;
-}
+}
diff --git a/src/trajectory.cpp b/src/trajectory.cpp
@@ -104,6 +104,9 @@ auto Trajectory::downsample(float const res, Index const lores, bool const shrin
       }
     }
   }
+  if (minSamp > maxSamp) {
+    Log::Fail("No valid trajectory points remain after downsampling");
+  }
   Index const dsSamples = maxSamp + 1 - minSamp;
   Log::Print(
     "Downsample res {} mm, factor {}, matrix {}, voxel-size {} mm, read-points {}-{}{}",
-Original file line number
+Diff line change
@@ Expand Up / @@ -25,4 +25,4 @@ int main_downsamp(args::Subparser &parser) @@
       writer.writeTensor(HD5::Keys::Noncartesian, ks2.dimensions(), ks2.data());
       return EXIT_SUCCESS;
-    }
+    }