Creating synthetic models

Create a grid, e.g., using the helper function of the crtomo.grid module:

import crtomo.grid as CRGrid
grid = CRGrid.crt_grid.create_surface_grid(
    nr_electrodes=30, spacing=1
)

Now create a tomo manager instance to manage parameter sets

import crtomo.tdManager as CRman
man = CRman.tdMan(grid=grid)

Then, create two new parameter sets, one for the magnitudes, one for the phase values:

pid_mag, pid_pha = man.add_homogeneous_model(magnitude=100, phase=-5)

Note

If you are only interested in the magnitude model, leave the phase at zero

Note

You can alternatively load a rho.dat file, as used by CRTomo using

pid_mag, pid_pha = man.load_rho_file('rho.dat')

Modify the one parmeter set at a time using one of the modifying functions:

man.parman.modify_area(
    pid_mag,
    xmin=1, xmax=5,
    zmin=-3, zmax=-2,
    value=10
)

man.parman.modify_area(
    pid_pha,
    xmin=1, xmax=5,
    zmin=-3, zmax=-2,
    value=-30
)

Parameter sets can be plotted using the plot manager:

result_mag = man.plot.plot_elements_to_ax(
    pid_mag,
    plot_colorbar=True,
)
fig = result_mag[0]
fig.savefig('model_mag.png', dpi=300)

result_pha = man.plot.plot_elements_to_ax(
    pid_pha,
    plot_colorbar=True,
)
fig = result_pha[0]
fig.savefig('model_pha.png', dpi=300)

(Source code)

(png, hires.png, pdf)

(png, hires.png, pdf)