EDI formats
EDI files have historically been the most common format for MT transfer functions. The metadata is quite flexible and can be verbose or scant in details about the transfer function. Needless to say there are many flavor of EDI files. The EDI module has been test with various flavors, if you find your’s is not read properly be sure to raise an issue with the “transfer function” label.
Here we will load a few different types of EDI files into TF. The files are common products from contracting companies.
Note: Spectra data can be readily converted. The spectra conversion is adopted from Anna Kelbert’s EMTF fortran code. This computes the inverse_signal_power
and residual_covariance
estimates for the errors to accurately estimate error upon rotation.
[48]:
from mt_metadata.transfer_functions import TF
from mt_metadata.transfer_functions.io.edi import EDI
Common output by Phoenix Geophysics
[3]:
from mt_metadata import TF_EDI_PHOENIX
[25]:
tf_object = TF(fn=TF_EDI_PHOENIX)
tf_object.read()
edi_object = EDI(fn=TF_EDI_PHOENIX)
[26]:
print(edi_object.Header)
print(edi_object.Info)
print(edi_object.Measurement)
print(edi_object.Data)
>HEAD
ACQBY=Phoenix
ACQDATE=2014-07-28T00:00:00+00:00
COORDINATE_SYSTEM=geographic
DATAID=14-IEB0537A
DATUM=WGS84
ELEV=158.000
EMPTY=1e+32
FILEBY=Phoenix
FILEDATE=2023-03-16
LAT=-22:49:25.40
LON=139:17:40.90
PROGDATE=2021-12-01
PROGNAME=mt_metadata
PROGVERS=0.2.2
STDVERS=SEG 1.0
UNITS=M
>INFO
CUTOFF: 0.00 COH: 35 % VAR: 25 %
Comp MTU box S/N Temp
Comp Chan# Sensor Azimuth
Ebat:12.3V Hbat:12.3V Rbat:11.9V
Elevation: 158 M. DECL: 0.000
Ex & Ey: MTU5A 2189 39 C
Ex1 1 100.0 M 0.0 DGtn
Ey1 2 100.0 M 90.0 DGtn
FILE: IEB0537A IEB0564M
Hx & Hy: MTU5A 2189 39 C
Hx1 3 COIL2318 0.0 DGtn
Hy1 4 COIL2319 90.0 DGtn
Hz1 5 COIL2320
Hz: MTU5A 2189 39 C
Lat 22:49:254S Long 139:17:409E
Lat 22:49.423 S Lng 139:17.681 E
MTU-Box Gains:E`s x 4 H`s x 4
MTU-DFT VERSION: TStoFT.38
MTU-RBS VERSION:R2012-0216-B22
MTU-Ref Serial Number: U-2779
Notch Filters set for 50 Hz.
PROCESSED FROM DFT TIME SERIES
RBS: 7 COH: 0.85 RHO VAR: 0.75
RHx2 6 COIL2485 0.0 DGtn
RHy2 7 COIL2487 90.0 DGtn
RUN INFORMATION
Reference Field: Remote H - Ref.
Reference Site: IEB0564M
Rx & Ry: MTU5A 2779 40 C
STATION 1
SYSTEM INFORMATION
Site Desc; BadR: 0 SatR: 54
Site Layout by:
Site Permitted by:
>=DEFINEMEAS
MAXCHAN=7
MAXRUN=999
MAXMEAS=7
REFLAT=-22.823722
REFLON=139.294694
REFELEV=158.0
REFTYPE=CART
UNITS=M
>HMEAS ID=5371.0537 CHTYPE=HX X=8.50 Y=8.50 Z=0.00 AZM=0.00 DIP=0.00 ACQCHAN=CH3
>HMEAS ID=5372.0537 CHTYPE=HY X=-8.50 Y=8.50 Z=0.00 AZM=90.00 DIP=0.00 ACQCHAN=CH4
>HMEAS ID=5373.0537 CHTYPE=HZ X=21.20 Y=-21.20 Z=0.00 AZM=0.00 DIP=0.00 ACQCHAN=CH5
>EMEAS ID=5374.0537 CHTYPE=EX X=-50.00 Y=-0.00 Z=0.00 X2=50.00 Y2=0.00 Z2=0.00 AZM=0.00 ACQCHAN=CH1
>EMEAS ID=5375.0537 CHTYPE=EY X=22.40 Y=-44.70 Z=0.00 X2=-22.40 Y2=44.70 Z2=0.00 AZM=116.62 ACQCHAN=CH2
>HMEAS ID=5376.0537 CHTYPE=RRHX X=8.50 Y=45008.50 Z=0.00 AZM=0.00 DIP=0.00 ACQCHAN=CH6
>HMEAS ID=5377.0537 CHTYPE=RRHY X=-8.50 Y=45008.50 Z=0.00 AZM=90.00 DIP=0.00 ACQCHAN=CH7
>=MTSECT
NFREQ=80
SECTID=14-IEB0537A
NCHAN=7
MAXBLKS=80
HX=5371.0537
HY=5372.0537
HZ=5373.0537
EX=5374.0537
EY=5375.0537
RRHX=5376.0537
RRHY=5377.0537
[6]:
print(tf_object)
Station: 14-IEB0537A
--------------------------------------------------
Survey: BOULIA
Project: IEB
Acquired by: Phoenix
Acquired date: 2014-07-28
Latitude: -22.824
Longitude: 139.295
Elevation: 158.000
Declination:
Value: 0.0
Model: WMM
Coordinate System: geographic
Impedance: True
Tipper: True
N Periods: 80
Period Range:
Min: 3.12500E-03 s
Max: 2.94118E+03 s
Frequency Range:
Min: 3.40000E-04 Hz
Max: 3.20000E+02 Hz
[7]:
tf_object.station_metadata
[7]:
{
"station": {
"acquired_by.author": "Phoenix",
"channels_recorded": [
"ex",
"ey",
"hx",
"hy",
"hz",
"rrhx",
"rrhy"
],
"comments": "CUTOFF: 0.00 COH: 35 % VAR: 25 %\nComp MTU box S/N Temp\nComp Chan# Sensor Azimuth\nEbat:12.3V Hbat:12.3V Rbat:11.9V\nElevation: 158 M. DECL: 0.000\nEx & Ey: MTU5A 2189 39 C\nEx1 1 100.0 M 0.0 DGtn\nEy1 2 100.0 M 90.0 DGtn\nFILE: IEB0537A IEB0564M\nHx & Hy: MTU5A 2189 39 C\nHx1 3 COIL2318 0.0 DGtn\nHy1 4 COIL2319 90.0 DGtn\nHz1 5 COIL2320\nHz: MTU5A 2189 39 C\nLat 22:49:254S Long 139:17:409E\nLat 22:49.423 S Lng 139:17.681 E\nMTU-Box Gains:E`s x 4 H`s x 4\nMTU-DFT VERSION: TStoFT.38\nMTU-RBS VERSION:R2012-0216-B22\nMTU-Ref Serial Number: U-2779\nNotch Filters set for 50 Hz.\nPROCESSED FROM DFT TIME SERIES\nRBS: 7 COH: 0.85 RHO VAR: 0.75\nRHx2 6 COIL2485 0.0 DGtn\nRHy2 7 COIL2487 90.0 DGtn\nRUN INFORMATION\nReference Field: Remote H - Ref.\nReference Site: IEB0564M\nRx & Ry: MTU5A 2779 40 C\nSTATION 1\nSYSTEM INFORMATION\nSite Desc; BadR: 0 SatR: 54\nSite Layout by:\nSite Permitted by:",
"data_type": "MT",
"geographic_name": null,
"id": "14-IEB0537A",
"location.datum": "WGS84",
"location.declination.model": "WMM",
"location.declination.value": 0.0,
"location.elevation": 158.0,
"location.latitude": -22.823722222222223,
"location.longitude": 139.29469444444445,
"orientation.method": null,
"orientation.reference_frame": "geographic",
"provenance.creation_time": "2014-08-01T00:00:00+00:00",
"provenance.software.author": null,
"provenance.software.name": "Phoenix",
"provenance.software.version": "MT-Editor Ver 0.99.2.106",
"provenance.submitter.author": "Phoenix",
"provenance.submitter.email": null,
"provenance.submitter.organization": null,
"release_license": "CC0-1.0",
"run_list": [
"14-IEB0537Aa"
],
"time_period.end": "2014-07-28T23:38:25+00:00",
"time_period.start": "2014-07-28T02:57:00+00:00",
"transfer_function.coordinate_system": "geopgraphic",
"transfer_function.id": "14-IEB0537A",
"transfer_function.processed_date": null,
"transfer_function.processing_parameters": [
"processing_parameters=[XPR Weighting=RHO Variance.]"
],
"transfer_function.remote_references": [],
"transfer_function.runs_processed": [
"14-IEB0537Aa"
],
"transfer_function.sign_convention": null,
"transfer_function.units": null
}
}
[8]:
tf_object.survey_metadata
[8]:
{
"survey": {
"acquired_by.author": "Phoenix",
"citation_dataset.doi": null,
"citation_journal.doi": null,
"datum": "WGS84",
"geographic_name": null,
"id": "BOULIA",
"name": null,
"northwest_corner.latitude": 0.0,
"northwest_corner.longitude": 0.0,
"project": "IEB",
"project_lead.email": null,
"project_lead.organization": null,
"release_license": "CC0-1.0",
"southeast_corner.latitude": 0.0,
"southeast_corner.longitude": 0.0,
"summary": null,
"time_period.end_date": "1980-01-01",
"time_period.start_date": "1980-01-01"
}
}
[9]:
tf_object.dataset
[9]:
<xarray.Dataset> Dimensions: (period: 80, output: 5, input: 5) Coordinates: * period (period) float64 0.003125 ... 2.941e+03 * output (output) <U2 'ex' 'ey' 'hx' 'hy' 'hz' * input (input) <U2 'ex' 'ey' 'hx' 'hy' 'hz' Data variables: transfer_function (period, output, input) complex128 (nan+na... transfer_function_error (period, output, input) float64 nan ... nan transfer_function_model_error (period, output, input) float64 nan ... nan inverse_signal_power (period, output, input) complex128 (nan+na... residual_covariance (period, output, input) complex128 (7.7123... Attributes: (12/14) survey: BOULIA project: IEB id: 14-IEB0537A name: None latitude: -22.823722222222223 longitude: 139.29469444444445 ... ... datum: WGS84 acquired_by: Phoenix start: 2014-07-28T02:57:00+00:00 end: 2014-07-28T23:38:25+00:00 runs_processed: ['14-IEB0537Aa'] coordinate_system: geographic
Common output by Metronix Geophysics
[27]:
from mt_metadata import TF_EDI_METRONIX
[28]:
tf_object = TF(fn=TF_EDI_METRONIX)
tf_object.read()
edi_object = EDI(fn=TF_EDI_METRONIX)
[29]:
print(edi_object.Header)
print(edi_object.Info)
print(edi_object.Measurement)
print(edi_object.Data)
>HEAD
ACQBY=Metronix
ACQDATE=2014-08-17T04:58:00+00:00
COORDINATE_SYSTEM=geographic
COUNTRY=Germany
DATAID=GEO858
DATUM=WGS84
ELEV=181.000
EMPTY=1e+32
ENDDATE=2014-08-17T20:03:00+00:00
FILEBY=Metronix
FILEDATE=2023-03-16
LAT=22:41:28.96
LON=139:42:18.14
PROGDATE=2021-12-01
PROGNAME=mt_metadata
PROGVERS=0.2.2
STATE=LX
STDVERS=SEG 1.0
UNITS=millivolts_per_kilometer_per_nanotesla
>INFO
>=DEFINEMEAS
MAXCHAN=9
MAXRUN=999
MAXMEAS=1000
REFLOC=Braunschweig
REFLAT=22.691378
REFLON=139.705040
REFELEV=181.0
REFTYPE=CART
UNITS=m
>EMEAS ID=1000.0001 CHTYPE=EX X=-50.00 Y=0.00 Z=0.00 X2=50.00 Y2=0.00 Z2=0.00 AZM=0.00 ACQCHAN=0.0
>EMEAS ID=1001.0001 CHTYPE=EY X=0.00 Y=-50.00 Z=0.00 X2=0.00 Y2=50.00 Z2=0.00 AZM=90.00 ACQCHAN=0.0
>HMEAS ID=1002.0001 CHTYPE=HX X=0.00 Y=0.00 Z=0.00 AZM=0.00 DIP=0.00 ACQCHAN=0.0
>HMEAS ID=1003.0001 CHTYPE=HY X=0.00 Y=0.00 Z=0.00 AZM=0.00 DIP=0.00 ACQCHAN=0.0
>HMEAS ID=1004.0001 CHTYPE=HZ X=0.00 Y=0.00 Z=0.00 AZM=0.00 DIP=0.00 ACQCHAN=0.0
>=MTSECT
NFREQ=73
SECTID=GEO858
NCHAN=0
MAXBLKS=999
EX=1000.0001
EY=1001.0001
HX=1002.0001
HY=1003.0001
HZ=1004.0001
RRHX=0
RRHY=0
[30]:
print(tf_object)
Station: GEO858
--------------------------------------------------
Survey: None
Project: None
Acquired by: Metronix
Acquired date: 2014-08-17
Latitude: 22.691
Longitude: 139.705
Elevation: 181.000
Declination:
Value: 0.0
Model: WMM
Coordinate System: geographic
Impedance: True
Tipper: True
N Periods: 73
Period Range:
Min: 5.15464E-03 s
Max: 1.44928E+03 s
Frequency Range:
Min: 6.90000E-04 Hz
Max: 1.94000E+02 Hz
[31]:
tf_object.station_metadata
[31]:
{
"station": {
"acquired_by.author": "Metronix",
"channels_recorded": [
"ex",
"ey",
"hx",
"hy",
"hz",
"rrhx",
"rrhy"
],
"data_type": "MT",
"geographic_name": null,
"id": "GEO858",
"location.datum": "WGS84",
"location.declination.model": "WMM",
"location.declination.value": 0.0,
"location.elevation": 181.0,
"location.latitude": 22.691378333333333,
"location.longitude": 139.70504,
"orientation.method": null,
"orientation.reference_frame": "geographic",
"provenance.creation_time": "2014-10-17T00:00:00+00:00",
"provenance.software.author": null,
"provenance.software.name": "Metronix",
"provenance.software.version": "Version 14 AUG 2014 SVN 1277 MINGW64",
"provenance.submitter.author": "Metronix",
"provenance.submitter.email": null,
"provenance.submitter.organization": null,
"release_license": "CC0-1.0",
"run_list": [
"GEO858a"
],
"time_period.end": "2014-08-17T20:03:00+00:00",
"time_period.start": "2014-08-17T04:58:00+00:00",
"transfer_function.coordinate_system": "geopgraphic",
"transfer_function.id": "GEO858",
"transfer_function.processed_date": null,
"transfer_function.processing_parameters": [],
"transfer_function.remote_references": [],
"transfer_function.runs_processed": [
"GEO858a"
],
"transfer_function.sign_convention": null,
"transfer_function.units": null
}
}
[32]:
tf_object.survey_metadata
[32]:
{
"survey": {
"acquired_by.author": "Metronix",
"citation_dataset.doi": null,
"citation_journal.doi": null,
"country": [
"Germany"
],
"datum": "WGS84",
"geographic_name": null,
"id": null,
"name": null,
"northwest_corner.latitude": 0.0,
"northwest_corner.longitude": 0.0,
"project": null,
"project_lead.email": null,
"project_lead.organization": null,
"release_license": "CC0-1.0",
"southeast_corner.latitude": 0.0,
"southeast_corner.longitude": 0.0,
"summary": null,
"time_period.end_date": "1980-01-01",
"time_period.start_date": "1980-01-01"
}
}
[33]:
tf_object.dataset
[33]:
<xarray.Dataset> Dimensions: (period: 73, output: 5, input: 5) Coordinates: * period (period) float64 0.005155 ... 1.449e+03 * output (output) <U2 'ex' 'ey' 'hx' 'hy' 'hz' * input (input) <U2 'ex' 'ey' 'hx' 'hy' 'hz' Data variables: transfer_function (period, output, input) complex128 (nan+na... transfer_function_error (period, output, input) float64 nan ... nan transfer_function_model_error (period, output, input) float64 nan ... nan inverse_signal_power (period, output, input) complex128 (nan+na... residual_covariance (period, output, input) complex128 0j ... 0j Attributes: (12/14) survey: None project: None id: GEO858 name: None latitude: 22.691378333333333 longitude: 139.70504 ... ... datum: WGS84 acquired_by: Metronix start: 2014-08-17T04:58:00+00:00 end: 2014-08-17T20:03:00+00:00 runs_processed: ['GEO858a'] coordinate_system: geographic
Common output by CGG
[34]:
from mt_metadata import TF_EDI_CGG
[35]:
tf_object = TF(fn=TF_EDI_CGG)
tf_object.read()
edi_object = EDI(fn=TF_EDI_CGG)
[36]:
print(edi_object.Header)
print(edi_object.Info)
print(edi_object.Measurement)
print(edi_object.Data)
>HEAD
ACQBY=GSC_CGG
ACQDATE=2014-06-05T00:00:00+00:00
COORDINATE_SYSTEM=geographic
DATAID=TEST01
DATUM=WGS84
ELEV=175.270
EMPTY=1e+32
FILEBY=None
FILEDATE=2023-03-16
LAT=-30:55:49.03
LOC=Australia
LON=127:13:45.23
PROGDATE=2021-12-01
PROGNAME=mt_metadata
PROGVERS=0.2.2
PROSPECT=EGC
STDVERS=SEG 1.0
UNITS=M
>INFO
H_SITE=E_SITE
PROCESSING PARAMETERS:
SITE INFO:
>=DEFINEMEAS
MAXCHAN=999
MAXRUN=999
MAXMEAS=999
REFLOC="TEST01"
REFLAT=-30.930285
REFLON=127.229230
REFELEV=175.27
REFTYPE=cartesian
UNITS=M
>HMEAS ID=1001.001 CHTYPE=HX X=0.00 Y=0.00 Z=0.00 AZM=0.00 DIP=0.00 ACQCHAN=0.0
>HMEAS ID=1002.001 CHTYPE=HY X=0.00 Y=0.00 Z=0.00 AZM=90.00 DIP=0.00 ACQCHAN=0.0
>HMEAS ID=1003.001 CHTYPE=HZ X=0.00 Y=0.00 Z=0.00 AZM=0.00 DIP=0.00 ACQCHAN=0.0
>EMEAS ID=1004.001 CHTYPE=EX X=0.00 Y=0.00 Z=0.00 X2=0.00 Y2=0.00 Z2=0.00 AZM=0.00 ACQCHAN=0.0
>EMEAS ID=1005.001 CHTYPE=EY X=0.00 Y=0.00 Z=0.00 X2=0.00 Y2=0.00 Z2=0.00 AZM=0.00 ACQCHAN=0.0
>HMEAS ID=1006.001 CHTYPE=RRHX X=0.00 Y=0.00 Z=0.00 AZM=0.00 DIP=0.00 ACQCHAN=0.0
>HMEAS ID=1007.001 CHTYPE=RRHY X=0.00 Y=0.00 Z=0.00 AZM=90.00 DIP=0.00 ACQCHAN=0.0
>=MTSECT
NFREQ=73
SECTID=None
NCHAN=0
MAXBLKS=999
EX=0
EY=0
HX=0
HY=0
HZ=0
RRHX=0
RRHY=0
[37]:
print(tf_object)
Station: TEST01
--------------------------------------------------
Survey: None
Project: EGC
Acquired by: GSC_CGG
Acquired date: 2014-06-05
Latitude: -30.930
Longitude: 127.229
Elevation: 175.270
Declination:
Value: 0.0
Model: WMM
Coordinate System: geographic
Impedance: True
Tipper: True
N Periods: 73
Period Range:
Min: 1.21153E-03 s
Max: 1.21153E+03 s
Frequency Range:
Min: 8.25404E-04 Hz
Max: 8.25404E+02 Hz
[38]:
tf_object.station_metadata
[38]:
{
"station": {
"acquired_by.author": "GSC_CGG",
"channels_recorded": [
"ex",
"ey",
"hx",
"hy",
"hz",
"rrhx",
"rrhy"
],
"comments": "H_SITE=E_SITE\nPROCESSING PARAMETERS:\nSITE INFO:",
"data_type": "MT",
"geographic_name": null,
"id": "TEST01",
"location.datum": "WGS84",
"location.declination.model": "WMM",
"location.declination.value": 0.0,
"location.elevation": 175.27,
"location.latitude": -30.930285,
"location.longitude": 127.22923,
"orientation.method": null,
"orientation.reference_frame": "geographic",
"provenance.creation_time": "2014-10-07T00:00:00+00:00",
"provenance.software.author": null,
"provenance.software.name": null,
"provenance.software.version": "Antlr3.Runtime:3.5.0.2;ContourEngine:1.0.41.8272;CoordinateSystemService:1.4.0.8439;DocumentCommon:1.4.0.8465;Fluent:2.1.0.0;GeoApi:1.7.0.0;hasp_net_windows:7.0.1.36183;Launcher:1.4.0.8471;MapDocument:1.4.0.8469;MTDocument:1.4.0.8459;MTDocumentDataProvider:1.4.0.8467;MTInversionCommon:1.4.0.8371;Ookii.Dialogs.Wpf:1.0.0.0;PlotElement:1.4.0.8466;PluginHost:1.4.0.8440;ProjNet:1.2.5085.21309;ShellEngine:1.4.0.8380;System.Windows.Interactivity:4.0.0.0;Utils:1.4.0.8449;Xceed.Wpf.AvalonDock:2.0.0.0;Xceed.Wpf.AvalonDock.Themes.Aero:2.0.0.0;Xceed.Wpf.Toolkit:1.9.0.0;",
"provenance.submitter.email": null,
"provenance.submitter.organization": null,
"release_license": "CC0-1.0",
"run_list": [
"TEST01a"
],
"time_period.end": "1980-01-01T00:00:00+00:00",
"time_period.start": "2014-06-05T00:00:00+00:00",
"transfer_function.coordinate_system": "geopgraphic",
"transfer_function.id": "TEST01",
"transfer_function.processed_date": null,
"transfer_function.processing_parameters": [
"processing_parameters=[NDec=1, NFFT=128, Ntype=1, RRType=None, RemoveLargeLines=true, RotMaxE=false]"
],
"transfer_function.remote_references": [],
"transfer_function.runs_processed": [
"TEST01a"
],
"transfer_function.sign_convention": null,
"transfer_function.units": null
}
}
[39]:
tf_object.survey_metadata
[39]:
{
"survey": {
"acquired_by.author": "GSC_CGG",
"citation_dataset.doi": null,
"citation_journal.doi": null,
"datum": "WGS84",
"geographic_name": "Australia",
"id": null,
"name": null,
"northwest_corner.latitude": 0.0,
"northwest_corner.longitude": 0.0,
"project": "EGC",
"project_lead.email": null,
"project_lead.organization": null,
"release_license": "CC0-1.0",
"southeast_corner.latitude": 0.0,
"southeast_corner.longitude": 0.0,
"summary": null,
"time_period.end_date": "1980-01-01",
"time_period.start_date": "1980-01-01"
}
}
[40]:
tf_object.dataset
[40]:
<xarray.Dataset> Dimensions: (period: 73, output: 5, input: 5) Coordinates: * period (period) float64 0.001212 ... 1.212e+03 * output (output) <U2 'ex' 'ey' 'hx' 'hy' 'hz' * input (input) <U2 'ex' 'ey' 'hx' 'hy' 'hz' Data variables: transfer_function (period, output, input) complex128 (nan+na... transfer_function_error (period, output, input) float64 nan ... nan transfer_function_model_error (period, output, input) float64 nan ... nan inverse_signal_power (period, output, input) complex128 (nan+na... residual_covariance (period, output, input) complex128 0j ... 0j Attributes: (12/14) survey: None project: EGC id: TEST01 name: None latitude: -30.930285 longitude: 127.22923 ... ... datum: WGS84 acquired_by: GSC_CGG start: 2014-06-05T00:00:00+00:00 end: 1980-01-01T00:00:00+00:00 runs_processed: ['TEST01a'] coordinate_system: geographic
Transfer Functions represented as Spectra
The spectra conversion is adopted from Anna Kelbert’s EMTF fortran code. This computes the inverse_signal_power
and residual_covariance
estimates for the errors to accurately estimate error upon rotation.
[41]:
from mt_metadata import TF_EDI_SPECTRA
[42]:
tf_object = TF(fn=TF_EDI_SPECTRA)
tf_object.read()
edi_object = EDI(fn=TF_EDI_SPECTRA)
[43]:
print(edi_object.Header)
print(edi_object.Info)
print(edi_object.Measurement)
print(edi_object.Data)
>HEAD
ACQBY=Quantec Consulting
ACQDATE=2004-07-03T00:00:00+00:00
COORDINATE_SYSTEM=geographic
COUNTRY=USA
COUNTY=Santa Fe
DATAID=SAGE_2005
DATUM=WGS84
ELEV=0.000
EMPTY=1e+32
ENDDATE=2004-07-06T00:00:00+00:00
FILEBY=Quantec Consulting
FILEDATE=2023-03-16
LAT=35:33:00.00
LON=-106:17:00.00
PROGDATE=2021-12-01
PROGNAME=mt_metadata
PROGVERS=0.2.2
STATE=New Mexico
STDVERS=1.0
UNITS=M
>INFO
MAXLINES=1000
>=DEFINEMEAS
MAXCHAN=7
MAXRUN=999
MAXMEAS=99999
REFLAT=35.550000
REFLON=-106.283333
REFELEV=0.0
REFTYPE=cartesian
UNITS=M
>HMEAS ID=11.001 CHTYPE=HX X=4858.00 Y=-3530.00 Z=0.00 AZM=107.00 DIP=0.00 ACQCHAN=0.0
>HMEAS ID=11.001 CHTYPE=RRHX X=4858.00 Y=-3530.00 Z=0.00 AZM=107.00 DIP=0.00 ACQCHAN=0.0
>HMEAS ID=12.001 CHTYPE=HY X=4858.00 Y=-3530.00 Z=0.00 AZM=-163.00 DIP=0.00 ACQCHAN=0.0
>HMEAS ID=12.001 CHTYPE=RRHY X=4858.00 Y=-3530.00 Z=0.00 AZM=-163.00 DIP=0.00 ACQCHAN=0.0
>HMEAS ID=13.001 CHTYPE=HZ X=4858.00 Y=-3530.00 Z=0.00 AZM=0.00 DIP=0.00 ACQCHAN=0.0
>EMEAS ID=14.001 CHTYPE=EX X=4872.00 Y=-3577.00 Z=0.00 X2=4843.00 Y2=-3482.00 Z2=0.00 AZM=106.98 ACQCHAN=0.0
>EMEAS ID=15.001 CHTYPE=EY X=4906.00 Y=-3515.00 Z=0.00 X2=4810.00 Y2=-3544.00 Z2=0.00 AZM=-163.19 ACQCHAN=0.0
>=MTSECT
NFREQ=33
SECTID=Ex
NCHAN=7
MAXBLKS=100
HX=11.001
HY=12.001
HZ=13.001
EX=14.001
EY=15.001
RRHX=11.001
RRHY=12.001
[44]:
print(tf_object)
Station: SAGE_2005
--------------------------------------------------
Survey: None
Project: None
Acquired by: Quantec Consulting
Acquired date: 2004-07-03
Latitude: 35.550
Longitude: -106.283
Elevation: 0.000
Declination:
Value: 0.0
Model: WMM
Coordinate System: geographic
Impedance: True
Tipper: True
N Periods: 33
Period Range:
Min: 4.19639E-03 s
Max: 2.09732E+02 s
Frequency Range:
Min: 4.76800E-03 Hz
Max: 2.38300E+02 Hz
[45]:
tf_object.station_metadata
[45]:
{
"station": {
"acquired_by.author": "Quantec Consulting",
"channels_recorded": [
"ex",
"ey",
"hx",
"hy",
"hz",
"rrhx",
"rrhy"
],
"comments": "MAXLINES=1000",
"data_type": "MT",
"geographic_name": null,
"id": "SAGE_2005",
"location.datum": "WGS84",
"location.declination.model": "WMM",
"location.declination.value": 0.0,
"location.elevation": 0.0,
"location.latitude": 35.55,
"location.longitude": -106.28333333333333,
"orientation.method": null,
"orientation.reference_frame": "geographic",
"provenance.creation_time": "2001-01-31T00:00:00+00:00",
"provenance.software.author": null,
"provenance.software.name": "Quantec Consulting",
"provenance.software.version": "GEOTOOLS_2.3",
"provenance.submitter.author": "Quantec Consulting",
"provenance.submitter.email": null,
"provenance.submitter.organization": null,
"release_license": "CC0-1.0",
"run_list": [
"SAGE_2005a"
],
"time_period.end": "2004-07-06T00:00:00+00:00",
"time_period.start": "2004-07-03T00:00:00+00:00",
"transfer_function.coordinate_system": "geopgraphic",
"transfer_function.id": "SAGE_2005",
"transfer_function.processed_date": null,
"transfer_function.processing_parameters": [],
"transfer_function.remote_references": [],
"transfer_function.runs_processed": [
"SAGE_2005a"
],
"transfer_function.sign_convention": null,
"transfer_function.units": null
}
}
[46]:
tf_object.survey_metadata
[46]:
{
"survey": {
"acquired_by.author": "Quantec Consulting",
"citation_dataset.doi": null,
"citation_journal.doi": null,
"country": [
"USA"
],
"datum": "WGS84",
"geographic_name": null,
"id": null,
"name": null,
"northwest_corner.latitude": 0.0,
"northwest_corner.longitude": 0.0,
"project": null,
"project_lead.email": null,
"project_lead.organization": null,
"release_license": "CC0-1.0",
"southeast_corner.latitude": 0.0,
"southeast_corner.longitude": 0.0,
"summary": null,
"time_period.end_date": "1980-01-01",
"time_period.start_date": "1980-01-01"
}
}
[47]:
tf_object.dataset
[47]:
<xarray.Dataset> Dimensions: (period: 33, output: 5, input: 5) Coordinates: * period (period) float64 0.004196 0.005952 ... 209.7 * output (output) <U2 'ex' 'ey' 'hx' 'hy' 'hz' * input (input) <U2 'ex' 'ey' 'hx' 'hy' 'hz' Data variables: transfer_function (period, output, input) complex128 (nan+na... transfer_function_error (period, output, input) float64 nan ... nan transfer_function_model_error (period, output, input) float64 nan ... nan inverse_signal_power (period, output, input) complex128 (nan+na... residual_covariance (period, output, input) complex128 (0.0062... Attributes: (12/14) survey: None project: None id: SAGE_2005 name: None latitude: 35.55 longitude: -106.28333333333333 ... ... datum: WGS84 acquired_by: Quantec Consulting start: 2004-07-03T00:00:00+00:00 end: 2004-07-06T00:00:00+00:00 runs_processed: ['SAGE_2005a'] coordinate_system: geographic
[ ]: