EIDORS: Electrical Impedance Tomography and Diffuse Optical Tomography Reconstruction Software |
EIDORS
(mirror) Main Documentation Tutorials − Image Reconst − Data Structures − Applications − FEM Modelling − GREIT − Old tutorials − Workshop Download Contrib Data GREIT Browse Docs Browse SVN News Mailing list (archive) FAQ Developer
|
Reconstruction and waveforms of lungsDataThis tutorial shows interesting lung images and waveforms which can be reconstructed from EIT data. This example uses data from a 10-day old spontaneously breathing neonate. In a prone position with the head to the left, most of the tidal volume flows to the right lung.Data are available Here. Data were recorded from a 10-day old spontaneously breathing neonate lying in the prone position with the head turned to the left, as documented in: S. Heinrich, H. Schiffmann, A. Frerichs, A. Klockgether-Radke, I. Frerichs, Body and head position effects on regional lung ventilation in infants: an electrical impedance tomography study. Intensive Care Med., 32:1392-1398, 2006. Download data%zn = 'if-neonate-spontaneous.zip'; %urlwrite(['http://eidors3d.sf.net/data_contrib/if-neonate-spontaneous/',zn],zn); zipfilecontents = unzip('../../data_contrib/if-neonate-spontaneous/if-neonate-spontaneous.zip'); Select reconstruction algorithmHere we use the GREIT algorithm on a neonate_16el geometry.fmdl = mdl_normalize(mk_library_model('neonate_16el'),1); [fmdl.stimulation,fmdl.meas_select] = mk_stim_patterns(16,1,'{ad}','{ad}'); opt.imgsz = [64 64]; opt.noise_figure = 0.5; imdl = mk_GREIT_model(mk_image(fmdl,1), 0.25, [], opt); Show reconstructed imagevv= eidors_readdata('P04P-1016.get'); vi=vv(:,45); vh=vv(:,61); imr = inv_solve(imdl,vh,vi); clf; axes('position',[0.05,0.5,0.25,0.45]); imr.calc_colours = struct('ref_level',0,'greylev',0.2,'backgnd',[1,1,1]); show_slices(imr); Indicate points to show waveformsyposns = [45 20 50]; xposns = [50 40 27]; ofs= [0,22,15]; % Show positions on image hold on; for i = 1:length(xposns) plot(xposns(i),yposns(i),'s','LineWidth',10); end; hold off; Show waveforms in the lungs and heart% Show plots imgs = calc_slices(inv_solve(imdl, vh, vv)); axes('position',[0.32,0.6,0.63,0.25]); taxis = (0:size(imgs,3)-1)/13; % frame rate = 13 hold all for i = 1:length(xposns); plot(taxis,ofs(i)+squeeze(imgs(yposns(i),xposns(i),:)),'LineWidth',2); end hold off set(gca,'yticklabel',[]); xlim([0 16]); print_convert neonate_intro04a.jpg Figure: Left: Recontructed tidal images, and Right: pixel waveforms in the (red) heart, right lung (yellow) and left lung (blue) |
Last Modified: $Date: 2018-03-20 17:30:05 -0400 (Tue, 20 Mar 2018) $ by $Author: aadler $