From 8dc1cebbe3ba72d9791fc9b244cb726b70bf95c5 Mon Sep 17 00:00:00 2001
From: Paul McCarthy <pauldmccarthy@gmail.com>
Date: Fri, 5 Jul 2024 17:51:47 +0100
Subject: [PATCH] ENH: accessors for zfstats
---
fsl/data/featimage.py | 21 +++++++++++++++++++++
1 file changed, 21 insertions(+)
diff --git a/fsl/data/featimage.py b/fsl/data/featimage.py
index b67bcbb7..1ac5da33 100644
--- a/fsl/data/featimage.py
+++ b/fsl/data/featimage.py
@@ -72,9 +72,11 @@ class FEATImage(fslimage.Image):
if featanalysis.hasStats(featDir):
design = featanalysis.loadDesign( featDir, settings)
names, cons = featanalysis.loadContrasts(featDir)
+ ftests = featanalysis.loadFTests( featDir)
else:
design = None
names, cons = [], []
+ ftests = []
fslimage.Image.__init__(self, path, **kwargs)
@@ -83,12 +85,14 @@ class FEATImage(fslimage.Image):
self.__design = design
self.__contrastNames = names
self.__contrasts = cons
+ self.__ftests = ftests
self.__settings = settings
self.__residuals = None
self.__pes = [None] * self.numEVs()
self.__copes = [None] * self.numContrasts()
self.__zstats = [None] * self.numContrasts()
+ self.__zfstats = [None] * self.numFTests()
self.__clustMasks = [None] * self.numContrasts()
if 'name' not in kwargs:
@@ -102,6 +106,7 @@ class FEATImage(fslimage.Image):
self.__pes = None
self.__copes = None
self.__zstats = None
+ self.__zfstats = None
self.__clustMasks = None
@@ -191,6 +196,11 @@ class FEATImage(fslimage.Image):
return len(self.__contrasts)
+ def numFTests(self):
+ """Returns the number of f-tests in the analysis."""
+ return len(self.__ftests)
+
+
def contrastNames(self):
"""Returns a list containing the name of each contrast in the analysis.
"""
@@ -274,6 +284,17 @@ class FEATImage(fslimage.Image):
return self.__zstats[con]
+ def getZFStats(self, ftest):
+ """Returns the Z statistic image for the given f-test (0-indexed). """
+
+ if self.__zfstats[ftest] is None:
+ zfile = featanalysis.getZFStatFile(self.__featDir, ftest)
+ self.__zstats[con] = fslimage.Image(
+ zfile,
+ name=f'{self.__analysisName}: zfstat{con + 1}')
+ return self.__zfstats[ftest]
+
+
def getClusterMask(self, con):
"""Returns the cluster mask image for the given contrast (0-indexed).
"""
--
GitLab