20 Sep
2016
20 Sep
'16
12:34 p.m.
Hi Chris,
The design of MSNoise is that you do the moving stack of 60 days
(=smoothing over the last 59 days) and compare that to a REF. This is
achieved with "mov stack"= "msnoise stack -m". If you want to make the
step stacking , i.e. only exporting the stack for each nth day, you need
to run the "msnoise stack -s" command, but it's really not well
documented and implemented, e.g. the dvv plots will probably fail.
The fact that the first 60 days are the same are a choice that goes back
to the orginal codes I implemented, assuming the stack should be the
same as long as the number of days of movstacks are not available. It's
a bad implementation and I really would like to change this in the future.
The "-i 9999" command is only useful if you want to re-stack data
computed more than 1 day ago.
Let me know how it goes
Thomas
On 17/09/2016 11:49, 蒋承孜 wrote:
Hello,
I came across a stack error when using MSNoise and could not find the solution in the
documentation.
Ihavefinishedthecomputationofthecross-correlationsformydataandnowIwouldliketo
do60-daystack with10-dayintervals. However, itreturned wrongresultswhenItriedtorun
'msnoisestack -s -i999'. Forinstance,
Iwouldliketostackcross-correlationfunctionscomputedfrom2009-09-10to2011-07-02.
Theresulton2009-09-10wascorrect, buttheremainingstacks were exactthesame with2009-09-10.
Theprocesswas pastedbelow. Iamnotsurewhatiswrong. HowcanI resolveit?
Thankyouinadvance,
Chris Chiang
Undergraduate Student
School of Geophysics and Information Technology
China University of Geosciences, Beijing
29 Xueyuan Road, Haidian District, Beijing, China
PROCESS:
(E:\Anaconda\Anaconda2) E:\Data\Mountain\filter0.1-0.2>msnoisestack -s -i999
LetsSTACK !
2016-09-1610:14:40 [DEBUG] Startingthestepstack
[]
2016-09-1610:14:40 [DEBUG] ProcessingYP_NE41:YP_NE64-ZZ-1
2016-09-1610:14:40 [DEBUG] Found650updateddays
2016-09-1610:14:40 [DEBUG] NewDataforYP_NE41:YP_NE64-ZZ-1
2016-09-1610:14:44 [DEBUG] YP_NE41_YP_NE642009-09-10 [2009-09-10 - 2009-11-08]
(60daystack)
e:\anaconda\anaconda2\lib\site-packages\obspy\io\sac\util.py:371: UserWarning:
Oldheaderhasinvalidreftime.
warnings.warn(msg)
2016-09-1610:14:44 [DEBUG] YP_NE41_YP_NE642009-11-09 [2009-09-10 - 2009-11-08]
(60daystack)
2016-09-1610:14:44 [DEBUG] YP_NE41_YP_NE642010-01-08 [2009-09-10 - 2009-11-08]
(60daystack)
2016-09-1610:14:45 [DEBUG] YP_NE41_YP_NE642010-03-09 [2009-09-10 - 2009-11-08]
(60daystack)
2016-09-1610:14:45 [DEBUG] YP_NE41_YP_NE642010-05-08 [2009-09-10 - 2009-11-08]
(60daystack)
2016-09-1610:14:45 [DEBUG] YP_NE41_YP_NE642010-07-07 [2009-09-10 - 2009-11-08]
(60daystack)
2016-09-1610:14:45 [DEBUG] YP_NE41_YP_NE642010-09-05 [2009-09-10 - 2009-11-08]
(60daystack)
2016-09-1610:14:46 [DEBUG] YP_NE41_YP_NE642010-11-04 [2009-09-10 - 2009-11-08]
(60daystack)
2016-09-1610:14:46 [DEBUG] YP_NE41_YP_NE642011-01-03 [2009-09-10 - 2009-11-08]
(60daystack)
2016-09-1610:14:46 [DEBUG] YP_NE41_YP_NE642011-03-04 [2009-09-10 - 2009-11-08]
(60daystack)
2016-09-1610:14:46 [DEBUG] YP_NE41_YP_NE642011-05-03 [2009-09-10 - 2009-11-08]
(60daystack)
2016-09-1610:14:46 [DEBUG] YP_NE41_YP_NE642011-07-02 [2009-09-10 - 2009-11-08]
(60daystack
_______________________________________________
MSNoise mailing list
MSNoise(a)mailman-as.oma.be
http://mailman-as.oma.be/mailman/listinfo/msnoise