Add read_nwb_sorting_analyzer: build a SortingAnalyzer from an NWB Units table

src/spikeinterface/core/basesorting.py

@@ -471,6 +471,23 @@ def get_last_spike_frame(self, segment_index: int | None = None) -> int:
             return 0
         return int(np.max(spikes_in_segment["sample_index"]))
 
+    def get_last_spike_time(self, segment_index: int | None = None) -> float:
+        """Get the time of the last spike in a segment across all units.
+
+        Parameters
+        ----------
+        segment_index : int or None, default: None
+            The segment index (required for multi-segment)
+
+        Returns
+        -------
+        float
+            The time of the last spike in seconds, or 0.0 if no spikes exist.
+        """
+        segment_index = self._check_segment_index(segment_index)
+        last_spike_frame = self.get_last_spike_frame(segment_index=segment_index)
+        return self.sample_index_to_time(last_spike_frame, segment_index=segment_index)
+
     def get_times(
         self,
         segment_index: int | None = None,

src/spikeinterface/extractors/__init__.py

@@ -5,6 +5,7 @@
 
 from .neoextractors import get_neo_num_blocks, get_neo_streams
 from .phykilosortextractors import read_kilosort_as_analyzer
+from .nwbextractors import read_nwb_sorting_analyzer
 
 from warnings import warn