pyLQT 1.0.dev1__py3-none-any.whl → 1.0.dev2__py3-none-any.whl

pyLQT/__init__.py

@@ -0,0 +1,3 @@
+from .pipeline import run_pyLQT
+
+__all__ = ["run_pyLQT"]

pyLQT/cli.py

@@ -0,0 +1,71 @@
+import argparse
+
+from pipeline import run_pyLQT
+
+
+def main():
+    """Run the pyLQT pipeline through the command-line interface."""
+
+    parser = argparse.ArgumentParser(
+        prog="pyLQT",
+        description="Run the pyLQT lesion quantification pipeline."
+    )
+
+    parser.add_argument(
+        "--lesion_cfg",
+        required=True,
+        help='[Path] Path to lesion config'
+    )
+    
+    parser.add_argument(
+        "--atlas_cfg",
+        default=None,
+        help='[Path, optional] Path to atlas config. If none provided, a default atlas is downloaded.'
+    )
+    
+    parser.add_argument(
+        "--tracts_cfg",
+        default=None,
+        help='[Path, optional] Path to tracts config. If none provided, default tracts are downloaded.'
+    )
+    
+    parser.add_argument(
+        "--out_folder",
+        required=True,
+        help='[Path] Output folder path'
+    )
+    
+    parser.add_argument(
+        "--registration_type",
+        choices=['affine', 'syn'],
+        default='syn',
+        help='[str, default=syn] Registration type, from lesion space and atlas space to tract space'
+    )
+    
+    parser.add_argument(
+        "--spared_threshold",
+        default=0.5,
+        help='[float (0,1), default=0.5] Defines threshold for spared parcel-to-parcel connectivity'
+    )
+    
+    parser.add_argument(
+        "--delete_tmp",
+        default=False,
+        help='[bool, default=False] Delete temporary folder'
+    )
+
+    args = parser.parse_args()
+
+    run_pyLQT(
+        lesion_cfg=args.lesion_cfg,
+        atlas_cfg=args.atlas_cfg,
+        tracts_cfg=args.tracts_cfg,
+        out_folder=args.out_folder,
+        registration_type=args.registration_type,
+        spared_threshold=args.spared_threshold,
+        delete_tmp=args.delete_tmp,
+    )
+
+
+if __name__ == "__main__":
+    main()

{lesion_metrics → pyLQT/lesion_metrics}/gray_matter/metrics.py

@@ -14,7 +14,7 @@ def get_parcel_damage(
     ):
     """Compute lesion damage for each parcel of an atlas.
 
-    The lesion is resampled to the atlas space, and the proportion of
+    The lesion is resampled to the atlas grid, and the proportion of
     each parcel overlapping with the lesion is computed. A parcel damage 
     map is created, where each parcel voxel is assigned the damage proportion 
     of its parcel.

{lesion_metrics → pyLQT/lesion_metrics}/utils/downloads.py

@@ -6,6 +6,7 @@ import shutil
 import gzip
 import nibabel as nib
 import numpy as np
+import pandas as pd
 
 from templateflow import api as tf
 from pathlib import Path
@@ -135,7 +136,72 @@ def uncompress_trk_gz(
 
     return out_path
 
-def default_download(out_folder: Path):
+def default_atlas_download(out_folder: Path):
+    atlas_cfg_folder = out_folder / "glasser360_config"
+    atlas_cfg_folder.mkdir(parents=True, exist_ok=True)
+    
+    atlas_path = atlas_cfg_folder / "glasser360_MNI152NLin2009cAsym.nii.gz"
+    atlas_labels_path = atlas_cfg_folder / "glasser360_labels.csv"
+    atlas_labels_path_txt = atlas_labels_path.with_suffix(".txt")
+    atlas_cfg_path = atlas_cfg_folder / "glasser360_config.json"
+    
+    if atlas_path.exists():
+        warnings.warn(
+            f"atlas already exists, using it: {atlas_path}"
+        )
+    else:
+        url_atlas = "https://ndownloader.figshare.com/files/43043653"
+        
+        response = requests.get(url_atlas, stream=True, timeout=60)
+        response.raise_for_status()
+        
+        with atlas_path.open("wb") as f:
+            for chunk in response.iter_content(chunk_size=1024*1024):
+                if chunk:
+                    f.write(chunk)
+
+    if atlas_labels_path.exists():
+        warnings.warn(
+            f"atlas labels already exist, using it: {atlas_labels_path}"
+        )
+    else:
+        url_labels =  "https://ndownloader.figshare.com/files/43041079"
+        
+        response = requests.get(url_labels, stream=True, timeout=60)
+        response.raise_for_status()
+        
+        with atlas_labels_path_txt.open("wb") as f:
+            for chunk in response.iter_content(chunk_size=1024*1024):
+                if chunk:
+                    f.write(chunk)
+
+        with open(atlas_labels_path_txt, "r") as f:
+            labels = [line.strip() for line in f if line.strip()]
+
+        atlas_labels = pd.DataFrame({
+            "Parcel": range(1, len(labels) + 1),
+            "Label": labels
+        })
+
+        atlas_labels.to_csv(atlas_labels_path, index=False)
+        atlas_labels_path_txt.unlink()
+
+    config = {
+        "name": "glasser360",
+        "atlas": "glasser360_MNI152NLin2009cAsym.nii.gz",
+        "atlas_labels": "glasser360_labels.csv",
+        "space": "MNI152NLin2009cAsym"
+    }
+            
+    atlas_cfg = Path(atlas_cfg_folder / "glasser360_config.json")
+            
+    with open(atlas_cfg, "w") as f:
+        json.dump(config, f, indent=4)
+        
+    return atlas_cfg
+    
+
+def default_tracts_download(out_folder: Path):
     """Download the default HCP1065 tractography dataset and create its configuration.
 
     The function downloads the HCP1065 average tractograms and their
@@ -162,7 +228,9 @@ def default_download(out_folder: Path):
         
     if tracts_folder.exists():
         warnings.warn(
-            f"Tracts folder already exists, using it: {tracts_folder}"
+            f"Tracts folder already exists, using it: {tracts_folder}",
+            UserWarning,
+            stacklevel=2
         )
     else:
         url = "https://github.com/data-others/atlas/releases/download/hcp1065/hcp1065_avg_tracts_trk.zip"
@@ -185,7 +253,9 @@ def default_download(out_folder: Path):
 
     if references_folder.exists():
         warnings.warn(
-            f"References folder already exists, using it: {references_folder}"
+            f"References folder already exists, using it: {references_folder}",
+            UserWarning,
+            stacklevel=2
         )
     else:
         url = "https://github.com/data-others/atlas/releases/download/hcp1065/hcp1065_avg_tracts_nifti.zip"
@@ -200,14 +270,13 @@ def default_download(out_folder: Path):
                     f.write(chunk)
 
         with zipfile.ZipFile(zip_path) as zip_file:
-            zip_file.extractall(zip_path.parent)
+            zip_file.extractall(tracts_cfg_folder)
 
         zip_path.unlink()
 
-        references_folder = zip_path.parent / "nifti"
-        references_folder.rename(
-            zip_path.parent / "hcp1065_avg_tracts"
-        )
+        extracted_niftis = tracts_cfg_folder / "nifti"
+        if extracted_niftis.exists():
+            extracted_niftis.rename(references_folder)
             
     tracts_name = "HCP1065"
     tracts_space = "MNI152NLin2009aAsym"
@@ -247,4 +316,8 @@ def get_ref_space_tf(space: str, res: int = 1):
     ref_space_volume[brain_mask.get_fdata() <= 0] = 0
     ref_space = nib.Nifti1Image(ref_space_volume, t1w_ref.affine, t1w_ref.header)
     
-    return ref_space
+    return ref_space
+
+if __name__=="__main__":
+    
+    default_atlas_download(Path("D:/internship_myspace_lab/pyLQT/test_download_atlas"))

{lesion_metrics → pyLQT/lesion_metrics}/utils/registration.py

@@ -26,14 +26,18 @@ def register_affine(
     AffineMap
         Affine mapping from the moving image space to the static image space.
     """
-    metric = MutualInformationMetric()
+    
+    moving_data = moving.get_fdata()
+    static_data = static.get_fdata()
+    
+    metric = MutualInformationMetric(nbins=32, sampling_proportion=0.2)
 
-    affreg = AffineRegistration(metric=metric)
+    affreg = AffineRegistration(metric=metric, level_iters=[500, 250, 100])
 
     transform = TranslationTransform3D()
     translation = affreg.optimize(
-        static.get_fdata(),
-        moving.get_fdata(),
+        static_data,
+        moving_data,
         transform,
         params0=None,
         static_grid2world=static.affine,
@@ -42,8 +46,8 @@ def register_affine(
 
     transform = RigidTransform3D()
     rigid = affreg.optimize(
-        static.get_fdata(),
-        moving.get_fdata(),
+        static_data,
+        moving_data,
         transform,
         starting_affine=translation.affine,
         params0=None,
@@ -53,8 +57,8 @@ def register_affine(
 
     transform = AffineTransform3D()
     mapping = affreg.optimize(
-        static.get_fdata(),
-        moving.get_fdata(),
+        static_data,
+        moving_data,
         transform,
         starting_affine=rigid.affine,
         params0=None,
@@ -93,7 +97,7 @@ def register_syn(
     """
     metric = CCMetric(dim=3)
 
-    syn = SymmetricDiffeomorphicRegistration(metric)
+    syn = SymmetricDiffeomorphicRegistration(metric, level_iters=[20, 10, 5])
 
     mapping = syn.optimize(static.get_fdata(),
                            moving.get_fdata(),

pyLQT/pipeline.py

@@ -0,0 +1,183 @@
+from lesion_metrics.white_matter.metrics import get_tracts_disconnections_severity, get_tracts_disconnections_map
+from lesion_metrics.white_matter.utils import get_streamlines_disconnection, get_aggregated_tracts, clip_tractogram_to_grid
+from lesion_metrics.network.metrics import get_sspl_increase, get_parcel_wise_disconnection_severity
+from lesion_metrics.gray_matter.metrics import get_parcel_damage
+from lesion_metrics.utils.configs import atlas_cfg_read, tracts_cfg_read, lesion_cfg_read
+from lesion_metrics.utils.downloads import default_tracts_download, default_atlas_download, create_tract_list, uncompress_trk_gz, get_ref_space_tf
+from lesion_metrics.utils.registration import register_affine, register_syn
+from pathlib import Path
+
+import nibabel as nib
+import shutil
+import tempfile
+
+def run_pyLQT(
+    lesion_cfg,
+    out_folder,
+    atlas_cfg=None,
+    tracts_cfg=None,
+    registration_type="syn",
+    spared_threshold=0.5,
+    delete_tmp=False,
+):
+    """Run the pyLQT lesion quantification pipeline.
+
+    Parameters
+    ----------
+    lesion_cfg : str or pathlib.Path
+        Path to the lesion configuration file.
+    out_folder : str or pathlib.Path
+        Output directory.
+    atlas_cfg : str or pathlib.Path, optional
+        Path to the atlas configuration file. If None, the default atlas
+        is downloaded.
+    tracts_cfg : str or pathlib.Path, optional
+        Path to the tractography configuration file. If None, the default
+        tractography dataset is downloaded.
+    registration_type : {"affine", "syn"}, default="syn"
+        Registration method used to map the lesion and atlas to tract space.
+    spared_threshold : float, default=0.5
+        Threshold in (0, 1) for spared parcel-to-parcel connectivity.
+    delete_tmp : bool, default=False
+        Whether to delete the temporary directory after processing.
+    """
+
+    if registration_type not in ("affine", "syn"):
+        raise ValueError("registration_type must be 'affine' or 'syn'")
+
+    spared_threshold = float(spared_threshold)
+    if not 0 < spared_threshold < 1:
+        raise ValueError("spared_threshold must be between 0 and 1")
+
+    out_folder = Path(out_folder)
+    out_folder.mkdir(parents=True, exist_ok=True)
+    tmp_folder = Path(tempfile.mkdtemp(prefix="LQT_tmp_", dir=out_folder))
+
+    lesion_cfg_path = Path(lesion_cfg)
+    atlas_cfg_path = (
+        default_atlas_download(out_folder)
+        if atlas_cfg is None
+        else Path(atlas_cfg)
+    )
+    tracts_cfg_path = (
+        default_tracts_download(out_folder)
+        if tracts_cfg is None
+        else Path(tracts_cfg)
+    )
+
+    lesion_cohort_name, lesion_folder, lesion_space = lesion_cfg_read(lesion_cfg_path)
+    atlas_name, atlas, atlas_labels, atlas_space = atlas_cfg_read(atlas_cfg_path)
+    tracts_name, tracts_folder, references_folder, tracts_space = (tracts_cfg_read(tracts_cfg_path))
+
+    if lesion_space.endswith((".nii.gz", ".nii")):
+        lesion_ref_space = nib.load(lesion_space)
+    else:
+        lesion_ref_space = get_ref_space_tf(lesion_space)
+
+    if atlas_space.endswith((".nii.gz", ".nii")):
+        atlas_ref_space = nib.load(atlas_space)
+    else:
+        atlas_ref_space = get_ref_space_tf(atlas_space)
+
+    if tracts_space.endswith((".nii.gz", ".nii")):
+        tracts_ref_space = nib.load(tracts_space)
+    else:
+        tracts_ref_space = get_ref_space_tf(tracts_space)
+
+    registrator = (register_affine if registration_type == "affine" else register_syn)
+
+    lesion_tracts_mapping = registrator(lesion_ref_space, tracts_ref_space)
+    atlas_tracts_mapping = registrator(atlas_ref_space, tracts_ref_space)
+
+    atlas_trans_volume = atlas_tracts_mapping.transform(atlas.get_fdata(), interpolation="nearest")
+    atlas_trans = nib.Nifti1Image(
+        atlas_trans_volume,
+        tracts_ref_space.affine,
+        tracts_ref_space.header,
+    )
+    nib.save(atlas_trans, tmp_folder / "atlas_trans.nii.gz")
+
+    tracts_pathlist_gz, references_pathlist = create_tract_list(
+        references_folder=references_folder,
+        tracts_folder=tracts_folder,
+    )
+
+    tracts_pathlist = []
+    for tract_path in tracts_pathlist_gz:
+        relative_path = tract_path.relative_to(tracts_folder)
+        out_path = tmp_folder / relative_path.with_suffix("")
+        out_path.parent.mkdir(parents=True, exist_ok=True)
+
+        tracts_pathlist.append(uncompress_trk_gz(tract_path, out_path))
+
+    aggregated_tracts_path = tmp_folder / "aggregated_tracts.trk"
+    aggregated_tracts = get_aggregated_tracts(
+        tracts_pathlist=tracts_pathlist,
+        out_path=aggregated_tracts_path,
+    )
+
+    for lesion_path in lesion_folder.glob("*.nii*"):
+        lesion_name = lesion_path.name.removesuffix(".nii.gz").removesuffix(".nii")
+        out_dir = out_folder / lesion_name
+        out_dir.mkdir(parents=True, exist_ok=True)
+
+        tmp_dir = tmp_folder / lesion_name
+        tmp_dir.mkdir(parents=True, exist_ok=True)
+
+        lesion = nib.load(lesion_path)
+        lesion_trans_volume = lesion_tracts_mapping.transform(lesion.get_fdata(), interpolation="linear")
+        lesion_trans = nib.Nifti1Image(
+            lesion_trans_volume,
+            tracts_ref_space.affine,
+            tracts_ref_space.header,
+        )
+        nib.save(lesion_trans, tmp_dir / "lesion_trans.nii.gz")
+
+        aggregated_tracts_disconnected_path = tmp_dir / "aggregated_disconnected.trk"
+        aggregated_tracts_disconnected, _ = get_streamlines_disconnection(
+            tract=aggregated_tracts,
+            lesion=lesion_trans,
+            out_path=aggregated_tracts_disconnected_path,
+        )
+        aggregated_tracts_disconnected = clip_tractogram_to_grid(aggregated_tracts_disconnected)
+
+        get_parcel_damage(
+            lesion=lesion_trans,
+            atlas=atlas_trans,
+            atlas_labels=atlas_labels,
+            parcel_column="Parcel",
+            out_path=out_dir / "parcel_damage",
+        )
+
+        get_tracts_disconnections_severity(
+            tracts_path=tracts_pathlist,
+            references_path=references_pathlist,
+            lesion=lesion_trans,
+            disc_tracts_out_path=tmp_dir / "wm_disconnections_severity",
+            disc_frame_out_path=out_dir / "wm_disconnections_severity",
+        )
+
+        get_tracts_disconnections_map(
+            tracts=aggregated_tracts,
+            tracts_disconnected=aggregated_tracts_disconnected,
+            out_path=out_dir / "wm_disconnections_map",
+        )
+
+        structural_connectivity, _, _, per_spared_connections = (
+            get_parcel_wise_disconnection_severity(
+                aggregated_tracts=aggregated_tracts,
+                aggregated_tracts_disconnected=aggregated_tracts_disconnected,
+                atlas=atlas_trans,
+                out_path=out_dir / "parcel_wise_disconnection_severity",
+            )
+        )
+
+        get_sspl_increase(
+            structural_connectivity=structural_connectivity,
+            per_spared_connections=per_spared_connections,
+            out_path=out_dir / "sspl_increase",
+            spared_threshold=spared_threshold,
+        )
+
+    if delete_tmp:
+        shutil.rmtree(tmp_folder)

pylqt-1.0.dev2.dist-info/METADATA

@@ -0,0 +1,136 @@
+Metadata-Version: 2.4
+Name: pyLQT
+Version: 1.0.dev2
+Summary: Python version of LQT from Griffis et al. 2021 (https://www.sciencedirect.com/science/article/pii/S2213158221000838)
+Author-email: Antonin HUDRY <hudryantoninm@gmail.com>
+License-Expression: GPL-3.0-only
+Keywords: LQT,lqt,lesion,brain,
+Classifier: Development Status :: 2 - Pre-Alpha
+Classifier: Intended Audience :: Science/Research
+Classifier: Intended Audience :: Healthcare Industry
+Classifier: Programming Language :: Python
+Requires-Python: >=3.10
+Description-Content-Type: text/markdown
+License-File: LICENSE
+Requires-Dist: numpy<2.0
+Requires-Dist: nibabel>5.0
+Requires-Dist: pandas>3.0
+Requires-Dist: scipy>1.15
+Requires-Dist: requests>2.0
+Requires-Dist: dipy>1.10
+Requires-Dist: templateflow
+Requires-Dist: networkx>3.6.1
+Dynamic: license-file
+
+# pyLQT
+
+pyLQT, standing for “Python Lesion Quantification Toolkit”, is a toolbox designed for analysis brain-lesioned patients scans with minimal resources, regrouping several metrics for lesion quantification. It is a sucessor to LQT [[1]](#ref1) which was originally developped on MATLAB and remastered on R [[2]](#ref2).
+
+The metrics are divided into three main domains:
+
+* **Grey matter:** lesion effects on grey-matter regions and parcels.
+* **White matter:** lesion effects on white-matter tracts and streamlines.
+* **Network:** lesion effects on structural or functional brain networks.
+
+## installation
+
+```bash
+pip install pyLQT
+```
+
+## Running pyLQT
+
+### Configuration files
+
+pyLQT uses three JSON configuration files to define the lesion cohort, atlas, and tractography dataset. Each file specifies the relevant data paths and reference space. For the simplest usage, only the lesion cohort configuration file is required, in which case the atlas and tractogram files are automatically downloaded.
+
+##### Lesion configuration
+
+Defines the lesion cohort, the folder containing the lesion NIfTI files, and their reference space.
+
+```json
+{
+    "cohort_name": "ExampleCohort",
+    "lesion_folder": "lesions",
+    "space": "MNI152NLin6Asym"
+}
+```
+
+##### Atlas configuration
+
+Specifies the atlas NIfTI file, its corresponding label CSV, and its reference space. The CSV must contain a `Parcel` column with parcel indices matching the atlas NIfTI and a `Label` column with the corresponding parcel names.
+
+```json
+{
+    "name": "glasser360",
+    "atlas": "glasser360.nii.gz",
+    "atlas_labels": "glasser360_labels.csv",
+    "space": "MNI152NLin2009cAsym"
+}
+```
+
+##### Tractography configuration
+
+Specifies the tractography dataset, the folder containing its `.trk` files, the folder containing its reference NIfTI files, and its reference space.
+
+```json
+{
+    "name": "HCP1065",
+    "tracts": "hcp1065_avg_tracts_trk",
+    "references": "hcp1065_avg_tracts",
+    "space": "MNI152NLin2009aAsym"
+}
+```
+
+**Notes:**
+- All paths are relative to the corresponding configuration file.
+- NIfTI files can be uncompressed (`.nii`) or gzipped (`.nii.gz`).
+- The `space` field can specify either a standard space supported by TemplateFlow (e.g., `MNI152NLin2009cAsym`) or a relative path to a reference NIfTI file (`.nii` or `.nii.gz`).
+- All paths and names above are examples and should be adapted to your dataset.
+
+### basic use
+
+Below is the simplest usage that only requires the [lesion cohort configuration file](#lesion-configuration). By default, the Glasser360 atlas and the HCP1065 tractogram are downloaded if not provided in the arguments.
+
+```bash
+pyLQT --lesion_cfg <path/to/lesion_config> --out_folder <path/to/outputs>
+```
+
+
+| Argument | Default | Options / Range | Description |
+|---|---|---|---|
+| `--lesion_cfg` | — | Path to a `.json` file | Path to the lesion cohort configuration file. |
+| `--atlas_cfg` | `None` | Path to a `.json` file or `None` | Path to the atlas configuration file. |
+| `--tracts_cfg` | `None` | Path to a `.json` file or `None` | Path to the tracts configuration file. |
+| `--out_folder` | — | Path to a directory | Path to the output folder. |
+| `--registration_type` | `syn` | `syn`, `affine` | Registration type used to transform the lesion and atlas from their respective spaces to tract space. |
+| `--spared_threshold` | `0.5` | Float between `0` and `1` | Threshold used to define spared parcel-to-parcel connectivity. |
+| `--delete_tmp` | `True` | `True`, `False` | Whether to delete temporary files generated during processing. |
+
+
+## Dependencies
+
+pyLQT relies on the following libraries:
+-    `numpy`
+-    `nibabel`
+-    `pandas`
+-    `scipy`
+-    `requests`
+-    `dipy`
+-    `templateflow`
+-    `networkx`
+
+
+## Links
+
+- [Source](https://github.com/AntoninHUDRY/pyLQT)
+- [documentation]() [TODO, deploy via github]
+- [PyPI](https://pypi.org/project/pyLQT/)
+
+## References
+
+<a id="ref1"></a>
+[1] Griffis et al. (2021). [Lesion Quantification Toolkit: A MATLAB software tool for estimating grey matter damage and white matter disconnections in patients with focal brain lesions](https://doi.org/10.1016/j.nicl.2021.102639). *NeuroImage: Clinical*.
+
+<a id="ref2"></a>
+[2] Dworkin Jordan (https://github.com/jdwor/LQT)

pylqt-1.0.dev2.dist-info/RECORD

@@ -0,0 +1,23 @@
+pyLQT/__init__.py,sha256=y4DvjD5OqKInZ-wJary448UmEBajl8gNbZv50dgabwI,58
+pyLQT/cli.py,sha256=aQHoOa8ZXhFOgJbJTRbAVCCeFLPGUc4qpgVnL1QLEfA,1881
+pyLQT/pipeline.py,sha256=PRS_yzaASATGOYz6GjNhYq733szgwF8k2vPlDqeTjsE,7313
+pyLQT/lesion_metrics/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+pyLQT/lesion_metrics/gray_matter/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+pyLQT/lesion_metrics/gray_matter/metrics.py,sha256=38R3etfmmKo3hFHkYh7Jw2EUrXbW7I0Jkl8So_-5JWA,2616
+pyLQT/lesion_metrics/gray_matter/utils.py,sha256=czxyAESP0neVCue424L0hon-5Tdbsl0pN0zJHWG5HiQ,840
+pyLQT/lesion_metrics/network/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+pyLQT/lesion_metrics/network/metrics.py,sha256=Ka3G8_dIq-8H3b8mvRR96QM7p_gR2yFLKUxQOLYMWM4,8281
+pyLQT/lesion_metrics/network/utils.py,sha256=fcakfKSMortNIE3TfJP8zwanzBiT1McsYnHcTgwbxE4,3283
+pyLQT/lesion_metrics/utils/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+pyLQT/lesion_metrics/utils/configs.py,sha256=FFur04Sh48IlxUTmz4kRkyLKWffi7bvzyw8B2GclpRU,4961
+pyLQT/lesion_metrics/utils/downloads.py,sha256=dXmB_bjt3uEPjJfGThkXw0lojidTTSm_Rx08wyJYafk,10567
+pyLQT/lesion_metrics/utils/registration.py,sha256=iqD9qIPLSxsmAPkGbnvyPZBAeVGfnG5_9P0UYZ0JGIg,3307
+pyLQT/lesion_metrics/white_matter/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+pyLQT/lesion_metrics/white_matter/metrics.py,sha256=sujuHry19XpVWFNo52Ntq7h3JkWEcgsZ1-rNFr7E940,6155
+pyLQT/lesion_metrics/white_matter/utils.py,sha256=oo492yzxsPACMvHEyYhoYVqgVRrv8mpqP7FOf5OU_Ms,7087
+pylqt-1.0.dev2.dist-info/licenses/LICENSE,sha256=IwGE9guuL-ryRPEKi6wFPI_zOhg7zDZbTYuHbSt_SAk,35823
+pylqt-1.0.dev2.dist-info/METADATA,sha256=KiPMLFsDYUD4RwjwyrciKoxH4SrAcMPevyL8aecYkek,5225
+pylqt-1.0.dev2.dist-info/WHEEL,sha256=mOdXR9C3wMVov7cSY8PY-XCzo3SmmWCHvGQr5csJoao,91
+pylqt-1.0.dev2.dist-info/entry_points.txt,sha256=GHGfSFJ0r77qGtkRMJeCuNj2YHm5aynNxdi3JXxyyCc,41
+pylqt-1.0.dev2.dist-info/top_level.txt,sha256=TwID1Q-SuDZ2HjBLKSwruNTXKrRGgAh16rMZ5pjUk4U,6
+pylqt-1.0.dev2.dist-info/RECORD,,

{pylqt-1.0.dev1.dist-info → pylqt-1.0.dev2.dist-info}/WHEEL

@@ -1,5 +1,5 @@
 Wheel-Version: 1.0
-Generator: setuptools (84.0.0)
+Generator: setuptools (84.2.1)
 Root-Is-Purelib: true
 Tag: py3-none-any
 

pylqt-1.0.dev2.dist-info/entry_points.txt

@@ -0,0 +1,2 @@
+[console_scripts]
+pyLQT = pyLQT.cli:main

pylqt-1.0.dev2.dist-info/top_level.txt

@@ -0,0 +1 @@
+pyLQT

pyLQT_cli.py

@@ -1,180 +0,0 @@
-from lesion_metrics.white_matter.metrics import get_tracts_disconnections_severity, get_tracts_disconnections_map
-from lesion_metrics.white_matter.utils import get_streamlines_disconnection, get_aggregated_tracts, clip_tractogram_to_grid
-from lesion_metrics.network.metrics import get_sspl_increase, get_parcel_wise_disconnection_severity
-from lesion_metrics.gray_matter.metrics import get_parcel_damage
-from lesion_metrics.utils.configs import atlas_cfg_read, tracts_cfg_read, lesion_cfg_read
-from lesion_metrics.utils.downloads import default_download, create_tract_list, uncompress_trk_gz, get_ref_space_tf
-from lesion_metrics.utils.registration import register_affine, register_syn
-from pathlib import Path
-
-import nibabel as nib
-import shutil
-import tempfile
-import argparse
-
-
-def run_pyLQT():
-    """Run the pyLQT lesion quantification pipeline.
-
-    Parses the command-line arguments, loads the lesion, atlas, and
-    tractography configurations, performs the required registrations,
-    and computes the different lesion quantification measures.
-    """
-    
-    parser = argparse.ArgumentParser()
-    
-    parser.add_argument(
-        "--lesion_cfg",
-        required=True,
-        help='[Path] Path to lesion config'
-    )
-    
-    parser.add_argument(
-        "--atlas_cfg",
-        required=True,
-        help='[Path] Path to atlas config'
-    )
-    
-    parser.add_argument(
-        "--tracts_cfg",
-        default=None,
-        help='[Path] Path to tracts config'
-    )
-    
-    parser.add_argument(
-        "--out_folder",
-        required=True,
-        help='[Path] Output folder path'
-    )
-    
-    parser.add_argument(
-        "--registration_type",
-        choices=['affine', 'syn'],
-        default='syn',
-        help='[str, default=syn] Registration type, from lesion space and atlas space to tract space'
-    )
-    
-    parser.add_argument(
-        "--spared_threshold",
-        default=0.5,
-        help='[float (0,1), default=0.5] Defines threshold for spared parcel-to-parcel connectivity'
-    )
-    
-    parser.add_argument(
-        "--delete_tmp",
-        default=False,
-        help='[bool, default=False] Delete temporary folder'
-    )
-
-    args = parser.parse_args()
-    
-    out_folder = Path(args.out_folder)
-    out_folder.mkdir(parents=True, exist_ok=True)
-    tmp_folder = Path(tempfile.mkdtemp(prefix="LQT_tmp_", dir=out_folder))
-    
-    lesion_cfg_path = Path(args.lesion_cfg)
-    atlas_cfg_path = Path(args.atlas_cfg)
-    if not args.tracts_cfg:
-        tracts_cfg_path = default_download(out_folder)
-    else: 
-        tracts_cfg_path = Path(args.tracts_cfg)
-    
-    lesion_cohort_name, lesion_folder, lesion_space = lesion_cfg_read(lesion_cfg_path)
-    atlas_name, atlas, atlas_labels, atlas_space = atlas_cfg_read(atlas_cfg_path)
-    tracts_name, tracts_folder, references_folder, tracts_space = tracts_cfg_read(tracts_cfg_path)
-    
-    if lesion_space.endswith(".nii.gz") or lesion_space.endswith(".nii"):
-        lesion_ref_space = nib.nifti1.load(lesion_space)
-    else: lesion_ref_space = get_ref_space_tf(lesion_space)
-    if atlas_space.endswith(".nii.gz") or atlas_space.endswith(".nii"):
-        atlas_ref_space = nib.nifti1.load(atlas_space)
-    else: atlas_ref_space = get_ref_space_tf(atlas_space)
-    if tracts_space.endswith(".nii.gz") or tracts_space.endswith(".nii"):
-        tracts_ref_space = nib.nifti1.load(tracts_space)
-    else: tracts_ref_space = get_ref_space_tf(tracts_space)
-    
-    if args.registration_type == "affine":
-        registrator = register_affine
-    else:
-        registrator = register_syn
-    lesion_tracts_mapping = registrator(lesion_ref_space, tracts_ref_space)
-    atlas_tracts_mapping = registrator(atlas_ref_space, tracts_ref_space)
-    
-    atlas_trans_volume = atlas_tracts_mapping.transform(atlas.get_fdata(), interpolation='nearest')
-    atlas_trans = nib.Nifti1Image(atlas_trans_volume, tracts_ref_space.affine, tracts_ref_space.header)
-    nib.nifti1.save(atlas_trans, tmp_folder / "atlas_trans.nii.gz")
-    
-    spared_threshold = float(args.spared_threshold)
-    tracts_pathlist_gz, references_pathlist = create_tract_list(references_folder=references_folder,
-                                                                tracts_folder=tracts_folder)
-    
-    tracts_pathlist = []
-    for tract_path in tracts_pathlist_gz:
-        relative_path = tract_path.relative_to(tracts_folder)
-        out_path = tmp_folder / relative_path.with_suffix("")
-
-        out_path.parent.mkdir(parents=True, exist_ok=True)
-
-        tracts_pathlist.append(
-            uncompress_trk_gz(tract_path, out_path)
-        )
-    
-    aggregated_tracts_path = tmp_folder / "aggregated_tracts.trk"
-    aggregated_tracts = get_aggregated_tracts(tracts_pathlist=tracts_pathlist,
-                                              out_path=aggregated_tracts_path)
-    
-    for lesion_path in lesion_folder.glob("*.nii*"):
-        
-        lesion_name = lesion_path.name.removesuffix(".nii.gz").removesuffix(".nii")
-        out_dir = out_folder / lesion_name
-        out_dir.mkdir(parents=True, exist_ok=True)
-        tmp_dir = tmp_folder / lesion_name
-        tmp_dir.mkdir(parents=True, exist_ok=True)
-        
-        lesion = nib.nifti1.load(lesion_path)
-        lesion_trans_volume = lesion_tracts_mapping.transform(lesion.get_fdata(), interpolation='linear')
-        lesion_trans = nib.Nifti1Image(lesion_trans_volume, tracts_ref_space.affine, tracts_ref_space.header)
-        nib.nifti1.save(lesion_trans, tmp_dir / "lesion_trans.nii.gz")
-
-        aggregated_tracts_disconnected_path = tmp_dir / f"aggregated_disconnected.trk"
-        aggregated_tracts_disconnected, _ = get_streamlines_disconnection(tract=aggregated_tracts,
-                                                                          lesion=lesion_trans,
-                                                                          out_path=aggregated_tracts_disconnected_path)
-        aggregated_tracts_disconnected = clip_tractogram_to_grid(aggregated_tracts_disconnected)
-        
-        get_parcel_damage(lesion=lesion_trans, 
-                          atlas=atlas_trans, 
-                          atlas_labels=atlas_labels,
-                          parcel_column="Parcel",
-                          out_path=out_dir / "parcel_damage")
-        
-        get_tracts_disconnections_severity(tracts_path=tracts_pathlist,
-                                           references_path=references_pathlist,
-                                           lesion=lesion_trans,
-                                           disc_tracts_out_path=tmp_dir / "wm_disconnections_severity",
-                                           disc_frame_out_path=out_dir / "wm_disconnections_severity")
-        
-        get_tracts_disconnections_map(tracts=aggregated_tracts,
-                                      tracts_disconnected=aggregated_tracts_disconnected,
-                                      out_path=out_dir / "wm_disconnections_map")
-        
-        structural_connectivity, _, _, per_spared_connections = (
-            get_parcel_wise_disconnection_severity(
-                aggregated_tracts=aggregated_tracts,
-                aggregated_tracts_disconnected=aggregated_tracts_disconnected,
-                atlas=atlas_trans,
-                out_path=out_dir / "parcel_wise_disconnection_severity",
-            )
-        )
-        
-        get_sspl_increase(structural_connectivity=structural_connectivity,
-                          per_spared_connections=per_spared_connections,
-                          out_path=out_dir / "sspl_increase",
-                          spared_threshold=spared_threshold)
-
-    if args.delete_tmp == True:
-        shutil.rmtree(tmp_folder)
-        
-        
-if __name__ == "__main__":
-    run_pyLQT()

pylqt-1.0.dev1.dist-info/METADATA

@@ -1,62 +0,0 @@
-Metadata-Version: 2.4
-Name: pyLQT
-Version: 1.0.dev1
-Summary: Python version of LQT from Griffis et al. 2021 (https://www.sciencedirect.com/science/article/pii/S2213158221000838)
-Author-email: Antonin HUDRY <hudryantoninm@gmail.com>
-License-Expression: GPL-3.0-only
-Keywords: LQT,lqt,lesion,brain,
-Classifier: Development Status :: 2 - Pre-Alpha
-Classifier: Intended Audience :: Science/Research
-Classifier: Intended Audience :: Healthcare Industry
-Classifier: Programming Language :: Python
-Requires-Python: >=3.10
-Description-Content-Type: text/markdown
-License-File: LICENSE
-Requires-Dist: numpy<2.0
-Requires-Dist: nibabel>5.0
-Requires-Dist: pandas>3.0
-Requires-Dist: scipy>1.15
-Requires-Dist: requests>2.0
-Requires-Dist: dipy>1.10
-Requires-Dist: templateflow
-Dynamic: license-file
-
-# pyLQT
-Python library to revamp the work on LQT, since it is no longer actively maintained. 
-
-# install
-
-Make sure to update your python environment with the `requirements.txt`: `pip install -r requirements.txt`
-
-# Running pyLQT
-
-The main pyLQT pipeline can be run using the `full_pipeline.py` file.
-
-### Usage
-
-```bash
-python full_pipeline.py \
-    --lesion_cfg <path/to/lesion_config> \
-    --atlas_cfg <path/to/atlas_config> \
-    --tracts_cfg <path/to/tracts_config> \
-    --out_folder <path/to/output> \
-    --spared_threshold 
-```
-
-| Argument             | Required | Default | Description                                                                                                                |
-| -------------------- | -------- | ------- | -------------------------------------------------------------------------------------------------------------------------- |
-| `--lesion_cfg`    | Yes      | —       | Path to the lesion configuration file (`.json`).                                   |
-| `--atlas_cfg`        | Yes      | —       | Path to the atlas configuration file (`.json`).                                                                                      |
-| `--tracts_cfg`           | No       | `None`  | Path to the tracts configuration file (`.json`).                        |
-| `--out_folder`       | Yes      | —       | Path to the output folder.                                                                                                 |
-| `--spared_threshold` | No       | `0.5`   | Threshold between 0 and 1 used to define spared parcel-to-parcel connectivity.                                             |
-| `--delete_tmp`       | No       | `True`  | Whether to delete the temporary files generated during processing.                                                         |
-
-### Configuration files
-
-The CLI version of pyLQT relies on 3 configuration files, one for the lesions, one for the atlas and one for the tracts:
-
-- Lesions: "cohort_name": name of the cohort of lesion used mainly for the naming of subsequent files; "lesion_folder": relative folder (should be next the config file) containing all the lesion niftis. The folder is scanned integrally so that all niftis in the folder and subfolders are processed; "space": world space in which the lesion are registered, can be a nifti directly given as relative path (with extensions) or a standard space name that will be fetched via templateflow (i.e. MNI152NLin2009cAsym).
-- Atlas: "name": name of the atlas used mainly for naming subsequent files; "atlas": relative path to nifti file corresponding to atlas; "atlas_labels": relative path to `.csv` file containing the labels of the atlas. Format should contain columns "Parcel" and "Label", corresponding to the parcel number inside the nifti and the name of the parcel respectively; "space": world space in which the lesion are registered, can be a nifti directly given as relative path (with extensions) or a standard space name that will be fetched via templateflow (i.e. MNI152NLin2009cAsym).
-
-Note: niftis can be given as raw or gzipped niftis (`.nii` or `.nii.gz`)

pylqt-1.0.dev1.dist-info/RECORD

@@ -1,21 +0,0 @@
-pyLQT_cli.py,sha256=FlTLPGcoq-xwPeGwVUfe_W_dg8yMpYV1X1jPBf-EefA,7937
-lesion_metrics/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-lesion_metrics/gray_matter/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-lesion_metrics/gray_matter/metrics.py,sha256=RzUS9glbin9Xn5ZldB7fGg2XiPS1P21e_7l02ZErk-I,2617
-lesion_metrics/gray_matter/utils.py,sha256=czxyAESP0neVCue424L0hon-5Tdbsl0pN0zJHWG5HiQ,840
-lesion_metrics/network/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-lesion_metrics/network/metrics.py,sha256=Ka3G8_dIq-8H3b8mvRR96QM7p_gR2yFLKUxQOLYMWM4,8281
-lesion_metrics/network/utils.py,sha256=fcakfKSMortNIE3TfJP8zwanzBiT1McsYnHcTgwbxE4,3283
-lesion_metrics/utils/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-lesion_metrics/utils/configs.py,sha256=FFur04Sh48IlxUTmz4kRkyLKWffi7bvzyw8B2GclpRU,4961
-lesion_metrics/utils/downloads.py,sha256=yu-rcuIDcgsDa7gxoLGyIewsk2DF2WGgCYe5bCOso0M,8031
-lesion_metrics/utils/registration.py,sha256=rB4vEMBrlshofOxaKrSqyE3Wvgq2_o3y5T3dCS_9sM0,3174
-lesion_metrics/white_matter/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-lesion_metrics/white_matter/metrics.py,sha256=sujuHry19XpVWFNo52Ntq7h3JkWEcgsZ1-rNFr7E940,6155
-lesion_metrics/white_matter/utils.py,sha256=oo492yzxsPACMvHEyYhoYVqgVRrv8mpqP7FOf5OU_Ms,7087
-pylqt-1.0.dev1.dist-info/licenses/LICENSE,sha256=IwGE9guuL-ryRPEKi6wFPI_zOhg7zDZbTYuHbSt_SAk,35823
-pylqt-1.0.dev1.dist-info/METADATA,sha256=gMSuT2w08jR1QSTH_C9TRsLVsJBEnCh98jSDH4_eCGs,3990
-pylqt-1.0.dev1.dist-info/WHEEL,sha256=YVMoNqKzERt-wjUZwJ33xBGAwnFl-4cqbYkTtWa4itE,91
-pylqt-1.0.dev1.dist-info/entry_points.txt,sha256=1mhG9fCXeDnL9VEBbsIN34GU90953e85325pG_wo6l8,46
-pylqt-1.0.dev1.dist-info/top_level.txt,sha256=KT_xuwXScSnV5QoZAXnCzeTOq23xQ_V9ESEh4jrfxBM,25
-pylqt-1.0.dev1.dist-info/RECORD,,

pylqt-1.0.dev1.dist-info/entry_points.txt

@@ -1,2 +0,0 @@
-[console_scripts]
-pyLQT = pyLQT_cli:run_pyLQT

pylqt-1.0.dev1.dist-info/top_level.txt

@@ -1,2 +0,0 @@
-lesion_metrics
-pyLQT_cli