Compare commits

..

5 Commits

Author SHA1 Message Date
3be655017f updated paths to header image file 2024-06-10 12:15:28 +02:00
228434397d add notebook collection from seismo-live.org on ObsPy
activate LFS
2024-06-10 11:47:40 +02:00
dbcb4a3d8a initial import of PPSD exercise from seismo live web site 2024-06-03 12:17:07 +02:00
d829006103 update README.md 2024-04-22 11:46:34 +02:00
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

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

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
}

Binary file not shown.

Binary file not shown.

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>

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
}

File diff suppressed because one or more lines are too long

File diff suppressed because one or more lines are too long

File diff suppressed because one or more lines are too long

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}

File diff suppressed because one or more lines are too long

1
ObsPy/06_FDSN.ipynb Normal file

File diff suppressed because one or more lines are too long

File diff suppressed because one or more lines are too long

File diff suppressed because one or more lines are too long

12315
ObsPy/data/2014.ndk Normal file

File diff suppressed because it is too large Load Diff

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

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

Binary file not shown.

1
ObsPy/data/IV_BOB.xml Normal file

File diff suppressed because one or more lines are too long

21877
ObsPy/data/all_stations.xml Normal file

File diff suppressed because it is too large Load Diff

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>

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

Binary file not shown.

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

Binary file not shown.

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

Binary file not shown.

BIN
ObsPy/data/mtcarmel_100hz.mseed (Stored with Git LFS) Normal 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

Binary file not shown.

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

Binary file not shown.

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

Binary file not shown.

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

Binary file not shown.

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

Binary file not shown.

3
ObsPy/images/Event.svg Normal file

File diff suppressed because one or more lines are too long

After

Width:  |  Height:  |  Size: 27 KiB

File diff suppressed because one or more lines are too long

After

Width:  |  Height:  |  Size: 26 KiB

File diff suppressed because one or more lines are too long

After

Width:  |  Height:  |  Size: 36 KiB

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

Binary file not shown.

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

Binary file not shown.

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

Binary file not shown.

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 ## License
The content of this repository is licensed under Creative Commons Attribution 4.0 International (CC-BY 4.0) ***unless otherwise stated*** (see below) The content of this repository is licensed under Creative Commons Attribution 4.0 International (CC-BY 4.0) ***unless otherwise stated*** (see below)

View File

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