seismology-RUB/DataAnalysis2022
10
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

Compare commits

merge into: seismology-RUB:main
Branches Tags
seismology-RUB:main
seismology-RUB:feature/2024
...
pull from: seismology-RUB:feature/2024
Branches Tags
seismology-RUB:feature/2024
seismology-RUB:main

5 Commits
main ... feature/20

Author SHA1 Message Date
Kasper D. Fischer
3be655017f updated paths to header image file 2024-06-10 12:15:28 +02:00
Kasper D. Fischer
228434397d add notebook collection from seismo-live.org on ObsPy 
activate LFS
 2024-06-10 11:47:40 +02:00
Kasper D. Fischer
dbcb4a3d8a initial import of PPSD exercise from seismo live web site 2024-06-03 12:17:07 +02:00
Kasper D. Fischer
d829006103 update README.md 2024-04-22 11:46:34 +02:00
Kasper D. Fischer
fc3b809875 Add data analysis code workspace configuration 
This commit adds a new file `dataanalaysis.code-workspace` which includes the necessary configuration for a data analysis project. The file specifies the root folder path and contains default settings.
 2024-04-22 11:35:38 +02:00
42 changed files with 36313 additions and 2 deletions
Show Stats Expand all files Collapse all files

3
.gitattributes vendored Normal file
View File

@ -0,0 +1,3 @@
*.png filter=lfs diff=lfs merge=lfs -text
*.mseed filter=lfs diff=lfs merge=lfs -text
*.sac filter=lfs diff=lfs merge=lfs -text

223
04-PPSD/1-probabilistic_power_spectral_densities.ipynb Normal file
View File

@ -0,0 +1,223 @@
{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"<div style='background-image: url(\"../images/header.svg\") ; padding: 0px ; background-size: cover ; border-radius: 5px ; height: 250px'>\n",
" <div style=\"float: right ; margin: 50px ; padding: 20px ; background: rgba(255 , 255 , 255 , 0.7) ; width: 50% ; height: 150px\">\n",
" <div style=\"position: relative ; top: 50% ; transform: translatey(-50%)\">\n",
" <div style=\"font-size: xx-large ; font-weight: 900 ; color: rgba(0 , 0 , 0 , 0.8) ; line-height: 100%\">Noise</div>\n",
" <div style=\"font-size: large ; padding-top: 20px ; color: rgba(0 , 0 , 0 , 0.5)\">Lab: Probabilistic Power Spectral Densities</div>\n",
" </div>\n",
" </div>\n",
"</div>\n",
"\n",
"\n",
"Seismo-Live: http://seismo-live.org\n",
"\n",
"##### Authors:\n",
"* Tobias Megies ([@megies](https://github.com/megies))\n",
"\n",
"---"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"%matplotlib inline\n",
"import matplotlib.pyplot as plt\n",
"plt.style.use(\"bmh\")\n",
"plt.rcParams['figure.figsize'] = 10, 6"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
" * read waveform data from file `data/GR.FUR..BHN.D.2015.361` (station `FUR`, [LMU geophysical observatory in Fürstenfeldbruck](https://www.geophysik.uni-muenchen.de/observatory/seismology))\n",
" * read corresponding station metadata from file `data/station_FUR.stationxml`\n",
" * print info on both waveforms and station metadata"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"from obspy import read, read_inventory\n",
"\n",
"st = read(\"data/GR.FUR..BHN.D.2015.361\")\n",
"inv = read_inventory(\"data/station_FUR.stationxml\")\n",
"\n",
"print(st)\n",
"print(inv)\n",
"inv.plot(projection=\"ortho\");"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
" * compute probabilistic power spectral densities using `PPSD` class from obspy.signal, see http://docs.obspy.org/tutorial/code_snippets/probabilistic_power_spectral_density.html (but use the inventory you read from StationXML as metadata)\n",
" * plot the processed `PPSD` (`plot()` method attached to `PPSD` object)"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"from obspy.signal import PPSD\n",
"\n",
"tr = st[0]\n",
"ppsd = PPSD(stats=tr.stats, metadata=inv)\n",
"\n",
"ppsd.add(tr)\n",
"ppsd.plot()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Since longer term stacks would need too much waveform data and take way too long to compute, we prepared one year continuous data preprocessed for a single channel of station `FUR` to play with..\n",
"\n",
" * load long term pre-computed PPSD from file `PPSD_FUR_HHN.npz` using `PPSD`'s `load_npz()` staticmethod (i.e. it is called directly from the class, not an instance object of the class)\n",
" * plot the PPSD (default is full time-range, depending on how much data and spread is in the data, adjust `max_percentage` option of `plot()` option) (might take a couple of minutes..!)\n",
" * do a cumulative plot (which is good to judge non-exceedance percentage dB thresholds)"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"from obspy.signal import PPSD\n",
"\n",
"ppsd = PPSD.load_npz(\"data/PPSD_FUR_HHN.npz\", allow_pickle=True)"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"ppsd.plot(max_percentage=10)\n",
"ppsd.plot(cumulative=True)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
" * do different stacks of the data using the [`calculate_histogram()` (see docs!)](http://docs.obspy.org/packages/autogen/obspy.signal.spectral_estimation.PPSD.calculate_histogram.html) method of `PPSD` and visualize them\n",
" * compare differences in different frequency bands qualitatively (anthropogenic vs. \"natural\" noise)..\n",
" * nighttime stack, daytime stack\n",
" * advanced exercise: Use the `callback` option and use some crazy custom callback function in `calculate_histogram()`, e.g. stack together all data from birthdays in your family.. or all German holidays + Sundays in the time span.. or from dates of some bands' concerts on a tour.. etc."
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"ppsd.calculate_histogram(time_of_weekday=[(-1, 0, 2), (-1, 22, 24)])\n",
"ppsd.plot(max_percentage=10)\n",
"ppsd.calculate_histogram(time_of_weekday=[(-1, 8, 16)])\n",
"ppsd.plot(max_percentage=10)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
" * do different stacks of the data using the [`calculate_histogram()` (see docs!)](http://docs.obspy.org/packages/autogen/obspy.signal.spectral_estimation.PPSD.calculate_histogram.html) method of `PPSD` and visualize them\n",
" * compare differences in different frequency bands qualitatively (anthropogenic vs. \"natural\" noise)..\n",
" * weekdays stack, weekend stack"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"ppsd.calculate_histogram(time_of_weekday=[(1, 0, 24), (2, 0, 24), (3, 0, 24), (4, 0, 24), (5, 0, 24)])\n",
"ppsd.plot(max_percentage=10)\n",
"ppsd.calculate_histogram(time_of_weekday=[(6, 0, 24), (7, 0, 24)])\n",
"ppsd.plot(max_percentage=10)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
" * do different stacks of the data using the [`calculate_histogram()` (see docs!)](http://docs.obspy.org/packages/autogen/obspy.signal.spectral_estimation.PPSD.calculate_histogram.html) method of `PPSD` and visualize them\n",
" * compare differences in different frequency bands qualitatively (anthropogenic vs. \"natural\" noise)..\n",
" * seasonal stacks (e.g. northern hemisphere autumn vs. spring/summer, ...)"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"ppsd.calculate_histogram(month=[10, 11, 12, 1])\n",
"ppsd.plot(max_percentage=10)\n",
"ppsd.calculate_histogram(month=[4, 5, 6, 7])\n",
"ppsd.plot(max_percentage=10)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
" * do different stacks of the data using the [`calculate_histogram()` (see docs!)](http://docs.obspy.org/packages/autogen/obspy.signal.spectral_estimation.PPSD.calculate_histogram.html) method of `PPSD` and visualize them\n",
" * compare differences in different frequency bands qualitatively (anthropogenic vs. \"natural\" noise)..\n",
" * stacks by specific month\n",
" * maybe even combine several of above restrictions.. (e.g. only nighttime on weekends)"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"jupytext": {
"encoding": "# -*- coding: utf-8 -*-"
},
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.12.3"
}
},
"nbformat": 4,
"nbformat_minor": 2
}

BIN
04-PPSD/data/GR.FUR..BHN.D.2015.361 Normal file
View File

Binary file not shown.

BIN
04-PPSD/data/PPSD_FUR_HHN.npz Normal file
View File

Binary file not shown.

232
04-PPSD/data/station_FUR.stationxml Normal file
View File

@ -0,0 +1,232 @@
<?xml version='1.0' encoding='UTF-8'?>
<FDSNStationXML xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns="http://www.fdsn.org/xml/station/1" schemaVersion="1.0" xsi:schemaLocation="http://www.fdsn.org/xml/station/1 http://www.fdsn.org/xml/station/fdsn-station-1.0.xsd">
<Source>Jane</Source>
<Sender>Jane</Sender>
<Module>JANE WEB SERVICE: fdsnws-station | Jane version: 0.0.0-gd6e1</Module>
<ModuleURI>http://jane/fdsnws/station/1/query?network=GR&amp;station=FUR&amp;starttime=2015-12-01&amp;endtime=2016-02-01&amp;level=response</ModuleURI>
<Created>2016-05-03T10:14:19+00:00</Created>
<Network startDate="2006-12-16T00:00:00.000000Z" code="GR">
<Description>GRSN</Description>
<SelectedNumberStations>1</SelectedNumberStations>
<TotalNumberStations>4</TotalNumberStations>
<Station startDate="2006-12-16T00:00:00.000000Z" code="FUR">
<Latitude plusError="0.0" minusError="0.0">48.162899</Latitude>
<Longitude plusError="0.0" minusError="0.0">11.2752</Longitude>
<Elevation plusError="0.0" minusError="0.0">565.0</Elevation>
<Site>
<Name>Fuerstenfeldbruck, Bavaria, GR-Net</Name>
</Site>
<CreationDate>2006-12-16T00:00:00.000</CreationDate>
<SelectedNumberChannels>12</SelectedNumberChannels>
<TotalNumberChannels>12</TotalNumberChannels>
<Channel locationCode="" code="HHN" startDate="2006-12-16T00:00:00.000">
<Latitude plusError="0.0" minusError="0.0">48.162899</Latitude>
<Longitude plusError="0.0" minusError="0.0">11.2752</Longitude>
<Elevation plusError="0.0" minusError="0.0">565.0</Elevation>
<Depth plusError="0.0" minusError="0.0">0.0</Depth>
<Azimuth plusError="0.0" minusError="0.0">0.0</Azimuth>
<Dip plusError="0.0" minusError="0.0">0.0</Dip>
<Type>TRIGGERED</Type>
<Type>GEOPHYSICAL</Type>
<SampleRate plusError="0.0" minusError="0.0">100.0</SampleRate>
<ClockDrift plusError="0.0" minusError="0.0">0.02</ClockDrift>
<CalibrationUnits>
<Name>A</Name>
<Description>Amperes</Description>
</CalibrationUnits>
<Sensor>
<Type>Streckeisen STS-2/N seismometer</Type>
</Sensor>
<Response>
<InstrumentSensitivity>
<Value>9.4368E8</Value>
<Frequency>0.02</Frequency>
<InputUnits>
<Name>M/S</Name>
<Description>Velocity in Meters per Second</Description>
</InputUnits>
<OutputUnits>
<Name>COUNTS</Name>
<Description>Digital Counts</Description>
</OutputUnits>
</InstrumentSensitivity>
<Stage number="1">
<PolesZeros>
<InputUnits>
<Name>M/S</Name>
<Description>Velocity in Meters per Second</Description>
</InputUnits>
<OutputUnits>
<Name>V</Name>
<Description>Volts</Description>
</OutputUnits>
<PzTransferFunctionType>LAPLACE (RADIANS/SECOND)</PzTransferFunctionType>
<NormalizationFactor>6.0077E7</NormalizationFactor>
<NormalizationFrequency plusError="0.0" minusError="0.0">1.0</NormalizationFrequency>
<Zero number="0">
<Real plusError="0.0" minusError="0.0">0.0</Real>
<Imaginary plusError="0.0" minusError="0.0">0.0</Imaginary>
</Zero>
<Zero number="1">
<Real plusError="0.0" minusError="0.0">0.0</Real>
<Imaginary plusError="0.0" minusError="0.0">0.0</Imaginary>
</Zero>
<Pole number="0">
<Real plusError="0.0" minusError="0.0">-0.037004</Real>
<Imaginary plusError="0.0" minusError="0.0">0.037016</Imaginary>
</Pole>
<Pole number="1">
<Real plusError="0.0" minusError="0.0">-0.037004</Real>
<Imaginary plusError="0.0" minusError="0.0">-0.037016</Imaginary>
</Pole>
<Pole number="2">
<Real plusError="0.0" minusError="0.0">-251.33</Real>
<Imaginary plusError="0.0" minusError="0.0">0.0</Imaginary>
</Pole>
<Pole number="3">
<Real plusError="0.0" minusError="0.0">-131.04</Real>
<Imaginary plusError="0.0" minusError="0.0">-467.29</Imaginary>
</Pole>
<Pole number="4">
<Real plusError="0.0" minusError="0.0">-131.04</Real>
<Imaginary plusError="0.0" minusError="0.0">467.29</Imaginary>
</Pole>
</PolesZeros>
<StageGain>
<Value>1500.0</Value>
<Frequency>0.02</Frequency>
</StageGain>
</Stage>
<Stage number="2">
<Coefficients>
<InputUnits>
<Name>V</Name>
<Description>Volts</Description>
</InputUnits>
<OutputUnits>
<Name>COUNTS</Name>
<Description>Digital Counts</Description>
</OutputUnits>
<CfTransferFunctionType>DIGITAL</CfTransferFunctionType>
</Coefficients>
<Decimation>
<InputSampleRate plusError="0.0" minusError="0.0">200.0</InputSampleRate>
<Factor>1</Factor>
<Offset>0</Offset>
<Delay plusError="0.0" minusError="0.0">0.0</Delay>
<Correction plusError="0.0" minusError="0.0">0.0</Correction>
</Decimation>
<StageGain>
<Value>629121.0</Value>
<Frequency>0.0</Frequency>
</StageGain>
</Stage>
</Response>
</Channel>
<Channel locationCode="" code="BHN" startDate="2006-12-16T00:00:00.000">
<Latitude plusError="0.0" minusError="0.0">48.162899</Latitude>
<Longitude plusError="0.0" minusError="0.0">11.2752</Longitude>
<Elevation plusError="0.0" minusError="0.0">565.0</Elevation>
<Depth plusError="0.0" minusError="0.0">0.0</Depth>
<Azimuth plusError="0.0" minusError="0.0">0.0</Azimuth>
<Dip plusError="0.0" minusError="0.0">0.0</Dip>
<Type>TRIGGERED</Type>
<Type>GEOPHYSICAL</Type>
<SampleRate plusError="0.0" minusError="0.0">20.0</SampleRate>
<ClockDrift plusError="0.0" minusError="0.0">0.02</ClockDrift>
<CalibrationUnits>
<Name>A</Name>
<Description>Amperes</Description>
</CalibrationUnits>
<Sensor>
<Type>Streckeisen STS-2/N seismometer</Type>
</Sensor>
<Response>
<InstrumentSensitivity>
<Value>9.4368E8</Value>
<Frequency>0.02</Frequency>
<InputUnits>
<Name>M/S</Name>
<Description>Velocity in Meters per Second</Description>
</InputUnits>
<OutputUnits>
<Name>COUNTS</Name>
<Description>Digital Counts</Description>
</OutputUnits>
</InstrumentSensitivity>
<Stage number="1">
<PolesZeros>
<InputUnits>
<Name>M/S</Name>
<Description>Velocity in Meters per Second</Description>
</InputUnits>
<OutputUnits>
<Name>V</Name>
<Description>Volts</Description>
</OutputUnits>
<PzTransferFunctionType>LAPLACE (RADIANS/SECOND)</PzTransferFunctionType>
<NormalizationFactor>6.0077E7</NormalizationFactor>
<NormalizationFrequency plusError="0.0" minusError="0.0">1.0</NormalizationFrequency>
<Zero number="0">
<Real plusError="0.0" minusError="0.0">0.0</Real>
<Imaginary plusError="0.0" minusError="0.0">0.0</Imaginary>
</Zero>
<Zero number="1">
<Real plusError="0.0" minusError="0.0">0.0</Real>
<Imaginary plusError="0.0" minusError="0.0">0.0</Imaginary>
</Zero>
<Pole number="0">
<Real plusError="0.0" minusError="0.0">-0.037004</Real>
<Imaginary plusError="0.0" minusError="0.0">0.037016</Imaginary>
</Pole>
<Pole number="1">
<Real plusError="0.0" minusError="0.0">-0.037004</Real>
<Imaginary plusError="0.0" minusError="0.0">-0.037016</Imaginary>
</Pole>
<Pole number="2">
<Real plusError="0.0" minusError="0.0">-251.33</Real>
<Imaginary plusError="0.0" minusError="0.0">0.0</Imaginary>
</Pole>
<Pole number="3">
<Real plusError="0.0" minusError="0.0">-131.04</Real>
<Imaginary plusError="0.0" minusError="0.0">-467.29</Imaginary>
</Pole>
<Pole number="4">
<Real plusError="0.0" minusError="0.0">-131.04</Real>
<Imaginary plusError="0.0" minusError="0.0">467.29</Imaginary>
</Pole>
</PolesZeros>
<StageGain>
<Value>1500.0</Value>
<Frequency>0.02</Frequency>
</StageGain>
</Stage>
<Stage number="2">
<Coefficients>
<InputUnits>
<Name>V</Name>
<Description>Volts</Description>
</InputUnits>
<OutputUnits>
<Name>COUNTS</Name>
<Description>Digital Counts</Description>
</OutputUnits>
<CfTransferFunctionType>DIGITAL</CfTransferFunctionType>
</Coefficients>
<Decimation>
<InputSampleRate plusError="0.0" minusError="0.0">200.0</InputSampleRate>
<Factor>1</Factor>
<Offset>0</Offset>
<Delay plusError="0.0" minusError="0.0">0.0</Delay>
<Correction plusError="0.0" minusError="0.0">0.0</Correction>
</Decimation>
<StageGain>
<Value>629121.0</Value>
<Frequency>0.0</Frequency>
</StageGain>
</Stage>
</Response>
</Channel>
</Station>
</Network>
</FDSNStationXML>

69
ObsPy/00_Introduction.ipynb Normal file
View File

@ -0,0 +1,69 @@
{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"<div style='background-image: url(\"../images/header.svg\") ; padding: 0px ; background-size: cover ; border-radius: 5px ; height: 250px'>\n",
" <div style=\"float: right ; margin: 50px ; padding: 20px ; background: rgba(255 , 255 , 255 , 0.7) ; width: 50% ; height: 150px\">\n",
" <div style=\"position: relative ; top: 50% ; transform: translatey(-50%)\">\n",
" <div style=\"font-size: xx-large ; font-weight: 900 ; color: rgba(0 , 0 , 0 , 0.8) ; line-height: 100%\">ObsPy Tutorial</div>\n",
" <div style=\"font-size: large ; padding-top: 20px ; color: rgba(0 , 0 , 0 , 0.5)\">Introduction/Index</div>\n",
" </div>\n",
" </div>\n",
"</div>"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"![](images/obspy_logo_full_524x179px.png)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"These notebooks are an introduction to ObsPy and will take about 2 days to fully complete. Some basic knowledge of scientific Python is assumed."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Timetable/Contents\n",
"\n",
"1. **[Introduction to File Formats and read/write in ObsPy](01_File_Formats.ipynb)** [45 Minutes] - *[[Solution]](01_File_Formats-with_solutions.ipynb)*\n",
"2. **[Handling Time](02_UTCDateTime.ipynb)** [35 Minutes] - *[[Solution]](02_UTCDateTime-with_solutions.ipynb)*\n",
"3. **[Handling Waveform Data: Stream & Trace Classes](03_waveform_data.ipynb)** [45 Minutes] - *[[Solution]](03_waveform_data-with_solutions.ipynb)*\n",
"4. **[Handling Station Metainformation](04_Station_metainformation.ipynb)** [15 Minutes] - *[[Solution]](04_Station_metainformation-with_solutions.ipynb)*\n",
"5. **[Handling Event Metainformation](05_Event_metadata.ipynb)** [25 Minutes] - *[[Solution]](05_Event_metadata-with_solutions.ipynb)*\n",
"6. **[Retrieving Data from Data Centers](06_FDSN.ipynb)** [30 Minutes] - *[[Solution]](06_FDSN-with_solutions.ipynb)*\n",
"7. **[Basic Downloading/Processing Exercise](07_Basic_Processing_Exercise.ipynb)** [60 - x Minutes] - *[[Solution]](07_Basic_Processing_Exercise-with_solutions.ipynb)*\n",
"8. **[Exercise: Event Detection and Magnitude Estimation in an Aftershock Series](08_Exercise__2008_MtCarmel_Earthquake_and_Aftershock_Series.ipynb)** [open end] - *[[Solution]](08_Exercise__2008_MtCarmel_Earthquake_and_Aftershock_Series-with_solutions.ipynb)*"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.10.14"
}
},
"nbformat": 4,
"nbformat_minor": 2
}

1
ObsPy/01_File_Formats.ipynb Normal file
View File

File diff suppressed because one or more lines are too long

1
ObsPy/02_UTCDateTime.ipynb Normal file
View File

File diff suppressed because one or more lines are too long

1
ObsPy/03_waveform_data.ipynb Normal file
View File

File diff suppressed because one or more lines are too long

1
ObsPy/04_Station_metainformation.ipynb Normal file
View File

@ -0,0 +1 @@
{"cells":[{"cell_type":"markdown","metadata":{},"source":["<div style='background-image: url(\"../images/header.svg\") ; padding: 0px ; background-size: cover ; border-radius: 5px ; height: 250px'>\n"," <div style=\"float: right ; margin: 50px ; padding: 20px ; background: rgba(255 , 255 , 255 , 0.7) ; width: 50% ; height: 150px\">\n"," <div style=\"position: relative ; top: 50% ; transform: translatey(-50%)\">\n"," <div style=\"font-size: xx-large ; font-weight: 900 ; color: rgba(0 , 0 , 0 , 0.8) ; line-height: 100%\">ObsPy Tutorial</div>\n"," <div style=\"font-size: large ; padding-top: 20px ; color: rgba(0 , 0 , 0 , 0.5)\">Handling Station Metadata</div>\n"," </div>\n"," </div>\n","</div>"]},{"cell_type":"markdown","metadata":{},"source":["Seismo-Live: http://seismo-live.org\n","\n","##### Authors:\n","* Lion Krischer ([@krischer](https://github.com/krischer))\n","* Tobias Megies ([@megies](https://github.com/megies))\n","---"]},{"cell_type":"markdown","metadata":{},"source":["![](images/obspy_logo_full_524x179px.png)"]},{"cell_type":"code","execution_count":null,"metadata":{},"outputs":[],"source":["%matplotlib inline\n","import matplotlib.pyplot as plt\n","plt.style.use('ggplot')\n","plt.rcParams['figure.figsize'] = 12, 8"]},{"cell_type":"markdown","metadata":{},"source":["- for station metadata, the de-facto standard of the future (replacing SEED/RESP) is [FDSN StationXML](http://www.fdsn.org/xml/station/)\n","- FDSN StationXML files can be read using **`read_inventory()`**"]},{"cell_type":"code","execution_count":null,"metadata":{},"outputs":[],"source":["from obspy import read_inventory\n","# real-world StationXML files often deviate from the official schema definition\n","# therefore file-format autodiscovery sometimes fails and we have to force the file format\n","inventory = read_inventory(\"./data/station_PFO.xml\", format=\"STATIONXML\")\n","print(type(inventory))"]},{"cell_type":"markdown","metadata":{},"source":["- the nested ObsPy Inventory class structure (Inventory/Station/Channel/Response/...) is closely modelled after FDSN StationXML\n","<img src=\"images/Inventory.svg\" width=90%>"]},{"cell_type":"code","execution_count":null,"metadata":{},"outputs":[],"source":["!head data/station_BFO.xml"]},{"cell_type":"code","execution_count":null,"metadata":{},"outputs":[],"source":["print(inventory)"]},{"cell_type":"code","execution_count":null,"metadata":{},"outputs":[],"source":["network = inventory[0]\n","print(network)"]},{"cell_type":"code","execution_count":null,"metadata":{},"outputs":[],"source":["station = network[0]\n","print(station)"]},{"cell_type":"code","execution_count":null,"metadata":{},"outputs":[],"source":["channel = station[0]\n","print(channel)"]},{"cell_type":"code","execution_count":null,"metadata":{},"outputs":[],"source":["print(channel.response)"]},{"cell_type":"code","execution_count":null,"metadata":{},"outputs":[],"source":["from obspy import read\n","st = read(\"./data/waveform_PFO.mseed\")\n","print(st)"]},{"cell_type":"code","execution_count":null,"metadata":{},"outputs":[],"source":["inv = read_inventory(\"./data/station_PFO.xml\", format=\"STATIONXML\")"]},{"cell_type":"code","execution_count":null,"metadata":{},"outputs":[],"source":["print(st[0].stats)"]},{"cell_type":"markdown","metadata":{},"source":["- the instrument response can be deconvolved from the waveform data using the convenience method **`Stream.remove_response()`**\n","- evalresp is used internally to calculate the instrument response"]},{"cell_type":"code","execution_count":null,"metadata":{},"outputs":[],"source":["st.plot()\n","st.remove_response(inventory=inv)\n","st.plot()"]},{"cell_type":"markdown","metadata":{},"source":["- several options can be used to specify details of the deconvolution (water level, frequency domain prefiltering), output units (velocity/displacement/acceleration), demeaning, tapering and to specify if any response stages should be omitted"]},{"cell_type":"code","execution_count":null,"metadata":{},"outputs":[],"source":["st = read(\"./data/waveform_PFO.mseed\")\n","st.remove_response(inventory=inv, water_level=60, pre_filt=(0.01, 0.02, 8, 10), output=\"DISP\")\n","st.plot()"]},{"cell_type":"markdown","metadata":{},"source":["- station metadata not present in StationXML yet but in Dataless SEED or RESP files can be used for instrument correction using the `.simulate()` method of Stream/Trace in a similar fashion"]}],"metadata":{"kernelspec":{"display_name":"Python 3","language":"python","name":"python3"}},"nbformat":4,"nbformat_minor":2}

1
ObsPy/05_Event_metadata.ipynb Normal file
View File

File diff suppressed because one or more lines are too long

1
ObsPy/06_FDSN.ipynb Normal file
View File

File diff suppressed because one or more lines are too long

1
ObsPy/07_Basic_Processing_Exercise.ipynb Normal file
View File

File diff suppressed because one or more lines are too long

1
ObsPy/08_Exercise__2008_MtCarmel_Earthquake_and_Aftershock_Series.ipynb Normal file
View File

File diff suppressed because one or more lines are too long

12315
ObsPy/data/2014.ndk Normal file
View File

File diff suppressed because it is too large Load Diff

13
ObsPy/data/GCMT_2014_04_01__Mw_8_1 Normal file
View File

@ -0,0 +1,13 @@
PDEW 2014 4 1 23 46 47.30 -19.6100 -70.7700 25.0 0.0 8.2 NEAR COAST OF NORTHERN C
event name: 201404012346A
time shift: 44.2000
half duration: 28.0000
latitude: -19.7000
longitude: -70.8100
depth: 21.6000
Mrr: 9.200000e+27
Mtt: -3.900000e+26
Mpp: -8.810000e+27
Mrt: 6.370000e+27
Mrp: -1.530000e+28
Mtp: 2.050000e+27

174
ObsPy/data/GCMT_2014_04_01__Mw_8_1.xml Normal file
View File

@ -0,0 +1,174 @@
<q:quakeml xsi:schemaLocation="http://quakeml.org/schema/xsd http://quakeml.org/schema/xsd/QuakeML-1.2.xsd" xmlns="http://quakeml.org/xmlns/bed/1.2" xmlns:q="http://quakeml.org/xmlns/quakeml/1.2" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance">
<eventParameters publicID="smi:service.iris.edu/fdsnws/event/1/query">
<event publicID="smi:service.iris.edu/fdsnws/event/1/query?eventid=4597319">
<description>
<text>NEAR COAST OF NORTHERN C</text>
<type>region name</type>
</description>
<focalMechanism publicID="smi:ds.iris.edu/spudservice/momenttensor/gcmtid/C201404012346A/quakeml#focalmechanism">
<momentTensor publicID="smi:ds.iris.edu/spudservice/momenttensor/gcmtid/C201404012346A/quakeml#momenttensor">
<dataUsed>
<waveType>body waves</waveType>
<stationCount>166</stationCount>
<componentCount>407</componentCount>
<shortestPeriod>50.0</shortestPeriod>
</dataUsed>
<dataUsed>
<waveType>surface waves</waveType>
<stationCount>139</stationCount>
<componentCount>219</componentCount>
<shortestPeriod>50.0</shortestPeriod>
</dataUsed>
<dataUsed>
<waveType>mantle waves</waveType>
<stationCount>170</stationCount>
<componentCount>483</componentCount>
<shortestPeriod>200.0</shortestPeriod>
</dataUsed>
<derivedOriginID>smi:www.iris.edu/spudservice/momenttensor/gcmtid/C201404012346A#cmtorigin</derivedOriginID>
<momentMagnitudeID>smi:ds.iris.edu/spudservice/momenttensor/gcmtid/C201404012346A/quakeml#magnitude</momentMagnitudeID>
<scalarMoment>
<value>1898000000000000000000</value>
</scalarMoment>
<tensor>
<Mrr>
<value>920000000000000000000</value>
<uncertainty>4000000000000000000</uncertainty>
</Mrr>
<Mtt>
<value>-39000000000000000000</value>
<uncertainty>2000000000000000000</uncertainty>
</Mtt>
<Mpp>
<value>-881000000000000000000</value>
<uncertainty>2000000000000000000</uncertainty>
</Mpp>
<Mrt>
<value>637000000000000000000</value>
<uncertainty>23000000000000000000</uncertainty>
</Mrt>
<Mrp>
<value>-1530000000000000000000</value>
<uncertainty>33000000000000000000</uncertainty>
</Mrp>
<Mtp>
<value>205000000000000000000</value>
<uncertainty>1000000000000000000</uncertainty>
</Mtp>
</tensor>
<sourceTimeFunction>
<type>triangle</type>
<duration>56.0</duration>
</sourceTimeFunction>
</momentTensor>
<triggeringOriginID>smi:www.iris.edu/spudservice/momenttensor/gcmtid/C201404012346A#reforigin</triggeringOriginID>
<nodalPlanes>
<nodalPlane1>
<strike>
<value>355</value>
</strike>
<dip>
<value>15</value>
</dip>
<rake>
<value>106</value>
</rake>
</nodalPlane1>
<nodalPlane2>
<strike>
<value>159</value>
</strike>
<dip>
<value>76</value>
</dip>
<rake>
<value>86</value>
</rake>
</nodalPlane2>
</nodalPlanes>
<principalAxes>
<tAxis>
<azimuth>
<value>63</value>
</azimuth>
<plunge>
<value>59</value>
</plunge>
<length>
<value>1903000000000000000000</value>
</length>
</tAxis>
<pAxis>
<azimuth>
<value>252</value>
</azimuth>
<plunge>
<value>30</value>
</plunge>
<length>
<value>-1892000000000000000000</value>
</length>
</pAxis>
<nAxis>
<azimuth>
<value>160</value>
</azimuth>
<plunge>
<value>4</value>
</plunge>
<length>
<value>-12000000000000000000</value>
</length>
</nAxis>
</principalAxes>
<creationInfo>
<agencyID>GCMT</agencyID>
<version>V10</version>
</creationInfo>
</focalMechanism>
<magnitude publicID="smi:ds.iris.edu/spudservice/momenttensor/gcmtid/C201404012346A/quakeml#magnitude">
<mag>
<value>8.1</value>
</mag>
<type>Mwc</type>
</magnitude>
<origin publicID="smi:www.iris.edu/spudservice/momenttensor/gcmtid/C201404012346A#reforigin">
<time>
<value>2014-04-01T23:46:47.300Z</value>
</time>
<longitude>
<value>-70.77</value>
</longitude>
<latitude>
<value>-19.61</value>
</latitude>
<depth>
<value>25000.0</value>
</depth>
</origin>
<origin publicID="smi:www.iris.edu/spudservice/momenttensor/gcmtid/C201404012346A#cmtorigin">
<time>
<value>2014-04-01T23:47:31.500Z</value>
</time>
<longitude>
<value>-70.81</value>
</longitude>
<latitude>
<value>-19.7</value>
</latitude>
<depth>
<value>21600.0</value>
</depth>
<timeFixed>false</timeFixed>
<epicenterFixed>false</epicenterFixed>
</origin>
<preferredOriginID>smi:www.iris.edu/spudservice/momenttensor/gcmtid/C201404012346A#cmtorigin</preferredOriginID>
<preferredFocalMechanismID>smi:ds.iris.edu/spudservice/momenttensor/gcmtid/C201404012346A/quakeml#focalmechanism</preferredFocalMechanismID>
<type>earthquake</type>
</event>
<creationInfo>
<creationTime>2015-05-07T09:58:12.311Z</creationTime>
<version>V10</version>
</creationInfo>
</eventParameters>
</q:quakeml>

BIN
ObsPy/data/IV_BOB.mseed (Stored with Git LFS) Normal file
View File

Binary file not shown.

1
ObsPy/data/IV_BOB.xml Normal file
View File

File diff suppressed because one or more lines are too long

21877
ObsPy/data/all_stations.xml Normal file
View File

File diff suppressed because it is too large Load Diff

36
ObsPy/data/event_tohoku_mainshock.xml Normal file
View File

@ -0,0 +1,36 @@
<?xml version="1.0" encoding="UTF-8"?><q:quakeml xmlns:q="http://quakeml.org/xmlns/quakeml/1.2" xmlns:iris="http://service.iris.edu/fdsnws/event/1/" xmlns="http://quakeml.org/xmlns/bed/1.2" xmlns:xsi="http://www.w3.org/2000/10/XMLSchema-instance" xsi:schemaLocation="http://quakeml.org/schema/xsd http://quakeml.org/schema/xsd/QuakeML-1.2.xsd"><eventParameters publicID="smi:service.iris.edu/fdsnws/event/1/query"><event publicID="smi:service.iris.edu/fdsnws/event/1/query?eventid=3279407">
<type>earthquake</type>
<description xmlns:iris="http://service.iris.edu/fdsnws/event/1/" iris:FEcode="228">
<type>Flinn-Engdahl region</type>
<text>NEAR EAST COAST OF HONSHU, JAPAN</text>
</description>
<preferredMagnitudeID>smi:service.iris.edu/fdsnws/event/1/query?magnitudeid=16642444</preferredMagnitudeID>
<preferredOriginID>smi:service.iris.edu/fdsnws/event/1/query?originid=9933375</preferredOriginID>
<origin xmlns:iris="http://service.iris.edu/fdsnws/event/1/" publicID="smi:service.iris.edu/fdsnws/event/1/query?originid=9933375" iris:contributor="ISC" iris:catalog="ISC" iris:contributorOriginId="02227159" iris:contributorEventId="16461282">
<time>
<value>2011-03-11T05:46:23.2000</value>
</time>
<creationInfo>
<author>ISC</author>
</creationInfo>
<latitude>
<value>38.2963</value>
</latitude>
<longitude>
<value>142.498</value>
</longitude>
<depth>
<value>19700.0</value>
</depth>
</origin>
<magnitude publicID="smi:service.iris.edu/fdsnws/event/1/query?magnitudeid=16642444">
<originID>smi:service.iris.edu/fdsnws/event/1/query?originid=9933383</originID>
<type>MW</type>
<mag>
<value>9.1</value>
</mag>
<creationInfo>
<author>GCMT</author>
</creationInfo>
</magnitude>
</event></eventParameters></q:quakeml>

184
ObsPy/data/event_tohoku_with_big_aftershocks.xml Normal file
View File

@ -0,0 +1,184 @@
<?xml version="1.0" encoding="UTF-8"?><q:quakeml xmlns:q="http://quakeml.org/xmlns/quakeml/1.2" xmlns:iris="http://service.iris.edu/fdsnws/event/1/" xmlns="http://quakeml.org/xmlns/bed/1.2" xmlns:xsi="http://www.w3.org/2000/10/XMLSchema-instance" xsi:schemaLocation="http://quakeml.org/schema/xsd http://quakeml.org/schema/xsd/QuakeML-1.2.xsd">
<eventParameters publicID="smi:service.iris.edu/fdsnws/event/1/query">
<event publicID="smi:service.iris.edu/fdsnws/event/1/query?eventid=3279407">
<type>earthquake</type>
<description xmlns:iris="http://service.iris.edu/fdsnws/event/1/" iris:FEcode="228">
<type>Flinn-Engdahl region</type>
<text>NEAR EAST COAST OF HONSHU, JAPAN</text>
</description>
<preferredMagnitudeID>smi:service.iris.edu/fdsnws/event/1/query?magnitudeid=16642444</preferredMagnitudeID>
<preferredOriginID>smi:service.iris.edu/fdsnws/event/1/query?originid=9933375</preferredOriginID>
<origin xmlns:iris="http://service.iris.edu/fdsnws/event/1/" publicID="smi:service.iris.edu/fdsnws/event/1/query?originid=9933375" iris:contributor="ISC" iris:catalog="ISC" iris:contributorOriginId="02227159" iris:contributorEventId="16461282">
<time>
<value>2011-03-11T05:46:23.2000</value>
</time>
<creationInfo>
<author>ISC</author>
</creationInfo>
<latitude>
<value>38.2963</value>
</latitude>
<longitude>
<value>142.498</value>
</longitude>
<depth>
<value>19700.0</value>
</depth>
</origin>
<magnitude publicID="smi:service.iris.edu/fdsnws/event/1/query?magnitudeid=16642444">
<originID>smi:service.iris.edu/fdsnws/event/1/query?originid=9933383</originID>
<type>MW</type>
<mag>
<value>9.1</value>
</mag>
<creationInfo>
<author>GCMT</author>
</creationInfo>
</magnitude>
</event>
<event publicID="smi:service.iris.edu/fdsnws/event/1/query?eventid=3279489">
<type>earthquake</type>
<description xmlns:iris="http://service.iris.edu/fdsnws/event/1/" iris:FEcode="228">
<type>Flinn-Engdahl region</type>
<text>NEAR EAST COAST OF HONSHU, JAPAN</text>
</description>
<preferredMagnitudeID>smi:service.iris.edu/fdsnws/event/1/query?magnitudeid=16642643</preferredMagnitudeID>
<preferredOriginID>smi:service.iris.edu/fdsnws/event/1/query?originid=9933511</preferredOriginID>
<origin xmlns:iris="http://service.iris.edu/fdsnws/event/1/" publicID="smi:service.iris.edu/fdsnws/event/1/query?originid=9933511" iris:contributor="ISC" iris:catalog="ISC" iris:contributorOriginId="02227176" iris:contributorEventId="2707261">
<time>
<value>2011-03-11T06:15:37.5700</value>
</time>
<creationInfo>
<author>ISC</author>
</creationInfo>
<latitude>
<value>36.2274</value>
</latitude>
<longitude>
<value>141.088</value>
</longitude>
<depth>
<value>25400.0</value>
</depth>
</origin>
<magnitude publicID="smi:service.iris.edu/fdsnws/event/1/query?magnitudeid=16642643">
<originID>smi:service.iris.edu/fdsnws/event/1/query?originid=9933514</originID>
<type>MW</type>
<mag>
<value>7.9</value>
</mag>
<creationInfo>
<author>GCMT</author>
</creationInfo>
</magnitude>
</event>
<event publicID="smi:service.iris.edu/fdsnws/event/1/query?eventid=3279409">
<type>earthquake</type>
<description xmlns:iris="http://service.iris.edu/fdsnws/event/1/" iris:FEcode="229">
<type>Flinn-Engdahl region</type>
<text>OFF EAST COAST OF HONSHU, JAPAN</text>
</description>
<preferredMagnitudeID>smi:service.iris.edu/fdsnws/event/1/query?magnitudeid=16642722</preferredMagnitudeID>
<preferredOriginID>smi:service.iris.edu/fdsnws/event/1/query?originid=9933564</preferredOriginID>
<origin xmlns:iris="http://service.iris.edu/fdsnws/event/1/" publicID="smi:service.iris.edu/fdsnws/event/1/query?originid=9933564" iris:contributor="ISC" iris:catalog="ISC" iris:contributorOriginId="02227183" iris:contributorEventId="16304073">
<time>
<value>2011-03-11T06:25:50.7400</value>
</time>
<creationInfo>
<author>ISC</author>
</creationInfo>
<latitude>
<value>38.051</value>
</latitude>
<longitude>
<value>144.6297</value>
</longitude>
<depth>
<value>19800.0</value>
</depth>
</origin>
<magnitude publicID="smi:service.iris.edu/fdsnws/event/1/query?magnitudeid=16642722">
<originID>smi:service.iris.edu/fdsnws/event/1/query?originid=9933568</originID>
<type>MW</type>
<mag>
<value>7.6</value>
</mag>
<creationInfo>
<author>GCMT</author>
</creationInfo>
</magnitude>
</event>
<event publicID="smi:service.iris.edu/fdsnws/event/1/query?eventid=3282643">
<type>earthquake</type>
<description xmlns:iris="http://service.iris.edu/fdsnws/event/1/" iris:FEcode="228">
<type>Flinn-Engdahl region</type>
<text>NEAR EAST COAST OF HONSHU, JAPAN</text>
</description>
<preferredMagnitudeID>smi:service.iris.edu/fdsnws/event/1/query?magnitudeid=16176888</preferredMagnitudeID>
<preferredOriginID>smi:service.iris.edu/fdsnws/event/1/query?originid=9535390</preferredOriginID>
<origin xmlns:iris="http://service.iris.edu/fdsnws/event/1/" publicID="smi:service.iris.edu/fdsnws/event/1/query?originid=9535390" iris:contributor="NEIC PDE-M" iris:catalog="NEIC PDE">
<time>
<value>2011-04-07T14:32:43.2900</value>
</time>
<creationInfo>
<author>NEIC</author>
</creationInfo>
<latitude>
<value>38.276</value>
</latitude>
<longitude>
<value>141.588</value>
</longitude>
<depth>
<value>42000.0</value>
</depth>
</origin>
<magnitude publicID="smi:service.iris.edu/fdsnws/event/1/query?magnitudeid=16176888">
<originID>smi:service.iris.edu/fdsnws/event/1/query?originid=9535390</originID>
<type>MW</type>
<mag>
<value>7.1</value>
</mag>
<creationInfo>
<author>UCMT</author>
</creationInfo>
</magnitude>
</event>
<event publicID="smi:service.iris.edu/fdsnws/event/1/query?eventid=3317566">
<type>earthquake</type>
<description xmlns:iris="http://service.iris.edu/fdsnws/event/1/" iris:FEcode="229">
<type>Flinn-Engdahl region</type>
<text>OFF EAST COAST OF HONSHU, JAPAN</text>
</description>
<preferredMagnitudeID>smi:service.iris.edu/fdsnws/event/1/query?magnitudeid=16183694</preferredMagnitudeID>
<preferredOriginID>smi:service.iris.edu/fdsnws/event/1/query?originid=9540887</preferredOriginID>
<origin xmlns:iris="http://service.iris.edu/fdsnws/event/1/" publicID="smi:service.iris.edu/fdsnws/event/1/query?originid=9540887" iris:contributor="NEIC PDE-M" iris:catalog="NEIC PDE">
<time>
<value>2011-07-10T00:57:10.8000</value>
</time>
<creationInfo>
<author>NEIC</author>
</creationInfo>
<latitude>
<value>38.034</value>
</latitude>
<longitude>
<value>143.264</value>
</longitude>
<depth>
<value>23000.0</value>
</depth>
</origin>
<magnitude publicID="smi:service.iris.edu/fdsnws/event/1/query?magnitudeid=16183694">
<originID>smi:service.iris.edu/fdsnws/event/1/query?originid=9540887</originID>
<type>MW</type>
<mag>
<value>7.0</value>
</mag>
<creationInfo>
<author>UCMT</author>
</creationInfo>
</magnitude>
</event>
</eventParameters>
</q:quakeml>

BIN
ObsPy/data/example.mseed (Stored with Git LFS) Normal file
View File

Binary file not shown.

BIN
ObsPy/data/example.sac (Stored with Git LFS) Normal file
View File

Binary file not shown.

BIN
ObsPy/data/mtcarmel.mseed (Stored with Git LFS) Normal file
View File

Binary file not shown.

BIN
ObsPy/data/mtcarmel_100hz.mseed (Stored with Git LFS) Normal file
View File

Binary file not shown.

587
ObsPy/data/station_BFO.xml Normal file
View File

@ -0,0 +1,587 @@
<?xml version="1.0" encoding="UTF-8"?>
<FDSNStationXML xmlns="http://www.fdsn.org/xml/station/1" schemaVersion="1.0">
<Source>SeisComP3</Source>
<Sender>ODC</Sender>
<Created>2015-02-20T10:41:48</Created>
<Network code="GR" startDate="1976-02-17T00:00:00" restrictedStatus="open">
<Description>German Regional Seismic Network, BGR Hannover</Description>
<Station code="BFO" startDate="1991-01-01T00:00:00" restrictedStatus="open">
<Latitude>48.3311</Latitude>
<Longitude>8.3303</Longitude>
<Elevation>589</Elevation>
<Site>
<Name>GRSN Station Black Forest Observatory</Name>
<Country>Germany</Country>
</Site>
<CreationDate>1991-01-01T00:00:00</CreationDate>
<Channel code="BHE" startDate="1991-01-01T00:00:00" endDate="2011-10-19T00:00:00" restrictedStatus="open" locationCode="">
<Latitude>48.3311</Latitude>
<Longitude>8.3303</Longitude>
<Elevation>589</Elevation>
<Depth>0</Depth>
<Azimuth>90</Azimuth>
<Dip>0</Dip>
<SampleRate>20</SampleRate>
<SampleRateRatio>
<NumberSamples>20</NumberSamples>
<NumberSeconds>1</NumberSeconds>
</SampleRateRatio>
<StorageFormat>Steim2</StorageFormat>
<ClockDrift>0.05</ClockDrift>
<Sensor resourceId="Sensor#20120723141113.907112.4138">
<Description>STS2</Description>
</Sensor>
<DataLogger resourceId="Datalogger#20120723141113.883481.3977"/>
<Response>
<InstrumentSensitivity>
<Value>598802400</Value>
<Frequency>1</Frequency>
<InputUnits>
<Name>m/s</Name>
</InputUnits>
<OutputUnits>
<Name/>
</OutputUnits>
</InstrumentSensitivity>
<Stage number="1">
<PolesZeros resourceId="ResponsePAZ#20120723141113.90714.4139" name="GR.BFO..BHE_1991.001_STS2">
<InputUnits>
<Name>M/S</Name>
</InputUnits>
<OutputUnits>
<Name>V</Name>
</OutputUnits>
<PzTransferFunctionType>LAPLACE (RADIANS/SECOND)</PzTransferFunctionType>
<NormalizationFactor>1</NormalizationFactor>
<NormalizationFrequency>1</NormalizationFrequency>
<Pole number="0">
<Real>-0.0367429</Real>
<Imaginary>0.036754</Imaginary>
</Pole>
<Pole number="1">
<Real>-0.0367429</Real>
<Imaginary>-0.036754</Imaginary>
</Pole>
<Zero number="2">
<Real>0</Real>
<Imaginary>0</Imaginary>
</Zero>
<Zero number="3">
<Real>0</Real>
<Imaginary>0</Imaginary>
</Zero>
</PolesZeros>
<StageGain>
<Value>598802400</Value>
<Frequency>1</Frequency>
</StageGain>
</Stage>
<Stage number="2">
<Coefficients>
<InputUnits>
<Name>V</Name>
</InputUnits>
<OutputUnits>
<Name>COUNTS</Name>
</OutputUnits>
<CfTransferFunctionType>DIGITAL</CfTransferFunctionType>
</Coefficients>
<Decimation>
<InputSampleRate>6.930355364e-310</InputSampleRate>
<Factor>1</Factor>
<Offset>0</Offset>
<Delay>0</Delay>
<Correction>0</Correction>
</Decimation>
<StageGain>
<Value>1</Value>
<Frequency>0</Frequency>
</StageGain>
</Stage>
</Response>
</Channel>
<Channel code="BHN" startDate="1991-01-01T00:00:00" endDate="2011-10-19T00:00:00" restrictedStatus="open" locationCode="">
<Latitude>48.3311</Latitude>
<Longitude>8.3303</Longitude>
<Elevation>589</Elevation>
<Depth>0</Depth>
<Azimuth>0</Azimuth>
<Dip>0</Dip>
<SampleRate>20</SampleRate>
<SampleRateRatio>
<NumberSamples>20</NumberSamples>
<NumberSeconds>1</NumberSeconds>
</SampleRateRatio>
<StorageFormat>Steim2</StorageFormat>
<ClockDrift>0.05</ClockDrift>
<Sensor resourceId="Sensor#20120723141113.906463.4134">
<Description>STS2</Description>
</Sensor>
<DataLogger resourceId="Datalogger#20120723141113.883481.3977"/>
<Response>
<InstrumentSensitivity>
<Value>598802400</Value>
<Frequency>1</Frequency>
<InputUnits>
<Name>m/s</Name>
</InputUnits>
<OutputUnits>
<Name/>
</OutputUnits>
</InstrumentSensitivity>
<Stage number="1">
<PolesZeros resourceId="ResponsePAZ#20120723141113.906491.4135" name="GR.BFO..BHN_1991.001_STS2">
<InputUnits>
<Name>M/S</Name>
</InputUnits>
<OutputUnits>
<Name>V</Name>
</OutputUnits>
<PzTransferFunctionType>LAPLACE (RADIANS/SECOND)</PzTransferFunctionType>
<NormalizationFactor>1</NormalizationFactor>
<NormalizationFrequency>1</NormalizationFrequency>
<Pole number="0">
<Real>-0.0367429</Real>
<Imaginary>0.036754</Imaginary>
</Pole>
<Pole number="1">
<Real>-0.0367429</Real>
<Imaginary>-0.036754</Imaginary>
</Pole>
<Zero number="2">
<Real>0</Real>
<Imaginary>0</Imaginary>
</Zero>
<Zero number="3">
<Real>0</Real>
<Imaginary>0</Imaginary>
</Zero>
</PolesZeros>
<StageGain>
<Value>598802400</Value>
<Frequency>1</Frequency>
</StageGain>
</Stage>
<Stage number="2">
<Coefficients>
<InputUnits>
<Name>V</Name>
</InputUnits>
<OutputUnits>
<Name>COUNTS</Name>
</OutputUnits>
<CfTransferFunctionType>DIGITAL</CfTransferFunctionType>
</Coefficients>
<Decimation>
<InputSampleRate>6.930355389e-310</InputSampleRate>
<Factor>1</Factor>
<Offset>0</Offset>
<Delay>0</Delay>
<Correction>0</Correction>
</Decimation>
<StageGain>
<Value>1</Value>
<Frequency>0</Frequency>
</StageGain>
</Stage>
</Response>
</Channel>
<Channel code="BHZ" startDate="1991-01-01T00:00:00" endDate="2011-10-19T00:00:00" restrictedStatus="open" locationCode="">
<Latitude>48.3311</Latitude>
<Longitude>8.3303</Longitude>
<Elevation>589</Elevation>
<Depth>0</Depth>
<Azimuth>0</Azimuth>
<Dip>-90</Dip>
<SampleRate>20</SampleRate>
<SampleRateRatio>
<NumberSamples>20</NumberSamples>
<NumberSeconds>1</NumberSeconds>
</SampleRateRatio>
<StorageFormat>Steim2</StorageFormat>
<ClockDrift>0.05</ClockDrift>
<Sensor resourceId="Sensor#20120723141113.905894.4130">
<Description>STS2</Description>
</Sensor>
<DataLogger resourceId="Datalogger#20120723141113.883481.3977"/>
<Response>
<InstrumentSensitivity>
<Value>598802400</Value>
<Frequency>1</Frequency>
<InputUnits>
<Name>m/s</Name>
</InputUnits>
<OutputUnits>
<Name/>
</OutputUnits>
</InstrumentSensitivity>
<Stage number="1">
<PolesZeros resourceId="ResponsePAZ#20120723141113.905922.4131" name="GR.BFO..BHZ_1991.001_STS2">
<InputUnits>
<Name>M/S</Name>
</InputUnits>
<OutputUnits>
<Name>V</Name>
</OutputUnits>
<PzTransferFunctionType>LAPLACE (RADIANS/SECOND)</PzTransferFunctionType>
<NormalizationFactor>1</NormalizationFactor>
<NormalizationFrequency>1</NormalizationFrequency>
<Pole number="0">
<Real>-0.0367429</Real>
<Imaginary>0.036754</Imaginary>
</Pole>
<Pole number="1">
<Real>-0.0367429</Real>
<Imaginary>-0.036754</Imaginary>
</Pole>
<Zero number="2">
<Real>0</Real>
<Imaginary>0</Imaginary>
</Zero>
<Zero number="3">
<Real>0</Real>
<Imaginary>0</Imaginary>
</Zero>
</PolesZeros>
<StageGain>
<Value>598802400</Value>
<Frequency>1</Frequency>
</StageGain>
</Stage>
<Stage number="2">
<Coefficients>
<InputUnits>
<Name>V</Name>
</InputUnits>
<OutputUnits>
<Name>COUNTS</Name>
</OutputUnits>
<CfTransferFunctionType>DIGITAL</CfTransferFunctionType>
</Coefficients>
<Decimation>
<InputSampleRate>6.930355466e-310</InputSampleRate>
<Factor>1</Factor>
<Offset>0</Offset>
<Delay>0</Delay>
<Correction>0</Correction>
</Decimation>
<StageGain>
<Value>1</Value>
<Frequency>0</Frequency>
</StageGain>
</Stage>
</Response>
</Channel>
</Station>
</Network>
<Network code="II" startDate="1980-01-01T00:00:00" restrictedStatus="open">
<Description>(GSN) Global Seismograph Network (IRIS/IDA)</Description>
<Station code="BFO" startDate="1996-05-29T00:00:00" restrictedStatus="open">
<Latitude>48.3319</Latitude>
<Longitude>8.3311</Longitude>
<Elevation>589</Elevation>
<Site>
<Name>Black Forest Observatory, Schiltach, Germany</Name>
<Country> Germany</Country>
</Site>
<CreationDate>1996-05-29T00:00:00</CreationDate>
<Channel code="BHE" startDate="2008-01-16T00:00:00" endDate="2012-06-21T14:29:59" restrictedStatus="open" locationCode="00">
<Latitude>48.3319</Latitude>
<Longitude>8.3311</Longitude>
<Elevation>589</Elevation>
<Depth>0</Depth>
<Azimuth>90</Azimuth>
<Dip>0</Dip>
<SampleRate>20</SampleRate>
<SampleRateRatio>
<NumberSamples>20</NumberSamples>
<NumberSeconds>1</NumberSeconds>
</SampleRateRatio>
<StorageFormat>Steim1</StorageFormat>
<ClockDrift>2.5e-06</ClockDrift>
<Sensor resourceId="Sensor#20121207101201.213013.1484">
<Type>VBB</Type>
<Description>STS-1H</Description>
<Manufacturer>Streckeisen</Manufacturer>
<Model>STS-1H</Model>
</Sensor>
<DataLogger resourceId="Datalogger#20121207101201.210619.1481">
<Description>BFO.2008.016.H00</Description>
</DataLogger>
<Response>
<InstrumentSensitivity>
<Value>3.1632e+10</Value>
<Frequency>0.05</Frequency>
<InputUnits>
<Name>M/S</Name>
</InputUnits>
<OutputUnits>
<Name/>
</OutputUnits>
</InstrumentSensitivity>
<Stage number="1">
<PolesZeros resourceId="ResponsePAZ#20121207101201.213237.1485" name="BFO.2008.016.HE00">
<InputUnits>
<Name>M/S</Name>
</InputUnits>
<OutputUnits>
<Name>V</Name>
</OutputUnits>
<PzTransferFunctionType>LAPLACE (HERTZ)</PzTransferFunctionType>
<NormalizationFactor>100.06</NormalizationFactor>
<NormalizationFrequency>0.05</NormalizationFrequency>
<Pole number="0">
<Real>-6.28</Real>
<Imaginary>7.78213</Imaginary>
</Pole>
<Pole number="1">
<Real>-6.28</Real>
<Imaginary>-7.78213</Imaginary>
</Pole>
<Pole number="2">
<Real>-0.002142</Real>
<Imaginary>0.001753</Imaginary>
</Pole>
<Pole number="3">
<Real>-0.002142</Real>
<Imaginary>-0.001753</Imaginary>
</Pole>
<Zero number="4">
<Real>0</Real>
<Imaginary>0</Imaginary>
</Zero>
<Zero number="5">
<Real>0</Real>
<Imaginary>0</Imaginary>
</Zero>
</PolesZeros>
<StageGain>
<Value>2539.63</Value>
<Frequency>0.05</Frequency>
</StageGain>
</Stage>
<Stage number="2">
<Coefficients>
<InputUnits>
<Name>V</Name>
</InputUnits>
<OutputUnits>
<Name>COUNTS</Name>
</OutputUnits>
<CfTransferFunctionType>DIGITAL</CfTransferFunctionType>
</Coefficients>
<Decimation>
<InputSampleRate>0</InputSampleRate>
<Factor>1</Factor>
<Offset>0</Offset>
<Delay>0</Delay>
<Correction>0</Correction>
</Decimation>
<StageGain>
<Value>700143</Value>
<Frequency>0</Frequency>
</StageGain>
</Stage>
</Response>
</Channel>
<Channel code="BHN" startDate="2008-01-16T00:00:00" endDate="2012-06-21T14:29:59" restrictedStatus="open" locationCode="00">
<Latitude>48.3319</Latitude>
<Longitude>8.3311</Longitude>
<Elevation>589</Elevation>
<Depth>0</Depth>
<Azimuth>0</Azimuth>
<Dip>0</Dip>
<SampleRate>20</SampleRate>
<SampleRateRatio>
<NumberSamples>20</NumberSamples>
<NumberSeconds>1</NumberSeconds>
</SampleRateRatio>
<StorageFormat>Steim1</StorageFormat>
<ClockDrift>2.5e-06</ClockDrift>
<Sensor resourceId="Sensor#20121207101201.23289.1507">
<Type>VBB</Type>
<Description>STS-1H</Description>
<Manufacturer>Streckeisen</Manufacturer>
<Model>STS-1H</Model>
</Sensor>
<DataLogger resourceId="Datalogger#20121207101201.210619.1481">
<Description>BFO.2008.016.H00</Description>
</DataLogger>
<Response>
<InstrumentSensitivity>
<Value>2.39545e+10</Value>
<Frequency>0.05</Frequency>
<InputUnits>
<Name>M/S</Name>
</InputUnits>
<OutputUnits>
<Name/>
</OutputUnits>
</InstrumentSensitivity>
<Stage number="1">
<PolesZeros resourceId="ResponsePAZ#20121207101201.233133.1508" name="BFO.2008.016.HN00">
<InputUnits>
<Name>M/S</Name>
</InputUnits>
<OutputUnits>
<Name>V</Name>
</OutputUnits>
<PzTransferFunctionType>LAPLACE (HERTZ)</PzTransferFunctionType>
<NormalizationFactor>100.076</NormalizationFactor>
<NormalizationFrequency>0.05</NormalizationFrequency>
<Pole number="0">
<Real>-6.29</Real>
<Imaginary>7.77405</Imaginary>
</Pole>
<Pole number="1">
<Real>-6.29</Real>
<Imaginary>-7.77405</Imaginary>
</Pole>
<Pole number="2">
<Real>-0.002157</Real>
<Imaginary>0.001657</Imaginary>
</Pole>
<Pole number="3">
<Real>-0.002157</Real>
<Imaginary>-0.001657</Imaginary>
</Pole>
<Zero number="4">
<Real>0</Real>
<Imaginary>0</Imaginary>
</Zero>
<Zero number="5">
<Real>0</Real>
<Imaginary>0</Imaginary>
</Zero>
</PolesZeros>
<StageGain>
<Value>2572.99</Value>
<Frequency>0.05</Frequency>
</StageGain>
</Stage>
<Stage number="2">
<Coefficients>
<InputUnits>
<Name>V</Name>
</InputUnits>
<OutputUnits>
<Name>COUNTS</Name>
</OutputUnits>
<CfTransferFunctionType>DIGITAL</CfTransferFunctionType>
</Coefficients>
<Decimation>
<InputSampleRate>0</InputSampleRate>
<Factor>1</Factor>
<Offset>0</Offset>
<Delay>0</Delay>
<Correction>0</Correction>
</Decimation>
<StageGain>
<Value>700143</Value>
<Frequency>0</Frequency>
</StageGain>
</Stage>
</Response>
</Channel>
<Channel code="BHZ" startDate="2008-01-16T00:00:00" endDate="2012-06-21T14:29:59" restrictedStatus="open" locationCode="00">
<Latitude>48.3319</Latitude>
<Longitude>8.3311</Longitude>
<Elevation>589</Elevation>
<Depth>0</Depth>
<Azimuth>0</Azimuth>
<Dip>-90</Dip>
<SampleRate>20</SampleRate>
<SampleRateRatio>
<NumberSamples>20</NumberSamples>
<NumberSeconds>1</NumberSeconds>
</SampleRateRatio>
<StorageFormat>Steim1</StorageFormat>
<ClockDrift>2.5e-06</ClockDrift>
<Sensor resourceId="Sensor#20121207101201.250929.1527">
<Type>VBB</Type>
<Description>STS-1V</Description>
<Manufacturer>Streckeisen</Manufacturer>
<Model>STS-1V</Model>
</Sensor>
<DataLogger resourceId="Datalogger#20121207101201.210619.1481">
<Description>BFO.2008.016.H00</Description>
</DataLogger>
<Response>
<InstrumentSensitivity>
<Value>2.27009e+10</Value>
<Frequency>0.05</Frequency>
<InputUnits>
<Name>M/S</Name>
</InputUnits>
<OutputUnits>
<Name/>
</OutputUnits>
</InstrumentSensitivity>
<Stage number="1">
<PolesZeros resourceId="ResponsePAZ#20121207101201.251156.1528" name="BFO.2008.016.HZ00">
<InputUnits>
<Name>M/S</Name>
</InputUnits>
<OutputUnits>
<Name>V</Name>
</OutputUnits>
<PzTransferFunctionType>LAPLACE (HERTZ)</PzTransferFunctionType>
<NormalizationFactor>96.1326</NormalizationFactor>
<NormalizationFrequency>0.05</NormalizationFrequency>
<Pole number="0">
<Real>-6.27451</Real>
<Imaginary>7.53309</Imaginary>
</Pole>
<Pole number="1">
<Real>-6.27451</Real>
<Imaginary>-7.53309</Imaginary>
</Pole>
<Pole number="2">
<Real>-0.002003</Real>
<Imaginary>0.001898</Imaginary>
</Pole>
<Pole number="3">
<Real>-0.002003</Real>
<Imaginary>-0.001898</Imaginary>
</Pole>
<Zero number="4">
<Real>0</Real>
<Imaginary>0</Imaginary>
</Zero>
<Zero number="5">
<Real>0</Real>
<Imaginary>0</Imaginary>
</Zero>
</PolesZeros>
<StageGain>
<Value>2419.66</Value>
<Frequency>0.05</Frequency>
</StageGain>
</Stage>
<Stage number="2">
<Coefficients>
<InputUnits>
<Name>V</Name>
</InputUnits>
<OutputUnits>
<Name>COUNTS</Name>
</OutputUnits>
<CfTransferFunctionType>DIGITAL</CfTransferFunctionType>
</Coefficients>
<Decimation>
<InputSampleRate>0</InputSampleRate>
<Factor>1</Factor>
<Offset>0</Offset>
<Delay>0</Delay>
<Correction>0</Correction>
</Decimation>
<StageGain>
<Value>700143</Value>
<Frequency>0</Frequency>
</StageGain>
</Stage>
</Response>
</Channel>
</Station>
</Network>
</FDSNStationXML>

467
ObsPy/data/station_PFO.xml Normal file
View File

@ -0,0 +1,467 @@
<?xml version="1.0" encoding="ISO-8859-1"?>
<FDSNStationXML xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns="http://www.fdsn.org/xml/station/1" schemaVersion="1.0" xsi:schemaLocation="http://www.fdsn.org/xml/station/1 http://www.fdsn.org/xml/station/fdsn-station-1.0.xsd">
<Source>IRIS-DMC</Source>
<Sender>IRIS-DMC</Sender>
<Module>IRIS WEB SERVICE: fdsnws-station | version: 1.0.11</Module>
<ModuleURI>http://service.iris.edu/fdsnws/station/1/query?network=II&amp;level=response&amp;station=PFO&amp;starttime=2011-03-11T05%3A46%3A13.000000&amp;endtime=2011-03-11T05%3A46%3A33.000000&amp;channel=BHZ</ModuleURI>
<Created>2014-02-17T15:43:45</Created>
<Network code="II" startDate="1986-01-01T00:00:00" endDate="2020-12-12T23:59:59" restrictedStatus="open">
<Description>Global Seismograph Network (GSN - IRIS/IDA)</Description>
<TotalNumberStations>50</TotalNumberStations>
<SelectedNumberStations>1</SelectedNumberStations>
<Station code="PFO" startDate="2006-07-13T00:00:00" endDate="2020-12-31T23:59:59" restrictedStatus="open" alternateNetworkCodes=".EARTHSCOPE,_ANSS,.UNRESTRICTED,_CARIBE-EWS,_GSN-BROADBAND,_US-REF,_GSN,_STS-1,_ANSS-BB,_FDSN,_FDSN-ALL,_US-ALL,_REALTIME">
<Latitude>33.6107</Latitude>
<Longitude>-116.4555</Longitude>
<Elevation>1280.0</Elevation>
<Site>
<Name>Pinon Flat, California, USA</Name>
</Site>
<CreationDate>1986-10-24T00:00:00</CreationDate>
<TotalNumberChannels>250</TotalNumberChannels>
<SelectedNumberChannels>2</SelectedNumberChannels>
<Channel locationCode="00" startDate="2010-07-30T18:50:00" restrictedStatus="open" endDate="2012-07-02T03:59:59" code="BHZ">
<Comment>
<Value>S/N #119005</Value>
</Comment>
<Latitude>33.6107</Latitude>
<Longitude>-116.4555</Longitude>
<Elevation>1280.0</Elevation>
<Depth>5.3</Depth>
<Azimuth>0.0</Azimuth>
<Dip>-90.0</Dip>
<Type>CONTINUOUS</Type>
<Type>GEOPHYSICAL</Type>
<SampleRate>20.0</SampleRate>
<ClockDrift>2.5E-6</ClockDrift>
<CalibrationUnits>
<Name>V</Name>
<Description>Volts</Description>
</CalibrationUnits>
<Sensor>
<Type>Streckeisen STS-1 Seismometer with Metrozet E300</Type>
</Sensor>
<Response>
<InstrumentSensitivity>
<Value>5.24814E9</Value>
<Frequency>0.05</Frequency>
<InputUnits>
<Name>M/S</Name>
<Description>Velocity in Meters Per Second</Description>
</InputUnits>
<OutputUnits>
<Name>COUNTS</Name>
<Description>Digital Counts</Description>
</OutputUnits>
</InstrumentSensitivity>
<Stage number="1">
<PolesZeros>
<InputUnits>
<Name>M/S</Name>
<Description>Velocity in Meters Per Second</Description>
</InputUnits>
<OutputUnits>
<Name>V</Name>
<Description>Volts</Description>
</OutputUnits>
<PzTransferFunctionType>LAPLACE (HERTZ)</PzTransferFunctionType>
<NormalizationFactor>11.5758</NormalizationFactor>
<NormalizationFrequency>.05</NormalizationFrequency>
<Zero number="0">
<Real>0</Real>
<Imaginary>0</Imaginary>
</Zero>
<Zero number="1">
<Real>0</Real>
<Imaginary>0</Imaginary>
</Zero>
<Zero number="2">
<Real>-12.5</Real>
<Imaginary>0</Imaginary>
</Zero>
<Zero number="3">
<Real>-.0242718</Real>
<Imaginary>0</Imaginary>
</Zero>
<Zero number="4">
<Real>-.0242718</Real>
<Imaginary>0</Imaginary>
</Zero>
<Pole number="5">
<Real>-.0019211</Real>
<Imaginary>.00194895</Imaginary>
</Pole>
<Pole number="6">
<Real>-.0019211</Real>
<Imaginary>-.00194895</Imaginary>
</Pole>
<Pole number="7">
<Real>-.0242315</Real>
<Imaginary>.00153484</Imaginary>
</Pole>
<Pole number="8">
<Real>-.0242315</Real>
<Imaginary>-.00153484</Imaginary>
</Pole>
<Pole number="9">
<Real>-7.691</Real>
<Imaginary>9.25817</Imaginary>
</Pole>
<Pole number="10">
<Real>-7.691</Real>
<Imaginary>-9.25817</Imaginary>
</Pole>
</PolesZeros>
<StageGain>
<Value>3314.4</Value>
<Frequency>.05</Frequency>
</StageGain>
</Stage>
<Stage number="2">
<PolesZeros>
<InputUnits>
<Name>V</Name>
<Description>Volts</Description>
</InputUnits>
<OutputUnits>
<Name>V</Name>
<Description>Volts</Description>
</OutputUnits>
<PzTransferFunctionType>LAPLACE (HERTZ)</PzTransferFunctionType>
<NormalizationFactor>1</NormalizationFactor>
<NormalizationFrequency>.05</NormalizationFrequency>
</PolesZeros>
<StageGain>
<Value>1</Value>
<Frequency>.05</Frequency>
</StageGain>
</Stage>
<Stage number="3">
<Coefficients>
<InputUnits>
<Name>V</Name>
<Description>Volts</Description>
</InputUnits>
<OutputUnits>
<Name>COUNTS</Name>
<Description>Digital Counts</Description>
</OutputUnits>
<CfTransferFunctionType>DIGITAL</CfTransferFunctionType>
</Coefficients>
<Decimation>
<InputSampleRate>20</InputSampleRate>
<Factor>1</Factor>
<Offset>0</Offset>
<Delay>0.0</Delay>
<Correction>0.0</Correction>
</Decimation>
<StageGain>
<Value>1583330</Value>
<Frequency>0</Frequency>
</StageGain>
</Stage>
<Stage number="4">
<Coefficients>
<InputUnits>
<Name>COUNTS</Name>
<Description>Digital Counts</Description>
</InputUnits>
<OutputUnits>
<Name>COUNTS</Name>
<Description>Digital Counts</Description>
</OutputUnits>
<CfTransferFunctionType>DIGITAL</CfTransferFunctionType>
<Numerator>-.0000000000000000365342</Numerator>
<Numerator>.0000000367488</Numerator>
<Numerator>-.00000042706</Numerator>
<Numerator>.00000114502</Numerator>
<Numerator>-.000000187594</Numerator>
<Numerator>-.000000337274</Numerator>
<Numerator>.00000278747</Numerator>
<Numerator>-.00000374403</Numerator>
<Numerator>.00000541172</Numerator>
<Numerator>.00000747336</Numerator>
<Numerator>-.000517759</Numerator>
<Numerator>.000210677</Numerator>
<Numerator>.0000463258</Numerator>
<Numerator>-.000608222</Numerator>
<Numerator>.00144175</Numerator>
<Numerator>-.00240627</Numerator>
<Numerator>.00322534</Numerator>
<Numerator>-.00350639</Numerator>
<Numerator>.00281441</Numerator>
<Numerator>-.000771971</Numerator>
<Numerator>-.00280512</Numerator>
<Numerator>.00777805</Numerator>
<Numerator>-.0135815</Numerator>
<Numerator>.0191765</Numerator>
<Numerator>-.0229704</Numerator>
<Numerator>.0240398</Numerator>
<Numerator>-.0220986</Numerator>
<Numerator>.00860734</Numerator>
<Numerator>.0117525</Numerator>
<Numerator>-.0447787</Numerator>
<Numerator>.0964923</Numerator>
<Numerator>-.191755</Numerator>
<Numerator>.527652</Numerator>
<Numerator>.724167</Numerator>
<Numerator>-.156905</Numerator>
<Numerator>.0442574</Numerator>
<Numerator>.00314168</Numerator>
<Numerator>-.0266714</Numerator>
<Numerator>.0361532</Numerator>
<Numerator>-.0385687</Numerator>
<Numerator>.0310842</Numerator>
<Numerator>-.0235259</Numerator>
<Numerator>.0153211</Numerator>
<Numerator>-.00740398</Numerator>
<Numerator>.00109645</Numerator>
<Numerator>.00309797</Numerator>
<Numerator>-.0051932</Numerator>
<Numerator>.00556131</Numerator>
<Numerator>-.0047611</Numerator>
<Numerator>.00338213</Numerator>
<Numerator>-.00192052</Numerator>
<Numerator>.000715218</Numerator>
<Numerator>.0000767719</Numerator>
<Numerator>-.000451897</Numerator>
<Numerator>.0005027</Numerator>
<Numerator>-.000565037</Numerator>
<Numerator>-.00005568</Numerator>
<Numerator>.0000157736</Numerator>
<Numerator>-.00000141985</Numerator>
<Numerator>.000000814909</Numerator>
<Numerator>.000000680795</Numerator>
<Numerator>-.00000125273</Numerator>
<Numerator>.00000152435</Numerator>
<Numerator>-.000000283336</Numerator>
<Numerator>-.0000000106384</Numerator>
<Numerator>.00000000125712</Numerator>
<Numerator>-.0000000000542954</Numerator>
</Coefficients>
<Decimation>
<InputSampleRate>20</InputSampleRate>
<Factor>1</Factor>
<Offset>0</Offset>
<Delay>1.6305</Delay>
<Correction>1.6305</Correction>
</Decimation>
<StageGain>
<Value>1</Value>
<Frequency>0</Frequency>
</StageGain>
</Stage>
</Response>
</Channel>
<Channel locationCode="10" startDate="2010-07-23T00:00:00" restrictedStatus="open" endDate="2013-10-30T23:59:59" code="BHZ">
<Comment>
<Value>S/N #484</Value>
</Comment>
<Latitude>33.6107</Latitude>
<Longitude>-116.4555</Longitude>
<Elevation>1280.0</Elevation>
<Depth>5.3</Depth>
<Azimuth>0.0</Azimuth>
<Dip>-90.0</Dip>
<Type>CONTINUOUS</Type>
<Type>GEOPHYSICAL</Type>
<SampleRate>40.0</SampleRate>
<ClockDrift>1.25E-6</ClockDrift>
<CalibrationUnits>
<Name>V</Name>
<Description>Volts</Description>
</CalibrationUnits>
<Sensor>
<Type>Nanometrics Trillium 240 Seismometer</Type>
</Sensor>
<Response>
<InstrumentSensitivity>
<Value>2.00625E9</Value>
<Frequency>0.05</Frequency>
<InputUnits>
<Name>M/S</Name>
<Description>Velocity in Meters Per Second</Description>
</InputUnits>
<OutputUnits>
<Name>COUNTS</Name>
<Description>Digital Counts</Description>
</OutputUnits>
</InstrumentSensitivity>
<Stage number="1">
<PolesZeros>
<InputUnits>
<Name>M/S</Name>
<Description>Velocity in Meters Per Second</Description>
</InputUnits>
<OutputUnits>
<Name>V</Name>
<Description>Volts</Description>
</OutputUnits>
<PzTransferFunctionType>LAPLACE (HERTZ)</PzTransferFunctionType>
<NormalizationFactor>9667.55</NormalizationFactor>
<NormalizationFrequency>.05</NormalizationFrequency>
<Zero number="0">
<Real>0</Real>
<Imaginary>0</Imaginary>
</Zero>
<Zero number="1">
<Real>0</Real>
<Imaginary>0</Imaginary>
</Zero>
<Zero number="2">
<Real>-14.588</Real>
<Imaginary>0</Imaginary>
</Zero>
<Zero number="3">
<Real>-25.4807</Real>
<Imaginary>0</Imaginary>
</Zero>
<Zero number="4">
<Real>-510.408</Real>
<Imaginary>0</Imaginary>
</Zero>
<Pole number="5">
<Real>-.00281527</Real>
<Imaginary>.00280021</Imaginary>
</Pole>
<Pole number="6">
<Real>-.00281527</Real>
<Imaginary>-.00280021</Imaginary>
</Pole>
<Pole number="7">
<Real>-16.994</Real>
<Imaginary>0</Imaginary>
</Pole>
<Pole number="8">
<Real>-30.558</Real>
<Imaginary>41.237</Imaginary>
</Pole>
<Pole number="9">
<Real>-30.558</Real>
<Imaginary>-41.237</Imaginary>
</Pole>
<Pole number="10">
<Real>-88.76</Real>
<Imaginary>181.91</Imaginary>
</Pole>
<Pole number="11">
<Real>-88.76</Real>
<Imaginary>-181.91</Imaginary>
</Pole>
</PolesZeros>
<StageGain>
<Value>1200</Value>
<Frequency>.05</Frequency>
</StageGain>
</Stage>
<Stage number="2">
<PolesZeros>
<InputUnits>
<Name>V</Name>
<Description>Volts</Description>
</InputUnits>
<OutputUnits>
<Name>V</Name>
<Description>Volts</Description>
</OutputUnits>
<PzTransferFunctionType>LAPLACE (HERTZ)</PzTransferFunctionType>
<NormalizationFactor>1</NormalizationFactor>
<NormalizationFrequency>.05</NormalizationFrequency>
</PolesZeros>
<StageGain>
<Value>1</Value>
<Frequency>.05</Frequency>
</StageGain>
</Stage>
<Stage number="3">
<Coefficients>
<InputUnits>
<Name>V</Name>
<Description>Volts</Description>
</InputUnits>
<OutputUnits>
<Name>COUNTS</Name>
<Description>Digital Counts</Description>
</OutputUnits>
<CfTransferFunctionType>DIGITAL</CfTransferFunctionType>
</Coefficients>
<Decimation>
<InputSampleRate>40</InputSampleRate>
<Factor>1</Factor>
<Offset>0</Offset>
<Delay>0.0</Delay>
<Correction>0.0</Correction>
</Decimation>
<StageGain>
<Value>1671840</Value>
<Frequency>0</Frequency>
</StageGain>
</Stage>
<Stage number="4">
<Coefficients>
<InputUnits>
<Name>COUNTS</Name>
<Description>Digital Counts</Description>
</InputUnits>
<OutputUnits>
<Name>COUNTS</Name>
<Description>Digital Counts</Description>
</OutputUnits>
<CfTransferFunctionType>DIGITAL</CfTransferFunctionType>
<Numerator>.000000000000418952</Numerator>
<Numerator>.000330318</Numerator>
<Numerator>.00102921</Numerator>
<Numerator>-.00314123</Numerator>
<Numerator>.000205709</Numerator>
<Numerator>.00152521</Numerator>
<Numerator>-.00623193</Numerator>
<Numerator>.0104801</Numerator>
<Numerator>-.0131202</Numerator>
<Numerator>.0107821</Numerator>
<Numerator>-.00144455</Numerator>
<Numerator>-.0158729</Numerator>
<Numerator>.0395074</Numerator>
<Numerator>-.0651036</Numerator>
<Numerator>.0853716</Numerator>
<Numerator>-.0891913</Numerator>
<Numerator>.0500619</Numerator>
<Numerator>.837233</Numerator>
<Numerator>.266723</Numerator>
<Numerator>-.166693</Numerator>
<Numerator>.095284</Numerator>
<Numerator>-.0509218</Numerator>
<Numerator>.0161458</Numerator>
<Numerator>.00706362</Numerator>
<Numerator>-.0183877</Numerator>
<Numerator>.0199414</Numerator>
<Numerator>-.0154895</Numerator>
<Numerator>.00852735</Numerator>
<Numerator>-.00255789</Numerator>
<Numerator>-.00181103</Numerator>
<Numerator>.00242649</Numerator>
<Numerator>-.00375769</Numerator>
<Numerator>.000467293</Numerator>
<Numerator>.000633072</Numerator>
<Numerator>-.00000156874</Numerator>
<Numerator>-.000012548</Numerator>
<Numerator>.000000321041</Numerator>
<Numerator>-.0000000263324</Numerator>
<Numerator>-.0000000509997</Numerator>
</Coefficients>
<Decimation>
<InputSampleRate>40</InputSampleRate>
<Factor>1</Factor>
<Offset>0</Offset>
<Delay>0.43045</Delay>
<Correction>0.43045</Correction>
</Decimation>
<StageGain>
<Value>1</Value>
<Frequency>0</Frequency>
</StageGain>
</Stage>
</Response>
</Channel>
</Station>
</Network>
</FDSNStationXML>

66
ObsPy/data/tohoku.xml Normal file
View File

@ -0,0 +1,66 @@
<?xml version="1.0" encoding="US-ASCII" standalone="yes"?>
<q:quakeml xmlns:q="http://quakeml.org/xmlns/quakeml/1.2" xmlns="http://quakeml.org/xmlns/bed/1.2" xmlns:ingv="http://webservices.rm.ingv.it/fdsnws/event/1">
<eventParameters publicID="smi:webservices.rm.ingv.it/fdsnws/event/1/query">
<event publicID="smi:webservices.rm.ingv.it/fdsnws/event/1/query?eventId=2322919">
<type>earthquake</type>
<description>
<type>region name</type>
<text>Near east coast of eastern Honshu, Japan</text>
</description>
<preferredMagnitudeID>smi:webservices.rm.ingv.it/fdsnws/event/1/query?magnitudeId=2122949</preferredMagnitudeID>
<preferredOriginID>smi:webservices.rm.ingv.it/fdsnws/event/1/query?originId=1180759</preferredOriginID>
<creationInfo>
<agencyID>INGV</agencyID>
<author>CSEM</author>
<creationTime>2013-01-05T23:02:38</creationTime>
</creationInfo>
<origin publicID="smi:webservices.rm.ingv.it/fdsnws/event/1/query?originId=1180759">
<evaluationMode>manual</evaluationMode>
<type>hypocenter</type>
<time>
<value>2011-03-11T05:46:23.000000</value>
</time>
<latitude>
<value>38.3</value>
</latitude>
<longitude>
<value>142.5</value>
</longitude>
<depth>
<value>21000</value>
</depth>
<depthType>from location</depthType>
<quality>
<associatedPhaseCount>88</associatedPhaseCount>
<usedPhaseCount>65</usedPhaseCount>
<minimumDistance>0.00000</minimumDistance>
<maximumDistance>0.00000</maximumDistance>
<associatedStationCount>82</associatedStationCount>
<usedStationCount>65</usedStationCount>
</quality>
<evaluationStatus>reviewed</evaluationStatus>
<methodID>smi:webservices.rm.ingv.it/fdsnws/event/1/query?methodId=2</methodID>
<earthModelID>smi:webservices.rm.ingv.it/fdsnws/event/1/query?earthModelId=1</earthModelID>
<creationInfo>
<agencyID>INGV</agencyID>
<author>CSEM</author>
<version>1000</version>
<creationTime>2013-01-05T23:02:39</creationTime>
</creationInfo>
</origin>
<magnitude publicID="smi:webservices.rm.ingv.it/fdsnws/event/1/query?magnitudeId=2122949">
<originID>smi:webservices.rm.ingv.it/fdsnws/event/1/query?originId=1180759</originID>
<type>Mw</type>
<mag>
<value>9.0</value>
</mag>
<methodID>smi:webservices.rm.ingv.it/fdsnws/event/1/query?methodId=2</methodID>
<creationInfo>
<agencyID>INGV</agencyID>
<author>CSEM</author>
<creationTime>2013-01-05T23:02:39</creationTime>
</creationInfo>
</magnitude>
</event>
</eventParameters>
</q:quakeml>

BIN
ObsPy/data/waveform_BFO_BHE.sac (Stored with Git LFS) Normal file
View File

Binary file not shown.

BIN
ObsPy/data/waveform_BFO_BHN.sac (Stored with Git LFS) Normal file
View File

Binary file not shown.

BIN
ObsPy/data/waveform_BFO_BHZ.sac (Stored with Git LFS) Normal file
View File

Binary file not shown.

BIN
ObsPy/data/waveform_PFO.mseed (Stored with Git LFS) Normal file
View File

Binary file not shown.

BIN
ObsPy/data/waveform_PFO_synthetics.mseed (Stored with Git LFS) Normal file
View File

Binary file not shown.

3
ObsPy/images/Event.svg Normal file
View File

File diff suppressed because one or more lines are too long
Side by Side

After

Width:  |  Height:  |  Size: 27 KiB

3
ObsPy/images/Inventory.svg Normal file
View File

File diff suppressed because one or more lines are too long
Side by Side

After

Width:  |  Height:  |  Size: 26 KiB

3
ObsPy/images/Stream_Trace.svg Normal file
View File

File diff suppressed because one or more lines are too long
Side by Side

After

Width:  |  Height:  |  Size: 36 KiB

BIN
ObsPy/images/cos_taper.png (Stored with Git LFS) Normal file
View File

Binary file not shown.

BIN
ObsPy/images/notebook_toolbar.png (Stored with Git LFS) Normal file
View File

Binary file not shown.

BIN
ObsPy/images/obspy_logo_full_524x179px.png (Stored with Git LFS) Normal file
View File

Binary file not shown.

4
README.md
View File

@ -1,6 +1,6 @@
# Seismological Data Analysis 2022
# Seismological Data Analysis 2024
Jupyter Notebooks for class "Seismological Data Analysis" (summer term 2022)
Jupyter Notebooks for class "Seismological Data Analysis" (summer term 2024)
## License
The content of this repository is licensed under Creative Commons Attribution 4.0 International (CC-BY 4.0) ***unless otherwise stated*** (see below)

8
dataanalaysis.code-workspace Normal file
View File

@ -0,0 +1,8 @@
{
"folders": [
{
"path": "."
}
],
"settings": {}
}