[daip] COMB OPCODE="POLC" not working.
R. Craig Walker
cwalker at nrao.edu
Tue Mar 19 22:28:04 EDT 2013
It appears that COMB with OPCODE='POLC' produces identical output to
what it gives with OPCODE='POLI'. The history indicates it should be
different, but the images don't show the corrections. I asked Rick
Perley if he had been using COMB with POLC and he immediately said
it doesn't work, confirming what I see. He had not reported it yet.
While preparing this message, I also had a problem with IMEAN in version
NEW claiming to have a file open error, but otherwise working - reported
below.
Here are the history records for two runs of COMB with POLI and
POLC.
3680 COMB RELEASE='31DEC13 '
3681 COMB INNAME='3C274_98A_43' INCLASS='QCL001'
3682 COMB INSEQ= 76 INDISK= 5
3683 COMB IN2NAME='3C274_98A_43' IN2CLASS='UCL001'
3684 COMB IN2SEQ= 76 IN2DISK= 5
3685 COMB OUTNAME='3C274_98A_43' OUTCLASS='POLI '
3686 COMB OUTSEQ= 76 OUTDISK= 5
3687 COMB USERID= 604
3688 COMB CTYPE='POLI' /Pol.int
3689 COMB BLC= 1 1 1 1 1 1 1 / Bottom left corner
3690 COMB TRC= 1024 1024 1 1 1 1 1 / Top right corner
3691 COMB A(1)= 1.0000E+00 A(2)= 0.0000E+00
3692 COMB / Undefined pixels magic-value blanked
3680 COMB RELEASE='31DEC13 '
3681 COMB INNAME='3C274_98A_43' INCLASS='QCL001'
3682 COMB INSEQ= 76 INDISK= 5
3683 COMB IN2NAME='3C274_98A_43' IN2CLASS='UCL001'
3684 COMB IN2SEQ= 76 IN2DISK= 5
3685 COMB OUTNAME='3C274_98A_43' OUTCLASS='POLC '
3686 COMB OUTSEQ= 76 OUTDISK= 5
3687 COMB USERID= 604
3688 COMB CTYPE='POLC' /Poli Cor
3689 COMB BLC= 1 1 1 1 1 1 1 / Bottom left corner
3690 COMB TRC= 1024 1024 1 1 1 1 1 / Top right corner
3691 COMB A(1)= 1.0000E+00 A(2)= 0.0000E+00
3692 COMB / Undefined pixels magic-value blanked
3693 COMB CBMULT(3)= 0.00000 / scaled to match image 1 units
3694 COMB CBMULT(4)= 0.00000 / scaled to match image 1 units
3695 COMB SIGMA= 8.5034E-05 /Input image 1 noise level
3696 COMB SIGMA= 7.9251E-05 /Input image 2 noise level
Here are the headers, including statistics from IMEAN:
POLI image:
AIPS 1: Image=3C274 (MA) Filename=3C274_98A_43.POLI . 76
AIPS 1: Telescope=VLBA Receiver=VLBA
AIPS 1: Observer=BW098A User #= 604
AIPS 1: Observ. date=12-JAN-2013 Map date=17-MAR-2013
AIPS 1: Minimum= 1.18302594E-07 Maximum= 3.35429818E-03 JY/BEAM
AIPS 1: ----------------------------------------------------------------
AIPS 1: Type Pixels Coord value at Pixel Coord incr Rotat
AIPS 1: RA---SIN 1024 12 30 49.423 312.00 -0.000035 0.00
AIPS 1: DEC--SIN 1024 12 23 28.044 513.00 0.000035 0.00
AIPS 1: FREQ 1 4.3151852E+10 1.00 1.9200000E+08 0.00
AIPS 1: STOKES 1 5.0000000E+00 1.00 1.0000000E+00 0.00
AIPS 1: ----------------------------------------------------------------
AIPS 1: Coordinate equinox 2000.00
AIPS 1: Map type=NORMAL Number of iterations= 888
AIPS 1: Conv size= 0.00045 X 0.00020 Position angle= -15.48
AIPS 1: Phase shifted in X -0.007 in Y 0.000
AIPS 1: Rest freq 0.000 Vel type: OPTICAL wrt YOU
AIPS 1: Alt ref. value 0.00000E+00 wrt pixel 0.00
AIPS 1: Maximum version number of extension files of type HI is 1
AIPS 1: Maximum version number of extension files of type CG is 1
AIPS 1: Keyword = 'CCFLUX ' value = 9.865344E-03
AIPS 1: Keyword = 'CCTOTAL ' value = 9.865341E-03
AIPS 1: Keyword = 'OLDRFQ ' value = 4.30560000D+10
AIPS 1: Keyword = 'MAXABSU ' value = -1.000000E+00
AIPS 1: Keyword = 'SOURNAM1' value = '3C274 '
AIPS 1: Keyword = 'SOURNAM2' value = ' '
AIPS 1: Keyword = 'ACTNOISE' value = 6.026580E-05
AIPS 1: Keyword = 'ACTMEAN ' value = 8.833273E-05
POLC image
AIPS 1: Image=3C274 (MA) Filename=3C274_98A_43.POLC . 76
AIPS 1: Telescope=VLBA Receiver=VLBA
AIPS 1: Observer=BW098A User #= 604
AIPS 1: Observ. date=12-JAN-2013 Map date=17-MAR-2013
AIPS 1: Minimum= 1.18302594E-07 Maximum= 3.35429818E-03 JY/BEAM
AIPS 1: ----------------------------------------------------------------
AIPS 1: Type Pixels Coord value at Pixel Coord incr Rotat
AIPS 1: RA---SIN 1024 12 30 49.423 312.00 -0.000035 0.00
AIPS 1: DEC--SIN 1024 12 23 28.044 513.00 0.000035 0.00
AIPS 1: FREQ 1 4.3151852E+10 1.00 1.9200000E+08 0.00
AIPS 1: STOKES 1 5.0000000E+00 1.00 1.0000000E+00 0.00
AIPS 1: ----------------------------------------------------------------
AIPS 1: Coordinate equinox 2000.00
AIPS 1: Map type=NORMAL Number of iterations= 888
AIPS 1: Conv size= 0.00045 X 0.00020 Position angle= -15.48
AIPS 1: Phase shifted in X -0.007 in Y 0.000
AIPS 1: Rest freq 0.000 Vel type: OPTICAL wrt YOU
AIPS 1: Alt ref. value 0.00000E+00 wrt pixel 0.00
AIPS 1: Maximum version number of extension files of type HI is 1
AIPS 1: Maximum version number of extension files of type CG is 1
AIPS 1: Keyword = 'CCFLUX ' value = 9.865344E-03
AIPS 1: Keyword = 'CCTOTAL ' value = 9.865341E-03
AIPS 1: Keyword = 'OLDRFQ ' value = 4.30560000D+10
AIPS 1: Keyword = 'MAXABSU ' value = -1.000000E+00
AIPS 1: Keyword = 'SOURNAM1' value = '3C274 '
AIPS 1: Keyword = 'SOURNAM2' value = ' '
AIPS 1: Keyword = 'ACTNOISE' value = 6.026580E-05
AIPS 1: Keyword = 'ACTMEAN ' value = 8.833273E-05
Note the Minimum, Maximum, ACTNOISE, and ACTMEAN are identical
to many significant figures. When I load them to the TV I see
no effect when one replaces the other. Also note that the Minimum
is positive for the POLC which it should not be if the Ricean
bias were being removed.
I have data processed in February 2012 where POLC appeared to
be working.
Here are my inputs:
AIPS 1: COMB: Task to combine in many ways two overlapping images
AIPS 1: Adverbs Values Comments
AIPS 1: ----------------------------------------------------------------
AIPS 1: INNAME '3C274_98A_43' First image name
AIPS 1: INCLASS 'QCL001' First image class
AIPS 1: INSEQ 76 First image seq. #
AIPS 1: INDISK 5 First image disk drive #
AIPS 1: IN2NAME '3C274_98A_43' Second image name
AIPS 1: IN2CLASS 'UCL001' Second image class
AIPS 1: IN2SEQ 76 Second image seq. #
AIPS 1: IN2DISK 5 Second image disk drive #
AIPS 1: IN3NAME ' ' First noise image name
AIPS 1: IN3CLASS ' ' First noise image class
AIPS 1: IN3SEQ 0 First noise image seq. #
AIPS 1: IN3DISK 0 First noise image disk #
AIPS 1: IN4NAME ' ' Second noise image name
AIPS 1: IN4CLASS ' ' Second noise image class
AIPS 1: IN4SEQ 0 Second noise image seq. #
AIPS 1: IN4DISK 0 Second noiseimage disk #
AIPS 1: DOALIGN 1 Should images be coincident?
AIPS 1: (See HELP.)
AIPS 1: OUTNAME '3C274_98A_43' Output image name
AIPS 1: OUTCLASS 'POLCT' Output image class
AIPS 1: OUTSEQ 76 Output image seq. #
AIPS 1: OUTDISK 5 Output image disk drive #
AIPS 1: BLC *all 0 Bottom left corner
AIPS 1: TRC *all 0 Top right corner
AIPS 1: OPCODE 'POLC' Algorithm type:
AIPS 1: 'SUM ','DIV ','SPIX','POLI',
AIPS 1: 'POLA','MULT','OPTD','CLIP'
AIPS 1: 'REAL','IMAG','MEAN','RM '
AIPS 1: 'POLC','SUMM'
AIPS 1: APARM *all 0 Parameters for algorithm:
AIPS 1: (1) - (4) scale and offset
AIPS 1: (8) > 0 => blank with 0.0
AIPS 1: (9) Map1 clip level
AIPS 1: (10) Map2 clip level
AIPS 1: see HELP COMB
AIPS 1: BPARM 8.503E-05 7.925E-05 Noise/control parameters:
AIPS 1: *rest 0 (1) Map1 noise level
AIPS 1: -1 use 3rd image
AIPS 1: 0 ignore noise
AIPS 1: (2) Map2 noise level
AIPS 1: -1 use 4th image
AIPS 1: 0 ignore noise
AIPS 1: (3) > 0 => output noise
AIPS 1: in addition
AIPS 1: (4) < 0.5 => clip w inputs
AIPS 1: > 1.5 => clip w S/N
AIPS 1: else => clip w noise
AIPS 1: (5) minimum ok abs(S/N) or
AIPS 1: maximum ok noise
AIPS 1: (6) max output noise
AIPS 1: 0 -> any
AIPS 1: see HELP COMB
AIPS 1: DOHIST -1 -2 => copy 1st HI only
AIPS 1: -3 => copy no HI files
The run above actually produced a new file (outcl='polct') using
version 'new' (31DEC12).
SECONDARY PROBLEM REPORT:
When I ran IMEAN in NEW on that file, I get:
>go imean
AIPS 1: Found in Version=NEW
AIPS 1: Waiting for returned adverbs
IMEAN1: Task IMEAN (release of 31DEC12) begins
IMEAN1: Initial guess for PIXSTD taken from ACTNOISE in header
IMEAN1: Guess is Mean= 8.8333E-05 Rms= 6.0266E-05
IMEAN1: returns adverbs to AIPS
IMEAN1: ZTXOP2: using translated file name =
IMEAN1: ZTXOP2:
IMEAN1: ZTXOP2: FORTRAN OPEN ERROR = 10
IMEAN1: ZERROR: ON FILE
IMEAN1: ZERROR: IN ZTXOP2 ERRNO = 10 (No child processes)
IMEAN1: ERROR ON OUTPUT FILE OPEN: ERR= 6
IMEAN1: Image= 3C274_98A_43.POLI . 76 5 xywind= 1 1 1024 1024
IMEAN1: Mean and rms found by fitting peak in histogram:
IMEAN1: Mean= 8.7791E-05 Rms= 6.0969E-05 **** from histogram
IMEAN1: Mean and rms found by including all data:
IMEAN1: Mean= 1.211E-04 Rms= 8.636E-05 JY/BEAM over 1048576. pixels
IMEAN1: Flux density = 1.5218E+00 Jy. beam area = 83.44 pixels
IMEAN1: Minimum= 1.1830E-07 at 929 642 1 1
IMEAN1: Skypos: RA 12 30 49.42190799 DEC 12 23 28.0481750
IMEAN1: Maximum= 3.3543E-03 at 314 516 1 1
IMEAN1: Skypos: RA 12 30 49.42337722 DEC 12 23 28.0437650
IMEAN1: Skypos: PPOL 43151.852 MHZ
IMEAN1: Purports to die of UNNATURAL causes
IMEAN1: localhos 31DEC12 NEW: Cpu= 0.0 Real= 0 IO= 4
AIPS 1: Resumes
No clue what the problem is with a file open. In TST with no
change of inputs, it is happy and gives the same answers.
BACK TO COMB REPORT
Here is the header for the POLCT file after IMEAN:
AIPS 1: Image=3C274 (MA) Filename=3C274_98A_43.POLCT . 76
AIPS 1: Telescope=VLBA Receiver=VLBA
AIPS 1: Observer=BW098A User #= 604
AIPS 1: Observ. date=12-JAN-2013 Map date=19-MAR-2013
AIPS 1: Minimum=-1.05711959E-04 Maximum= 3.35333752E-03 JY/BEAM
AIPS 1: ----------------------------------------------------------------
AIPS 1: Type Pixels Coord value at Pixel Coord incr Rotat
AIPS 1: RA---SIN 1024 12 30 49.423 312.00 -0.000035 0.00
AIPS 1: DEC--SIN 1024 12 23 28.044 513.00 0.000035 0.00
AIPS 1: FREQ 1 4.3151852E+10 1.00 1.9200000E+08 0.00
AIPS 1: STOKES 1 5.0000000E+00 1.00 1.0000000E+00 0.00
AIPS 1: ----------------------------------------------------------------
AIPS 1: Coordinate equinox 2000.00
AIPS 1: Map type=NORMAL Number of iterations= 888
AIPS 1: Conv size= 0.00045 X 0.00020 Position angle= -15.48
AIPS 1: Phase shifted in X -0.007 in Y 0.000
AIPS 1: Rest freq 0.000 Vel type: OPTICAL wrt YOU
AIPS 1: Alt ref. value 0.00000E+00 wrt pixel 0.00
AIPS 1: Maximum version number of extension files of type HI is 1
AIPS 1: Maximum version number of extension files of type CG is 1
AIPS 1: Keyword = 'CCFLUX ' value = 9.865344E-03
AIPS 1: Keyword = 'CCTOTAL ' value = 9.865341E-03
AIPS 1: Keyword = 'OLDRFQ ' value = 4.30560000D+10
AIPS 1: Keyword = 'MAXABSU ' value = -1.000000E+00
AIPS 1: Keyword = 'SOURNAM1' value = '3C274 '
AIPS 1: Keyword = 'SOURNAM2' value = ' '
AIPS 1: Keyword = 'ACTNOISE' value = 4.457602E-05
AIPS 1: Keyword = 'ACTMEAN ' value = -2.481640E-05
So it is different from POLI and has the expected negative
minimum. It would seem that the problem was introduced
in TST.
Cheers,
Craig
---------------------------------------------------------------------
R. Craig Walker Array Operations Center
cwalker at nrao.edu National Radio Astronomy Observatory
Phone 575 835 7247 P. O. Box O
Fax 575 835 7027 Socorro NM 87801 USA
--------------------------------------------------------------web
More information about the Daip
mailing list