import crtomo.grid as CRGrid
grid = CRGrid.crt_grid.create_surface_grid(
nr_electrodes=30, spacing=1
)
import crtomo.tdManager as CRman
man = CRman.tdMan(grid=grid)
pid_mag, pid_pha = man.add_homogeneous_model(magnitude=100, phase=-5)
man.parman.modify_area(pid_mag, 1, 5, -5, -2, value=10)
r = man.plotman.plot_elements_to_ax(
pid_mag,
title='Forward model',
plot_colorbar=True,
cblabel=r'$|\rho|~[\Omega m]$',
)
fig = r[0]
fig.savefig('model_mag.png', dpi=300)
r = man.plotman.plot_elements_to_ax(
pid_pha,
title='Forward model',
plot_colorbar=True,
cblabel=r'$\phi~[mrad]$',
)
fig = r[0]
fig.savefig('model_mag.png', dpi=300)
man.configs.gen_dipole_dipole(skipc=0)
rmag_rpha = man.measurements()
K = man.configs.compute_K_factors(spacing=1)
rhoa = rmag_rpha[:, 0] * K
cid_rho = man.configs.add_measurements(rhoa)
r = man.configs.plot_pseudosection_type2(
cid_rho,
cblabel=r'$\rho_a~[\Omega m]$',
)
import numpy as np
maglog = np.log10(rmag_rpha[:, 0])
cid = man.configs.add_measurements(maglog)
mid_mag, mid_pha = man.assignments['measurements']
r = man.configs.plot_pseudosection_type2(
mid_mag,
cblabel=r'$|Z|~[\Omega]$',
)
fig = r[0]
fig.savefig('pseudosection_mag.png', dpi=300)
r = man.configs.plot_pseudosection_type2(
maglog,
cblabel=r'$log_{10}(|Z|~[\Omega])$',
)
man.invert(output_directory=None, catch_output=False)
eps_data = man._read_eps_ctr('inversion/')
r = man.plot_eps_data(eps_data)
man.plot_eps_data_hist(eps_data)