Merge pull request #2170 from effigies/fix/repetition_time_mismatch

fix: Account for NIfTI time units when checking RepetitionTime

Author: rwblair

src/schema/meta/expression_tests.yaml

@@ -74,6 +74,14 @@
 # General expressions
 - expression: 1 + 2
   result: 3
+- expression: 1 - 2
+  result: -1
+- expression: 3 * 4
+  result: 12
+- expression: 3 / 2
+  result: 1.5
+- expression: 3 % 2
+  result: 1
 - expression: '"cat" + "dog"'
   result: 'catdog'
 - expression: match('string', '.*')

src/schema/rules/checks/func.yaml

@@ -40,9 +40,17 @@ RepetitionTimeMismatch:
     - type(sidecar.RepetitionTime) != "null"
     - type(nifti_header) != "null"
   checks:
-    # Implement millisecond rounding via AND
-    - sidecar.RepetitionTime - nifti_header.pixdim[4] < 0.001
-    - sidecar.RepetitionTime - nifti_header.pixdim[4] > -0.001
+    # Use index(...) % 3 to treat unknown units like seconds
+    - |
+      nifti_header.pixdim[4]
+        * 10 ** (-3 * (index(["sec", "msec", "usec", "unknown"], nifti_header.xyzt_units.t) % 3))
+      - sidecar.RepetitionTime
+      > -0.001
+    - |
+      nifti_header.pixdim[4]
+        * 10 ** (-3 * (index(["sec", "msec", "usec", "unknown"], nifti_header.xyzt_units.t) % 3))
+      - sidecar.RepetitionTime
+      < 0.001
 
 # 54
 BoldNot4d:

tools/schemacode/src/bidsschematools/expressions.py

@@ -62,7 +62,7 @@ def parse(expression_string: str) -> "ASTNode":
 notOp = Literal("!")
 compOp = one_of(("==", "!=", "<", "<=", ">", ">=", "in"))
 addOp = one_of(("+", "-"))
-mulOp = one_of(("*", "/"))
+mulOp = one_of(("*", "/", "%"))
 expOp = Literal("**")
 
 lpar, rpar = Suppress("("), Suppress(")")