Synopsis
Identifies groups of pixels above/below threshold
Syntax
dmimgblob infile outfile threshold [srconly] [clobber] [verbose]
Description
`dmimgblob' will identify groups of contigious pixels that are above (and optionally below) some user selected threshold. The output image is a mask where the pixels in the mask respresent which unique set the pixel belongs to. Pixels above threshold are given values 1 to N (where N is the number of blobs found). This can be used as a crude source detection.
When the 'srconly' parameter is set to no (default), then pixels which fall below the threshold are also grouped. These (presumably) background pixels are given group numbers in the output mask from -1 to -M (where M is the number of background blobs).
Any Null/NaN valued pixels or pixels outside the data-subspace are assigned a mask value = 0.
Note: This tool uses a recursive algorithm that may exhaust the stack for particularly large datasets. If you get a 'Segmentation violation' please try increase the 'stacksize' (see 'man limit'); (t)csh: limit stacksize unlimited; (ba)sh: ulimit -s unlimited.
Examples
Example 1
dmimgblob infile=img.fits outfile=mask.fits threshold=10
The file 'mask.fits' will contain a map of where the input image, img.fits, was >= 10 (in whatever units the input image is). Contigious pixel < 10 are labeled with negative values in the mask image.
Example 2
dmimgblob infile=img.fits outfile=mask.fits threshold=10 srconly=yes
Similar to above except that only values >= the threshold are included (as are Null/NaN or data-subpsace excluded values with a mask value = 0).
Parameters
name | type | ftype | def | min | max | reqd |
---|---|---|---|---|---|---|
infile | file | input | yes | |||
outfile | file | output | yes | |||
threshold | real | yes | ||||
srconly | boolean | no | ||||
clobber | boolean | no | ||||
verbose | integer | 0 | 0 | 5 |
Detailed Parameter Descriptions
Parameter=infile (file required filetype=input)
The input image.
The input 2D image.
Parameter=outfile (file required filetype=output)
The output file name
The output mask file (long data-type). Each unique blob of pixels in the input image above threshold are identified with unique values in the output mask file.
Parameter=threshold (real required)
Threshold to apply to image
Input image values >= this value are masked w/ positive values in the output mask. Value < this value are optionally identified with negative values in the output.
Parameter=srconly (boolean default=no)
Include values only >= threshold?
If set to 'yes' then only values >= the threshold are included. If set to 'no' then values < are also included.
Parameter=clobber (boolean default=no)
Remove output if it exists?
Used to specify whether or not to clobber existing file that has the same name as the specified output file
Parameter=verbose (integer default=0 min=0 max=5)
The tool chatter level
Verbose can be from 0 to 5, generating different amounts of debugging output.
Changes in CIAO 4.15
-
Updated the tool to recognize real-valued NULLs and to recognize valid data ranges.
Bugs
There are no known bugs for this tool.
See Also
- dm
- dmfiltering
- tools
- dmappend, dmcontour, dmellipse, dmfilth, dmimg2jpg, dmimgadapt, dmimgblob, dmimgcalc, dmimgdist, dmimgfilt, dmimghist, dmimghull, dmimglasso, dmimgpick, dmimgpm, dmimgproject, dmimgreproject, dmimgthresh, dmmaskbin, dmmaskfill, dmnautilus, dmregrid, dmregrid2, dmstat, evalpos, imgmoment, mean_energy_map, pileup_map