How to draw a perfect circle and distance calculation

1 回表示 (過去 30 日間)
amj
amj 2018 年 3 月 9 日
コメント済み: amj 2018 年 3 月 11 日
Hi,
I got 2 issues:
  1. to draw a perfect circle
  2. to calculate a perfect distance
This is my code:
load fisheriris
rng(1); % For reproducibility
n = size(meas,1);
idx = randsample(n,1)
X = meas(~ismember(1:n,idx),3:4); % Training data
Y = meas(idx,3:4)
MdlKDT = KDTreeSearcher(X) % not used by mahanalobis
MdlES = ExhaustiveSearcher(X) % Prepare an exhaustive nearest neighbors searcher.
r = 0.3; % Search radius
IdxKDT = rangesearch(MdlKDT,Y,r)
IdxES = rangesearch(MdlES,Y,r)
[IdxKDT IdxES]
cellfun(@isequal,IdxKDT,IdxES)
figure;
plot(X(:,1),X(:,2),'.k');
hold on;
plot(Y(:,1),Y(:,2),'*r');
Mdl = KDTreeSearcher(X)
[n,d] = knnsearch(Mdl,Y,'k',10);
ctr = Y - d(end);
diameter = 2*d(end);
% Draw a circle around the 10 nearest neighbors.
h = rectangle('position',[ctr,diameter,diameter],...
'curvature',[1 1]);
h.LineStyle = ':';
If i used the fisheriris data, it is works well. But if i'm using Faithful Geyser data (the data as in the attachment), the circle becomes weird. The distance calculation as shown in second picture with yellow color also effected where the data points do not lie within a circle.. This is the code to apply the faithful data:
faithfuldat = xlsread('faithful.csv');
fData = faithfuldat(:,:);
[n,dim]=size(fData);
idx = randsample(n,1)
X = fData(~ismember(1:n,idx),:); % Training data
Y = fData(idx,:)
I really appreciate if anyone could advice me how to improve the circle and distance.

採用された回答

Image Analyst
Image Analyst 2018 年 3 月 9 日
編集済み: Image Analyst 2018 年 3 月 9 日
To draw a circle, you can use rectangle(), viscircles(), or the FAQ: http://matlab.wikia.com/wiki/FAQ#How_do_I_create_a_circle.3F
Add
axis equal;
to your code after you plot to make your circles look circular on the screen.
  11 件のコメント
amj
amj 2018 年 3 月 10 日
Dear image analsyt,if so, i would try tomorrow. i have been trying since his morning. I'm still could not fix it
amj
amj 2018 年 3 月 11 日
Dear Image Analyst,
I modified the code:
xx = (range(diff(X(:,1))))
yy = (range(diff(X(:,2))))
ratio = (xx / yy) /sqrt(2)
for j = 1:length(Y)
c = Y(j,:)
Yx = c(1);
Yy = c(2);
pos = [c-r (2*r) (2*r)/ratio];
rectangle('Position',pos,'Curvature',[1 1])
end
axis([1 5.5 45 95])
Unluckily, the central points are not lie in center part. If you could advice an improvement, i would be grateful

サインインしてコメントする。

その他の回答 (0 件)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by