python库读取freesurfer的皮质分割解剖学统计
freesurfer-stats的Python项目详细描述
python库读取freesurfer的皮质分割解剖学统计 主题/stats/[rl]h.aparc.*.stats
自由曲面https://surfer.nmr.mgh.harvard.edu/
安装
pip3 install --user freesurfer-stats
发布遵循语义版本控制<;https://semver.org/>;方案。
用法
>>>fromfreesurfer_statsimportCorticalParcellationStats>>>stats=CorticalParcellationStats.read('tests/subjects/fabian/stats/lh.aparc.DKTatlas.stats')>>>stats.headers['subjectname']'fabian'>>>stats.headers['CreationTime'].isoformat()'2019-05-09T21:05:54+00:00'>>>stats.headers['cvs_version']'Id: mris_anatomical_stats.c,v 1.79 2016/03/14 15:15:34 greve Exp'>>>stats.headers['cmdline'][:64]'mris_anatomical_stats -th3 -mgz -cortex ../label/lh.cortex.label'>>>stats.hemisphere'left'>>>stats.whole_brain_measurements['estimated_total_intracranial_volume_mm^3']01.670487e+06Name:estimated_total_intracranial_volume_mm^3,dtype:float64>>>stats.whole_brain_measurements['white_surface_total_area_mm^2']098553Name:white_surface_total_area_mm^2,dtype:int64>>>stats.structural_measurements[['structure_name','surface_area_mm^2',...'gray_matter_volume_mm^3']].head()structure_namesurface_area_mm^2gray_matter_volume_mm^30caudalanteriorcingulate147242581caudalmiddlefrontal303982392cuneus259767223entorhinal49923794fusiform30799064
加载多个统计文件
>>>importglob,pandas>>>fromfreesurfer_statsimportCorticalParcellationStats>>>defload_whole_brain_measurements(stats_path)->pandas.DataFrame:...stats=CorticalParcellationStats.read(stats_path)...stats.whole_brain_measurements['subject']=stats.headers['subjectname']...stats.whole_brain_measurements['source_basename']=os.path.basename(stats_path)...stats.whole_brain_measurements['hemisphere']=stats.hemisphere...returnstats.whole_brain_measurements...>>>whole_brain_measurements=pandas.concat(...map(load_whole_brain_measurements,glob.glob('tests/subjects/fabian/stats/*h.aparc*.stats')),...sort=False)>>>whole_brain_measurements.reset_index(drop=True,inplace=True)>>>whole_brain_measurements[['subject','source_basename','hemisphere',...'white_surface_total_area_mm^2','pial_surface_total_area_mm^2']]subjectsource_basenamehemispherewhite_surface_total_area_mm^2pial_surface_total_area_mm^20fabianlh.aparc.DKTatlas.statsleft98553.0NaN1fabianrh.aparc.statsright99468.9NaN2fabianrh.aparc.a2009s.statsright99494.9NaN3fabianrh.aparc.DKTatlas.statsright99494.9NaN4fabianlh.aparc.statsleft98536.5NaN5fabianlh.aparc.pial.statsleftNaN118601.06fabianrh.aparc.pial.statsrightNaN121260.07fabianlh.aparc.a2009s.statsleft98553.0NaN>>>