SAGA 9.4.2 | Tool Library Documentation

Resection (Terrestrial)


Description

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 point numberFFocal Length in mmDefault: 0.000000
Pixel Width (um)floating point numberWPixel 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 point numberppXPrincipal Point Offset in XDefault: 0.000000
Principal Point Offset in Y (pixels)floating point numberppYPrincipal Point Offset in YDefault: 0.000000
Provide Radial Distortion Parameters?booleanGIVE_DISTORTIONSDo you wish to provide the Radial Distortion Parameters?Default: 0
K1floating point numberK1K1Default: 0.000000
K2floating point numberK2K2Default: 0.000000
K3floating point numberK3K3Default: 0.000000
Xfloating point numberXcApproximate CoordinateDefault: 0.000000
Yfloating point numberYcApproximate CoordinateDefault: 0.000000
Zfloating point numberZcApproximate CoordinateDefault: 0.000000
Xfloating point numberXtApproximate CoordinateDefault: 0.000000
Yfloating point numberYtApproximate CoordinateDefault: 0.000000
Zfloating point numberZtApproximate CoordinateDefault: 0.000000
Output Text Filefile pathOUTPUT FILEThe file to write the Calculation Results to.-

Command Line


Usage: saga_cmd imagery_photogrammetry 0 [-POINTS ] [-F ] [-W ] [-EST_OFFSETS ] [-ppX ] [-ppY ] [-GIVE_DISTORTIONS ] [-K1 ] [-K2 ] [-K3 ] [-Xc ] [-Yc ] [-Zc ] [-Xt ] [-Yt ] [-Zt ] [-OUTPUT_FILE ]
  -POINTS:          	Measured Points (PC)
	point cloud, input
  -F:            	Focal Length (mm)
	floating point number
	Default: 0.000000
  -W:            	Pixel Width (um)
	floating point number
	Default: 0.000000
  -EST_OFFSETS:     	Estimate Principal Point Offsets?
	boolean
	Default: 0
  -ppX:          	Principal Point Offset in X (pixels)
	floating point number
	Default: 0.000000
  -ppY:          	Principal Point Offset in Y (pixels)
	floating point number
	Default: 0.000000
  -GIVE_DISTORTIONS:	Provide Radial Distortion Parameters?
	boolean
	Default: 0
  -K1:           	K1
	floating point number
	Default: 0.000000
  -K2:           	K2
	floating point number
	Default: 0.000000
  -K3:           	K3
	floating point number
	Default: 0.000000
  -Xc:           	X
	floating point number
	Default: 0.000000
  -Yc:           	Y
	floating point number
	Default: 0.000000
  -Zc:           	Z
	floating point number
	Default: 0.000000
  -Xt:           	X
	floating point number
	Default: 0.000000
  -Yt:           	Y
	floating point number
	Default: 0.000000
  -Zt:           	Z
	floating point number
	Default: 0.000000
  -OUTPUT_FILE:     	Output Text File
	file path