drops some data unexpectedly if cover has overlapping sids.
The following workaround produces the desired result.
sids_ = pystare.to_compressed_range(sids)
sids_ = pystare.expand_intervals(sids_,int(max_res),multi_resolution=False)
p = p[p.stare_intersects(sids_)]
This workaround is not ideal as the range code can be fragile.
It is unfortunate that the multi_resolution capability is broken at present, leading to more sids to compare. See bug pystare.expand_intervals issue 122.
p is a STAREDataFrame
drops some data unexpectedly if cover has overlapping sids.
The following workaround produces the desired result.
This workaround is not ideal as the range code can be fragile.
It is unfortunate that the multi_resolution capability is broken at present, leading to more sids to compare. See bug pystare.expand_intervals issue 122.