profoundSegimExtend: Find Additional Image Segments

Description

This function takes a target segmentation map, and finds additional sources on a target image. This might be useful where you want to do matched aperture photometry, but there is the possibility that the target image might have genuine additional sources that are not in your current segmentation map. This might cause issues with the sky estimation etc (even if you do not actually care about getting photometry for these additional sources).

Usage

profoundSegimExtend(image = NULL, segim = NULL, mask = segim, ...)

Arguments

image

Numeric matrix; required, the image we want to analyse. Note, image NAs are treated as masked pixels.

segim

Integer matrix; required, the segmentation map of the image. This matrix *must* be the same dimensions as image.

mask

Boolean matrix; optional, parts of the image to mask out (i.e. ignore), where 1 means mask out and 0 means use for analysis. If provided, this matrix *must* be the same dimensions as image.

…

Further arguments to be passed to profoundProFound.

Value

Integer matrix; the segmentation map matched pixel by pixel to image and segim. Newly identified segments are appended in number to the input segim, so if the maximum previous segment ID was 10, the new sources would start from 11 etc.

Details

This function is mostly useful if you are chaining together a script of detection images based on certain bands that are then run on data which might have very different (and previously undetected) astrophysical sources.

Examples

Run this code

# NOT RUN {
image=readFITS(system.file("extdata", 'VIKING/mystery_VIKING_Z.fits', package="ProFound"))

profound=profoundProFound(image, magzero=30, plot=TRUE)

#Now we remove the brightest central source:
profound$segim[profound$segim==1]=0

#And now we find it again:
segim_extend=profoundSegimExtend(image=image, segim=profound$segim, skycut=2)
profoundSegimPlot(image=image$imDat, segim=segim_extend)
# }