-
Notifications
You must be signed in to change notification settings - Fork 21
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Can't read multiple SAC files using read_data
#88
Comments
This is being caused by #87 but s = read_data("sac", "./*.SAC")
SeisData with 1 channels (1 shown)
ID: C.FR..BH
NAME:
LOC: 63.7469 N, -68.5451 E, 25.0 m
FS: 40.0
GAIN: 1.0
RESP: a0 1.0, f0 1.0, 0z, 0p
UNITS:
SRC: /home/timclements/Downloads/SACTE…
MISC: 0 entries
NOTES: 1 entries
T: 2019-09-19T07:06:33 (2 gaps)
X: -6.722e+04
-6.701e+04
...
-2.850e+04
(nx = 432000)
C: 0 open, 0 total there are two gaps. We can check out the data gaps with timing information s.t
1-element Array{Array{Int64,2},1}:
[1 1568876793300000; 144001 -3600000000; 288001 -3600000000; 432000 0] each data segment has 144,000 samples or 1 hour of data. We can check that with some helper functions: SeisIO.u2d.(SeisIO.t_win(s.t[1],s.fs[1]) .* SeisIO.μs)
3×2 Array{Dates.DateTime,2}:
2019-09-19T07:06:33.3 2019-09-19T08:06:33.275
2019-09-19T07:06:33.3 2019-09-19T08:06:33.275
2019-09-19T07:06:33.3 2019-09-19T08:06:33.275 Since all three channels have the same name and start/end times, calling ungap(s)
SeisData with 1 channels (1 shown)
ID: C.FR..BH
NAME:
LOC: 63.7469 N, -68.5451 E, 25.0 m
FS: 40.0
GAIN: 1.0
RESP: a0 1.0, f0 1.0, 0z, 0p
UNITS:
SRC: /home/timclements/Downloads/SACTE…
MISC: 0 entries
NOTES: 2 entries
T: 2019-09-19T07:06:33 (0 gaps)
X: -3.430e+04
-3.494e+04
...
-2.850e+04
(nx = 144000)
C: 0 open, 0 total which it does. Solving #87 will solve this issue. |
Duplicate of issue #87; leaving open while I investigate. If my guess is right, the fix is easy, but I need to do some further checking before I push code live (and I'm still looking for a suitable Travis-CI replacement, whose greed is ultimately the cause of these delays)... |
I'm unable to read multiple SAC channels into a single SeisData instance using
read_data
. Tried on v1.2.1 and also latest master. The files are the same as uploaded for #87 .The text was updated successfully, but these errors were encountered: