Description:
Performs a
mixed cross-self-correlation calculation on an input NDF image file. The resulting correlation
image/plot is stored to disk.
The cross-self-correlated image may be used to find faint diffuse objects for a given scale
length.
The circular exponential profile template used is of a size that optimises the search for galaxies of the
scale length requested by the user.
To reduce the influence of bright objects or cosmic rays; the user may elect to employ a cut out pixel
count value where any pixel found to be above that value is ignored. The cutout value is determined
by the user inputting a global background count value (available via HISTPEAK), the background
count standard deviation and the number of standard deviations above sky level at which the cutout
should occur.
The user is required to enter a value for the scale length of of the object(s) of interest and also the
image pixel size.
The method assumes some sort of symmetry is present in the objects detected but appears to work
well on a wide range of image types.
A border is present in the output image which is of the same width as the radius of the template. All
pixels within this border are assigned the value bad.
The correlation is optimised by making the template size 1.8x that of the galaxy scale length required.
This factor was determined from simulations by Phillipps and Davies at Cardiff.
Usage:
SELFCW IN OUT SCALE PSIZE BACK USEALL MULT [SIGMA] [NSIGMA]
Parameters:
BACK = _REAL (Read)
The modal pixel count value found in the input NDF. Units
counts.
IN = _NDF (Read)
The name of the NDF that is to be examined.
NSIGMA =
_REAL(Read)
The number of standard deviations above the sky level count value at which the
pixel count cutoff occurs.
MULT = _REAL (Read)
A multiplying factor used to modify
the output range.
OUT = _NDF (Write)
The name of the NDF that will be created.
PSIZE = _REAL (Read)
The size of each pixel in arc seconds. If the image contains a SKY
co-ordinate frame this value will be determined automatically.
SCALE = _REAL (Read)
The scale length of the galaxies being searched for. Units arc seconds.
SIGMA = _REAL
(Read)
The standard deviation of the pixel background count within the input NDF. Should
be determined using a routine such as HISTPEAK which ignores outliers.
USEALL =
_LOGICAL (Read)
Used to indicate whether a pixel count threshold is to be a applied when
calculating the self-correlation.
Examples:
selfcw in=p2 out=scp2 scale=15. psize=0.96
back=1000.2
useall=true
A self-correlation image, optimised for galaxies of a 15 arc second scale length, is
generated using image P2 as the input source image and SCP2 as the output image. The pixel size on
the image is .96 arc second and the background count value for the source image is 1000.2 selfcw
in=lsbg1 out=lsbg2 scale=8. back=444. useall=false
sigma=12. nsigma=4.
A self-correlation image, optimised for galaxies of a 8 arc second scale
length, is generated using image P2 as the input source image and SCP2 as the output image. The
background count value for the source image is 444 and the pixel size in arc seconds will be
determined from a SKY frame in the image’s WCS component if possible.
All pixels with a count value greater than
444.12.x4.
are excluded from the correlation calculations.
Notes:
It is assumed that the x and y axis pixels sizes
are the same size.
To establish the statistical significance of a detection, this application should be used in conjunction
with MIXUP to allow noise equivalent images to be generated and correlated thereby establishing a 3
sigma limit.