[MSNoise] Compute MWCS

[Thomas Lecocq]

Dear Weiyu,

It looks like your filters are really narrow band, which will make very 
few points available for the linear regression. I suspect this is the 
reason for the warnings & the NaN outputs.

Cheers,

Thomas

Le 10/06/2019 à 10:50, weiyu a écrit :
> Hi MSNoise group,
> I recently ran into a problem when compute MWCS. I have 9 filters in my project, when I ran “msnoise compute-mwcs”,it worked successfully, but when calculated filter 8 and filter 9 something happened, the parameters were “0.04 0.05 0.0667 0.1” and “0.025 0.033 0.05 0.1” respectively. The code reported some Warnings:
>
> 2019-06-10 15:39:13 [INFO] We will recompute all MWCS based on the new REF for SX.YDH:SX.YSG
> 2019-06-10 15:39:13 [DEBUG] Ready to compute
> 2019-06-10 15:39:14 [INFO] There are MWCS jobs for some days to recompute for CQ.CHK:CQ.CHS
> 2019-06-10 15:39:14 [DEBUG] Processing MWCS for: CQ_CHK_CQ_CHS.ZZ.08 - 2016-01-01 - 05 days
> /data1/dell/Software/anaconda3/envs/msnoise/lib/python3.6/site-packages/scipy/optimize/minpack.py:779: OptimizeWarning: Covariance of the parameters could not be estimated
>    category=OptimizeWarning)
> /data1/dell/Software/anaconda3/envs/msnoise/lib/python3.6/site-packages/msnoise/move2obspy.py:349: RuntimeWarning: invalid value encountered in double_scalars
>    e = np.sum((phi - m * v) ** 2) / (np.size(v) - 1)
> /data1/dell/Software/anaconda3/envs/msnoise/lib/python3.6/site-packages/msnoise/move2obspy.py:349: RuntimeWarning: divide by zero encountered in double_scalars
>    e = np.sum((phi - m * v) ** 2) / (np.size(v) - 1)
> 2019-06-10 15:39:16 [DEBUG] Processing MWCS for: CQ_CHK_CQ_CHS.ZZ.08 - 2016-01-01 - 15 days
> 2019-06-10 15:39:16 [DEBUG] Processing MWCS for: CQ_CHK_CQ_CHS.ZZ.08 - 2016-01-01 - 30 days
> 2019-06-10 15:39:17 [DEBUG] Processing MWCS for: CQ_CHK_CQ_CHS.ZZ.08 - 2016-01-01 - 45 days
>
> With these warnings, the code still worked, and then new warnings occured:
>
> 2019-06-10 15:42:33 [DEBUG] Processing MWCS for: CQ_CHK_CQ_CHS.ZZ.08 - 2016-04-25 - 15 days
> 2019-06-10 15:42:33 [DEBUG] Processing MWCS for: CQ_CHK_CQ_CHS.ZZ.08 - 2016-04-25 - 30 days
> /data1/dell/Software/anaconda3/envs/msnoise/lib/python3.6/site-packages/msnoise/move2obspy.py:327: RuntimeWarning: divide by zero encountered in true_divide
>    w = 1.0 / (1.0 / (coh[index_range] ** 2) - 1.0)
> /data1/dell/Software/anaconda3/envs/msnoise/lib/python3.6/site-packages/msnoise/move2obspy.py:327: RuntimeWarning: invalid value encountered in true_divide
>    w = 1.0 / (1.0 / (coh[index_range] ** 2) - 1.0)
> 2019-06-10 15:42:34 [DEBUG] Processing MWCS for: CQ_CHK_CQ_CHS.ZZ.08 - 2016-04-25 - 45 days
> 2019-06-10 15:42:34 [DEBUG] Processing MWCS for: CQ_CHK_CQ_CHS.ZZ.08 - 2016-04-25 - 60 days
>
> And the code was still running, the output was as follows:
>
> -3.300000000000000000e+02 -7.905952284653005258e+00 nan 9.937050757172795068e-01
> -3.250000000000000000e+02 -8.138657988098142226e+00 nan 9.976592601521057091e-01
> -3.200000000000000000e+02 -7.632117540664275168e+00 nan 9.894723420074960085e-01
> -3.150000000000000000e+02 -4.285274200620775353e+00 nan 8.972202676215104145e-01
> -3.100000000000000000e+02 -1.489405237146531613e+00 nan 9.666564114380370443e-01
> -3.050000000000000000e+02 -1.821892139101924268e+00 nan 9.604215787607996990e-01
> -3.000000000000000000e+02 8.124339913942842273e+00 nan 9.809556426874297408e-01
> -2.950000000000000000e+02 1.436912500259247771e+00 nan 9.710551233524897796e-01
> -2.900000000000000000e+02 6.558446922684406655e-01 nan 9.932296151167433473e-01
> -2.850000000000000000e+02 5.593336746118490188e-01 inf 9.946785904215598162e-01
> -2.800000000000000000e+02 1.623208176830735067e-02 inf 9.997050009388863678e-01
> -2.750000000000000000e+02 -8.591600654862702413e+00 inf 9.924982060141236406e-01
> -2.700000000000000000e+02 5.866839658532229329e+00 nan 9.421224113061066330e-01
> -2.650000000000000000e+02 4.197498248311506863e+00 nan 9.094520333511715293e-01
> -2.600000000000000000e+02 6.883655670067920340e+00 nan 9.634636129273534966e-01
> -2.550000000000000000e+02 6.902610666221385749e+00 nan 9.689572985482811651e-01
> -2.500000000000000000e+02 -2.224539350609665789e+00 nan 9.149593016370045628e-01
> -2.450000000000000000e+02 -2.569921624096909163e+00 nan 9.193301735266846686e-01
> -2.400000000000000000e+02 -1.787990836057278932e+00 inf 9.549244033726577152e-01
> -2.350000000000000000e+02 -7.489977845697655123e-01 nan 9.923064807292454326e-01
>
> I don’t know why this problem only occured at filter 8 and 9.
> looking forward to know what’s going on
>
> very very thanks!
> weiyu
>
>
> _______________________________________________
> MSNoise mailing list
> MSNoise at mailman-as.oma.be
> http://mailman-as.oma.be/mailman/listinfo/msnoise
>