[daip] Can't run IMAGR on data
Amy Mioduszewski
amiodusz at nrao.edu
Tue Dec 6 13:01:10 EST 2011
Hi Andy,
You have your cellsize set to 1 asec, since this is EVN data my guess is that
much too big. Set it smaller (0.01?) and will hopefully fix things. Sometimes
the error messages have nothing to do with the actual problem.
Amy
Andy Biggs wrote:
> 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
>
> _______________________________________________
> Daip mailing list
> Daip at listmgr.cv.nrao.edu
> http://listmgr.cv.nrao.edu/mailman/listinfo/daip
>
More information about the Daip
mailing list