partition_skymodel_by_flux

partition_skymodel_by_flux(sc, model, flux_threshold=-inf)

Partition skymodel according to flux

Bright skycomponents are put into a SkyModel as a list, and weak skycomponents are inserted into SkyModel as an image.

Parameters:

• sc – List of skycomponents

• model – Model image

• flux_threshold –

Returns:

SkyModel

For example:

fluxes = numpy.linspace(0, 1.0, 11)
sc = [create_skycomponent(direction=phasecentre, flux=numpy.array([[f]]), frequency=frequency,
                          polarisation_frame=PolarisationFrame('stokesI')) for f in fluxes]

sm = partition_skymodel_by_flux(sc, model, flux_threshold=0.31)
assert len(sm.components) == 7, len(sm.components)