AggregateAndPlot crash plotting showers #42
Description
When testing a fix for issue 40 i noticed that AggregateAndPlot will crash for the December-2018 data as shown below. Data are successfully being loaded, but it looks like the lower and upper bounds of something are +/- inf. I'll see if i can track it down.
(wmpl311) mmcintyre@mark-mcintyre-vm:~/data/agptest$ python -m wmpl.Trajectory.AggregateAndPlot -o ./plots /srv/meteor/rms/gmn/extracted_data/trajectories/2018/201812/
Writing summary file...
Trajectory summary saved to: ./plots/trajectory_summary.txt
Plotting showers...
Processing shower: 4 GEM
/home/mmcintyre/src/WesternMeteorPyLib/wmpl/Trajectory/AggregateAndPlot.py:1268: RuntimeWarning: divide by zero encountered in log10
ax_mass.hist(np.log10(mass_data_filt), bins=len(mass_data_filt), cumulative=-1,
Traceback (most recent call last):
File "<frozen runpy>", line 198, in _run_module_as_main
File "<frozen runpy>", line 88, in _run_code
File "/home/mmcintyre/src/WesternMeteorPyLib/wmpl/Trajectory/AggregateAndPlot.py", line 2386, in <module>
generateShowerPlots(output_dir, traj_list, min_members=30, P_0m=P_0M, max_radiant_err=0.5)
File "/home/mmcintyre/src/WesternMeteorPyLib/wmpl/Trajectory/AggregateAndPlot.py", line 1268, in generateShowerPlots
ax_mass.hist(np.log10(mass_data_filt), bins=len(mass_data_filt), cumulative=-1,
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/home/mmcintyre/miniconda3/envs/wmpl311/lib/python3.11/site-packages/matplotlib/__init__.py", line 1465, in inner
return func(ax, *map(sanitize_sequence, args), **kwargs)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/home/mmcintyre/miniconda3/envs/wmpl311/lib/python3.11/site-packages/matplotlib/axes/_axes.py", line 6862, in hist
m, bins = np.histogram(x[i], bins, weights=w[i], **hist_kwargs)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/home/mmcintyre/miniconda3/envs/wmpl311/lib/python3.11/site-packages/numpy/lib/histograms.py", line 780, in histogram
bin_edges, uniform_bins = _get_bin_edges(a, bins, range, weights)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/home/mmcintyre/miniconda3/envs/wmpl311/lib/python3.11/site-packages/numpy/lib/histograms.py", line 426, in _get_bin_edges
first_edge, last_edge = _get_outer_edges(a, range)
^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/home/mmcintyre/miniconda3/envs/wmpl311/lib/python3.11/site-packages/numpy/lib/histograms.py", line 315, in _get_outer_edges
raise ValueError(
ValueError: supplied range of [-inf, -inf] is not finite