问题描述
当前,我正在处理一些tiff文件,这些文件包含有关降水的信息,并且我想将它们与预先存在的数据框组合在一起。 此数据框具有提供纬度和经度的列,最终我希望能够将tiff文件的降水值与数据框合并以得到正确的纬度和经度。 如何找到与正确的经度和纬度相对应的降水值?
解决方法
有一个来自GDAL Python绑定存储库的示例,val_at_coord.py
返回给定地理位置的栅格值:
https://github.com/OSGeo/gdal/blob/master/gdal/swig/python/samples/val_at_coord.py
用法:
val_at_coord.py [-display_xy] [longitude latitude | -coordtype=georef X Y] filename
默认情况下,前两个参数应该是 经度,纬度顺序。如果在之前指定了coordtype = georef 接下来的2个值将被解释为 数据集空间参考系统。