17 Sep
2016
17 Sep
'16
12:49 p.m.
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
From cz_chiang(a)cugb.edu.cn Sun Sep 18 02:37:29 2016
Return-Path: <cz_chiang(a)cugb.edu.cn>
X-Original-To: msnoise(a)mailman-as.oma.be
Delivered-To: msnoise(a)mailman-as.oma.be
Received: from cugb.edu.cn (mail.cugb.edu.cn [202.204.105.6])
by mailman-as.oma.be (Postfix) with ESMTP id 7BE08200C2
for <msnoise(a)mailman-as.oma.be>be>; Sun, 18 Sep 2016 02:37:28 +0000 (UTC)
Received: by ajax-webmail-app2 (Coremail) ; Sun, 18 Sep 2016 10:43:10 +0800
(GMT+08:00)
Date: Sun, 18 Sep 2016 10:43:10 +0800 (GMT+08:00)
From: =?GBK?B?va+z0NfO?= <cz_chiang(a)cugb.edu.cn>
To: msnoise(a)mailman-as.oma.be
Message-ID: <772507fd.59e6.1573b2d0b3d.Coremail.cz_chiang(a)cugb.edu.cn>
MIME-Version: 1.0
X-Originating-IP: [121.194.89.182]
X-Priority: 3
X-Mailer: Coremail Webmail Server Version 4.0.8 dev build
20150107(58648.7033.6860) Copyright (c) 2002-2016 www.mailtech.cn cug
X-SendMailWithSms: false
X-CM-TRANSID: DAMKCgAnioo+_91XIf5yAQ--.15638W
X-CM-SenderInfo: qrqriiqurrjjo6fxwuoohg3hdfq/1tbiAQICCVPfY5cb6AAAsD
X-Coremail-Antispam: 1Ur529EdanIXcx71UUUUU7IcSsGvfJ3iIAIbVAYjsxI4VWxJw
CS07vEb4IE77IF4wCS07vE1I0E4x80FVAKz4kxMIAIbVAFxVCaYxvI4VCIwcAKzIAtYxBI
daVFxhVjvjDU=
Content-Type: text/plain; charset=GBK
Content-Transfer-Encoding: 7bit
X-Content-Filtered-By: Mailman/MimeDel 2.1.12
Subject: [MSNoise] MSNoise stack error
X-BeenThere: msnoise(a)mailman-as.oma.be
X-Mailman-Version: 2.1.12
Precedence: list
List-Id: Python Package for Monitoring Seismic Velocity Changes using Ambient
Seismic Noise <msnoise.mailman-as.oma.be>
List-Unsubscribe: <http://mailman-as.oma.be/mailman/options/msnoise>,
<mailto:msnoise-request@mailman-as.oma.be?subject=unsubscribe>
List-Archive: <http://mailman-as.oma.be/mailman/private/msnoise/>
List-Post: <mailto:msnoise@mailman-as.oma.be>
List-Help: <mailto:msnoise-request@mailman-as.oma.be?subject=help>
List-Subscribe: <http://mailman-as.oma.be/mailman/listinfo/msnoise>,
<mailto:msnoise-request@mailman-as.oma.be?subject=subscribe>
X-List-Received-Date: Sun, 18 Sep 2016 02:37:29 -0000
Hello,
Sorry for the terrible format of the last e-mail I sent. I came across a stack error when
using MSNoise and could not find the solution in the documentation. I have finished the
computation of cross-correlations for my data. Now I would like to do 60-day stack with
10-day intervals. However, It returned wrong results when I tried to run 'msnoise
stack -s -i 999'. The process is pasted below. The result on 2009-09-10 is correct,
but the remaining stacks are exact the same with 2009-09-10. I have no idea what is wrong.
How can I resolve it?
Thank you in advance,
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>msnoise stack -s -i 999
Lets STACK !
2016-09-16 10:14:40 [DEBUG] Starting the step stack
[]
2016-09-16 10:14:40 [DEBUG] Processing YP_NE41:YP_NE64-ZZ-1
2016-09-16 10:14:40 [DEBUG] Found 650 updated days
2016-09-16 10:14:40 [DEBUG] New Data for YP_NE41:YP_NE64-ZZ-1
2016-09-16 10:14:44 [DEBUG] YP_NE41_YP_NE64 2009-09-10 [2009-09-10 - 2009-11-08] (60 day
stack)
e:\anaconda\anaconda2\lib\site-packages\obspy\io\sac\util.py:371: UserWarning: Old header
has invalid reftime.
warnings.warn(msg)
2016-09-16 10:14:44 [DEBUG] YP_NE41_YP_NE64 2009-11-09 [2009-09-10 - 2009-11-08] (60 day
stack)
2016-09-16 10:14:44 [DEBUG] YP_NE41_YP_NE64 2010-01-08 [2009-09-10 - 2009-11-08] (60 day
stack)
2016-09-16 10:14:45 [DEBUG] YP_NE41_YP_NE64 2010-03-09 [2009-09-10 - 2009-11-08] (60 day
stack)
2016-09-16 10:14:45 [DEBUG] YP_NE41_YP_NE64 2010-05-08 [2009-09-10 - 2009-11-08] (60 day
stack)
2016-09-16 10:14:45 [DEBUG] YP_NE41_YP_NE64 2010-07-07 [2009-09-10 - 2009-11-08] (60 day
stack)
2016-09-16 10:14:45 [DEBUG] YP_NE41_YP_NE64 2010-09-05 [2009-09-10 - 2009-11-08] (60 day
stack)
2016-09-16 10:14:46 [DEBUG] YP_NE41_YP_NE64 2010-11-04 [2009-09-10 - 2009-11-08] (60 day
stack)
2016-09-16 10:14:46 [DEBUG] YP_NE41_YP_NE64 2011-01-03 [2009-09-10 - 2009-11-08] (60 day
stack)
2016-09-16 10:14:46 [DEBUG] YP_NE41_YP_NE64 2011-03-04 [2009-09-10 - 2009-11-08] (60 day
stack)
2016-09-16 10:14:46 [DEBUG] YP_NE41_YP_NE64 2011-05-03 [2009-09-10 - 2009-11-08] (60 day
stack)
2016-09-16 10:14:46 [DEBUG] YP_NE41_YP_NE64 2011-07-02 [2009-09-10 - 2009-11-08] (60 day
stack)
From hsy119(a)mail.ustc.edu.cn Tue Sep 20 02:10:23 2016
Return-Path: <hsy119(a)mail.ustc.edu.cn>
X-Original-To: msnoise(a)mailman-as.oma.be
Delivered-To: msnoise(a)mailman-as.oma.be
Received: from ustc.edu.cn (smtp2.ustc.edu.cn [202.38.64.46])
by mailman-as.oma.be (Postfix) with ESMTP id 2A38F2003C
for <msnoise(a)mailman-as.oma.be>be>; Tue, 20 Sep 2016 02:10:21 +0000 (UTC)
Received: by ajax-webmail-newmailweb.ustc.edu.cn (Coremail) ; Tue, 20 Sep
2016 10:10:18 +0800 (GMT+08:00)
X-CM-HeaderCharset: UTF-8
X-Originating-IP: [114.214.191.239]
Date: Tue, 20 Sep 2016 10:10:18 +0800 (GMT+08:00)
From: hsy119(a)mail.ustc.edu.cn
To: msnoise(a)mailman-as.oma.be
X-Priority: 3
X-Mailer: Coremail Webmail Server Version XT3.0.6a_preview build
20150605(69773.7394) Copyright (c) 2002-2016 www.mailtech.cn ustc-xl
X-SendMailWithSms: false
MIME-Version: 1.0
Message-ID: <477c397d.1cc5c.157455bacbd.Coremail.hsy119(a)mail.ustc.edu.cn>
X-CM-TRANSID: LkAmygDXh_yKmuBXgleTCA--.37168W
X-CM-SenderInfo: 5kv1ii2z6ptxnoox23vfohv3gofq/1tbiAQYEAFQhmNIukQAAsD
X-Coremail-Antispam: 1Ur529EdanIXcx71UUUUU7IcSsGvfJ3iIAIbVAYjsxI4VWxJw
CS07vEb4IE77IF4wCS07vE1I0E4x80FVAKz4kxMIAIbVAFxVCaYxvI4VCIwcAKzIAtYxBI
daVFxhVjvjDU=
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: base64
X-Content-Filtered-By: Mailman/MimeDel 2.1.12
Subject: [MSNoise] Error when plot dtt figure
X-BeenThere: msnoise(a)mailman-as.oma.be
X-Mailman-Version: 2.1.12
Precedence: list
List-Id: Python Package for Monitoring Seismic Velocity Changes using Ambient
Seismic Noise <msnoise.mailman-as.oma.be>
List-Unsubscribe: <http://mailman-as.oma.be/mailman/options/msnoise>,
<mailto:msnoise-request@mailman-as.oma.be?subject=unsubscribe>
List-Archive: <http://mailman-as.oma.be/mailman/private/msnoise/>
List-Post: <mailto:msnoise@mailman-as.oma.be>
List-Help: <mailto:msnoise-request@mailman-as.oma.be?subject=help>
List-Subscribe: <http://mailman-as.oma.be/mailman/listinfo/msnoise>,
<mailto:msnoise-request@mailman-as.oma.be?subject=subscribe>
X-List-Received-Date: Tue, 20 Sep 2016 02:10:23 -0000
SGksCgpTb3JyeSB0byBib3RoZXIgYWdhaW4uCkknbSB0cnlpbmcgdG8gcGxvdCBkdHQgZmlndXJl
IGJ1dCBJIG1lZXQgYW4gZXJyb3IuIEl0J3Mgc3RyYW5nZSBiZWNhdXNlIHRoZXJlJ3Mgbm8gcHJv
YmxlbSB3aGVuIEkgcGxvdCBhIGNlcnRhaW4gc3RhdGlvbiBwYWlyLCB3aGF0ZXZlciB0aGUgZGF0
YSBpcy4gQnV0IGZvciB0aGUgcmVzZXQgcGFpcnMsIHRoZXJlJ3MgYSBidWcgYXMgZm9sbG93cy4K
CgogICAgICAgICBEYXRlICAgICAgICAgQSAgICAgICAgRUEgICAgICAgIEVNICAgICAgIEVNMCAg
ICAgICAgIE0gICAgICAgIE0wICBcCjEgIDIwMTUtMDktMDEgIDAuMDEwODYyICAwLjAxMDMwNiAg
MC4wMDA3MDggIDAuMDAwNDQ1IC0wLjAwMDE5OCAgMC4wMDAzODMgIAoKICAgICAgICAgICAgIFBh
aXJzIAoxICBHMl8xNTkyX0cyXzE2NTAgClRyYWNlYmFjayAobW9zdCByZWNlbnQgY2FsbCBsYXN0
KToKICBGaWxlICIvcHVibGljL3NvZnR3YXJlL2FuYWNvbmRhMy9iaW4vbXNub2lzZSIsIGxpbmUg
MTEsIGluIDxtb2R1bGU+CiAgICBzeXMuZXhpdChydW4oKSkKICBGaWxlICIvcHVibGljL3NvZnR3
YXJlL2FuYWNvbmRhMy9saWIvcHl0aG9uMy41L3NpdGUtcGFja2FnZXMvbXNub2lzZS9zY3JpcHRz
L21zbm9pc2UucHkiLCBsaW5lIDYxNCwgaW4gcnVuCiAgICBjbGkob2JqPXt9KQogIEZpbGUgIi9w
dWJsaWMvc29mdHdhcmUvYW5hY29uZGEzL2xpYi9weXRob24zLjUvc2l0ZS1wYWNrYWdlcy9jbGlj
ay9jb3JlLnB5IiwgbGluZSA3MTYsIGluIF9fY2FsbF9fCiAgICByZXR1cm4gc2VsZi5tYWluKCph
cmdzLCAqKmt3YXJncykKICBGaWxlICIvcHVibGljL3NvZnR3YXJlL2FuYWNvbmRhMy9saWIvcHl0
aG9uMy41L3NpdGUtcGFja2FnZXMvY2xpY2svY29yZS5weSIsIGxpbmUgNjk2LCBpbiBtYWluCiAg
ICBydiA9IHNlbGYuaW52b2tlKGN0eCkKICBGaWxlICIvcHVibGljL3NvZnR3YXJlL2FuYWNvbmRh
My9saWIvcHl0aG9uMy41L3NpdGUtcGFja2FnZXMvY2xpY2svY29yZS5weSIsIGxpbmUgMTA2MCwg
aW4gaW52b2tlCiAgICByZXR1cm4gX3Byb2Nlc3NfcmVzdWx0KHN1Yl9jdHguY29tbWFuZC5pbnZv
a2Uoc3ViX2N0eCkpCiAgRmlsZSAiL3B1YmxpYy9zb2Z0d2FyZS9hbmFjb25kYTMvbGliL3B5dGhv
bjMuNS9zaXRlLXBhY2thZ2VzL2NsaWNrL2NvcmUucHkiLCBsaW5lIDEwNjAsIGluIGludm9rZQog
ICAgcmV0dXJuIF9wcm9jZXNzX3Jlc3VsdChzdWJfY3R4LmNvbW1hbmQuaW52b2tlKHN1Yl9jdHgp
KQogIEZpbGUgIi9wdWJsaWMvc29mdHdhcmUvYW5hY29uZGEzL2xpYi9weXRob24zLjUvc2l0ZS1w
YWNrYWdlcy9jbGljay9jb3JlLnB5IiwgbGluZSA4ODksIGluIGludm9rZQogICAgcmV0dXJuIGN0
eC5pbnZva2Uoc2VsZi5jYWxsYmFjaywgKipjdHgucGFyYW1zKQogIEZpbGUgIi9wdWJsaWMvc29m
dHdhcmUvYW5hY29uZGEzL2xpYi9weXRob24zLjUvc2l0ZS1wYWNrYWdlcy9jbGljay9jb3JlLnB5
IiwgbGluZSA1MzQsIGluIGludm9rZQogICAgcmV0dXJuIGNhbGxiYWNrKCphcmdzLCAqKmt3YXJn
cykKICBGaWxlICIvcHVibGljL3NvZnR3YXJlL2FuYWNvbmRhMy9saWIvcHl0aG9uMy41L3NpdGUt
cGFja2FnZXMvY2xpY2svZGVjb3JhdG9ycy5weSIsIGxpbmUgMTcsIGluIG5ld19mdW5jCiAgICBy
ZXR1cm4gZihnZXRfY3VycmVudF9jb250ZXh0KCksICphcmdzLCAqKmt3YXJncykKICBGaWxlICIv
cHVibGljL3NvZnR3YXJlL2FuYWNvbmRhMy9saWIvcHl0aG9uMy41L3NpdGUtcGFja2FnZXMvbXNu
b2lzZS9zY3JpcHRzL21zbm9pc2UucHkiLCBsaW5lIDU1OSwgaW4gZHR0CiAgICBtYWluKHN0YTEs
IHN0YTIsIGZpbHRlcmlkLCBjb21wLCBkYXksIG1vdl9zdGFjaywgc2hvdywgb3V0ZmlsZSkKICBG
aWxlICIvcHVibGljL3NvZnR3YXJlL2FuYWNvbmRhMy9saWIvcHl0aG9uMy41L3NpdGUtcGFja2Fn
ZXMvbXNub2lzZS9wbG90cy9kdHQucHkiLCBsaW5lIDkxLCBpbiBtYWluCiAgICBwbHQucGxvdCh4
bGluZU0wLCB5bGluZU0wLCAncicsIGxhYmVsPSdNMCcpCiAgRmlsZSAiL3B1YmxpYy9zb2Z0d2Fy
ZS9hbmFjb25kYTMvbGliL3B5dGhvbjMuNS9zaXRlLXBhY2thZ2VzL21hdHBsb3RsaWIvcHlwbG90
LnB5IiwgbGluZSAzMTU0LCBpbiBwbG90CiAgICByZXQgPSBheC5wbG90KCphcmdzLCAqKmt3YXJn
cykKICBGaWxlICIvcHVibGljL3NvZnR3YXJlL2FuYWNvbmRhMy9saWIvcHl0aG9uMy41L3NpdGUt
cGFja2FnZXMvbWF0cGxvdGxpYi9fX2luaXRfXy5weSIsIGxpbmUgMTgxMiwgaW4gaW5uZXIKICAg
IHJldHVybiBmdW5jKGF4LCAqYXJncywgKiprd2FyZ3MpCiAgRmlsZSAiL3B1YmxpYy9zb2Z0d2Fy
ZS9hbmFjb25kYTMvbGliL3B5dGhvbjMuNS9zaXRlLXBhY2thZ2VzL21hdHBsb3RsaWIvYXhlcy9f
YXhlcy5weSIsIGxpbmUgMTQyNCwgaW4gcGxvdAogICAgZm9yIGxpbmUgaW4gc2VsZi5fZ2V0X2xp
bmVzKCphcmdzLCAqKmt3YXJncyk6CiAgRmlsZSAiL3B1YmxpYy9zb2Z0d2FyZS9hbmFjb25kYTMv
bGliL3B5dGhvbjMuNS9zaXRlLXBhY2thZ2VzL21hdHBsb3RsaWIvYXhlcy9fYmFzZS5weSIsIGxp
bmUgMzg2LCBpbiBfZ3JhYl9uZXh0X2FyZ3MKICAgIGZvciBzZWcgaW4gc2VsZi5fcGxvdF9hcmdz
KHJlbWFpbmluZywga3dhcmdzKToKICBGaWxlICIvcHVibGljL3NvZnR3YXJlL2FuYWNvbmRhMy9s
aWIvcHl0aG9uMy41L3NpdGUtcGFja2FnZXMvbWF0cGxvdGxpYi9heGVzL19iYXNlLnB5IiwgbGlu
ZSAzNjAsIGluIF9wbG90X2FyZ3MKICAgIHkgPSBfY2hlY2tfMWQodHVwWy0xXSkKICBGaWxlICIv
cHVibGljL3NvZnR3YXJlL2FuYWNvbmRhMy9saWIvcHl0aG9uMy41L3NpdGUtcGFja2FnZXMvbWF0
cGxvdGxpYi9jYm9vay5weSIsIGxpbmUgMjIwNywgaW4gX2NoZWNrXzFkCiAgICByZXR1cm4gbnAu
YXRsZWFzdF8xZCh4KQogIEZpbGUgIi9wdWJsaWMvc29mdHdhcmUvYW5hY29uZGEzL2xpYi9weXRo
b24zLjUvc2l0ZS1wYWNrYWdlcy9udW1weS9jb3JlL3NoYXBlX2Jhc2UucHkiLCBsaW5lIDUwLCBp
biBhdGxlYXN0XzFkCiAgICBhcnkgPSBhc2FueWFycmF5KGFyeSkKICBGaWxlICIvcHVibGljL3Nv
ZnR3YXJlL2FuYWNvbmRhMy9saWIvcHl0aG9uMy41L3NpdGUtcGFja2FnZXMvbnVtcHkvY29yZS9u
dW1lcmljLnB5IiwgbGluZSA1MzMsIGluIGFzYW55YXJyYXkKICAgIHJldHVybiBhcnJheShhLCBk
dHlwZSwgY29weT1GYWxzZSwgb3JkZXI9b3JkZXIsIHN1Ym9rPVRydWUpCiAgRmlsZSAiL3B1Ymxp
Yy9zb2Z0d2FyZS9hbmFjb25kYTMvbGliL3B5dGhvbjMuNS9zaXRlLXBhY2thZ2VzL3BhbmRhcy9j
b3JlL3Nlcmllcy5weSIsIGxpbmUgNTgzLCBpbiBfX2dldGl0ZW1fXwogICAgcmVzdWx0ID0gc2Vs
Zi5pbmRleC5nZXRfdmFsdWUoc2VsZiwga2V5KQogIEZpbGUgIi9wdWJsaWMvc29mdHdhcmUvYW5h
Y29uZGEzL2xpYi9weXRob24zLjUvc2l0ZS1wYWNrYWdlcy9wYW5kYXMvaW5kZXhlcy9iYXNlLnB5
IiwgbGluZSAxOTgwLCBpbiBnZXRfdmFsdWUKICAgIHR6PWdldGF0dHIoc2VyaWVzLmR0eXBlLCAn
dHonLCBOb25lKSkKICBGaWxlICJwYW5kYXMvaW5kZXgucHl4IiwgbGluZSAxMDMsIGluIHBhbmRh
cy5pbmRleC5JbmRleEVuZ2luZS5nZXRfdmFsdWUgKHBhbmRhcy9pbmRleC5jOjMzMzIpCiAgRmls
ZSAicGFuZGFzL2luZGV4LnB5eCIsIGxpbmUgMTExLCBpbiBwYW5kYXMuaW5kZXguSW5kZXhFbmdp
bmUuZ2V0X3ZhbHVlIChwYW5kYXMvaW5kZXguYzozMDM1KQogIEZpbGUgInBhbmRhcy9pbmRleC5w
eXgiLCBsaW5lIDE1OSwgaW4gcGFuZGFzLmluZGV4LkluZGV4RW5naW5lLmdldF9sb2MgKHBhbmRh
cy9pbmRleC5jOjQwMTgpCiAgRmlsZSAicGFuZGFzL2hhc2h0YWJsZS5weXgiLCBsaW5lIDMwMywg
aW4gcGFuZGFzLmhhc2h0YWJsZS5JbnQ2NEhhc2hUYWJsZS5nZXRfaXRlbSAocGFuZGFzL2hhc2h0
YWJsZS5jOjY2MTApCiAgRmlsZSAicGFuZGFzL2hhc2h0YWJsZS5weXgiLCBsaW5lIDMwOSwgaW4g
cGFuZGFzLmhhc2h0YWJsZS5JbnQ2NEhhc2hUYWJsZS5nZXRfaXRlbSAocGFuZGFzL2hhc2h0YWJs
ZS5jOjY1NTQpCgpUaGVyZSdzIG5vIGVycm9yIHdoZW4gSSBjb21wdXRlIGNjLCBtd2NzLCBkdHQs
b25seSBzb21lIHdhcm5pbmcuClNvIGhvdyBjYW4gSSBmaXggdGhpcz8KVGhhbmtzIGZvciB5b3Vy
IGhlbHAKClNodXllCgoK
From Thomas.Lecocq(a)seismology.be Tue Sep 20 10:30:15
2016
Return-Path: <Thomas.Lecocq(a)seismology.be>
X-Original-To: msnoise(a)mailman-as.oma.be
Delivered-To: msnoise(a)mailman-as.oma.be
Received: from 9.mo4.mail-out.ovh.net (9.mo4.mail-out.ovh.net [46.105.40.176])
by mailman-as.oma.be (Postfix) with ESMTP id 9D8E6205E6
for <msnoise(a)mailman-as.oma.be>be>; Tue, 20 Sep 2016 10:30:15 +0000 (UTC)
Received: from player691.ha.ovh.net (b9.ovh.net [213.186.33.59])
by mo4.mail-out.ovh.net (Postfix) with ESMTP id 6C2D6FFAC0B
for <msnoise(a)mailman-as.oma.be>be>; Tue, 20 Sep 2016 12:30:15 +0200 (CEST)
Received: from [10.10.50.104] (unknown [83.240.218.218])
(Authenticated sender: tom(a)asktom.be)
by player691.ha.ovh.net (Postfix) with ESMTPSA id 3DF52260095
for <msnoise(a)mailman-as.oma.be>be>; Tue, 20 Sep 2016 12:30:14 +0200 (CEST)
To: msnoise(a)mailman-as.oma.be
References: <477c397d.1cc5c.157455bacbd.Coremail.hsy119(a)mail.ustc.edu.cn>
From: Thomas Lecocq <Thomas.Lecocq(a)seismology.be>
Message-ID: <04b37855-7e20-f38a-2525-28ebfabbbf53(a)seismology.be>
Date: Tue, 20 Sep 2016 11:30:05 +0100
User-Agent: Mozilla/5.0 (Windows NT 10.0; WOW64; rv:45.0) Gecko/20100101
Thunderbird/45.2.0
MIME-Version: 1.0
In-Reply-To: <477c397d.1cc5c.157455bacbd.Coremail.hsy119(a)mail.ustc.edu.cn>
Content-Type: text/plain; charset=windows-1252; format=flowed
Content-Transfer-Encoding: 7bit
X-Ovh-Tracer-Id: 1132373832982925598
X-VR-SPAMSTATE: OK
X-VR-SPAMSCORE: 0
X-VR-SPAMCAUSE:
gggruggvucftvghtrhhoucdtuddrfeeluddrkeeigddtleculddtuddrfeeltddrtddtmdcutefuodetggdotefrodftvfcurfhrohhfihhlvgemucfqggfjnecuuegrihhlohhuthemuceftddtnecu
Subject: Re: [MSNoise] Error when plot dtt figure
X-BeenThere: msnoise(a)mailman-as.oma.be
X-Mailman-Version: 2.1.12
Precedence: list
List-Id: Python Package for Monitoring Seismic Velocity Changes using Ambient
Seismic Noise <msnoise.mailman-as.oma.be>
List-Unsubscribe: <http://mailman-as.oma.be/mailman/options/msnoise>,
<mailto:msnoise-request@mailman-as.oma.be?subject=unsubscribe>
List-Archive: <http://mailman-as.oma.be/mailman/private/msnoise/>
List-Post: <mailto:msnoise@mailman-as.oma.be>
List-Help: <mailto:msnoise-request@mailman-as.oma.be?subject=help>
List-Subscribe: <http://mailman-as.oma.be/mailman/listinfo/msnoise>,
<mailto:msnoise-request@mailman-as.oma.be?subject=subscribe>
X-List-Received-Date: Tue, 20 Sep 2016 10:30:15 -0000
Hi Shuye,
I don't understand what you mean with "the reset pairs" ? What is the
command you used to start this plot ?
Thomas
On 20/09/2016 03:10, hsy119(a)mail.ustc.edu.cn wrote:
Hi,
Sorry to bother again.
I'm trying to plot dtt figure but I meet an error. It's strange because
there's no problem when I plot a certain station pair, whatever the data is. But for
the reset pairs, there's a bug as follows.
Date A EA EM EM0 M M0 \
1 2015-09-01 0.010862 0.010306 0.000708 0.000445 -0.000198 0.000383
Pairs
1 G2_1592_G2_1650
Traceback (most recent call last):
File "/public/software/anaconda3/bin/msnoise", line 11, in <module>
sys.exit(run())
File
"/public/software/anaconda3/lib/python3.5/site-packages/msnoise/scripts/msnoise.py",
line 614, in run
cli(obj={})
File "/public/software/anaconda3/lib/python3.5/site-packages/click/core.py",
line 716, in __call__
return self.main(*args, **kwargs)
File "/public/software/anaconda3/lib/python3.5/site-packages/click/core.py",
line 696, in main
rv = self.invoke(ctx)
File "/public/software/anaconda3/lib/python3.5/site-packages/click/core.py",
line 1060, in invoke
return _process_result(sub_ctx.command.invoke(sub_ctx))
File "/public/software/anaconda3/lib/python3.5/site-packages/click/core.py",
line 1060, in invoke
return _process_result(sub_ctx.command.invoke(sub_ctx))
File "/public/software/anaconda3/lib/python3.5/site-packages/click/core.py",
line 889, in invoke
return ctx.invoke(self.callback, **ctx.params)
File "/public/software/anaconda3/lib/python3.5/site-packages/click/core.py",
line 534, in invoke
return callback(*args, **kwargs)
File
"/public/software/anaconda3/lib/python3.5/site-packages/click/decorators.py",
line 17, in new_func
return f(get_current_context(), *args, **kwargs)
File
"/public/software/anaconda3/lib/python3.5/site-packages/msnoise/scripts/msnoise.py",
line 559, in dtt
main(sta1, sta2, filterid, comp, day, mov_stack, show, outfile)
File
"/public/software/anaconda3/lib/python3.5/site-packages/msnoise/plots/dtt.py",
line 91, in main
plt.plot(xlineM0, ylineM0, 'r', label='M0')
File
"/public/software/anaconda3/lib/python3.5/site-packages/matplotlib/pyplot.py",
line 3154, in plot
ret = ax.plot(*args, **kwargs)
File
"/public/software/anaconda3/lib/python3.5/site-packages/matplotlib/__init__.py",
line 1812, in inner
return func(ax, *args, **kwargs)
File
"/public/software/anaconda3/lib/python3.5/site-packages/matplotlib/axes/_axes.py",
line 1424, in plot
for line in self._get_lines(*args, **kwargs):
File
"/public/software/anaconda3/lib/python3.5/site-packages/matplotlib/axes/_base.py",
line 386, in _grab_next_args
for seg in self._plot_args(remaining, kwargs):
File
"/public/software/anaconda3/lib/python3.5/site-packages/matplotlib/axes/_base.py",
line 360, in _plot_args
y = _check_1d(tup[-1])
File
"/public/software/anaconda3/lib/python3.5/site-packages/matplotlib/cbook.py",
line 2207, in _check_1d
return np.atleast_1d(x)
File
"/public/software/anaconda3/lib/python3.5/site-packages/numpy/core/shape_base.py",
line 50, in atleast_1d
ary = asanyarray(ary)
File
"/public/software/anaconda3/lib/python3.5/site-packages/numpy/core/numeric.py",
line 533, in asanyarray
return array(a, dtype, copy=False, order=order, subok=True)
File
"/public/software/anaconda3/lib/python3.5/site-packages/pandas/core/series.py",
line 583, in __getitem__
result = self.index.get_value(self, key)
File
"/public/software/anaconda3/lib/python3.5/site-packages/pandas/indexes/base.py",
line 1980, in get_value
tz=getattr(series.dtype, 'tz', None))
File "pandas/index.pyx", line 103, in pandas.index.IndexEngine.get_value
(pandas/index.c:3332)
File "pandas/index.pyx", line 111, in pandas.index.IndexEngine.get_value
(pandas/index.c:3035)
File "pandas/index.pyx", line 159, in pandas.index.IndexEngine.get_loc
(pandas/index.c:4018)
File "pandas/hashtable.pyx", line 303, in
pandas.hashtable.Int64HashTable.get_item (pandas/hashtable.c:6610)
File "pandas/hashtable.pyx", line 309, in
pandas.hashtable.Int64HashTable.get_item (pandas/hashtable.c:6554)
There's no error when I compute cc, mwcs, dtt,only some warning.
So how can I fix this?
Thanks for your help
Shuye
_______________________________________________
MSNoise mailing list
MSNoise(a)mailman-as.oma.be
http://mailman-as.oma.be/mailman/listinfo/msnoise
20 Sep
20 Sep
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
3147
days inactive
3150
days old
1 comments
2 participants
participants (2)
-
Thomas Lecocq -
蒋承孜