def approximate_rpc_as_projective(rpc_model,
col_range,
lin_range,
alt_range,
verbose=False):
"""
Returns a least-square approximation of the RPC functions as a projection
matrix. The approximation is optimized on a sampling of the 3D region
defined by the altitudes in alt_range and the image tile defined by
col_range and lin_range.
"""
### step 1: generate cartesian coordinates of 3d points used to fit the
### best projection matrix
# get mesh points and convert them to geodetic then to geocentric
# coordinates
cols, lins, alts = generate_point_mesh(col_range, lin_range, alt_range)
lons, lats, alts = rpc_model.direct_estimate(cols, lins, alts)
x, y, z = geographiclib.geodetic_to_geocentric(lats, lons, alts)
### step 2: estimate the camera projection matrix from corresponding
# 3-space and image entities
world_points = np.vstack([x, y, z]).T
image_points = np.vstack([cols, lins]).T
P = estimation.camera_matrix(world_points, image_points)
### step 3: for debug, test the approximation error
if verbose:
# compute the projection error made by the computed matrix P, on the
# used learning points
colPROJ = np.zeros(len(x))
linPROJ = np.zeros(len(x))
for i in range(len(x)):
v = np.dot(P, [[x[i]], [y[i]], [z[i]], [1]])
colPROJ[i] = v[0] / v[2]
linPROJ[i] = v[1] / v[2]
print('approximate_rpc_as_projective: (min, max, mean)')
print_distance_between_vectors(cols, colPROJ, 'cols')
print_distance_between_vectors(lins, linPROJ, 'rows')
return P
def approximate_rpc_as_projective(rpc_model, col_range, lin_range, alt_range,
verbose=False):
"""
Returns a least-square approximation of the RPC functions as a projection
matrix. The approximation is optimized on a sampling of the 3D region
defined by the altitudes in alt_range and the image tile defined by
col_range and lin_range.
"""
### step 1: generate cartesian coordinates of 3d points used to fit the
### best projection matrix
# get mesh points and convert them to geodetic then to geocentric
# coordinates
cols, lins, alts = generate_point_mesh(col_range, lin_range, alt_range)
lons, lats, alts = rpc_model.direct_estimate(cols, lins, alts)
x, y, z = geographiclib.geodetic_to_geocentric(lats, lons, alts)
### step 2: estimate the camera projection matrix from corresponding
# 3-space and image entities
world_points = np.vstack([x, y, z]).T
image_points = np.vstack([cols, lins]).T
P = estimation.camera_matrix(world_points, image_points)
### step 3: for debug, test the approximation error
if verbose:
# compute the projection error made by the computed matrix P, on the
# used learning points
colPROJ = np.zeros(len(x))
linPROJ = np.zeros(len(x))
for i in range(len(x)):
v = np.dot(P, [[x[i]],[y[i]],[z[i]],[1]])
colPROJ[i] = v[0]/v[2]
linPROJ[i] = v[1]/v[2]
print('approximate_rpc_as_projective: (min, max, mean)')
print_distance_between_vectors(cols, colPROJ, 'cols')
print_distance_between_vectors(lins, linPROJ, 'rows')
return P