How can I find all points inside a ellipse

48 views (last 30 days)
azdoud youssef on 29 Dec 2017
Edited: Deepu Kumar on 30 Jul 2021
I have the bounds of an ellipse in x- and y-coordinates. My end goal is to get a matrix of all x-values and a corresponding matrix of all y-values for the points in/out that ellipse. If possible, I would also like to color in the specified ellipse on a 2d plot. Does anyone know which commands/functions I can use to do this? Most of what I have seen simply confirms whether input values lie in/out an ellipse. Thank you!
such that a sample of the ellipse is represented as :
xCenter = 15;
yCenter = 10;
theta = 0 : 0.01 : 2*pi;
x = xRadius * cos(theta) + xCenter;
y = yRadius * sin(theta) + yCenter;
plot(x, y, 'LineWidth', 3);
Walter Roberson on 29 Dec 2017
There are an infinite number of points inside (or on) an ellipse. Even if you confine yourself to discrete x values, there would be an infinite number of y values, corresponding to the infinite density of the chord. To do anything useful, you need to either work with a symbolic formula or else discretize both x and y.

Star Strider on 29 Dec 2017
‘My end goal is to get a matrix of all x-values and a corresponding matrix of all y-values for the points in/out that ellipse.’
There would be an infinite number of them, since you have defined an infinite 2D space. Your computer (or any computer) likely cannot handle it.
‘If possible, I would also like to color in the specified ellipse on a 2d plot’
That’s more tractable. See the patch (link) and related functions (such as area and fill, linked to at the end of that page) to create your plot. Remember to specify axis equal to avoid distorting your ellipse plot.
Star Strider on 30 Dec 2017
If you want to know the points inside and outside the ellipse, and there are a specific number of them, use the inpolygon (link) function. I know that you mentioned that you apparently do not want to use it, however it is the most efficient option.
Otherwise, if you want the area of the ellipse, you can use the trapz function twice (once in each dimension), or create an anonymous function representation for your ellipse and use integral2, or for a vector valued function, integral once in each dimension.

azdoud youssef on 30 Dec 2017
clear ,close all;
xCenter = 15;
yCenter = 10;
theta = 0 : 0.01 : 2*pi;
x = xRadius * cos(theta) + xCenter;
y = yRadius * sin(theta) + yCenter;
xq=rand(1,3000)*100;
yq=rand(1,3000)*100;
in = inpolygon(xq,yq,x,y);
hold on
plot(xq(in), yq(in), 'ro');
plot(xq(~in), yq(~in), 'b*');
plot(x, y, 'LineWidth', 1);
axis square;
xlim([0 20]);
ylim([0 20]);
grid on;
% Draw lines along the axes
hold on;
'LineWidth', 4, 'Color', [1,0,0]);
'LineWidth', 4, 'Color', [1,0,0]);
% Enlarge figure to full screen.
set(gcf, 'Units', 'Normalized', 'OuterPosition', [0 0 1 1]);
Waqar Khan on 18 Mar 2021
@azdoud youssef how did you get the values from ellipse, Could you please share the code or give me suggestions.

Deepu Kumar on 30 Jul 2021
Edited: Deepu Kumar on 30 Jul 2021
I think this is what you were asking.
In place of xCenter and yCenter, you can choose what ever the values that you want.
xCenter = 0;
yCenter = 0;
theta = 0 : (2*pi)/100 : 2*pi;
x(i,j) = xx * cos(theta(j)) + xCenter;
y(i,j) = yy * sin(theta(j)) + yCenter;
end
end
plot(x, y, 'LineWidth', 3);