'MSoilR100AnoR1Z-13.75W15H2.5L5S0Noise0',
'MSoilR1000AnoR1Z-13.75W15H2.5L5S0Noise0'
]
x_axis = 'y'
for fi in filenames:
print(fi)
x_raw, y_raw, z_raw, U_raw, maxdepth, shape_raw, p1, p2, SimName, ano_prop = MALM.load_MALM_sens3d(
filename='./loadmalm/' + fi + '.pkl')
# pEXP.plot_field(x_raw, y_raw,)
shape = (200, 200)
xp, yp, U = gridder.interp(x_raw, y_raw, U_raw, shape)
parameters = para.set_par(shape=shape, max_elevation=abs(maxdepth))
interp = True
scaled = parameters[0]
SI = parameters[1]
zp, qorder, nlay = parameters[2:5]
minAlt_ridge, maxAlt_ridge = parameters[5:7]
#%%
# ridges analysis parameters
nlay = 25
max_elevation = 30
minAlt_ridge = max_elevation * 0.05
maxAlt_ridge = max_elevation * 0.65
interp = True
smooth = False
data_struct = grav.load_grav_fatiando(
name='loadgrav/za3000_zb3500_l500_ofs0_dens1200')
xp, yp, zp, U = data_struct['xyzg']
shape = data_struct['shape']
model = data_struct['model']
dens = data_struct['density']
# scaled, SI, zp, qorder, nlay, minAlt_ridge, maxAlt_ridge = para.set_par(shape=shape,max_elevation=max_elevation)
x1, x2, y1, y2, z1, z2 = np.array(model[0].get_bounds())
p1 = [min(yp), 0]
p2 = [max(yp), 0]
max_elevation = z2 * 1.2
scaled, SI, zp, qorder, nlay, minAlt_ridge, maxAlt_ridge = para.set_par(
shape=shape, max_elevation=max_elevation)
interp = True
qorder = 0
#%%
# Plot the data
pEXP.plot_line(xp, yp, U, p1, p2, interp=interp)
#%%
# Pad the edges of grids (if necessary)
xp, yp, U, shape = dEXP.pad_edges(xp, yp, U, shape, pad_type=0) # reflexion=5
p1 = [min(yp), 0]
p2 = [max(yp), 0]
x_axis = 'y'
shape, max_elevation = dataset[4:6]
dict_data = dataset[7]
dict_data['AnoBool']
xA = (dict_data['HZ'][0][0] + dict_data['HZ'][0][1]) / 2
x1 = dict_data['HZ'][0][2]
y1 = dict_data['HZ'][0][0]
y2 = dict_data['HZ'][0][1]
z1 = dict_data['HZ'][1]
z2 = z1 - dict_data['HZ'][2]
p = dataset[6] # line points
# set imaging pseudo-inversion parameters
parameters = para.set_par(shape=shape, max_elevation=max_elevation)
scaled = parameters[0]
SI = parameters[1]
zp, qorder, nlay = parameters[2:5]
minAlt_ridge, maxAlt_ridge = parameters[5:7]
max_elevation = 30
# nlay = 50
# xp, yp, zp = coord_xyz_int
xp, yp, zp = coord_xyz
# len(xp)
Uini = Uload[0] # U_raw, Ucor, U_int, Ucor_int
p1, p2 = p