Add numpy 1.26.x support and Python 3.12/3.13 compatibility

.github/workflows/github-actions-ci.yml

@@ -26,7 +26,7 @@ jobs:
     name: Lint
     runs-on: "ubuntu-latest"
     steps:
-      - uses: actions/checkout@v2
+      - uses: actions/checkout@v4
       - name: flake8 linting
         run: |
           pip install flake8
@@ -43,15 +43,24 @@ jobs:
     runs-on: ${{ matrix.os }}
     strategy:
       matrix:
-        os: ["macos-latest", "windows-latest", "ubuntu-latest"]
-        python-version: ["3.8", "3.9", "3.10", "3.11"]
+        # macos-15-intel for x86_64, macos-latest for ARM64
+        os: ["macos-15-intel", "macos-latest", "windows-latest", "ubuntu-latest"]
+        python-version: ["3.10", "3.11", "3.12", "3.13"]
+        exclude:
+          # Exclude Windows Python 3.13 because numpy<2 lacks support: conda-forge
+          # has no numpy 1.x builds, so conda falls back to an experimental
+          # MINGW-W64 build that segfaults on import; remove numpy<2 to enable it.
+          - os: windows-latest
+            python-version: "3.13"
+      fail-fast: false
+    defaults:
       fail-fast: false
     defaults:
       run:
         shell: bash -l {0}
     steps:
-      - uses: actions/checkout@v2
-      - uses: conda-incubator/setup-miniconda@v2
+      - uses: actions/checkout@v4
+      - uses: conda-incubator/setup-miniconda@v3
         with:
           auto-update-conda: true
           python-version: ${{ matrix.python-version }}
@@ -71,9 +80,9 @@ jobs:
     runs-on: ["self-hosted"]
     strategy:
       matrix:
-        image: ["allensdk_local_py38:latest"]
+        image: ["allensdk_local_py310:latest"]
     steps:
-      - uses: actions/checkout@v2
+      - uses: actions/checkout@v4
       - name: run test in docker
         run: |
           docker run \

.github/workflows/nightly.yml

@@ -10,10 +10,10 @@ jobs:
     runs-on: ["self-hosted"]
     strategy:
       matrix:
-        image: ["allensdk_local_py38:latest"]
+        image: ["allensdk_local_py310:latest"]
         branch: ["master", "rc/**"]
     steps:
-      - uses: actions/checkout@v2
+      - uses: actions/checkout@v4
         with:
           ref: ${{ matrix.branch }}
       - name: run test in docker

.github/workflows/notebook_runner.yml

@@ -9,11 +9,11 @@ jobs:
     strategy:
       matrix:
         os: ["ubuntu-latest"]
-        python-version: ["3.8"]
+        python-version: ["3.10"]
       fail-fast: false
     runs-on: ${{ matrix.os }}
     steps:
-      - uses: actions/checkout@v3
+      - uses: actions/checkout@v4
       - name: Set up Python
         uses: actions/setup-python@v4
         with:

allensdk/brain_observatory/behavior/data_objects/running_speed/running_acquisition.py

@@ -129,7 +129,7 @@ def from_nwb(
         cls,
         nwbfile: NWBFile
     ) -> "RunningAcquisition":
-        running_module = nwbfile.modules['running']
+        running_module = nwbfile.processing['running']
         dx_interface = running_module.get_data_interface('dx')
 
         dx = dx_interface.data

allensdk/brain_observatory/behavior/data_objects/running_speed/running_speed.py

@@ -188,7 +188,7 @@ def from_nwb(
         nwbfile: NWBFile,
         filtered=True
     ) -> "RunningSpeed":
-        running_module = nwbfile.modules['running']
+        running_module = nwbfile.processing['running']
         interface_name = 'speed' if filtered else 'speed_unfiltered'
         running_interface = running_module.get_data_interface(interface_name)
 

allensdk/brain_observatory/behavior/data_objects/stimuli/templates.py

@@ -1,4 +1,5 @@
 import os
+import warnings
 import numpy as np
 from pathlib import Path
 from typing import Optional, List, Dict
@@ -184,17 +185,44 @@ def to_nwb(self, nwbfile: NWBFile,
 
             nwbfile.add_stimulus_template(visual_stimulus_image_series)
 
+        # DEPRECATED (11/2025): _add_image_index_to_nwb is no longer called.
+        # The IndexSeries using indexed_timeseries is deprecated in pynwb 2.x
+        # and would require significant refactoring to use indexed_images instead.
+        # Since this NWB writing code will not be used for new data generation,
+        # we skip this step rather than refactor. The stimulus template data
+        # is still written above; only the IndexSeries linking is skipped.
+        # See: https://pynwb.readthedocs.io/en/stable/pynwb.image.html#pynwb.image.IndexSeries
         if 'image_index' in stimulus_presentations.value \
                 and self._image_template_key is not None:
-            nwbfile = self._add_image_index_to_nwb(
-                nwbfile=nwbfile, presentations=stimulus_presentations)
+            warnings.warn(
+                "As of 11/2025, Templates.to_nwb() no longer adds the image "
+                "index (IndexSeries) to the NWB file. The IndexSeries "
+                "indexed_timeseries field is deprecated in pynwb 2.x. "
+                "The stimulus template data is still written.",
+                UserWarning,
+                stacklevel=2
+            )
 
         return nwbfile
 
     def _add_image_index_to_nwb(
             self, nwbfile: NWBFile, presentations: Presentations):
         """Adds the image index and start_time for all stimulus templates
-        to NWB"""
+        to NWB
+
+        .. deprecated:: 2.16.3
+            This method is deprecated as of 11/2025. The IndexSeries
+            indexed_timeseries field is deprecated in pynwb 2.x. This method
+            is no longer called from to_nwb() and will be removed in a future
+            release.
+        """
+        warnings.warn(
+            "_add_image_index_to_nwb is deprecated and will be removed in a "
+            "future release. The IndexSeries indexed_timeseries field is "
+            "deprecated in pynwb 2.x.",
+            DeprecationWarning,
+            stacklevel=2
+        )
         stimulus_templates = self.value[self._image_template_key]
         presentations = presentations.value
 
@@ -208,7 +236,7 @@ def _add_image_index_to_nwb(
         image_index = IndexSeries(
             name=nwb_template.name,
             data=stimulus_index['image_index'].values,
-            unit='None',
+            unit='N/A',
             indexed_timeseries=nwb_template,
             timestamps=stimulus_index['start_time'].values)
         nwbfile.add_stimulus(image_index)

allensdk/brain_observatory/ecephys/ecephys_project_api/http_engine.py

@@ -93,10 +93,10 @@ def write_bytes(path: str, stream: Iterable[bytes]):
 class AsyncHttpEngine(HttpEngine):
 
     def __init__(
-        self, 
-        scheme: str, 
-        host: str, 
-        session: Optional[aiohttp.ClientSession] = None, 
+        self,
+        scheme: str,
+        host: str,
+        session: Optional[aiohttp.ClientSession] = None,
         **kwargs
     ):
         """ Simple tool for making asynchronous streaming http requests.
@@ -105,11 +105,11 @@ def __init__(
         ----------
         scheme :
             e.g "http" or "https"
-        host : 
+        host :
             will be used as the base for request urls
-        session : 
-            If provided, this preconstructed session will be used rather than 
-            a new one. Keep in mind that AsyncHttpEngine closes its session 
+        session :
+            If provided, this preconstructed session will be used rather than
+            a new one. Keep in mind that AsyncHttpEngine closes its session
             when it is garbage collected!
         **kwargs :
             Will be passed to parent.
@@ -119,14 +119,25 @@ def __init__(
         super(AsyncHttpEngine, self).__init__(scheme, host, **kwargs)
 
         if session:
-            self.session = session
+            self._session = session
+            self._owns_session = False
             warnings.warn(
-                "Recieved preconstructed session, ignoring timeout parameter."
+                "Received preconstructed session, ignoring timeout parameter."
             )
         else:
-            self.session = aiohttp.ClientSession(
+            # Defer session creation until actually needed in an async context
+            # (aiohttp 3.9+ requires ClientSession to be created within an event loop)
+            self._session = None
+            self._owns_session = True
+
+    @property
+    def session(self) -> aiohttp.ClientSession:
+        """Lazily create the aiohttp session when first accessed."""
+        if self._session is None:
+            self._session = aiohttp.ClientSession(
                 timeout=aiohttp.client.ClientTimeout(self.timeout)
             )
+        return self._session
 
     async def _stream_coroutine(
         self, 
@@ -169,10 +180,10 @@ def stream(
         return functools.partial(self._stream_coroutine, route)
 
     def __del__(self):
-        if hasattr(self, "session"):
+        if hasattr(self, "_session") and self._session is not None:
             nest_asyncio.apply()
             loop = asyncio.get_event_loop()
-            loop.run_until_complete(self.session.close())
+            loop.run_until_complete(self._session.close())
 
     @staticmethod
     def write_bytes(

allensdk/brain_observatory/ecephys/ecephys_project_cache.py

@@ -737,9 +737,12 @@ def get_grouped_uniques(this, other, foreign_key, field_key, unique_key, inplace
     if not inplace:
         this = this.copy()
 
-    uniques = other.groupby(foreign_key)\
-        .apply(lambda grp: pd.DataFrame(grp)[field_key].unique())
-    this[unique_key] = 0
+    # Select only the field_key column before apply to avoid FutureWarning about
+    # grouping columns in pandas 2.2+
+    uniques = other.groupby(foreign_key)[field_key]\
+        .apply(lambda x: x.unique())
+    # Use object dtype to allow storing arrays of strings
+    this[unique_key] = None
     this.loc[uniques.index.values, unique_key] = uniques.values
 
     if not inplace:

allensdk/brain_observatory/gaze_mapping/_gaze_mapper.py

@@ -289,7 +289,9 @@ def pupil_position_on_monitor_in_degrees(self,
 
         mag = np.linalg.norm(self.monitor.position)
         meridian = np.degrees(np.arctan(x / mag))
-        elevation = np.degrees(np.arctan(y / np.linalg.norm([x, mag], axis=0)))
+        # Use np.vstack to create a homogeneous 2D array (numpy 1.24+ compatibility)
+        elevation = np.degrees(np.arctan(y / np.linalg.norm(
+            np.vstack([x, np.full_like(x, mag, dtype=float)]), axis=0)))
 
         angles = np.vstack([meridian, elevation]).T
 

allensdk/brain_observatory/nwb/nwb_api.py

@@ -59,9 +59,9 @@ def get_running_speed(self, lowpass=True) -> RunningSpeed:
         """
 
         interface_name = 'speed' if lowpass else 'speed_unfiltered'
-        values = self.nwbfile.modules['running'].get_data_interface(
+        values = self.nwbfile.processing['running'].get_data_interface(
             interface_name).data[:]
-        timestamps = self.nwbfile.modules['running'].get_data_interface(
+        timestamps = self.nwbfile.processing['running'].get_data_interface(
             interface_name).timestamps[:]
 
         return RunningSpeed(
@@ -87,7 +87,7 @@ def get_image(self, name, module, image_api=None) -> sitk.Image:
         if image_api is None:
             image_api = ImageApi
 
-        nwb_img = self.nwbfile.modules[module].get_data_interface(
+        nwb_img = self.nwbfile.processing[module].get_data_interface(
             'images')[name]
         data = nwb_img.data
         resolution = nwb_img.resolution  # px/cm

allensdk/brain_observatory/receptive_field_analysis/chisquarerf.py

@@ -313,8 +313,16 @@ def interpolate_RF(rf_map, deg_per_pnt):
         1,
     )
 
-    interpolated = si.interp2d(x_coor, y_coor, rf_map)
-    interpolated = interpolated(x_interpolated, y_interpolated)
+    # interp2d was removed in scipy 1.14, use RectBivariateSpline instead
+    try:
+        interpolator = si.interp2d(x_coor, y_coor, rf_map)
+        interpolated = interpolator(x_interpolated, y_interpolated)
+    except NotImplementedError:
+        # RectBivariateSpline uses (row, col) ordering vs interp2d's (x, y),
+        # so arguments are swapped: (y_coor, x_coor) and (y_interpolated, x_interpolated)
+        # Use kx=ky=1 (linear) to match interp2d default behavior
+        interpolator = si.RectBivariateSpline(y_coor, x_coor, rf_map, kx=1, ky=1)
+        interpolated = interpolator(y_interpolated, x_interpolated)
 
     return interpolated
 

allensdk/config/app/application_config.py

@@ -351,8 +351,9 @@ def apply_configuration_from_file(self, config_file_path):
             cfg_string = self.from_json_file(config_file_path)
             try:
                 config.readfp(io.BytesIO(cfg_string))
-            except (NameError, TypeError):
-                config.read_string(cfg_string)  # Python 3
+            except (NameError, TypeError, AttributeError):
+                # readfp was removed in Python 3.12
+                config.read_string(cfg_string)
         else:
             config.read(config_file_path)