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
|
Using Netgen to create elliptical EIDORS ModelsEIDORS can use Netgen to create sophisticated 2D and 3D modelsHere are some examples of the varity of models which can be generated using the function: ng_mk_ellip_models.
Simple 3D ellipse. Major, minor axes = [1.5, 0.8]. No electrodesfmdl= ng_mk_ellip_models([1,1.5,0.8],[0],[]); show_fem(fmdl); Simple 2D cylinder. Axes = [1.5,0.8]. Refined to 0.1fmdl= ng_mk_ellip_models([0,1.5,0.8,0.1],[],[]); show_fem(fmdl); 3D cylinder. Axes = [1.5,0.8]. 2 planes of 8 elecs with radius 0.1fmdl= ng_mk_ellip_models([1,1.2,0.8],[8,0.3,0.7],[0.1]); show_fem(fmdl); 3D cylinder. Axes= [1.3,1] = 1. 7 rect elecs with no refinementfmdl= ng_mk_ellip_models([3,1.3],[7,1],[0.2,0.3]); show_fem(fmdl); 2D cylinder. Axes = [1.2,0.8]. 11 rect elecs with refinement. Rotated 45 degreesfmdl= ng_mk_ellip_models([0,1.2,0.8],[11],[0.2,0,0.05]); th = 45* [2*pi/360]; fmdl.nodes = fmdl.nodes*[cos(th),sin(th);-sin(th),cos(th)]; show_fem(fmdl); 2D cylinder. elecs at 0, 90 and 120 degreesfmdl= ng_mk_ellip_models([0,1.2,0.8],[0;90;120],[0.2,0,0.03]); show_fem(fmdl); 3D cylinder. Various elecs at 0, 30, 60, 90 in planesel_pos = [0,0.5;30,1;60,1.5;90,2.0]; el_sz = [0.2,0,0.1;0.1,0,0.05;0.2,0.2,0.02;0.2,0.4,0.5]; fmdl= ng_mk_ellip_models([3,0.8,1.2],el_pos,el_sz); show_fem(fmdl); Simple 3D cylinder with a ballextra={'ball','solid ball = sphere(0.5,0.5,1;0.4);'}; fmdl= ng_mk_ellip_models([2,1.2,0.8],[8,1],[0.1],extra); img= mk_image(fmdl, 1); img.elem_data(fmdl.mat_idx{2}) = 2; show_fem(img); 3D cylinder with a two separate ballsb1 = 'solid ball1= sphere(0.5, 0.5,1;0.2);'; b2 = 'solid ball2= sphere(0.5,-0.5,1;0.2);'; extra = {'ball1','ball2',[b1,b2]}; [fmdl,mat_idx]= ng_mk_ellip_models([2,1.2,0.8],[8,1],[0.1],extra); img= mk_image(fmdl, 1); img.elem_data(mat_idx{2}) = 2; img.elem_data(mat_idx{3}) = 0.5; show_fem(img); |
Last Modified: $Date: 2017-02-28 13:12:08 -0500 (Tue, 28 Feb 2017) $ by $Author: aadler $