M33 combined 6cm+20cm catalog description (2016 July 31)

Text catalog: catalog_2016Jul31.txt.gz
FITS catalog: catalog_2016Jul31.fits
Islands (segmentation) image: m33-fixbeam-islands-g.fits.gz
Detection image: m33-fixbeam-combined-filter.fits.gz

The 6cm+20cm detection image was created by averaging the entire 6cm and 20cm bandpasses. The 20cm high and low frequency images were used, and 14 of the 16 125-MHz channels from the 6cm observations were used. (6cm channels 1 and 2 were contaminated by interference.) The 6cm and 20cm bands are equally weighted in the combined image. That's not optimal for a typical extragalactic source (which has a power law spectrum nu**alpha with alpha=-0.7) but is a reasonably "spectrum-neutral" weighting for the source detection. The resulting source lists do not appear to be a strong function of the weighting.

For this version of the catalog (as opposed to the 2015Feb03 version), the images at all frequencies were convolved to a fixed round beam with FWHM 5.9 arcsec. That reduces the resolution of the final image, but it improves the consistency of the multi-resolution filtering step. So I think this should produce more accurate flux densities and (especially) spectral indices than the previous catalog.

A multi-resolution (median pyramid) image processing algorithm was used both for source detection and to subtract the background from the image to make the source detection threshold more uniform. The images without background subtraction can be compared to the original summed image to see the effect of the background subtraction.

The island map (segmentation map) shows where sources where detected. The numbers in the island map correspond to the "Island" field in the catalog. Islands were detected by thresholding the multi-resolution image. Islands detected at more than one resolution are merged. It is possible for several sources detected separately at high resolution to be embedded in one larger source at low resolution; in that case the morphology and sizes of the high-resolution islands are used to determine how to break the large island into pieces. Padding is added around the islands to give enough context for an accurate fit of the fluxes, positions and sizes. Rarely the padding causes convergence problems, in which case it is reduced to get a better fit. The FITS file with the island numbers is available for download above.

The pixels in each island are fitted with elliptical Gaussian source models to determine sizes, positions, and peak and integrated flux densities. For most islands a constant background offset is also fitted (although that is omitted when it causes convergence problems). The data are fitted simultaneously in 16 frequency bands (2 20cm bands and 14 6cm bands). The Gaussian is constrained to have the same position and source size in all bands, with the flux and background varying between bands. The fitted size quoted in the catalog is includes the round 5.9 arcsec FWHM beam size. For an unresolved source the fitted size will be approximately equal to 5.9 arcsec. Note that the fitted size is not constrained to be larger than this value; sources with sizes smaller than the beam PSF will have integrated fluxes smaller than their peak fluxes.

In addition to the Gaussian fits, the catalog also includes integrated and peak fluxes for the island. The peak is just the brightest pixel in the island, and the integrated flux is the sum of all the flux within the island with a PSF-dependent correction factor (close to unity) for flux that spills outside the island. Usually the island fluxes are are noisier than the Gaussian fits but they work better when the source is not well modeled as a Gaussian. The image with the Gaussian source models from the catalog subtracted can be used to identify sources where the Gaussian fits are poor.

The images with islands overlaid show the positions of the sources. Red islands have source entries in the catalog; green islands either are too faint (less than 4-sigma) or failed in fitting, usually because the region was crowded.

When islands touch, the Gaussian fits to all the islands in the group are carried out simultaneously, so the effects of overlapping sources is properly included. That is an improvement over the previous version of the catalog. The resulting fits in crowded regions are much improved over the older catalog.

Catalog column descriptions

The 2016Jul31 M33 catalog from 6cm+20cm data includes 2024 sources detected at 4-sigma or greater in island or in spectral index fit. It is available for download in text or FITS format. The text catalog has comments at the top that describe the contents. Here is a description of the columns:

NumberColumnDescription
1-3RA2000RA (hh mm ss.sss)
4-6Dec2000Declination (hh mm ss.sss)
7RA_deg(deg) RA
8Dec_deg(deg) Declination
9X(pixels) X pixel position of fitted peak
10Y(pixels) Y pixel position of fitted peak
11Fint(uJy) Flux at mean frequency nu0 = 3.36125 GHz
12FiErr(uJy) Error on Fint
13SpindSpectral index, Fnu = Fint*(nu/nu0)**Spind
14SpErrError on SpErr
15chisqChi-square for power-law fit to flux values in 4 bands
16major(arcsec) Major axis FWHM (fitted, includes Gaussian beam with FWHM = 5.9 arcsec)
17minor(arcsec) Minor axis FWHM (fitted, includes Gaussian beam with FWHM = 5.9 arcsec)
18PA(deg) Position angle of major axis
19IslandIsland number (integer)
20NpixNumber of pixels in island
21Isl_Fp(uJy) Peak pixel in island
22I_Fpsig(uJy) Noise in Isl_Fp
23Isl_Fi(uJy) Integrated flux in island
24I_Fisig(uJy) Noise in Isl_Fi
251.38GHz(uJy) Integrated flux in 20cm low-frequency band (1.3775 GHz)
261.81GHz(uJy) Integrated flux in 20cm high-frequency band (1.8135 GHz)
274.68Ghz(uJy) Integrated flux in 6cm low-frequency band (4.6790 GHz)
285.58Ghz(uJy) Integrated flux in 6cm high-frequency band (5.5750 GHz)
29sig1.38(uJy) Noise in F(1.38GHz)
30sig1.81(uJy) Noise in F(1.81GHz)
31sig4.68(uJy) Noise in F(4.68GHz)
32sig5.58(uJy) Noise in F(5.58GHz)
33Bkgd(uJy) Constant background in island fit

Richard L. White, rlw@stsci.edu
2016 July 31