import pyeit.mesh.plot as mplot from pyeit.eit.fem import Forward from pyeit.eit.pde import pdeprtni from pyeit.eit.utils import eit_scan_lines # build tetrahedron # 3D tetrahedron must have a bbox bbox = [[-1, -1, -1], [1, 1, 1]] # save calling convention as distmesh 2D ms, el_pos = mesh.create(h0=0.15, bbox=bbox) no2xy = ms['node'] el2no = ms['element'] # report the status of the 2D mesh quality.stats(no2xy, el2no) """ 1. FEM forward simulations """ # setup EIT scan conditions el_dist, step = 7, 1 ex_mat = eit_scan_lines(16, el_dist) # calculate simulated data fwd = Forward(ms, el_pos) # in python, index start from 0 ex_line = ex_mat[1].ravel() # change alpha anomaly = [{'x': 0.40, 'y': 0.40, 'z': 0.0, 'd': 0.30, 'alpha': 100.0}] ms_test = mesh.set_alpha(ms, anomaly=anomaly, background=1.0) tri_perm = ms_test['alpha']
import pyeit.mesh.plot as mplot from pyeit.eit.fem import Forward from pyeit.eit.interp2d import sim2pts from pyeit.eit.utils import eit_scan_lines # build tetrahedron # 3D tetrahedron must have a bbox bbox = [[-1, -1, -1], [1, 1, 1]] # save calling convention as distmesh 2D mesh_obj, el_pos = mesh.create(h0=0.15, bbox=bbox) pts = mesh_obj['node'] tri = mesh_obj['element'] # report the status of the 2D mesh quality.stats(pts, tri) """ 1. FEM forward simulations """ # setup EIT scan conditions el_dist, step = 7, 1 ex_mat = eit_scan_lines(16, el_dist) # calculate simulated data fwd = Forward(mesh_obj, el_pos) # in python, index start from 0 ex_line = ex_mat[1].ravel() # change alpha anomaly = [{'x': 0.40, 'y': 0.40, 'z': 0.0, 'd': 0.30, 'perm': 100.0}] mesh_new = mesh.set_perm(mesh_obj, anomaly=anomaly, background=1.0)