EIDORS: Electrical Impedance Tomography and Diffuse Optical Tomography Reconstruction Software

Hosted by

Build simple 2D models with distmesh

The simplest way to build distmesh models is to use the function mk_common_models with the parameter like b2d0c.
−b is the mesh density (a is least, f is largest)
−2d asks for a 2D distmesh model
−0 is the level of electrode mesh refinement (0−4)
−c is a circular model

Create coarse mesh (416 elems)

% Simple distmesh model $Id: simple_2d_model01.m 2767 2011-07-14 16:32:47Z bgrychtol $

% Create small model with little mesh refinement
imdm=mk_common_model('a2d1c',16);
subplot(121)
show_fem(imdm.fwd_model);
axis on; axis equal; axis([-1.1 1.1 -1.1 1.1]);

% close in view
subplot(122)
show_fem(imdm.fwd_model);
axis equal; axis([-.1 .1 .85 1.05]);

print_convert simple_2d_model01a.png '-density 125'

Figure: Circular model generated by distmesh (416 elements) Left: complete model Right: magnified view of electrode (3 nodes per electrode)

Create fine mesh (3602 elems)

% Simple distmesh model $Id: simple_2d_model02.m 2767 2011-07-14 16:32:47Z bgrychtol $

% Create small model with little mesh refinement
imdm=mk_common_model('d2d4c',16);
subplot(121)
show_fem(imdm.fwd_model);
axis on; axis equal; axis([-1.1 1.1 -1.1 1.1]);

% close in view
subplot(122)
show_fem(imdm.fwd_model);
axis equal; axis([-.1 .1 .85 1.05]);

print_convert simple_2d_model02a.png '-density 125'

Figure: Circular model generated by distmesh (3602 elements) Left: complete model Left: magnified view of electrode (3 nodes per electrode)