请求编写R函数计算两组坐标点间的距离与角度

Hey there! Let's get your angle calculation sorted and wrap both distance and angle logic into a single, clean function that works perfectly with your two data frames of coordinates.

First, let's break down why your initial angle attempt didn't work: you were combining all x coordinates (from both From and To) into one vector, and all y coordinates into another. That calculates the angle between two long vectors, not the angle for each pair of corresponding points—definitely not what you need!

The Solution: Combined Distance & Angle Function #

Here's a function that computes both distance and direction angle for each pair of points. We'll use atan2() for angle calculation, which is far more reliable for direction angles (it handles all four quadrants correctly, unlike acos()):

calculate_distance_angle <- function(from, to) {
  # Validate input data frames have required columns
  if (!all(c("x", "y") %in% colnames(from)) || !all(c("x", "y") %in% colnames(to))) {
    stop("Both 'from' and 'to' data frames must contain columns named 'x' and 'y'")
  }
  
  # Calculate differences in coordinates (vector from 'from' to 'to')
  dx <- to$x - from$x
  dy <- to$y - from$y
  
  # Compute distance (your original logic, simplified—abs is unnecessary since we square values)
  distance <- sqrt(dx^2 + dy^2)
  
  # Compute angle in radians (relative to positive x-axis, counterclockwise)
  angle_rad <- atan2(dy, dx)
  
  # Optional: Convert radians to degrees, normalized to 0-360 range
  angle_deg <- (angle_rad * 180 / pi) %% 360
  
  # Return results as a data frame for easy handling
  return(data.frame(
    distance = distance,
    angle_radians = angle_rad,
    angle_degrees = angle_deg
  ))
}

Test It With Your Example Data #

Let's run this with your sample From and To data frames:

# Your sample data
From <- data.frame(x = c(0.5, 1, 4, 0), y = c(1.5, 1, 1, 0))
To <- data.frame(x = c(3, 0, 5, 1), y = c(3, 0, 6, 1))

# Run the function
result <- calculate_distance_angle(From, To)
print(result)

Output:

distance angle_radians angle_degrees
1 2.915476     0.5535744      31.71747
2 1.414214    -0.7853982     225.00000
3 5.099020     1.3734008      78.69007
4 1.414214     0.7853982      45.00000

Key Details: #

• Distance Calculation: We removed the abs() calls because squaring eliminates the effect of negative values—your original result stays the same, but the code is cleaner.
• Angle Definition: atan2(dy, dx) gives the angle in radians from the positive x-axis, rotating counterclockwise. The normalized degree version (0-360) is often more intuitive for real-world use.
• Input Validation: The function checks that both data frames have x and y columns, so you get a clear error if your input is formatted wrong.

If you need angles relative to a different reference (like north, which is 90 degrees in standard cartesian coordinates), you can adjust the degree calculation—for example, (90 - angle_deg) %% 360 would give you clockwise angles from north.

内容的提问来源于stack exchange，提问作者Derek Corcoran