profoundSegimExtend: Find Additional Image Segments
In ProFound: Photometry Tools

Description Usage Arguments Details Value Author(s) See Also Examples

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).

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

`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`.

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.

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.

Aaron Robotham

profoundProFound, profoundMakeSegim

## 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)

## End(Not run)