如何遍历shapefile(.shp)并获得每个坐标的高度？

import shapefile
from osgeo import gdal
import numpy as np

# 打开shapefile文件
sf = shapefile.Reader('your_shapefile.shp')

# 获取要素和记录
shapes = sf.shapes()
records = sf.records()

# 遍历要素和记录
for shape, record in zip(shapes, records):
    # 获取要素的几何形状
    points = shape.points
    
    # 遍历坐标点
    for point in points:
        # 获取坐标点的经度和纬度
        lon, lat = point
        
        # 获取坐标点的高度
        # 这里使用示例数据，实际情况需要根据数据源获取高度信息
        elevation = get_elevation(lon, lat)
        
        # 处理高度数据，例如打印每个坐标点的高度
        print(f"坐标点({lon}, {lat})的高度为：{elevation}")

# 获取坐标点的高度信息（示例函数）
def get_elevation(lon, lat):
    # 使用gdal库读取高程数据
    dataset = gdal.Open('your_elevation_data.tif')
    band = dataset.GetRasterBand(1)
    
    # 根据经纬度获取像素值
    gt = dataset.GetGeoTransform()
    x = int((lon - gt[0]) / gt[1])
    y = int((lat - gt[3]) / gt[5])
    elevation = band.ReadAsArray(x, y, 1, 1)[0][0]
    
    return elevation