把线状矢量的点坐标读取为numpy数组，并对相邻的点计算距离，转为米单位，当相邻点之间的距离大于100米时，则在相邻点之间加点

好的，我会根据您提供的线状矢量的点坐标，将其读取为numpy数组，并计算相邻点之间的距离。如果距离大于100米，则在相邻点之间添加新的点。

import numpy as np
from math import radians, sin, cos, sqrt, atan2

# 请提供线状矢量的点坐标（经纬度）
points = [(0, 0), (1, 1), (2, 3), (4, 4)]

# 将点坐标读取为numpy数组
point_array = np.array(points)

# 定义haversine函数来计算两个经纬度点之间的距离
def haversine(lat1, lon1, lat2, lon2):
    R = 6371.0  # 地球平均半径（单位：km）

    # 将经纬度转换为弧度
    lat1_rad = radians(lat1)
    lon1_rad = radians(lon1)
    lat2_rad = radians(lat2)
    lon2_rad = radians(lon2)

    # 计算差值
    dlat = lat2_rad - lat1_rad
    dlon = lon2_rad - lon1_rad

    # 应用haversine公式计算距离
    a = sin(dlat / 2)**2 + cos(lat1_rad) * cos(lat2_rad) * sin(dlon / 2)**2
    c = 2 * atan2(sqrt(a), sqrt(1 - a))

    distance = R * c * 1000  # 转换为米单位
    return distance

# 计算相邻点之间的距离并添加新的点
new_points = [points[0]]
for i in range(len(points) - 1):
    lat1, lon1 = points[i]
    lat2, lon2 = points[i+1]
    distance = haversine(lat1, lon1, lat2, lon2)

    if distance > 100:
        num_points = int(distance / 100)  # 计算需要添加的点的数量
        lat_diff = (lat2 - lat1) / (num_points + 1)
        lon_diff = (lon2 - lon1) / (num_points + 1)

        for j in range(1, num_points + 1):
            new_lat = lat1 + j * lat_diff
            new_lon = lon1 + j * lon_diff
            new_points.append((new_lat, new_lon))

    new_points.append((lat2, lon2))

# 将新的点坐标读取为numpy数组
new_point_array = np.array(new_points)

# 计算相邻点之间的距离
distances = []
for i in range(len(new_points) - 1):
    lat1, lon1 = new_points[i]
    lat2, lon2 = new_points[i+1]
    distance = haversine(lat1, lon1, lat2, lon2)
    distances.append(distance)

print("原始点坐标数组：")
print(point_array)

print("加点后的点坐标数组：")
print(new_point_array)

print("相邻点之间的距离（单位：米）：")
print(distances)

请将您的线状矢量的点坐标替换到 `points` 变量中，然后运行上述代码即可得到相邻点之间的距离，并在距离大于100米时添加新的点。