Tool Resection (Terrestrial)

Single Image Spatial Resection (Terrestrial): from at least 3 image points with known coordinates, the Cardan angles of the image orientation and the coordinates of the perspective center are calculated by a least-squares adjustment. The Cardan angles refer to the following transformation between image coordinates (x_image) and global coordinates (X_Global):
x_image = R_1(omega) * R_2(kappa) * R_3(alpha) * X_Global
Here R_1, R_2, R_3 denote rotation matrices of a right-handed (passive) coordinate transformation.

The inputs consist of a point cloud containing the identical points with their pixel indices as additional attributes. The origin of pixels is in the lower left corner of the image. The interior orientation parameters of the camera is to be provided. These include Focal Length (mm), Pixel Size (um), Principal Point Offsets (pixels) and optionally the Radial Distortion Parameters. The distortion model being used is as follows:
x_d = x_u (1 - dR)
y_d = y_u (1 - dR)
where, dR = K1 * r_u ^ 2 + K2 * r_u ^ 4 + K3 * r_u ^ 6,
r_u ^ 2 = x_u ^ 2 + y_u ^ 2,
x_u, y_u are the undistorted (corrected) image coordinates in mm,
x_d, y_d are the distorted (observed) image coordinates in mm,
K1 is in [mm ^ -2], K2 is in [mm ^ -4], K3 is in [mm ^ -6].
Approximate coordinates in [m] for the Projection Center and the center of image are also to be provided.

Following the adjustment the results are written to a text file. The main contents of the text file are: for each iteration the Sum of Squared Residuals, the A-Posteriori Standard Deviation (Sigma Naught) and the Condition of Normal Matrix, and the final estimated exterior orientation parameters (Xc, Yc, Zc, Omega, Kappa, Alpha).

Optionally the Principal Point Offsets can be estimated. This requires at least 4 image points as input.

[Warning] The tool is dedicated to the terrestrial case which means the viewing direction must be sufficiently different from the vertical direction.
The estimation of Principal Point Offsets is not recommended if the condition of the Normal Matrix is less than 10 ^ -7.

• Author: HfT Stuttgart (c) 2013
• Menu: Imagery|Photogrammetry

Parameters

Command-line

Usage: saga_cmd imagery_photogrammetry 0 [-POINTS <str>] [-F <double>] [-W <double>] [-EST_OFFSETS <str>] [-ppX <double>] [-ppY <double>] [-GIVE_DISTORTIONS <str>] [-K1 <double>] [-K2 <double>] [-K3 <double>] [-Xc <double>] [-Yc <double>] [-Zc <double>] [-Xt <double>] [-Yt <double>] [-Zt <double>] [-OUTPUT FILE <str>]
  -POINTS:<str>          	Measured Points (PC)
	Point Cloud, input
  -F:<double>            	Focal Length (mm)
	Floating point
	Default: 0.000000
  -W:<double>            	Pixel Width (um)
	Floating point
	Default: 0.000000
  -EST_OFFSETS:<str>     	Estimate Principal Point Offsets?
	Boolean
	Default: 0
  -ppX:<double>          	Principal Point Offset in X (pixels)
	Floating point
	Default: 0.000000
  -ppY:<double>          	Principal Point Offset in Y (pixels)
	Floating point
	Default: 0.000000
  -GIVE_DISTORTIONS:<str>	Provide Radial Distortion Parameters?
	Boolean
	Default: 0
  -K1:<double>           	K1
	Floating point
	Default: 0.000000
  -K2:<double>           	K2
	Floating point
	Default: 0.000000
  -K3:<double>           	K3
	Floating point
	Default: 0.000000
  -Xc:<double>           	X
	Floating point
	Default: 0.000000
  -Yc:<double>           	Y
	Floating point
	Default: 0.000000
  -Zc:<double>           	Z
	Floating point
	Default: 0.000000
  -Xt:<double>           	X
	Floating point
	Default: 0.000000
  -Yt:<double>           	Y
	Floating point
	Default: 0.000000
  -Zt:<double>           	Z
	Floating point
	Default: 0.000000
  -OUTPUT FILE:<str>     	Output Text File
	File path