[daip] Can't run IMAGR on data
Andy Biggs
abiggs at eso.org
Tue Dec 6 10:54:04 EST 2011
Hi. I've just finished the flagging and calibration of an EVN dataset - everything went very smoothly. As soon as I try and image the data though I get the following error:
localh> IMAGR1: Task IMAGR (release of 31DEC11) begins
localh> IMAGR1: Doing no flagging this time
localh> IMAGR1: Create HDF259 .IMAGR . 1 (UV) on disk 1 cno 2
localh> IMAGR1: Beginning channel 1 through 16 with 8 IFs
localh> IMAGR1: IMACPY: NO DATA SELECTED
localh> IMAGR1: IMACPY: ERROR COPYING Input UVdata
localh> IMAGR1: TO UVdata work object
localh> IMAGR1: Blanking channel 1 image plane
localh> IMAGR1: OBRGET: KEYWORD CNO NOT FOUND
localh> IMAGR1: KEYWORD =CNO
localh> IMAGR1: OBRGET: PROBLEM WITH OBJECT: CLEAN field number 0001
localh> IMAGR1: IMABLK: ERROR BLANKING CLEAN field number 0001
localh> IMAGR1: NO DATA FOUND
localh> IMAGR1: Deleting UV work file:
localh> IMAGR1: Destroyed UV image file: catno= 2 disk= 1
localh> IMAGR1: Purports to die of UNNATURAL causes
localh> IMAGR1: localhos 31DEC11 TST: Cpu= 4.1 Real= 48 IO= 972
I get the same error whether I use the multi-source file with calibration applied, or a split file with the calibration already applied. I also have no problems with any other tasks - plotting in POSSM, VPLOT and WIPER works fine. I have also tried this on two computers, a Linux box and a MacbookPro, with the same results.
My inputs to IMAGR seem to be very boring:
AIPS 1: IMAGR: Wide field imaging/Clean task
AIPS 1: Adverbs Values Comments
AIPS 1: ----------------------------------------------------------------
AIPS 1: INNAME 'HDF259' Input UV data (name)
AIPS 1: INCLASS 'SPLIT' Input UV data (class)
AIPS 1: INSEQ 1 Input UV data (seq. #)
AIPS 1: INDISK 1 Input UV data disk drive #
AIPS 1: SOURCES 'HDF259' Source name
AIPS 1: *rest ' '
AIPS 1: QUAL -1 Calibrator qualifier -1=>all
AIPS 1: CALCODE ' ' Calibrator code ' '=>all
AIPS 1: TIMERANG *all 0 Time range to use
AIPS 1: SELBAND -1 Bandwidth to select (kHz)
AIPS 1: SELFREQ -1 Frequency to select (MHz)
AIPS 1: FREQID -1 Freq. ID to select.
AIPS 1: SUBARRAY 0 Sub-array, 0=>all
AIPS 1: ANTENNAS *all 0 Antennas to plot
AIPS 1: BASELINE *all 0 Baselines with ANTENNAS
AIPS 1: DOCALIB -1 > 0 calibrate data & weights
AIPS 1: > 99 do NOT calibrate weights
AIPS 1: GAINUSE 0 CL (or SN) table to apply
AIPS 1: DOPOL -1 If >0.5 correct polarization.
AIPS 1: PDVER 0 PD table to apply (DOPOL>0)
AIPS 1: BLVER -1 BL table to apply.
AIPS 1: FLAGVER 0 Flag table version
AIPS 1: DOBAND -1 If >0.5 apply bandpass cal.
AIPS 1: Method used depends on value
AIPS 1: of DOBAND (see HELP file).
AIPS 1: BPVER 0 Bandpass table version
AIPS 1: SMOOTH *all 0 Smoothing function. See
AIPS 1: HELP SMOOTH for details.
AIPS 1: STOKES 'I' Stokes parameters (see HELP)
AIPS 1: BCHAN 1 Low freq. channel 0 for cont.
AIPS 1: ECHAN 0 Highest freq channel
AIPS 1: CHANNEL 0 Restart channel number
AIPS 1: NCHAV 128 Number of chan. to average.
AIPS 1: CHINC 0 Channel incr. between maps.
AIPS 1: BIF 1 First IF in average.
AIPS 1: EIF 8 Last IF in average.
AIPS 1: OUTNAME ' ' Output image name (name)
AIPS 1: OUTDISK 1 Output image disk drive #
AIPS 1: OUTSEQ 0 Output seq. no.
AIPS 1: OUTVER 0 CC ver. no (Continuum only)
AIPS 1: *** SET OUTVER ON RESTARTS
AIPS 1: IN2NAME ' ' UV work file name
AIPS 1: IN2CLASS ' ' UV work file class
AIPS 1: IN2SEQ 0 UV work file seq
AIPS 1: *** SET TO KEEP WORK FILE
AIPS 1: IN2DISK 1 UV work file disk
AIPS 1: CELLSIZE 1 1 (X,Y) size of grid in asec
AIPS 1: IMSIZE 512 512 Minimum image size
AIPS 1: NFIELD 1 Number of fields (max 4096)
AIPS 1: DO3DIMAG -1 > 0 => use different tangent
AIPS 1: points for each field
AIPS 1: FLDSIZE *all 0 Clean size of each field.
AIPS 1: RASHIFT *all 0 RA shift per field (asec)
AIPS 1: DECSHIFT *all 0 DEC shift per field (asec)
AIPS 1: UVTAPER 0 0 (U,V) Gaussian taper
AIPS 1: units are kilo-lambda
AIPS 1: UVRANGE 0 0 Min & max baseline (klambda)
AIPS 1: GUARD 0 0 x,y guard band fractional
AIPS 1: radius
AIPS 1: ROTATE 0 Rotate image CCW from N by
AIPS 1: ROTATE degrees
AIPS 1: ZEROSP *all 0 0-spacing fluxes and weights
AIPS 1: SEE HELP!!
AIPS 1: UVWTFN 'NA' UV dist. weight function
AIPS 1: UVSIZE 0 0 Array size for doing uniform
AIPS 1: weights. 0 -> actual field
AIPS 1: size.
AIPS 1: ROBUST 0 Robustness power: -5 -> pure
AIPS 1: uniform weights, 5 => natural
AIPS 1: UVBOX 0 Additional rows and columns
AIPS 1: used in weighting.
AIPS 1: UVBXFN 1 Box function type when UVBOX
AIPS 1: > 0. 0 -> 1 round pill box.
AIPS 1: XTYPE 5 Conv. function type in x
AIPS 1: default spheroidal
AIPS 1: YTYPE 5 Conv. function type in y
AIPS 1: default spheroidal
AIPS 1: XPARM *all 0 Conv. function parms for x
AIPS 1: YPARM *all 0 Conv. function parms for y
AIPS 1: NITER 10 Maximum # of Clean components
AIPS 1: BCOMP *all 0 Begin at BCOMP component
AIPS 1: Specify for each field.
AIPS 1: ALLOKAY 0 For restart: > 0 => beams
AIPS 1: okay, > 1 => work file too
AIPS 1: NBOXES 0 Number of boxes for Clean
AIPS 1: NB: field 1 only.
AIPS 1: CLBOX *all 0 Four coordinates for each box
AIPS 1: BOXFILE *all ' ' Input file of field params
AIPS 1: and Clean boxes; ' ' => use
AIPS 1: FLDSIZE, RASHIFT, DECSHIFT,
AIPS 1: NBOXES, CLBOX only.
AIPS 1: OBOXFILE *all ' ' Output file for final Clean
AIPS 1: boxes
AIPS 1: GAIN 0.1 Clean loop gain
AIPS 1: FLUX 0 Minimum Clean component (Jy)
AIPS 1: MINPATCH 51 Min. BEAM half-width in AP.
AIPS 1: BMAJ 0 FWHM(asec) major axis Clean
AIPS 1: restoring beam.
AIPS 1: BMIN 0 FWHM(asec) minor axis Clean
AIPS 1: restoring beam.
AIPS 1: BPA 0 Clean beam position angle
AIPS 1: OVERLAP 0 1 => restore components to
AIPS 1: overlapped fields, >=2=>
AIPS 1: expect overlaps in Cleaning
AIPS 1: ONEBEAM -1 > 0 use only 1 dirty beam
AIPS 1: per scale in multi-facet
AIPS 1: Cleans
AIPS 1: OVRSWTCH 0 Not 0 => switch from OVERLAP
AIPS 1: >= 2 to OVERLAP 1 - see HELP
AIPS 1: PHAT 0 Prussian hat height.
AIPS 1: FACTOR 0 Speedup factor see HELP
AIPS 1: CMETHOD ' ' Modeling method:
AIPS 1: 'DFT','GRID',' '
AIPS 1: IMAGRPRM *all 0 Task enrichment parameters
AIPS 1: (1) Antenna diameter (m)
AIPS 1: (2) Source Spectral index
AIPS 1: (3) Frequency scaling factor
AIPS 1: (4) > 0 -> SDI Clean factor
AIPS 1: (5) >0 => scale residuals
AIPS 1: (6) Half-width in x of box
AIPS 1: (7) Half-width in y of box
AIPS 1: (8) Filter components whose
AIPS 1: neighborhood is weaker than
AIPS 1: IMAGRPRM(8) Jy. 0 -> don't
AIPS 1: (9) Radius in pixels for the
AIPS 1: IMAGRPRM(8) test.
AIPS 1: (10) multiplier of image size
AIPS 1: to get beam size: 0 => 2;
AIPS 1: 2, 1, 0.5 0.25 supported
AIPS 1: (11-16) Multi-scale controls
AIPS 1: (17) spectral index radius
AIPS 1: 0 -> no correction
AIPS 1: (18) Limit grids (see help)
AIPS 1: (19) Dynamic range limit
AIPS 1: (20) Retry factor (see help)
AIPS 1: IM2PARM *all 0 Yet more parameters:
AIPS 1: (1) Auto boxes: allowed #
AIPS 1: (2) : island level
AIPS 1: (3) : peak required
AIPS 1: (4) : limit wrt max
AIPS 1: (5) : extend boxes
AIPS 1: (6) : edge skip
AIPS 1: (7) reset boxes for next chan
AIPS 1: (11) baseline-dependent avg
AIPS 1: max time in sec
AIPS 1: (12) field size 0 -> infinite
AIPS 1: NGAUSS 0 Number of scales to use
AIPS 1: WGAUSS *all 0 Scales in arc sec >= 0
AIPS 1: FGAUSS *all 0 Minimum flux for each resol.
AIPS 1: MAXPIXEL 0 Maximum pixels searched in
AIPS 1: each major cycle.
AIPS 1: IN3NAME ' ' Spectral index image name
AIPS 1: IN3CLASS ' ' Spectral index image class
AIPS 1: IN3SEQ 0 Spectral index image sequence
AIPS 1: number
AIPS 1: IN3DISK 0 Spectral index image disk
AIPS 1: IN4NAME ' ' Spectral curvature name
AIPS 1: IN4CLASS ' ' Spectral curvature class
AIPS 1: IN4SEQ 0 Spectral curvature sequence
AIPS 1: number
AIPS 1: IN4DISK 0 Spectral curvature disk
AIPS 1: FQTOL -1 Frequency tolerance in kHz
AIPS 1: (primary beam & spec index)
AIPS 1: DOTV 1 Display residuals on TV ?
AIPS 1: Start with field = DOTV
AIPS 1: GRCHAN 0 Graphics channel of boundary
AIPS 1: BADDISK *all 0 Disks to avoid for scratch.
And here is the file header:
AIPS 1: Image=HDF259 (UV) Filename=HDF259 .SPLIT . 1
AIPS 1: Telescope=EVN Receiver=VLBA
AIPS 1: Observer=EC033 User #= 5000
AIPS 1: Observ. date=13-JUN-2011 Map date=06-DEC-2011
AIPS 1: # visibilities 326254 Sort order TB
AIPS 1: Rand axes: UU-L-SIN VV-L-SIN WW-L-SIN TIME1 BASELINE
AIPS 1: INTTIM CORR-ID
AIPS 1: ----------------------------------------------------------------
AIPS 1: Type Pixels Coord value at Pixel Coord incr Rotat
AIPS 1: COMPLEX 3 0.0000000E+00 1.00 1.0000000E+00 0.00
AIPS 1: STOKES 2 -1.0000000E+00 1.00 -1.0000000E+00 0.00
AIPS 1: FREQ 16 1.6265212E+09 0.53 5.0000000E+05 0.00
AIPS 1: IF 8 1.0000000E+00 1.00 1.0000000E+00 0.00
AIPS 1: RA 1 12 37 13.870 1.00 3600.000 0.00
AIPS 1: DEC 1 62 18 26.500 1.00 3600.000 0.00
AIPS 1: ----------------------------------------------------------------
AIPS 1: Coordinate equinox 2000.00
AIPS 1: Rest freq 0.000 Vel type: OPTICAL wrt YOU
AIPS 1: Alt ref. value 0.00000E+00 wrt pixel 0.53
AIPS 1: Maximum version number of extension files of type HI is 1
AIPS 1: Maximum version number of extension files of type NX is 1
AIPS 1: Maximum version number of extension files of type FQ is 1
AIPS 1: Maximum version number of extension files of type AN is 1
AIPS 1: Keyword = 'OLDRFQ ' value = 1.62652125D+09
As usual, I am running a pretty up-to-date binary installation of 31DEC11.
Thanks.
Andy
More information about the Daip
mailing list