SAGA-GIS Tool Library Documentation (v9.1.2)

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.

Parameters

 NameTypeIdentifierDescriptionConstraints
InputMeasured Points (PC)Point Cloud, inputPOINTSList of Measured Points as PC-
OptionsFocal Length (mm)Floating pointFFocal Length in mmDefault: 0.000000
Pixel Width (um)Floating pointWPixel Width in micro metersDefault: 0.000000
Estimate Principal Point Offsets?BooleanEST_OFFSETSDo you wish to estimate Principal Point Offsets?Default: 0
Principal Point Offset in X (pixels)Floating pointppXPrincipal Point Offset in XDefault: 0.000000
Principal Point Offset in Y (pixels)Floating pointppYPrincipal Point Offset in YDefault: 0.000000
Provide Radial Distortion Parameters?BooleanGIVE_DISTORTIONSDo you wish to provide the Radial Distortion Parameters?Default: 0
K1Floating pointK1K1Default: 0.000000
K2Floating pointK2K2Default: 0.000000
K3Floating pointK3K3Default: 0.000000
XFloating pointXcApproximate CoordinateDefault: 0.000000
YFloating pointYcApproximate CoordinateDefault: 0.000000
ZFloating pointZcApproximate CoordinateDefault: 0.000000
XFloating pointXtApproximate CoordinateDefault: 0.000000
YFloating pointYtApproximate CoordinateDefault: 0.000000
ZFloating pointZtApproximate CoordinateDefault: 0.000000
Output Text FileFile pathOUTPUT FILEThe file to write the Calculation Results to.-

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