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

Question

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

Answer 1

rasterio当前不像GDAL那样支持使用RPC从世界坐标转换为图像坐标。请在此处查看有关此主题的讨论：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)