如何根据经纬度把h5文件对应点坐标得到

要根据经纬度获取h5文件中的像元坐标，你可以使用一些地理信息处理库来实现。下面是一个示例代码，使用了`pyproj`库来进行坐标转换和投影操作：

import pandas as pd
import h5py
from pyproj import Proj, transform

# 读取Excel表格中的经纬度数据
df = pd.read_excel('data.xlsx')  # 假设Excel文件名为data.xlsx，需要与代码在同一目录下
latitudes = df['纬度列名']
longitudes = df['经度列名']

# 定义投影坐标系
in_proj = Proj(init='epsg:4326')  # 输入经纬度坐标系，这里使用WGS84坐标系
out_proj = Proj(init='epsg:XXXX')  # 输出坐标系，根据你的具体数据情况指定，例如UTM投影坐标系

# 打开h5文件
h5_file = h5py.File('data.h5', 'r')  # 假设h5文件名为data.h5，需要与代码在同一目录下

# 提取像元值
pixel_values = []
for lat, lon in zip(latitudes, longitudes):
    # 经纬度转换为目标坐标系的坐标
    x, y = transform(in_proj, out_proj, lon, lat)

    # 将坐标四舍五入为整数，表示像元坐标
    pixel_coords = int(round(x)), int(round(y))

    # 从h5文件中读取像元值
    pixel_value = h5_file['dataset_name'][pixel_coords[1], pixel_coords[0]]  # dataset_name是h5文件中的数据集名字

    pixel_values.append(pixel_value)

# 关闭h5文件
h5_file.close()

# 打印提取到的像元值
print(pixel_values)

上述代码中的`'epsg:XXXX'`需要根据你的实际数据情况指定输出坐标系，例如使用UTM投影坐标系。你可以根据需要修改该参数。

希望这能帮助到你！如果还有其他问题，请随时提问。