https://doi.org/10.5201/ipol.2018.194
Tip revision: 3fe71958c866901c02af48041516988f9d2ace04 authored by Software Heritage on 06 December 2017, 00:00:00 UTC
ipol: Deposit 1325 in collection ipol
ipol: Deposit 1325 in collection ipol
Tip revision: 3fe7195
isAdjacent.m
function [r, t, st] = isAdjacent(seg1,seg2,rho,theta)
%% [r, t, st] = isAdjacent(seg1,seg2,rho,theta)
% For input output description follow GG.m script and comments
%
% Written by Boshra Rajaei 3/11/16
% b.rajaei@sadjad.ac.ir
%
% Copyright (c) 2016-2017 boshra rajaei <b.rajaei@sadjad.ac.ir>
%
% This program is a free software: you can redistribute it and/or modify
% it under the terms of the GNU Affero General Public License as
% published by the Free Software Foundation, either version 3 of the
% License, or (at your option) any later version.
%
% This program is distributed in the hope that it will be useful,
% but WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
% GNU Affero General Public License for more details.
%
% You should have received a copy of the GNU Affero General Public License
% along with this program. If not, see <http://www.gnu.org/licenses/>.
%
%%
line1(1,:) = seg1.strtip;
line1(2,:) = seg1.endtip;
line2(1,:) = seg2.strtip;
line2(2,:) = seg2.endtip;
angle_thr = 5;
dist_thr = 5;
r = -1;t=-1;st=-1;
flag = 0;
for i=1:2 %the two tips of each line
for j=1:2
if(i==j) continue; end
st = (i~=1); %i=1,j=2 : st = 0 i=2,j=1 : st = 1 %(i-1)*2+j;%
dist = norm(line1(i,:)-line2(j,:));
if(dist>rho) continue; end
p12 = line1(i,:); p21 = line2(j,:);
p11 = line1; p11(i,:)=[]; p22 = line2; p22(j,:)=[];
p1 = p12 - p11;
p2 = p22 - p21;
ang = radtodeg(acos(p1*p2'/(norm(p1)*norm(p2)))); %line i , j
p3 = p21 - p12;
ang1 = radtodeg(acos(p1*p3'/(norm(p1)*norm(p3)))); %line i and middle line
ang2 = radtodeg(acos(p3*p2'/(norm(p3)*norm(p2)))); %line j and middle line
if(ang>theta) continue; end
if(dist>dist_thr && (ang1>ang+angle_thr || ang2>ang+angle_thr)) continue; end
if((dist<=rho && ang<=theta/3) || (dist<=2*rho/3 && ang<=2*theta/3) || (dist<=rho/3 && ang<=theta))
r = dist;
t = ang;
flag = 1;
end
end
if(flag)
if(r<=rho/3) r = rho/3;
elseif(r<=2*rho/3) r = 2*rho/3;
else r = rho;
end
if(t<=3) t=3;
elseif(t<=theta/3) t = theta/3;
elseif(t<=2*theta/3) t = 2*theta/3;
else t = theta;
end
return;
end
end