Notebook created by Martín Villanueva - martin.villanueva@usm.cl - LIRAE/UTFSM - October 2017.
In [1]:
import numpy as np
import matplotlib.pyplot as mpl
%matplotlib inline
# just for presentation purposes
import warnings
warnings.filterwarnings("ignore")
In [5]:
%%bash
pip install pycupid
Collecting pycupid
Using cached pycupid-0.1.4-cp36-cp36m-macosx_10_6_intel.macosx_10_9_intel.macosx_10_9_x86_64.macosx_10_10_intel.macosx_10_10_x86_64.whl
Requirement already satisfied: numpy>=1.11.2 in /Users/martin/anaconda3/lib/python3.6/site-packages (from pycupid)
Installing collected packages: pycupid
Successfully installed pycupid-0.1.4
In [2]:
from astropy.io import fits
hdulist = fits.open('FITS/ALMA00000034.fits')
hdulist.info()
print('\nTarget Object: '+hdulist[0].header['OBJECT'])
Filename: FITS/ALMA00000034.fits
No. Name Ver Type Cards Dimensions Format
0 PRIMARY 1 PrimaryHDU 731 (100, 100, 1, 1) float32
Target Object: Orion KL
In [3]:
orion = hdulist[0]
tmp = orion.copy()
orion.header
Out[3]:
SIMPLE = T /Standard FITS
BITPIX = -32 /Floating point (32 bit)
NAXIS = 4
NAXIS1 = 100
NAXIS2 = 100
NAXIS3 = 1
NAXIS4 = 1
BSCALE = 1.000000000000E+00 /PHYSICAL = PIXEL*BSCALE + BZERO
BZERO = 0.000000000000E+00
BMAJ = 5.160116155942E-04
BMIN = 3.801464703348E-04
BPA = 1.758169860840E+02
BTYPE = 'Intensity'
OBJECT = 'Orion KL'
BUNIT = 'JY/BEAM ' /Brightness (pixel) unit
EQUINOX = 2.000000000000E+03
RADESYS = 'FK5 '
LONPOLE = 1.800000000000E+02
LATPOLE = -5.376388888889E+00
PC001001= 1.000000000000E+00
PC002001= 0.000000000000E+00
PC003001= 0.000000000000E+00
PC004001= 0.000000000000E+00
PC001002= 0.000000000000E+00
PC002002= 1.000000000000E+00
PC003002= 0.000000000000E+00
PC004002= 0.000000000000E+00
PC001003= 0.000000000000E+00
PC002003= 0.000000000000E+00
PC003003= 1.000000000000E+00
PC004003= 0.000000000000E+00
PC001004= 0.000000000000E+00
PC002004= 0.000000000000E+00
PC003004= 0.000000000000E+00
PC004004= 1.000000000000E+00
CTYPE1 = 'RA---SIN'
CRVAL1 = 8.380979166667E+01
CDELT1 = -1.111111111111E-04
CRPIX1 = 5.100000000000E+01
CUNIT1 = 'deg '
CTYPE2 = 'DEC--SIN'
CRVAL2 = -5.376388888889E+00
CDELT2 = 1.111111111111E-04
CRPIX2 = 5.100000000000E+01
CUNIT2 = 'deg '
CTYPE3 = 'FREQ '
CRVAL3 = 2.309497185811E+11
CDELT3 = 1.513671875000E+07
CRPIX3 = 1.000000000000E+00
CUNIT3 = 'Hz '
CTYPE4 = 'STOKES '
CRVAL4 = 1.000000000000E+00
CDELT4 = 1.000000000000E+00
CRPIX4 = 1.000000000000E+00
CUNIT4 = ' '
PV2_1 = 0.000000000000E+00
PV2_2 = 0.000000000000E+00
RESTFRQ = 2.309497185811E+11 /Rest Frequency (Hz)
SPECSYS = 'TOPOCENT' /Spectral reference frame
ALTRVAL = -0.000000000000E+00 /Alternate frequency reference value
ALTRPIX = 1.000000000000E+00 /Alternate frequency reference pixel
VELREF = 259 /1 LSR, 2 HEL, 3 OBS, +256 Radio
COMMENT casacore non-standard usage: 4 LSD, 5 GEO, 6 SOU, 7 GAL
TELESCOP= 'ALMA '
OBSERVER= 'srandall'
DATE-OBS= '2012-01-20T01:30:59.808001'
TIMESYS = 'UTC '
OBSRA = 8.380979166667E+01
OBSDEC = -5.376388888889E+00
OBSGEO-X= 2.225142180269E+06
OBSGEO-Y= -5.440307370349E+06
OBSGEO-Z= -2.481029851874E+06
DATE = '2012-04-06T06:43:52.749000' /Date FITS file was written
ORIGIN = 'CASA 3.4.0 (stable r19148)'
HISTORY CASA START LOGTABLE
HISTORY 2012-02-17T14:47:02 INFO SRCCODE='imager::clean()'
HISTORY 2012-02-17T14:47:02 INFO SRCCODE='imager::clean()'
HISTORY Sun Jan 22 01:51:31 2012 HISTORY flagdata [] taskname = flagdata
HISTORY Sun Jan 22 01:51:32 2012 HISTORY flagdata [] vis = "uid___A002_X37
HISTORY >b127_X2f.ms"
HISTORY Sun Jan 22 01:51:32 2012 HISTORY flagdata [] mode = "manualflag"
HISTORY Sun Jan 22 01:59:49 2012 HISTORY flagdata [] taskname = flagdata
HISTORY Sun Jan 22 01:59:49 2012 HISTORY flagdata [] vis = "uid___A002_X37
HISTORY >b127_X2f.ms"
HISTORY Sun Jan 22 01:59:49 2012 HISTORY flagdata [] mode = "manualflag"
HISTORY Sun Jan 22 02:01:39 2012 HISTORY flagdata [] taskname = flagdata
HISTORY Sun Jan 22 02:01:39 2012 HISTORY flagdata [] vis = "uid___A002_X37
HISTORY >b127_X2f.ms"
HISTORY Sun Jan 22 02:01:39 2012 HISTORY flagdata [] mode = "manualflag"
HISTORY Sun Jan 22 02:12:33 2012 HISTORY calibrater::setdata [Beginning selectv
HISTORY >is--(MSSelection version)-------, chanmode=none nchan=1 start=0 step=1
HISTORY >Start='0km/s' mStep='0km/s' msSelect='']
HISTORY Sun Jan 22 02:18:28 2012 HISTORY calibrater::correct [] Beginning corre
HISTORY >ct---------------------------
HISTORY Sun Jan 22 02:25:20 2012 HISTORY applycal [] taskname=applycal
HISTORY Sun Jan 22 02:25:20 2012 HISTORY applycal [] version: 3.3.0 3.3.0 rev.
HISTORY > 16856 Thu 2011/11/03 18:24:40 UTC
HISTORY Sun Jan 22 02:25:20 2012 HISTORY applycal [] vis = "uid___A002_
HISTORY >X37b127_X2f.ms"
HISTORY Sun Jan 22 02:25:20 2012 HISTORY applycal [] field = "1,2,3"
HISTORY Sun Jan 22 02:25:20 2012 HISTORY applycal [] spw = "33,35,37,39
HISTORY >,65,67,69,71,73,75,77,79,81,83,85,87,89,91,93,95"
HISTORY Sun Jan 22 02:25:20 2012 HISTORY applycal [] intent = ""
HISTORY Sun Jan 22 02:25:20 2012 HISTORY applycal [] selectdata = True
HISTORY Sun Jan 22 02:25:20 2012 HISTORY applycal [] timerange = ""
HISTORY Sun Jan 22 02:25:20 2012 HISTORY applycal [] uvrange = ""
HISTORY Sun Jan 22 02:25:20 2012 HISTORY applycal [] antenna = ""
HISTORY Sun Jan 22 02:25:20 2012 HISTORY applycal [] scan = ""
HISTORY Sun Jan 22 02:25:20 2012 HISTORY applycal [] observation = ""
HISTORY Sun Jan 22 02:25:20 2012 HISTORY applycal [] msselect = ""
HISTORY Sun Jan 22 02:25:20 2012 HISTORY applycal [] gaintable = ['uid___A002
HISTORY >_X37b127_X2f.ms.wvr.smooth']
HISTORY Sun Jan 22 02:25:20 2012 HISTORY applycal [] gainfield = ['1,2,3']
HISTORY Sun Jan 22 02:25:20 2012 HISTORY applycal [] interp = ['linear']
HISTORY Sun Jan 22 02:25:20 2012 HISTORY applycal [] spwmap = []
HISTORY Sun Jan 22 02:25:20 2012 HISTORY applycal [] gaincurve = False
HISTORY Sun Jan 22 02:25:20 2012 HISTORY applycal [] opacity = []
HISTORY Sun Jan 22 02:25:20 2012 HISTORY applycal [] parang = False
HISTORY Sun Jan 22 02:25:20 2012 HISTORY applycal [] calwt = False
HISTORY Sun Jan 22 02:25:20 2012 HISTORY applycal [] flagbackup = False
HISTORY Sun Jan 22 23:37:55 2012 HISTORY ms::split() [fieldids=2,3 spwids=33,35
HISTORY >,37,39,65,67,69,71,73,75,77,79,81,83,85,87,89,91,93,95 step=[1] which='
HISTORY >orrected'] dataset1 split from /data1/home/randall/sciver/uid___A002_X3
HISTORY >b127_X2f.ms
HISTORY Sun Jan 22 23:37:55 2012 HISTORY split [] taskname=split
HISTORY Sun Jan 22 23:37:55 2012 HISTORY split [] version: 3.3.0 3.3.0 rev. 168
HISTORY >56 Thu 2011/11/03 18:24:40 UTC
HISTORY Sun Jan 22 23:37:55 2012 HISTORY split [] vis = "uid___A002_X37
HISTORY >b127_X2f.ms"
HISTORY Sun Jan 22 23:37:55 2012 HISTORY split [] outputvis = "dataset1"
HISTORY Sun Jan 22 23:37:55 2012 HISTORY split [] datacolumn = "corrected"
HISTORY Sun Jan 22 23:37:55 2012 HISTORY split [] field = "2,3"
HISTORY Sun Jan 22 23:37:55 2012 HISTORY split [] spw = "33,35,37,39,65
HISTORY >,67,69,71,73,75,77,79,81,83,85,87,89,91,93,95"
HISTORY Sun Jan 22 23:37:55 2012 HISTORY split [] width = 1
HISTORY Sun Jan 22 23:37:55 2012 HISTORY split [] antenna = ""
HISTORY Sun Jan 22 23:37:55 2012 HISTORY split [] timebin = "0.0s"
HISTORY Sun Jan 22 23:37:55 2012 HISTORY split [] timerange = ""
HISTORY Sun Jan 22 23:37:55 2012 HISTORY split [] scan = ""
HISTORY Sun Jan 22 23:37:55 2012 HISTORY split [] intent = ""
HISTORY Sun Jan 22 23:37:55 2012 HISTORY split [] array = ""
HISTORY Sun Jan 22 23:37:55 2012 HISTORY split [] uvrange = ""
HISTORY Sun Jan 22 23:37:55 2012 HISTORY split [] correlation = ""
HISTORY Sun Jan 22 23:37:55 2012 HISTORY split [] observation = ""
HISTORY Sun Jan 22 23:37:55 2012 HISTORY split [] combine = ""
HISTORY Sun Jan 22 23:37:55 2012 HISTORY split [] keepflags = True
HISTORY Mon Jan 23 00:36:56 2012 HISTORY flagdata [] taskname = flagdata
HISTORY Mon Jan 23 00:36:56 2012 HISTORY flagdata [] vis = "dataset1"
HISTORY Mon Jan 23 00:36:56 2012 HISTORY flagdata [] mode = "shadow"
HISTORY Mon Jan 23 00:48:32 2012 HISTORY calibrater::setdata [Beginning selectv
HISTORY >is--(MSSelection version)-------, chanmode=none nchan=1 start=0 step=1
HISTORY >Start='0km/s' mStep='0km/s' msSelect='']
HISTORY Mon Jan 23 00:48:33 2012 HISTORY calibrater::solve [] Beginning solve--
HISTORY >---------------------------
HISTORY Mon Jan 23 00:49:11 2012 HISTORY calibrater::setdata [Beginning selectv
HISTORY >is--(MSSelection version)-------, chanmode=none nchan=1 start=0 step=1
HISTORY >Start='0km/s' mStep='0km/s' msSelect='']
HISTORY Mon Jan 23 00:49:12 2012 HISTORY calibrater::solve [] Beginning solve--
HISTORY >---------------------------
HISTORY Mon Jan 23 01:17:52 2012 HISTORY calibrater::setdata [Beginning selectv
HISTORY >is--(MSSelection version)-------, chanmode=none nchan=1 start=0 step=1
HISTORY >Start='0km/s' mStep='0km/s' msSelect='']
HISTORY Mon Jan 23 01:17:53 2012 HISTORY calibrater::solve [] Beginning solve--
HISTORY >---------------------------
HISTORY Mon Jan 23 01:22:00 2012 HISTORY calibrater::setdata [Beginning selectv
HISTORY >is--(MSSelection version)-------, chanmode=none nchan=1 start=0 step=1
HISTORY >Start='0km/s' mStep='0km/s' msSelect='']
HISTORY Mon Jan 23 01:22:01 2012 HISTORY calibrater::solve [] Beginning solve--
HISTORY >---------------------------
HISTORY Mon Jan 23 02:07:39 2012 HISTORY setjy [] taskname=setjy
HISTORY Mon Jan 23 02:07:39 2012 HISTORY setjy [] version: 3.3.0 3.3.0 rev. 168
HISTORY >56 Thu 2011/11/03 18:24:40 UTC
HISTORY Mon Jan 23 02:07:39 2012 HISTORY setjy [] vis = "dataset1"
HISTORY Mon Jan 23 02:07:39 2012 HISTORY setjy [] field = "0"
HISTORY Mon Jan 23 02:07:39 2012 HISTORY setjy [] spw = ""
HISTORY Mon Jan 23 02:07:39 2012 HISTORY setjy [] selectdata = False
HISTORY Mon Jan 23 02:07:39 2012 HISTORY setjy [] timerange = ""
HISTORY Mon Jan 23 02:07:39 2012 HISTORY setjy [] scan = ""
HISTORY Mon Jan 23 02:07:39 2012 HISTORY setjy [] observation = ""
HISTORY Mon Jan 23 02:07:39 2012 HISTORY setjy [] modimage = ""
HISTORY Mon Jan 23 02:07:39 2012 HISTORY setjy [] listmodimages = False
HISTORY Mon Jan 23 02:07:39 2012 HISTORY setjy [] scalebychan = False
HISTORY Mon Jan 23 02:07:39 2012 HISTORY setjy [] fluxdensity = [1.5, 0, 0, 0]
HISTORY Mon Jan 23 02:07:39 2012 HISTORY setjy [] spix = 0
HISTORY Mon Jan 23 02:07:39 2012 HISTORY setjy [] reffreq = "1GHz"
HISTORY Mon Jan 23 02:07:39 2012 HISTORY setjy [] standard = "Perley-Butler
HISTORY > 2010"
HISTORY Mon Jan 23 02:09:15 2012 HISTORY imager::data selection [] Selected 163
HISTORY >2 out of 114240 rows.
HISTORY Mon Jan 23 02:18:44 2012 HISTORY calibrater::setdata [Beginning selectv
HISTORY >is--(MSSelection version)-------, chanmode=none nchan=1 start=0 step=1
HISTORY >Start='0km/s' mStep='0km/s' msSelect='']
HISTORY Mon Jan 23 02:18:45 2012 HISTORY calibrater::solve [] Beginning solve--
HISTORY >---------------------------
HISTORY Mon Jan 23 02:31:09 2012 HISTORY calibrater::setdata [Beginning selectv
HISTORY >is--(MSSelection version)-------, chanmode=none nchan=1 start=0 step=1
HISTORY >Start='0km/s' mStep='0km/s' msSelect='']
HISTORY Mon Jan 23 02:41:23 2012 HISTORY calibrater::fluxscale() [] Table g3_co
HISTORY >rr is of type: T Jones
HISTORY Mon Jan 23 02:57:53 2012 HISTORY calibrater::setdata [Beginning selectv
HISTORY >is--(MSSelection version)-------, chanmode=none nchan=1 start=0 step=1
HISTORY >Start='0km/s' mStep='0km/s' msSelect='']
HISTORY Mon Jan 23 02:57:53 2012 HISTORY calibrater::solve [] Beginning solve--
HISTORY >---------------------------
HISTORY Mon Jan 23 02:59:16 2012 HISTORY calibrater::setdata [Beginning selectv
HISTORY >is--(MSSelection version)-------, chanmode=none nchan=1 start=0 step=1
HISTORY >Start='0km/s' mStep='0km/s' msSelect='']
HISTORY Mon Jan 23 02:59:17 2012 HISTORY calibrater::solve [] Beginning solve--
HISTORY >---------------------------
HISTORY Mon Jan 23 03:10:37 2012 HISTORY calibrater::setdata [Beginning selectv
HISTORY >is--(MSSelection version)-------, chanmode=none nchan=1 start=0 step=1
HISTORY >Start='0km/s' mStep='0km/s' msSelect='']
HISTORY Mon Jan 23 03:10:38 2012 HISTORY calibrater::solve [] Beginning solve--
HISTORY >---------------------------
HISTORY Mon Jan 23 03:11:58 2012 HISTORY calibrater::setdata [Beginning selectv
HISTORY >is--(MSSelection version)-------, chanmode=none nchan=1 start=0 step=1
HISTORY >Start='0km/s' mStep='0km/s' msSelect='']
HISTORY Mon Jan 23 03:11:59 2012 HISTORY calibrater::correct [] Beginning corre
HISTORY >ct---------------------------
HISTORY Mon Jan 23 03:13:44 2012 HISTORY applycal [] taskname=applycal
HISTORY Mon Jan 23 03:13:44 2012 HISTORY applycal [] version: 3.3.0 3.3.0 rev.
HISTORY > 16856 Thu 2011/11/03 18:24:40 UTC
HISTORY Mon Jan 23 03:13:44 2012 HISTORY applycal [] vis = "dataset1"
HISTORY Mon Jan 23 03:13:44 2012 HISTORY applycal [] field = "0"
HISTORY Mon Jan 23 03:13:44 2012 HISTORY applycal [] spw = ""
HISTORY Mon Jan 23 03:13:44 2012 HISTORY applycal [] intent = ""
HISTORY Mon Jan 23 03:13:44 2012 HISTORY applycal [] selectdata = True
HISTORY Mon Jan 23 03:13:44 2012 HISTORY applycal [] timerange = ""
HISTORY Mon Jan 23 03:13:44 2012 HISTORY applycal [] uvrange = ""
HISTORY Mon Jan 23 03:13:44 2012 HISTORY applycal [] antenna = ""
HISTORY Mon Jan 23 03:13:44 2012 HISTORY applycal [] scan = ""
HISTORY Mon Jan 23 03:13:44 2012 HISTORY applycal [] observation = ""
HISTORY Mon Jan 23 03:13:44 2012 HISTORY applycal [] msselect = ""
HISTORY Mon Jan 23 03:13:44 2012 HISTORY applycal [] gaintable = ['b1_corr_sm
HISTORY >ooth', 'g2_corr', 'g3_corr']
HISTORY Mon Jan 23 03:13:44 2012 HISTORY applycal [] gainfield = ['']
HISTORY Mon Jan 23 03:13:44 2012 HISTORY applycal [] interp = ['linear']
HISTORY Mon Jan 23 03:13:44 2012 HISTORY applycal [] spwmap = []
HISTORY Mon Jan 23 03:13:44 2012 HISTORY applycal [] gaincurve = False
HISTORY Mon Jan 23 03:13:44 2012 HISTORY applycal [] opacity = []
HISTORY Mon Jan 23 03:13:44 2012 HISTORY applycal [] parang = False
HISTORY Mon Jan 23 03:13:44 2012 HISTORY applycal [] calwt = False
HISTORY Mon Jan 23 03:13:44 2012 HISTORY applycal [] flagbackup = False
HISTORY Mon Jan 23 03:22:59 2012 HISTORY calibrater::setdata [Beginning selectv
HISTORY >is--(MSSelection version)-------, chanmode=none nchan=1 start=0 step=1
HISTORY >Start='0km/s' mStep='0km/s' msSelect='']
HISTORY Mon Jan 23 03:23:00 2012 HISTORY calibrater::correct [] Beginning corre
HISTORY >ct---------------------------
HISTORY Mon Jan 23 03:27:44 2012 HISTORY applycal [] taskname=applycal
HISTORY Mon Jan 23 03:27:44 2012 HISTORY applycal [] version: 3.3.0 3.3.0 rev.
HISTORY > 16856 Thu 2011/11/03 18:24:40 UTC
HISTORY Mon Jan 23 03:27:44 2012 HISTORY applycal [] vis = "dataset1"
HISTORY Mon Jan 23 03:27:44 2012 HISTORY applycal [] field = "1"
HISTORY Mon Jan 23 03:27:44 2012 HISTORY applycal [] spw = ""
HISTORY Mon Jan 23 03:27:44 2012 HISTORY applycal [] intent = ""
HISTORY Mon Jan 23 03:27:44 2012 HISTORY applycal [] selectdata = True
HISTORY Mon Jan 23 03:27:44 2012 HISTORY applycal [] timerange = ""
HISTORY Mon Jan 23 03:27:44 2012 HISTORY applycal [] uvrange = ""
HISTORY Mon Jan 23 03:27:44 2012 HISTORY applycal [] antenna = ""
HISTORY Mon Jan 23 03:27:44 2012 HISTORY applycal [] scan = ""
HISTORY Mon Jan 23 03:27:44 2012 HISTORY applycal [] observation = ""
HISTORY Mon Jan 23 03:27:44 2012 HISTORY applycal [] msselect = ""
HISTORY Mon Jan 23 03:27:44 2012 HISTORY applycal [] gaintable = ['b1_corr',
HISTORY > 'g3_corr', 'g5_corr']
HISTORY Mon Jan 23 03:27:44 2012 HISTORY applycal [] gainfield = ['', '0', '0
HISTORY >']
HISTORY Mon Jan 23 03:27:44 2012 HISTORY applycal [] interp = ['linear']
HISTORY Mon Jan 23 03:27:44 2012 HISTORY applycal [] spwmap = []
HISTORY Mon Jan 23 03:27:44 2012 HISTORY applycal [] gaincurve = False
HISTORY Mon Jan 23 03:27:44 2012 HISTORY applycal [] opacity = []
HISTORY Mon Jan 23 03:27:44 2012 HISTORY applycal [] parang = False
HISTORY Mon Jan 23 03:27:44 2012 HISTORY applycal [] calwt = False
HISTORY Mon Jan 23 03:27:44 2012 HISTORY applycal [] flagbackup = False
HISTORY Mon Jan 23 03:28:10 2012 HISTORY calibrater::setdata [Beginning selectv
HISTORY >is--(MSSelection version)-------, chanmode=none nchan=1 start=0 step=1
HISTORY >Start='0km/s' mStep='0km/s' msSelect='']
HISTORY Mon Jan 23 03:28:11 2012 HISTORY calibrater::correct [] Beginning corre
HISTORY >ct---------------------------
HISTORY Mon Jan 23 03:30:14 2012 HISTORY applycal [] taskname=applycal
HISTORY Mon Jan 23 03:30:14 2012 HISTORY applycal [] version: 3.3.0 3.3.0 rev.
HISTORY > 16856 Thu 2011/11/03 18:24:40 UTC
HISTORY Mon Jan 23 03:30:14 2012 HISTORY applycal [] vis = "dataset1"
HISTORY Mon Jan 23 03:30:14 2012 HISTORY applycal [] field = "1"
HISTORY Mon Jan 23 03:30:14 2012 HISTORY applycal [] spw = ""
HISTORY Mon Jan 23 03:30:14 2012 HISTORY applycal [] intent = ""
HISTORY Mon Jan 23 03:30:14 2012 HISTORY applycal [] selectdata = True
HISTORY Mon Jan 23 03:30:14 2012 HISTORY applycal [] timerange = ""
HISTORY Mon Jan 23 03:30:14 2012 HISTORY applycal [] uvrange = ""
HISTORY Mon Jan 23 03:30:14 2012 HISTORY applycal [] antenna = ""
HISTORY Mon Jan 23 03:30:14 2012 HISTORY applycal [] scan = ""
HISTORY Mon Jan 23 03:30:14 2012 HISTORY applycal [] observation = ""
HISTORY Mon Jan 23 03:30:14 2012 HISTORY applycal [] msselect = ""
HISTORY Mon Jan 23 03:30:14 2012 HISTORY applycal [] gaintable = ['b1_corr_sm
HISTORY >ooth', 'g3_corr', 'g5_corr']
HISTORY Mon Jan 23 03:30:14 2012 HISTORY applycal [] gainfield = ['', '0', '0
HISTORY >']
HISTORY Mon Jan 23 03:30:14 2012 HISTORY applycal [] interp = ['linear']
HISTORY Mon Jan 23 03:30:14 2012 HISTORY applycal [] spwmap = []
HISTORY Mon Jan 23 03:30:14 2012 HISTORY applycal [] gaincurve = False
HISTORY Mon Jan 23 03:30:14 2012 HISTORY applycal [] opacity = []
HISTORY Mon Jan 23 03:30:14 2012 HISTORY applycal [] parang = False
HISTORY Mon Jan 23 03:30:14 2012 HISTORY applycal [] calwt = False
HISTORY Mon Jan 23 03:30:14 2012 HISTORY applycal [] flagbackup = False
HISTORY Wed Feb 15 18:25:46 2012 HISTORY setjy [] taskname=setjy
HISTORY Wed Feb 15 18:25:46 2012 HISTORY setjy [] version: 3.3.0 3.3.0 rev. 168
HISTORY >56 Thu 2011/11/03 18:24:40 UTC
HISTORY Wed Feb 15 18:25:46 2012 HISTORY setjy [] vis = "dataset1"
HISTORY Wed Feb 15 18:25:46 2012 HISTORY setjy [] field = "0"
HISTORY Wed Feb 15 18:25:46 2012 HISTORY setjy [] spw = ""
HISTORY Wed Feb 15 18:25:46 2012 HISTORY setjy [] selectdata = False
HISTORY Wed Feb 15 18:25:46 2012 HISTORY setjy [] timerange = ""
HISTORY Wed Feb 15 18:25:46 2012 HISTORY setjy [] scan = ""
HISTORY Wed Feb 15 18:25:46 2012 HISTORY setjy [] observation = ""
HISTORY Wed Feb 15 18:25:46 2012 HISTORY setjy [] modimage = ""
HISTORY Wed Feb 15 18:25:46 2012 HISTORY setjy [] listmodimages = False
HISTORY Wed Feb 15 18:25:46 2012 HISTORY setjy [] scalebychan = False
HISTORY Wed Feb 15 18:25:47 2012 HISTORY setjy [] fluxdensity = [1.399999999999
HISTORY >9999, 0, 0, 0]
HISTORY Wed Feb 15 18:25:47 2012 HISTORY setjy [] spix = 0
HISTORY Wed Feb 15 18:25:47 2012 HISTORY setjy [] reffreq = "1GHz"
HISTORY Wed Feb 15 18:25:47 2012 HISTORY setjy [] standard = "Perley-Butler
HISTORY > 2010"
HISTORY Wed Feb 15 18:27:25 2012 HISTORY imager::data selection [] Selected 163
HISTORY >2 out of 114240 rows.
HISTORY Wed Feb 15 18:35:01 2012 HISTORY flagdata [] taskname = flagdata
HISTORY Wed Feb 15 18:35:01 2012 HISTORY flagdata [] vis = "dataset1"
HISTORY Wed Feb 15 18:35:01 2012 HISTORY flagdata [] mode = "shadow"
HISTORY Wed Feb 15 18:37:43 2012 HISTORY flagdata [] taskname = flagdata
HISTORY Wed Feb 15 18:37:43 2012 HISTORY flagdata [] vis = "dataset1"
HISTORY Wed Feb 15 18:37:43 2012 HISTORY flagdata [] mode = "manualflag"
HISTORY Wed Feb 15 18:39:16 2012 HISTORY calibrater::setdata [Beginning selectv
HISTORY >is--(MSSelection version)-------, chanmode=none nchan=1 start=0 step=1
HISTORY >Start='0km/s' mStep='0km/s' msSelect='']
HISTORY Wed Feb 15 18:39:18 2012 HISTORY calibrater::solve [] Beginning solve--
HISTORY >---------------------------
HISTORY Wed Feb 15 18:40:35 2012 HISTORY calibrater::setdata [Beginning selectv
HISTORY >is--(MSSelection version)-------, chanmode=none nchan=1 start=0 step=1
HISTORY >Start='0km/s' mStep='0km/s' msSelect='']
HISTORY Wed Feb 15 18:40:36 2012 HISTORY calibrater::solve [] Beginning solve--
HISTORY >---------------------------
HISTORY Wed Feb 15 18:42:12 2012 HISTORY calibrater::setdata [Beginning selectv
HISTORY >is--(MSSelection version)-------, chanmode=none nchan=1 start=0 step=1
HISTORY >Start='0km/s' mStep='0km/s' msSelect='']
HISTORY Wed Feb 15 18:42:13 2012 HISTORY calibrater::solve [] Beginning solve--
HISTORY >---------------------------
HISTORY Wed Feb 15 18:50:10 2012 HISTORY calibrater::setdata [Beginning selectv
HISTORY >is--(MSSelection version)-------, chanmode=none nchan=1 start=0 step=1
HISTORY >Start='0km/s' mStep='0km/s' msSelect='']
HISTORY Wed Feb 15 18:50:12 2012 HISTORY calibrater::correct [] Beginning corre
HISTORY >ct---------------------------
HISTORY Wed Feb 15 18:51:36 2012 HISTORY applycal [] taskname=applycal
HISTORY Wed Feb 15 18:51:36 2012 HISTORY applycal [] version: 3.3.0 3.3.0 rev.
HISTORY > 16856 Thu 2011/11/03 18:24:40 UTC
HISTORY Wed Feb 15 18:51:36 2012 HISTORY applycal [] vis = "dataset1"
HISTORY Wed Feb 15 18:51:36 2012 HISTORY applycal [] field = "0"
HISTORY Wed Feb 15 18:51:36 2012 HISTORY applycal [] spw = ""
HISTORY Wed Feb 15 18:51:36 2012 HISTORY applycal [] intent = ""
HISTORY Wed Feb 15 18:51:36 2012 HISTORY applycal [] selectdata = True
HISTORY Wed Feb 15 18:51:36 2012 HISTORY applycal [] timerange = ""
HISTORY Wed Feb 15 18:51:36 2012 HISTORY applycal [] uvrange = ""
HISTORY Wed Feb 15 18:51:36 2012 HISTORY applycal [] antenna = ""
HISTORY Wed Feb 15 18:51:36 2012 HISTORY applycal [] scan = ""
HISTORY Wed Feb 15 18:51:37 2012 HISTORY applycal [] observation = ""
HISTORY Wed Feb 15 18:51:37 2012 HISTORY applycal [] msselect = ""
HISTORY Wed Feb 15 18:51:37 2012 HISTORY applycal [] gaintable = ['b1_corr_sm
HISTORY >ooth', 'g2_corr', 'g3_corr']
HISTORY Wed Feb 15 18:51:37 2012 HISTORY applycal [] gainfield = ['']
HISTORY Wed Feb 15 18:51:37 2012 HISTORY applycal [] interp = ['linear']
HISTORY Wed Feb 15 18:51:37 2012 HISTORY applycal [] spwmap = []
HISTORY Wed Feb 15 18:51:37 2012 HISTORY applycal [] gaincurve = False
HISTORY Wed Feb 15 18:51:37 2012 HISTORY applycal [] opacity = []
HISTORY Wed Feb 15 18:51:37 2012 HISTORY applycal [] parang = False
HISTORY Wed Feb 15 18:51:37 2012 HISTORY applycal [] calwt = False
HISTORY Wed Feb 15 18:51:37 2012 HISTORY applycal [] flagbackup = False
HISTORY Wed Feb 15 18:57:39 2012 HISTORY calibrater::setdata [Beginning selectv
HISTORY >is--(MSSelection version)-------, chanmode=none nchan=1 start=0 step=1
HISTORY >Start='0km/s' mStep='0km/s' msSelect='']
HISTORY Wed Feb 15 18:57:41 2012 HISTORY calibrater::correct [] Beginning corre
HISTORY >ct---------------------------
HISTORY Wed Feb 15 18:58:44 2012 HISTORY calibrater::setdata [Beginning selectv
HISTORY >is--(MSSelection version)-------, chanmode=none nchan=1 start=0 step=1
HISTORY >Start='0km/s' mStep='0km/s' msSelect='']
HISTORY Wed Feb 15 18:58:45 2012 HISTORY calibrater::correct [] Beginning corre
HISTORY >ct---------------------------
HISTORY Wed Feb 15 19:02:23 2012 HISTORY applycal [] taskname=applycal
HISTORY Wed Feb 15 19:02:23 2012 HISTORY applycal [] version: 3.3.0 3.3.0 rev.
HISTORY > 16856 Thu 2011/11/03 18:24:40 UTC
HISTORY Wed Feb 15 19:02:23 2012 HISTORY applycal [] vis = "dataset1"
HISTORY Wed Feb 15 19:02:23 2012 HISTORY applycal [] field = "1"
HISTORY Wed Feb 15 19:02:23 2012 HISTORY applycal [] spw = ""
HISTORY Wed Feb 15 19:02:23 2012 HISTORY applycal [] intent = ""
HISTORY Wed Feb 15 19:02:23 2012 HISTORY applycal [] selectdata = True
HISTORY Wed Feb 15 19:02:23 2012 HISTORY applycal [] timerange = ""
HISTORY Wed Feb 15 19:02:23 2012 HISTORY applycal [] uvrange = ""
HISTORY Wed Feb 15 19:02:23 2012 HISTORY applycal [] antenna = ""
HISTORY Wed Feb 15 19:02:23 2012 HISTORY applycal [] scan = ""
HISTORY Wed Feb 15 19:02:23 2012 HISTORY applycal [] observation = ""
HISTORY Wed Feb 15 19:02:23 2012 HISTORY applycal [] msselect = ""
HISTORY Wed Feb 15 19:02:23 2012 HISTORY applycal [] gaintable = ['b1_corr_sm
HISTORY >ooth', 'g3_corr', 'g5_corr']
HISTORY Wed Feb 15 19:02:23 2012 HISTORY applycal [] gainfield = ['', '0', '0
HISTORY >']
HISTORY Wed Feb 15 19:02:23 2012 HISTORY applycal [] interp = ['linear']
HISTORY Wed Feb 15 19:02:23 2012 HISTORY applycal [] spwmap = []
HISTORY Wed Feb 15 19:02:23 2012 HISTORY applycal [] gaincurve = False
HISTORY Wed Feb 15 19:02:23 2012 HISTORY applycal [] opacity = []
HISTORY Wed Feb 15 19:02:23 2012 HISTORY applycal [] parang = False
HISTORY Wed Feb 15 19:02:23 2012 HISTORY applycal [] calwt = False
HISTORY Wed Feb 15 19:02:23 2012 HISTORY applycal [] flagbackup = False
HISTORY Thu Feb 16 13:05:52 2012 HISTORY concat [] taskname=concat
HISTORY Thu Feb 16 13:05:52 2012 HISTORY concat [] vis = "comb_data"
HISTORY Thu Feb 16 13:05:52 2012 HISTORY concat [] concatvis = "dataset2"
HISTORY Thu Feb 16 13:05:52 2012 HISTORY concat [] freqtol = "500kHz"
HISTORY Thu Feb 16 13:05:52 2012 HISTORY concat [] dirtol = ""
HISTORY Thu Feb 16 13:15:26 2012 HISTORY concat [] taskname=concat
HISTORY Thu Feb 16 13:15:26 2012 HISTORY concat [] vis = "comb_data"
HISTORY Thu Feb 16 13:15:26 2012 HISTORY concat [] concatvis = "dataset3"
HISTORY Thu Feb 16 13:15:27 2012 HISTORY concat [] freqtol = "500kHz"
HISTORY Thu Feb 16 13:15:27 2012 HISTORY concat [] dirtol = ""
HISTORY Thu Feb 16 13:47:23 2012 HISTORY ms::split() [fieldids= spwids=* step=[
HISTORY >1] which='corrected'] comb_corr split from /mnt/sco-red/data1/home/sran
HISTORY >all/sciver/comb_data
HISTORY Thu Feb 16 13:47:23 2012 HISTORY split [] taskname=split
HISTORY Thu Feb 16 13:47:23 2012 HISTORY split [] version: 3.3.0 3.3.0 rev. 168
HISTORY >56 Thu 2011/11/03 18:24:40 UTC
HISTORY Thu Feb 16 13:47:23 2012 HISTORY split [] vis = "comb_data"
HISTORY Thu Feb 16 13:47:23 2012 HISTORY split [] outputvis = "comb_corr"
HISTORY Thu Feb 16 13:47:23 2012 HISTORY split [] datacolumn = "corrected"
HISTORY Thu Feb 16 13:47:23 2012 HISTORY split [] field = ""
HISTORY Thu Feb 16 13:47:23 2012 HISTORY split [] spw = ""
HISTORY Thu Feb 16 13:47:23 2012 HISTORY split [] width = 1
HISTORY Thu Feb 16 13:47:23 2012 HISTORY split [] antenna = ""
HISTORY Thu Feb 16 13:47:23 2012 HISTORY split [] timebin = "0.0s"
HISTORY Thu Feb 16 13:47:23 2012 HISTORY split [] timerange = ""
HISTORY Thu Feb 16 13:47:23 2012 HISTORY split [] scan = ""
HISTORY Thu Feb 16 13:47:23 2012 HISTORY split [] intent = ""
HISTORY Thu Feb 16 13:47:23 2012 HISTORY split [] array = ""
HISTORY Thu Feb 16 13:47:23 2012 HISTORY split [] uvrange = ""
HISTORY Thu Feb 16 13:47:23 2012 HISTORY split [] correlation = ""
HISTORY Thu Feb 16 13:47:24 2012 HISTORY split [] observation = ""
HISTORY Thu Feb 16 13:47:24 2012 HISTORY split [] combine = ""
HISTORY Thu Feb 16 13:47:24 2012 HISTORY split [] keepflags = True
HISTORY Fri Feb 17 14:28:09 2012 HISTORY ms::split() [fieldids=1 spwids=4:1850~
HISTORY >1880 step=[1] which='data'] cont231 split from /mnt/sco-red/data1/home/
HISTORY >randall/sciver/comb_corr
HISTORY Fri Feb 17 14:28:09 2012 HISTORY split [] taskname=split
HISTORY Fri Feb 17 14:28:09 2012 HISTORY split [] version: 3.3.0 3.3.0 rev. 168
HISTORY >56 Thu 2011/11/03 18:24:40 UTC
HISTORY Fri Feb 17 14:28:09 2012 HISTORY split [] vis = "comb_corr"
HISTORY Fri Feb 17 14:28:09 2012 HISTORY split [] outputvis = "cont231"
HISTORY Fri Feb 17 14:28:09 2012 HISTORY split [] datacolumn = "data"
HISTORY Fri Feb 17 14:28:09 2012 HISTORY split [] field = "1"
HISTORY Fri Feb 17 14:28:09 2012 HISTORY split [] spw = "4:1850~1880"
HISTORY Fri Feb 17 14:28:09 2012 HISTORY split [] width = 1
HISTORY Fri Feb 17 14:28:09 2012 HISTORY split [] antenna = ""
HISTORY Fri Feb 17 14:28:09 2012 HISTORY split [] timebin = "0.0s"
HISTORY Fri Feb 17 14:28:09 2012 HISTORY split [] timerange = ""
HISTORY Fri Feb 17 14:28:09 2012 HISTORY split [] scan = ""
HISTORY Fri Feb 17 14:28:09 2012 HISTORY split [] intent = ""
HISTORY Fri Feb 17 14:28:09 2012 HISTORY split [] array = ""
HISTORY Fri Feb 17 14:28:09 2012 HISTORY split [] uvrange = ""
HISTORY Fri Feb 17 14:28:09 2012 HISTORY split [] correlation = ""
HISTORY Fri Feb 17 14:28:09 2012 HISTORY split [] observation = ""
HISTORY Fri Feb 17 14:28:09 2012 HISTORY split [] combine = ""
HISTORY Fri Feb 17 14:28:09 2012 HISTORY split [] keepflags = True
HISTORY Fri Feb 17 14:29:38 2012 HISTORY imager::data selection [mode=none ncha
HISTORY >n=[-1] start=[0] step=[1] mstart='Radialvelocity: 0' mstep='Radialveloc
HISTORY >ty: 0' spectralwindowids=[] fieldids=[] msselect=]
HISTORY Fri Feb 17 14:29:38 2012 HISTORY imager::defineimage() [Defining image
HISTORY > properties:nx=100 ny=100 cellx='0.4arcsec' celly='0.4arcsec' stokes=I'
HISTORY > mode=MFS nchan=-1 start=0 step=1 spwids=[-1] fieldid=0 facets=1 frame=
HISTORY > distance='0', phaseCenter='field-0 ' mStart='Radialvelocity: 0' qStep
HISTORY >'0 '' mFreqStart='Frequency: 0]
HISTORY Fri Feb 17 14:29:38 2012 HISTORY imager::weight() [] Weighting MS: Imag
HISTORY >ing weights will be changed
HISTORY Fri Feb 17 14:29:38 2012 HISTORY imager::weight() [] Natural weighting
HISTORY Fri Feb 17 14:29:39 2012 HISTORY imager::clean() [] Using multifield Ho
HISTORY >gbom clean
HISTORY Fri Feb 17 14:29:39 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 0, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:29:39 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:29:40 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:29:40 2012 HISTORY imager::clean() [] Fitted beam used in
HISTORY > restoration: 1.85764 by 1.36853 (arcsec) at pa 175.817 (deg)
HISTORY Fri Feb 17 14:30:11 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 100, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:30:11 2012 HISTORY imager::clean() [] Starting deconvolut
HISTORY >ion
HISTORY Fri Feb 17 14:30:14 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:30:14 2012 HISTORY imager::clean() [] Beam used in restor
HISTORY >ation: 1.85764 by 1.36853 (arcsec) at pa 175.817 (deg)
HISTORY Fri Feb 17 14:30:20 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 100, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:30:20 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:30:21 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:30:41 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 100, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:30:41 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:30:42 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:31:29 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 100, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:31:29 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:31:31 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:31:43 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 100, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:31:43 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:31:45 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:32:23 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 100, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:32:23 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:32:24 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:32:45 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 100, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:32:45 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:32:48 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:33:00 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 100, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:33:00 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:33:01 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:33:08 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 100, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:33:08 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:33:09 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:33:29 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 100, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:33:29 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:33:32 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:33:38 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 100, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:33:38 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:33:39 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:33:43 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 100, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:33:43 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:33:44 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:33:49 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 100, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:33:49 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:33:50 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:33:57 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 500, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:33:57 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:33:57 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:34:01 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 500, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:34:01 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:34:02 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:34:15 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 1000, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:34:15 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:34:16 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:35:52 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 1000, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:35:52 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:35:55 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:36:15 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 1000, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:36:15 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:36:16 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:36:37 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 1000, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:36:37 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:36:38 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:37:26 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 1000, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:37:26 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:37:30 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:37:52 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 1000, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:37:52 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:37:53 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:38:42 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 1000, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:38:42 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:38:48 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:39:11 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 1000, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:39:11 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:39:12 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:39:55 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 1000, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:39:55 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:40:00 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:40:10 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 1000, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:40:10 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:40:11 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:40:57 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 1000, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:40:57 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:41:03 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:41:21 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 1000, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:41:21 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:41:22 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:41:32 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 1000, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:41:32 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:41:33 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:42:37 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 1000, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:42:37 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:42:42 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:43:04 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 1000, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:43:04 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:43:08 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:43:19 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 1000, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:43:19 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:43:20 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:43:26 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 1000, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:43:26 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:43:26 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:44:03 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 1000, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:44:03 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:44:06 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:44:13 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 1000, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:44:13 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:44:14 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:44:17 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 1000, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:44:17 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:44:18 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:44:24 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 1000, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:44:24 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:44:25 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:46:26 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 1000, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:46:26 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:46:29 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:46:36 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 1000, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:46:36 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:46:37 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:46:51 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 1000, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:46:51 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:46:52 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY Fri Feb 17 14:47:01 2012 HISTORY imager::clean() [] Clean gain = 0.1, N
HISTORY >iter = 1000, Threshold = 1.5 mJy
HISTORY Fri Feb 17 14:47:01 2012 HISTORY imager::clean() [] Continuing deconvol
HISTORY >ution
HISTORY Fri Feb 17 14:47:02 2012 HISTORY imager::clean() [] Threshhold not reac
HISTORY >hed yet.
HISTORY CASA END LOGTABLE
FRAMEID = 'ALMA00000034_00_00'
In [4]:
import aplpy
fig = aplpy.FITSFigure(orion)
fig.show_colorscale(cmap='gist_heat')
fig.add_colorbar()
fig.show_grid()
fig.show_contour(colors='white')
INFO: Setting slices=[0, 0] [aplpy.core]
INFO: Auto-setting vmin to -1.677e-01 [aplpy.core]
INFO: Auto-setting vmax to 1.093e+00 [aplpy.core]
In [5]:
# CUPID Clumping algorithms
from pycupid import clumpfind, fellwalker, gaussclumps, reinhold
from ipywidgets import interact, fixed, IntSlider, FloatSlider
In [6]:
# data as a NumPy array (Don't run this twice)
data = orion.data.sum(axis=(0,1))
data = np.flip(data,0)
In [7]:
def run_cupid(data, algorithm=fellwalker, rms=0.1, config=None):
# CUPID clumping algorithm
out = algorithm(data, rms, config)
out[out==out.min()] = 0.
# Plotting the result with APLpy
fig = mpl.figure(figsize=(15, 7))
# fig1
fig1 = aplpy.FITSFigure(orion, figure=fig, subplot=[0.1,0.1,0.35,0.8], auto_refresh=True)
fig1.show_colorscale(cmap='gist_heat')
fig1.show_grid()
fig1.show_contour(colors='white')
# fig2
tmp.data[0,0,:,:] = np.flip(out.T, 0)
fig2 = aplpy.FITSFigure(tmp, figure=fig, subplot=[0.53,0.1,0.35,0.8], auto_refresh=True)
fig2.show_colorscale(cmap='magma')
fig2.show_grid()
fig.canvas.draw()
mpl.show()
In [8]:
algorithms = {'gaussclumps':gaussclumps, 'clumpfind':clumpfind,
'fellwalker':fellwalker, 'reinhold':reinhold}
interact(run_cupid,
data = fixed(data),
algorithm = algorithms,
rms = FloatSlider(min=data.min(), max=data.max()*0.2, step=0.01, value=0.05),
config = fixed(None)
);
In [9]:
def run_fellwalker(data, RMS=0.05, CLEANITER=None, FLATSLOPE=None, MAXBAD=None, MINDIP=None,
MINHEIGHT=None, MINPIX=None, MAXJUMP=None):
# CUPID clumping algorithm
config = dict()
if CLEANITER is not None: config['CLEANITER']=CLEANITER
if FLATSLOPE is not None: config['FLATSLOPE']=FLATSLOPE
if MAXBAD is not None: config['MAXBAD']=MAXBAD
if MINDIP is not None: config['MINDIP']=MINDIP
if MINHEIGHT is not None: config['MINHEIFHT']=MINHEIGHT
if MINPIX is not None: config['MINPIX']=MINPIX
if MAXJUMP is not None: config['MAXJUMP']=MAXJUMP
out = fellwalker(data, RMS, config)
out[out==out.min()] = 0.
# Plotting the result with APLpy
fig = mpl.figure(figsize=(15, 7))
# fig1
fig1 = aplpy.FITSFigure(orion, figure=fig, subplot=[0.1,0.1,0.35,0.8], auto_refresh=True)
fig1.show_colorscale(cmap='gist_heat')
fig1.show_grid()
fig1.show_contour(colors='white')
# fig2
tmp.data[0,0,:,:] = np.flip(out.T, 0)
fig2 = aplpy.FITSFigure(tmp, figure=fig, subplot=[0.53,0.1,0.35,0.8], auto_refresh=True)
fig2.show_colorscale(cmap='magma')
fig2.show_grid()
fig.canvas.draw()
mpl.show()
In [11]:
interact(run_fellwalker,
data = fixed(data),
RMS = FloatSlider(min=data.min(), max=data.max()*0.2, step=0.01, value=0.05),
CLEANITER = IntSlider(min=1, max=10, step=1, value=1),
FLATSLOPE = FloatSlider(min=0, max=data.max()*0.25, step=0.01, value=0.05),
MAXBAD = FloatSlider(min=0., max=1., step=0.01, value=0.05),
MINDIP = FloatSlider(min=0., max=data.max()*0.25, step=0.01, value=0.05),
MINHEIGHT = FloatSlider(min=0., max=data.max()*0.25, step=0.01, value=0.05),
MINPIX = IntSlider(min=0, max=20, step=1, value=7),
MAXJUMP = IntSlider(min=1, max=10, step=1, value=4)
);
In [ ]:
Content source: ChileanVirtualObservatory/pycupid
Similar notebooks: