PNTREFN

Description

PNTREFN refines tie points (TP) and ground control points (GCP) in a CATALYST Professional OrthoEngine project file.

Parameters

Name	Type	Caption	Length	Value range
OEPROJI *	String	Input OrthoEngine project file	1 - 192
REFNMODE *	String	Refinement mode	4 - 192	PROG \| DIRECT Default: PROG
GCPRES	Float	Maximum GCP residual	0 - 1	-5 - 1000
TPRES	Float	Maximum TP residual	0 - 1	-5 - 1000
OEPROJO *	String	Output OrthoEngine project file	1 - 192
POINTOPT *	String	Options for rejected points	3 - 192	INACT \| DEL Default: INACT
MINTP	Integer	Minimum number of TPs	0 - 1
REFNOPT	String	Options for rejected points	0 - 192
IMSTAT	Float	Output PNTREFN statistics	0 - 48
REPORT	String	Report mode	0 - 192	Quick links

* Required parameter

Parameter descriptions

OEPROJI

The path and file name of the OrthoEngine project file (*.prj) containing the points you want to refine.

REFNMODE

The mode of refinement to use: progressive or direct.

In progressive (PROG) mode, points greater than 1,000, 500, 100, 20, and 10 pixels, respectively, are removed iteratively before applying the maximum GCP residual, maximum TP residual, or both, as specified for GCPRES and TPRES, respectively.

In direct (DIRECT) mode, the maximum residual specified is applied exclusively.

GCPRES

The maximum residual for x and y per GCP.

For example, to refine the OrthoEngine project until all GCPs have a residual of two or less in both the x and y direction, enter 2.

If the entered value is negative, it uses outlier method to remove points. Any points within an image have residuals from the mean value greater than -GCPRES * standard deviation of the image would be removed.

TPRES

The maximum residual for x and y per TP.

For example, to refine the OrthoEngine project until all TPs have a residual of two or less in both the x and y direction, enter 2.

If the entered value is negative, it uses outlier method to remove points. Any points within an image have residuals from the mean value greater than -TPRES * standard deviation of the image would be removed.

OEPROJO

The path and file name of the output OrthoEngine project file to which to write the refined TPs, GCPs, or both.

You can specify your input project as your output project; however, by doings so, the input file will be overwritten with the new information.

POINTOPT

The disposition of the rejected points.

That is, on completion of refinement, you can elect to either render the rejected points inactive (INACT) or delete (DEL) them.

MINTP

The minimum number of TPs per image.

Images with fewer than the minimum number of TPs before and after refinement will be rendered inactive.

REFNOPT

Refinement options. Current supported options are UPDATE_ANGLES, NO_CHECK and NO_MESSAGE.

UPDATE_ANGLES: Update input omega, phi and kappa values with the computed kappa values
NO_CHECK: To turn off checking for removed points.
NO_MESSAGE: To disable messages from bundle adjustment.
MINGCP=#: To keep the minimum number of GCPs per image. Value can be an integer (such as MINGCP=10) or an integer with % (such as MINGCP=10%).
MAXGCPRES=#: This option is used together with MINGCP. It allows the user to set the maximum residual for each image together with MINGCP to avoid large blunders. This program will continue to remove GCPs from an image which have residuals greater than MAXGCPRES even when it has reached the minimum number of GCPs.

IMSTAT

On completion of processing, PNTREFN updates IMSTAT with information that may be useful in a scripting environment, as follows:

first IMSTAT array element: Original number of GCPs.
second IMSTAT array element: Original number of TPs.
third IMSTAT array element: Initial TP x-RMS.
fourth IMSTAT array element: Initial TP y-RMS.
fifth IMSTAT array element: Initial maximum TP x-residual.
sixth IMSTAT array element: Initial maximum TP y-residual.
seventh IMSTAT array element: Refined TP x-RMS.
eighth IMSTAT array element: Refined TP y-RMS.
ninth IMSTAT array element: Refined maximum TP x- residual.
tenth IMSTAT array element: Refined maximum TP y-residual.
eleventh IMSTAT array element: Initial GCP x-RMS.
twelfth IMSTAT array element: Initial GCP y-RMS.
thirteenth IMSTAT array element: Initial maximum GCP x-residual.
fourteenth IMSTAT array element: Initial maximum GCP y-residual.
fifteenth IMSTAT array element: Refined GCP x-RMS.
sixteenth IMSTAT array element: Refined GCP y-RMS.
seventeenth IMSTAT array element: Refined maximum GCP x-residual.
eighteenth IMSTAT array element: Refined maximum GCP y-residual.

REPORT

Specifies where to direct the generated report.

Available options are:

TERM: generates a report on the terminal (default)
DISK: generates a report on file "IMPRPT.LST"
OFF: turns off report generation
<filename>: appends a report to the specified file

Details

PNTREFN refines TPs and GCPs in an OrthoEngine project file.

The algorithm removes points in either a progressive or direct mode according to the specified residual. In progressive mode, large residuals are removed in the following steps:

Remove residuals greater than 1,000 pixels
Remove residuals greater than 500 pixels
Remove residuals greater than 100 pixels
Remove residuals greater than 20 pixels
Remove residuals greater than 10 pixels
Remove residuals according to the specified residual

In direct mode, the specified residual is applied exclusively.

Processing is complete when the specified residual is reached.

Example

In the following example, TPs and GCPs are designated inactive until the x and y residual for each meets the specified GCP and TP residual.

EASI>OEPROJI	=	"airphoto_model.prj"
EASI>REFNMODE =	"PROG"  ! Progressive mode
EASI>GCPRES	=	1  ! Refine until remaining GCPs have an x-residual and a y-residual of 1
EASI>TPRES	=	2  ! Refine until remaining TPs have an x-residual and a y-residual of 2
EASI>OEPROJO	=	"refined.prj"
EASI>POINTOPT =      "INACT"        ! Designate refined points inactive
EASI>MINTP	=	10 ! Designate images with fewer than 10 TPs inactive

EASI>R PNTREFN

Environments	PYTHON :: EASI
Quick links	Description :: Parameters :: Parameter descriptions :: Details :: Example