feather_compare¶
- uvcombine.feather_compare(hires, lores, SAS, LAS, lowresfwhm, highresextnum=0, lowresextnum=0, beam_divide_lores=True, lowpassfilterSD=False, min_beam_fraction=0.1, plot_min_beam_fraction=0.001, doplot=True, return_samples=False, weights=None)[source]¶
Compare the single-dish and interferometer data over the region where they should agree
- Parameters:
- highresfitsfilestr
The high-resolution FITS file
- lowresfitsfilestr
The low-resolution (single-dish) FITS file
- SAS
astropy.units.Quantity The smallest angular scale to plot
- LAS
astropy.units.Quantity The largest angular scale to plot (probably the LAS of the high resolution data)
- lowresfwhm
astropy.units.Quantity The full-width-half-max of the single-dish (low-resolution) beam; or the scale at which you want to try to match the low/high resolution data
- highresextnumint, optional
Select the HDU when passing a multi-HDU FITS file for the high-resolution data.
- lowresextnumint, optional
Select the HDU when passing a multi-HDU FITS file for the low-resolution data.
- beam_divide_lores: bool
Divide the low-resolution data by the beam weight before plotting? (should do this: otherwise, you are plotting beam-downweighted data)
- min_beam_fractionfloat
The minimum fraction of the beam to include; values below this fraction will be discarded when deconvolving
- plot_min_beam_fractionfloat
Like min_beam_fraction, but used only for plotting
- doplotbool
If true, make plots. Otherwise will just return the results.
- return_samplesbool, optional
Return the samples in the overlap region. This includes: the angular scale at each point, the ratio, the high-res values, and the low-res values.
- weights
ndarray, optional Provide an array of weights with the spatial shape of the high-res data. This is useful when either of the data have emission at the map edge, which will lead to ringing in the Fourier transform. A weights array can be provided to smoothly taper the edges of each map to avoid this issue.
- Returns:
- statsdict
Statistics on the ratio of the high-resolution FFT data to the low-resolution FFT data over the range SAS < x < LAS. Sigma-clipped stats are included.