我现在有北京市石景山区近1万个点位的坐标和数据我想按照石景山各街道筛分出来帮忙写一段脚本可以自动将这些点位归到所属街道并应用到GIS展示中

以下是一个Python脚本，可以根据点位坐标和街道边界数据自动将点位归属到所属街道，并生成新的数据文件用于GIS展示。

import geopandas as gpd
from shapely.geometry import Point

# 读取点位数据
points = gpd.read_file('points.geojson')

# 读取街道边界数据
streets = gpd.read_file('streets.geojson')

# 创建空的街道列
points['street'] = ''

# 对每个点位，判断它属于哪个街道
for i, point in points.iterrows():
    for j, street in streets.iterrows():
        if Point(point['geometry']).within(street['geometry']):
            points.loc[i, 'street'] = street['name']
            break

# 保存新的数据文件
points.to_file('points_with_street.geojson', driver='GeoJSON')

请将脚本中的points.geojson和streets.geojson替换为您自己的点位和街道边界数据文件路径。在运行脚本后，它将自动为每个点位添加一个street列，并将其归属到所属街道。最终生成的points_with_street.geojson文件可以用于GIS展示