基于RPC的将大地测量转换为图像坐标系

Question

我正在寻找实现从大地坐标系到基于RPC的图像坐标系的正向转换的rasterio函数。在这种情况下，谁能帮助我？

Answer 1

rasterio当前不支持使用RPC将世界坐标转换为图像坐标，就像GDAL一样。请参阅此处有关此主题的讨论：https://github.com/mapbox/rasterio/issues/410

但是，我可以推荐一个Python库，该库实现与世界/图像坐标之间的向量化转换：https://github.com/cmla/rpcm

此库的自述文件仅显示了该工具的命令行用法，但是您可以直接在python中调用它：

import numpy as np
from rpcm import rpc_from_geotiff

rpc = rpc_from_geotiff("path_to_your_file.tiff")

lon = np.array([..., ...])
lat = np.array([..., ...])
x, y = rpc.projection(lon, lat)