DASPy-toolbox 1.2.0__tar.gz → 1.2.1__tar.gz

{daspy_toolbox-1.2.0 → daspy_toolbox-1.2.1}/DASPy_toolbox.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: DASPy-toolbox
-Version: 1.2.0
+Version: 1.2.1
 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

{daspy_toolbox-1.2.0 → daspy_toolbox-1.2.1}/DASPy_toolbox.egg-info/SOURCES.txt

@@ -23,6 +23,7 @@ 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/util.py

{daspy_toolbox-1.2.0 → daspy_toolbox-1.2.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: DASPy-toolbox
-Version: 1.2.0
+Version: 1.2.1
 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

{daspy_toolbox-1.2.0 → daspy_toolbox-1.2.1}/daspy/advanced_tools/denoising.py

@@ -1,6 +1,6 @@
 # Purpose: Remove noise from data
 # Author: Minzhe Hu, Zefeng Li
-# Date: 2024.5.13
+# Date: 2025.9.18
 # Email: hmz2018@mail.ustc.edu.cn
 import numpy as np
 from copy import deepcopy
@@ -57,7 +57,7 @@ def common_mode_noise_removal(data, method='median'):
         common = np.mean(data, 0)
 
     xx = np.sum(common ** 2)
-    data_dn = np.zeros((nch, nt))
+    data_dn = np.zeros((nch, nt), dtype=data.dtype)
     for i in range(nch):
         xc = np.sum(common * data[i])
         data_dn[i] = data[i] - xc / xx * common

{daspy_toolbox-1.2.0 → daspy_toolbox-1.2.1}/daspy/basic_tools/preprocessing.py

@@ -1,9 +1,10 @@
 # Purpose: Some preprocess methods
 # Author: Minzhe Hu
-# Date: 2025.6.25
+# Date: 2025.10.30
 # Email: hmz2018@mail.ustc.edu.cn
 import warnings
 import numpy as np
+from numpy.fft import rfft, irfft, rfftfreq
 from scipy.signal import detrend
 from scipy.signal.windows import tukey
 from daspy.basic_tools.filter import lowpass_cheby_2
@@ -320,7 +321,7 @@ def padding(data, dn, reverse=False):
         return data_pd
 
 
-def time_integration(data, fs, c=0):
+def time_integration(data, fs, domain='time', c=0):
     """
     Integrate DAS data in time.
 
@@ -329,10 +330,19 @@ def time_integration(data, fs, c=0):
     :param c: float. A constant added to the result.
     :return: Integrated data.
     """
-    return np.cumsum(data, axis=1) / fs + c
-
-
-def time_differential(data, fs, prepend=0):
+    if domain == 'time':
+        return np.cumsum(data, axis=1) / fs + c
+    elif domain in ['frequency', 'freq']:
+        nsp = data.shape[1]
+        freqs = rfftfreq(nsp, d=1/fs)
+        spectrum = rfft(data, axis=1)
+        H = np.zeros_like(freqs, dtype=complex)
+        nonzero = freqs != 0
+        H[nonzero] = 1 / (1j * 2 * np.pi * freqs[nonzero])
+        return np.real(irfft(spectrum * H))
+
+
+def time_differential(data, fs, domain='time', prepend=0):
     """
     Differentiate DAS data in time.
 
@@ -342,9 +352,16 @@ def time_differential(data, fs, prepend=0):
         performing the difference. 
     :return: Differentiated data.
     """
-    if prepend == 'mean':
-        prepend = np.mean(data, axis=1).reshape((-1, 1))
-    return np.diff(data, axis=1, prepend=prepend) * fs
+    if domain == 'time':
+        if prepend == 'mean':
+            prepend = np.mean(data, axis=1).reshape((-1, 1))
+        return np.diff(data, axis=1, prepend=prepend) * fs
+    elif domain in ['frequency', 'freq']:
+        nsp = data.shape[1]
+        freqs = rfftfreq(nsp, d=1./fs)
+        spectrum = rfft(data, axis=1)
+        H = 1j * 2 * np.pi * freqs
+        return np.real(irfft(spectrum * H))
 
 
 def distance_integration(data, dx, c=0):

{daspy_toolbox-1.2.0 → daspy_toolbox-1.2.1}/daspy/core/collection.py

@@ -1,6 +1,6 @@
 # Purpose: Module for handling Collection objects.
 # Author: Minzhe Hu
-# Date: 2025.9.15
+# Date: 2025.9.18
 # Email: hmz2018@mail.ustc.edu.cn
 import os
 import warnings
@@ -277,10 +277,11 @@ class Collection(object):
             if len(flist) == 1:
                 sec = read(flist[0], tmin=tmin, tmax=tmax, **kwargs)
             else:
-                sec = read(flist[0], tmin=tmin, **kwargs)
+                sec = read(flist[0], tmin=tmin-tolerance, **kwargs)
                 for f in flist[1:-1]:
                     sec += read(f, **kwargs)
-                sec += read(flist[-1], tmax=tmax, **kwargs)
+                sec += read(flist[-1], tmax=tmax+tolerance, **kwargs)
+                sec.trimming(tmin=tmin, tmax=tmax)
             return sec
         else:
             self.flist = flist
@@ -326,8 +327,8 @@ class Collection(object):
                 kwargs_list[j]['zi'] = 0
 
     def process(self, operations, savepath='./processed', merge=1,
-                suffix='_pro', ftype=None, dtype=None, save_operations=False,
-                tolerance=0.5, **read_kwargs):
+                suffix='_pro', ftype=None, dtype=None, file_format='auto',
+                save_operations=False, tolerance=0.5, **read_kwargs):
         """
         :param operations: list or None. Each element of operations list
             should be [str of method name, dict of kwargs]. None for read
@@ -339,6 +340,8 @@ class Collection(object):
         :param ftype: None or str. File format for saving. None for automatic
             detection, or 'pkl', 'pickle', 'tdms', 'h5', 'hdf5', 'segy', 'sgy',
             'npy'.
+        :param file_format: Format in which the file is saved.
+        :type file_format: str or function
         :param dtype: str. The data type of the saved data.
         :parma save_operations: bool. If True, save the operations to
             method_list.pkl and kwargs_list.pkl in savepath.
@@ -369,7 +372,8 @@ class Collection(object):
                     warnings.warn(f'{f} is an empty file. Continuous data is '
                                   'interrupted here.')
                     if m > 0:
-                        sec_merge.save(filepath, dtype=dtype)
+                        sec_merge.save(filepath, file_format=file_format,
+                                       dtype=dtype)
                         m = 0
                     self._kwargs_initialization(method_list, kwargs_list)
                     continue
@@ -377,7 +381,8 @@ class Collection(object):
                     sec = read(f, ftype=self.ftype, **read_kwargs)
                     if sec.data.size == 0:
                         if m > 0:
-                            sec_merge.save(filepath, dtype=dtype)
+                            sec_merge.save(filepath, ftype=ftype,
+                                           file_format=file_format, dtype=dtype)
                             m = 0
                         self._kwargs_initialization(method_list, kwargs_list)
                         continue
@@ -385,23 +390,24 @@ class Collection(object):
                     warnings.warn(f'Error reading {f}: {e}. Continuous data is '
                                   'interrupted here.')
                     if m > 0:
-                        sec_merge.save(filepath, dtype=dtype)
+                        sec_merge.save(filepath, ftype=ftype,
+                                       file_format=file_format, dtype=dtype)
                         m = 0
                     self._kwargs_initialization(method_list, kwargs_list)
                     continue
                 for j, method in enumerate(method_list):
                     if method in ['taper', 'cosine_taper']:
                         if not ((i==0 and kwargs_list[j]['side'] != 'right') or
-                                (i == len(self) - 1 and kwargs_list[j]['side'] !=
-                                'left')):
+                                (i == len(self) - 1 and kwargs_list[j]['side']
+                                 != 'left')):
                             continue
                     out = getattr(sec, method)(**kwargs_list[j])
                     if method == 'time_integration':
                         kwargs_list[j]['c'] = sec.data[:, -1].copy()
                     elif method == 'time_differential':
                         kwargs_list[j]['prepend'] = sec.data[:, -1].copy()
-                    elif method in ['bandpass', 'bandstop', 'lowpass', 'highpass',
-                                    'lowpass_cheby_2']:
+                    elif method in ['bandpass', 'bandstop', 'lowpass',
+                                    'highpass', 'lowpass_cheby_2']:
                         kwargs_list[j]['zi'] = out
                 
                 if m == 0:
@@ -415,7 +421,8 @@ class Collection(object):
                     warnings.warn(f'The start time of {f} does not correspond '
                                   'to the end time of the previous file. '
                                   'Continuous data is interrupted here.')
-                    sec_merge.save(filepath, dtype=dtype)
+                    sec_merge.save(filepath, ftype=ftype,
+                                   file_format=file_format, dtype=dtype)
                     sec_merge = sec
                     f0, f1 = os.path.splitext(os.path.basename(f))
                     f1 = f1 if ftype is None else ftype
@@ -423,10 +430,12 @@ class Collection(object):
                     m = 0
                 m += 1
                 if m == merge:
-                    sec_merge.save(filepath, dtype=dtype)
+                    sec_merge.save(filepath, ftype=ftype,
+                                   file_format=file_format, dtype=dtype)
                     m = 0
             if m > 0:
-                sec_merge.save(filepath, dtype=dtype)
+                sec_merge.save(filepath, ftype=ftype, file_format=file_format,
+                               dtype=dtype)
         except KeyboardInterrupt as e:
             with open(method_file, 'wb') as f:
                 pickle.dump(method_list, f)
@@ -452,25 +461,16 @@ class Collection(object):
 def _create_cascade_method(method_name):
     def cascade_method(self, savepath='./processed', merge=1,
                        suffix=f'_{method_name}', ftype=None, dtype=None,
-                       save_operations=False, **kwargs):
+                       file_format='auto', save_operations=False, tolerance=0.5,
+                       **kwargs):
         """
         Automatically generated method for {method_name}.
         Applies the {method_name} operation to the data and saves the result.
-
-        :param savepath: str. Path to save processed files.
-        :param merge: int or str. int for merge several processed files into 1.
-            'all' for merge all files.
-        :param suffix: str. Suffix for processed files.
-        :param ftype: None or str. None for automatic detection, or 'pkl',
-            'pickle', 'tdms', 'h5', 'hdf5', 'segy', 'sgy', 'npy'.
-        :param dtype: str. The data type of the saved data.
-        :parma save_operations: bool. If True, save the operations to
-            method_list.pkl and kwargs_list.pkl in savepath.
-        :param kwargs: dict. Parameters for the {method_name} operation.
         """
         operations = [[method_name, kwargs]]
         self.process(operations, savepath=savepath, merge=merge, suffix=suffix,
-                     ftype=ftype, dtype=dtype, save_operations=save_operations)
+                     ftype=ftype, dtype=dtype, file_format=file_format,
+                     save_operations=save_operations, tolerance=tolerance)
     return cascade_method
 
 

{daspy_toolbox-1.2.0 → daspy_toolbox-1.2.1}/daspy/core/dasdatetime.py

@@ -84,7 +84,7 @@ class DASDateTime(datetime):
             match = re.match(r'(.*)(UTC|GMT)([+-]?\d{1,2})(.*)', date_string,
                              re.IGNORECASE)
             if match:
-                dt1, tz_prefix, offset, dt2 = match.groups()
+                dt1, _, offset, dt2 = match.groups()
                 offset_hours = int(offset)
                 tz = timezone(timedelta(hours=offset_hours))
                 return cls.strptime(dt1 + dt2, format.replace('%Z', '')).\

daspy_toolbox-1.2.1/daspy/core/make_example.py

@@ -0,0 +1,32 @@
+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.2.0 → daspy_toolbox-1.2.1}/daspy/core/read.py

@@ -1,6 +1,6 @@
 # Purpose: Module for reading DAS data.
 # Author: Minzhe Hu
-# Date: 2025.9.15
+# Date: 2025.10.30
 # Email: hmz2018@mail.ustc.edu.cn
 # Partially modified from
 # https://github.com/RobbinLuo/das-toolkit/blob/main/DasTools/DasPrep.py
@@ -188,7 +188,8 @@ def _read_h5(fname, headonly=False, file_format='auto', chmin=None, chmax=None,
                 start_time = DASDateTime.fromisoformat(h5_file[data_key].
                                                        attrs['start_time'])
             else:
-                start_time = DASDateTime.fromtimestamp(float(header[100])).utc()
+                start_time = DASDateTime.fromtimestamp(float(header[100])).\
+                    utc()
                 transpose = True
             metadata = {'dx': header[1],
                         'fs': header[6] / header[15] / header[98],
@@ -216,13 +217,9 @@ def _read_h5(fname, headonly=False, file_format='auto', chmin=None, chmax=None,
                 fs = float(attrs['FreqRes'])
             except KeyError:
                 try:
-                    fs = (attrs['PulseRateFreq'][0] /
-                            attrs['SamplingRes'][0]) / 1000
+                    fs = float(1000 / attrs['Spacing'][1])
                 except KeyError:
-                    try:
-                        fs = 1000 / attrs['Spacing'][1]
-                    except KeyError:
-                        fs = attrs['SamplingRate'][0]
+                    fs = attrs['PulseRateFreq'][0]
             time = h5_file[f'{group}/Source1/time']
             if len(time.shape) == 2: # Febus A1-R
                 start_time = DASDateTime.fromtimestamp(time[0, 0]).utc()
@@ -260,11 +257,9 @@ def _read_h5(fname, headonly=False, file_format='auto', chmin=None, chmax=None,
                              'Smart Earth ZD-DAS', 'Unknown']:
             dataset = h5_file['Acquisition/Raw[0]/RawData/']
             attrs = h5_file['Acquisition'].attrs
-            try:
+            if 'NumberOfLoci' in attrs.keys():
                 if dataset.shape[0] != attrs['NumberOfLoci']:
                     transpose = True
-            except KeyError:
-                pass
 
             try:
                 fs = h5_file['Acquisition/Raw[0]'].attrs['OutputDataRate']
@@ -282,7 +277,7 @@ def _read_h5(fname, headonly=False, file_format='auto', chmin=None, chmax=None,
                         stime = h5_file['Acquisition/Raw[0]/RawDataTime/'][0]
                     except KeyError:
                         stime = 0
-                        
+
             if isinstance(stime, bytes):
                 stime = stime.decode('ascii')
 
@@ -415,11 +410,14 @@ def _read_h5(fname, headonly=False, file_format='auto', chmin=None, chmax=None,
         elif transpose:
             if len(dataset.shape) == 3:
                 fs = int(metadata['fs'])
-                j0, k0 = divmod(sj.start, fs)
-                j1, k1 = divmod(sj.stop, fs)
-                k1 = (j1 - j0) * fs + k1
+                fs_b = attrs.get('BlockRate', [1000])[0] / 1e3
+                nsp_b = round(fs/fs_b)
+                half_ol = round((dataset.shape[1] - nsp_b) / 2)
+                j0, k0 = divmod(sj.start, nsp_b)
+                j1, k1 = divmod(sj.stop, nsp_b)
+                k1 = (j1 - j0) * nsp_b + k1
                 j1 += 1
-                data = dataset[j0:j1, :, si]
+                data = dataset[j0:j1, half_ol:half_ol+nsp_b, si]
                 data = data.reshape((-1, data.shape[-1]))[k0:k1, :].T
             else:
                 data = dataset[sj, si].T
@@ -444,9 +442,9 @@ def _read_tdms(fname, headonly=False, file_format='auto', chmin=None,
                 version = float(properties['iDASVersion'][:3])
                 if version < 2.3:
                     file_format = 'Silixa iDAS'
-                elif 2.3 <= version < 2.6:
+                elif 2.3 <= version < 2.7:
                     file_format = 'Silixa iDAS-v2'
-                elif version >= 2.6:
+                elif version >= 2.7:
                     file_format = 'Silixa iDAS-v3'
             elif group_name == ['DAS']:
                 key = 'DAS'
@@ -483,12 +481,11 @@ def _read_tdms(fname, headonly=False, file_format='auto', chmin=None,
                         if isinstance(properties[time_key], str):
                             start_time = DASDateTime.fromisoformat(
                                 properties[time_key])
+                            break
                         elif isinstance(properties[time_key], np.datetime64):
                             start_time = DASDateTime.from_datetime(
                                 properties[time_key].item())
-                        else:
-                            continue
-                        break
+                            break
                 metadata['start_time'] = start_time
             if 'GaugeLength' in properties.keys():
                 metadata['gauge_length'] = properties['GaugeLength']

{daspy_toolbox-1.2.0 → daspy_toolbox-1.2.1}/daspy/core/section.py

@@ -1,6 +1,6 @@
 # Purpose: Module for handling Section objects.
 # Author: Minzhe Hu
-# Date: 2025.9.15
+# Date: 2025.10.31
 # Email: hmz2018@mail.ustc.edu.cn
 import warnings
 import os
@@ -150,15 +150,15 @@ class Section(object):
         else:
             raise TypeError('The input should be Section or np.ndarray.')
 
-        if len(data) != self.nch:
-            if len(data[0]) == self.nch:
-                data = data.T
-            else:
-                raise ValueError('These two Sections have different number of '
-                                 'channels, please check.')
-        if out.data is None:
+        if (out.data is None) or (out.nch * out.nsp == 0):
             out.data = data
         else:
+            if len(data) != self.nch:
+                if len(data[0]) == self.nch:
+                    data = data.T
+                else:
+                    raise ValueError('These two Sections have different number '
+                                     'of channels, please check.')
             out.data = np.hstack((out.data, data))
 
         return out
@@ -548,7 +548,7 @@ class Section(object):
             self.start_distance += channel[0] * self.dx
             return self
         else:
-            return data
+            return data * self.scale
 
     def plot(self, xmode='distance', tmode='origin', obj='waveform',
              kwargs_pro={}, **kwargs):
@@ -645,7 +645,10 @@ class Section(object):
             if tmode == 'origin':
                 if hasattr(self, 'origin_time'):
                     kwargs['t0'] -= self.origin_time
-                    kwargs.setdefault('ylabel', 'Times(s) after occurance')
+                    if ('transpose' in kwargs.keys()) and kwargs['transpose']:
+                        kwargs.setdefault('xlabel', 'Times(s) after occurance')
+                    else:
+                        kwargs.setdefault('ylabel', 'Times(s) after occurance')
                 else:
                     tmode == 'start'
             if tmode == 'start':
@@ -765,18 +768,26 @@ class Section(object):
         self.data = cosine_taper(self.data, p=p, side=side)
         return self
 
-    def downsampling(self, xint=None, tint=None, stack=True,
+    def downsampling(self, xint=None, tint=None, fs=None, dx=None, stack=True,
                      lowpass_filter=True):
         """
         Downsample DAS data.
 
         :param xint: int. Spatial downsampling factor.
         :param tint: int. Time downsampling factor.
+        :param fs: float. Target sampling rate after downsampling. It is used
+            if tint is None.
+        :param dx: float. Target channel interval after downsampling. It is
+            used if xint is None.
         :param stack: bool. If True, stacking will replace decimation.
         :param lowpass_filter: bool. Lowpass cheby2 filter before time
             downsampling or not.
         :return: Downsampled data.
         """
+        if xint is None and dx is not None:
+            xint = round(dx / self.dx)
+        if tint is None and fs is not None:
+            tint = round(self.fs / fs)
         self.data = downsampling(self.data, xint=xint, tint=tint, stack=stack,
                                  lowpass_filter=lowpass_filter)
         if xint and xint > 1:
@@ -851,25 +862,26 @@ class Section(object):
                     return self
         warnings.warn('Unable to convert data type.')
 
-    def time_integration(self, c=0):
+    def time_integration(self, domain='time', c=0):
         """
         Integrate DAS data in time.
 
         :param c: float. A constant added to the result.
         """
-        self.data = time_integration(self.data, self.fs, c=c)
+        self.data = time_integration(self.data, self.fs, domain=domain, c=c)
         if hasattr(self, 'data_type'):
             self._time_int_dif_attr(mode=1)
         return self
 
-    def time_differential(self, prepend=0):
+    def time_differential(self, domain='time', prepend=0):
         """
         Differentiate DAS data in time.
 
-        :param prepend: 'mean' or values to prepend to `data` along axis prior to
-            performing the difference. 
+        :param prepend: 'mean' or values to prepend to `data` along axis prior
+            to performing the difference. 
         """
-        self.data = time_differential(self.data, self.fs, prepend=prepend)
+        self.data = time_differential(self.data, self.fs, domain=domain,
+                                      prepend=prepend)
         if hasattr(self, 'data_type'):
             self._time_int_dif_attr(mode=-1)
         return self
@@ -901,7 +913,8 @@ class Section(object):
             None.
         """
         if zi is None:
-            self.data = bandpass(self.data, self.fs, freqmin, freqmax, **kwargs)
+            self.data = bandpass(self.data, self.fs, freqmin, freqmax,
+                                 **kwargs)
             return self
         else:
             self.data, zf = bandpass(self.data, self.fs, freqmin, freqmax,

{daspy_toolbox-1.2.0 → daspy_toolbox-1.2.1}/daspy/core/write.py

@@ -1,6 +1,6 @@
 # Purpose: Module for writing DAS data.
 # Author: Minzhe Hu
-# Date: 2025.9.15
+# Date: 2025.9.18
 # Email: hmz2018@mail.ustc.edu.cn
 import os
 import warnings
@@ -168,7 +168,6 @@ def _write_h5(sec, fname, raw_fname=None, file_format='auto'):
                 start_time = sec.start_time.utc() if \
                     isinstance(sec.start_time, datetime) else \
                     DASDateTime.fromtimestamp(sec.start_time)
-
                 if file_format in ['OptaSense ODH4+', 'OptaSense QuantX',
                                    'Unknown']:
                     h5_file.get('Acquisition/Raw[0]/').create_dataset('RawData',
@@ -186,6 +185,7 @@ def _write_h5(sec, fname, raw_fname=None, file_format='auto'):
                         attrs['PartStartTime'] = stime_str
                     h5_file['Acquisition'].attrs['MeasurementStartTime'] = \
                         stime_str
+                h5_file['Acquisition'].attrs['NumberOfLoci'] = sec.nch
                 stimestamp = start_time.timestamp()
                 datatime = (np.arange(stimestamp, stimestamp + sec.nsp /
                     sec.fs, 1 / sec.fs) * 1e6).astype(int)
@@ -434,9 +434,14 @@ def _write_h5(sec, fname, raw_fname=None, file_format='auto'):
             elif file_format in ['OptaSense ODH4+', 'OptaSense QuantX',
                                 'Silixa iDAS-MG', 'Sintela Onyx v1.0',
                                 'Smart Earth ZD-DAS', 'Unknown']:
+                if file_format in ['Silixa iDAS-MG', 'Sintela Onyx v1.0',
+                                   'Smart Earth ZD-DAS']:
+                    _update_h5_dataset(h5_file, 'Acquisition/Raw[0]/',
+                                       'RawData', sec.data.T)
+                else:
+                    _update_h5_dataset(h5_file, 'Acquisition/Raw[0]/',
+                                       'RawData', sec.data)
                 h5_file['Acquisition'].attrs['NumberOfLoci'] = sec.nch
-                _update_h5_dataset(h5_file, 'Acquisition/Raw[0]/', 'RawData',
-                                    sec.data)
                 if isinstance(sec.start_time, datetime):
                     if isinstance(h5_file['Acquisition/Raw[0]/RawData'].
                                     attrs['PartStartTime'], bytes):

{daspy_toolbox-1.2.0 → daspy_toolbox-1.2.1}/setup.py

@@ -2,14 +2,14 @@ from setuptools import setup, find_packages
 
 
 setup(
-    name='DASPy-toolbox', version='1.2.0',
+    name='DASPy-toolbox', version='1.2.1',
     description=(
         '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.'
         ),
-    long_description=open('README.md').read(),
+    long_description=open('README.md', encoding="utf-8").read(),
     author='Minzhe Hu, Zefeng Li',
     author_email='hmz2018@mail.ustc.edu.cn',
     maintainer='Minzhe Hu',
@@ -17,7 +17,7 @@ setup(
     license='MIT License',
     url='https://github.com/HMZ-03/DASPy',
     packages=find_packages(),
-    entry_points={  
+    entry_points={
         'console_scripts': [
             'daspy = daspy.main:main',
             ]