README for data in Chapter 2 in Developing parametric finite element models of the residual limb – prosthetic socket system towards clinical application


Author: J. W. Steer
Date: 19 June 2019

This contains instructions for opening the .pkl file containing the data behind the figures. The .pkl file is strutured using dictionaries

Opening the data file
---------------------

The data file is encoded using the pickle package in python 3.6. The data is all stored within numpy arrays. The authors recommend downloading the anaconda distribution of python. 

"""
import pickle
import numpy
with open('data.pkl', 'rb') as f:
	Data = pickle.load(f)
"""

Figure 4.4
--------

Socket Rectification maps

Participant = ['A', 'B']
Socket = ['TSB', 'PTB', 'KBM']


Example
Data['Figure 4.4']['A']['TSB']
Returns a Nx4 numpy array of the individuals corresponding to vertices and rectification 

Figure 4.5
--------

Socket Rectification PDE

Participant = ['A', 'B']
Socket = ['TSB', 'PTB', 'KBM']


Example
Data['Figure 4.5']['A']['TSB']
Returns a numpy array of the PDE for each socket

Figure 4.6
--------

Difference in soft tissue depth

Participant = ['A', 'B']
Socket = ['TSB', 'PTB', 'KBM']


Example
Data['Figure 4.6']['A']['TSB']
Returns a Nx4 numpy array of the vertices and soft tissue depth


Figure 4.7
--------

Load magnitude for sockets

Participant = ['A', 'B']
Socket = ['TSB', 'PTB', 'KBM']


Example
Data['Figure 4.7']['A']['TSB']
Returns a Nx2 numpy array of the load and angle off-vertical 



Figure 4.8/9
--------

Experimental pressure data with each session

Participant = ['A', 'B']
Socket = ['TSB', 'PTB', 'KBM']


Example
Data['Figure 4.8']['A']['TSB']
Returns a Nx7 numpy array of the session number, 'A-P', 'A-M', 'A-D', 'L-P', 'L-M' and 'A-D' 

Figure 4.10
--------

Example pressure distribution 

Participant = ['A', 'B']
Socket = ['TSB', 'PTB', 'KBM']


Example
Data['Figure 4.10']['A']['TSB']
Returns a Nx2 numpy array of anterior and lateral pressure sensor


Figure 4.11
--------

Predicted pressure distribution

Participant = ['A', 'B']
Socket = ['TSB', 'PTB', 'KBM']


Example
Data['Figure 4.11']['A']['TSB']
Returns a Nx4 numpy array of the vertices and pressure values 


Figure 4.12
--------

Predicted pressure at key locations 

Participant = ['A', 'B']
Socket = ['TSB', 'PTB', 'KBM']


Example
Data['Figure 4.12']['A']['TSB']
Returns a Nx4 numpy array at the Residuum Tip, Fibula Head, Tibial Tuberosity and Patella Tendon