diff --git a/src/cai_lab_to_nwb/zaki_2024/utils/__init__.py b/src/cai_lab_to_nwb/zaki_2024/utils/__init__.py
index a9b7d56..38262bf 100644
--- a/src/cai_lab_to_nwb/zaki_2024/utils/__init__.py
+++ b/src/cai_lab_to_nwb/zaki_2024/utils/__init__.py
@@ -1,2 +1 @@
-from .motion_correction import add_motion_correction, load_motion_correction_data
 from .cell_registration import add_cell_registration, get_global_ids_from_csv
diff --git a/src/cai_lab_to_nwb/zaki_2024/utils/motion_correction.py b/src/cai_lab_to_nwb/zaki_2024/utils/motion_correction.py
deleted file mode 100644
index b09af0f..0000000
--- a/src/cai_lab_to_nwb/zaki_2024/utils/motion_correction.py
+++ /dev/null
@@ -1,80 +0,0 @@
-from typing import List
-
-import numpy as np
-from neuroconv.tools import get_module
-from pydantic.types import PathType
-from pynwb import NWBFile, TimeSeries
-from roiextractors.extraction_tools import DtypeType
-
-
-def load_motion_correction_data(folder_path: PathType) -> np.ndarray:
-    """Read the xy shifts in the 'motion' field from the zarr object.
-
-    Parameters
-    ----------
-    folder_path: PathType
-        Path to minian output.
-
-    Returns
-    -------
-    motion_correction: numpy.ndarray
-        The first column is the x shifts. The second column is the y shifts.
-    """
-    import zarr
-    import warnings
-
-    if "/motion.zarr" not in zarr.open(folder_path, mode="r"):
-        warnings.warn(f"Group '/motion.zarr' not found in the Zarr store.", UserWarning)
-        return None
-    else:
-        dataset = zarr.open(str(folder_path) + "/motion.zarr", "r")
-        # from zarr field motion.zarr/shift_dim we can verify that the two column refer respectively to
-        # ['height','width'] --> ['y','x']. Following best practice we swap the two columns
-        motion_correction = dataset["motion"][:, [1, 0]]
-        return motion_correction
-
-
-def add_motion_correction(
-    nwbfile: NWBFile,
-    motion_correction_series: np.ndarray,
-    one_photon_series_name: str,
-) -> None:
-    """Add motion correction data to the NWBFile.
-
-    The x, y shifts for the imaging data (identified by 'one_photon_series_name' are added to the NWBFile as a TimeSeries.
-    The series is added to the 'ophys' processing module.
-
-    Parameters
-    ----------
-    nwbfile: NWBFile
-        The NWBFile where the motion correction time series will be added to.
-    motion_correction_series: numpy.ndarray
-        The x, y shifts for the imaging data.
-    one_photon_series_name: str
-        The name of the one photon series in the NWBFile.
-    """
-
-    assert (
-        one_photon_series_name in nwbfile.acquisition
-    ), f"The one photon series '{one_photon_series_name}' does not exist in the NWBFile."
-    name_suffix = one_photon_series_name.replace("OnePhotonSeries", "")
-    motion_correction_time_series_name = "MotionCorrectionSeries" + name_suffix
-    ophys = get_module(nwbfile, "ophys")
-    if motion_correction_time_series_name in ophys.data_interfaces:
-        raise ValueError(
-            f"The motion correction time series '{motion_correction_time_series_name}' already exists in the NWBFile."
-        )
-
-    one_photon_series = nwbfile.acquisition[one_photon_series_name]
-    num_frames = one_photon_series.data.maxshape[0]
-    assert (
-        num_frames == motion_correction_series.shape[0]
-    ), f"The number of frames for motion correction ({motion_correction_series.shape[0]}) does not match the number of frames ({num_frames}) from the {one_photon_series_name} imaging data."
-    xy_translation = TimeSeries(
-        name="MotionCorrectionSeries" + name_suffix,
-        description=f"The x, y shifts for the {one_photon_series_name} imaging data.",
-        data=motion_correction_series,
-        unit="px",
-        timestamps=one_photon_series.timestamps,
-    )
-    ophys.add(xy_translation)