Comments (1)
Hi Stefania,
For completeness sake I provide my answers also here on Github:
To answer your questions:
-
when you add models to the folder .mdt/.../components/user/composite_models it will be picked up next time you start MDT. If you are using the GUI you have to restart it, if you are using iPython you have to reload the MDT library, if you are using command line it should just pick it up. Typically, if it doesn't pick it up there is something wrong with the model.
-
MDT was originally based on the premise of multi-compartment models for modelling diffusion MRI signal. As such, the basic components in MDT are "compartments" which, when defined, can be freely mixed into (multi-compartment) models. The compartments were initially items such as "Tensor", "Ball", "Cylinder", "Stick" etc., following the definitions from Basser, Zhang and other early dMRI modellers. Later scientific work experimented with T1 and T2 decay in the dMRI models. Adding this to MDT was easy by defining T1 and T2 decays as additional "compartments", and mixing these in with the multi-compartment models. This works, although at the expense of naming clarity, since T1 and T2 decays are not compartments as defined in multi-compartment dMRI models.
To give an example, the ActiveAx model is defined as such:
model_expression = '''
S0 * ((Weight(w_csf) * Ball) +
(Weight(w_ic) * CylinderGPD) +
(Weight(w_ec) * Zeppelin))
'''
Meaning we have the Ball, CylinderGPD and Zeppelin compartments (these are well-defined dMRI compartments) mixed together using weights and S0 for the unweighted signal.
If you want to add T1 decay to that, you could change it to:
model_expression = '''
S0 * ExpT1DecTM_simple * ((Weight(w_csf) * Ball) +
(Weight(w_ic) * CylinderGPD) +
(Weight(w_ec) * Zeppelin))
'''
And voila, you have T1 decay in the ActiveAx model.
Have a look at the file ".mdt//components/standard/compartment_models/relaxometry.py for an overview of available decay models.
The decay constants are adding as an additional free parameter. These can be fixed using provided methods. Please see the documentation for that.
- The parameter transform is used by MDT to scale parameters just prior to optimization, and scale back again for presentation. To the users, MDT uses SI units and you should provide the T1 in seconds. To the optimization routine, MDT scales it such that the ranges of all parameters are approximately equivalent. Optimization routines struggle if one parameter is defined between 0.01 and 0.10 and another between 3000 and 6000 (S0 for example). Scaling helps to put everything in perspective.
Best,
Robbert
from mdt.
Related Issues (20)
- Test fails HOT 2
- interpretation NODDI output HOT 6
- Array problems HOT 1
- RuntimeError: clEnqueueNDRangeKernel failed: OUT_OF_RESOURCES HOT 2
- option --patience HOT 2
- Access tu more GPU nodes HOT 3
- Threading when using CPUs HOT 4
- Difference between CHARMED models HOT 1
- Confusing error message "Double operator in model listing." HOT 1
- Open CL errors HOT 14
- Noise estimation HOT 8
- Singularity runs, but gives semi-random results HOT 8
- Initialization HOT 2
- Value Error cannot find model 'NODDI' HOT 3
- AttributeError: module 'pyopencl.cltypes' has no attribute '4' HOT 5
- fixing parameters in fit process
- Numpy errors HOT 1
- Issue with invalid array size HOT 2
- mdt-list-models is empty, unable to select expected models for fit
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from mdt.