das2numpy 0.0.1__py3-none-any.whl

das2numpy-0.0.1.dist-info/METADATA

@@ -0,0 +1,89 @@
+Metadata-Version: 2.1
+Name: das2numpy
+Version: 0.0.1
+Summary: A simple and universal package for loading large amounts of distributed acoustic sensing (DAS) data.
+Author-email: Erik Genthe <erik.genthe@desy.de>
+Project-URL: Homepage, https://git.physnet.uni-hamburg.de/wave/das2numpy
+Classifier: Programming Language :: Python :: 3
+Classifier: License :: OSI Approved :: GNU General Public License v3 (GPLv3)
+Classifier: Operating System :: OS Independent
+Requires-Python: >=3.8
+Description-Content-Type: text/markdown
+License-File: LICENSE
+
+# Module for loading Distributed Acoustic Sensing (DAS) data. SILIXA / OPTASENSE
+
+Python: If you want to get started quickly, have a look at the [examples](#example-python-script).
+
+
+## Install dependencies
+If you want to use pip for installing, you can just execute *install_dependencies.sh*.
+Otherwise, have a look into install_dependencies.sh and install the listed packages yourself.
+
+
+## Use as python module
+### API
+
+
+#### Recommended: simplest interface
+```python
+def load_array(t_start:datetime, t_end:datetime, channel_start:int, channel_end:int) -> NP.ndarray:
+```
+
+```
+Loads data and returns it as a numpy array. 
+Args:
+    t_start (datetime): datetime object which defines the start of the data to load.
+    t_end (datetime): datetime object which defines the end of the data to load.
+    channel_start (int): The starting index of sensor in the data (inclusive).
+    channel_end (int): The ending index of sensors in the data (exclusive).
+Returns:
+    A 2d-numpy-array containing the data.
+    The first axis corresponds to the time, the second, to the channel
+ ```
+
+
+#### More detailed interface
+```python
+def load_array(t_start:datetime, t_end:datetime, t_step:int, channel_start:int, channel_end:int, channel_step:int) -> NP.ndarray:
+```
+
+``` Loading data into numpy array.
+Returns nothing, the data can be accessed by accessing the data field of this instance.
+Warning: using a different value then 1 for t_step or channel_step can result in a high cpu-usage.
+        Consider using multithreaded=True in the constructor and a high amount of workers if needed.
+Constraints: 
+    t_start has to be less or equal t_end, 
+    same for channel_start and channel_end.
+    t_step and channel_step have to be greater then 0
+Args:
+    t_start (datetime): datetime object which defines the start of the data to load.
+    t_end (datetime): datetime object which defines the end of the data to load.
+    t_step (int): If you, for example only want to load the data of every fourth timestep use t_end=4
+    channel_start (int): The starting index of sensor in the data (inclusive).
+    channel_end (int): The ending index of sensors in the data (exclusive).
+    channel_step (int): Like t_step, but for the sensor position.
+Returns:
+    A 2d-numpy-array containing the data.
+    The first axis corresponds to the time, the second, to the channel
+```
+
+### Lower level interfaces
+There are also lower level interfaces in the module.
+For example, the above interfaces also exist with POSIX timestamps in milliseconds instead of datetime objects. These timestamps have exactly the same resolution as the time axis of the resulting array.
+
+
+### Example python-script
+
+[example.py](example.py)
+
+
+## Use as command line interface
+
+Example call (make sure that the current working directory is not inside idas2numpy):
+```python -m idas2numpy "SILIXA" ~/iDAS/work/2024-05-10-desy/ 2024-05-10T10:01:00 2024-05-10T10:02:00 10 0 1000 10 default
+```
+
+For more information:
+```python -m idas2numpy -h
+```

das2numpy-0.0.1.dist-info/RECORD

@@ -0,0 +1,16 @@
+example.py,sha256=PHPB62mrEP3AT0ykTxr89IOoXY_U2-ZUJ1Kovf2g9Qs,1175
+das2numpy/__init__.py,sha256=L9A_nHlQC_jpcQpTfFpHPlulpik0gr3jzZUbxETj8kM,785
+das2numpy/__main__.py,sha256=UhJ7jyS4iT39aWaNTspNHaROKaIkY7pPssVcqAwVmdE,2832
+das2numpy/chunk.py,sha256=FgtA1Sd8oYOydcFBoWMm8qg3Pg9yYhyopohG_6fgxas,10951
+das2numpy/filefinder.py,sha256=dOoN2nNpjKRpdV_1fWqFFgnwJWmphDgpENv2hcKsU9A,4466
+das2numpy/test.py,sha256=UZ93sWuPoIYxpjHW7FBzWUZP4wi9LKGIZ48NOKXrZvo,5812
+das2numpy/utils.py,sha256=07wnRfwPvZFmZftxcsf4utTfBMIyAQFox1MzH4i0jU4,4512
+das2numpy/setups/light_tdms_reader.py,sha256=0rWl56Sb-4h_6j5CmpJqUtPrW53nhdEE4N02mpqRVIQ,20909
+das2numpy/setups/optasense_b35idefix.py,sha256=MxYxSBNy6vZMMu6Is53o0xy6d_STatPqIKlhV2ukMwI,3329
+das2numpy/setups/optasense_b35idefix_fast.py,sha256=qanSrfXRnE6txdA8XAah-UUE819t2td97rFaBt6z31k,4087
+das2numpy/setups/silixa.py,sha256=MH1rC6DxsKt3qHFJlnr79rpoaEFKukCleS2ll_jthJI,2679
+das2numpy-0.0.1.dist-info/LICENSE,sha256=va7uzkBxuvq1TlSHpdrjOEqKO8RAciSGbcR9Z55LoyU,35132
+das2numpy-0.0.1.dist-info/METADATA,sha256=bn4euK2-BhM0KveqzQsu4YeCgj69jE5-Ji_OTtBlN1E,3538
+das2numpy-0.0.1.dist-info/WHEEL,sha256=Z4pYXqR_rTB7OWNDYFOm1qRk0RX6GFP2o8LgvP453Hk,91
+das2numpy-0.0.1.dist-info/top_level.txt,sha256=uJ4_X_a22JACdRAOHOcCqCtm8ciFMlWZ9om2zFQi9h4,18
+das2numpy-0.0.1.dist-info/RECORD,,

das2numpy-0.0.1.dist-info/WHEEL

@@ -0,0 +1,5 @@
+Wheel-Version: 1.0
+Generator: setuptools (70.3.0)
+Root-Is-Purelib: true
+Tag: py3-none-any
+

das2numpy-0.0.1.dist-info/top_level.txt

@@ -0,0 +1,2 @@
+das2numpy
+example

example.py

@@ -0,0 +1,39 @@
+import sys
+from datetime import datetime
+import matplotlib.pyplot as PP
+from das2numpy import loader, utils
+
+
+print("Load data to numpy-array")
+t_start = datetime(2024,  5, 11, 1, 0, 0)
+t_end   = datetime(2024,  5, 11, 1, 0, 5)
+channel_start = 0
+channel_end = -1
+loader = loader("/pnfs/desy.de/m/project/iDAS/work/2024-05-11-desy", "SILIXA", 1)
+data = loader.load_array(t_start, t_end, channel_start, channel_end)
+
+print("Reduce data by binning (mean averaging)")
+bin_factors = (10, 10)
+data = utils.bin(data, bin_factors) # Reduce time sampling and spatial sampling by a factor of 10 by averaging each.
+sampling_hz = 1000.0 / bin_factors[0]
+channel_spacing = 1.0 * bin_factors[1]
+
+print("Create plot with pyplot")
+PP.title(f"{t_start.isoformat()}")
+PP.imshow(
+    data,
+    cmap = "seismic",
+    aspect = "auto",
+    interpolation = "nearest",
+    origin = "lower",
+    vmin = -2e-6,
+    vmax = +2e-6,
+    extent = (
+        channel_start, channel_start + (data.shape[1] * channel_spacing), 
+        data.shape[0] / sampling_hz, 0
+    )
+)
+PP.xlabel("Position [m]")
+PP.ylabel("Time [s]")
+PP.colorbar(label="Strain-rate [$\\frac{m}{m \\cdot s}$]")
+PP.savefig("waterfall.png")