Remove lines outside plotted circle

Question

Xen 2017-3-20

1
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/330817-remove-lines-outside-plotted-circle

回答： Riya 2025-6-3

Hi people.

I am plotting a circle using the rectangle function and then I am plotting several lines in the same plot using the line function; the lines are within the circle but also extend out of it (of course). How can I keep only the part of the lines within the circle? Or maybe plot the lines initially with a method that somehow confines them within the circle.

I am sure this is easily done as an image analysis problem, but I was wondering if it can be done with plots.

Thanks.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Riya 2025-6-3

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/330817-remove-lines-outside-plotted-circle#answer_1565801

在 MATLAB Online 中打开

Hi Xen,

I understand that you are plotting a circle using the “rectangle” function and adding lines that go beyond the circle. You want to show only the parts of the lines that are inside the circle.

This can be achieved using basic geometry and parametric line equations. The approach is to detect intersection points of the lines with the circle and then redraw only the segments that lie inside.

Kindly refer to the following sample code for the same:

% Define the circle parameters 
center = [0, 0]; 
radius = 5; 
% Draw the circle 
theta = linspace(0, 2*pi, 500); 
x_circle = center(1) + radius * cos(theta); 
y_circle = center(2) + radius * sin(theta); 
plot(x_circle, y_circle, 'k', 'LineWidth', 1.5); 
axis equal; hold on; 
% Define line segments [x1 y1 x2 y2] 
lines = [ 
    -6, -1, 6, 1; 
    -4, 4, 4, -4; 
    -10, 0, 10, 0 
    ]; 
% Loop through each line 
for i = 1:size(lines, 1) 
    x1 = lines(i, 1); y1 = lines(i, 2); 
    x2 = lines(i, 3); y2 = lines(i, 4); 
    
    dx = x2 - x1; 
    dy = y2 - y1; 
    
    a = dx^2 + dy^2; 
    b = 2 * (dx * (x1 - center(1)) + dy * (y1 - center(2))); 
    c = (x1 - center(1))^2 + (y1 - center(2))^2 - radius^2; 
    
    discriminant = b^2 - 4 * a * c; 
    
    if discriminant < 0 
        % Check if both points lie inside the circle 
        if norm([x1 y1] - center) < radius && norm([x2 y2] - center) < radius 
            plot([x1 x2], [y1 y2], 'r'); 
        end 
    else 
        % Compute the intersection points 
        t1 = (-b + sqrt(discriminant)) / (2 * a); 
        t2 = (-b - sqrt(discriminant)) / (2 * a); 
        
        t_min = max(0, min(t1, t2)); 
        t_max = min(1, max(t1, t2)); 
        
        if t_max >= t_min 
            xi1 = x1 + dx * t_min; 
            yi1 = y1 + dy * t_min; 
            xi2 = x1 + dx * t_max; 
            yi2 = y1 + dy * t_max; 
            
            plot([xi1 xi2], [yi1 yi2], 'r', 'LineWidth', 2); 
        end 
    end 
end 
title('Clipped Lines Inside Circle'); 

This code reshapes the plotting so that only the segment of the line within the circle is shown. This works by calculating the intersection points using the quadratic formula and then restricting the line to the portion that lies inside the circle.

Here is the output of the above sample code for your reference: