EIDORS: Electrical Impedance Tomography and Diffuse Optical Tomography Reconstruction Software

Hosted by

Simulate a moving 2D target on a fixed mesh

Step 1: Create fine mesh (6076 elems)

% Simulate movement $Id: simulate_move1_01.m 1535 2008-07-26 15:36:27Z aadler $
imdm=mk_common_model('e2d4c',16);
smdl= imdm.fwd_model; % simulation model

Step 2: Simulate movement

% Simulate movement $Id: simulate_move1_02.m 2770 2011-07-14 16:54:10Z bgrychtol $
[vh,vi,xyr_pt]= simulate_2d_movement( 20, smdl,[.75,.05]);

% Show model and simulated targets
show_fem(smdl);
theta= linspace(0,2*pi,50); xr= cos(theta); yr= sin(theta);
hold on;
for i=1:length(xyr_pt)
    hh= plot(xyr_pt(3,i)*xr+ xyr_pt(1,i), ...
             xyr_pt(3,i)*yr+ xyr_pt(2,i));
    set(hh,'LineWidth',3,'Color',[0,0,1]);
    text(xyr_pt(1,i),xyr_pt(2,i),sprintf('%d',i), ...
        'HorizontalAlignment','center','FontSize',8, ...
        'Color',[0,0,1],'FontWeight','bold');
end
hold off;
axis([-1.1 1.1 -1.1 1.1]); axis equal
print_convert simulate_move1_02a.png '-density 100'

Figure: Simulation mesh and simulated target positions

Step 3: Reconstruct images

% Reconstruct $Id: simulate_move1_03.m 2770 2011-07-14 16:54:10Z bgrychtol $

imdl= mk_common_model('c2c2',16');
img= inv_solve(imdl,vh,vi);

subplot(121)
show_fem(img)
axis image

subplot(122)
show_slices(img)

print_convert simulate_move1_03a.png '-density 125'

Figure: Reconstructed images Left FEM of first image Right Reconstructed image sequence