基于RPC的大地坐标到图像坐标系的转换

Question

我正在寻找基于 RPC 从大地坐标系到图像坐标系的正向变换的光栅 function。 任何人都可以在这种情况下帮助我吗？

Answer 1

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

但是，我可以推荐一个 python 库，它实现与世界/图像坐标的矢量化变换： https://github.com/cmla/rpcm

这个库的 README 只显示了该工具的命令行用法，但您可以直接在 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)