基於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)