Solving issue #13122: Aiming to achieve better sensitivity of head m…

…ovement detection for chpi extraction (#13123)

Co-authored-by: pre-commit-ci[bot] <66853113+pre-commit-ci[bot]@users.noreply.github.com>

doc/changes/devel/13123.bugfix.rst

@@ -0,0 +1 @@
+Fixed CTF cHPI coil location detection in :func:`mne.chpi.extract_chpi_locs_ctf` by :newcontrib:`Shresth Keshari`.

doc/changes/names.inc

@@ -277,6 +277,7 @@
 .. _Senwen Deng: https://snwn.de
 .. _Seyed Yahya Shirazi: https://neuromechanist.github.io
 .. _Sheraz Khan: https://github.com/SherazKhan
+.. _Shresth Keshari: https://github.com/shresth-keshari
 .. _Shristi Baral: https://github.com/shristibaral
 .. _Silvia Cotroneo: https://github.com/sfc-neuro
 .. _Simeon Wong: https://github.com/dtxe

mne/chpi.py

@@ -234,7 +234,7 @@ def extract_chpi_locs_ctf(raw, verbose=None):
 
     # find indices where chpi locations change
     indices = [0]
-    indices.extend(np.where(np.all(np.diff(chpi_data, axis=1), axis=0))[0] + 1)
+    indices.extend(np.where(np.any(np.diff(chpi_data, axis=1), axis=0))[0] + 1)
     # data in channels are in ctf device coordinates (cm)
     rrs = chpi_data[:, indices].T.reshape(len(indices), 3, 3)  # m
     # map to mne device coords