在MATLAB中找到两个未知方程曲线之间的距离

Question

我目前正在进行MATLAB项目，我需要找到两块膝盖骨头之间的距离。

使用chenvese程序进行活动轮廓并设法绘制骨骼的曲线后，我需要找到它们之间的空间。

我使用interp1函数尝试了各种程序，但我的曲线没有相同的大小/长度。所以我试着用findobj提取曲线方程，但它没有让我到处都是。

%*********************** Chen Vese ***********************************
m=zeros(size(IBinaire,1),size(IBinaire,2));
m(300:900,400:1200)=1; %200,1200
seg = chenvese(IBinaire,m,500,0.2,'chan');

%*********************** Courbes de contours *************************
c = contour(seg);
s = getcontourlines(c);
plot(s(1).x,-s(1).y, 'b', s(2).x,-s(2).y,'g', s(3).x,-s(3).y,'r')

h = findobj('type', 'line', 'marker', '-and', 'b', [1 0]);
xx1 = get(h, 'XData')
yy1 = get(h, 'YData')


h = findobj('type', 'line', 'marker', '-and', 'g', [1 0]);
xx2 = get(h, 'XData')
yy2 = get(h, 'YData')

%****************Distance entre lignes ************************
z1=xx1+1i*yy1;
z2=xx2+1i*yy2;
i=interligne(z1,z2); 
plot(i)

interligne是一个计算两条曲线之间距离的程序，但它不起作用，告诉我曲线大小不一样......

Answer 1

以下是一些示例代码，如何使用interp1来解决它：

% just make some test data
s = [
    struct('x',linspace(0,10,20)','y',1+(3-linspace(0,10,20)').^2);
    struct('x',linspace(0,10,30)','y',-(4-linspace(0,10,30)').^2);
    ];

% the logic is here
f = @(i,p) interp1(1:numel(s(i).x), [s(i).x,s(i).y],p);

d = @(X) norm(f(1,X(1))-f(2,X(2)));

X = fminsearch(d,[5,5]) % use a reasonable start point

% just some visualization
res = [f(1,X(1));f(2,X(2))]
clf
line(s(1).x,s(1).y,'Color','b')
line(s(2).x,s(2).y,'Color','g')
line(res(:,1),res(:,2),'Color','r','Marker','x')
axis([0,10,-5,5])

Answer 2

在这种情况下，您可能正在寻找线上任意两点之间的最小距离：

%example data line A
A=rand(6,1)+i*rand(6,1);
%example data line B
B=rand(4,1)+i*rand(4,1);
AA=repmat(A,size(B.'));
BB=repmat(B.',size(A));
pwdist=sqrt(real(AA(:)-BB(:)).^2+imag(AA(:)-BB(:)).^2)
[d,f]=min(pwdist);
[a,b]=ind2sub(size(AA),f);

使用pdist2也是如此，但需要统计工具箱。如果可用，请使用它。它可能更快。