计算两条线之间的角度而不需要计算斜率？ （爪哇）

Question

我有两条线：L1 和 L2。我想计算两条线之间的角度。 L1 有点：{(x1, y1), (x2, y2)}，L2 有点：{(x3, y3), (x4, y4)}。

如何计算这两条线之间形成的角度，而无需计算斜率？我当前遇到的问题是，有时我有水平线（沿 x 轴的线），并且以下公式失败（除以零例外）：

arctan((m1 - m2) / (1 - (m1 * m2)))

其中 m1 和 m2 分别是线 1 和线 2 的斜率。是否有一个公式/算法可以计算两条线之间的角度而不会出现被零除的异常？任何帮助将不胜感激。

这是我的代码片段：

// Calculates the angle formed between two lines
public static double angleBetween2Lines(Line2D line1, Line2D line2)
{
    double slope1 = line1.getY1() - line1.getY2() / line1.getX1() - line1.getX2();
    double slope2 = line2.getY1() - line2.getY2() / line2.getX1() - line2.getX2();
    double angle = Math.atan((slope1 - slope2) / (1 - (slope1 * slope2)));
    return angle;
}

谢谢。

Answer 1

atan2 函数减轻了处理atan 的痛苦。

它被声明为 double atan2(double y, double x) 并将直角坐标 (x,y) 从极坐标转换为角度 theta坐标 (r,theta)

所以我将你的代码重写为

public static double angleBetween2Lines(Line2D line1, Line2D line2)
{
    double angle1 = Math.atan2(line1.getY1() - line1.getY2(),
                               line1.getX1() - line1.getX2());
    double angle2 = Math.atan2(line2.getY1() - line2.getY2(),
                               line2.getX1() - line2.getX2());
    return angle1-angle2;
}

Answer 2

在这种情况下，点积可能更有用。 在这里您可以找到一个 Java 几何包，它提供了一些有用的帮助器。下面是他们用于确定两个 3-d 点之间角度的计算。希望它能让您开始：

public static double computeAngle (double[] p0, double[] p1, double[] p2)
{
  double[] v0 = Geometry.createVector (p0, p1);
  double[] v1 = Geometry.createVector (p0, p2);

  double dotProduct = Geometry.computeDotProduct (v0, v1);

  double length1 = Geometry.length (v0);
  double length2 = Geometry.length (v1);

  double denominator = length1 * length2;

  double product = denominator != 0.0 ? dotProduct / denominator : 0.0;

  double angle = Math.acos (product);

  return angle;
}

祝您好运！

Answer 3

dx1 = x2-x1;
dy1 = y2-y1;
dx2 = x4-x3;
dy2 = y4-y3;

d = dx1*dx2 + dy1*dy2;   // dot product of the 2 vectors
l2 = (dx1*dx1+dy1*dy1)*(dx2*dx2+dy2*dy2) // product of the squared lengths

angle = acos(d/sqrt(l2));

两个向量的点积等于角度的余弦乘以两个向量的长度。这将计算点积，除以向量的长度，并使用反余弦函数来恢复角度。

Answer 4

也许我的 Android 坐标系统方法对某人有用（使用 Android PointF 类来存储点）

/**
 * Calculate angle between two lines with two given points
 *
 * @param A1 First point first line
 * @param A2 Second point first line
 * @param B1 First point second line
 * @param B2 Second point second line
 * @return Angle between two lines in degrees
 */

public static float angleBetween2Lines(PointF A1, PointF A2, PointF B1, PointF B2) {
    float angle1 = (float) Math.atan2(A2.y - A1.y, A1.x - A2.x);
    float angle2 = (float) Math.atan2(B2.y - B1.y, B1.x - B2.x);
    float calculatedAngle = (float) Math.toDegrees(angle1 - angle2);
    if (calculatedAngle < 0) calculatedAngle += 360;
    return calculatedAngle;
}

它为任何象限返回以度为单位的正值：0 <= x < 360

你可以在这里查看我的实用类

Answer 5

获取角度的公式为 tan a = (slope1-slope2)/(1+slope1*slope2)

您正在使用：

tan a = (slope1 - slope2) / (1 - slope1 * slope2)

所以它应该是：

double angle = Math.atan((slope1 - slope2) / (1 + slope1 * slope2));

Answer 6

首先，你确定括号的顺序正确吗？我认为（可能是错误的）应该是这样的：

   double slope1 = (line1.getY1() - line1.getY2()) / (line1.getX1() - line1.getX2());
   double slope2 = (line2.getY1() - line2.getY2()) / (line2.getX1() - line2.getX2());

其次，您可以对除零的 div 做两件事：您可以捕获异常并处理它

double angle;
try
{
    angle = Math.atan((slope1 - slope2) / (1 - (slope1 * slope2)));
catch (DivideByZeroException dbze)
{
    //Do something about it!
}

……或者您可以检查您的除数是否不为零 在尝试操作之前。

if ((1 - (slope1 * slope2))==0)
{
    return /*something meaningful to avoid the div by zero*/
}
else 
{
    double angle = Math.atan((slope1 - slope2) / (1 - (slope1 * slope2)));
    return angle;
 }

Answer 7

检查这个Python代码：

import math
def angle(x1,y1,x2,y2,x3,y3):

  if (x1==x2==x3 or y1==y2==y3):
    return 180
  else:
    dx1 = x2-x1
    dy1 = y2-y1
    dx2 = x3-x2
    dy2 = y3-y2
    if x1==x2:
      a1=90
    else:
      m1=dy1/dx1
      a1=math.degrees(math.atan(m1))
    if x2==x3:
      a2=90
    else:
      m2=dy2/dx2
      a2=math.degrees(math.atan(m2))
    angle = abs(a2-a1)
    return angle

print angle(0,4,0,0,9,-6)

Answer 8

dx1=x2-x1 ; dy1=y2-y1 ; dx2=x4-x3 ;dy2=y4-y3.

Angle(L1,L2)=pi()/2*((1+sign(dx1))* (1-sign(dy1^2))-(1+sign(dx2))*(1-sign(dy2^2)))
           +pi()/4*((2+sign(dx1))*sign(dy1)-(2+sign(dx2))*sign(dy2))
           +sign(dx1*dy1)*atan((abs(dx1)-abs(dy1))/(abs(dx1)+abs(dy1)))
           -sign(dx2*dy2)*atan((abs(dx2)-abs(dy2))/(abs(dx2)+abs(dy2)))