FRD2P

Description

Applies a reduction to pole filter to a potential field image. The filter is applied in the frequency domain. A 2-D Fast Fourier transformation (FFT) is first applied to the image. After filtering, the image is transformed back to the spatial domain.

Parameters

Name	Type	Caption	Length	Value range
FILI *	String	Input file name	1 - 192
FILO *	String	Output file name	1 - 192
DBIC *	Integer	Input potential field channel	1 - 1
DBOC *	Integer	Output filtered image channel	1 - 1
DBIW	Integer	Raster input window	0 - 4	Xoffset, Yoffset, Xsize, Ysize
DBOW	Integer	Raster output window	0 - 4	Xoffset, Yoffset, Xsize, Ysize
SRCANG *	Float	Source field direction angles	2 - 2	0 - 180 Default: 30
REGANG *	Float	Regional field direction angles	2 - 2
FILB	String	Buffer file name	0 - 192
REPORT	String	Report mode	0 - 192	Quick links
MONITOR	String	Monitor mode	0 - 3	ON, OFF Default: ON

* Required parameter

Parameter descriptions

FILI

Specifies the name of the PCIDSK file from which image data is read.

FILO

Specifies the name of the PCIDSK file to receive the transformed image data.

The output file can be the same as the input file. The specified output file must exist before running FRD2P.

DBIC

Specifies the input channels containing the potential fields to use in the transformation.

DBOC

Specifies the output channel to receive the filtered image data.

DBIW

Specifies the raster window (Xoffset, Yoffset, Xsize, Ysize) read from the input image. If this parameter is not specified, the entire image is processed by default.

Xoffset, Yoffset define the upper-left starting pixel coordinates of the window. Xsize is the number of pixels that define the window width. Ysize is the number of lines that define the window height.

DBOW

Specifies the raster window (Xoffset, Yoffset, Xsize, Ysize) to be output. If this parameter is not specified, the entire layer is output by default.

Xoffset, Yoffset define the upper-left starting pixel coordinates of the window. Xsize is the number of pixels that define the window width. Ysize is the number of lines that define the window height.

SRCANG

Specifies, in degrees, the inclination and azimuth angles of the source magnetic field. Default values for both angles are 30.0.

REGANG

Specifies, in degrees, the inclination and azimuth angles of the regional magnetic field. Default values for both angles are 30.0.

FILB

Optionally specifies the name of the temporary buffer file.

Because this function uses large amounts of RAM, a temporary buffer file is automatically created for large images. This parameter records the name of the buffer file created by FCONT. Typically, no user interaction is required for this parameter.

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

MONITOR

The program progress can be monitored by printing the percentage of processing completed. A system parameter, MONITOR, controls this activity.

Available options are:

• ON: Turns the monitor ON (default)
• OFF: Turns the monitor OFF (recommended if running in batch or background mode)

Details

FRD2P applies the reduction to pole filter to a potential field image. The filter is applied in the frequency domain. The image is transformed using a 2-D Fast Fourier transformation (FFT), then transformed back to the spatial domain after filtering.

Because this function does not require the entire image to be loaded into memory, there is no limit on the size of input image. If necessary, the function automatically creates a buffer file, if FILB is not specified. The buffer file name is recorded in FILB for efficiently reusing the buffer in another run.

This function is most efficient if the input window has dimensions of a power of 2; otherwise, FCONT must pad extra rows and columns to force the image dimensions a power of 2.

Example

EASI>fili	=	'magfield.pix'	! input file
EASI>filo	=	'magfield.pix'	! output file
EASI>dbic	=	1	! potential fields channel
EASI>dboc	=	2	! output filtered channel
EASI>dbiw	=		! process entire image
EASI>dbow	=		! output entire image
EASI>srcang	=	90,0	! source field direction angles
EASI>regang	=	45,60	! regional field direction angles
EASI>filb	=		! default, no buffer file specified

EASI>run FRD2P

Algorithm

The input image is first transformed to frequency domain using 2-D FFT. The dimensions of the transformed image are power of 2 and are at least as large as the input image dimensions. After applying the filter, the frequency image is transformed back to the spatial domain and truncated to the input image size.

The reduction to pole filter has the following form:

        f^2 / [c1 f + i(a1 u + b1 v)] [c2 f + i(a2 u + b2 v)]

The filter gain at u=v=0 is undefined. A value of 1.0 is inserted. The resolution of u, v are given as:

        Delta_u = 1/(SizeU * Delta_x)
        Delta_v = 1/(SizeV * Delta_y)