junifer 0.0.6.dev317__tar.gz → 0.0.6.dev330__tar.gz

{junifer-0.0.6.dev317 → junifer-0.0.6.dev330}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: junifer
-Version: 0.0.6.dev317
+Version: 0.0.6.dev330
 Summary: JUelich NeuroImaging FEature extractoR
 Author-email: Fede Raimondo <f.raimondo@fz-juelich.de>, Synchon Mandal <s.mandal@fz-juelich.de>
 Maintainer-email: Fede Raimondo <f.raimondo@fz-juelich.de>, Synchon Mandal <s.mandal@fz-juelich.de>

junifer-0.0.6.dev330/docs/changes/newsfragments/409.bugfix

@@ -0,0 +1 @@
+Use correct file output suffices for AFNI-based markers by `Synchon Mandal`_

junifer-0.0.6.dev330/docs/changes/newsfragments/409.enh

@@ -0,0 +1 @@
+Ease asset storage for AFNI-based markers by `Synchon Mandal`_

junifer-0.0.6.dev330/docs/changes/newsfragments/412.change

@@ -0,0 +1 @@
+Introduce new CLI parameter ``--verbose-datalad`` to control ``datalad``'s logging handler by `Fede Raimondo`_

junifer-0.0.6.dev330/docs/changes/newsfragments/412.feature

@@ -0,0 +1 @@
+Add support for coloured :class:`logging.Formatter` to use with ``junifer``'s :class:`logging.StreamHandler` based on environment by `Fede Raimondo`_

{junifer-0.0.6.dev317 → junifer-0.0.6.dev330}/junifer/_version.py

@@ -12,5 +12,5 @@ __version__: str
 __version_tuple__: VERSION_TUPLE
 version_tuple: VERSION_TUPLE
 
-__version__ = version = '0.0.6.dev317'
-__version_tuple__ = version_tuple = (0, 0, 6, 'dev317')
+__version__ = version = '0.0.6.dev330'
+__version_tuple__ = version_tuple = (0, 0, 6, 'dev330')

{junifer-0.0.6.dev317 → junifer-0.0.6.dev330}/junifer/api/queue_context/gnu_parallel_local_adapter.py

@@ -37,6 +37,9 @@ class GnuParallelLocalAdapter(QueueContextAdapter):
         virtual environment of any kind (default None).
     verbose : str, optional
         The level of verbosity (default "info").
+    verbose_datalad : str or None, optional
+        The level of verbosity for datalad. If None, will be the same
+        as ``verbose`` (default None).
     submit : bool, optional
         Whether to submit the jobs (default False).
 
@@ -65,6 +68,7 @@ class GnuParallelLocalAdapter(QueueContextAdapter):
         pre_collect: Optional[str] = None,
         env: Optional[dict[str, str]] = None,
         verbose: str = "info",
+        verbose_datalad: Optional[str] = None,
         submit: bool = False,
     ) -> None:
         """Initialize the class."""
@@ -76,6 +80,7 @@ class GnuParallelLocalAdapter(QueueContextAdapter):
         self._pre_collect = pre_collect
         self._check_env(env)
         self._verbose = verbose
+        self._verbose_datalad = verbose_datalad
         self._submit = submit
 
         self._log_dir = self._job_dir / "logs"
@@ -155,6 +160,11 @@ class GnuParallelLocalAdapter(QueueContextAdapter):
 
     def run(self) -> str:
         """Return run commands."""
+        verbose_args = f"--verbose {self._verbose}"
+        if self._verbose_datalad:
+            verbose_args = (
+                f"{verbose_args} --verbose-datalad {self._verbose_datalad}"
+            )
         return (
             f"#!/usr/bin/env {self._shell}\n\n"
             "# This script is auto-generated by junifer.\n\n"
@@ -169,7 +179,7 @@ class GnuParallelLocalAdapter(QueueContextAdapter):
             f"{self._job_dir.resolve()!s}/{self._executable} "
             f"{self._arguments} run "
             f"{self._yaml_config_path.resolve()!s} "
-            f"--verbose {self._verbose} "
+            f"{verbose_args} "
             f"--element"
         )
 
@@ -184,6 +194,11 @@ class GnuParallelLocalAdapter(QueueContextAdapter):
 
     def collect(self) -> str:
         """Return collect commands."""
+        verbose_args = f"--verbose {self._verbose}"
+        if self._verbose_datalad:
+            verbose_args = (
+                f"{verbose_args} --verbose-datalad {self._verbose_datalad}"
+            )
         return (
             f"#!/usr/bin/env {self._shell}\n\n"
             "# This script is auto-generated by junifer.\n\n"
@@ -193,7 +208,7 @@ class GnuParallelLocalAdapter(QueueContextAdapter):
             f"{self._job_dir.resolve()!s}/{self._executable} "
             f"{self._arguments} collect "
             f"{self._yaml_config_path.resolve()!s} "
-            f"--verbose {self._verbose}"
+            f"{verbose_args}"
         )
 
     def prepare(self) -> None:

{junifer-0.0.6.dev317 → junifer-0.0.6.dev330}/junifer/api/queue_context/htcondor_adapter.py

@@ -37,6 +37,9 @@ class HTCondorAdapter(QueueContextAdapter):
         virtual environment of any kind (default None).
     verbose : str, optional
         The level of verbosity (default "info").
+    verbose_datalad : str or None, optional
+        The level of verbosity for datalad. If None, will be the same
+        as ``verbose`` (default None).
     cpus : int, optional
         The number of CPU cores to use (default 1).
     mem : str, optional
@@ -83,6 +86,7 @@ class HTCondorAdapter(QueueContextAdapter):
         pre_collect: Optional[str] = None,
         env: Optional[dict[str, str]] = None,
         verbose: str = "info",
+        verbose_datalad: Optional[str] = None,
         cpus: int = 1,
         mem: str = "8G",
         disk: str = "1G",
@@ -99,6 +103,7 @@ class HTCondorAdapter(QueueContextAdapter):
         self._pre_collect = pre_collect
         self._check_env(env)
         self._verbose = verbose
+        self._verbose_datalad = verbose_datalad
         self._cpus = cpus
         self._mem = mem
         self._disk = disk
@@ -201,10 +206,15 @@ class HTCondorAdapter(QueueContextAdapter):
 
     def run(self) -> str:
         """Return run commands."""
+        verbose_args = f"--verbose {self._verbose} "
+        if self._verbose_datalad is not None:
+            verbose_args = (
+                f"{verbose_args} --verbose-datalad {self._verbose_datalad} "
+            )
         junifer_run_args = (
             "run "
             f"{self._yaml_config_path.resolve()!s} "
-            f"--verbose {self._verbose} "
+            f"{verbose_args}"
             "--element $(element)"
         )
         log_dir_prefix = (
@@ -246,10 +256,16 @@ class HTCondorAdapter(QueueContextAdapter):
 
     def collect(self) -> str:
         """Return collect commands."""
+        verbose_args = f"--verbose {self._verbose} "
+        if self._verbose_datalad is not None:
+            verbose_args = (
+                f"{verbose_args} --verbose-datalad {self._verbose_datalad} "
+            )
+
         junifer_collect_args = (
             "collect "
             f"{self._yaml_config_path.resolve()!s} "
-            f"--verbose {self._verbose}"
+            f"{verbose_args}"
         )
         log_dir_prefix = f"{self._log_dir.resolve()!s}/junifer_collect"
         fixed = (

{junifer-0.0.6.dev317 → junifer-0.0.6.dev330}/junifer/cli/cli.py

@@ -45,6 +45,32 @@ __all__ = [
 ]
 
 
+def _validate_optional_verbose(
+    ctx: click.Context, param: str, value: Optional[str]
+):
+    """Validate optional verbose option.
+
+    Parameters
+    ----------
+    ctx : click.Context
+        The context of the command.
+    param : str
+        The parameter to validate.
+    value : str
+        The value to validate.
+
+    Returns
+    -------
+    str or int or None
+        The validated value.
+
+    """
+    if value is None:
+        return value
+    else:
+        return _validate_verbose(ctx, param, value)
+
+
 def _validate_verbose(
     ctx: click.Context, param: str, value: str
 ) -> Union[str, int]:
@@ -105,8 +131,17 @@ def cli() -> None:  # pragma: no cover
     callback=_validate_verbose,
     default="info",
 )
+@click.option(
+    "--verbose-datalad",
+    type=click.UNPROCESSED,
+    callback=_validate_optional_verbose,
+    default=None,
+)
 def run(
-    filepath: click.Path, element: tuple[str], verbose: Union[str, int]
+    filepath: click.Path,
+    element: tuple[str],
+    verbose: Union[str, int],
+    verbose_datalad: Optional[Union[str, int]],
 ) -> None:
     """Run feature extraction.
 
@@ -120,10 +155,12 @@ def run(
         The element(s) to operate on.
     verbose : click.Choice
         The verbosity level: warning, info or debug (default "info").
+    verbose_datalad : click.Choice or None
+        The verbosity level for datalad: warning, info or debug (default None).
 
     """
     # Setup logging
-    configure_logging(level=verbose)
+    configure_logging(level=verbose, level_datalad=verbose_datalad)
     # TODO(synchon): add validation
     # Parse YAML
     config = parse_yaml(filepath)
@@ -167,7 +204,17 @@ def run(
     callback=_validate_verbose,
     default="info",
 )
-def collect(filepath: click.Path, verbose: Union[str, int]) -> None:
+@click.option(
+    "--verbose-datalad",
+    type=click.UNPROCESSED,
+    callback=_validate_optional_verbose,
+    default=None,
+)
+def collect(
+    filepath: click.Path,
+    verbose: Union[str, int],
+    verbose_datalad: Union[str, int, None],
+) -> None:
     """Collect extracted features.
 
     \f
@@ -178,10 +225,12 @@ def collect(filepath: click.Path, verbose: Union[str, int]) -> None:
         The filepath to the configuration file.
     verbose : click.Choice
         The verbosity level: warning, info or debug (default "info").
+    verbose_datalad : click.Choice or None
+        The verbosity level for datalad: warning, info or debug (default None).
 
     """
     # Setup logging
-    configure_logging(level=verbose)
+    configure_logging(level=verbose, level_datalad=verbose_datalad)
     # TODO: add validation
     # Parse YAML
     config = parse_yaml(filepath)
@@ -208,12 +257,19 @@ def collect(filepath: click.Path, verbose: Union[str, int]) -> None:
     callback=_validate_verbose,
     default="info",
 )
+@click.option(
+    "--verbose-datalad",
+    type=click.UNPROCESSED,
+    callback=_validate_optional_verbose,
+    default=None,
+)
 def queue(
     filepath: click.Path,
     element: tuple[str],
     overwrite: bool,
     submit: bool,
     verbose: Union[str, int],
+    verbose_datalad: Union[str, int, None],
 ) -> None:
     """Queue feature extraction.
 
@@ -231,6 +287,8 @@ def queue(
         Whether to submit the job.
     verbose : click.Choice
         The verbosity level: warning, info or debug (default "info").
+    verbose_datalad : click.Choice or None
+        The verbosity level for datalad: warning, info or debug (default None).
 
     Raises
     ------
@@ -239,7 +297,7 @@ def queue(
 
     """
     # Setup logging
-    configure_logging(level=verbose)
+    configure_logging(level=verbose, level_datalad=verbose_datalad)
     # TODO: add validation
     # Parse YAML
     config = parse_yaml(filepath)  # type: ignore
@@ -365,9 +423,16 @@ def selftest(subpkg: str) -> None:
     callback=_validate_verbose,
     default="info",
 )
+@click.option(
+    "--verbose-datalad",
+    type=click.UNPROCESSED,
+    callback=_validate_optional_verbose,
+    default=None,
+)
 def reset(
     filepath: click.Path,
     verbose: Union[str, int],
+    verbose_datalad: Union[str, int, None],
 ) -> None:
     """Reset generated assets.
 
@@ -379,10 +444,12 @@ def reset(
         The filepath to the configuration file.
     verbose : click.Choice
         The verbosity level: warning, info or debug (default "info").
+    verbose_datalad : click.Choice or None
+        The verbosity level for datalad: warning, info or debug (default None).
 
     """
     # Setup logging
-    configure_logging(level=verbose)
+    configure_logging(level=verbose, level_datalad=verbose_datalad)
     # Parse YAML
     config = parse_yaml(filepath)
     # Perform operation
@@ -409,11 +476,18 @@ def reset(
     callback=_validate_verbose,
     default="info",
 )
+@click.option(
+    "--verbose-datalad",
+    type=click.UNPROCESSED,
+    callback=_validate_optional_verbose,
+    default=None,
+)
 def list_elements(
     filepath: click.Path,
     element: tuple[str],
     output_file: Optional[click.Path],
     verbose: Union[str, int],
+    verbose_datalad: Union[str, int, None],
 ) -> None:
     """List elements of a dataset.
 
@@ -430,10 +504,12 @@ def list_elements(
         stdout is not performed.
     verbose : click.Choice
         The verbosity level: warning, info or debug (default "info").
+    verbose_datalad : click.Choice or None
+        The verbosity level for datalad: warning, info or debug (default None
 
     """
     # Setup logging
-    configure_logging(level=verbose)
+    configure_logging(level=verbose, level_datalad=verbose_datalad)
     # Parse YAML
     config = parse_yaml(filepath)
     # Fetch datagrabber

{junifer-0.0.6.dev317 → junifer-0.0.6.dev330}/junifer/markers/falff/_afni_falff.py

@@ -50,7 +50,7 @@ class AFNIALFF(metaclass=Singleton):
     @lru_cache(maxsize=None, typed=True)
     def compute(
         self,
-        data: "Nifti1Image",
+        input_path: Path,
         highpass: float,
         lowpass: float,
         tr: Optional[float],
@@ -59,8 +59,8 @@ class AFNIALFF(metaclass=Singleton):
 
         Parameters
         ----------
-        data : 4D Niimg-like object
-            Images to process.
+        input_path : pathlib.Path
+            Path to the input data.
         highpass : positive float
             Highpass cutoff frequency.
         lowpass : positive float
@@ -82,19 +82,17 @@ class AFNIALFF(metaclass=Singleton):
         """
         logger.debug("Creating cache for ALFF computation via AFNI")
 
-        # Create component-scoped tempdir
-        tempdir = WorkDirManager().get_tempdir(prefix="afni_alff+falff")
-
-        # Save target data to a component-scoped tempfile
-        nifti_in_file_path = tempdir / "input.nii"  # needs to be .nii
-        nib.save(data, nifti_in_file_path)
+        # Create element-scoped tempdir
+        element_tempdir = WorkDirManager().get_element_tempdir(
+            prefix="afni_lff"
+        )
 
         # Set 3dRSFC command
-        alff_falff_out_path_prefix = tempdir / "alff_falff"
+        lff_out_path_prefix = element_tempdir / "output"
         bp_cmd = [
             "3dRSFC",
-            f"-prefix {alff_falff_out_path_prefix.resolve()}",
-            f"-input {nifti_in_file_path.resolve()}",
+            f"-prefix {lff_out_path_prefix.resolve()}",
+            f"-input {input_path.resolve()}",
             f"-band {highpass} {lowpass}",
             "-no_rsfa -nosat -nodetrend",
         ]
@@ -104,49 +102,48 @@ class AFNIALFF(metaclass=Singleton):
         # Call 3dRSFC
         run_ext_cmd(name="3dRSFC", cmd=bp_cmd)
 
-        # Create element-scoped tempdir so that the ALFF and fALFF maps are
-        # available later as nibabel stores file path reference for
-        # loading on computation
-        element_tempdir = WorkDirManager().get_element_tempdir(
-            prefix="afni_alff_falff"
-        )
-
+        # Read header to get output suffix
+        niimg = nib.load(input_path)
+        header = niimg.header
+        sform_code = header.get_sform(coded=True)[1]
+        if sform_code == 4:
+            output_suffix = "tlrc"
+        else:
+            output_suffix = "orig"
+        # Set params suffix
         params_suffix = f"_{highpass}_{lowpass}_{tr}"
 
         # Convert alff afni to nifti
-        alff_afni_to_nifti_out_path = (
-            element_tempdir / f"alff{params_suffix}_output.nii"
+        alff_nifti_out_path = (
+            element_tempdir / f"output_alff{params_suffix}.nii"
         )  # needs to be .nii
         convert_alff_cmd = [
             "3dAFNItoNIFTI",
-            f"-prefix {alff_afni_to_nifti_out_path.resolve()}",
-            f"{alff_falff_out_path_prefix}_ALFF+orig.BRIK",
+            f"-prefix {alff_nifti_out_path.resolve()}",
+            f"{lff_out_path_prefix}_ALFF+{output_suffix}.BRIK",
         ]
         # Call 3dAFNItoNIFTI
         run_ext_cmd(name="3dAFNItoNIFTI", cmd=convert_alff_cmd)
 
         # Convert falff afni to nifti
-        falff_afni_to_nifti_out_path = (
-            element_tempdir / f"falff{params_suffix}_output.nii"
+        falff_nifti_out_path = (
+            element_tempdir / f"output_falff{params_suffix}.nii"
         )  # needs to be .nii
         convert_falff_cmd = [
             "3dAFNItoNIFTI",
-            f"-prefix {falff_afni_to_nifti_out_path.resolve()}",
-            f"{alff_falff_out_path_prefix}_fALFF+orig.BRIK",
+            f"-prefix {falff_nifti_out_path.resolve()}",
+            f"{lff_out_path_prefix}_fALFF+{output_suffix}.BRIK",
         ]
         # Call 3dAFNItoNIFTI
         run_ext_cmd(name="3dAFNItoNIFTI", cmd=convert_falff_cmd)
 
         # Load nifti
-        alff_data = nib.load(alff_afni_to_nifti_out_path)
-        falff_data = nib.load(falff_afni_to_nifti_out_path)
-
-        # Delete tempdir
-        WorkDirManager().delete_tempdir(tempdir)
+        alff_data = nib.load(alff_nifti_out_path)
+        falff_data = nib.load(falff_nifti_out_path)
 
         return (
             alff_data,
             falff_data,
-            alff_afni_to_nifti_out_path,
-            falff_afni_to_nifti_out_path,
-        )  # type: ignore
+            alff_nifti_out_path,
+            falff_nifti_out_path,
+        )

{junifer-0.0.6.dev317 → junifer-0.0.6.dev330}/junifer/markers/falff/_junifer_falff.py

@@ -47,7 +47,7 @@ class JuniferALFF(metaclass=Singleton):
     @lru_cache(maxsize=None, typed=True)
     def compute(
         self,
-        data: "Nifti1Image",
+        input_path: Path,
         highpass: float,
         lowpass: float,
         tr: Optional[float],
@@ -56,8 +56,8 @@ class JuniferALFF(metaclass=Singleton):
 
         Parameters
         ----------
-        data : 4D Niimg-like object
-            Images to process.
+        input_path : pathlib.Path
+            Path to the input data.
         highpass : positive float
             Highpass cutoff frequency.
         lowpass : positive float
@@ -80,9 +80,10 @@ class JuniferALFF(metaclass=Singleton):
         logger.debug("Creating cache for ALFF computation via junifer")
 
         # Get scan data
-        niimg_data = data.get_fdata().copy()
+        niimg = nib.load(input_path)
+        niimg_data = niimg.get_fdata().copy()
         if tr is None:
-            tr = float(data.header["pixdim"][4])  # type: ignore
+            tr = float(niimg.header["pixdim"][4])  # type: ignore
             logger.info(f"`tr` not provided, using `tr` from header: {tr}")
 
         # Bandpass the data within the lowpass and highpass cutoff freqs
@@ -120,19 +121,17 @@ class JuniferALFF(metaclass=Singleton):
         # Calculate ALFF
         alff = numerator / np.sqrt(niimg_data.shape[-1])
         alff_data = nimg.new_img_like(
-            ref_niimg=data,
+            ref_niimg=niimg,
             data=alff,
         )
         falff_data = nimg.new_img_like(
-            ref_niimg=data,
+            ref_niimg=niimg,
             data=falff,
         )
 
-        # Create element-scoped tempdir so that the ALFF and fALFF maps are
-        # available later as nibabel stores file path reference for
-        # loading on computation
+        # Create element-scoped tempdir
         element_tempdir = WorkDirManager().get_element_tempdir(
-            prefix="junifer_alff+falff"
+            prefix="junifer_lff"
         )
         output_alff_path = element_tempdir / "output_alff.nii.gz"
         output_falff_path = element_tempdir / "output_falff.nii.gz"

{junifer-0.0.6.dev317 → junifer-0.0.6.dev330}/junifer/markers/falff/falff_base.py

@@ -146,7 +146,7 @@ class ALFFBase(BaseMarker):
             estimator = JuniferALFF()
         # Compute ALFF + fALFF
         alff, falff, alff_path, falff_path = estimator.compute(  # type: ignore
-            data=input_data["data"],
+            input_path=input_data["path"],
             highpass=self.highpass,
             lowpass=self.lowpass,
             tr=self.tr,

{junifer-0.0.6.dev317 → junifer-0.0.6.dev330}/junifer/markers/falff/tests/test_falff_spheres.py

@@ -69,7 +69,9 @@ def test_ALFFSpheres(caplog: pytest.LogCaptureFixture, tmp_path: Path) -> None:
             # Fit transform marker on data
             output = marker.fit_transform(element_data)
 
-            assert "Creating cache" in caplog.text
+            # Tests for ALFFParcels run before this with the same data and that
+            # should create the cache
+            assert "Calculating ALFF and fALFF" in caplog.text
 
             # Get BOLD output
             assert "BOLD" in output

{junifer-0.0.6.dev317 → junifer-0.0.6.dev330}/junifer/markers/reho/_afni_reho.py

@@ -50,7 +50,7 @@ class AFNIReHo(metaclass=Singleton):
     @lru_cache(maxsize=None, typed=True)
     def compute(
         self,
-        data: "Nifti1Image",
+        input_path: Path,
         nneigh: int = 27,
         neigh_rad: Optional[float] = None,
         neigh_x: Optional[float] = None,
@@ -65,8 +65,8 @@ class AFNIReHo(metaclass=Singleton):
 
         Parameters
         ----------
-        data : 4D Niimg-like object
-            Images to process.
+        input_path : pathlib.Path
+            Path to the input data.
         nneigh : {7, 19, 27}, optional
             Number of voxels in the neighbourhood, inclusive. Can be:
 
@@ -128,19 +128,17 @@ class AFNIReHo(metaclass=Singleton):
         """
         logger.debug("Creating cache for ReHo computation via AFNI")
 
-        # Create component-scoped tempdir
-        tempdir = WorkDirManager().get_tempdir(prefix="afni_reho")
-
-        # Save target data to a component-scoped tempfile
-        nifti_in_file_path = tempdir / "input.nii"  # needs to be .nii
-        nib.save(data, nifti_in_file_path)
+        # Create element-scoped tempdir
+        element_tempdir = WorkDirManager().get_element_tempdir(
+            prefix="afni_reho"
+        )
 
         # Set 3dReHo command
-        reho_out_path_prefix = tempdir / "reho"
+        reho_out_path_prefix = element_tempdir / "output"
         reho_cmd = [
             "3dReHo",
             f"-prefix {reho_out_path_prefix.resolve()}",
-            f"-inset {nifti_in_file_path.resolve()}",
+            f"-inset {input_path.resolve()}",
         ]
         # Check ellipsoidal / cuboidal volume arguments
         if neigh_rad:
@@ -164,28 +162,28 @@ class AFNIReHo(metaclass=Singleton):
         # Call 3dReHo
         run_ext_cmd(name="3dReHo", cmd=reho_cmd)
 
-        # Create element-scoped tempdir so that the ReHo map is
-        # available later as nibabel stores file path reference for
-        # loading on computation
-        element_tempdir = WorkDirManager().get_element_tempdir(
-            prefix="afni_reho"
-        )
+        # Read header to get output suffix
+        niimg = nib.load(input_path)
+        header = niimg.header
+        sform_code = header.get_sform(coded=True)[1]
+        if sform_code == 4:
+            output_suffix = "tlrc"
+        else:
+            output_suffix = "orig"
+
         # Convert afni to nifti
-        reho_afni_to_nifti_out_path = (
+        reho_nifti_out_path = (
             element_tempdir / "output.nii"  # needs to be .nii
         )
         convert_cmd = [
             "3dAFNItoNIFTI",
-            f"-prefix {reho_afni_to_nifti_out_path.resolve()}",
-            f"{reho_out_path_prefix}+orig.BRIK",
+            f"-prefix {reho_nifti_out_path.resolve()}",
+            f"{reho_out_path_prefix}+{output_suffix}.BRIK",
         ]
         # Call 3dAFNItoNIFTI
         run_ext_cmd(name="3dAFNItoNIFTI", cmd=convert_cmd)
 
         # Load nifti
-        output_data = nib.load(reho_afni_to_nifti_out_path)
-
-        # Delete tempdir
-        WorkDirManager().delete_tempdir(tempdir)
+        output_data = nib.load(reho_nifti_out_path)
 
-        return output_data, reho_afni_to_nifti_out_path  # type: ignore
+        return output_data, reho_nifti_out_path