This repositry contains the Matlab m-files to run our LD 3D TFM algorithm. The first part of the package includes the FIDVC algorithm, which is how the 3D displacement fields are calculated. The second part of the code package includes converting those displacement fields into 3D surface tractions as described in Toyjanova, Hannen et al., Soft Matter, 2014.
Most of the Matlab file structures are very similar to our LD 3D TFM algorithm. The current code structure is set up to compute viscoelastic tractions based on the constitutive model depicted in our Soft Matter paper (Toyjanova, Hannen et al., Soft Matter, 2014), but the user can simply tune the code to fit the specifics of their mechanical viscoelastic model. Feel free to contact us with questions or with help on implementing a user-specific model into our 3D VTFM code structure.
To run you need a compatible C compiler. Please see (http://www.mathworks.com/support/compilers/R2015a/index.html)
- To check if the 3D image stack have the required speckle pattern and intensity values for correlation please use our DVC simulator.
- The 3D image stack need to be saved in a 3 dimensional matrix (intensity values are stored at x, y and z position) in vol.mat* files.
- We recommend that the input image stack at each dimension should have at least 1.5 times of the subset size as the number of pixels. The default subset size is 128x128x64, so we recommend that the minimum input volume size should be 192x192x96.
- The size of the input image stack should be divisible by 0.5 times the size of the subset.
- Make sure that the main files and the supplemental m files (from file exchange) are added to the path in Matlab.
- Download and save the example volume data in the example folder.
- Run the exampleRunFile.m file
FIDVC in the 3D-VTFM requires a 3D stack to be read in, which depending on the volume size can require a large amount of RAM in Matlab.
- Main files
- calculateNormals.m
- calculateSurfaceUi.m
- findSurface.m
- fun3VDTFM.m
- removeOutliers.m
- Supplement m files from the MATLAB file exchange:
- gridfit.m
- inpaint_nans.m
- inpaint_nans3.m
- Example Run files
- exampleRunFile.m
- example volume data (vol00.mat, vol01.mat).
What are the requirements for the input 3D image stack?
Please refer to input 3D Image Stack Requirements.
If used please cite: Toyjanova J., Hannen, E., Bar-Kochba E., Darling, E.M., Henann, D.L., and Franck, C., (2014) 3D Viscoelastic Traction Force Microscopy. Soft Matter doi: 10.1039/c4sm01271b
@article{toyjanova20143d,
title={3D Viscoelastic traction force microscopy},
author={Toyjanova, Jennet and Hannen, Erin and Bar-Kochba, Eyal and Darling, Eric M and Henann, David L and Franck, Christian},
journal={Soft matter},
volume={10},
number={40},
pages={8095--8106},
year={2014},
publisher={Royal Society of Chemistry}
}
For questions, please first refer to FAQ and Questions/Issues. Add a new question if similar issue hasn't been reported. We shall help you at the earliest. The author's contact information can be found at Franck Lab.