Merge pull request #1310 from hornauerp/master

Fixed bug in maxwellrawio.py that shuffled channels
NeuralEnsemble · Jul 27, 2023 · 479c86f · 479c86f
2 parents 972901f + d97fe0c
commit 479c86f
Showing 1 changed file with 4 additions and 3 deletions.
diff --git a/neo/rawio/maxwellrawio.py b/neo/rawio/maxwellrawio.py
@@ -183,9 +183,10 @@ def _get_analogsignal_chunk(self, block_index, seg_index, i_start, i_stop,
             if np.array(channel_indexes).size > 1 and np.any(np.diff(channel_indexes) < 0):
                 # get around h5py constraint that it does not allow datasets
                 # to be indexed out of order
-                sorted_channel_indexes = np.sort(channel_indexes)
-                resorted_indexes = np.array(
-                    [list(channel_indexes).index(ch) for ch in sorted_channel_indexes])
+                order_f = np.argsort(channel_indexes) 
+                sorted_channel_indexes = channel_indexes[order_f] 
+                # use argsort again on order_f to obtain resorted_indexes 
+                resorted_indexes = np.argsort(order_f)
 
         try:
             if resorted_indexes is None: