feather(highres, lowres, imagename='', sdfactor=1.0, effdishdiam=- 1.0, lowpassfiltersd=False)¶
Combine two images using their Fourier transforms
imagename (string=’’) - Name of output feathered image
highres (string) - Name of high resolution (interferometer) image
lowres (string) - Name of low resolution (single dish) image
sdfactor (double=1.0) - Scale factor to apply to Single Dish image
effdishdiam (double=-1.0) - New effective SingleDish diameter to use in m
lowpassfiltersd (bool=False) - Filter out the high spatial frequencies of the SD image
feather can be used as one method of combining single-dish and interferometric images after they have been separately made. The algorithm converts each image to the gridded visibility plane, combines them, and reconverts them into a combined image. Each image must include a well-defined beam shape (clean beam) in order for feathering to work well; these could be a ‘clean beam’ for interferometric images, and a ‘primary-beam’ for a single-dish image. The two images must have the same flux density normalization scale.
More information about feather can be found in the Image Combination section of the CASAdocs.
Feathering is a simple method for combining two images with different spatial resolution. The processing steps are:
Regrid the low-resolution image to a temporary copy matching the resolution of the high-resolution image.
Transform each image to the spatial-frequency plane (gridded).
Scale the low-resolution image (uv-grid) by the ratio of the volumes of the two ‘clean beams’ (high-res/low-res).
Add to this, the uv-grid of the high-resolution image, scaled by (1-\(\omega\) t) where ‘\(\omega\) t’ is the Fourier transform of the ‘clean beam’ defined in the low-resolution image.
Transform back to the image plane.
One commonly used option for feather is sdfactor, which can be used to adjust the flux scale of the SD image in the case that the fluxes don’t match where there is overlap between the spatial-frequencies in the SD image and the interferometer image.
NOTE: The tclean task allows another method of combining single-dish and interferometric data. The single-dish image can be used as a starting model for the interferometric image-reconstruction. If there is some overlap between the spatial-frequencies contained in the single-dish image and the interferometer sampling function, then such a starting model will help constrain the solutions on the short-baselines of the interferometric data.
Name of output feathered image. Default is none; example: imagename=’orion_combined.im’.
Name of high resolution (interferometer) image. Default is none; example: highres=’orion_vla.im’. This image is often a clean image obtained from synthesis observations.
Name of low resolution (single dish) image. Default is none; example: lowres=’orion_gbt.im’. This image is often a image from a single-dish observations or a clean image obtained from lower resolution synthesis observations.
Value by which to scale the Single Dish image. Default is 1.0. Basically modifying the flux scale of the SD image.
New effective SingleDish diameter to use in meters [m]. One can only reduce the dish effective dish diameter in feathering. Default is -1.0 which means leave as is.
If True the high spatial frequency in the SD image is rejected. Any data outside the maximum uv distance that the SD has illuminated is filtered out.
Creating a image called ‘M100_Feather_CO.image’ from an ALMA interferometric cube, ‘M100_combine_CO_cube.image.subim’, and a single dish ALMA total power image, ‘M100_TP_CO_cube.regrid.subim.depb’. The inputs have been appropriately cleaned, regridded, and cropped beforehand.
Creating an image called ‘feather.im’ by combining the cleaned, synthesis image, ‘synth.im’ and the SD image, ‘single_dish.im’ while increasing the flux scale of the SD image by setting sdfactor = 1.2.
feather(imagename ='feather.im', highres ='synth.im', lowres ='single_dish.im'sdfactor = 1.2)
No additional development details
- Parameter Details
Detailed descriptions of each function parameter
imagename (string='')- Name of output feathered imageDefault: noneExample: imagename=’orion_combined.im’
highres (string)- Name of high resolution (interferometer) imageDefault: noneExample: imagename=’orion_vla.im’
lowres (string)- Name of low resolution (single dish) imageDefault: noneExample: imagename=’orion_gbt.im’
sdfactor (double=1.0)- Value by which to scale the Single Dish image.Default: 1.0Basically modifying the flux scale of the SD image
effdishdiam (double=-1.0)- New effective SingleDish diameter to use in mDefault: -1.0 (leave as is)Obviously one can only reduce the disheffective dish diameter in feathering.
lowpassfiltersd (bool=False)- Filter out the high spatial frequencies of the SD imageDefault: FalseIf True the high spatial frequency in the SDimage is rejected.Any data outside the maximum uv distancethat the SD has illuminated is filteredout.