feat: add dataset-level registration QC summary with error decomposition

autoafids/workflow/Snakefile

@@ -382,6 +382,13 @@ rule all:
         )
         if config["LEAD_DBS_DIR"] or config["FMRIPREP_DIR"]
         else [],
+        regqc_summary=[
+            os.path.join(
+                root, "dataset", "regqc", "dataset_desc-reg_qc_summary.html"
+            )
+        ]
+        if config["LEAD_DBS_DIR"] or config["FMRIPREP_DIR"]
+        else [],
         afidspred=inputs[config["modality"]].expand(
             bids(
                 root=root,

autoafids/workflow/envs/regqc.yaml

@@ -2,7 +2,7 @@
 name: regqc
 channels: [khanlab,conda-forge, defaults]
 dependencies:
-  - python=3.9 
+  - python=3.10
   - plotly
   - simpleitk
   - scipy

autoafids/workflow/rules/regqc.smk

@@ -183,3 +183,44 @@ rule regqc:
         "../envs/regqc.yaml"
     script:
         "../scripts/regqc.py"
+
+
+rule regqc_summary:
+    input:
+        csvs=inputs[config["modality"]].expand(
+            bids(
+                root=root,
+                datatype="regqc",
+                desc="reg",
+                suffix="qc.csv",
+                **inputs[config["modality"]].wildcards,
+            )
+        ),
+        htmls=inputs[config["modality"]].expand(
+            bids(
+                root=root,
+                datatype="regqc",
+                desc="reg",
+                suffix="qc.html",
+                **inputs[config["modality"]].wildcards,
+            )
+        ),
+        fcsvs=inputs[config["modality"]].expand(
+            bids(
+                root=root,
+                datatype="regqc",
+                desc="reg",
+                suffix="afids.fcsv",
+                **inputs[config["modality"]].wildcards,
+            )
+        ),
+    output:
+        summary_html=os.path.join(
+            root, "dataset", "regqc", "dataset_desc-reg_qc_summary.html"
+        ),
+    params:
+        gt_fcsv=lambda wildcards: get_ref_paths()[1],
+    conda:
+        "../envs/regqc.yaml"
+    script:
+        "../scripts/regqc_summary.py"