Advanced imaging and self-calibration

If you have followed the previous sections you should have ended up with a calibrated data set for all your sources (all of them stored in separated SPLIT files). Although the imaging process can be perfectly done within AIPS, as explained previously, Difmap allows you to perform a more dynamic imaging/calibration of the data. Its approach usually produces a easier way to perform the last and more refined steps in the calibration of an EVN dataset (basically flagging, imaging with different weightings, and self-calibration). For this reason, in this section we will describe how to improve the final EVN images by working with the data in both AIPS and Difmap and which steps will produce a better calibration on the target source. Although the tutorial is focused for phase-referencing observations, it can also be applied for other techniques.

This guide assumes that you have a small knowledge of how Difmap works. Otherwise, we recommend to follow a more detailed tutorial, as the one done by Doug Gobeille at Brandeis University. An exhaustive manual is also provided by ESO in this link (PS file), with a detailed explanation of all the tasks available in Difmap.

In the previous tutorial we have produced a standard calibration of your data by running fringe on your strongest sources, bandpass calibrating the data, and extrapolating the solutions to your target source. However, sometimes the obtained solutions are not optimal and some stations can exhibit a systematic offset in amplitude due to a wrong calibration (we remind you that in the EVN we have a wide variety of stations, with different capabilities, sizes, and sensitivities).

Because of that, we can improve the calibration if we self-calibrate the data using one strong source that has been observed regularly during the whole observation. In phase-referencing sessions this one should be the phase calibrator. Given that this source should be strong enough, close to our target source, and observed every few minutes, we can get much better solutions self-calibrating it and then extrapolating these solutions again to our target source. In those cases where your target source is bright enough, we can directly conduct this self-calibration directly on it, but in that case note that the absolute position of your source is lost (no astrometric studies can be performed).

The mentioned procedure will be carried out as follows:

  • Load the phase-calibrator data into Difmap and inspect the data
  • Produce an initial image of the source
  • Self-calibrate the data based on this image
  • Repeat this loop until we obtain accurate solutions
  • Save and import the data again into AIPS
  • Create a SN table containing the self-calibration solutions and extrapolate them to the target source
  • Finally image the target source

Load and inspect the data

During the last steps of the previous tutorial we exported the SPLIT files to the hard disk with FITTP. Hence we have the phase-calibrator data (calibrated) as a UVFITS file. We can then open Difmap and load these data into it (contrary to AIPS, Difmap works with the UVFITS files on the fly):

obs filename-phasecal-split.uvfits
select PI
Be careful!
Difmap loads all the data into memory. This produces that if the data size is too large and your resources are low, all the data cannot be loaded. For such reason you have the select command to select a subset of your data. In the general case, the EVN data files are typically small enough to not produce issues.

select allows you to select which polarization is going to be loaded (I, Q, U, RR,...) and optionally the number of subbands to be loaded (otherwise all of them). By default Difmap only produces I stokes if both, RR and LL are available. We can avoid this behavior selecting select PI (meaning pseudo-I). In this case Difmap will compute I stokes and will keep the data (RR or LL) for those stations with only one polarization.

In Difmap there are several tasks to plot your data and interact with them. The most common ones are:

  • radplot: shows the amplitudes/phases versus uv-distance.
  • vplot: amplitudes/phases versus time per baseline and IF. Specially useful to flag RFI that is only present in some baselines/stations at particular times.
  • uvplot: uv-coverage.
  • projplot: same as radplot but considering the position angle, and you can rotate to see the data from different projected angles.
We note that if you press a number between 1 and 3 you obtain a plot of amplitudes, phases, or amplitudes and phases, respectively. Additionally you can flag data by clicking with the mouse (it flags the closest point) or by making a flagging area (clicking the C character and moving the mouse).

Imaging with Difmap

Once we have a dataset free of RFI


Correcting the target source