Once the sq.pp.filter_cells() PR is merged, we should be able to resue basically all code to just discard points that don't fall into a seg mask.