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
>
>
>
|