missing muselet sources in segmentation map
Occasionally SExtractor can create a catalog with sources that do not have a corresponding ID in the segmentation map.
This can produce the following error
[DEBUG] writing sources to: /media/Elements/OST/muselet/lines
Created:1257/2795Traceback (most recent call last):
File "<ipython-input-31-41668c07050a>", line 1, in <module>
cont,sing,raw=muselet('ost_subcube_randomNoise.fits', step=3)
File "/home/ppxag5/anaconda3/envs/py/lib/python3.7/site-packages/mpdaf/sdetect/muselet.py",
line 1596, in muselet
n_cpu=n_cpu)
File "/home/ppxag5/anaconda3/envs/py/lib/python3.7/site-packages/mpdaf/sdetect/muselet.py",
line 1491, in step3
write_line_sources(cat_lines, dir_, cube, ima_size, n_cpu=n_cpu)
File "/home/ppxag5/anaconda3/envs/py/lib/python3.7/site-packages/mpdaf/sdetect/muselet.py",
line 1355, in write_line_sources
write_line_source_single(row, dir_, cube, ima_size)
File "/home/ppxag5/anaconda3/envs/py/lib/python3.7/site-packages/mpdaf/sdetect/muselet.py",
line 1191, in write_line_source_single
size = get_mask_minsize(im_seg, [dec, ra]) + 2. #pad by 1 arcsec border
File "/home/ppxag5/anaconda3/envs/py/lib/python3.7/site-packages/mpdaf/sdetect/muselet.py",
line 1134, in get_mask_minsize
max_dist = np.max(np.abs(obj_coord_pix), axis=0)
File "/home/ppxag5/anaconda3/envs/py/lib/python3.7/site-packages/numpy/core/fromnumeric.py",
line 2505, in amax
initial=initial)
File "/home/ppxag5/anaconda3/envs/py/lib/python3.7/site-packages/numpy/core/fromnumeric.py",
line 86, in _wrapreduction
return ufunc.reduce(obj, axis, dtype, out, **passkwargs)