Is anyone familiar with mosacing techniques in AIPS, Miriad, or both? I have 21cm VLA spectral line data, taken mosaic-style, and am having a heck of a time trying to stitch the pointings together. Here's what I have done so far:
Reduced in the standard manner with AIPS up until the 'SPLIT' step
Wrote out each pointing with FITTP
Loaded into Miriad
In Miriad, 'invert', 'mossdi', and 'restor', mostly with the defaults for parameters
The best way I can describe the results is like a bad patchwork quilt...does anyone know of a good cookbook or guide for putting together pointings into a mosaic? I know AIPS has some tasks for this (VTESS/UTESS/LTESS) but I haven't had success with them. Thanks!
Katie
Monday, August 11, 2008
Mosaicing: AIPS/Miriad/whatever works!
Tuesday, July 22, 2008
How does CALIB make use of clean components?
A question from Sam Tun....
I've been working on some self-cal programs for the Owens Valley Solar Array, and I am having similar problems as these. I would want a self-caled map that gives me a good approximation to my CLEAN map fluxes, but no go on that. So, my question is, does anyone know what is done to the clean components to get them to return appropriate visibilities to feed into self-cal? I believe that in AIPS you feed the CLEAN map as a model into calib, but does anyone know the details of what goes on in there (or where to find them)? Gracias.
Total Flux in an Image Region?
A question from an anonymous reader...
But I have this issue: How to find the total flux in a given region? tvstat can make regions, but it shows only mean and rms. Is there any task that will show total flux under a selection?
Wednesday, July 16, 2008
How does DBCON deal with flagging?
Does anybody out there understand how DBCON flags your data? The help file is being less than transparent. DBCON seems to copy over all the FG tables from the first image you give it to the DBCON'd file, but does not appear to copy FG tables from the second image.
There is this note in the help file:
"Also, any CL, FG, TY, WX, IM, MC, PC, AT, CT, OB, or GC tables with version=1 will have their source numbers translated and appended to the end of the corresponding table (if any) from the first file."
So, let's say I have 13 flag tables for Image1 and 6 flag tables for Image2 (which happens a lot with the new crazy FG table creation scheme). Does this mean FG #1 for Image2 gets appended to FG tables 1-13 from Image1? And then in the futere, I will be apply a combination of FG#1 to Image1 and FG #13 to Image2? That's kinda dumb...
Wednesday, July 9, 2008
How to Know When You Can Self-Cal and What Solint To Use
Have you ever done a self-calibration run to find out the self cal is actually making your images worse, not better? Have you ever guessed at what SOLINT to use while self calibrating? I know I have!
People always say you should evaluate the signal-to-noise of your data before self-calibrating, but I never understood what this meant until today! There is a simple equation to find out if you can self-cal and if the SOLINT you are considering might be too short....
First image your data and clean it pretty good. Afterwards you can look in the image header and not the total cleaned flux. This is your 'Signal'.
Second of all, you want to calculate the noise in your data, per baseline per SOLINT. First, measure the rms noise in your image, in Jy/beam (sig_image). Next, calculate the number of baselines in your data (N_base where N_base = ((N_ant * N_ant-1) / 2) and N_ant is the number of operational antennae). Finally, figure out how much time-on-source went into making your image, in minutes (TOS). The noise of your data per baseline for a given SOLINT (in minutes) is then:
Noise = sig_image * sqrt(N_base) * sqrt(TOS/SOLINT)
Now, compare the 'Signal' with the 'Noise'. For 'P' self cal, you want the Signal to be at least 5 times greater than the Noise. If it's not, then increase your SOLINT. For A&P self cal, you probably want a signal-to-noise of 10-20, at least.
Note: If you are doing multi-facet imaging (at lower frequencies), you want to use the total flux in your data-- that is the sum from all facets, The image headers tell you this as 'CCTOTAL'.
Self Calibration and DBCON
Let's say you have two data sets that you want to DBCON together. Should you self calibrate them individually? or after DBCONing?
The tip I got is you should self-cal each data set individually as best you can, then DBCON. However, after DBCONing, there might be some small gain offsets between the two data sets. Do a final A&P self cal on the DBCONed data set with a really long SOLINT, so that you basically have one SN solution per data set. This will ensure the original data sets are as consistent with one another as possible.
Two More Ways to Identify Bad Data
Crystal also said a good way to identify RFI is to look at Stokes V, either in TVFLG or UVPLT or whatever. RFI is usually polarized and will pop out as unusually high baselines/channels/times.
Stokes V will usually mimic your amplitude structure (as a function of baseline length). So, for example, if your source has high amplitudes at short baselines, Stokes V should also be a bit higher on these short baselines (so don't get confused and flag them).
And another way to find that bad data-- Data weights. Plot your data in UVPLT with BPARM = 0 13. The weights for all the data should cluster-- if there are any points that are anomalously low or high, flag them! Crystal says she likes to use WIPER. Obviously, be more concerned about data with really high weights that really low weights-- because this data will affect your overall data set more, as it is up-weighted!
Subscribe to:
Posts (Atom)
If you have a question about AIPS that CANNOT be answered by reading the