ip_data_filename = dir_path + "ip_data.xyz"
########################################################
# Load Data and Topography
# ------------------------
#
# Here we load the observed data and topography.
#
#
topo_xyz = np.loadtxt(str(topo_filename))
dc_data = read_dcip_xyz(
dc_data_filename,
"volt",
data_header="V/A",
uncertainties_header="UNCERT",
is_surface_data=False,
)
ip_data = read_dcip_xyz(
ip_data_filename,
"apparent_chargeability",
data_header="APP_CHG",
uncertainties_header="UNCERT",
is_surface_data=False,
)
##########################################################
# Plot Observed DC Data in Pseudosection
# --------------------------------------
# containing 3D DC resistivity data.
#
# Load 3D topography
topo_xyz = np.loadtxt(str(topo_filename))
# Load 3D data. Here, the data are loaded from an XYZ formatted data file.
# The user must supply the proper headers for the function to identify the
# correct column. Using the 'dict_headers' keyword argument, we can load and
# organize additional columns in the data file as a dictionary.
data_3d, out_dict = read_dcip_xyz(
data_filename,
"volt",
a_headers=["XA", "YA", "ZA"],
b_headers=["XB", "YB", "ZB"],
m_headers=["XM", "YM", "ZM"],
n_headers=["XN", "YN", "ZN"],
data_header="V/A",
uncertainties_header="UNCERT",
dict_headers=["LINEID"],
)
#######################################################################
# Plot 3D Pseudosection
# ---------------------
#
# Here we demonstrate how 3D DC resistivity data can be represented on a 3D
# pseudosection plot. To use this utility, you must have Python's *plotly*
# package. Here, we represent the data as apparent conductivities.
#
dc_data_filename = dir_path + "dc_data.xyz"
########################################################
# Load Data and Topography
# ------------------------
#
# Here we load the observed data and topography.
#
#
topo_xyz = np.loadtxt(str(topo_filename))
dc_data = read_dcip_xyz(
dc_data_filename,
"volt",
data_header="V/A",
uncertainties_header="UNCERT",
is_surface_data=False,
)
##########################################################
# Plot Observed Data in Pseudosection
# -----------------------------------
#
# Here we plot the observed DC and IP data in 3D pseudosections.
# To use this utility, you must have Python's *plotly* package.
# Here, we represent the DC data as apparent conductivities
# and the IP data as apparent chargeabilities.
#
topo_filename = dir_path + "topo_xyz.txt"
dc_data_filename = dir_path + "dc_data.xyz"
########################################################
# Load Data and Topography
# ------------------------
#
# Here we load the observed data and topography.
#
#
topo_xyz = np.loadtxt(str(topo_filename))
dc_data = read_dcip_xyz(dc_data_filename,
'volt',
data_header='V/A',
uncertainties_header='UNCERT',
is_surface_data=False)
##########################################################
# Plot Observed Data in Pseudosection
# -----------------------------------
#
# Here we plot the observed DC and IP data in 3D pseudosections.
# To use this utility, you must have Python's *plotly* package.
# Here, we represent the DC data as apparent conductivities
# and the IP data as apparent chargeabilities.
#
# Convert predicted data to apparent conductivities
apparent_conductivity = 1 / apparent_resistivity_from_voltage(
dc_data_filename = dir_path + "dc_data.xyz"
ip_data_filename = dir_path + "ip_data.xyz"
########################################################
# Load Data and Topography
# ------------------------
#
# Here we load the observed data and topography.
#
#
topo_xyz = np.loadtxt(str(topo_filename))
dc_data = read_dcip_xyz(dc_data_filename,
'volt',
data_header='V/A',
uncertainties_header='UNCERT',
is_surface_data=False)
ip_data = read_dcip_xyz(ip_data_filename,
'apparent_chargeability',
data_header='APP_CHG',
uncertainties_header='UNCERT',
is_surface_data=False)
##########################################################
# Plot Observed DC Data in Pseudosection
# --------------------------------------
#
# Here we plot the observed DC data in 3D pseudosection.
# To use this utility, you must have Python's *plotly* package.