DASPy-toolbox 1.1.1__tar.gz → 1.1.2__tar.gz

{DASPy-toolbox-1.1.1 → daspy_toolbox-1.1.2}/DASPy_toolbox.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
-Metadata-Version: 2.1
+Metadata-Version: 2.2
 Name: DASPy-toolbox
-Version: 1.1.1
+Version: 1.1.2
 Summary: DASPy is an open-source project dedicated to provide a python package for DAS (Distributed Acoustic Sensing) data processing, which comprises classic seismic data processing techniques and Specialized algorithms for DAS applications.
 Home-page: https://github.com/HMZ-03/DASPy
 Author: Minzhe Hu, Zefeng Li
@@ -22,14 +22,24 @@ Requires-Dist: h5py
 Requires-Dist: segyio
 Requires-Dist: nptdms
 Requires-Dist: tqdm
-
-<img src="./website/USTC.svg" height="170" />&emsp;<img src="./website/DAMS.png" height="150" />
-
-# DASPy
+Dynamic: author
+Dynamic: author-email
+Dynamic: classifier
+Dynamic: description
+Dynamic: home-page
+Dynamic: license
+Dynamic: maintainer
+Dynamic: maintainer-email
+Dynamic: requires-dist
+Dynamic: requires-python
+Dynamic: summary
+
+<img src="https://raw.github.com/hmz-03/daspy/main/website/logo.png" height="200" />
 
 [![Supported Python versions](https://img.shields.io/badge/python-3.9%20|%203.10%20|%203.11%20|%203.12-blue)](https://pypi.org/project/DASPy-toolbox/)
 [![License](https://img.shields.io/pypi/l/daspy-toolbox.svg)](https://opensource.org/license/mit)
 [![PyPI Version](https://img.shields.io/pypi/v/daspy-toolbox.svg)](https://pypi.org/project/DASPy-toolbox/)
+
 [![DOI](https://img.shields.io/badge/DOI-10.1785/0220240124-blue.svg)](https://doi.org/10.1785/0220240124)
 [![PyPI Downloads](https://img.shields.io/pypi/dm/daspy-toolbox.svg?label=pypi)](https://pypi.org/project/DASPy-toolbox/)
 [![Conda Downloads](https://img.shields.io/conda/dn/conda-forge/daspy-toolbox?label=conda)](https://anaconda.org/conda-forge/daspy-toolbox)

daspy_toolbox-1.1.2/DASPy_toolbox.egg-info/SOURCES.txt

@@ -0,0 +1,28 @@
+LICENSE
+README.md
+setup.py
+DASPy_toolbox.egg-info/PKG-INFO
+DASPy_toolbox.egg-info/SOURCES.txt
+DASPy_toolbox.egg-info/dependency_links.txt
+DASPy_toolbox.egg-info/entry_points.txt
+DASPy_toolbox.egg-info/requires.txt
+DASPy_toolbox.egg-info/top_level.txt
+daspy/__init__.py
+daspy/advanced_tools/__init__.py
+daspy/advanced_tools/channel.py
+daspy/advanced_tools/decomposition.py
+daspy/advanced_tools/denoising.py
+daspy/advanced_tools/fdct.py
+daspy/advanced_tools/strain2vel.py
+daspy/basic_tools/__init__.py
+daspy/basic_tools/filter.py
+daspy/basic_tools/freqattributes.py
+daspy/basic_tools/preprocessing.py
+daspy/basic_tools/visualization.py
+daspy/core/__init__.py
+daspy/core/collection.py
+daspy/core/dasdatetime.py
+daspy/core/example.pkl
+daspy/core/read.py
+daspy/core/section.py
+daspy/core/write.py

{DASPy-toolbox-1.1.1 → daspy_toolbox-1.1.2}/PKG-INFO

@@ -1,6 +1,6 @@
-Metadata-Version: 2.1
+Metadata-Version: 2.2
 Name: DASPy-toolbox
-Version: 1.1.1
+Version: 1.1.2
 Summary: DASPy is an open-source project dedicated to provide a python package for DAS (Distributed Acoustic Sensing) data processing, which comprises classic seismic data processing techniques and Specialized algorithms for DAS applications.
 Home-page: https://github.com/HMZ-03/DASPy
 Author: Minzhe Hu, Zefeng Li
@@ -22,14 +22,24 @@ Requires-Dist: h5py
 Requires-Dist: segyio
 Requires-Dist: nptdms
 Requires-Dist: tqdm
-
-<img src="./website/USTC.svg" height="170" />&emsp;<img src="./website/DAMS.png" height="150" />
-
-# DASPy
+Dynamic: author
+Dynamic: author-email
+Dynamic: classifier
+Dynamic: description
+Dynamic: home-page
+Dynamic: license
+Dynamic: maintainer
+Dynamic: maintainer-email
+Dynamic: requires-dist
+Dynamic: requires-python
+Dynamic: summary
+
+<img src="https://raw.github.com/hmz-03/daspy/main/website/logo.png" height="200" />
 
 [![Supported Python versions](https://img.shields.io/badge/python-3.9%20|%203.10%20|%203.11%20|%203.12-blue)](https://pypi.org/project/DASPy-toolbox/)
 [![License](https://img.shields.io/pypi/l/daspy-toolbox.svg)](https://opensource.org/license/mit)
 [![PyPI Version](https://img.shields.io/pypi/v/daspy-toolbox.svg)](https://pypi.org/project/DASPy-toolbox/)
+
 [![DOI](https://img.shields.io/badge/DOI-10.1785/0220240124-blue.svg)](https://doi.org/10.1785/0220240124)
 [![PyPI Downloads](https://img.shields.io/pypi/dm/daspy-toolbox.svg?label=pypi)](https://pypi.org/project/DASPy-toolbox/)
 [![Conda Downloads](https://img.shields.io/conda/dn/conda-forge/daspy-toolbox?label=conda)](https://anaconda.org/conda-forge/daspy-toolbox)

{DASPy-toolbox-1.1.1 → daspy_toolbox-1.1.2}/README.md

@@ -1,10 +1,9 @@
-<img src="./website/USTC.svg" height="170" />&emsp;<img src="./website/DAMS.png" height="150" />
-
-# DASPy
+<img src="https://raw.github.com/hmz-03/daspy/main/website/logo.png" height="200" />
 
 [![Supported Python versions](https://img.shields.io/badge/python-3.9%20|%203.10%20|%203.11%20|%203.12-blue)](https://pypi.org/project/DASPy-toolbox/)
 [![License](https://img.shields.io/pypi/l/daspy-toolbox.svg)](https://opensource.org/license/mit)
 [![PyPI Version](https://img.shields.io/pypi/v/daspy-toolbox.svg)](https://pypi.org/project/DASPy-toolbox/)
+
 [![DOI](https://img.shields.io/badge/DOI-10.1785/0220240124-blue.svg)](https://doi.org/10.1785/0220240124)
 [![PyPI Downloads](https://img.shields.io/pypi/dm/daspy-toolbox.svg?label=pypi)](https://pypi.org/project/DASPy-toolbox/)
 [![Conda Downloads](https://img.shields.io/conda/dn/conda-forge/daspy-toolbox?label=conda)](https://anaconda.org/conda-forge/daspy-toolbox)

{DASPy-toolbox-1.1.1 → daspy_toolbox-1.1.2}/daspy/advanced_tools/channel.py

@@ -84,7 +84,7 @@ def channel_checking(data, deg=10, thresh=5, continuity=True, adjacent=2,
         bad_chn = np.argwhere(deviation < -thresh * mad).ravel().tolist()
     elif mode == 'high':
         bad_chn = np.argwhere(deviation > thresh * mad).ravel().tolist()
-    elif mode == 'high':
+    elif mode == 'both':
         bad_chn = np.argwhere(deviation < -thresh * mad).ravel().tolist() + \
                 np.argwhere(deviation > thresh * mad).ravel().tolist()
     good_chn = list(set(range(nch)) - set(bad_chn))

{DASPy-toolbox-1.1.1 → daspy_toolbox-1.1.2}/daspy/basic_tools/visualization.py

@@ -1,13 +1,13 @@
 # Purpose: Plot data
 # Author: Minzhe Hu
-# Date: 2024.11.8
+# Date: 2025.1.6
 # Email: hmz2018@mail.ustc.edu.cn
 import numpy as np
 import matplotlib.pyplot as plt
 from collections.abc import Sequence
 
 
-def plot(data: np.ndarray, dx=None, fs=None, ax=None, obj='waveform', dpi=150,
+def plot(data: np.ndarray, dx=None, fs=None, ax=None, obj='waveform', dpi=300,
          title=None, transpose=False, t0=0, x0=0, pick=None, f=None, k=None,
          t=None, c=None, cmap=None, vmin=None, vmax=None, dB=False,
          xmode='distance', tmode='time', xlim=None, ylim=None, xlog=False,
@@ -29,7 +29,9 @@ def plot(data: np.ndarray, dx=None, fs=None, ax=None, obj='waveform', dpi=150,
     :param title: str. The title of this axes.
     :param transpose: bool. Transpose the figure or not.
     :param t0, x0: The beginning of time and space.
-    :param pick: Sequence of picked phases. Required if obj=='phasepick'.
+    :param pick: dictionary of sequence of picked phases. Key should be 'P' for
+        P phase, 'S' for S phase and 'N' for unknown phase type. Required if
+        obj=='phasepick'.
     :param f: Sequence of frequency. Required if obj is one of 'spectrum',
         'spectrogram', 'fk' or 'dispersion'.
     :param k: Wavenumber sequence. Required if obj=='fk'.
@@ -82,14 +84,18 @@ def plot(data: np.ndarray, dx=None, fs=None, ax=None, obj='waveform', dpi=150,
             extent = [x0 * 1e-3, (x0 + nch * dx) * 1e-3, t0 + nt / fs, t0]
 
         if obj == 'phasepick' and len(pick):
-            pck = np.array(pick).astype(float)
-            if xmode.lower() == 'distance':
-                pck[:, 0] = (x0 + pck[:, 0] * dx) * 1e-3
-            elif xmode.lower() == 'channel':
-                pck[:, 0] = x0 + pck[:, 0]
-            if tmode.lower() == 'sampling':
-                pck[:, 1] = pck[:, 1] / fs
-            ax.scatter(pck[:,0], t0 + pck[:,1], marker=',', s=0.1, c='black')
+            pick_color = {'P': 'r', 'S': 'b', 'N': 'k'}
+            for phase, pck in pick.items():
+                if len(pck):
+                    pck = np.array(pck).astype(float)
+                    if xmode.lower() == 'distance':
+                        pck[:, 0] = (x0 + pck[:, 0] * dx) * 1e-3
+                    elif xmode.lower() == 'channel':
+                        pck[:, 0] = x0 + pck[:, 0]
+                    if tmode.lower() == 'sampling':
+                        pck[:, 1] = pck[:, 1] / fs
+                    ax.scatter(pck[:,0], t0 + pck[:,1], marker=',', s=0.1,
+                            c=pick_color[phase])
 
     elif obj in ['spectrum', 'spectrogram', 'fk', 'dispersion']:
         if isinstance(data[0,0], (complex, np.complex64)):

{DASPy-toolbox-1.1.1 → daspy_toolbox-1.1.2}/daspy/core/collection.py

@@ -1,6 +1,6 @@
 # Purpose: Module for handling Collection objects.
 # Author: Minzhe Hu
-# Date: 2025.1.1
+# Date: 2025.1.23
 # Email: hmz2018@mail.ustc.edu.cn
 import os
 import warnings
@@ -63,7 +63,7 @@ class Collection(object):
                     sec.gauge_length = None
                 ftime.append(sec.start_time)
                 metadata_list.append((sec.nch, sec.nt, sec.dx, sec.fs,
-                                      sec.gauge_length))
+                                      sec.gauge_length, sec.duration))
 
             if len(set(metadata_list)) > 1:
                 warnings.warn('More than one kind of setting detected.')
@@ -71,6 +71,8 @@ class Collection(object):
             for i, key in enumerate(['nch', 'nt', 'dx', 'fs', 'gauge_length']):
                 if not hasattr(self, key):
                     setattr(self, key, metadata[i])
+                if flength is None:
+                    flength = metadata[-1]
             self.ftime = ftime
         elif meta_from_file:
             i = int(len(self.flist) > 1)
@@ -211,6 +213,9 @@ class Collection(object):
                 flist.append(self.flist[i])
                 ftime.append(self.ftime[i])
 
+        if len(flist) == 0:
+            warnings.warn('Out of collection time range.')
+            return None
         if readsec:
             sec = read(flist[0], **kwargs)
             for f in flist[1:]:
@@ -218,10 +223,9 @@ class Collection(object):
             sec.trimming(tmin=stime, tmax=etime)
             return sec
         else:
-            coll = self.copy()
-            coll.flist = flist
-            coll.ftime = ftime
-            return coll
+            self.flist = flist
+            self.ftime = ftime
+            return self
 
     def _optimize_for_continuity(self, operations):
         method_list = []
@@ -231,8 +235,8 @@ class Collection(object):
         for opera in operations:
             method, kwargs = opera
             if method == 'downsampling':
-                if hasattr(kwargs, 'lowpass_filter') and not\
-                        kwargs['lowpass_filter']:
+                if (hasattr(kwargs, 'lowpass_filter') and not\
+                    kwargs['lowpass_filter']) or not hasattr(kwargs, 'tint'):
                     method_list.append('downsampling')
                     kwargs_list.append(kwargs)
                 else:
@@ -271,6 +275,16 @@ class Collection(object):
             merge = len(self)
         for i in tqdm(range(0, len(self))):
             f = self[i]
+            if os.path.getsize(f) == 0:
+                for j, method in enumerate(method_list):
+                    if method == 'time_integration':
+                        kwargs_list[j]['c'] = 0
+                    elif method == 'time_differential':
+                        kwargs_list[j]['prepend'] = 0
+                    elif method in ['bandpass', 'bandstop', 'lowpass',
+                                    'highpass', 'lowpass_cheby_2']:
+                        kwargs_list[j]['zi'] = 0
+                continue
             sec = read(f, ftype=self.ftype, **read_kwargs)
             for j, method in enumerate(method_list):
                 if method in ['taper', 'cosine_taper']:

{DASPy-toolbox-1.1.1 → daspy_toolbox-1.1.2}/daspy/core/dasdatetime.py

@@ -1,6 +1,6 @@
 # Purpose: Module for handling DASDateTime objects.
 # Author: Minzhe Hu
-# Date: 2024.11.18
+# Date: 2025.2.11
 # Email: hmz2018@mail.ustc.edu.cn
 import time
 from typing import Iterable
@@ -77,4 +77,26 @@ class DASDateTime(datetime):
 
     def to_obspy_UTCDateTime(self):
         from obspy import UTCDateTime
-        return UTCDateTime(UTCDateTime(self.to_datetime()))
+        return UTCDateTime(UTCDateTime(self.to_datetime()))
+
+    @classmethod
+    def strptime(cls, date_string, format):
+        """
+        string, format -> new datetime parsed from a string
+        (like time.strptime()).
+        """
+        from _strptime import _strptime
+        tt, fraction, gmtoff_fraction = _strptime(date_string, format)
+        tzname, gmtoff = tt[-2:]
+        args = tt[:6] + (fraction,)
+        if gmtoff is not None:
+            tzdelta = timedelta(seconds=gmtoff, microseconds=gmtoff_fraction)
+            if tzname:
+                tz = timezone(tzdelta, tzname)
+            else:
+                tz = timezone(tzdelta)
+            args += (tz,)
+        elif tt[-3] == 0:
+            args += (utc,)
+
+        return cls(*args)

{DASPy-toolbox-1.1.1 → daspy_toolbox-1.1.2}/daspy/core/section.py

@@ -1,6 +1,6 @@
 # Purpose: Module for handling Section objects.
 # Author: Minzhe Hu
-# Date: 2024.12.12
+# Date: 2025.2.25
 # Email: hmz2018@mail.ustc.edu.cn
 import warnings
 import os
@@ -517,7 +517,7 @@ class Section(object):
             it will be set to parameter obj.
         """
         if 'data' not in kwargs.keys():
-            if obj == 'waveform':
+            if obj == 'waveform' or obj == 'phasepick':
                 data = deepcopy(self.data)
             elif obj == 'spectrum':
                 data, f = self.spectrum(**kwargs_pro)
@@ -537,8 +537,13 @@ class Section(object):
 
         if 'ax' not in kwargs.keys() and 'title' not in kwargs.keys():
             kwargs['title'] = obj
-            if hasattr(self, 'data_type'):
-                kwargs['title'] += f' ({self.data_type})'
+
+        if obj == 'phasepick':
+            for phase, pck in kwargs['pick'].items():
+                if len(pck) != 0:
+                    pck = np.array(pck)
+                    pck[:, 0] -= self.start_channel
+                    kwargs['pick'][phase] = pck
 
         if xmode == 'channel':
             kwargs.setdefault('x0', self.start_channel)
@@ -611,8 +616,8 @@ class Section(object):
         """
         Normalize for each individual channel using Z-score method.
 
-        :param method: str. Method for normalization, should be one of 'max' or
-            'z-score'.
+        :param method: str. Method for normalization, should be one of 'max',
+            'z-score', or 'one-bit'.
         """
         self.data = normalization(self.data, method=method)
         if hasattr(self, 'data_type'):
@@ -1053,11 +1058,13 @@ class Section(object):
         :return: self or Good channels and bad channels.
         """
         good_chn, bad_chn = channel_checking(self.data, **kwargs)
+        good_chn = good_chn + self.start_channel
+        bad_chn = bad_chn + self.start_channel
         if use:
             self.channel_data(good_chn, replace=True)
             return self
         else:
-            return good_chn + self.start_channel, bad_chn + self.start_channel
+            return good_chn, bad_chn
 
     def turning_points(self, data_type='default', **kwargs):
         """

{DASPy-toolbox-1.1.1 → daspy_toolbox-1.1.2}/setup.py

@@ -2,7 +2,7 @@ from setuptools import setup, find_packages
 
 
 setup(
-    name='DASPy-toolbox', version='1.1.1',
+    name='DASPy-toolbox', version='1.1.2',
     description=(
         'DASPy is an open-source project dedicated to provide a python package '
         'for DAS (Distributed Acoustic Sensing) data processing, which '

DASPy-toolbox-1.1.1/DASPy_toolbox.egg-info/SOURCES.txt

@@ -1,53 +0,0 @@
-LICENSE
-README.md
-setup.py
-DASPy_toolbox.egg-info/PKG-INFO
-DASPy_toolbox.egg-info/SOURCES.txt
-DASPy_toolbox.egg-info/dependency_links.txt
-DASPy_toolbox.egg-info/entry_points.txt
-DASPy_toolbox.egg-info/requires.txt
-DASPy_toolbox.egg-info/top_level.txt
-daspy/__init__.py
-daspy/advanced_tools/__init__.py
-daspy/advanced_tools/channel.py
-daspy/advanced_tools/decomposition.py
-daspy/advanced_tools/denoising.py
-daspy/advanced_tools/fdct.py
-daspy/advanced_tools/strain2vel.py
-daspy/basic_tools/__init__.py
-daspy/basic_tools/filter.py
-daspy/basic_tools/freqattributes.py
-daspy/basic_tools/preprocessing.py
-daspy/basic_tools/visualization.py
-daspy/core/__init__.py
-daspy/core/collection.py
-daspy/core/dasdatetime.py
-daspy/core/example.pkl
-daspy/core/make_example.py
-daspy/core/read.py
-daspy/core/section.py
-daspy/core/write.py
-daspy/seismic_detection/__init__.py
-daspy/seismic_detection/calc_travel_time.py
-daspy/seismic_detection/core.py
-daspy/seismic_detection/detection.py
-daspy/seismic_detection/location.py
-daspy/seismic_detection/magnitude.py
-daspy/seismic_detection/phase_picking.py
-daspy/seismic_detection/gamma/__init__.py
-daspy/seismic_detection/gamma/_base.py
-daspy/seismic_detection/gamma/_bayesian_mixture.py
-daspy/seismic_detection/gamma/_gaussian_mixture.py
-daspy/seismic_detection/gamma/app.py
-daspy/seismic_detection/gamma/seismic_ops.py
-daspy/seismic_detection/gamma/utils.py
-daspy/structure_imaging/__init__.py
-daspy/structure_imaging/ambient_noise.py
-daspy/structure_imaging/core.py
-daspy/structure_imaging/dispersion.py
-daspy/structure_imaging/fault_zone.py
-daspy/structure_imaging/inversion.py
-daspy/traffic_monitoring/JamDetection.py
-daspy/traffic_monitoring/SpeedMeasurement.py
-daspy/traffic_monitoring/VehicleDetection.py
-daspy/traffic_monitoring/__init__.py

DASPy-toolbox-1.1.1/daspy/core/make_example.py

@@ -1,32 +0,0 @@
-import h5py
-import numpy as np
-from daspy import read, Section, DASDateTime
-from daspy.core.dasdatetime import utc
-
-
-origin_time = DASDateTime(2016, 3, 21, 7, 37, 10, 535000, tzinfo=utc)
-
-# read DAS data
-dx, fs = 1, 1000
-sec = Section(np.zeros((8721, 0)), dx, fs, data_type='Strain rate')
-filename = ['PoroTomo_iDAS16043_160321073651.h5',
-            'PoroTomo_iDAS16043_160321073721.h5',
-            'PoroTomo_iDAS16043_160321073751.h5',
-            'PoroTomo_iDAS16043_160321073821.h5']
-
-first = True
-for fn in filename:
-    with h5py.File(fn,'r') as fp:
-        if first:
-            sec.start_time = DASDateTime.fromtimestamp(fp['t'][0]).astimezone(utc)
-            sec.start_channel = fp['channel'][0]
-            first = False
-        sec += fp['das'][()].T
-
-sec.origin_time = origin_time
-sec.trimming(mode=0, xmin=2500, xmax=3000)
-sec.trimming(tmin=origin_time, tmax=origin_time+90)
-sec.downsampling(tint=10)
-sec.trimming(tmin=origin_time+20, tmax=origin_time+70)
-
-sec.save('example.py')

DASPy-toolbox-1.1.1/daspy/seismic_detection/__init__.py

@@ -1 +0,0 @@
-from daspy.seismic_detection.core import DetectSection

DASPy-toolbox-1.1.1/daspy/seismic_detection/calc_travel_time.py

@@ -1,23 +0,0 @@
-import sys
-import os
-import multiprocessing
-import numpy as np
-from obspy.taup.tau_model import TauModel
-from obspy.taup import TauPyModel
-from daspy.seismic_detection.location import calc_travel_time
-
-
-cores = multiprocessing.cpu_count() - 1
-pool = multiprocessing.Pool(processes=cores)
-
-
-filepath, elo, ela, edep, savefile = sys.argv[1:6]
-
-elo, ela, edep = map(float, (elo, ela, edep))
-model = TauPyModel()
-model.model = TauModel.from_file(os.path.join(filepath, 'model.npz'))
-geometry = np.load(os.path.join(filepath, 'geometry.npy'))
-tasks = [(model, elo, ela, edep, *geo) for geo in geometry]
-delay = pool.starmap(calc_travel_time, tasks)
-delay = np.array(delay)
-np.save(os.path.join(filepath, savefile), delay)

DASPy-toolbox-1.1.1/daspy/seismic_detection/core.py

@@ -1,186 +0,0 @@
-import numpy as np
-from collections.abc import Iterable
-from tqdm import tqdm
-import daspy.basic_tools.visualization as DBV
-from daspy.core import Section, DASDateTime
-from daspy.seismic_detection.phase_picking import *
-
-
-def plot(data, obj='waveform', title=None, cmap=None, vmax=None, vmin=None,
-         **kwargs):
-    if obj == 'pickmap':
-        cmap = 'hot' if cmap is None else cmap
-        vmax = np.percentile(abs(data), 99) if vmax is None else vmax
-        vmin = np.percentile(abs(data), 80) if vmin is None else vmin
-        title = obj
-        obj = 'waveform'
-
-    return DBV.plot(data, title=title, cmap=cmap, vmin=vmin, vmax=vmax,
-                    **kwargs)
-
-class DetectSection(Section):
-
-    def __str__(self):
-        describe = ''
-        n = max(map(len, self.__dict__.keys()))
-        for key, value in self.__dict__.items():
-            if key in ['data', 'geometry', 'pickmap']:
-                describe = '{}: shape{}\n'.format(key.rjust(n), value.shape) \
-                    + describe
-            elif key in ['dx', 'start_distance', 'gauge_length']:
-                describe += '{}: {} m\n'.format(key.rjust(n), value)
-            elif key == 'fs':
-                describe += '{}: {} Hz\n'.format(key.rjust(n), value)
-            elif key == 'start_time':
-                if isinstance(value, DASDateTime):
-                    describe += '{}: {}\n'.format(key.rjust(n), value)
-                else:
-                    describe += '{}: {} s\n'.format(key.rjust(n), value)
-            elif key == 'pick':
-                describe += '{}: {}\n'.format(key.rjust(n), len(value))
-            else:
-                describe += '{}: {}\n'.format(key.rjust(n), value)
-        return describe
-    
-    __repr__ = __str__
-
-    def plot(self, obj='waveform', kwargs_pro={}, **kwargs):
-        if obj == 'pickmap':
-            if 'data' not in kwargs.keys():
-                if not hasattr(self, 'pickmap') or len(kwargs_pro):
-                    self.calc_pickmap(**kwargs_pro)
-                kwargs['data']= self.pickmap
-            kwargs.setdefault('cmap', 'hot')
-            kwargs.setdefault('vmax', np.percentile(abs(kwargs['data']), 99.9))
-            kwargs.setdefault('vmin', np.percentile(abs(kwargs['data']), 95))
-            kwargs.setdefault('title', obj)
-            obj = 'waveform'
-        elif obj == 'phasepick':
-            if 'data' not in kwargs.keys():
-                kwargs['data']= self.data
-            if 'pick' not in kwargs.keys():
-                if not hasattr(self, 'pick') or len(kwargs_pro):
-                    self.map_picking(**kwargs_pro)
-                kwargs['pick'] = self.pick
-                if len(kwargs['pick']):
-                    kwargs['pick'][:, 0] -= self.start_channel
-
-        super().plot(obj=obj, kwargs_pro=kwargs_pro, **kwargs)
-
-    def calc_pickmap(self, method='sta_lta', **kwargs):
-        if isinstance(method, str):
-            method = [method]
-        self.pickmap = np.zeros_like(self.data)
-        for m in method:
-            if m == 'sta_lta':
-                self.pickmap += sta_lta_map(self.data, self.fs, **kwargs)
-            elif m == 'kurto':
-                self.pickmap += kurto_map(self.data, self.fs, **kwargs)
-            elif m == 'skew':
-                self.pickmap += skew_map(self.data, self.fs, **kwargs)
-
-        return self.pickmap
-
-    def map_picking(self, thres1=5, thres2=5, choose_max=False, min_dt=0,
-                    **kwargs):
-        if not hasattr(self, 'pickmap') or len(kwargs):
-            self.calc_pickmap(**kwargs)
-
-        pick, amp = map_picking(self.pickmap, thres1=thres1, thres2=thres2,
-                                choose_max=choose_max,
-                                min_dsp=min_dt*self.fs).astype(float)
-        if len(pick):
-            pick[:, 0] += self.start_channel
-            pick[:, 1] = pick[:, 1] / self.fs
-
-        self.pick = pick
-        self.pick_amp = amp
-        return pick
-
-    def symmetry_detection(self, win=5000):
-            sec = self.copy()
-            sec.normalization()
-            win_ch = round(win / self.dx)
-            cc = np.zeros(self.nch)
-            for sch in range(self.nch):
-                win_ch_use = min(win_ch, sch, sec.nch-sch-1)
-                if win_ch_use > 0:
-                    panel1 = sec.data[sch-win_ch_use:sch]
-                    panel1 = panel1[::-1]
-                    panel2 = sec.data[sch+1:sch+win_ch_use+1]
-                    cc[sch] = np.sum(panel1 * panel2) / win_ch
-            self.signal_channel = self.start_channel + np.argmax(cc)
-            return cc / sec.nt
-
-
-    def trimming_pick(self, pick, t_win=0.5):
-        sp = round(t_win * self.fs)
-        if isinstance(pick, Iterable):
-            if len(pick) != self.nch:
-                raise ValueError('Length of pick should equals to channel '
-                                 'number.')
-            if isinstance(pick[0], DASDateTime):
-                pick = [p - self.start_time for p in pick]
-            pick = np.round(np.array(pick) * self.fs).astype(int)
-            data = np.zeros((self.nch, 2 * sp))
-            for i, p in enumerate(pick):
-                data[i] = self.data[i, p-sp:p+sp]
-            self.data = data
-            self.start_time = -t_win
-        else:
-            self.trimming(tmin=pick-t_win, tmax=pick+t_win)
-        return self
-
-
-    def relative_arrival(self, pick, t_win=0.5, dch=30, dt_max=0.1, thresh=0,
-                         positive=True, weight_pow=1):
-        dsp_max = round(dt_max * self.fs)
-        sec_win1 = self.copy()
-        sec_win2 = self.copy()
-        sec_win1.trimming_pick(pick, t_win=t_win)
-        sec_win2.trimming_pick(pick, t_win=t_win+dt_max)
-        G = []
-        d = []
-        for i in tqdm(range(self.nch)):
-            for j in range(1, dch+1):
-                if i+j >= self.nch:
-                    continue
-                cc = xcorr(sec_win2.data[i], sec_win1.data[i+j])
-                if not positive:
-                    cc = np.abs(cc)
-                w = max(cc)
-                if w < thresh:
-                    continue
-                w = w ** weight_pow
-                d.append((np.argmax(cc) - dsp_max) * w)
-                Gl = np.zeros(self.nch-1)
-                Gl[i:i+j] = w
-                G.append(Gl)
-
-        G = np.array(G)
-        idx = np.where(np.sum(G, axis=0) > 0)[0]
-        print(len(idx))
-        G = np.mat(G[:, idx])
-        print(G.shape)
-        d = np.mat(d).T
-        try:
-            m = (G.T * G).I * G.T * d
-        except:
-            m = G.I * d
-            print('xxx')
-        time_shift_difference = np.zeros(self.nch - 1)
-        time_shift_difference[idx] = np.array(m)[:,0]
-        time_shift = np.cumsum(time_shift_difference) / self.fs
-        time_shift = np.insert(time_shift, 0, 0)
-        if isinstance(pick, Iterable):
-            if isinstance(pick[0], DASDateTime):
-                pick = [p - self.start_time for p in pick]
-
-        return - time_shift + np.array(pick)
-
-
-    @classmethod
-    def from_section(clc, raw_sec):
-        raw_dict = raw_sec.__dict__
-        data = raw_dict.pop('data')
-        return clc(data, **raw_dict)