# calculating angle between line of best fit and x axis

조회 수: 3(최근 30일)
C.G. 2021년 10월 7일
댓글: Star Strider 2021년 10월 7일
I have a set of coordinate points from a picture. I have used polyfit to find the line of best fit of these, but I want to find the angle between the line of best fit and the x axis, so the angle the pile of grains is sitting at. Can somebody tell me how I can do this?
%find the circles and get the coordinates
imshow(I)
[X,Y]=meshgrid(1:size(I));
Rmin = 10;
Rmax = 20;
[cd, rd] = imfindcircles(I,[Rmin Rmax], 'ObjectPolarity', 'dark', 'Sensitivity', 0.94, 'Method','TwoStage');
viscircles(cd, rd, 'Color', 'b');
subplot(1,3,1)
imshow(I)
hold on
scatter(cd(:,1), cd(:,2),'r*')
cd = cd(cd(:,2) <315, :)
subplot(1,3,2)
imshow(I)
hold on
scatter(cd(:,1), cd(:,2),'r*')
%plot line of best fit
p = polyfit(cd(:,1),cd(:,2),1);
Bfit = polyval(p,cd(:,1));
subplot(1,3,3)
imshow(I)
hold on
scatter(cd(:,1),cd(:,2))
hold on
plot(cd(:,1),Bfit,'-r')

댓글을 달려면 로그인하십시오.

### 채택된 답변

Star Strider 2021년 10월 7일
The slope is simply the tangent, so —
slope_angle = atand(p(1)) % Slope Angle (°)
Here, that comes out to be about 14.6°. The slope appears negative in the image plot because of the axes orientation being different in the images. Plotting it with both axes in their normal orientation (added at the end) shows it to be positive.
%find the circles and get the coordinates
imshow(I)
[X,Y]=meshgrid(1:size(I));
Rmin = 10;
Rmax = 20;
[cd, rd] = imfindcircles(I,[Rmin Rmax], 'ObjectPolarity', 'dark', 'Sensitivity', 0.94, 'Method','TwoStage');
viscircles(cd, rd, 'Color', 'b');
subplot(1,3,1)
imshow(I)
hold on
scatter(cd(:,1), cd(:,2),'r*')
cd = cd(cd(:,2) <315, :)
cd = 17×2
124.1955 220.2184 68.5603 219.1154 274.1195 266.7293 194.0057 236.0971 322.3071 281.7616 225.4457 250.7806 455.6494 314.4799 376.8527 296.7009 81.8340 218.4861 255.0560 263.2335
subplot(1,3,2)
imshow(I)
hold on
scatter(cd(:,1), cd(:,2),'r*')
%plot line of best fit
p = polyfit(cd(:,1),cd(:,2),1)
p = 1×2
0.2602 197.6515
Bfit = polyval(p,cd(:,1));
subplot(1,3,3)
imshow(I)
hold on
scatter(cd(:,1),cd(:,2))
hold on
plot(cd(:,1),Bfit,'-r')
slope_angle = atand(p(1)) % Slope Angle (°)
slope_angle = 14.5834
figure
plot(cd(:,1),cd(:,2),'.')
hold on
plot(cd(:,1),Bfit,'-r')
hold off
axis('equal')
title(sprintf('Slope = %.1f°',slope_angle))
.
##### 댓글 수: 2표시숨기기 이전 댓글 수: 1
Star Strider 2021년 10월 7일
As always, my pleasure!
.

댓글을 달려면 로그인하십시오.

### 추가 답변(1개)

Alan Stevens 2021년 10월 7일
편집: Alan Stevens 2021년 10월 7일
Take the arctangent of the gradient of the straight line.
##### 댓글 수: 1표시숨기기 없음
C.G. 2021년 10월 7일
how do i get the gradient from my code above?

댓글을 달려면 로그인하십시오.

### Community Treasure Hunt

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

Start Hunting!

Translated by