Note you should not regrid the mosaiced image. In regridding, the information that describes the mosaicing process is lost, which will make life difficult for you.
Typical inputs to regrid would be
REGRID | |
in=lmc.sd | Input single dish image. |
tin=lmc.mosaic | The mosaic giving the coordinate system |
that we want to regrid to. | |
out=lmc.sd_regrid | The output, regridded, single-dish image. |
For interpolation to produce a faithful output, the grid of the input cannot be too coarse - the number of pixels across the single-dish beam should be at least 3. Although this is usually readily satisfied on the spatial axes, spectral line users should be more wary: the number of channels across the spectral axis response function is usually very small. In addition, whereas the spatial resolution of the single dish will be poorer than the mosaic, this is not necessarily true for the spectral axis. So, in addition to regridding, some spectral smearing might need to be performed to reduce the spectral resolution of the single-dish image to that of the mosaic. Unfortunately, Miriad is poorly equipped to solve these problems; these issues are left as a difficult exercise for the user.
A common approach to estimating flux calibration factor is to compare the single-dish and mosaic images in that annulus of spatial frequencies where both are sensitive. For example, a Parkes observation will sample spatial frequencies from 0 to near 64 meters, and a ATCA mosaic (assuming the shortest ATCA spacing is observed) will be sensitive from about 20 meters upwards. As the reliability of the data near the extremes is suspect, spatial frequencies from 25 to 40 meters should be modestly reliably represented in both an ATCA mosaic and a Parkes image. Comparing in the overlap annulus works very well when there is a simple source which is well represented in this region (e.g. a dominant point source).
The tasks to perform the linear and non-linear combination both have parameters for setting and deducing the flux calibration factor.
While it is possible to imagine algorithms that could deduce the beamwidth based on the spatial frequency overlap between the single-dish and mosaic data, in practise this is not possible (the overlap region is not wide nor is the data reliable enough). In estimating beamwidth parameters, the beamwidth and flux calibration factor are highly coupled; an error in the beamwidth has much the same effect in the overlap region as an error in the flux scale. Thus you cannot determine the beamwidth and flux scale simultaneously.
The linear combination method assumes the single-dish point-spread function is a gaussian form, whereas the non-linear method takes a image dataset as the point-spread function (the point-spread function need not be symmetric). If a point-spread function is not readily available for the non-linear method, a gaussian dataset can be produced.
Miriad manager