def run(args): plot_scatter_files(args.data, mask_file=args.mask_file, masked_opacity=args.masked_opacity, volume=args.volume, thresholds=args.thresholds, limits=args.limits, x_jitter=args.x_jitter, y_jitter=args.y_jitter, uniq=args.unique_points, include_stats=args.stats) if args.output_img: pl.savefig(args.output_img) pl.show()
def test_plot_scatter_files_with_styles(kw): fns = glob(pjoin(pymvpa_dataroot, *('haxby2001/sub001/anatomy/lowres00*.nii.gz'.split('/')))) figs = plot_scatter_files(fns, **kw) # all of those produce 1 figure since we have only 2 files assert(figs) if kw.get('style', None) in {'full', 'pair2', 'upper_triang'}: # TODO: unify API assert(isinstance(figs[0], Figure)) else: assert(isinstance(figs, Figure)) # now let's give 3 files figs = plot_scatter_files(fns + fns[-1:], **kw) assert_equal(len(figs), 2) pl.close('all')
def test_plot_scatter_files_mask(): fns = glob(pjoin(pymvpa_dataroot, *('haxby2001/sub001/anatomy/lowres00*.nii.gz'.split('/')))) # figure out reasonable threshold data = _get_data(fns[0]) min_, max_ = np.min(data), np.max(data) plot_scatter_files(fns, mask_file=fns[0], mask_thresholds=(min_ + max_)/2., masked_opacity=0.5) # now fancier ones thrs = [min_ + (max_ - min_) * 0.2, min_ + (max_ - min_) * 0.8] plot_scatter_files(fns, mask_file=fns[0], mask_thresholds=thrs, masked_opacity=0.5) plot_scatter_files(fns, mask_file=fns[0], mask_thresholds=thrs[::-1], masked_opacity=0.5) pl.close('all')