Hi David,
I've put them here, so let me know when you've picked them up, as its
taking up my home space :
http://www.mrao.cam.ac.uk/~buckle/GBS/Orion/
cheers,
Jane
David Berry wrote:
> 2009/5/11 Jane Buckle <[log in to unmask]>:
>
>> Dear All,
>>
>> I am seeing startlingly different behaviour when between makecube and
>> wcsmosaic (see the two attached images, which were created from the same
>> data using these two routines) when I use them to create cubes from four
>> files that have different spectral extents, although they are repeated
>> observations of the same source. I have changed the way that I clean up the
>> raw datafiles, and so one pair of files has the original spectral extent,
>> while the other pair have had about 30% clipped from the ends of the
>> spectra. The actual makecube command that I use is
>> "makecube in=/mraosw/data4/harp/mjlsg11/raw_cleaned/a20081111_00036_01_0001
>> out=temp autogrid=false detectors=! genvar=tsys outcat=! inweight=true
>> pixsize=\[7.2,7.2\] sparse=false trim=true weights=false spread=nearest
>> params=\[1.0,1.33\] system=tracking usedetpos=false"
>>
>> If I use makecube, then mfittrend with the following commands (from
>> hislist), I get pictureA:
>>
>> 1: 2009 May 11 15:32:21.519 - MAKECUBE (SMURF V0.4.1)
>>
>> Parameters: ALIGNSYS=FALSE AUTOGRID=FALSE BADMASK='OR'
>> CROTA=-55.000000000042 DETECTORS=!
>> FBL=[1.5135963227979,0.00037905821943821]
>> FBR=[1.508831094304,0.0071842303514043]
>> FLBND=[1.508806786544,0.00035475123779951,-424.20370675136]
>> FUBND=[1.5193107862253,0.01119288100032,442.75726064658]
>> FTL=[1.5192864791151,0.0043632445261451]
>> FTR=[1.5145213640061,0.011168574035431] GENVAR='tsys' [log in to unmask]
>> INWEIGHT=TRUE LBND=[-118,-99] LBOUND=[-118,-99,-512] NTILE=1
>> NPOLBIN=1 OUT=@temp OUTCAT=! PARAMS=[1,1.33] PIXSIZE=[7.2,7.2]
>> MSG_FILTER='NORM' REF=! REFLAT='0:19:53' REFLON='5:47:00.0'
>> SPARSE=FALSE SPECBOUNDS='-423.7804 442.3339' TILEDIMS=! TRIM=TRUE
>> WEIGHTS=FALSE SPREAD='nearest' SYSTEM='tracking' UBND=[120,100]
>> UBOUND=[120,100,511] USEDETPOS=FALSE
>> Software: /soft/star/bin/smurf/makecube
>>
>> 2: 2009 May 11 15:32:35.914 - MFITTREND (KAPPA 1.9-11)
>>
>> Parameters: AUTO=TRUE AXIS=@3 CLIP=[2,2,2.5,3] FITTYPE='Polynomial'
>> IN=@temp MASK=! METHOD='Single' MODIFYIN=FALSE NUMBIN=32 ORDER=1
>> OUT=@ngc2071_12co_b RMSCLIP=! SUBTRACT=TRUE TITLE=! VARIANCE=TRUE
>> Software: /soft/star/bin/kappa/mfittrend
>>
>> If I use makecube and mfittrend on the individual files, and then wcsmosaic
>> (hislist below) to put the final cubes together, I get pictureB.
>>
>> 1: 2009 May 11 14:44:22.095 - MAKECUBE (SMURF V0.4.1)
>>
>> Parameters: ALIGNSYS=FALSE AUTOGRID=FALSE BADMASK='OR'
>> CROTA=-55.000000000042 DETECTORS=!
>> FBL=[1.5137192681308,0.00050775633035526]
>> FBR=[1.5089740600772,0.0072843421205392]
>> FLBND=[1.5089497522708,0.00048344931397717,-424.20370675136]
>> FUBND=[1.5190906087172,0.011052734630646,442.75726064658]
>> FTL=[1.5190663015624,0.0042516937035154]
>> FTR=[1.5143211976496,0.011028427618018] GENVAR='tsys'
>> IN=@/mraosw/data4/harp/mjlsg11/raw_cleaned/a20081111_00036_01_000?
>> INWEIGHT=TRUE LBND=[-117,-94] LBOUND=[-117,-94,-512] NTILE=1
>> NPOLBIN=1 OUT=@temp OUTCAT=! PARAMS=[1,1.33] PIXSIZE=[7.2,7.2]
>> MSG_FILTER='NORM' REF=! REFLAT='0:19:53' REFLON='5:47:00.0'
>> SPARSE=FALSE SPECBOUNDS='-423.7804 442.3339' TILEDIMS=! TRIM=TRUE
>> WEIGHTS=FALSE SPREAD='nearest' SYSTEM='tracking' UBND=[120,93]
>> UBOUND=[120,93,511] USEDETPOS=FALSE
>> Software: /star/bin/smurf/makecube
>>
>> 2: 2009 May 11 14:44:36.318 - MFITTREND (KAPPA 1.9-11)
>>
>> Parameters: AUTO=TRUE AXIS=@3 CLIP=[2,2,2.5,3] FITTYPE='Polynomial'
>> IN=@temp MASK=! METHOD='Single' MODIFYIN=FALSE NUMBIN=32 ORDER=1
>> OUT=@mc_20081111_00036_01 RMSCLIP=! SUBTRACT=TRUE TITLE=!
>> VARIANCE=TRUE
>> Software: /star/bin/kappa/mfittrend
>>
>> 3: 2009 May 11 14:48:29.082 - WCSMOSAIC (KAPPA 1.9-11)
>>
>> Parameters: ACC=0.05
>> FLBND=[1.5089097074157,0.00048344931397717,-424.20370675136]
>> FUBND=[1.5195481106442,0.011430254274311,442.75726064658]
>> GENVAR=FALSE ILEVEL=2 [log in to unmask] LBND=!
>> LBOUND=[-117,-99,-512] MAXPIX=1000 METHOD='nearest'
>> OUT=@test_ngc2071_12co REF=! UBND=! UBOUND=[122,104,511]
>> VARIANCE=TRUE WLIM=1E-10
>> Software: /star/bin/kappa/wcsmosaic
>>
>>
>> Both of these are run on the same 4 files:
>>
>>
>
> Hi Jane.
>
>
>
>> /mraosw/data4/harp/mjlsg11/raw_cleaned/a20081111_00036_01_0001
>> /mraosw/data4/harp/mjlsg11/raw_cleaned/a20081111_00037_01_0001
>> /mraosw/data4/harp/mjlsg11/raw_cleaned/a20090131_00028_01_0001
>> /mraosw/data4/harp/mjlsg11/raw_cleaned/a20090131_00029_01_0001
>>
>
> Could you put these four files somewhere where I/we can get at them?
>
> Thanks,
>
> David
>
>
>
>
>> the files from 20081111 have a spectral extent, from ndftrace:
>> No. of dimensions: 3
>> Dimension size(s): 1024 x 16 x 5488
>> Pixel bounds : 1:1024, 1:16, 1:5488
>> Total pixels : 89915392
>> ........
>> Frame title : "3-d compound coordinate system"
>> Domain : DSBSPECTRUM-SPACEINDEX-TIME
>> First pixel centre : -423.7804, 1, 2008-11-11 12:45:27
>>
>> while the files from 20090131 have:
>> No. of dimensions: 3
>> Dimension size(s): 716 x 12 x 5805
>> Pixel bounds : 155:870, 1:12, 1:5805
>> Total pixels : 49876560
>> ........
>> Frame title : "3-d compound coordinate system"
>> Domain : DSBSPECTRUM-SPACEINDEX-TIME
>> First pixel centre : -293.3614, 1, 2009-01-31 08:26:41
>>
>>
>> Any clues about this? What this looks like visually is the detector gains
>> problem, but it clearly is not. Is it possible that part of the 'gains'
>> issue was that the spectral dimension was not being identically written for
>> all files, and makecube needs files that have identical spectral extents?
>>
>> Anyway, I hope for enlightenment, and a method of reducing files of varying
>> spectral extents.
>>
>>
>> Cheers,
>>
>> Jane
>>
>> --
>>
>> Dr Jane V. Buckle
>> Astrophysics Group, Cavendish Lab, J J Thomson Avenue, Cambridge CB3 0HE
>> Tel: +44 (0)1223 337298
>>
>>
>>
>>
--
Dr Jane V. Buckle
Astrophysics Group, Cavendish Lab, J J Thomson Avenue, Cambridge CB3 0HE
Tel: +44 (0)1223 337298
|