import os
    import distutils.util
    from pathlib import Path
    import numpy as np
    import matplotlib.pyplot as plt
    from true_vs_apparent.common.plot_utils import init_graphing, make_interactive, mean_sd_plot, style_axes
    from true_vs_apparent.common.database import create_db, BiplaneViconSubject, pre_fetch
    from true_vs_apparent.common.analysis_utils import prepare_db, extract_sub_rot, read_ludewig_data
    from true_vs_apparent.common.json_utils import get_params
    from true_vs_apparent.common.arg_parser import mod_arg_parser
    import logging
    from logging.config import fileConfig

    config_dir = Path(
        mod_arg_parser("Compare our dataset against Ludewig's", __package__,
                       __file__))
    params = get_params(config_dir / 'parameters.json')

    # if the database variables have been retained (i.e. we are re-running the script) then skip retrieving data from
    # disk
    if not bool(distutils.util.strtobool(os.getenv('VARS_RETAINED', 'False'))):
        # ready db
        db = create_db(params.biplane_vicon_db_dir,
                       BiplaneViconSubject,
                       include_anthro=True)
        if params.excluded_trials:
            db = db[~db['Trial_Name'].str.contains('|'.join(params.
                                                            excluded_trials))]
        db['Trial'].apply(pre_fetch)

    # logging
    import matplotlib.ticker as ticker
    import spm1d
    from true_vs_apparent.common.plot_utils import (init_graphing,
                                                    mean_sd_plot, extract_sig,
                                                    style_axes, sig_filter,
                                                    output_spm_p)
    from true_vs_apparent.common.database import create_db, BiplaneViconSubject, pre_fetch
    from true_vs_apparent.common.analysis_utils import prepare_db, extract_sub_rot_norm, sub_rot_at_max_elev
    from true_vs_apparent.common.json_utils import get_params
    from true_vs_apparent.common.arg_parser import mod_arg_parser
    import logging
    from logging.config import fileConfig

    config_dir = Path(
        mod_arg_parser(
            "Compare yx'y'' and zx'y'' GH axial rotation against true axial rotation ",
            __package__, __file__))
    params = get_params(config_dir / 'parameters.json')

    if not bool(distutils.util.strtobool(os.getenv('VARS_RETAINED', 'False'))):
        # ready db
        db = create_db(params.biplane_vicon_db_dir,
                       BiplaneViconSubject,
                       include_anthro=True)
        db['age_group'] = db['Age'].map(lambda age: '<35'
                                        if age < 40 else '>45')
        exc_trials = [
            "O45_003_CA_t01", "O45_003_SA_t02", "O45_003_FE_t02",
            "U35_010_FE_t01"
        ]
        db = db[~db['Trial_Name'].str.contains('|'.join(exc_trials))]
示例#3
0
    import distutils.util
    from pathlib import Path
    import numpy as np
    import matplotlib.pyplot as plt
    import spm1d
    from true_vs_apparent.common.plot_utils import init_graphing, make_interactive, mean_sd_plot, spm_plot, style_axes
    from true_vs_apparent.common.database import create_db, BiplaneViconSubject, pre_fetch
    from true_vs_apparent.common.analysis_utils import prepare_db, extract_sub_rot, extract_sub_rot_norm
    from true_vs_apparent.common.json_utils import get_params
    from true_vs_apparent.common.arg_parser import mod_arg_parser
    import logging
    from logging.config import fileConfig

    config_dir = Path(
        mod_arg_parser(
            "Compare yx'y'', zx'y'' and true GH axial rotation against zero",
            __package__, __file__))
    params = get_params(config_dir / 'parameters.json')

    if not bool(distutils.util.strtobool(os.getenv('VARS_RETAINED', 'False'))):
        # ready db
        db = create_db(params.biplane_vicon_db_dir,
                       BiplaneViconSubject,
                       include_anthro=True)
        db['age_group'] = db['Age'].map(lambda age: '<35'
                                        if age < 40 else '>45')
        if params.excluded_trials:
            db = db[~db['Trial_Name'].str.contains('|'.join(params.
                                                            excluded_trials))]
        db['Trial'].apply(pre_fetch)
    import os
    import distutils.util
    from pathlib import Path
    import numpy as np
    import matplotlib.pyplot as plt
    from true_vs_apparent.common.plot_utils import (init_graphing, mean_sd_plot, style_axes)
    from true_vs_apparent.common.database import create_db, BiplaneViconSubject, pre_fetch
    from true_vs_apparent.common.analysis_utils import prepare_db, extract_sub_rot_norm
    from true_vs_apparent.common.json_utils import get_params
    from true_vs_apparent.common.arg_parser import mod_arg_parser
    import logging
    from logging.config import fileConfig

    config_dir = Path(mod_arg_parser("Compute true axial rotation based on difference scapula anatomical "
                                     "coordinate systems",
                                     __package__, __file__))
    params = get_params(config_dir / 'parameters.json')

    if not bool(distutils.util.strtobool(os.getenv('VARS_RETAINED', 'False'))):
        # ready db
        db = create_db(params.biplane_vicon_db_dir, BiplaneViconSubject, include_anthro=True)
        db['age_group'] = db['Age'].map(lambda age: '<35' if age < 40 else '>45')
        exc_trials = ["O45_003_CA_t01", "O45_003_SA_t02", "O45_003_FE_t02", "U35_010_FE_t01"]
        db = db[~db['Trial_Name'].str.contains('|'.join(exc_trials))]
        db['Trial'].apply(pre_fetch)

    # relevant parameters
    output_path = Path(params.output_dir)

    # logging
示例#5
0
    import matplotlib.pyplot as plt
    import matplotlib.ticker as ticker
    import spm1d
    from true_vs_apparent.common.plot_utils import (init_graphing,
                                                    mean_sd_plot, extract_sig,
                                                    style_axes, sig_filter)
    from true_vs_apparent.common.database import create_db, BiplaneViconSubject, pre_fetch
    from true_vs_apparent.common.analysis_utils import prepare_db, extract_sub_rot_norm, sub_rot_at_max_elev
    from true_vs_apparent.common.json_utils import get_params
    from true_vs_apparent.common.arg_parser import mod_arg_parser
    import logging
    from logging.config import fileConfig

    config_dir = Path(
        mod_arg_parser(
            "Compare GH yx'y'' and zx'y'' GH axial rotation against true axial rotation "
            "normalized by 25 HT elev", __package__, __file__))
    params = get_params(config_dir / 'parameters.json')

    if not bool(distutils.util.strtobool(os.getenv('VARS_RETAINED', 'False'))):
        # ready db
        db = create_db(params.biplane_vicon_db_dir,
                       BiplaneViconSubject,
                       include_anthro=True)
        db['age_group'] = db['Age'].map(lambda age: '<35'
                                        if age < 40 else '>45')
        exc_trials = [
            "O45_003_CA_t01", "O45_003_SA_t02", "O45_003_FE_t02",
            "U35_010_FE_t01"
        ]
        db = db[~db['Trial_Name'].str.contains('|'.join(exc_trials))]
示例#6
0
    import os
    import distutils.util
    from pathlib import Path
    import numpy as np
    import matplotlib.pyplot as plt
    import spm1d
    from true_vs_apparent.common.plot_utils import init_graphing, make_interactive, mean_sd_plot, style_axes
    from true_vs_apparent.common.database import create_db, BiplaneViconSubject, pre_fetch
    from true_vs_apparent.common.analysis_utils import prepare_db, extract_sub_rot_norm, sub_rot_at_max_elev
    from true_vs_apparent.common.json_utils import get_params
    from true_vs_apparent.common.arg_parser import mod_arg_parser
    import logging
    from logging.config import fileConfig

    config_dir = Path(mod_arg_parser("Visualize how accounting for plane of elevation affects the ISB decomposition",
                                     __package__, __file__))
    params = get_params(config_dir / 'parameters.json')

    if not bool(distutils.util.strtobool(os.getenv('VARS_RETAINED', 'False'))):
        # ready db
        db = create_db(params.biplane_vicon_db_dir, BiplaneViconSubject, include_anthro=True)
        db['age_group'] = db['Age'].map(lambda age: '<35' if age < 40 else '>45')
        if params.excluded_trials:
            db = db[~db['Trial_Name'].str.contains('|'.join(params.excluded_trials))]
        db['Trial'].apply(pre_fetch)

    # relevant parameters
    output_path = Path(params.output_dir)

    # logging
    fileConfig(config_dir / 'logging.ini', disable_existing_loggers=False)
                                                    make_interactive,
                                                    mean_sd_plot, extract_sig,
                                                    style_axes, sig_filter,
                                                    output_spm_p)
    from true_vs_apparent.common.database import create_db, BiplaneViconSubject, pre_fetch
    from true_vs_apparent.common.analysis_utils import (
        prepare_db, extract_sub_rot_norm, sub_rot_at_max_elev,
        subj_name_to_number, extract_sub_rot, sub_rot_at_min_elev)
    from true_vs_apparent.common.json_utils import get_params
    from true_vs_apparent.common.arg_parser import mod_arg_parser
    import logging
    from logging.config import fileConfig

    config_dir = Path(
        mod_arg_parser(
            "Compare GH yx'y'', xz'y'' and true axial rotation between planes of elevation",
            __package__, __file__))
    params = get_params(config_dir / 'parameters.json')

    if not bool(distutils.util.strtobool(os.getenv('VARS_RETAINED', 'False'))):
        # ready db
        db = create_db(params.biplane_vicon_db_dir,
                       BiplaneViconSubject,
                       include_anthro=True)
        db['age_group'] = db['Age'].map(lambda age: '<35'
                                        if age < 40 else '>45')
        exc_trials = [
            "O45_003_CA_t01", "O45_003_SA_t02", "O45_003_FE_t02",
            "U35_010_FE_t01"
        ]
        db = db[~db['Trial_Name'].str.contains('|'.join(exc_trials))]
    if __package__ is None:
        print('Use -m option to run this library module as a script.')

    import os
    import distutils.util
    from pathlib import Path
    from matplotlib.backends.backend_pdf import PdfPages
    from true_vs_apparent.common.analysis_utils import get_trajs
    from true_vs_apparent.common.up_down import analyze_up_down
    from true_vs_apparent.common.database import create_db, BiplaneViconSubject, pre_fetch
    from true_vs_apparent.common.json_utils import get_params
    from true_vs_apparent.common.arg_parser import mod_arg_parser
    from logging.config import fileConfig

    config_dir = Path(
        mod_arg_parser('Identify up down sections per trial', __package__,
                       __file__))
    params = get_params(config_dir / 'parameters.json')
    if not bool(distutils.util.strtobool(os.getenv('VARS_RETAINED', 'False'))):
        # ready db
        db = create_db(params.biplane_vicon_db_dir, BiplaneViconSubject)
        db = db[~db['Trial_Name'].str.contains('|'.join(params.excluded_trials)
                                               )]
        db['Trial'].apply(pre_fetch)

    # relevant parameters
    output_path = Path(params.output_dir)

    # logging
    fileConfig(config_dir / 'logging.ini', disable_existing_loggers=False)
    log = logging.getLogger(params.logger_name)
示例#9
0
    import distutils.util
    from pathlib import Path
    import numpy as np
    import matplotlib.pyplot as plt
    import spm1d
    from true_vs_apparent.common.plot_utils import init_graphing, make_interactive, mean_sd_plot, spm_plot, style_axes
    from true_vs_apparent.common.database import create_db, BiplaneViconSubject, pre_fetch
    from true_vs_apparent.common.analysis_utils import prepare_db, extract_sub_rot, extract_sub_rot_norm
    from true_vs_apparent.common.json_utils import get_params
    from true_vs_apparent.common.arg_parser import mod_arg_parser
    import logging
    from logging.config import fileConfig

    config_dir = Path(
        mod_arg_parser(
            "Compare yx'y'', zx'y'' and true GH axial rotation for differences between "
            "genders", __package__, __file__))
    params = get_params(config_dir / 'parameters.json')

    if not bool(distutils.util.strtobool(os.getenv('VARS_RETAINED', 'False'))):
        # ready db
        db = create_db(params.biplane_vicon_db_dir,
                       BiplaneViconSubject,
                       include_anthro=True)
        db['age_group'] = db['Age'].map(lambda age: '<35'
                                        if age < 40 else '>45')
        if params.excluded_trials:
            db = db[~db['Trial_Name'].str.contains('|'.join(params.
                                                            excluded_trials))]
        db['Trial'].apply(pre_fetch)
示例#10
0
    import os
    import distutils.util
    from pathlib import Path
    import numpy as np
    import matplotlib.pyplot as plt
    import spm1d
    from true_vs_apparent.common.plot_utils import init_graphing, make_interactive, mean_sd_plot, spm_plot, style_axes
    from true_vs_apparent.common.database import create_db, BiplaneViconSubject, pre_fetch
    from true_vs_apparent.common.analysis_utils import prepare_db, extract_sub_rot, extract_sub_rot_norm
    from true_vs_apparent.common.json_utils import get_params
    from true_vs_apparent.common.arg_parser import mod_arg_parser
    import logging
    from logging.config import fileConfig

    config_dir = Path(
        mod_arg_parser('Compare using GC vs AC when analyzing trials',
                       __package__, __file__))
    params = get_params(config_dir / 'parameters.json')

    # if the database variables have been retained (i.e. we are re-running the script) then skip retrieving data from
    # disk
    if not bool(distutils.util.strtobool(os.getenv('VARS_RETAINED', 'False'))):
        # ready db
        db = create_db(params.biplane_vicon_db_dir,
                       BiplaneViconSubject,
                       include_anthro=True)
        if params.excluded_trials:
            db = db[~db['Trial_Name'].str.contains('|'.join(params.
                                                            excluded_trials))]
        db['Trial'].apply(pre_fetch)

    # relevant parameters
if __name__ == '__main__':
    if __package__ is None:
        print('Use -m option to run this library module as a script.')

    import os
    import distutils.util
    from pathlib import Path
    from true_vs_apparent.common.database import create_db, BiplaneViconSubject, pre_fetch
    from true_vs_apparent.common.analysis_utils import prepare_db
    from true_vs_apparent.common.analysis_er_utils import ready_er_db
    from true_vs_apparent.common.json_utils import get_params
    from true_vs_apparent.common.arg_parser import mod_arg_parser
    import logging
    from logging.config import fileConfig

    config_dir = Path(mod_arg_parser('Compare different methods of computing true GH and HT axial rotation',
                                     __package__, __file__))
    params = get_params(config_dir / 'parameters.json')

    if not bool(distutils.util.strtobool(os.getenv('VARS_RETAINED', 'False'))):
        # ready db
        db = create_db(params.biplane_vicon_db_dir, BiplaneViconSubject, include_anthro=True)
        db['age_group'] = db['Age'].map(lambda age: '<35' if age < 40 else '>45')
        if params.excluded_trials:
            db = db[~db['Trial_Name'].str.contains('|'.join(params.excluded_trials))]
        db['Trial'].apply(pre_fetch)

    # relevant parameters
    output_path = Path(params.output_dir)

    # logging
    fileConfig(config_dir / 'logging.ini', disable_existing_loggers=False)
示例#12
0
    import matplotlib.pyplot as plt
    import spm1d
    from true_vs_apparent.common.plot_utils import init_graphing, make_interactive, mean_sd_plot, spm_plot, style_axes
    from true_vs_apparent.common.database import create_db, BiplaneViconSubject, pre_fetch
    from true_vs_apparent.common.analysis_utils import (prepare_db,
                                                        extract_sub_rot,
                                                        extract_sub_rot_norm,
                                                        subj_name_to_number)
    from true_vs_apparent.common.json_utils import get_params
    from true_vs_apparent.common.arg_parser import mod_arg_parser
    import logging
    from logging.config import fileConfig

    config_dir = Path(
        mod_arg_parser(
            "Compare yx'y'' and zx'y'' GH axial rotation between planes of elevation",
            __package__, __file__))
    params = get_params(config_dir / 'parameters.json')

    if not bool(distutils.util.strtobool(os.getenv('VARS_RETAINED', 'False'))):
        # ready db
        db = create_db(params.biplane_vicon_db_dir,
                       BiplaneViconSubject,
                       include_anthro=True)
        db['age_group'] = db['Age'].map(lambda age: '<35'
                                        if age < 40 else '>45')
        if params.excluded_trials:
            db = db[~db['Trial_Name'].str.contains('|'.join(params.
                                                            excluded_trials))]
        db['Trial'].apply(pre_fetch)
    import distutils.util
    from pathlib import Path
    import numpy as np
    import matplotlib.pyplot as plt
    import spm1d
    from true_vs_apparent.common.plot_utils import init_graphing, make_interactive, mean_sd_plot, spm_plot, style_axes
    from true_vs_apparent.common.database import create_db, BiplaneViconSubject, pre_fetch
    from true_vs_apparent.common.analysis_utils import prepare_db, extract_sub_rot, extract_sub_rot_norm
    from true_vs_apparent.common.json_utils import get_params
    from true_vs_apparent.common.arg_parser import mod_arg_parser
    import logging
    from logging.config import fileConfig

    config_dir = Path(
        mod_arg_parser(
            "Compare yx'y'', zx'y'' and true GH axial rotation for differences between arm "
            "raising and lowering", __package__, __file__))
    params = get_params(config_dir / 'parameters.json')

    if not bool(distutils.util.strtobool(os.getenv('VARS_RETAINED', 'False'))):
        # ready db
        db = create_db(params.biplane_vicon_db_dir,
                       BiplaneViconSubject,
                       include_anthro=True)
        db['age_group'] = db['Age'].map(lambda age: '<35'
                                        if age < 40 else '>45')
        if params.excluded_trials:
            db = db[~db['Trial_Name'].str.contains('|'.join(params.
                                                            excluded_trials))]
        db['Trial'].apply(pre_fetch)
    from pathlib import Path
    import numpy as np
    import matplotlib.pyplot as plt
    from matplotlib.backends.backend_pdf import PdfPages
    from true_vs_apparent.common import plot_utils
    from true_vs_apparent.common.database import create_db, BiplaneViconSubject, pre_fetch
    from true_vs_apparent.common.json_utils import get_params
    from true_vs_apparent.common.arg_parser import mod_arg_parser
    from true_vs_apparent.common.analysis_er_utils import ready_er_db, plot_directives
    import matplotlib.ticker as plticker
    import logging
    from logging.config import fileConfig

    config_dir = Path(
        mod_arg_parser(
            "Create an overview of axial rotation for external rotation trials",
            __package__, __file__))
    params = get_params(config_dir / 'parameters.json')

    if not bool(distutils.util.strtobool(os.getenv('VARS_RETAINED', 'False'))):
        # ready db
        db = create_db(params.biplane_vicon_db_dir,
                       BiplaneViconSubject,
                       include_anthro=True)
        db['age_group'] = db['Age'].map(lambda age: '<35'
                                        if age < 40 else '>45')
        if params.excluded_trials:
            db = db[~db['Trial_Name'].str.contains('|'.join(params.
                                                            excluded_trials))]
        db['Trial'].apply(pre_fetch)
if __name__ == '__main__':
    if __package__ is None:
        print('Use -m option to run this library module as a script.')

    import os
    import distutils.util
    from pathlib import Path
    from true_vs_apparent.common import plot_utils
    from true_vs_apparent.common.database import create_db, BiplaneViconSubject, pre_fetch
    from true_vs_apparent.common.analysis_utils import prepare_db
    from true_vs_apparent.common.json_utils import get_params
    from true_vs_apparent.common.arg_parser import mod_arg_parser
    from logging.config import fileConfig

    config_dir = Path(
        mod_arg_parser('Overview of elevations trials', __package__, __file__))
    params = get_params(config_dir / 'parameters.json')
    if not bool(distutils.util.strtobool(os.getenv('VARS_RETAINED', 'False'))):
        # ready db
        db = create_db(params.biplane_vicon_db_dir,
                       BiplaneViconSubject,
                       include_anthro=True)
        db['age_group'] = db['Age'].map(lambda age: '<35'
                                        if age < 40 else '>45')
        if params.excluded_trials:
            db = db[~db['Trial_Name'].str.contains('|'.join(params.
                                                            excluded_trials))]
        db['Trial'].apply(pre_fetch)

    # relevant parameters
    output_path = Path(params.output_dir)
    import matplotlib.pyplot as plt
    import spm1d
    from true_vs_apparent.common.plot_utils import init_graphing, make_interactive, mean_sd_plot, spm_plot, style_axes
    from true_vs_apparent.common.database import create_db, BiplaneViconSubject, pre_fetch
    from true_vs_apparent.common.analysis_utils import (prepare_db,
                                                        extract_sub_rot,
                                                        extract_sub_rot_norm,
                                                        subj_name_to_number)
    from true_vs_apparent.common.json_utils import get_params
    from true_vs_apparent.common.arg_parser import mod_arg_parser
    import logging
    from logging.config import fileConfig

    config_dir = Path(
        mod_arg_parser(
            "Compare yx'y'' and zx'y'' GH, HT, ST (HT-GH) axial rotation between planes of "
            "elevation.", __package__, __file__))
    params = get_params(config_dir / 'parameters.json')

    if not bool(distutils.util.strtobool(os.getenv('VARS_RETAINED', 'False'))):
        # ready db
        db = create_db(params.biplane_vicon_db_dir,
                       BiplaneViconSubject,
                       include_anthro=True)
        db['age_group'] = db['Age'].map(lambda age: '<35'
                                        if age < 40 else '>45')
        if params.excluded_trials:
            db = db[~db['Trial_Name'].str.contains('|'.join(params.
                                                            excluded_trials))]
        db['Trial'].apply(pre_fetch)
    import matplotlib.ticker as ticker
    import spm1d
    from true_vs_apparent.common.plot_utils import (init_graphing,
                                                    mean_sd_plot, extract_sig,
                                                    style_axes, sig_filter,
                                                    output_spm_p)
    from true_vs_apparent.common.database import create_db, BiplaneViconSubject, pre_fetch
    from true_vs_apparent.common.analysis_utils import prepare_db, extract_sub_rot_norm, sub_rot_at_max_elev
    from true_vs_apparent.common.json_utils import get_params
    from true_vs_apparent.common.arg_parser import mod_arg_parser
    import logging
    from logging.config import fileConfig

    config_dir = Path(
        mod_arg_parser(
            "Compare GH yx'y'' and zx'y'' axial rotation against true axial rotation "
            "acromial angle", __package__, __file__))
    params = get_params(config_dir / 'parameters.json')

    if not bool(distutils.util.strtobool(os.getenv('VARS_RETAINED', 'False'))):
        # ready db
        db = create_db(params.biplane_vicon_db_dir,
                       BiplaneViconSubject,
                       include_anthro=True)
        db['age_group'] = db['Age'].map(lambda age: '<35'
                                        if age < 40 else '>45')
        exc_trials = [
            "O45_003_CA_t01", "O45_003_SA_t02", "O45_003_FE_t02",
            "U35_010_FE_t01", "O45_002_CA_t02", "O45_002_SA_t03",
            "O45_001_ERaR_t01", "O45_002_ERaR_t02"
        ]

if __name__ == '__main__':
    if __package__ is None:
        print('Use -m option to run this library module as a script.')

    import numpy as np
    from pathlib import Path
    from true_vs_apparent.common.database import create_db, BiplaneViconSubject
    from true_vs_apparent.common.analysis_utils import prepare_db
    from true_vs_apparent.common.json_utils import get_params
    from true_vs_apparent.common.arg_parser import mod_arg_parser
    from logging.config import fileConfig
    import logging

    config_dir = Path(mod_arg_parser('Apparent - true axial rotation = spherical area', __package__, __file__))
    params = get_params(config_dir / 'parameters.json')
    db = create_db(params.biplane_vicon_db_dir, BiplaneViconSubject)

    # logging
    fileConfig(config_dir / 'logging.ini', disable_existing_loggers=False)
    log = logging.getLogger(params.logger_name)

    # prepare database
    db_elev = db.loc[db['Trial_Name'].str.contains('_CA_|_SA_|_FE_')].copy()
    prepare_db(db_elev, params.torso_def, params.scap_lateral, params.dtheta_fine, params.dtheta_coarse,
               [params.min_elev, params.max_elev], should_fill=False, should_clean=False)

    if 'isb' in params.decomp_method:
        fnc_path = isb_path
        fnc_axial_orient_diff = isb_add_axial_diff
    import matplotlib.ticker as plticker
    from true_vs_apparent.common import plot_utils
    from true_vs_apparent.common.database import create_db, BiplaneViconSubject, pre_fetch
    from true_vs_apparent.common.json_utils import get_params
    from true_vs_apparent.common.arg_parser import mod_arg_parser
    from true_vs_apparent.common.analysis_er_utils import ready_er_db
    from true_vs_apparent.common.plot_utils import (mean_sd_plot,
                                                    make_interactive,
                                                    style_axes, extract_sig,
                                                    sig_filter, output_spm_p)
    import logging
    from logging.config import fileConfig

    config_dir = Path(
        mod_arg_parser(
            "Compare GH yx'y'', yx'y'' normalized, and xz'y'' against true axial rotation for "
            "external rotation trials", __package__, __file__))
    params = get_params(config_dir / 'parameters.json')

    if not bool(distutils.util.strtobool(os.getenv('VARS_RETAINED', 'False'))):
        # ready db
        db = create_db(params.biplane_vicon_db_dir,
                       BiplaneViconSubject,
                       include_anthro=True)
        db['age_group'] = db['Age'].map(lambda age: '<35'
                                        if age < 40 else '>45')
        exc_trials = [
            "O45_003_CA_t01", "O45_003_SA_t02", "O45_003_FE_t02",
            "U35_010_FE_t01"
        ]
        db = db[~db['Trial_Name'].str.contains('|'.join(exc_trials))]
示例#20
0
        print('Use -m option to run this library module as a script.')

    import numpy as np
    import matplotlib.pyplot as plt
    from pathlib import Path
    from true_vs_apparent.common.database import create_db, BiplaneViconSubject
    from true_vs_apparent.common.analysis_utils import prepare_db
    from true_vs_apparent.analysis.up_down_identify import extract_up_down_min_max
    from true_vs_apparent.common import plot_utils
    from true_vs_apparent.common.json_utils import get_params
    from true_vs_apparent.common.arg_parser import mod_arg_parser
    from logging.config import fileConfig
    import logging

    config_dir = Path(
        mod_arg_parser('Check U35_002_SA_t01 filling', __package__, __file__))
    params = get_params(config_dir / 'parameters.json')
    db = create_db(params.biplane_vicon_db_dir, BiplaneViconSubject)
    db = db.loc[['U35_002_SA_t01']]

    # logging
    fileConfig(config_dir / 'logging.ini', disable_existing_loggers=False)
    log = logging.getLogger(params.logger_name)

    # compute min and max ht elevation for each subject
    prepare_db(db,
               params.torso_def,
               params.scap_lateral,
               params.dtheta_fine,
               params.dtheta_coarse, [params.min_elev, params.max_elev],
               should_clean=False)
示例#21
0
    import distutils.util
    from pathlib import Path
    import matplotlib.pyplot as plt
    import matplotlib.ticker as ticker
    from true_vs_apparent.common.plot_utils import (init_graphing, style_axes,
                                                    make_interactive)
    from true_vs_apparent.common.database import create_db, BiplaneViconSubject, pre_fetch
    from true_vs_apparent.common.analysis_utils import prepare_db
    from true_vs_apparent.common.json_utils import get_params
    from true_vs_apparent.common.arg_parser import mod_arg_parser
    from true_vs_apparent.analysis.overview_elev_trials import summary_plotter, ind_plotter
    import logging
    from logging.config import fileConfig

    config_dir = Path(
        mod_arg_parser("Plot HT yx'y'' axial rotation for individual subjects",
                       __package__, __file__))
    params = get_params(config_dir / 'parameters.json')

    if not bool(distutils.util.strtobool(os.getenv('VARS_RETAINED', 'False'))):
        # ready db
        db = create_db(params.biplane_vicon_db_dir,
                       BiplaneViconSubject,
                       include_anthro=True)
        db['age_group'] = db['Age'].map(lambda age: '<35'
                                        if age < 40 else '>45')
        if params.excluded_trials:
            db = db[~db['Trial_Name'].str.contains('|'.join(params.
                                                            excluded_trials))]
        db['Trial'].apply(pre_fetch)

    # relevant parameters
示例#22
0
if __name__ == '__main__':
    if __package__ is None:
        print('Use -m option to run this library module as a script.')

    import os
    import distutils.util
    from pathlib import Path
    from true_vs_apparent.common.database import create_db, BiplaneViconSubject, pre_fetch
    from true_vs_apparent.common.analysis_utils import prepare_db
    from true_vs_apparent.common.json_utils import get_params
    from true_vs_apparent.common.arg_parser import mod_arg_parser
    import logging
    from logging.config import fileConfig

    config_dir = Path(
        mod_arg_parser('Determine the drift due to linearization', __package__,
                       __file__))
    params = get_params(config_dir / 'parameters.json')

    if not bool(distutils.util.strtobool(os.getenv('VARS_RETAINED', 'False'))):
        # ready db
        db = create_db(params.biplane_vicon_db_dir,
                       BiplaneViconSubject,
                       include_anthro=True)
        db['age_group'] = db['Age'].map(lambda age: '<35'
                                        if age < 40 else '>45')
        if params.excluded_trials:
            db = db[~db['Trial_Name'].str.contains('|'.join(params.
                                                            excluded_trials))]
        db['Trial'].apply(pre_fetch)

    # relevant parameters