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

Question

I am looking for the rasterio function that does forward transformation from geodetic coordinate system to image coordinate system based on RPC. Can anyone help me within this context?

Answer 1

rasterio currently doesn't support using RPCs to transform from world coordinates to image coordinates, like GDAL does. See the discussion on this topic here: https://github.com/mapbox/rasterio/issues/410

However, I can recommend a python library that implements a vectorized transform to/from world/image coordinates: https://github.com/cmla/rpcm

The README of this library only shows the command-line usage of the tool, but you can call it directly in 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)