cell_size=dem_cell_size,
elev_cell=elev_cell,
elev_max=1200.,
dem_rotation_angle=30)
#write new data file
d_obj.center_stations(new_model_fn)
d_obj.change_data_elevation(new_model_fn)
n_dfn = d_obj.write_data_file(save_path=sv_path,
fn_basename='um_data_ef03_tec.dat',
fill=False,
compute_error=False,
elevation=True)
# write covariance file
cov = modem.Covariance()
cov.smoothing_east = 0.4
cov.smoothing_north = 0.4
cov.smoothing_z = 0.4
cov.save_path = sv_path
cov.write_covariance_file(model_fn=new_model_fn)
# make paraview files if needed
m_obj.write_vtk_file()
d_obj.write_vtk_station_file()
# plot the model and make sure everything looks good
mm = modem.ModelManipulator(model_fn=new_model_fn,
data_fn=n_dfn,
depth_index=45)
m_obj.model_fn = None
m_obj.write_model_file(**{
'save_path': sv_path,
'model_fn_basename': 'mont_topography.rho'
})
#write new data file
d_obj.center_stations(m_obj.model_fn)
d_obj.change_data_elevation(m_obj.model_fn)
n_dfn = d_obj.write_data_file(fn_basename=r"mont_data_topo_c.dat",
compute_error=False,
fill=False,
elevation=True)
# write covariance file
cov = modem.Covariance()
cov.smoothing_east = 0.4
cov.smoothing_north = 0.4
cov.smoothing_z = 0.4
cov.save_path = sv_path
cov.write_covariance_file(model_fn=m_obj.model_fn)
#m_obj.write_vtk_file()
#d_obj.write_vtk_station_file()
mm = modem.ModelManipulator(model_fn=m_obj.model_fn,
data_fn=n_dfn,
depth_index=26,
xlimits=(-10, 10),
ylimits=(-10, 10))
new_model_fn = m_obj.add_topography_to_model(dem_fn,
model_fn,
model_center=model_center,
cell_size=dem_cell_size,
elev_cell=elev_cell,
elev_max=1260.)
#write new data file
n_dfn = d_obj.change_data_elevation(d_obj.data_fn,
new_model_fn,
new_data_fn=os.path.join(
sv_path, 'gz_data_topo.dat'))
# write covariance file
cov = modem.Covariance()
cov.smoothing_east = 0.4
cov.smoothing_north = 0.4
cov.smoothing_z = 0.4
cov.save_path = sv_path
cov.write_covariance_file(model_fn=new_model_fn)
# make paraview files if needed
m_obj.write_vtk_file()
d_obj.write_vtk_station_file()
# plot the model and make sure everything looks good
mm = modem.ModelManipulator(model_fn=new_model_fn,
data_fn=n_dfn,
depth_index=42,
xlimits=(-200, 10),
ylimits=(-200, 200))