else: data_min = data[data.nonzero()].min() return [data_min, data_max] print 'getting range' ranges = np.array([get_range(pfile) for pfile in pfiles]) dmin = ranges.min() dmax = ranges.max() print 'min=%.4f; max=%.4f' % (dmin, dmax) ### ### # Surface Viz for i, pfile in enumerate(pfiles): print pfile oprefix = path.join(odir, "C4_fstats_global_easythresh_surface_scan%i" % (i + 1)) sw = SurfWrap(name=factor, infile=pfile, cbar="red-yellow", outprefix=oprefix) sw.min = dmin sw.max = dmax sw.run(compilation="box") sw.montage(compilation="stick") ###
from os import path from surfwrap import SurfWrap obase = "/home2/data/Projects/CWAS/adhd200_rerun/cwas" comparisons = ["adhdc_vs_adhdi", "tdc_vs_adhdc", "tdc_vs_adhdi"] # GLOBAL! preproc = "global" for comparison in comparisons: print "group comparison: %s" % comparison # Input pfile distdir = path.join(obase, "%s_rois_random_k3200" % preproc) mdmrdir = path.join(distdir, "%s_gender+age+iq+mean_FD.mdmr" % comparison) pfile = path.join(mdmrdir, "log_pvals_%s.nii.gz" % comparison) # Output oprefix odir = path.join(mdmrdir, "images") if not path.exists(odir): os.mkdir(odir) oprefix = path.join(odir, path.basename(pfile).replace(".nii.gz", "")) # Surface Viz sw = SurfWrap(name=comparison, infile=pfile, cbar="red-yellow", outprefix=oprefix) sw.min = 1.3 sw.run() # import pysurfer # pysurfer = reload(pysurfer) # SurfWrap = pysurfer.SurfWrap
data_max = data.max() if data_max == 0: data_min = data_max else: data_min = data[data.nonzero()].min() return [data_min, data_max] print "getting range" ranges = np.array([get_range(pfile) for pfile in pfiles]) dmin = ranges[ranges != 0].min() dmax = ranges[ranges != 0].max() print "min=%.4f; max=%.4f" % (dmin, dmax) ### ### # Surface Viz for i, pfile in enumerate(pfiles): print pfile oprefix = path.join(odir, "C2_global_perms_surface_scan%i" % (i + 1)) sw = SurfWrap(name=factor, infile=pfile, cbar="red-yellow", outprefix=oprefix) sw.min = dmin sw.max = dmax sw.run(compilation="box") sw.montage(compilation="stick") ###
obase = "/home2/data/Projects/CWAS/adhd200_rerun/cwas" comparisons = ["adhdc_vs_adhdi", "tdc_vs_adhdc", "tdc_vs_adhdi"] # GLOBAL! preproc = "global" for comparison in comparisons: print "group comparison: %s" % comparison # Input pfile distdir = path.join(obase, "%s_rois_random_k3200" % preproc) mdmrdir = path.join(distdir, "%s_gender+age+iq+mean_FD.mdmr" % comparison) pfile = path.join(mdmrdir, "log_pvals_%s.nii.gz" % comparison) # Output oprefix odir = path.join(mdmrdir, "images") if not path.exists(odir): os.mkdir(odir) oprefix = path.join(odir, path.basename(pfile).replace(".nii.gz", "")) # Surface Viz sw = SurfWrap(name=comparison, infile=pfile, cbar="red-yellow", outprefix=oprefix) sw.min = 1.3 sw.run() # import pysurfer # pysurfer = reload(pysurfer) # SurfWrap = pysurfer.SurfWrap