junifer 0.0.6.dev317__py3-none-any.whl → 0.0.6.dev324__py3-none-any.whl

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.dev324'
+__version_tuple__ = version_tuple = (0, 0, 6, 'dev324')

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/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/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/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/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

junifer/markers/reho/_junifer_reho.py

@@ -48,15 +48,15 @@ class JuniferReHo(metaclass=Singleton):
     @lru_cache(maxsize=None, typed=True)
     def compute(
         self,
-        data: "Nifti1Image",
+        input_path: Path,
         nneigh: int = 27,
     ) -> tuple["Nifti1Image", Path]:
         """Compute ReHo map.
 
         Parameters
         ----------
-        data : 4D Niimg-like object
-            Images to process.
+        input_path : pathlib.Path
+            Path to the input data.
         nneigh : {7, 19, 27, 125}, optional
             Number of voxels in the neighbourhood, inclusive. Can be:
 
@@ -89,7 +89,8 @@ class JuniferReHo(metaclass=Singleton):
         logger.debug("Creating cache for ReHo computation via junifer")
 
         # Get scan data
-        niimg_data = data.get_fdata()
+        niimg = nib.load(input_path)
+        niimg_data = niimg.get_fdata().copy()
         # Get scan dimensions
         n_x, n_y, n_z, _ = niimg_data.shape
 
@@ -119,7 +120,7 @@ class JuniferReHo(metaclass=Singleton):
         # after #299 is merged
         # Calculate whole brain mask
         mni152_whole_brain_mask = nmask.compute_brain_mask(
-            target_img=data,
+            target_img=niimg,
             threshold=0.5,
             mask_type="whole-brain",
         )
@@ -227,7 +228,7 @@ class JuniferReHo(metaclass=Singleton):
 
         # Create new image like target image
         output_data = nimg.new_img_like(
-            ref_niimg=data,
+            ref_niimg=niimg,
             data=reho_map,
             copy_header=False,
         )

junifer/markers/reho/reho_base.py

@@ -125,7 +125,7 @@ class ReHoBase(BaseMarker):
             estimator = JuniferReHo()
         # Compute reho
         reho_map, reho_map_path = estimator.compute(  # type: ignore
-            data=input_data["data"],
+            input_path=input_data["path"],
             **reho_params,
         )
 

{junifer-0.0.6.dev317.dist-info → junifer-0.0.6.dev324.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: junifer
-Version: 0.0.6.dev317
+Version: 0.0.6.dev324
 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.dev317.dist-info → junifer-0.0.6.dev324.dist-info}/RECORD

@@ -1,6 +1,6 @@
 junifer/__init__.py,sha256=2McgH1yNue6Z1V26-uN_mfMjbTcx4CLhym-DMBl5xA4,266
 junifer/__init__.pyi,sha256=SsTvgq2Dod6UqJN96GH1lCphH6hJQQurEJHGNhHjGUI,508
-junifer/_version.py,sha256=ysdlII3ZUlFBiHIzSCvasFtc1CMEpYMsiJEwLoFxjKQ,428
+junifer/_version.py,sha256=Mek99XWm9XI8ealUwe-651JH0dQXcrlgJtYJdIejdA4,428
 junifer/conftest.py,sha256=PWYkkRDU8ly2lYwv7VBKMHje4et6HX7Yey3Md_I2KbA,613
 junifer/py.typed,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 junifer/stats.py,sha256=e9aaagMGtgpRfW3Wdpz9ocpnYld1IWylCDcjFUgX9Mk,6225
@@ -208,14 +208,14 @@ junifer/markers/complexity/tests/test_sample_entropy.py,sha256=rfbiguVq7CUwYIvYB
 junifer/markers/complexity/tests/test_weighted_perm_entropy.py,sha256=yDWKEaUbxrnrG6J2NlktLfwSBre5OuXd63kEof7t8PM,2373
 junifer/markers/falff/__init__.py,sha256=qxdx_3FsVrn7h3gtbocK0ZmvqZwPQZGKuVkPm31ejNM,217
 junifer/markers/falff/__init__.pyi,sha256=X-q2zBjUX0imQ37yN2Cg5gKfDvq8sh_9y2hRH4g5ufY,120
-junifer/markers/falff/_afni_falff.py,sha256=D_hhT-rfWXeIe-Pwh4pK2O7N9Q07VyMOUHeAz6Ooudo,4518
-junifer/markers/falff/_junifer_falff.py,sha256=TnZ6T8u0f0YitJ7KRk4gMcC1tv6UW8m4VVcvcjGRsL0,4450
-junifer/markers/falff/falff_base.py,sha256=XZcBfXA74tWrW--uCALkRLl1tB7toR3_HfzPJUblz8g,4909
+junifer/markers/falff/_afni_falff.py,sha256=PYkSOFMyaHoGYDvmBjKLW1ALyWBe7yI36JBqZ71ji2c,4223
+junifer/markers/falff/_junifer_falff.py,sha256=1PsavcopVjPtfmPZsnNi5ynl2GTfnCx9qjuKDb7YejE,4347
+junifer/markers/falff/falff_base.py,sha256=WtJTMRn_Vmv9RZaaeLeyZpCsQc8QgB3UqIPWsUI3Lh4,4915
 junifer/markers/falff/falff_parcels.py,sha256=sSb6QLaJKpL0GCTRWW3RnpOZCoy1f9lDLgJ0I_W_LlM,6017
 junifer/markers/falff/falff_spheres.py,sha256=GrakJYPB01y9BNBXM8WzWaae0mC-S06txiycvfBGcj0,6656
 junifer/markers/falff/py.typed,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 junifer/markers/falff/tests/test_falff_parcels.py,sha256=Z3n1i8dkYbdXgouUjfIif9yLv5MubBEdrtAA-a6kRcc,4349
-junifer/markers/falff/tests/test_falff_spheres.py,sha256=-VLEvFaF8CMCN_7FLYCSfP7MMjy-gm1Zgu13je5Pku8,4373
+junifer/markers/falff/tests/test_falff_spheres.py,sha256=PGsxFjMxsH8HxIHVtdcQcX8suFa18ma_12Almknzn88,4503
 junifer/markers/functional_connectivity/__init__.py,sha256=dGTn69eS7a3rylMQh_wKlO28UmYGjsoDEGu4q5sgQFA,230
 junifer/markers/functional_connectivity/__init__.pyi,sha256=qfw6WVyE65u-5NZNi0xPa8zZVtkRfFvwyl4jHH2Xl00,539
 junifer/markers/functional_connectivity/crossparcellation_functional_connectivity.py,sha256=uLdVGywmL7qrzloh1YBL4g4tPiamA47MgHF2DQH0JTU,5733
@@ -233,10 +233,10 @@ junifer/markers/functional_connectivity/tests/test_functional_connectivity_parce
 junifer/markers/functional_connectivity/tests/test_functional_connectivity_spheres.py,sha256=A9OtFdndiSGOcPHH-QLPh6qoiD03A6KjM_emwxAlPg0,4145
 junifer/markers/reho/__init__.py,sha256=WZf4A0XaRThjl8SlFOhvTLUfhTHp5koLxZgowsgTSAE,211
 junifer/markers/reho/__init__.pyi,sha256=_aFb-Ry_EP2OMU6xRL4GlfuDpSl_egHllL-fz7vXjcE,118
-junifer/markers/reho/_afni_reho.py,sha256=YBqNYipZO8EFM4Jmek_A36zr9n4nTR0RVDPpxjJKLxM,6466
-junifer/markers/reho/_junifer_reho.py,sha256=86oBH8UtWsHJJGz2-uRXAjdGNHvaSV6Xu1v7-9AJqLs,9340
+junifer/markers/reho/_afni_reho.py,sha256=SOWR5y9AYKfw1wj2Z4Wy7ckMUVTmeS378bayvPPVqqo,6225
+junifer/markers/reho/_junifer_reho.py,sha256=14ObaRa2-0JzcoYLJyhnx4bgPCpTdciDXmrn9-gPv20,9387
 junifer/markers/reho/py.typed,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-junifer/markers/reho/reho_base.py,sha256=Amji7hNtkT9D-teJ2tTzDsVHKRXETDwhICDKa8ZGehY,4071
+junifer/markers/reho/reho_base.py,sha256=Q88TbhIM4rQWdeQPLwwxwZ9DrR8l09orD1rdTkSYDtc,4077
 junifer/markers/reho/reho_parcels.py,sha256=UE1ia3uqbmTcZMc_FI625xVPLxBYvwpfrcvhekopbkI,6392
 junifer/markers/reho/reho_spheres.py,sha256=FCC2qncC85Kd82hg-MOu4T7NAKEkXHUaCcwC9taau9Y,6996
 junifer/markers/reho/tests/test_reho_parcels.py,sha256=bRtDi91qRcRYaRqqQjuSU6NuNz-KwLVCoTYo-e5VmsI,4075
@@ -341,10 +341,10 @@ junifer/utils/tests/test_config.py,sha256=7ltIXuwb_W4Mv_1dxQWyiyM10XgUAfsWKV6D_i
 junifer/utils/tests/test_fs.py,sha256=WQS7cKlKEZ742CIuiOYYpueeAhY9PqlastfDVpVVtvE,923
 junifer/utils/tests/test_helpers.py,sha256=k5qqfxK8dFyuewTJyR1Qn6-nFaYNuVr0ysc18bfPjyU,929
 junifer/utils/tests/test_logging.py,sha256=duO4ou365hxwa_kwihFtKPLaL6LC5XHiyhOijrrngbA,8009
-junifer-0.0.6.dev317.dist-info/AUTHORS.rst,sha256=rmULKpchpSol4ExWFdm-qu4fkpSZPYqIESVJBZtGb6E,163
-junifer-0.0.6.dev317.dist-info/LICENSE.md,sha256=MqCnOBu8uXsEOzRZWh9EBVfVz-kE9NkXcLCrtGXo2yU,34354
-junifer-0.0.6.dev317.dist-info/METADATA,sha256=aQiIDdr_Ov-v0da-VOALNQiJe1FIGM-Xz32tTTqlW00,8429
-junifer-0.0.6.dev317.dist-info/WHEEL,sha256=PZUExdf71Ui_so67QXpySuHtCi3-J3wvF4ORK6k_S8U,91
-junifer-0.0.6.dev317.dist-info/entry_points.txt,sha256=6O8ru0BP-SP7YMUZiizFNoaZ2HvJpadO2G7nKk4PwjI,48
-junifer-0.0.6.dev317.dist-info/top_level.txt,sha256=4bAq1R2QFQ4b3hohjys2JBvxrl0GKk5LNFzYvz9VGcA,8
-junifer-0.0.6.dev317.dist-info/RECORD,,
+junifer-0.0.6.dev324.dist-info/AUTHORS.rst,sha256=rmULKpchpSol4ExWFdm-qu4fkpSZPYqIESVJBZtGb6E,163
+junifer-0.0.6.dev324.dist-info/LICENSE.md,sha256=MqCnOBu8uXsEOzRZWh9EBVfVz-kE9NkXcLCrtGXo2yU,34354
+junifer-0.0.6.dev324.dist-info/METADATA,sha256=a_1zgab7_sp4jn5_jM6ZrKnJXL1InFZpK9zj39EUtVY,8429
+junifer-0.0.6.dev324.dist-info/WHEEL,sha256=PZUExdf71Ui_so67QXpySuHtCi3-J3wvF4ORK6k_S8U,91
+junifer-0.0.6.dev324.dist-info/entry_points.txt,sha256=6O8ru0BP-SP7YMUZiizFNoaZ2HvJpadO2G7nKk4PwjI,48
+junifer-0.0.6.dev324.dist-info/top_level.txt,sha256=4bAq1R2QFQ4b3hohjys2JBvxrl0GKk5LNFzYvz9VGcA,8
+junifer-0.0.6.dev324.dist-info/RECORD,,