Skip to main content


Showing posts from October, 2011

Querying with a Python script

I have written a Python script/module to query the spectral line catalog compilation. It does not rely on any SPLAT-VO protocol at all. It queries the database through the normal web interface, which makes it possible to query with the same parameters as through the web interface. Yay! The instructions below is for Ubuntu Linux, but is probably easily adapted to other platforms.

EDIT2:I have rewritten the module and the updated version is now part of astropy.astroquery ( However, the CLI interface is stripped out of this version. Furthermore, according to, the php interface will be replace by a java onlye sometime this year (, see Milestone no 4.) which would mean that I have to update it.  

The code is hosted on my trunk (bazaar versioning):

EDIT:I have moved my code to github (!)…

Interoperating GILDAS and MIRIAD, attempt 1

EDIT: There is a guide for interoperating GILDAS and MIRIAD here :

I am now trying to work with data from the Submillimeter Array (SMA), calibrated in IDL with the MIR package, continuum subtracted in MIRIAD, exported to UVFITS and read in to GILDAS for imaging. Yes, it is a mess.

It kind of works. Except the RMS noise somehow is 10-15% higher for the imaged data in GILDAS than for the same uv-data and same RMS channels in MIRIAD. Why? I do not know, but I have now resorted to work only in MIRIAD with the MIR calibrated data. I will conduct my normal analysis there, and output a FITS file of the cleaned data to work with my own scripts for creating various maps and spectra.

The best would be if I could write/use some Python modules for fitting Gaussians in the uv-plane as well, so I could finally break free of the dependence of MIRIAD/GILDAS/CASA/whatever package that is used at telescope X.