以编程方式找到光线和平面的交点？

import numpy as np

def find_intersection(ray_origin, ray_direction, plane_point, plane_normal):
    # 计算光线的参数方程
    t = np.linspace(0, 1, 100)  # 光线的参数范围
    ray_points = ray_origin + t[:, np.newaxis] * ray_direction

    # 计算交点的参数
    denom = np.dot(ray_direction, plane_normal)
    numer = np.dot(plane_point - ray_origin, plane_normal)
    intersection_param = numer / denom

    # 计算交点坐标
    intersection_point = ray_origin + intersection_param * ray_direction

    # 判断交点是否存在
    if intersection_param >= 0 and intersection_param <= 1:
        return intersection_point
    else:
        return None

# 示例用法
ray_origin = np.array([0, 0, 0])
ray_direction = np.array([1, 1, 1])
plane_point = np.array([0, 0, 1])
plane_normal = np.array([0, 0, 1])

intersection = find_intersection(ray_origin, ray_direction, plane_point, plane_normal)
if intersection is not None:
    print("交点坐标：", intersection)
else:
    print("交点不存在")