U f{j@sJddlZddlmZddlmZddlmZmZm Z m Z e dddZ dS) N) Polynomial)BaseRaw)_validate_typewarnloggerverbosec Csddlm}t|tdt|tdt|t}t|t}|jdksR|j|jkrjt d|jd|jt |dkr~t d t j ||dd }|jd d }|j\}} td |d| |||\} } d| dd| d} | dks| dkrt | d| dkrt | dn t| dtd| } tdd| dd|jddd|jd| ddd|dd} |dkrt| d||d ||8}n$t| d|| d ||7}td!|jd"| }|j|d#d$|j|jd"|jd"<W5QRX|jd|jd}dt|dd%} |dkrft| d|d|jdn*|dkrt| d|d|jdd S)&a)Realign two simultaneous recordings. Due to clock drift, recordings at a given same sample rate made by two separate devices simultaneously can become out of sync over time. This function uses event times captured by both acquisition devices to resample ``other`` to match ``raw``. Parameters ---------- raw : instance of Raw The first raw instance. other : instance of Raw The second raw instance. It will be resampled to match ``raw``. t_raw : array-like, shape (n_events,) The times of shared events in ``raw`` relative to ``raw.times[0]`` (0). Typically these could be events on some TTL channel like ``find_events(raw)[:, 0] - raw.first_event``. t_other : array-like, shape (n_events,) The times of shared events in ``other`` relative to ``other.times[0]``. %(verbose)s Notes ----- This function operates inplace. It will: 1. Estimate the zero-order (start offset) and first-order (clock drift) correction. 2. Crop the start of ``raw`` or ``other``, depending on which started recording first. 3. Resample ``other`` to match ``raw`` based on the clock drift. 4. Crop the end of ``raw`` or ``other``, depending on which stopped recording first (and the clock drift rate). This function is primarily designed to work on recordings made at the same sample rate, but it can also operate on recordings made at different sample rates to resample and deal with clock drift simultaneously. .. versionadded:: 0.22 r)statsrawotherz=t_raw and t_other must be 1D with the same shape, got shapes z and z5Fewer than 20 times passed, results may be unreliable)xydeg)r )domainzZero order coefficient: z First order coefficient: z!Linear correlation computed as R=z0.3fz and p=z0.2eg?z, cannot resample safelygư>z, results may be unreliableiz Drift rate: z0.1fu μs/sec (total drift over rz sec recording: z ms)z Cropping z sec from the start of NzResampling otherZsfreqT)rz sec from the end of )Zscipyr rrnparrayfloatndimshape ValueErrorlenrrfitconvertZcoefrinfoZpearsonrabstimescropZ load_dataZresampleZ_unlock)r r Zt_rawZt_otherrr ZpolyZ convertedZzero_ordZ first_ordrpmsgZdr_ms_sZ sfreq_newdeltar$X/var/www/html/bwave/python/venv/lib/python3.8/site-packages/mne/preprocessing/realign.py realign_raw sZ)           0       r&)N) numpyrZnumpy.polynomial.polynomialriorutilsrrrrr&r$r$r$r%s