[MSNoise] plots

Thomas Lecocq thomas.lecocq at seismology.be
Tue Nov 11 12:11:48 UTC 2014


Josiah,

A possible reason for "missing" data in the final output would be that 
either the dt/t computation did not manage to select enough reliable 
points for computing the weighted linear regression, an thus dt/t is not 
calculated for this/these dates. You can check in the DTT folder if all 
daily files are present.

Thomas

Le 31/10/2014 00:35, Josiah Ensing a écrit :
> Thank you. When I used MSNoise 1.2.5, the second issue I wrote about
> (mismatching x axis) is not a problem anymore. The first issue, is still a
> problem. But this time the 10days , moving window plot is plotted, and the
> 1 and 3 days moving windows a re plotted, but not for the full time scale,
> and the 5 days moving window is blank.
>
> 2014-10-31 4:03 GMT+13:00 Thomas Lecocq <thomas.lecocq at oma.be>:
>
>> Hi Josiah,
>>
>> your email was blocked because it contained no subject...
>>
>> Could you test using 1.2.5 ? you just have to replace the python files in
>> the directory and all should work exactly the same...
>>
>> Thomas
>>
>> Le 24/10/2014 00:49, Josiah Ensing a écrit :
>>
>>> Using MSNoise 1.24 there are two python scripts that plot; they are
>>> 07plot_dtt.py and s07plot_dtt.py. I get errors for either.
>>>
>>> Using s07plot_dtt.py yields attachment s07plot_dtt.png. There are two
>>> problems:
>>>
>>> 1)  I entered mov_stack values 1,3,5, and 10 days into the configurator
>>> but
>>> the graphical output includes only 1, 3, and 5 days moving window graphs,
>>> and the only graph with any data plotted onto it is the 5 days moving
>>> window.
>>>
>>> 2) the graph x axes automatically include a much larger time period than
>>> covered by the data plotted. This includes dates back to 2013 and forward
>>> into 2016, while all the data I have is from 2014. I am unsure if this
>>> problem is related to some of my settings in the configurator. My settings
>>> that I can see might be related are: enddate=2100-01-01, ref_begin=-100,
>>> ref_end=0, startdate=2014-05-04.
>>>
>>> Using 07plot_dtt.py yields attachment image 07.plot_dtt.png  the
>>> following:
>>>
>>>   runfile('/home/devora/MSNoise-master/07.plot_dtt.py',
>>>>> wdir=r'/home/devora/MSNoise-master')
>>> UMD has deleted: database_tools, msnoise_table_def
>>>
>>> loading 1 days
>>>
>>> loading 3 days
>>>
>>> Traceback (most recent call last):
>>>
>>> File "<stdin>", line 1, in <module>
>>>
>>> File
>>> "/usr/lib/python2.7/site-packages/spyderlib/widgets/
>>> externalshell/sitecustomize.py",
>>> line 540, in runfile
>>>
>>> execfile(filename, namespace)
>>>
>>> File "/home/devora/MSNoise-master/07.plot_dtt.py", line 128, in <module>
>>>
>>> plt.fill_between(ALL.index,ALL[dttname]-ALL[errname],ALL[
>>> dttname]+ALL[errname],lw=1,color='red',zorder=-1,alpha=0.3)
>>>
>>> File "/usr/lib64/python2.7/site-packages/matplotlib/pyplot.py", line
>>> 2757,
>>> in fill_between
>>>
>>> interpolate=interpolate, **kwargs)
>>>
>>> File "/usr/lib64/python2.7/site-packages/matplotlib/axes.py", line 6988,
>>> in
>>> fill_between
>>>
>>> x = ma.masked_invalid(self.convert_xunits(x))
>>>
>>> File "/usr/lib64/python2.7/site-packages/numpy/ma/core.py", line 2239, in
>>> masked_invalid
>>>
>>> condition = ~(np.isfinite(a))
>>>
>>> TypeError: ufunc 'isfinite' not supported for the input types, and the
>>> inputs could not be safely coerced to any supported types according to the
>>> casting rule ''safe''
>>>
>>>
>>> Has anyone else encountered these problems before? Does anyone have any
>>> suggestions for what I might do to overcome them?
>>>
>>>
>>> Thankyou
>>>
>>>
>>> Josiah
>>>
>>>
>>> _______________________________________________
>>> MSNoise mailing list
>>> MSNoise at mailman-as.oma.be
>>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>>
>> _______________________________________________
>> MSNoise mailing list
>> MSNoise at mailman-as.oma.be
>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>
>
>
>
>
> _______________________________________________
> MSNoise mailing list
> MSNoise at mailman-as.oma.be
> http://mailman-as.oma.be/mailman/listinfo/msnoise



More information about the MSNoise mailing list