在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

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

以下是一些示例代码，说明如何使用

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])

---

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

%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

也可以实现同样的功能，但需要统计工具箱。如果可用，请使用它。它可能更快。

我认为使用

interp1

你是正确的-但是使用单独的维度。如果你有统计工具箱，你也许可以使用

pdist2

。实际上，我更关注的是相同x值和每个x值的蓝色/绿色曲线之间的距离。我有pdist2，我正在尝试。我正在尝试适应我的问题，但我不知道如何告诉matlab使用我的蓝色/红色曲线。我使用了

findobj

的结果，但它给了我一个空白结果。请尝试

findall

。