Does anyone know how to use the matlab to calculate the minimu distance between a point outside oval and the oval surface？

22 次查看（过去 30 天）

huazai2020 2020-6-28

1
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/555925-does-anyone-know-how-to-use-the-matlab-to-calculate-the-minimu-distance-between-a-point-outside-oval

编辑： huazai2020 2020-7-18

Does anyone know how to use the matlab to calculate the minimu distance between a point outside oval and the oval surface？

3 个评论
显示 1更早的评论隐藏 1更早的评论

Image Analyst 2020-6-28

编辑：Image Analyst 2020-6-28

bwdist() perhaps?

huazai2020 2020-6-29

The oual can be both given in the form of an equation or data points. How to use the bwdist(), could you share me the code?

请先登录，再进行评论。

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

采纳的回答

Matt J 2020-6-30

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/555925-does-anyone-know-how-to-use-the-matlab-to-calculate-the-minimu-distance-between-a-point-outside-oval#answer_458782

编辑：Matt J 2020-6-30

在 MATLAB Online 中打开

You can use trustregprob from the File Exchange

https://www.mathworks.com/matlabcentral/fileexchange/53191-quadratic-minimization-with-norm-constraint?focused=6974391&tab=function

For example, consider the following ellipse and external point y,

A=[2 1;1,2]; %Ellipse equation matrix
y=[1;0.5]; %External point
z=[0.584808315593597 ; 0.201052451754066]; %Closest point
ellipsefun=@(p,q) sum([p;q].*(A*[p;q]))-1 ;
hold on
 fimplicit(ellipsefun, [-1 1.3  -1 1.3])
 plot(y(1),y(2),'rx',z(1),z(2),'bo');
 axis equal
 
hold off

I found the closest point z using the code below,

Aroot=chol(A);
L=Aroot\eye(2);
z=Aroot\trustregprob(L.'*L, L.'*y,1);

and the minimum distance is just,

>> distance=norm(z-y)
distance =
    0.5116

20 个评论
显示 18更早的评论隐藏 18更早的评论

huazai2020 2020-7-7

What is FEX AND H1? Are they these?

Find the projection of point P in R^n on the ellipsoid

E = { x = x0 + U*(z.*radii) : |z| = 1 }, where U is orthogonal matrix of the orientation of E, radii are the axis lengths, and x0 is the center.

Or on generalized conic E = { x : x'*A*x + b'*x + c = 0 }.

The projection is the minimization problem:

min | x - P | (or max | x - P|) for x in E.

Method: solve the Euler Lagrange equation with respect to the Lagrange multiplier, which can be written as polynomial equation (from an idea by Roger Stafford)

Bruno Luong 2020-7-8

FEX short word for File Exchange. Click on it I put the link on my message.

H1: https://www.mathworks.com/help/matlab/matlab_prog/add-help-for-your-program.html

请先登录，再进行评论。

更多回答（2 个）

Image Analyst 2020-7-4

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/555925-does-anyone-know-how-to-use-the-matlab-to-calculate-the-minimu-distance-between-a-point-outside-oval#answer_460961

在 MATLAB Online 中打开

"The oual can be both given in the form of an equation or data points." <== if you have data points (xb, yb) on the boundary of an ellipse/oval, you can use sqrt() to find the distances to some other point (xp, yp). Then use in() to find the minimum distance.

distances = sqrt(xb-xp).^2 + (yb-yp).^2);
minDistance = min(distances)

If you have an image then you need to get the boundary points first:

boundaries = bwboundaries(binaryImageOfEllipse);
boundaries = boundaries{1}; % Extract double array from cell array.
xb = boundaries(:, 2);
yb = boundaries(:, 1);

10 个评论
显示 8更早的评论隐藏 8更早的评论

Image Analyst 2020-7-18

在 MATLAB Online 中打开

No, you shouldn't use an inner loop because your minDistance will just be the distance between (x1,y1) and xwoint(Num2) and ywoint(Num2) instead of the whole array. Do it like this:

clc;
clear all;
A = xlsread('data1.xlsx');
xpoint=A(:,1);
ypoint=A(:,2);
Num1=length(xpoint);
xwoint=A(:,3);
ywoint=A(:,4);
Num2=length(xwoint);
for k=1:1:Num1
	x1 = xpoint(k);
	y1 = ypoint(k);
	distances = sqrt((xwoint-x1).^2 + (ywoint-y1).^2);
	minDistance(k) = min(distances)*1000;
end
plot(minDistance, 'b-');
% xlswrite('thick.xlsx',minDistance,'sheet1');