如何计算三维空间中两点相对于相机方向的夹角

import numpy as np

def calculate_angle(camera_pos, camera_dir, point1, point2):
    # Step 3: Calculate relative direction vectors
    dir_vector1 = point1 - camera_pos
    dir_vector2 = point2 - camera_pos

    # Step 4: Normalize direction vectors
    dir_vector1 /= np.linalg.norm(dir_vector1)
    dir_vector2 /= np.linalg.norm(dir_vector2)

    # Step 5: Calculate dot product
    dot_product = np.dot(dir_vector1, dir_vector2)

    # Step 6: Convert to angle in degrees
    angle = np.arccos(dot_product) * 180 / np.pi

    return angle

# Example usage
camera_pos = np.array([0, 0, 0])
camera_dir = np.array([0, 0, 1])
point1 = np.array([1, 1, 1])
point2 = np.array([-1, -1, -1])

angle = calculate_angle(camera_pos, camera_dir, point1, point2)
print("Angle between point1 and point2 relative to camera direction:", angle)