To receive notifications about scheduled maintenance, please subscribe to the mailing-list gitlab-operations@sympa.ethz.ch. You can subscribe to the mailing-list at https://sympa.ethz.ch

Commit ffe7438e authored by Feliks Kiszkurno's avatar Feliks Kiszkurno
Browse files

changed create data

parent 824c8ed5
......@@ -25,7 +25,7 @@ create_new_data = True # set to True if you need to reassign the classes
invert_existing_data = False # invert existing measurements
create_new_data_only = True # set to False in order to run ML classifications
reassign_classes = False; class_type = 'norm'
param_path = os.path.abspath(os.path.join(os.getcwd()) + '/' + 'TestDefinitions/sen_study.csv')
param_path = os.path.abspath(os.path.join(os.getcwd()) + '/' + 'TestDefinitions/sen_study2.csv')
test_definitions.init(path=param_path)
# Load existing data instead of creating new one.
......
......@@ -27,7 +27,7 @@ def plot_sen_corr(y_pred, y_answer, sen, clf_name, test_name, batch_name, *, tra
match_plot = data_sort['MATCH'].to_numpy()
fig, ax = plt.subplots(1, 1)
ax.plot(sen_plot, match_plot)
ax.scatter(sen_plot, match_plot)
fig.suptitle('Sensitivity vs correctness of classification')
ax.set_xlabel('Sensitivity [normalized]')
ax.set_ylabel('Correctness of classification [T/F]')
......
......@@ -36,10 +36,10 @@ def plot_and_save_pg(test_name, plot_title, ert_manager, input_model, result_ful
#cb[0].locator = tick_locator
#cb[0].update_ticks()
limits = pg.utils.interperc(ert_manager.inv.model, trimval=7.0)
limits = pg.utils.interperc(ert_manager.inv.model, trimval=10.0)
print(limits)
im1 = pg.show(ert_manager.paraDomain, ert_manager.inv.model, c_min=limits[0], c_max=limits[1], showMesh=True, ax=ax[1],
im1 = pg.show(ert_manager.paraDomain, ert_manager.inv.model, cMin=limits[0], cMax=limits[1], showMesh=True, ax=ax[1],
label='Resistivity \u03A9 *m')
ax[1].set_title('Result')
ax[1] = slopestabilitytools.set_labels(ax[1])
......
......@@ -20,10 +20,10 @@ import settings
def create_data(test_name, test_config, max_depth, *, lambda_param=20, z_weight=0.6):
#world_boundary_v = [-200, 0] # [NW edge] relatively to the middle
#world_boundary_h = [200, -100] # [SE edge]
world_boundary_v = [-1000, 0] # [right, left border] relatively to the middle
world_boundary_h = [1000, -100] # [top, bottom border]
world_boundary_v = [-200, 0] # [NW edge] relatively to the middle
world_boundary_h = [200, -100] # [SE edge]
#world_boundary_v = [-1000, 0] # [right, left border] relatively to the middle
#world_boundary_h = [1000, -100] # [top, bottom border]
test_results = {}
......@@ -43,8 +43,8 @@ def create_data(test_name, test_config, max_depth, *, lambda_param=20, z_weight=
slopestabilitytools.save_plot(fig_geometry, test_name, '_1_geometry')
measurement_scheme = ert.createERTData(
#elecs=np.linspace(start=-32, stop=32, num=44),
elecs=np.linspace(start=-50, stop=50, num=100),
elecs=np.linspace(start=-32, stop=32, num=44),
#elecs=np.linspace(start=-50, stop=50, num=100),
schemeName='dd')
for electrode in measurement_scheme.sensors():
......@@ -78,7 +78,7 @@ def create_data(test_name, test_config, max_depth, *, lambda_param=20, z_weight=
k0 = pg.physics.ert.createGeometricFactors(data)
#inversion_Domain = mt.createMesh(mt.createRectangle([-35, 0], [35, -25], quality=34, area=1))
#inversion_mesh = pg.meshtools.appendTriangleBoundary(inversion_Domain, marker=0, xbound=30, ybound=30)
model_inverted = ert_manager.invert(data=data, lam=lambda_param, paraDX=0.25, paraMaxCellSize=2, zWeight=z_weight, paraDepth=15,
model_inverted = ert_manager.invert(data=data, lam=lambda_param, paraDX=0.25, paraMaxCellSize=2, zWeight=z_weight,# paraDepth=15,
quality=34, zPower=0.4)
result_full = ert_manager.inv.model
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment