写一段Python代码 从nc文件中提取固定经纬度的 U、V值

为了从nc文件中提取特定经纬度的U、V值，你可以按照以下步骤操作：

1. 定义一个函数来匹配经纬度并获取对应的U和V值：

def extract_values_at(lat, lon, nc_data):
    """
    根据给定的经纬度从nc_data中提取U和V值

    参数:
    lat (float): 目标纬度
    lon (float): 目标经度
    nc_data (xr.Dataset): 使用xr.open_dataset()打开的nc数据

    返回:
    tuple: 包含U和V值的元组，如果找不到对应位置，则返回None
    """
    try:
        lat_idx, lon_idx = matching(lat, lon, nc_data)
        u_val = nc_data['u'][lat_idx, lon_idx]
        v_val = nc_data['v'][lat_idx, lon_idx]
        return u_val, v_val
    except Exception as e:
        print(f"无法提取({lat}, {lon})的U/V值: {e}")
        return None, None

2. 遍历站点列表，调用上述函数提取每个站点的U、V值：

# 假设你已经有了站点列表，包含经纬度和对应的站点代号
stations_data = [(a1, b1, c1, 'station_1'), ...]

for station in stations_data:
    lat, lon, _, _ = station  # 取出纬度和经度
    u_val, v_val = extract_values_at(lat, lon, nc_dataset)  # 传递到函数中
    if u_val is not None and v_val is not None:
        # 处理提取到的U和V值，如存储在一个新的字典或数据结构中
        result_dict[station[-1]] = {'U': u_val, 'V': v_val}

这里假设`nc_dataset`已经通过`xr.open_dataset()`打开了nc文件。