function [pts,status] = WP_int_line_sphere(line,sphere)
%

pts = [];
status = 0;
a = sphere(1);
b = sphere(2);
c = sphere(3);
radius = sphere(4);

[d,cp] = cv_dist_pt_line([a,b,c],line);
if d>radius
    return
end

dir = line(2,:) - line(1,:);
e = norm([a,b,c] - cp);
f = sqrt(radius^2-e^2);
pts = [cp-f*dir;cp+f*dir];
if norm(pts(1,:)-pts(2,:))>eps
    status = 1;
else
    status = 2;
end
