vpdistance.m

来自「使用支持向量机（svm）方法进行图像的单视度量(single view metr」· M 代码 · 共 12 行

function [d,alpha]=vpdistance(VP,Lines)
if(VP(3)==0)
    a1=atan2(Lines(:,2),Lines(:,1));
    a2=atan2(VP(2),VP(1))+pi/2;
    alpha=abs(a1-repmat(a2,size(a1,1),1));
    alpha=abs((alpha<pi/2).*alpha+(alpha>=pi/2).*(pi-alpha));
    alpha=abs((alpha<pi/2).*alpha+(alpha>=pi/2).*(pi-alpha));
    d=0;
else
    d=abs(((VP(1)*Lines(:,1)+VP(2)*Lines(:,2))/VP(3)+Lines(:,3))./sqrt(Lines(:,1).^2+Lines(:,2).^2));
    alpha=0;
end