Outputting the number of pixels in an island

[nmmaddox]

I see it is possible to output the total flux contained within an island. Is it possible to also output the number of pixels in the island? It doesn't need to be any sort of shape, just the number of pixels.

[darafferty]

Not directly, unfortunately. But it is stored in the Image object. So, for example, if you use PyBDSF within a Python script like so:

import bdsf
img = bdsf.process_image('image.fits')

then you can access the island size (number of pixels) for an island with:

isl = img.islands[island_id]
isl_npix = isl.size_active

If you use the interactive pybdsf shell, the Image object is hidden but can be accessed as _img.

[nmmaddox]

OK, thanks very much for this.