Comment trouver les 4 points les plus éloignés de ce rectangle

Question

J'ai tracé une liste de coins que j'ai obtenus en utilisant la détection Harris Corner.

Maintenant, je dois trouver le 4 favorise les points qui représenteront les coins du rectangle

Je sais que je peux obtenir les deux coins en diagonale supérieure et inférieure en utilisant

max(C); 
    min(C); 

Où C est une matrice n lignes avec une colonne pour x et y comme

x y 
    0 1 
    2 3 
    4 5 
    6 6 

Mais comment puis-je obtenir les deux autres coins?

Je pensais que je pouvais tourner la matrice et l'utilisation min et max à nouveau, mais bien sûr, que me renvoie une énorme matrice colonne de n (et je veux une matrice colonne 2)

Je me sens comme la réponse est évidente , mais je suis découpage :(

Answer 1

Je ne sais pas comment est limitée la méthode suivante, mais cela a fonctionné pour un exemple similaire à la vôtre:

% detect possible corners 
points = detectHarrisFeatures(BW); 
C = points.Location; 
% compute par-wise distances between all points 
D = pdist2(C,C); 
p = zeros(1,4); 
% compute maximum distance to find first pair 
[m,idx1] = max(D,[],2); 
[~,idx2] = max(m); 
idx1 = idx1(idx2); 
p(1:2) = [idx1, idx2]; 
% add first pair distance to distance matrix so the next pair will be 
% distant from this pair as well, and compute max distance again 
D = bsxfun(@plus,D,sum(D([idx1 idx2],:),1)); 
[m,idx1] = max(D,[],2); 
[~,idx2] = max(m); 
idx1 = idx1(idx2); 
p(3:4) = [idx1, idx2]; 
% plot 
imshow(BW); 
hold on; 
plot(C(:,1),C(:,2),'g.'); 
plot(C(p,1),C(p,2),'rx','LineWidth',2); 

une autre option est d'utiliser les fonctions FEX comme Polygon simplification et Decimate Polygon, et de définir le nombre de vertex désiré à 4.

Answer 2

Une autre solution de rechange, si vous ne voulez pas utiliser la détection d'angle:

im=zeros(100); 

im(40:70,30:80)=1; 
im=imrotate(im,rand*100); 

[x,y]=find(im); 
x=x+2*randn(size(x)); 
y=y+2*randn(size(y)); 
X=[x(:),y(:)]; 

d=ceil(pdist2(X,X)*10)/10; 

[a,b]=find(d==max(d(:))); 

xm=x(a); 
ym=y(a); 

figure,plot(x,y,'ks') 
hold on, plot(xm,ym,'ro','MarkerSize',12,'MArkerFaceColor','r') 
axis image