** 10 June 2002
on ev6
ls -ld /*/*/*F22*I* /*/*/*/*F22*I* | grep -v .fits> fred
rm fred2
awk '{print "du -s" ,$9,">>fred2"}' fred > fredx
source fredx
awk '{print $1}' fred2 > fred3
paste fred3 fred | awk '$1>1000' >relfiles
clean
(edit by hand)
(cutting the Decoupe stuff out: it is all stacked images)
995786 drwxr-xr-x 2 mellier crap 1536 Apr 28 18:14 /evsix2/mellier/DecoupeMarch2001/F22P5-I/
995786 drwxr-xr-x 2 mellier crap 1536 Apr 28 18:14 /i2/mellier/DecoupeMarch2001/F22P5-I/
1014370 drwxr-xr-x 2 mellier crap 1536 May 1 11:42 /evsix2/mellier/DecoupeMarch2001/F22P7-I/
1014370 drwxr-xr-x 2 mellier core 1536 May 1 11:42 /i2/mellier/DecoupeMarch2001/F22P7-I/
1013546 drwxr-xr-x 2 mellier crap 1536 May 1 12:06 /evsix2/mellier/DecoupeMarch2001/F22P8-I/
1013546 drwxr-xr-x 2 mellier crap 1536 May 1 12:06 /i2/mellier/DecoupeMarch2001/F22P8-I/
1012330 drwxr-xr-x 2 mellier crap 1536 May 1 12:31 /evsix2/mellier/DecoupeMarch2001/F22P9-I/
1012330 drwxr-xr-x 2 mellier crap 1536 May 1 12:31 /i2/mellier/DecoupeMarch2001/F22P9-I/
1019898 drwxr-xr-x 2 mellier crap 1536 May 2 09:06 /evsix2/mellier/DecoupeMarch2001/F22P10-I/
1019898 drwxr-xr-x 2 mellier crap 1536 May 2 09:06 /i2/mellier/DecoupeMarch2001/F22P10-I/
1011306 drwxr-xr-x 2 mellier core 1536 May 1 10:52 /i2/mellier/DecoupeMarch2001/F22P11-I/
1011306 drwxr-xr-x 2 mellier core 1536 May 1 10:52 /evsix2/mellier/DecoupeMarch2001/F22P11-I/
1012010 drwxr-xr-x 2 mellier crap 1536 May 1 11:17 /i2/mellier/DecoupeMarch2001/F22P12-I/
1012010 drwxr-xr-x 2 mellier crap 1536 May 1 11:17 /evsix2/mellier/DecoupeMarch2001/F22P12-I/
The good stuff:
5331340 drwxr-xr-x 3 mellier core 512 Jan 16 2000 /samplix3/mellier/F22P3-I-May99/
5360250 drwxr-xr-x 3 dantel system 512 Oct 25 2000 /samplix1/dantel/cfh12k_sept_oct_2000/F22P3b-I-Sep00/
5343722 drwxr-xr-x 3 dantel system 512 Oct 25 2000 /samplix1/dantel/cfh12k_sept_oct_2000/F22P4b-I-Sep00/
6008898 drwxr-xr-x 4 mellier core 512 Jan 14 2000 /samplix3/mellier/F22P4-I-Nov99/
5350106 drwxr-xr-x 3 dantel system 512 Oct 25 2000 /samplix1/dantel/cfh12k_sept_oct_2000/F22P5b-I-Sep00/
7519582 drwxr-xr-x 3 dantel system 512 Oct 25 2000 /samplix1/dantel/cfh12k_sept_oct_2000/F22P6b-I-Oct00/
3936970 drwxr-xr-x 3 mellier system 512 Jun 21 2000 /samplix3/mellier/F22P6b-I-Apr00-2/
7509296 drwxr-xr-x 3 dantel system 512 Oct 25 2000 /samplix1/dantel/cfh12k_sept_oct_2000/F22P7b-I-Oct00/
4724364 drwxr-xr-x 3 mellier system 512 Jun 21 2000 /samplix3/mellier/F22P11b-I-Apr00-2/
Note: mellier and dantel files *not* necessarily the same.
P1 Lost due to moon
P2 Lost due to moon
P3 /samplix3/mellier/F22P3-I-May99/notcoadded/*MDFG.fits.gz (5*12)
/samplix1/dantel/cfh12k_sept_oct_2000/F22P3b-I-Sep00/notcoadded/*MOBDFG.fits.gz (5*12)
P4 /samplix3/mellier/F22P4-I-Nov99/notcoadded/*MDFG.fits.gz (5*12)
/samplix1/dantel/cfh12k_sept_oct_2000/F22P4b-I-Sep00/notcoadded/*MOBDFG.fits.gz (5*12)
P5 /samplix1/dantel/cfh12k_sept_oct_2000/F22P5b-I-Sep00/notcoadded/*MOBDFG.fits.gz (5*12)
P6 /samplix3/mellier/F22P6b-I-Apr00-2/notcoadded/flatnotcoadded/*MODFOfl.fits (5*12)
/samplix1/dantel/cfh12k_sept_oct_2000/F22P6b-I-Oct00/notcoadded/*MOBDFG.fits.gz (7*12)
P7 /samplix1/dantel/cfh12k_sept_oct_2000/F22P7b-I-Oct00/notcoadded/*MOBDFG.fits.gz (7*12)
P8 nada
P9 nada
P10 nada
P11 /samplix3/mellier/F22P11b-I-Apr00-2/notcoadded/flatnotcoadded/*MODFOfl.fits (6*12)
P12 nada
New note added: 22 June 2002:
Somwhere in here you should run /samplix1/gwyn/findim.pl
In theory:
P1:
1.485930/0720/0.9/May99(01)[Lost: Moon]
2.485931/0720/0.9/May99(01)[Lost: Moon]
3.485932/0720/1.0/May99(01)[Lost: Moon]
4.485933/0720/0.9/May99(01)[Lost: Moon]
5.485934/0720/0.9/May99(01)[Lost: Moon]
6.486004/0720/1.2/May99(01)[Lost: Moon]
7.486005/0720/1.2/May99(01)[Lost: Moon]
P2:
1.486006/0720/1.0/May99(02)[Lost: Moon]
2.486007/0720/0.9/May99(02)[Lost: Moon]
3.486008/0720/0.9/May99(02)[Lost: Moon]
4.486009/0720/0.9/May99(02)[Lost: Moon]
5.486090/0720/0.8/May99(02)[Lost: Moon]
6.486091/0720/0.7/May99(02)[Lost: Moon]
7.486092/0720/0.8/May99(02)[Lost: Moon]
P3:
1.557282/0720/0.80/Nov00(03b)[ignore]
2.550361/0720/0.60/Sept00(03b)
3.550362/0720/0.60/Sept00(03b)
4.550363/0720/0.60/Sept00(03b)
5.550364/0720/0.75/Sept00(03b)[Lost: elongated]
6.550366/0720/0.65/Sept00(03b)
7.486093/0720/0.7/May99(03)[Lost: Moon]
8.486094/0720/0.6/May99(03)[Lost: Moon]
9.486370/0720/1.1/May99(03)
10.486371/0720/1.1/May99(03)
11.486372/0720/1.1/May99(03)
12.486373/0720/1.1/May99(03)
13.486374/0720/1.1/May99(03)
P4:
1.550367/0720/0.75/Sept00(04b)
2.550368/0720/0.75/Sept00(04b)
3.550369/0720/0.80/Sept00(04b)
4.550373/0720/0.60/Sept00(04b)
5.550374/0720/0.60/Sept00(04b)
6.508265/0720/0.7/Nov99(04)
7.508266/0720/0.8/Nov99(04)
8.508267/0720/0.8/Nov99(04)
9.508268/0720/0.8/Nov99(04)
10.508269/0720/0.8/Nov99(04)
P5:
1.550375/0720/0.55/Sept00(05b)
2.550376/0720/0.55/Sept00(05b)
3.550377/0720/0.60/Sept00(05b)
4.550378/0720/0.60/Sept00(05b)
5.550379/0720/0.65/Sept00(05b)
P6:
1.531288/0720/0.8/Apr00-2(06b)
2.531289/0720/0.7/Apr00-2(06b)
3.531290/0720/0.8/Apr00-2(06b)
4.531182/0720/1.0/Apr00-2(06b)
5.531183/0720/1.0/Apr00-2(06b)
6.551371/0720/0.7/Oct00-1(06b)
7.551372/0720/0.6/Oct00-1(06b)
8.551373/0720/0.6/Oct00-1(06b)
9.551374/0720/0.6/Oct00-1(06b)
10.551375/0720/0.7/Oct00-1(06b)
11.551376/0720/0.7/Oct00-1(06b)
12.551377/0720/0.8/Oct00-1(06b)
P7:
1.551379/0720/0.6/Oct00-1(07b)
2.551380/0720/0.6/Oct00-1(07b)
3.551381/0720/0.65/Oct00-1(07b)
4.551382/0720/0.65/Oct00-1(07b)
5.551383/0720/0.65/Oct00-1(07b)
6.551384/0720/0.6/Oct00-1(07b)
7.551385/0720/0.8/Oct00-1(07b)[Double images]
P8:
1.552797/0720/0.6/Oct00-3(8b)[Cirrus]
2.552798/0720/0.6/Oct00-3(8b)[Cirrus]
3.552799/0720/0.5/Oct00-3(8b)[Cirrus]
4.552800/0720/0.5/Oct00-3(8b)[Cirrus]
5.552801/0720/0.5/Oct00-3(8b)[Cirrus]
6.552855/0720/0.8/Oct00-3(8b)[Photometric]
P9:
1.552802/0720/0.5/Oct00-3(09b)
2.552803/0720/0.5/Oct00-3(09b)
3.552804/0720/0.5/Oct00-3(09b)
4.552805/0720/0.5/Oct00-3(09b)
5.552806/0720/0.5/Oct00-3(09b)
6.552856/0720/1.0/Oct00-3(09b)[Photometric]
P10:
1.552850/0720/0.8/Oct00-3(10b)
2.552851/0720/0.6/Oct00-3(10b)
3.552852/0720/0.6/Oct00-3(10b)
4.552853/0720/0.7/Oct00-3(10b)
5.552854/0720/0.8/Oct00-3(10b)
P11:
1.557283/0720/0.9/Nov00(11b)[ignore]
2.552953/0720/0.6/Oct00-3(11b)
3.552954/0720/0.6/Oct00-3(11b)
4.552858/0720/0.8/Oct00-3(11b)[Photometric]
5.552859/0720/0.8/Oct00-3(11b)[Photometric]
6.552860/0720/0.8/Oct00-3(11b)[Photometric]
7.531630/0720/1.0/Apr00-2(11b)
8.531631/0720/0.9/Apr00-2(11b)
9.531632/0720/0.8/Apr00-2(11b)
10.531633/0720/0.8/Apr00-2(11b)
11.531477/0720/1.3/Apr00-2(11b)[Lost: Clouds]
12.531478/0720/1.5/Apr00-2(11b)[Lost: Clouds]
P12:
1.552955/0720/0.6/Oct00-3(12b)
2.552956/0720/0.7/Oct00-3(12b)
3.552957/0720/0.7/Oct00-3(12b)
4.552958/0720/0.8/Oct00-3(12b)
5.552959/0720/1.0/Oct00-3(12b)
/samplix1/dantel/DATAprocessed-avril2000-2.DEF/astrometry/F22P11b-I-Apr00-2:
/samplix1/dantel/DATAprocessed-avril2000-2.DEF/astrometry/F22P6b-I-Apr00-2:
/samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p11b_I/notcoadded/flatnotcoadded:
/samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p11b_I/notcoadded:
/samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p6b_I/notcoadded/cut:
/samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p6b_I/notcoadded/flatnotcoadded:
/samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p6b_I/notcoadded:
/samplix1/dantel/DATAprocessed-avril2000-2.DEF/images_MODF:
/samplix1/dantel/cfh12k_sept_oct_2000/F22P3b-I-Sep00/notcoadded/flatnotcoadded:
/samplix1/dantel/cfh12k_sept_oct_2000/F22P3b-I-Sep00/notcoadded:
/samplix1/dantel/cfh12k_sept_oct_2000/F22P4b-I-Sep00/notcoadded/flatnotcoadded:
/samplix1/dantel/cfh12k_sept_oct_2000/F22P4b-I-Sep00/notcoadded:
/samplix1/dantel/cfh12k_sept_oct_2000/F22P5b-I-Sep00/notcoadded/flatnotcoadded:
/samplix1/dantel/cfh12k_sept_oct_2000/F22P5b-I-Sep00/notcoadded:
/samplix1/dantel/cfh12k_sept_oct_2000/F22P6b-I-Oct00/notcoadded/flatnotcoadded:
/samplix1/dantel/cfh12k_sept_oct_2000/F22P6b-I-Oct00/notcoadded:
/samplix1/dantel/cfh12k_sept_oct_2000/F22P7b-I-Oct00/notcoadded/flatnotcoadded/bad:
/samplix1/dantel/cfh12k_sept_oct_2000/F22P7b-I-Oct00/notcoadded/flatnotcoadded:
/samplix1/dantel/cfh12k_sept_oct_2000/F22P7b-I-Oct00/notcoadded:
/samplix1/dantel/cfh12k_sept_oct_2000/astrometry/oct2000/F22P6b-I-Oct00:
/samplix1/dantel/cfh12k_sept_oct_2000/astrometry/oct2000/F22P7b-I-Oct00:
/samplix1/dantel/cfh12k_sept_oct_2000/astrometry/sept2000/F22P3b-I-Sep00:
/samplix1/dantel/cfh12k_sept_oct_2000/astrometry/sept2000/F22P4b-I-Sep00:
/samplix1/dantel/cfh12k_sept_oct_2000/astrometry/sept2000/F22P5b-I-Sep00:
/samplix1/foucaud/Oct2000_3/Results/f22/f22p10b_I:
/samplix1/foucaud/Oct2000_3/Results/f22/f22p11b_I:
/samplix1/foucaud/Oct2000_3/Results/f22/f22p12b_I:
/samplix1/foucaud/Oct2000_3/Results/f22/f22p8b_I:
/samplix1/foucaud/Oct2000_3/Results/f22/f22p9b_I:
/samplix3/dantel/cfh12k_oct2000-MOBDF:
/samplix3/dantel/cfh12k_sep2000-MOBDF:
/samplix3/dantel/coadd-oct2000_2_3/F22/F22P10b-I-Oct00_3:
/samplix3/dantel/coadd-oct2000_2_3/F22/F22P11b-I-Oct00_3:
/samplix3/dantel/coadd-oct2000_2_3/F22/F22P12b-I-Oct00_3:
/samplix3/dantel/coadd-oct2000_2_3/F22/F22P8b-I-Oct00_3:
/samplix3/dantel/coadd-oct2000_2_3/F22/F22P9b-I-Oct00_3:
/samplix3/mellier/F22P11b-I-Apr00-2/notcoadded/flatnotcoadded:
/samplix3/mellier/F22P3-I-May99/notcoadded/flatnotcoadded:
/samplix3/mellier/F22P3-I-May99/notcoadded:
/samplix3/mellier/F22P4-I-Nov99/notcoadded/flatnotcoadded:
/samplix3/mellier/F22P4-I-Nov99/notcoadded:
/samplix3/mellier/F22P6b-I-Apr00-2/notcoadded/flatnotcoadded:
/u2/morin/terapix/dbterapix/v2/src/dbtest/tpx/data:
(less the *fl.fits)
/samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p11b_I/notcoadded:
/samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p6b_I/notcoadded/cut:
/samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p6b_I/notcoadded:
/samplix1/dantel/DATAprocessed-avril2000-2.DEF/images_MODF:
/samplix1/dantel/cfh12k_sept_oct_2000/F22P3b-I-Sep00/notcoadded:
/samplix1/dantel/cfh12k_sept_oct_2000/F22P4b-I-Sep00/notcoadded:
/samplix1/dantel/cfh12k_sept_oct_2000/F22P5b-I-Sep00/notcoadded:
/samplix1/dantel/cfh12k_sept_oct_2000/F22P6b-I-Oct00/notcoadded:
/samplix1/dantel/cfh12k_sept_oct_2000/F22P7b-I-Oct00/notcoadded:
/samplix1/foucaud/Oct2000_3/Results/f22/f22p10b_I:
/samplix1/foucaud/Oct2000_3/Results/f22/f22p11b_I:
/samplix1/foucaud/Oct2000_3/Results/f22/f22p12b_I:
/samplix1/foucaud/Oct2000_3/Results/f22/f22p8b_I:
/samplix1/foucaud/Oct2000_3/Results/f22/f22p9b_I:
/samplix3/dantel/cfh12k_oct2000-MOBDF:
/samplix3/dantel/cfh12k_sep2000-MOBDF:
/samplix3/mellier/F22P3-I-May99/notcoadded:
/samplix3/mellier/F22P4-I-Nov99/notcoadded:
/u2/morin/terapix/dbterapix/v2/src/dbtest/tpx/data:
** 11 June 2002
Let's go with
/samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p6b_I/notcoadded/cut( = cropped)
/samplix1/dantel/DATAprocessed-avril2000-2.DEF/images_MODF (includes other images)
/samplix3/dantel/cfh12k_oct2000-MOBDF (other images)
/samplix3/dantel/cfh12k_sep2000-MOBDF (other images)
/u2/morin/terapix/dbterapix/v2/src/dbtest/tpx/data (No good: just headers)
/samplix1/dantel/cfh12k_sept_oct_2000/F22P3b-I-Sep00/notcoadded
/samplix3/mellier/F22P3-I-May99/notcoadded
/samplix1/dantel/cfh12k_sept_oct_2000/F22P4b-I-Sep00/notcoadded
/samplix3/mellier/F22P4-I-Nov99/notcoadded
/samplix1/dantel/cfh12k_sept_oct_2000/F22P5b-I-Sep00/notcoadded
/samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p6b_I/notcoadded
/samplix1/dantel/cfh12k_sept_oct_2000/F22P6b-I-Oct00/notcoadded
/samplix1/dantel/cfh12k_sept_oct_2000/F22P7b-I-Oct00/notcoadded
/samplix1/foucaud/Oct2000_3/Results/f22/f22p8b_I
/samplix1/foucaud/Oct2000_3/Results/f22/f22p9b_I
/samplix1/foucaud/Oct2000_3/Results/f22/f22p10b_I
/samplix1/foucaud/Oct2000_3/Results/f22/f22p11b_I
/samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p11b_I/notcoadded
/samplix1/foucaud/Oct2000_3/Results/f22/f22p12b_I
(Fav's)
/samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p6b_I/notcoadded
/samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p11b_I/notcoadded
All available in
/samplix1/dantel/DATAprocessed-avril2000-2.DEF/images_MODF
According to Seb, don't use the
Fav's...
But according to Mirelle Dantel, do use them...
Final list:
cat >F22I.dirlist
/samplix1/dantel/cfh12k_sept_oct_2000/F22P3b-I-Sep00/notcoadded
/samplix3/mellier/F22P3-I-May99/notcoadded
/samplix1/dantel/cfh12k_sept_oct_2000/F22P4b-I-Sep00/notcoadded
/samplix3/mellier/F22P4-I-Nov99/notcoadded
/samplix1/dantel/cfh12k_sept_oct_2000/F22P5b-I-Sep00/notcoadded
/samplix1/dantel/cfh12k_sept_oct_2000/F22P6b-I-Oct00/notcoadded
/samplix1/dantel/cfh12k_sept_oct_2000/F22P7b-I-Oct00/notcoadded
/samplix1/foucaud/Oct2000_3/Results/f22/f22p8b_I
/samplix1/foucaud/Oct2000_3/Results/f22/f22p9b_I
/samplix1/foucaud/Oct2000_3/Results/f22/f22p10b_I
/samplix1/foucaud/Oct2000_3/Results/f22/f22p11b_I
/samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p11b_I/notcoadded
/samplix1/foucaud/Oct2000_3/Results/f22/f22p12b_I
/samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p6b_I/notcoadded
/samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p11b_I/notcoadded
/samplix1/dantel/DATAprocessed-avril2000-2.DEF/images_MODF
run: /samplix/gwyn/track/linkall.pl
Put the links in F22I.links
552959o11oct003p12.fits
^ ^^ ^ ^ ^^^
| || | | |||____ pointing number: 2 digits 03 to 12 (should be 01 to 16)
| || | | ||_____ the letter 'p'
| || | | |______ date: 1 digit for sub-run number
| || | |________ date: 2 digits for year
| || |___________ date:
| ||_____________ CCD number: 2 digits 00 to 11
| |______________ the letter 'o'
|____________________ the exposure number: 6 digits
3 letters for month, 2 digits for year 1 digit
switch to pix
cd /raid
mkdir gwyn
cd gwyn
mkdir F22
cd F22
mkdir F22I
cd F22I
mkdir data
cd data
ncftp -u gwyn samplix.iap.fr
cd /samplix1/gwyn/F22/F22I
get *
bye
gzip -d *gz
** 13 June 2002
Redid most the above with the Fav images.
in iraf:
cd F22I/data
!ls *.fits | perl -nle 'print $_."\[1:2048,1:4096\]"' >listin
!ls *.fits | perl -nle 's/\.fits/C.fits/ ; print $_' >listout
imcopy @listin @listout
** 14 June 2002
instead of all that:
rm imlist*
stack.pl
stackall
** 17 June 2002
ls *N.fits | perl -nle '($a=$_)=~s/N.fits/.fits/; rename $_,"../$a"'
mv mask* ..
mv imlist* ..
cd ..
sex.pl
** 19 June 2002
Bug in linkall.pl means I have to do a redo from start.
On ev6:
cd /samplix1/gwyn/track
linkall.pl
On pix1:
cd /raid/gwyn/F22/F22I/data
ncftp -u gwyn samplix.iap.fr
cd /samplix1/gwyn/F22/F22I
get *
bye
gzip -d *gz
mv * back
ln -s back/*.fits .
delhead CONSWV *.fits
rm imlist*
imlist.pl
stackall
2551.540u 197.060s 53:57.22 84.9% 0+0k 0+0io 245pf+0w
rm *p??.fits
ls *N.fits | perl -nle '($a=$_)=~s/N.fits/.fits/; rename $_,"../$a"'
cd ..
cp data/mask*.box .
imlist.pl
sex.pl >sex.log
(both imlist.pl and sex.pl are now pix1-based, not lasa28-based)
** 20 June 2002
maskch.pl
on lasa28:
scat -r -4000 -n 1000000 -d -c ua2 334.5 .33333 | awk '{print $2,$3,$4,$5}' >USNO
also using getout.f to get the outlines
this is checked with sm.1 on Morgan:
/home/gwyn/red/F22/F22I/inithead
Check all the masks with the help of Wayne & Schuster (1963)
00 images: most have a big box on the upper left that
should not be included in the final images, but could be
used in the astrometry part.
Also there are a lot very faint baddish columns
531182o04apr002p06.fits shows a flatfield gradient
0
486370o04may990p03.fits (which looks like the normal CCD00)
has the faint bad cols and the big bad square.
CCD05 problems with tvm...
486370o06may990p03.fits streaky stars
486370o07may990p03.fits strange effects: examine image without boxes
486370o10may990p03.fits lots of extensions
CCD00 leave as is: DONE
the blurry area in left side. might vary
from exposure to exposure
CCD01 leave as is: DONE
gaps in the bad columns in one exposures
CCD02 No problems DONE
CCD03 No problems DONE
CCD04 Leave as is DONE Blurry area in May99
CCD05 homogenize all but 990 DONE
may990 completely different
CCD06 All but may990 : homogenize: DONE
may990 fine as is.
CCD07 homogenize all by cating boxes DONE
CCD08 homogenize all: DONE (could have been ignored. The correction is the
strip along the bottom.
CCD09 homogenize all: DONE
CCD10 No problems DONE
CCD11 may990: added marginal bad column: keep it DONE
** June 21 2002
Script for dealing with CCD09
disp 486370o09may990p03.fits[700:850,3900:4096] 1 erase+ xcen=.25 ycen=.3333
tvmark 1 '' commands=add09may990 col=205
tvmark 1 '' commands=add09may990 col=204
disp 531182o09apr002p06.fits[700:850,3900:4096] 1 erase- xcen=.50 ycen=.3333
tvmark 1 '' commands=add09may990 col=205
tvmark 1 '' commands=add09apr002 col=204
disp 551371o09oct001p06.fits[700:850,3900:4096] 1 erase- xcen=.75 ycen=.3333
tvmark 1 '' commands=add09may990 col=205
tvmark 1 '' commands=add09oct001 col=204
disp 508265o09nov990p04.fits[700:850,3900:4096] 1 erase- xcen=.25 ycen=.6666
tvmark 1 '' commands=add09may990 col=205
tvmark 1 '' commands=add09nov990 col=204
disp 550361o09sep000p03.fits[700:850,3900:4096] 1 erase- xcen=.50 ycen=.6666
tvmark 1 '' commands=add09may990 col=205
tvmark 1 '' commands=add09sep000 col=204
disp 552797o09oct003p08.fits[700:850,3900:4096] 1 erase- xcen=.75 ycen=.6666
tvmark 1 '' commands=add09may990 col=205
tvmark 1 '' commands=add09oct003 col=204
Script for dealing with CCD08
disp 486370o08may990p03.fits[1:100,1:4096] 1 erase+ xcen=.15 ycen=.5 ymag=.2
tvmark 1 '' commands=mask08may990.tvm col=204
tvmark 1 '' commands=add08may990 col=205
disp 531182o08apr002p06.fits[1:100,1:4096] 1 erase- xcen=.30 ycen=.5 ymag=.2
tvmark 1 '' commands=mask08may990.tvm col=204
tvmark 1 '' commands=add08apr002 col=205
disp 551371o08oct001p06.fits[1:100,1:4096] 1 erase- xcen=.45 ycen=.5 ymag=.2
tvmark 1 '' commands=mask08may990.tvm col=204
tvmark 1 '' commands=add08oct001 col=205
disp 508265o08nov990p04.fits[1:100,1:4096] 1 erase- xcen=.60 ycen=.5 ymag=.2
tvmark 1 '' commands=mask08may990.tvm col=204
tvmark 1 '' commands=add08nov990 col=205
disp 550361o08sep000p03.fits[1:100,1:4096] 1 erase- xcen=.75 ycen=.5 ymag=.2
tvmark 1 '' commands=mask08may990.tvm col=204
tvmark 1 '' commands=add08sep000 col=205
disp 552797o08oct003p08.fits[1:100,1:4096] 1 erase- xcen=.90 ycen=.5 ymag=.2
tvmark 1 '' commands=mask08may990.tvm col=204
tvmark 1 '' commands=add08oct003 col=205
script for dealing with CCD07
disp 486370o07may990p03.fits[1550:1730,100:200] 1 erase+ xcen=.25 ycen=.3333
tvmark 1 '' commands=fred col=205
tvmark 1 '' commands=add07may990 col=204
disp 531182o07apr002p06.fits[1550:1730,100:200] 1 erase- xcen=.50 ycen=.3333
tvmark 1 '' commands=fred col=205
tvmark 1 '' commands=add07apr002 col=204
disp 551371o07oct001p06.fits[1550:1730,100:200] 1 erase- xcen=.75 ycen=.3333
tvmark 1 '' commands=fred col=205
tvmark 1 '' commands=add07oct001 col=204
disp 508265o07nov990p04.fits[1550:1730,100:200] 1 erase- xcen=.25 ycen=.6666
tvmark 1 '' commands=fred col=205
tvmark 1 '' commands=add07nov990 col=204
disp 550361o07sep000p03.fits[1550:1730,100:200] 1 erase- xcen=.50 ycen=.6666
tvmark 1 '' commands=fred col=205
tvmark 1 '' commands=add07sep000 col=204
disp 552797o07oct003p08.fits[1550:1730,100:200] 1 erase- xcen=.75 ycen=.6666
tvmark 1 '' commands=fred col=205
tvmark 1 '' commands=add07oct003 col=204
script for dealing with CCD01
disp 486370o00may990p03.fits[140:575,1400:4096] 1 erase+ xcen=.15 ycen=.5
tvmark 1 '' commands=mask00may990.tvm col=204
tvmark 1 '' commands=add00may990 col=205
disp 531182o00apr002p06.fits[140:575,1400:4096] 1 erase- xcen=.30 ycen=.5
tvmark 1 '' commands=mask00apr002.tvm col=204
tvmark 1 '' commands=add00apr002 col=205
disp 551371o00oct001p06.fits[140:575,1400:4096] 1 erase- xcen=.45 ycen=.5
tvmark 1 '' commands=mask00oct001.tvm col=204
tvmark 1 '' commands=add00oct001 col=205
disp 508265o00nov990p04.fits[140:575,1400:4096] 1 erase- xcen=.60 ycen=.5
tvmark 1 '' commands=mask00nov990.tvm col=204
tvmark 1 '' commands=add00nov990 col=205
disp 550361o00sep000p03.fits[140:575,1400:4096] 1 erase- xcen=.75 ycen=.5
tvmark 1 '' commands=mask00sep000.tvm col=204
tvmark 1 '' commands=add00sep000 col=205
disp 552797o00oct003p08.fits[140:575,1400:4096] 1 erase- xcen=.90 ycen=.5
tvmark 1 '' commands=mask00oct003.tvm col=204
tvmark 1 '' commands=add00oct003 col=205
look at the .gif's of the masks:
maskapr002.gif
maskmay990.gif
masknov990.gif
maskoct001.gif
maskoct003.gif
masksep000.gif
maskpic: takes the mask*tvm files (as specified
by imlistall) and makes
cute .pgm files of the masks.
mog.pl: takes these pgm and assembles them
into a big gif, one per epoch.
back to pix1
cd /raid/gwyn/F22/F22I/data
ncftp -u gwyn lasa28.astrsp-mrs.fr
cd /home/gwyn/red/F22/F22I/maskcheck/
get add*
bye
imlist.pl
stackall >stacklog
ls *N.fits | perl -nle '($a=$_)=~s/N.fits/.fits/; rename $_,"../$a"'
cd ..
cp data/mask*.box .
imlist.pl
sex.pl >sex.log
dofit1.pl dofit1.log
** 22 June 2002
531478o??apr002p11.fits
is completely blank it shall be discarded
But it should never have been read:
Ditto with the following:
There was a problem with linkall.pl
It has been fixed.
sep000 03 5.550364/0720/0.75/Sept00(03b)[Lost: elongated]
oct001 07 7.551385/0720/0.8/Oct00-1(07b)[Double images]
oct003 08 1.552797/0720/0.6/Oct00-3(8b)[Cirrus]
oct003 08 2.552798/0720/0.6/Oct00-3(8b)[Cirrus]
oct003 08 3.552799/0720/0.5/Oct00-3(8b)[Cirrus]
oct003 08 4.552800/0720/0.5/Oct00-3(8b)[Cirrus]
oct003 08 5.552801/0720/0.5/Oct00-3(8b)[Cirrus]
apr002 11 11.531477/0720/1.3/Apr00-2(11b)[Lost: Clouds]
apr002 11 12.531478/0720/1.5/Apr00-2(11b)[Lost: Clouds]
These images have been removed...
The process has been restarted though:
backup of linkall.pl is on Morgan:
/home/gwyn/red/scripts/linkall.pl
run:
rm 550364*
rm 551385*
rm 531477*
rm 531478*
Keeping the cirrus images (reluctantly)
on ev6
cd /samplix1/gwyn/track
linkall.pl
On pix1:
cd /raid/gwyn/F22/F22I/data
ncftp -u gwyn samplix.iap.fr
cd /samplix1/gwyn/F22/F22I
get *
bye
gzip -d *gz
mv * back
ln -s back/*.fits .
delhead CONSWV *.fits
rm imlist*
imlist.pl
stackall
ls *N.fits | perl -nle '($a=$_)=~s/N.fits/.fits/; rename $_,"../$a"'
cd ..
** 23 June 2002
cd /raid/gwyn/F22/F22I
cp data/mask*box
imlist.pl
sex.pl >sex.log
dofit1.pl >dofit1.log
check results...
checkdo.pl dofit1.log
worst centering errors:
486371o06may990p03.fits 75 -1 65 53 1.255
486373o06may990p03.fits 76 -1 64 48 1.241
486374o06may990p03.fits 72 -1 60 49 1.233
486372o06may990p03.fits 76 -1 65 52 1.217
486370o06may990p03.fits 75 -1 64 51 1.216
552954o06oct003p11.fits 75 -1 66 67 1.044
fewest initial matches:
550376o05sep000p05.fits 50 -1 42 43 0.645
550375o05sep000p05.fits 51 -1 43 45 0.717
550377o05sep000p05.fits 52 -1 43 45 0.694
550361o09sep000p03.fits 52 -1 30 30 0.570
550362o09sep000p03.fits 53 -1 32 32 0.557
550362o05sep000p03.fits 53 -1 50 50 0.781
fewest final matches:
550361o09sep000p03.fits 52 -1 30 30 0.570
550362o09sep000p03.fits 53 -1 32 32 0.557
552955o02oct003p12.fits 54 -1 33 33 0.674
550377o09sep000p05.fits 55 -1 33 34 0.725
550376o09sep000p05.fits 55 -1 33 34 0.654
552956o09oct003p12.fits 57 -1 35 35 0.529
Not too shabby.
Results of inspecting images:
486370o00may990p03.fits multiple dots
486371o00may990p03.fits multiple dots
486372o00may990p03.fits multiple dots
486373o00may990p03.fits multiple dots
486374o00may990p03.fits multiple dots
generally, the swath of crap is of variable width
551372o00oct001p06.fits streak
551383o00oct001p06.fits streak
552854o00oct003p10.fits streak
552958o00oct003p12.fits faint streak
551372o01oct001p06.fits streak
551383o01oct001p07.fits streak
552854o01oct003p10.fits streak
551372o02oct001p06.fits streak
551383o02oct001p06.fits streak
552854o02oct003p10.fits streak
552956o02oct003p12.fits faint streak
531288o01apr002p06.fits faint streak
551372o03oct001p06.fits streak
531288o04apr002p06.fits bad cosmic ray
550361o04sep000p03.fits streak
550362o04sep000p03.fits thread?
550363o04sep000p03.fits thread?
550366o04sep000p03.fits bad cosmic ray
550378o04sep000p05.fits scattered light
552806o04oct003p09.fits bad cosmic ray
552850o04oct003p10.fits tumbling satellite
552851o04oct003p10.fits bumpy cosmic rays
552856o04oct003p09.fits bad cosmic ray
552957o04oct003p12.fits big streak
486370o05may990p03.fits thread?
550379o05sep000p05.fits internal reflection
551372o05oct001p06.fits streak
551383o05oct001p07.fits streak
552801o05oct003p08.fits streak
552853o05oct003p10.fits streak
552954o05oct003p11.fits streak
486370o06may990p03.fits stars streaked
486371o06may990p03.fits stars streaked
486372o06may990p03.fits stars streaked
486373o06may990p03.fits stars streaked
486374o06may990p03.fits stars streaked
531183o06apr002p06.fits scattered light
531290o06apr002p06.fits scattered light
552802o06oct003p09.fits bright star just off frame
552803o06oct003p09.fits bright star just off frame
552804o06oct003p09.fits bright star just off frame
552805o06oct003p09.fits bright star just off frame
552806o06oct003p09.fits bright star just off frame
552855o06oct003p08.fits bright star just off frame
552856o06oct003p09.fits bright star just off frame
486371o07may990p03.fits slightly trailed images
486373o07may990p03.fits slightly trailed images
some other very slightly trailed image in may...
531183o07apr002p06.fits two "asteroids"
552802o07oct003p09.fits very bright star
552803o07oct003p09.fits very bright star
552804o07oct003p09.fits very bright star
552805o07oct003p09.fits very bright star
552806o07oct003p09.fits very bright star
552856o07oct003p09.fits very bright star
486370o08may990p03.fits very slightly trailed images
486371o08may990p03.fits very slightly trailed images
486372o08may990p03.fits slightly trailed images
486373o08may990p03.fits slightly trailed images
486374o08may990p03.fits slightly trailed images
550368o08sep000p04.fits bright star off to one side.
552806o08oct003p09.fits two horizontal lines
486370o09may990p03.fits streak on lower right corner
486371o09may990p03.fits very slightly trailed images
486373o09may990p03.fits slightly trailed images
531630o09apr002p11.fits bright star
531631o09apr002p11.fits bright star
531632o09apr002p11.fits bright star
531633o09apr002p11.fits bright star
550367o09sep000p04.fits bright star off edge
550368o09sep000p04.fits bright star off edge
550369o09sep000p04.fits bright star off edge
550373o09sep000p04.fits bright star off edge
550374o09sep000p04.fits bright star off edge
552956o09oct003p12.fits faint streak
486370o10may990o03.fits bright streak
550367o10sep000p04.fits scattered light at bottom
550368o10sep000p04.fits scattered light at bottom
550369o10sep000p04.fits scattered light at bottom
550373o10sep000p04.fits scattered light at bottom
550374o10sep000p04.fits scattered light at bottom
552957o10oct003p12.fits bright streak
may990o11 some more bits might need masking
486370o11may990o03.fits bright streak
*o11apr002p11.fits don't look good
531632o11apr002p11.fits doubled bright star
550ls 361o11sep000p03.fits faint streaks
550369o11sep000p04.fits ghostly doubled stars
550375o11sep000p05.fits ghostly residual stars
551372o11oct001p06.fits streak
551373o11oct001p06.fits bright streak
552954o11oct003p11.fits streak
And looking at the results of the first fit:
508269o03nov990p04.fits corner star unmatched
551371o06oct001p06.fits missing stars in bottom left
551372o06oct001p06.fits a few stars in bottom left
551373o06oct001p06.fits missing stars in bottom left
551374o06oct001p06.fits missing stars in bottom left
551375o06oct001p06.fits missing stars in bottom left
551376o06oct001p06.fits missing stars in bottom left
551377o06oct001p06.fits missing stars in bottom left
** 24 June 2002
Photometry:
I am assuming:
Y-date Header-date Table-date Zero-point
may990 1999-05-13 May 99 25.802
nov990 1999-11-09 Nov 99 25.894
apr002 2000-04-30 Apr00-2 26.3386
2000-05-01
2000-05-03
sep000 2000-09-22 Sept+Oct00 26.2906
oct001 2000-10-03 Sept+Oct00 26.2906
oct003 2000-10-23 Oct00-3 26.0078
2000-10-24
2000-10-25
dofit2.pl is being maintained on pix1
with this version of getzp
sub getzp {
my $head=shift @_;
my $expkey ='EXPTIME';
my $airkey ='AIRMASS';
my $datekey ='DATE-OBS';
my $colour='I';
my %fact= qw (B 0.20 V 0.12 R 0.10 I 0.08);
my %zpoint = qw {
I1999-05-13 25.802
I1999-11-09 25.894
I2000-04-30 26.3386
I2000-05-01 26.3386
I2000-05-03 26.3386
I2000-09-22 26.2906
I2000-10-03 26.2906
I2000-10-23 26.0078
I2000-10-24 26.0078
I2000-10-25 26.0078
};
open HEAD,$head;
while (
) {
chomp;
($key,$value)= split /=/;
($value) = split /\//,$value;
$key=~s/\s//g;
$value=~s/\s//g;
$value=~s/\'//g;
if ($key eq $expkey ) {$exptime = $value}
if ($key eq $airkey ) {$airmass = $value}
if ($key eq $datekey) {$date = $value}
}
close HEAD;
$fact=$fact{$colour};
unless (defined($zpoint{$colour.$date})) {
print "ERROR: couldn't find zero-point\n";
print "$head $colour $date\n";
die;
}
$zpoint =$zpoint{$colour.$date} + 2.5*&log10($exptime) - $airmass*$fact;
}
** 25 June 2002
Seb says:
Oct 00-2 No calibrations because of Supernovae dudes
Apr-00-2 first reduction run of Seb: Probably not good.
** 28 June 2002
checked gradients accross images.
Get the minibackground images from SExtractor.
use
/home/virmos/gwyn/scripts/back.pl
on pix1
Run backch.f on lasa28.astrsp-mrs.fr
on all the images.
backch.f prints out
image_name
median
peak-to-peak
horizontal_gradient
vertical_gradient
diagonal_gradient (sw-ne)
diagonal_gradient (se-nw)
The gradients should not be more than 2%
The peak to peak is actually affected by bright
stars more than anything else.
All images have less than 2% gradients.
Some have high medians:
531633o??apr002p11.fits 40000+
It shall be tossed...
** 1 July 2002
Checking the gain:
run /home/virmos/gwyn/scripts/noise.pl
(it runs SExtractor, throws away the catlogs
and outputs a object-less images with a mb.fits suffix)
run /home/virmos/gwyn/pg/stddev.f
ls *mb.fits | perl -nle 'system "stddev $_"' | grep -v read > stddev.log
This gives you the std. dev. by examining 10x10 pixel
sub-rasters, tossing those contaminated by objects or
bad columns, and averaging all the std-devs.
Using a small sub-raster means even if there is
fringing, you are probably either on a fringe or completely
off it. The average will change slightly, but the std-dev won't.
Copy back to lasa28
~/red/F22/F22I/backcheck
run sm.1 and sm.2
to get:
I seem to have deleted all the add-files
I don't think they got included in the last stackall.
Redo from start:
cp
stackall >stackall.log
ls *N.fits | perl -nle '($a=$_)=~s/N.fits/.fits/; rename $_,"../$a"'
cd ..
cd /raid/gwyn/F22/F22I
cp data/mask* .
imlist.pl
sex.pl >sex.log
Problems with one image.
Turns out sex.pl got cut off in mid stream. Not enough
stars --> not a good fit.
dofit1.pl dofit1.log
mv *06may* toss
dofit2.pl dofit2nomay06.log
crossref
dofit3.pl dofit3nomay06.log
** 9 July 2002
I'm supposed to be on vacation
set
PHOT_APERTURES 20 # MAG_APER aperture diameter(s) in pixels
sex.pl
dofit.pl USNO 1 dofit1.log
checkdo.pl dofit1.log
ls *.fits | perl -nle 'system "uphead $_ fead"'
dofit.pl USNO 2 dofit2.log
checkdo.pl dofit2.log
550376o03sep000p05.fits 74 2 52 65 0.727
508268o00nov990p04.fits 72 2 45 54 0.709
550377o03sep000p05.fits 75 2 49 63 0.701
508268o05nov990p04.fits 70 2 55 62 0.698
550376o04sep000p05.fits 57 2 42 52 0.690
486373o10may990p03.fits 61 2 46 56 0.689
fewest initial matches:
550366o05sep000p03.fits 49 2 42 48 0.574
550377o05sep000p05.fits 50 2 38 45 0.508
550376o05sep000p05.fits 50 2 38 42 0.498
550375o05sep000p05.fits 51 2 38 45 0.566
550362o09sep000p03.fits 51 2 30 32 0.390
550362o05sep000p03.fits 51 2 43 48 0.502
fewest final matches:
550361o09sep000p03.fits 52 2 30 32 0.393
552955o02oct003p12.fits 52 2 30 32 0.506
550376o09sep000p05.fits 55 2 28 32 0.415
550362o09sep000p03.fits 51 2 30 32 0.390
550377o09sep000p05.fits 55 2 29 33 0.409
552957o02oct003p12.fits 53 2 31 33 0.440
No real worries.
crossref
dofit.pl master 2 500 dofit3.log
cd astrocheck
ln -s ../*fead .
ln -s ../master
ln -s ../rdall.used
\ls *.fead | perl -nle 'system "getout $_"'
ls *out >outlist
cat *out | sort -n | head > outrel
cat *out | sort -nr | head >> outrel
cat *out | sort -nr +1 | head >> outrel
cat *out | sort -n +1 | head >> outrel
** 10 July 2002
Added a feature to stackall.f
Every pixel less than the median (less 2000)
is replaced with the saturation value: 65535
This will hopefully avoid the problem of solarized
saturated stars being identified as two stars.
stackall
ls *N.fits | perl -nle '($a=$_)=~s/N.fits/.fits/; rename $_,"../$a"'
cp mask*box ..
cd ..
mv *o06may*.fits toss
mv 531633o??apr002p11.fits toss
imlist.pl
sex.pl>sex.log
** 12 July 2002
It looks like the Oct00-3 I zero point is off
sub getzp {
my $head=shift @_;
my $expkey ='EXPTIME';
my $airkey ='AIRMASS';
my $datekey ='DATE-OBS';
my $colour='I';
my %fact= qw (B 0.20 V 0.12 R 0.10 I 0.08);
my %zpoint = qw {
I1999-05-13 25.802
I1999-11-09 25.894
I2000-04-30 26.3386
I2000-05-01 26.3386
I2000-05-03 26.3386
I2000-09-22 26.2906
I2000-10-03 26.2906
I2000-10-23 26.260
I2000-10-24 26.260
I2000-10-25 26.260
};
open HEAD,$head;
while () {
chomp;
($key,$value)= split /=/;
($value) = split /\//,$value;
$key=~s/\s//g;
$value=~s/\s//g;
$value=~s/\'//g;
if ($key eq $expkey ) {$exptime = $value}
if ($key eq $airkey ) {$airmass = $value}
if ($key eq $datekey) {$date = $value}
}
close HEAD;
$fact=$fact{$colour};
unless (defined($zpoint{$colour.$date})) {
print "ERROR: couldn't find zero-point\n";
print "$head $colour $date\n";
die;
}
$zpoint =$zpoint{$colour.$date} + 2.5*&log10($exptime) - $airmass*$fact;
}
Back to astrometry.
sex.pl >sex.log (did this before)
dofit.pl USNO 1 100 dofit1.log
checkdo.pl dofit1.log
dofit.pl USNO 2 100 dofit2.log
checkdo.pl dofit2.log
ls *.fits | perl -nle 'system "uphead $_ fead"'
dofit.pl USNO 2 100 dofit2.log
/raid/gwyn/F22/F22I #52 > checkdo.pl dofit2.log
worst centering errors:
551373o06oct001p06.fits 74 2 21 57 1.322
550366o10sep000p03.fits 67 2 50 64 0.730
550376o03sep000p05.fits 73 2 51 64 0.727
508267o09nov990p04.fits 73 2 52 65 0.723
550379o04sep000p05.fits 59 2 46 55 0.709
508268o00nov990p04.fits 72 2 45 54 0.709
fewest initial matches:
550366o05sep000p03.fits 49 2 42 48 0.574
550377o05sep000p05.fits 50 2 38 45 0.508
550376o05sep000p05.fits 50 2 38 42 0.498
550362o05sep000p03.fits 51 2 43 48 0.502
550375o05sep000p05.fits 51 2 38 45 0.566
550361o05sep000p03.fits 52 2 41 49 0.590
fewest final matches:
550376o01sep000p05.fits 58 2 32 34 0.568
550377o01sep000p05.fits 60 2 34 36 0.569
550375o01sep000p05.fits 60 2 36 38 0.536
552957o03oct003p12.fits 52 2 36 39 0.541
550376o05sep000p05.fits 50 2 38 42 0.498
550361o01sep000p03.fits 59 2 40 42 0.528
crossref rdall 2 8 .5 master
loaded 576169 objects from 859 files
sorting them... done
cross correlating, iteration: 1
found 449487 new matches
tossed 2350 bad objects
93916 reference objects
cross correlating, iteration: 2
found 1 new matches
tossed 10 bad objects
93906 reference objects
cross correlating, iteration: 3
found 0 new matches
tossed 0 bad objects
93906 reference objects
93906
dofit.pl master 2 500 dofit3.log
checkdo.pl dofit3.log
worst centering errors:
531632o02apr002p11.fits 433 2 321 414 0.163
531631o02apr002p11.fits 430 2 324 413 0.161
531630o02apr002p11.fits 381 2 298 367 0.159
531290o00apr002p06.fits 373 2 270 337 0.155
531288o00apr002p06.fits 364 2 257 323 0.151
531182o00apr002p06.fits 377 2 257 333 0.147
fewest initial matches:
531632o05apr002p11.fits 281 2 204 268 0.081
531632o04apr002p11.fits 282 2 201 265 0.048
531631o04apr002p11.fits 284 2 204 269 0.051
531631o05apr002p11.fits 285 2 207 270 0.080
531630o05apr002p11.fits 285 2 208 271 0.084
531630o04apr002p11.fits 288 2 210 275 0.055
fewest final matches:
531632o04apr002p11.fits 282 2 201 265 0.048
531289o10apr002p06.fits 296 2 230 267 0.068
531632o05apr002p11.fits 281 2 204 268 0.081
531631o04apr002p11.fits 284 2 204 269 0.051
531631o05apr002p11.fits 285 2 207 270 0.080
531630o05apr002p11.fits 285 2 208 271 0.084
crossref rdall 2 8 .5 master2
loaded 576169 objects from 859 files
sorting them... done
cross correlating, iteration: 1
found 449670 new matches
tossed 2300 bad objects
93932 reference objects
cross correlating, iteration: 2
found 2 new matches
tossed 4 bad objects
93931 reference objects
cross correlating, iteration: 3
found 0 new matches
tossed 0 bad objects
93931 reference objects
93931
dofit.pl master2 2 500 dofit4.log
checkdo.pl dofit4.log
worst centering errors:
531631o02apr002p11.fits 430 2 323 414 0.122
531630o02apr002p11.fits 381 2 300 367 0.119
531290o07apr002p06.fits 454 2 300 409 0.114
531288o00apr002p06.fits 368 2 258 324 0.113
531632o02apr002p11.fits 433 2 320 415 0.113
531182o05apr002p06.fits 399 2 275 348 0.112
fewest initial matches:
531632o05apr002p11.fits 281 2 204 267 0.070
531631o05apr002p11.fits 284 2 206 270 0.066
531632o04apr002p11.fits 284 2 201 267 0.047
531630o05apr002p11.fits 285 2 209 271 0.070
531631o04apr002p11.fits 286 2 204 269 0.047
531630o04apr002p11.fits 290 2 210 275 0.050
fewest final matches:
531632o05apr002p11.fits 281 2 204 267 0.070
531632o04apr002p11.fits 284 2 201 267 0.047
531289o10apr002p06.fits 297 2 231 268 0.060
531631o04apr002p11.fits 286 2 204 269 0.047
531631o05apr002p11.fits 284 2 206 270 0.066
531630o05apr002p11.fits 285 2 209 271 0.070
tar -cvf phot.tar *rd *fead (at least for the nonce)
crossref rdall 2 8 .5 master3
loaded 576169 objects from 859 files
sorting them... done
cross correlating, iteration: 1
found 449681 new matches
tossed 2294 bad objects
93932 reference objects
cross correlating, iteration: 2
found 4 new matches
tossed 2 bad objects
93930 reference objects
cross correlating, iteration: 3
found 0 new matches
tossed 0 bad objects
93930 reference objects
93930
dofit.pl master2 2 500 dofit5.log
worst centering errors:
531631o02apr002p11.fits 430 2 323 414 0.122
531630o02apr002p11.fits 381 2 300 367 0.119
531290o07apr002p06.fits 454 2 300 409 0.114
531288o00apr002p06.fits 368 2 258 324 0.113
531632o02apr002p11.fits 433 2 320 415 0.113
531182o05apr002p06.fits 399 2 275 348 0.112
fewest initial matches:
531632o05apr002p11.fits 281 2 204 267 0.070
531631o05apr002p11.fits 284 2 206 270 0.066
531632o04apr002p11.fits 284 2 201 267 0.047
531630o05apr002p11.fits 285 2 209 271 0.070
531631o04apr002p11.fits 286 2 204 269 0.047
531630o04apr002p11.fits 290 2 210 275 0.050
fewest final matches:
531632o05apr002p11.fits 281 2 204 267 0.070
531632o04apr002p11.fits 284 2 201 267 0.047
531289o10apr002p06.fits 297 2 231 268 0.060
531631o04apr002p11.fits 286 2 204 269 0.047
531631o05apr002p11.fits 284 2 206 270 0.066
531630o05apr002p11.fits 285 2 209 271 0.070
** 13 July 2002
The zero-point for Oct 00-3 is 26.260+/-0.006
bring *.fead and *.rd over to lasa28
/home/gwyn/red/F22/F22I/phot
Run erd.pl
erdall
Exposure.erd zero-point UTC-date Hawaii
calculated from (DATE-OBS) Standard
the .fead file Time
486370may990p03.erd 32.8002912379574 1999-05-13 03:42
486371may990p03.erd 32.8133312379574 1999-05-13 03:56
486372may990p03.erd 32.8225312379574 1999-05-13 04:09
486373may990p03.erd 32.8272512379574 1999-05-13 04:23
486374may990p03.erd 32.8334112379574 1999-05-13 04:36
508265nov990p04.erd 32.9521312379574 1999-11-09 19:10
etc.
run ~/pg/magmatch.f
(which now uses the median, instead of the average)
run ~/pg/connph.pl which
calculates zero-points for every image
based on every photometric image.
It outputs the average.
The graphical results are in:
photcalF22I.gif
photcalF22I.ps
which were outputed by sm.5
Lower left: layout of the 22hr pointings, outlined in red labelled in
black. Note the CCD06 of May 1999 has been excluded.
Upper left: layout of the 22hr pointings. "Photometric" exposures are
in heavy black and labelled "0". Exposures that "touch" (more than 150
common objects) photometric exposures are in red and labelled "1".
Exposures that "touch" level 1 exposures are outlined in blue and
labelled "2".
Big right figure: For every exposure, I calculated the zero-point you
would get if you assumed it was photometric, using the airmass and
exposure time from the header and the zero-point for the run from
Yannick's web page. (Except for Oct-00-3 where I used my own
zero-point)
I selected some exposures as really being photometric. These ones are
shown in black. For each "photometric" exposure, I derived zero-points
for all the other exposures, by comparing magnitudes in the overlap
regions. The difference between these magnitudes is plotted. The
exposures are labelled on the left and the UTC date and HST time are
shown on the right. The colour-coding is the same as the upper left
panel.
sm.6 produces magmatchdiff.ps
which shows the difference between
Average (3-sigma clipped)
Average (2-sigma clipped)
Median
Only the worst case scenarios (for Noverlaps>150) are shown
Not much of a difference really.
I checked a few bad cases. In every case, the median
is better. They are often cases where CCD's with
internal zero-point offsets are being compared.
The zero-point shifts internal to each pointing are
shown in interncheck.ps generated with sm.2
sm.1 is for comparing the results of matchex...
and has the macro comp in it.
There are night to night variations in the zero-points.
Look at ~/red/stan/*oct003*d.ps
SA92 2000-10-23 26.259
SA92 2000-10-24 26.238
SA95 2000-10-24 26.254 ->26.257
SA101 2000-10-24 26.280
SA101 2000-10-25 26.309
** 14 July 2002
photometry finalized
final astrometry checked using
lasa28:/home/gwyn/ag/rdmatch.f
which matches all the .rd files together
matching everything with r<2" and deltamag=.5
It runs on rdall (image zero-point)
mkrdall.pl makes rdall, working from
the .fead files.
Examine the worst case scenarios:
lasa28:/home/gwyn/red/F22/F22I/dofit5check/sm.1
which plots up various parameter from rdall.sum
name name #overlap #bad .95% fraction_bad
552856o05oct003p09.rd 486370o02may990p03.rd 303 10 0.232 0.033
486372o10may990p03.rd 486370o10may990p03.rd 251 8 0.225 0.032
486370o10may990p03.rd 486372o10may990p03.rd 250 8 0.236 0.032
551376o06oct001p06.rd 531182o00apr002p06.rd 373 12 0.255 0.032
551377o06oct001p06.rd 531182o00apr002p06.rd 362 17 0.285 0.047
550369o05sep000p04.rd 550373o05sep000p04.rd 384 13 0.228 0.034
550378o05sep000p05.rd 550375o05sep000p05.rd 473 15 0.181 0.032
531182o00apr002p06.rd 551376o06oct001p06.rd 380 12 0.255 0.032
bad here means worse than 0.3" (or 1.5 pixels)
and
lasa28:/home/gwyn/red/F22/F22I/dofit5check/sm.2
which shows the bad region
green <0.3"
blue <0.6"
red >0.6"
lasa28:/home/gwyn/ag/applyzp.pl
corrects to AB magnitude, then puts the zeropoint into the headers.
copy the headers back to pix1
(after doing some checking on the astrometry)
cd /raid/gwyn/F22/F22I/
mkdir swarp
mkdir stacks
cd swarp
mkdir vmem
mkdir resamp
doswarp.pl
** 15 July 2002
For future reference, swarp resamples all the images it can lay
it's hands on, regardless if they are relevant.
The resampled images come out 1x1 pixel...
The speed is not greatly affected.
SWarp fails (on pix1) if the center is specified
and there are more than 5 images.
** 16 July 2002
Successfully SWarped
*oct003p08 nov990p04 oct003p09
F22t2..
sex -c bigimage.sex F22It2.fits -CATALOG_NAME F22It2.cat -WEIGHT_IMAGE F22It2.weight.fits
swarping *sep000*
Saturation value is about 1000, not 40000
[1] doswarp.pl F22I.list &
/raid/gwyn/F22/F22I/swarp #39 >
> *Error*: cannot open for reading /raid/gwyn/F22/F22I/swarp/resamp//551382o09oct001p07.resamp.weight.fits
** 17 July 2002
F22It2.fits 0.830 square degrees
F22It3.fits 0.965 square degrees
sex -c bigimage.sex F22It3.fits -CATALOG_NAME F22It3.cat -WEIGHT_IMAGE F22It3.weight.fits
Redo SExtractor because I used the wrong threshold.
More data for P1, P2 taken on
20-21 September 2000
setting gain to 0.04
(that is the average of the FSCALES time the average gain of the
original images) 0.033 * 1.2
SATUR_LEVEL 1200. # level (in ADUs) at which arises saturation
bigimage.sex looks like this:
# Default configuration file for SExtractor V1.2b14 - > 2.0
# EB 23/07/98
# (*) indicates parameters which can be omitted from this config file.
#-------------------------------- Catalog -----------------------
CATALOG_NAME F22It3.cat # name of the output catalog
CATALOG_TYPE ASCII_HEAD # "NONE","ASCII_HEAD","ASCII","FITS_1.0"
# or "FITS_LDAC"
PARAMETERS_NAME header.param # name of the file
# containing catalog contents
#------------------------------- Extraction -------------------------
DETECT_TYPE CCD # "CCD" or "PHOTO" (*)
DETECT_MINAREA 3 # minimum number of pixels above threshold
DETECT_THRESH 1 # or , in mag.arcsec-2
ANALYSIS_THRESH 1 # or , in mag.arcsec-2
THRESH_TYPE RELATIVE
FILTER Y
FILTER_NAME gauss_3.0_7x7.conv # name of the file containing the filter
DEBLEND_NTHRESH 32 # Number of deblending sub-thresholds
DEBLEND_MINCONT 0.002 # Minimum contrast parameter for deblending
CLEAN Y
CLEAN_PARAM 1.0 # Cleaning efficiency
MASK_TYPE CORRECT # type of detection MASKing: can be one of
# "NONE", "BLANK" or "CORRECT"
#------------------------------ Photometry -----------------------------------
PHOT_APERTURES 15.,25. # MAG_APER aperture diameter(s) in pixels
PHOT_FLUXFRAC 0.2, 0.5, 0.8 # Fraction of FLUX_AUTO defining FLUX_RADIUS
PHOT_AUTOPARAMS 2.5, 3.5 # MAG_AUTO parameters: ,
SATUR_LEVEL 1000. # level (in ADUs) at which arises saturation
PHOT_AUTOAPERS 15.0, 15.0 # MAG_AUTO minimum apertures: estimation,photometry
MAG_ZEROPOINT 30.00 # magnitude zero-point
MAG_GAMMA 4.0 # gamma of emulsion (for photographic scans)
GAIN 0.04 # detector gain in e-/ADU.
PIXEL_SCALE 0. # size of pixel in arcsec (0=use FITS WCS info).
#------------------------- Star/Galaxy Separation -------------------
SEEING_FWHM 0.9 # stellar FWHM in arcsec
STARNNW_NAME default.nnw # Neural-Network_Weight table filename
#------------------------------ Background ---------------------------
BACK_SIZE 512 # Background mesh: or ,
BACK_FILTERSIZE 9 # Background filter: or ,
BACKPHOTO_TYPE LOCAL # can be "GLOBAL" or "LOCAL" (*)
BACKPHOTO_THICK 30 # thickness of the background LOCAL annulus (*)
BACK_TYPE AUTO
BACK_VALUE 0
#------------------------------ Check Image -----------------------
CHECKIMAGE_TYPE NONE # can be one of "NONE", "IDENTICAL",
# "BACKGROUND", "-BACKGROUND","BACKGROUND_RMS",
# "MINIBACKGROUND", "MINIBACK_RMS",
# "FILTERED", "OBJECTS", "-OBJECTS",
# "SEGMENTATION", or "APERTURES",
CHECKIMAGE_NAME check.fits # Filename for the check-image (*)
#--------------------- Memory (change with caution!) --------------------
MEMORY_OBJSTACK 2000 # number of objects in stack
MEMORY_PIXSTACK 5000000 # number of pixels in stack
MEMORY_BUFSIZE 512 # number of lines in buffer
#----------------------------- Miscellaneous --------------------------
VERBOSE_TYPE NORMAL # can be "QUIET", "NORMAL" or "FULL" (*)
#------------------------------- New Stuff ----------------------------
WEIGHT_TYPE MAP_WEIGHT
#WEIGHT_TYPE BACKGROUND
While header.param looks like this:
NUMBER
X_IMAGE
Y_IMAGE
ALPHA_J2000
DELTA_J2000
A_IMAGE
ERRA_IMAGE
B_IMAGE
ERRB_IMAGE
THETA_IMAGE
ERRTHETA_IMAGE
X2_IMAGE
ERRX2_IMAGE
Y2_IMAGE
ERRY2_IMAGE
XY_IMAGE
ERRXY_IMAGE
A_WORLD
ERRA_WORLD
B_WORLD
ERRB_WORLD
THETA_J2000
ERRTHETA_J2000
X2_WORLD
ERRX2_WORLD
Y2_WORLD
ERRY2_WORLD
XY_WORLD
ERRXY_WORLD
ISOAREA_IMAGE
THRESHOLD
MU_THRESHOLD
BACKGROUND
FLUX_MAX
MU_MAX
KRON_RADIUS
FLUX_RADIUS(3)
MAG_ISO
MAGERR_ISO
FLUX_ISO
FLUXERR_ISO
MAG_ISOCOR
MAGERR_ISOCOR
FLUX_ISOCOR
FLUXERR_ISOCOR
MAG_AUTO
MAGERR_AUTO
FLUX_AUTO
FLUXERR_AUTO
MAG_BEST
MAGERR_BEST
FLUX_BEST
FLUXERR_BEST
MAG_APER(2)
MAGERR_APER(2)
FLUX_APER(2)
FLUXERR_APER(2)
FLAGS
CLASS_STAR
Exciting stuff, eh?
The mapix job finished.
FTP'ing to lasa28
cut at: RA=334.531 and Dec= 0.409371
into four images
see:
/home/gwyn/red/F22/F22Icut4.ps
generated with /home/gwyn/red/F22/F22I/swarpcheck/sm.1
This works out to:
20006.784 3381.134 334.531 -.4
20006.713 12119.821 334.531 0.1
20006.584 27848.442 334.531 1
for x and
1988.075 17527.351 333.5 0.409371
19464.935 17526.346 334.5 0.409371
32572.268 17527.022 335.25 0.409371
for y
call it (20000,17500)
NAXIS1 = 37573 / Number of pixels along this axis
NAXIS2 = 34187 / Number of pixels along this axis
what we want is
imcopy F22I[1:20000,1:17500]
imcopy F22I[1:20000,17501:34187]
imcopy F22I[20001:37573,1:17500]
imcopy F22I[20001:37573,17501:34187]
From Seb:
voila comment on utilise toFITS2d:
/u3/hjmcc/bin/toFITS2d -x 9990 -X 15000 -y 3709 -Y 10000
F10I-center+P11_12_new2.fits >! F10I-center+P11_12_new2_sub.fits
ceci va produire une sous-image de 15000*10000 dont le coin bas a gauche a
pour coordonnees (9990,3709) dans l'ancienne image ...
cheers
seb
/u3/hjmcc/bin/toFITS2d -x 9990 -X 15000 -y 3709 -Y 10000
F10I-center+P11_12_new2.fits >! F10I-center+P11_12_new2_sub.fits
/u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 1 -Y 17500 F22I.fits >! F22Isw.fits
/u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 1 -Y 17500 F22I.fits >! F22Ise.fits
/u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 17500 -Y 16686 F22I.fits >! F22Inw.fits
/u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 17500 -Y 16686 F22I.fits >! F22Ine.fits
/u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 1 -Y 17500 F22I.weight.fits >! F22Isw.weight.fits
/u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 1 -Y 17500 F22I.weight.fits >! F22Ise.weight.fits
/u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 17500 -Y 16686 F22I.weight.fits >! F22Inw.weight.fits
/u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 17500 -Y 16686 F22I.weight.fits >! F22Ine.weight.fits
The north fiel
The north fields don't work.
There is a problem somewhere 72% into the image.
Actually, some of the image is missing.
** 19 July 2002
Working on the quadrants that work.
sex -c bigimage.sex F22Ise.fits -CATALOG_NAME F22Ise.cat -WEIGHT_IMAGE F22Ise.weight.fits
sex -c bigimage.sex F22Isw.fits -CATALOG_NAME F22Isw.cat -WEIGHT_IMAGE F22Isw.weight.fits
exsex F22Ise.cat X_IMAGE Y_IMAGE ALPHA_j2000 DELTA_j2000 MAG_AUTO MAGERR_AUTO 'flux_radius(2)' flags | grep -v 99.0000 > F22Ise.xyrdm
exsex F22Isw.cat X_IMAGE Y_IMAGE ALPHA_j2000 DELTA_j2000 MAG_AUTO MAGERR_AUTO 'flux_radius(2)' flags | grep -v 99.0000 > F22Isw.xyrdm
swarp (F22Itry2)
^[[1A^[[1M> Co-adding line: 29360 / 34187
^[[1A^[[1M> Co-adding line: 29376 / 34187
^[[1A^[[1M> Co-adding line: 29392 / 34187
^[[1A^[[1M> Co-adding line: 29408 / 34187
^[[1Afloating exception
Sigh.
The part that was written is the same as before
imstat F22Itry1/stacks/F22Isw.fits
# IMAGE NPIX MEAN MODE MIDPT MIN MAX
F22Itry1/stacks/F22Isw.fits 350000000 3.394 -1.397 -0.9151 -350.8 3163.
cl> e
imstat F22Itry2/stacks/F22Isw.fits
# IMAGE NPIX MEAN MODE MIDPT MIN MAX
F22Itry2/stacks/F22Isw.fits 350000000 3.394 -1.397 -0.9151 -350.8 3163.
cl>
Recompiling with -ieee at Henry's suggestion.
I don't think it will work.
Regions in DS9:
"phyiscal" does the correct thing on binned images
(
image;box(2807,3512.25,3,10.5,0) # color=red
^ ^ ^ ^ ^
| | | | |____ God knows
| | | |_________ length of the box sides
| | |___________ length of the box top and bottom
| |________________ center of box
|_____________________ center of box
** 20 July 2002
Looking for new images in the central part of the F22...
/evsix /evsix1/cfh12k_mai_1999/485930o:
/sampl /samplix1/DATA/cfh12k_oct2_2000/552568o:
485930 485930o/ bad
485931 485931o.phu bad
485932 485932o1.pl bad
485933 485933o/ bad
485934 485934o/ bad
486004 486004o.phu bad
486005 486005o.phu bad
486006 486006o.phu bad
486007 486007o.phu bad
486008 486008o.phu bad
486009 486009o.phu bad
486090 486090o.phu bad
486091 486091o.phu bad
486092 486092o.phu bad
486093 486093o.phu bad
486176 486176o.phu check it out no fits files
486177 486177o.phu check it out no fits files
486370 486370o10MDFGOfl.fits used
486372 486372o/ used
486373 486373o/ used
486374 486374o10MDFGOfl.fits used
486663 486663o/ B
486664 486664o/ B
486665 486665o/ B
486830 486830o/ B
486831 486831o/ B
486832 486832o/ B
486833 486833o/ B
508265 508265o00MDFG.fits.gz used
508267 508267o00MDFG.fits.gz used
508268 508268o00MDFG.fits.gz used
508269 508269o00MDFG.fits.gz used
508364 508364o10MDFOfl.fits B
508365 508365o00MDF.fits.gz B
508496 508496o00MDF.fits.gz B
508497 508497o10MDFOfl.fits.gz B
508498 508498o00MDF.fits.gz B
508594 508594o10MDFOfl.fits V
508595 508595o00MDF.fits.gz V
508596 508596o00MDF.fits.gz V
508597 508597o00MDF.fits.gz V
508670 508670o00MDF.fits.gz V
552568 552568o/ unreduced
552569 552569o09_hdr.fits unreduced
552570 552570o00.fits unreduced
552571 552571o.phu unreduced
552572 552572o.phu unreduced
552573 552573o00_hdr.fits unreduced
552574 552574o00.fits unreduced
552576 552576o.phu unreduced
552577 552577o/ unreduced
552578 552578o00_hdr.fits unreduced
552579 552579o00.fits unreduced
552632 552632o/ unreduced
552633 552633o/ unreduced
552634 552634o00_hdr.fits unreduced
552635 552635o00.fits unreduced
552637 552637o.phu unreduced
552638 552638o/ unreduced
552639 552639o00_hdr.fits unreduced
552640 552640o.phu unreduced
552641 552641o.phu unreduced
552642 552642o/ unreduced
552643 552643o00_hdr.fits unreduced unreduced
552644 552644o00.fits unreduced
552960 552960o11_hdr.fits unreduced
552962 552962o/ R
552963 552963o09_hdr.fits R
552964 552964o00.fits R
553047 553047f.phu unreduced
553049 553049o/ R
553050 553050d00_hdr.fits unreduced
553051 553051d00.fits unreduced
553052 553052d.phu unreduced
553053 553053d/ unreduced
553054 553054d00_hdr.fits unreduced
553055 553055d.phu unreduced
587773 587773~2.jpg R
594985 594985~1.jpg no fits files
603105 603105~1.hdr V
** 21 July 2002
The latest stack worked! Maybe there is something to this
-ieee stuff after all. I doubt it.
/u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 1 -Y 17500 F22I.fits >! F22Isw.fits
/u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 1 -Y 17500 F22I.fits >! F22Ise.fits
/u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 17500 -Y 16686 F22I.fits >! F22Inw.fits
/u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 17500 -Y 16686 F22I.fits >! F22Ine.fits
/u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 1 -Y 17500 F22I.weight.fits >! F22Isw.weight.fits
/u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 1 -Y 17500 F22I.weight.fits >! F22Ise.weight.fits
/u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 17500 -Y 16686 F22I.weight.fits >! F22Inw.weight.fits
/u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 17500 -Y 16686 F22I.weight.fits >! F22Ine.weight.fits
boxes:
F22Iseall.reg all the boxes, plus doubles, plus targets
F22Iseb4.reg don't know
F22Isebox.reg just the boxes, saturated stars and bad regions
F22Isestars.reg early star boxes
fred.reg
USNO.reg USNO catalog
tvm.pl does tvmark like stuff.
Running imstat to see if everything is hunky-dory
IMAGE NPIX MEAN MODE MIDPT MIN MAX
F22Isw.fits 350000000 3.394 -1.397 -0.9151 -350.8 3163.
It is. Compare with above.
SExtraction is happening in
pix1:/raid/gwyn/F22/F22Imapix
The catalogs from try1 (which failed 72% through
the final write) are in pix1:/raid/gwyn/F22/F22Imapix/try1
The good stuff is in pix1:/raid/gwyn/F22/F22Imapix
sex -c bigimage.sex F22Ise.fits -CATALOG_NAME F22Ise.cat -WEIGHT_IMAGE F22Ise.weight.fits
sex -c bigimage.sex F22Ine.fits -CATALOG_NAME F22Ine.cat -WEIGHT_IMAGE F22Ine.weight.fits
sex -c bigimage.sex F22Inw.fits -CATALOG_NAME F22Inw.cat -WEIGHT_IMAGE F22Inw.weight.fits
sex -c bigimage.sex F22Isw.fits -CATALOG_NAME F22Isw.cat -WEIGHT_IMAGE F22Isw.weight.fits
exsex F22Ise.cat X_IMAGE Y_IMAGE ALPHA_j2000 DELTA_j2000 MAG_AUTO MAGERR_AUTO 'flux_radius(2)' flags | grep -v 99.0000 > F22Ise.xyrdm
exsex F22Isw.cat X_IMAGE Y_IMAGE ALPHA_j2000 DELTA_j2000 MAG_AUTO MAGERR_AUTO 'flux_radius(2)' flags | grep -v 99.0000 > F22Isw.xyrdm
exsex F22Ine.cat X_IMAGE Y_IMAGE ALPHA_j2000 DELTA_j2000 MAG_AUTO MAGERR_AUTO 'flux_radius(2)' flags | grep -v 99.0000 > F22Ine.xyrdm
exsex F22Inw.cat X_IMAGE Y_IMAGE ALPHA_j2000 DELTA_j2000 MAG_AUTO MAGERR_AUTO 'flux_radius(2)' flags | grep -v 99.0000 > F22Inw.xyrdm
Can't display F22I.fits on mapix
** 3 August 2002
Updated F22I.imlist in /samplix1/gwyn/track
to include the P13 images.
** 9 August 2002
Rexamine areas:
22:21:31 -00:32:00
** 13 August 2002
(working from Victoria)
updated F22I.dirlist: first line has P13...
/samplix1/foucaud/Oct2000_2/Results/f22/
/samplix1/dantel/cfh12k_sept_oct_2000/F22P3b-I-Sep00/notcoadded
/samplix3/mellier/F22P3-I-May99/notcoadded
/samplix1/dantel/cfh12k_sept_oct_2000/F22P4b-I-Sep00/notcoadded
/samplix3/mellier/F22P4-I-Nov99/notcoadded
/samplix1/dantel/cfh12k_sept_oct_2000/F22P5b-I-Sep00/notcoadded
/samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p6b_I/notcoadded
/samplix1/dantel/cfh12k_sept_oct_2000/F22P6b-I-Oct00/notcoadded
/samplix1/dantel/cfh12k_sept_oct_2000/F22P7b-I-Oct00/notcoadded
/samplix1/foucaud/Oct2000_3/Results/f22/f22p8b_I
/samplix1/foucaud/Oct2000_3/Results/f22/f22p9b_I
/samplix1/foucaud/Oct2000_3/Results/f22/f22p10b_I
/samplix1/foucaud/Oct2000_3/Results/f22/f22p11b_I
/samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p11b_I/notcoadded
/samplix1/foucaud/Oct2000_3/Results/f22/f22p12b_I
/samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p6b_I/notcoadded
/samplix1/dantel/DATAprocessed-avril2000-2.DEF/f22p11b_I/notcoadded
Tried to run linkall.pl
didn't work because samplix1 isn't online
modified linkall.pl to do just P13
(I have now modified it back)
On pix1:
moved /raid/gwyn/F22/F22I to /raid/gwyn/F22/F22Iv1
working in F22I
mkdir data
cd data
cp /raid/gwyn/F22/F22Iv1/data/*.fits .
ncftp -u gwyn samplix.iap.fr
cd /samplix1/gwyn/F22/F22IP13
get *.fits
bye
cp /raid/gwyn/F22/F22Iv1/data/add* .
cp add00oct001 add00oct002
cp add01oct001 add01oct002
cp add02oct001 add02oct002
cp add03oct001 add03oct002
cp add04oct001 add04oct002
cp add05oct001 add05oct002
cp add06oct001 add06oct002
cp add07oct001 add07oct002
cp add08oct001 add08oct002
cp add09oct001 add09oct002
cp add10oct001 add10oct002
cp add11oct001 add11oct002
imlist.pl
stackall >stackall.log
ls *N.fits | perl -nle '($a=$_)=~s/N.fits/.fits/; rename $_,"../$a"'
cp mask*box ..
cd ..
mkdir toss
mv *o06may*.fits toss
mv 531633o??apr002p11.fits toss
imlist.pl
sex.pl>sex.log
(get the USNO catalog...)
** 13 August 2002
Using 26.26 for the oct002 zero-point:
my %zpoint = qw {
I1999-05-13 25.802
I1999-11-09 25.894
I2000-04-30 26.3386
I2000-05-01 26.3386
I2000-05-03 26.3386
I2000-09-22 26.2906
I2000-10-03 26.2906
I2000-10-20 26.26
I2000-10-21 26.26
I2000-10-23 26.260
I2000-10-24 26.260
I2000-10-25 26.260
};
dofit.pl USNO 1 100 dofit1.log
checkdo.pl dofit1.log
dofit.pl USNO 2 100 dofit2.log
problems with some images:
Flatfielding problems:
look at:
552632o10oct002p13.fits
(for example)
** 22 August 2002
On pix1:
cd /raid/gwyn/F22/F22I/data/
imcopy the [1:2048,1:4096] section of all the images
to backoct002/
(use a ls *.fits | awk '{print imcopy $1....}' >fred.cl
style of thing
in Linux:
cd backoct002
back.pl
tar -cvf mb.tar *mb.fits
switch to lasa28:
cd /home/gwyn/red/F22/F22I/backcheck/mbdir
ncftp -u virmos pix1.iap.fr
cd /raid/gwyn/F22/F22I/data/backoct002
get mb.tar
quit
tar -xvf mb.tar
ls *mb.fits | perl -nle 'system "backch $_"' > backch.res
results:
Yup. They're wrong.
examine
/home/gwyn/red/F22/F22I/backcheck/mbdir/backch.res
(morgan or lasa28)
Things like:
552642o01oct002p13mb.fits 23635. 15.26 0.53 13.67 11.56 12.23 bad
552641o01oct002p13mb.fits 23103. 15.46 0.27 13.75 11.70 12.27 bad
552643o01oct002p13mb.fits 24005. 15.60 0.53 13.75 11.54 12.49 bad
552644o01oct002p13mb.fits 24335. 16.62 0.88 14.39 11.64 13.21 bad
552641o06oct002p13mb.fits 24123. 21.48 9.25 1.09 2.86 8.36 bad
552642o06oct002p13mb.fits 24665. 20.81 9.25 1.25 2.70 8.44 bad
552644o06oct002p13mb.fits 25684. 21.51 9.29 1.66 2.15 9.15 bad
= no good at all
** 23 August 2002
/home/gwyn/mask (on Morgan) will contain a bunch of
programs for dealing with the masks.
Or at least one:
dedub.pl will take F22Iseall.reg
and F22Ise.xyrdm
and make fred.reg which has all the boxes of
F22Iseall.reg (no doubles) and all the usuable
stars of F22Ise.xyrdm (flag<4 and mag<22.5)
as red boxes and blue circular points
F22Isefin.reg is the definitive version.
** 24 August 2002
F22I??.xyrdm (x,y,RA,Dec,Mag_auto,Magerr_auto,flux_radius)
fluxradius in pixels...
The general masking procedure
- make some boxes with ds9 regions
- also make them using ~/mask/mksat.pl
- save them to F22Iswall.reg
- merge them with the relevant object list fro F22Isw.xyrdm
using ~/mask/dedub.pl. Output is fred.reg
- plot them again with ds9
- add more boxes to fred.reg using a distinctive colour
- add these new boxes (by hand) to F22Iswall.reg
** 25 August 2002
Rexamine areas:
22:21:31 -00:32:00
22:15:37 +00:11:53 (sw)
Big circle-ish area...
exsex F22Isw.cat X_IMAGE Y_IMAGE ALPHA_j2000 DELTA_j2000 MAG_APER MAGERR_APER 'flux_radius(2)' flags | grep -v 99.0000 > F22Isw.xyrdma
exsex F22Isw.cat ALPHA_j2000 DELTA_j2000 MAG_APER MAG_AUTO >F22Isw.xymm
There are a number of problems....
Which will be full resolved at a later date.
~/mask/dedub.f takes over from ~/mask/dedub.pl
It will also take the "lines" which it blankets with
boxes
grep cyan fred.reg > fred ; cat F22Ineall.reg fred >fred2 ; \mv fred2 F22Ineall.reg ; dedub
** 26 August 2002
awk '$5<22.5 && $8<4 {print $3,$4,$5,$6}' F22Isw.xyrdm >fredsw
awk '$5<22.5 && $8<4 {print $3,$4,$5,$6}' F22Ise.xyrdm >fredse
awk '$5<22.5 && $8<4 {print $3,$4,$5,$6}' F22Inw.xyrdm >frednw
awk '$5<22.5 && $8<4 {print $3,$4,$5,$6}' F22Ine.xyrdm >fredne
cat fredsw fredse frednw fredne > F22I.rdm
conmask.pl takes the F22I??all.reg files and converts them
from DS9 format physical coordinates to pairs of RA,Dec (2000).
It uses the F22I??.head and xy2skypv to do the transforming.
Using the F22I??b4.fits headers (and multiplying the coordinates
by 4) doesn't work as well.
Compare using ds9 to do the transformation (x,y)->(RA,Dec)
and chmask.pl to make convert from DS9 format to corners
appmask takes up to three aruguments:
box list (assumes ra1 dec1 ra2 dec2 in decimal degrees)
inital object list (ra,dec in columns 3,4)
output object list (same format as inital, but only
the objects outside the boxes)
appmask F22Ise.box F22Ise.xyrdm
appmask F22Isw.box F22Isw.xyrdm
appmask F22Ine.box F22Ine.xyrdm
appmask F22Inw.box F22Inw.xyrdm
appmask F22I.box F22I.xyrdm
(Should be obvious at this point that the .box files contain the masks...)
There are patches of badness, but what can you do?
New data from Seb
/samplix1/foucaud/Oct2000_2/Results/f22_2
Now linked correctly in:
/samplix1/gwyn/F22/F22IP13
on pix1:
cd /raid/gwyn/F22/F22I/data/
mkdir backoct0002bis
cd backoct0002bis
ncftp -u gwyn samplix.iap.fr
cd /samplix1/gwyn/F22/F22IP13
get *
bye
ls *.fits | awk '{print "imcopy " $1 "[1:2048,1:4096] " $1}' >fred.cl
(in IRAF)
cl backch.res
** 27 August 2002
Decision:
include
552573o??oct002p13
552574o??oct002p13
552576o??oct002p13
552577o??oct002p13
552578o??oct002p13
552579o??oct002p13
except for 00
Others have large flat field gradients...
ncftp -u gwyn samplix
cd /samplix1/gwyn/F22/F22IP13
get 552573*.fits
get 552574*.fits
get 552576*.fits
get 552577*.fits
get 552578*.fits
get 552579*.fits
bye
imlist.pl
stackall >stackall.log
ls *N.fits | perl -nle '($a=$_)=~s/N.fits/.fits/; rename $_,"../$a"'
cp mask*box ..
cd ..
mkdir toss
imlist.pl
sex.pl>sex.log
mv *o06may*.* toss
mv 531633o??apr002p11.* toss
mv 55257?o00oct002p13.* toss
** 28 August 2002
dofit.pl USNO 1 100 dofit1.log
checkdo.pl dofit1.log
dofit.pl USNO 2 100 dofit2.log
checkdo.pl dofit2.log
Looks good...
ls *.fits | perl -nle 'system "uphead $_ fead"'
dofit.pl USNO 2 100 dofit2.log
checkdo.pl dofit2.log
crossref
(checking...
megamatch rdall
get rdall.sum
plot number of overlapping galaxies vs. relative centering error
do this on Morgan: /home/gwyn/red/F22/F22I/dofit6check
sm.1
)
dofit.pl master 2 500 dofit3.log
checkdo.pl dofit3.log
Woah!
checkdo.pl dofit3.log
worst centering errors:
552576o10oct002p13.fits 421 2 56 158 1.446
531632o02apr002p11.fits 434 2 322 415 0.166
531631o02apr002p11.fits 431 2 325 414 0.164
531630o02apr002p11.fits 382 2 299 368 0.160
531182o05apr002p06.fits 399 2 275 349 0.144
531183o07apr002p06.fits 456 2 311 413 0.143
fewest initial matches:
531632o04apr002p11.fits 285 2 203 269 0.059
531632o05apr002p11.fits 287 2 208 273 0.090
531631o04apr002p11.fits 290 2 211 275 0.073
531631o05apr002p11.fits 294 2 214 277 0.093
531630o04apr002p11.fits 294 2 216 280 0.072
531630o05apr002p11.fits 297 2 218 281 0.094
fewest final matches:
552576o10oct002p13.fits 421 2 56 158 1.446
531632o04apr002p11.fits 285 2 203 269 0.059
531632o05apr002p11.fits 287 2 208 273 0.090
550361o05sep000p03.fits 316 2 203 274 0.080
531631o04apr002p11.fits 290 2 211 275 0.073
To fix:
Examine image: looks OK
plot detections: also looks OK
dofit 552576o10oct002p13.fits .cex master fead 2 y 500
WCS:
CRVALs 334.4591667 0.4119167
CRPIXs -2161.23 4218.18
CD: 5.667E-05 -7.500E-07 1.177E-06 5.676E-05
order: 2
PV1_: -1.15E-03 1.01E+00 -4.45E-03 0.00E+00 -2.18E-02 2.16E-02 -1.26E-03
PV2_: 1.16E-03 1.01E+00 -4.27E-03 0.00E+00 3.58E-02 -2.03E-02 3.23E-03
read 500 lines from 552576o10oct002p13.cex
read 96044 lines from master
retreived 497 stars from memory
In tol: 4.0 sec Out tol: 0.0 sec found 421 matches from 500 to 497
Doing linear part
In tol: 4.0 sec Out tol: 8.0 sec found 315 matches from 500 to 497
retreived 495 stars from memory
In tol: 4.0 sec Out tol: 8.0 sec found 315 matches from 500 to 495
In tol: 4.0 sec Out tol: 8.0 sec found 315 matches from 500 to 495
In tol: 4.0 sec Out tol: 8.0 sec found 315 matches from 500 to 495
Doing 2nd order
In tol: 4.0 sec Out tol: 8.0 sec found 316 matches from 500 to 495
In tol: 1.0 sec Out tol: 8.0 sec found 313 matches from 500 to 495
In tol: 1.0 sec Out tol: 8.0 sec found 312 matches from 500 to 495
In tol: 1.0 sec Out tol: 8.0 sec found 56 matches from 500 to 495
In tol: 2.0 sec Out tol: 4.0 sec found 158 matches from 500 to 495
552576o10oct002p13.fits 421 2 56 158 1.446
WCS:
CRVALs 334.4591667 0.4119167
CRPIXs -2161.06 4215.46
CD: 5.666E-05 -7.827E-07 1.150E-06 5.677E-05
order: 2
PV1_: -2.91E-03 1.01E+00 -3.83E-02 0.00E+00 -2.78E-02 1.65E-02 -1.38E-01
PV2_: -6.78E-04 9.79E-01 -2.75E-03 0.00E+00 -9.89E-02 -2.37E-02 -2.59E-04
Probably too many .cex detections
Fix it with:
dofit 552576o10oct002p13.fits .cex master fead 2 n 400
552576o10oct002p13.fits 332 2 257 318 0.060
We now resume our regualar programming:
crossref
dofit.pl master 2 500 dofit4.log
checkdo.pl dofit4.log
megamatch
(check...)
crossref
dofit.pl master 2 500 dofit5.log
checkdo.pl dofit5.log
megamatch
(checking...
megamatch rdall
get rdall.sum
plot number of overlapping galaxies vs. relative centering error
do this on Morgan: /home/gwyn/red/F22/F22I/dofit6check
sm.1
)
** 30 August 2002
Briefly re-ran
sex.pl
dofit.pl master 2 500 dofit6.log
with PHOT_APER 15
Then I switched it back again.
Having problems matching results from before.
Magnitudes have larger errors it seems.
Switching the magnitude threshold lower...
DETECT_MINAREA 5 # minimum number of pixels above threshold
DETECT_THRESH 5 # or , in mag.arcsec-2
ANALYSIS_THRESH 5 # or , in mag.arcsec-2
becomes
DETECT_MINAREA 3 # minimum number of pixels above threshold
DETECT_THRESH 3 # or , in mag.arcsec-2
ANALYSIS_THRESH 3 # or , in mag.arcsec-2
** 3 September 2002
To do:
From Henry:
- Figure out flags
- done
- square edges
- yes
- polygon code
- done
From Vincent:
- delivery mechanism: Put it somewhere and tell him
- images: F22I.fits and the subrasters
- catalogs: with the # thingies
- Henry flags: Two columns: ID, flag
- Stephen flags
- putting stuff on vcom: I have an account and /vcom/gwyn
From Stephen:
- mask edges: DONE
- document everything:
- star/galaxy separation
- bright star masking
- surface area
- new flags
- check catalog parameters parameters: Done
Bottom: -00:31:45 =~ 1122 -0.529166666
Top: +01:18:00 =~ 33095 1.3
West: 22:13:48 =~ 1115 333.45
East: 22:21:52 =~ 36360 335.46666666
Which becomes:
# Region file format: DS9 version 3.0
global color=red font="helvetica 10 normal" select=1 edit=1 move=1 delete=1 include=1 fixed=0 source
fk5;box(334.4583333,0.38541666,2.0166666,1.8291666666,-0)
cat>corners
333.45 -0.529166666
335.46666666 -0.529166666
335.46666666 1.3
333.45 1.3
/home/gwyn/red/F22/F22I/swarpcheck #8 > sky2xypv F22I.head corners
1112.108 1121.920 333.45 -0.529166666
36361.445 1121.937 335.46666666 -0.529166666
36357.678 33095.594 335.46666666 1.3
1115.896 33095.634 333.45 1.3
seeing.f
** 4 September 2002
Dicking with the boxes, just to be sure...
on Morgan
/home/gwyn/red/F22/F22I/swarpcheck #17 > conmask.pl
se
sw
ne
nw
cat F22I??.box > F22I.box
Had another look at F22Iseb4.fits
New and improved version of F22Iseall.fits
/home/gwyn/red/F22/F22I/swarpcheck/sm.o
produces F22Iboxes.gif
black boxes should be retained.
blue boxes are outside the boundaries and should be ditched.
Two more boxes to mask the seams...
RA : 334:31:40 +/- 5"
Dec: 0:24:28 +/- 5"
334.5291666666 334:31:45 22:18:07.00
334.5347222222 333:31:55 22:17:08.33
0:24:23.0 0.4063888888
0:24:33.0 0.4091666666
The boxes are:
F22Icr.box
333.45 0.4063888888 335.46666666 0.4091666666
334.5291666666 -0.529166666 334.5347222222 1.3
** OK, we were extracting the wrong -ing half light radius with exsex:
This has been fixed on Morgan
exsex F22Ise.cat X_IMAGE Y_IMAGE ALPHA_j2000 DELTA_j2000 MAG_AUTO MAGERR_AUTO 'flux_radius(2)' flags | grep -v 99.0000 > F22Ise.xyrdm
exsex F22Isw.cat X_IMAGE Y_IMAGE ALPHA_j2000 DELTA_j2000 MAG_AUTO MAGERR_AUTO 'flux_radius(2)' flags | grep -v 99.0000 > F22Isw.xyrdm
exsex F22Ine.cat X_IMAGE Y_IMAGE ALPHA_j2000 DELTA_j2000 MAG_AUTO MAGERR_AUTO 'flux_radius(2)' flags | grep -v 99.0000 > F22Ine.xyrdm
exsex F22Inw.cat X_IMAGE Y_IMAGE ALPHA_j2000 DELTA_j2000 MAG_AUTO MAGERR_AUTO 'flux_radius(2)' flags | grep -v 99.0000 > F22Inw.xyrdm
cat F22I??.xyrdm >F22I.xyrdm
appmask F22Ise.box F22Ise.xyrdm
appmask F22Isw.box F22Isw.xyrdm
appmask F22Ine.box F22Ine.xyrdm
appmask F22Inw.box F22Inw.xyrdm
appmask F22I.box F22I.xyrdm
The area of the F22I is:
3.349333 (dx=.01 degrees)
3.248901 (dx=.001 degrees)
** 9 September 2002
Belated documenation of star-galaxy separation
/home/gwyn/mask/seeing (on Morgan)
/home/gwyn/red/F22/F22I/swarpcheck #20 > seeing
read 107540 lines fromF22I.xyrdmcl
333.440521 335.479919
-0.543314993 1.31409574
nailing done the last squares
1.70000005 3.20000005
0 1.70000005 1.88750005
1 1.88750005 2.07500005
2 2.07500005 2.26250005
3 2.26250005 2.45000005
4 2.45000005 2.63750005
5 2.63750005 2.82500005
6 2.82500005 3.01250005
7 3.01250005 3.20000005
8 3.20000005 3.38750005
9 3.38750005 3.57500005
** 10 September 2002
Using lasa28:/home/gwyn/mask/area2.f
We get:
3.3426 sq deg at a resolution of 20000
3.3382 sq deg at a resolution of 10000
3.3293 sq deg at a resolution of 5000
3.3110 sq deg at a resolution of 2500
2500 3.3110
5000 3.3293
10000 3.3382
15000 3.3408
20000 3.3426
Running maskfits.f
(vcom:/vcom/gwyn/mask or lasa28:/home/gwyn/mask)
F22Ine.fits
WCS:
CRVALs 334.4612061 0.3846313
CRPIXs -1213.00 -406.00
CD: 5.722E-05 0.000E+00 0.000E+00 5.722E-05
reading F22Ine.fits done
oneing array
done
reading F22I.box done read 773 boxes
-18887.892-16378.080 16357.678 15595.594
doing right
doing top
writing F22IneMask.fits done
area (sq deg): 0.7365
F22Ise.fits
WCS:
CRVALs 334.4612061 0.3846313
CRPIXs -1213.00 17093.00
CD: 5.722E-05 0.000E+00 0.000E+00 5.722E-05
reading F22Ise.fits done
oneing array
done
reading F22I.box done read 468 boxes
-18887.892 1120.920 16357.678 33094.594
doing right
doing bottom
writing F22IseMask.fits done
area (sq deg): 0.7834
F22Isw.fits
WCS:
CRVALs 334.4612061 0.3846313
CRPIXs 18786.00 17093.00
CD: 5.722E-05 0.000E+00 0.000E+00 5.722E-05
reading F22Isw.fits done
oneing array
done
reading F22I.box done read 1205 boxes
1111.108 1120.920 36356.678 33094.594
doing left
doing bottom
writing F22IswMask.fits done
area (sq deg): 0.9182
F22Inw.fits
WCS:
CRVALs 334.4612061 0.3846313
CRPIXs 18786.00 -406.00
CD: 5.722E-05 0.000E+00 0.000E+00 5.722E-05
reading F22Inw.fits done
oneing array
done
reading F22I.box done read 1031 boxes
1111.108-16378.080 36356.678 15595.594
doing left
doing top
writing F22InwMask.fits done
area (sq deg): 0.9071
F22Ine area (sq deg): 0.7365
F22Ise area (sq deg): 0.7834
F22Isw area (sq deg): 0.9182
F22Inw area (sq deg): 0.9071
3.3452 sq deg total (adopted)
** 11 September 2002
Belated description of software...
area.f calculates the area by making small steps
in RA,Dec and counting the ones that that are
outside all the masks. It goes through every point
and checks to see if it is in every box. SLOW
ares2.f (faster) sets up a big grid, sets them all to 1
goes through every box and sets all points inside
the box to 0.
maskfits.f reads the 4 quartered .fits images, and makes
a mask image which is 0 inside a mask (or outside
the borders) and 1 outside a mask and inside the
borders. It also calculates the surface area.
The princples of star galaxy separation are documented at
http://www.astrsp-mrs.fr/gwyn/cat/starsep.html
seeing.f:(morgan:/home/gwyn/mask/seeing.f)
Given a file with ra,dec,magnitude,half-light-radius
Figures out typical stellar half-light-radius
by finding the peak of the half-light-radius
distribution in a given magnitude range.
Figures out the width of the locus by following
it down until falls below 20% of the peak of
the histogram.
Run it with ns=40 (40x40 boxes) to get a high resolution seeing map
of the F22I field.
Run it with ns=20 to get good star-galaxy separation
You need more counts in your histogram to get the
drop off than you do to get the peak
reads:
F22I.xyrdmcl
writes:
radecseeing (Ra,Dec,mode of half-light radius)
seeing.fort (binary file for sm.j)
histout (histogram (x,y,seeingblock), for reading by sm.n)
maskout (mag, half-light, useable, seeing blck, for reading by sm.n)
F22I.zones (ra1,ra2,dec1,dec2,hmax,hon,hoff)
objects with positions in the box (ra1:ra2,dec1:dec2)
with magntiudes in the range 19erdmall
#magmatch erdmall
erd.pl
magmatch erdall
ls *o09*fead | perl -nle 'system "getoutex $_"'
ls *exp >expall
sort -n *exp | head -1 >fred1
sort -n *exp | tail -1 >fred2
sort -n +1 *exp | head -1 >fred3
sort -n +1 *exp | head -1 >fred4
cat fred1 fred2 fred3 fred4 > exprel
connph.pl
(works now)
sm.5
(works)
/raid/gwyn/F22/F22I #21 > dofit.pl master 2 500 dofit7.log
update dofit on pix1 so it puts in the AstroGwyn header cards.
and on pix1:/raid/gwyn/F22/F22I
dofit.pl master 2 500 dofit7.log
tar -cvf feadrd.tar *.fead *.rd
gzip feadrd.tar
on lasa28:/home/gwyn/pg/F22I/
ncftp -u gwyn pix1.iap.fr
cd /raid/gwyn/F22/F22I
get feadrd.tar.gz
bye
gzip -dc feadrd.tar.gz | tar -xvf -
erd.pl
magmatch erdall (mag range 1850
Connectivity level
Photometric image 1 2 3 4 5 6 7 8 9
531182apr002p06 29 19 11 3 5 10
531183apr002p06 30 18 11 3 5 10
552850oct003p10 9 8 7 16 16 15 6
552851oct003p10 9 8 7 16 16 15 6
552852oct003p10 9 8 7 16 16 15 6
552853oct003p10 9 8 7 16 16 15 6
552854oct003p10 9 8 7 16 16 15 6
552855oct003p08 22 24 21 10
552856oct003p09 19 27 25 6
552858oct003p11 17 7 16 16 15 6
552859oct003p11 16 8 16 16 15 6
552860oct003p11 17 7 16 16 15 6
18150
Doesn't work
1860
Connectivity level
Photometric image 1 2 3 4 5 6 7 8 9
531182apr002p06 26 18 15 3 5 10
531183apr002p06 21 23 15 3 5 10
552850oct003p10 9 8 5 12 14 20 9
552851oct003p10 9 8 5 12 14 20 9
552852oct003p10 9 8 5 12 14 20 9
552853oct003p10 9 8 5 12 14 20 9
552854oct003p10 9 8 5 12 14 20 9
552855oct003p08 21 24 15 7 10
552856oct003p09 16 20 20 21
552858oct003p11 17 5 12 14 20 9
552859oct003p11 16 6 12 14 20 9
552860oct003p11 17 5 12 14 20 9
1850
Connectivity level
Photometric image 1 2 3 4 5 6 7 8 9
531182apr002p06 29 19 11 3 5 10
531183apr002p06 30 18 11 3 5 10
552850oct003p10 9 8 7 16 16 15 6
552851oct003p10 9 8 7 16 16 15 6
552852oct003p10 9 8 7 16 16 15 6
552853oct003p10 9 8 7 16 16 15 6
552854oct003p10 9 8 7 16 16 15 6
552855oct003p08 22 24 21 10
552856oct003p09 19 27 25 6
552858oct003p11 17 7 16 16 15 6
552859oct003p11 16 8 16 16 15 6
552860oct003p11 17 7 16 16 15 6
** 2 October 2002
Forget the above. In fact, there was a bug in
magmatch.f (which has been fixed)
connph.pl
sm.5 (just to check)
examine ~/log/F22Ip13zp.gif
applyzp.pl
ls *.fead | perl -nle '($a=$_)=~s/fead/head/; system "cp $_ $a"'
tar -cvf head.tar *head
move to pix1
cd /raid/gwyn/F22/F22I
tar -xvf head.tar
final astrometry checked using
lasa28:/home/gwyn/ag/rdmatch.f
which matches all the .rd files together
matching everything with r<2" and deltamag=.5
It runs on rdall (image zero-point)
mkrdall.pl makes rdall, working from
the .fead files.
** 3 October 2002
rerun crossref and dofit one more time...
on pix1:/raid/gwyn/F22/F22I
crossref
dofit.pl master 2 500 dofitX.log
mkdir rdcheck
cd rdcheck/
cp ../*.fead .
cp ../*.cex .
ls *.cex | perl -nle '($r=$_)=~s/cex/rd/; ($f=$_)=~s/cex/fead/; system "xy2skypv $f $_ $r"'
cp ../rdall .
rdmatch rdall
awk '$3>200 && $6>0.03' rdall.sum
/raid/gwyn/F22/F22I #46 > awk '$3>200 && $6>0.03' rdall.sum
552856o05oct003p09.rd 486371o02may990p03.rd 354 11 0.229 0.031
552574o03oct002p13.rd 552573o03oct002p13.rd 413 16 0.224 0.039
552576o03oct002p13.rd 552573o03oct002p13.rd 410 14 0.229 0.034
552577o03oct002p13.rd 552573o03oct002p13.rd 379 18 0.249 0.047
552578o03oct002p13.rd 552573o03oct002p13.rd 391 14 0.228 0.036
552579o03oct002p13.rd 552573o03oct002p13.rd 414 17 0.269 0.041
551384o05oct001p07.rd 551379o05oct001p07.rd 317 10 0.136 0.032
552576o05oct002p13.rd 552573o05oct002p13.rd 247 10 0.226 0.040
552577o05oct002p13.rd 552573o05oct002p13.rd 230 8 0.225 0.035
552578o05oct002p13.rd 552573o05oct002p13.rd 265 12 0.249 0.045
552579o05oct002p13.rd 552573o05oct002p13.rd 257 11 0.238 0.043
552576o05oct002p13.rd 552574o05oct002p13.rd 257 8 0.224 0.031
552573o05oct002p13.rd 552576o05oct002p13.rd 247 10 0.226 0.040
552574o05oct002p13.rd 552576o05oct002p13.rd 255 8 0.223 0.031
552578o05oct002p13.rd 552576o05oct002p13.rd 267 10 0.239 0.037
552578o05oct002p13.rd 552577o05oct002p13.rd 259 8 0.226 0.031
552573o05oct002p13.rd 552578o05oct002p13.rd 265 12 0.249 0.045
552576o05oct002p13.rd 552578o05oct002p13.rd 269 9 0.237 0.033
552577o05oct002p13.rd 552578o05oct002p13.rd 261 9 0.227 0.034
552573o05oct002p13.rd 552579o05oct002p13.rd 254 9 0.229 0.035
486370o02may990p03.rd 552805o05oct003p09.rd 276 13 0.250 0.047
486371o02may990p03.rd 552805o05oct003p09.rd 321 14 0.259 0.044
486373o02may990p03.rd 552805o05oct003p09.rd 323 10 0.221 0.031
486374o02may990p03.rd 552805o05oct003p09.rd 333 12 0.227 0.036
552803o05oct003p09.rd 552805o05oct003p09.rd 406 14 0.160 0.034
552856o05oct003p09.rd 552805o05oct003p09.rd 388 17 0.263 0.044
486370o02may990p03.rd 552856o05oct003p09.rd 314 10 0.232 0.032
552956o05oct003p12.rd 552959o05oct003p12.rd 401 13 0.220 0.032
550368o05sep000p04.rd 550373o05sep000p04.rd 399 14 0.224 0.035
550369o05sep000p04.rd 550373o05sep000p04.rd 388 15 0.236 0.039
552579o06oct002p13.rd 552574o06oct002p13.rd 286 9 0.226 0.031
552574o06oct002p13.rd 552579o06oct002p13.rd 286 9 0.226 0.031
552855o06oct003p08.rd 552798o06oct003p08.rd 477 19 0.273 0.040
552798o06oct003p08.rd 552855o06oct003p08.rd 467 18 0.273 0.039
486371o10may990p03.rd 486370o10may990p03.rd 259 9 0.220 0.035
486372o10may990p03.rd 486370o10may990p03.rd 252 8 0.222 0.032
486370o10may990p03.rd 486371o10may990p03.rd 259 9 0.220 0.035
486370o10may990p03.rd 486372o10may990p03.rd 251 8 0.222 0.032
looks good. let update the zero-points
# move to lasa28
cd /raid/gwyn/F22/F22I
tar -cvf - *.fead | gzip > fead.tar.gz
on lasa28:
cd /home/gwyn/pg/F22I
ncftp -u virmos pix1.iap.fr
cd /raid/gwyn/F22/F22I
get fead.tar.gz
bye
gzip -dc fead.tar.gz | tar -xvf -
# apply the zero-points.
applyzp.pl
ls *.fead | perl -nle '($a=$_)=~s/fead/head/; system "cp $_ $a"'
tar -cvf - *.head| gzip > head.tar.gz
#move to pix1
ncftp -u virmos pix1.iap.fr
cd /raid/gwyn/F22/F22I
put head.tar.gz
gzip -dc head.tar.gz | tar -xvf -
# get ready to swarp
telnet samplix.iap.fr
cd /samplix1/gwyn/F22
mkdir swarpF22I
cd swarpF22I
ncftp -u virmos pix1.iap.fr
cd /raid/gwyn/F22/F22I
get *.head
get *.fits
bye
mkdir masks
cd masks
ncftp -u virmos pix1.iap.fr
cd /raid/gwyn/F22/F22I/data/
get mask*.fits
bye
cd ..
ls *.fits >F22I.swarplist
mkdir stacks
mkdir resamp
mkdir vmem
get F22I.swarp
it looks like this: (note this may not be the same as was used for the first
image! mapix died, taking the original with it)
# Default configuration file for SWarp 1.33# EB 30/04/2002
#
#----------------------------------- Output -----------------------------------
IMAGEOUT_NAME F22Isw.fits # Output filename
WEIGHTOUT_NAME F22Isw.weight.fits # Output weight-map filename
HEADER_ONLY N # Only a header as an output file (Y/N)?
HEADER_SUFFIX .head # Filename extension for additional headers
#------------------------------- Input Weights --------------------------------
WEIGHT_TYPE MAP_WEIGHT # "BACKGROUND","MAP_RMS","MAP_VARIANCE"
# or "MAP_WEIGHT"
WEIGHT_SUFFIX .weight.fits # Suffix to use for weight-maps
WEIGHT_IMAGE # Weightmap filename if suffix not used
# (all or for each weight-map)
#------------------------------- Co-addition ----------------------------------
COMBINE Y # Combine resampled images (Y/N)?
COMBINE_TYPE MEDIAN # "MEDIAN","AVERAGE","MIN","MAX"
# or "WEIGHTED"
#-------------------------------- Astrometry ----------------------------------
CELESTIAL_TYPE NATIVE # "NATIVE", "EQUATORIAL", "GALACTIC",
# or "ECLIPTIC"
PROJECTION_TYPE TAN # Any WCS projection code or "NONE"
CENTER_TYPE ALL # "MANUAL", "ALL" or "MOST"
CENTER 00:00:00.0, +00:00:00.0 # Coordinates of the image center
#CENTER_TYPE MANUAL
#CENTER 22:20:00.0 +00:00:00.0
PIXELSCALE_TYPE MANUAL # "MANUAL", "FIT", "MIN", "MAX"
# or "MEDIAN"
PIXEL_SCALE 0.206 # Pixel scale
IMAGE_SIZE 0 # Image size (0 = AUTOMATIC)
#IMAGE_SIZE 1000,1000
#-------------------------------- Resampling ----------------------------------
RESAMPLE Y # Resample input images (Y/N)?
RESAMPLE_DIR /raid/gwyn/F22/F22I/swarp/resamp/ # Directory path for resampled images
RESAMPLE_SUFFIX .resamp.fits # filename extension for resampled images
RESAMPLING_TYPE LANCZOS3 # "NEAREST","BILINEAR","LANCZOS2",
# "LANCZOS3" or "LANCZOS4" (1 per axis)
OVERSAMPLING 0 # Oversampling in each dimension
# (0 = automatic)
INTERPOLATE N # Interpolate bad input pixels (Y/N)?
# (all or for each image)
FSCALASTRO_TYPE FIXED # "NONE" or "FIXED"
FSCALE_KEYWORD FSCALE # FITS keyword for the multiplicative
# factor to apply to each input image
FSCALE_DEFAULT 1.0 # Default FSCALE value if not in header
GAIN_KEYWORD GAIN # FITS keyword for eff. gain (e-/ADU)
GAIN_DEFAULT 0.0 # Default gain if no FITS keyword found
# 0 = infty (all or for each image)
#--------------------------- Background subtraction ---------------------------
SUBTRACT_BACK Y # Subtraction sky background (Y/N)?
# (all or for each image)
BACK_TYPE AUTO # "AUTO" or "MANUAL"
# (all or for each image)
BACK_DEFAULT 0.0 # Default background value in MANUAL
# (all or for each image)
BACK_SIZE 128 # Background mesh size (pixels)
# (all or for each image)
BACK_FILTERSIZE 3 # Background map filter range (meshes)
# (all or for each image)
#------------------------- Virtual memory management --------------------------
VMEM_DIR /raid/gwyn/F22/F22I/swarp/vmem/ # Directory path for swap files
VMEM_MAX 2047 # Maximum amount of virtual memory (MB)
MEM_MAX 512 # Maximum amount of usable RAM (MB)
#------------------------------ Miscellaneous ---------------------------------
DELETE_TMPFILES N # Delete temporary resampled FITS files
# (Y/N)?
WRITE_FILEINFO N # Write information about each input
# file in the output image header?
VERBOSE_TYPE FULL # "QUIET","NORMAL" or "FULL"
NTHREADS 1 # Number of simultaneous threads for
# the SMP version of SWarp
The actual command to be run is:
doswarp.pl F22I.swarplist F22Ihole.fits
which generates a swarp command
swarp 486370o00may990p03.fits 486370o01may990p03.fits 486370o02may990p03.fits 486370o03may990p03.fits 486370o04may990p03.fits 486370o05may990p03.fits 486370o07may990p03.fits 486370o08may990p03.fits 486370o09may990p03.fits 486370o10may990p03.fits 486370o11may990p03.fits 486371o00may990p03.fits 486371o01may990p03.fits 486371o02may990p03.fits 486371o03may990p03.fits 486371o04may990p03.fits 486371o05may990p03.fits 486371o07may990p03.fits 486371o08may990p03.fits 486371o09may990p03.fits 486371o10may990p03.fits 486371o11may990p03.fits 486372o00may990p03.fits 486372o01may990p03.fits 486372o02may990p03.fits 486372o03may990p03.fits 486372o04may990p03.fits 486372o05may990p03.fits 486372o07may990p03.fits 486372o08may990p03.fits 486372o09may990p03.fits 486372o10may990p03.fits 486372o11may990p03.fits 486373o00may990p03.fits 486373o01may990p03.fits 486373o02may990p03.fits 486373o03may990p03.fits 486373o04may990p03.fits 486373o05may990p03.fits 486373o07may990p03.fits 486373o08may990p03.fits 486373o09may990p03.fits 486373o10may990p03.fits 486373o11may990p03.fits 486374o00may990p03.fits 486374o01may990p03.fits 486374o02may990p03.fits 486374o03may990p03.fits 486374o04may990p03.fits 486374o05may990p03.fits 486374o07may990p03.fits 486374o08may990p03.fits 486374o09may990p03.fits 486374o10may990p03.fits 486374o11may990p03.fits 508265o00nov990p04.fits 508265o01nov990p04.fits 508265o02nov990p04.fits 508265o03nov990p04.fits 508265o04nov990p04.fits 508265o05nov990p04.fits 508265o06nov990p04.fits 508265o07nov990p04.fits 508265o08nov990p04.fits 508265o09nov990p04.fits 508265o10nov990p04.fits 508265o11nov990p04.fits 508266o00nov990p04.fits 508266o01nov990p04.fits 508266o02nov990p04.fits 508266o03nov990p04.fits 508266o04nov990p04.fits 508266o05nov990p04.fits 508266o06nov990p04.fits 508266o07nov990p04.fits 508266o08nov990p04.fits 508266o09nov990p04.fits 508266o10nov990p04.fits 508266o11nov990p04.fits 508267o00nov990p04.fits 508267o01nov990p04.fits 508267o02nov990p04.fits 508267o03nov990p04.fits 508267o04nov990p04.fits 508267o05nov990p04.fits 508267o06nov990p04.fits 508267o07nov990p04.fits 508267o08nov990p04.fits 508267o09nov990p04.fits 508267o10nov990p04.fits 508267o11nov990p04.fits 508268o00nov990p04.fits 508268o01nov990p04.fits 508268o02nov990p04.fits 508268o03nov990p04.fits 508268o04nov990p04.fits 508268o05nov990p04.fits 508268o06nov990p04.fits 508268o07nov990p04.fits 508268o08nov990p04.fits 508268o09nov990p04.fits 508268o10nov990p04.fits 508268o11nov990p04.fits 508269o00nov990p04.fits 508269o01nov990p04.fits 508269o02nov990p04.fits 508269o03nov990p04.fits 508269o04nov990p04.fits 508269o05nov990p04.fits 508269o06nov990p04.fits 508269o07nov990p04.fits 508269o08nov990p04.fits 508269o09nov990p04.fits 508269o10nov990p04.fits 508269o11nov990p04.fits 531182o00apr002p06.fits 531182o01apr002p06.fits 531182o02apr002p06.fits 531182o03apr002p06.fits 531182o04apr002p06.fits 531182o05apr002p06.fits 531182o06apr002p06.fits 531182o07apr002p06.fits 531182o08apr002p06.fits 531182o09apr002p06.fits 531182o10apr002p06.fits 531182o11apr002p06.fits 531183o00apr002p06.fits 531183o01apr002p06.fits 531183o02apr002p06.fits 531183o03apr002p06.fits 531183o04apr002p06.fits 531183o05apr002p06.fits 531183o06apr002p06.fits 531183o07apr002p06.fits 531183o08apr002p06.fits 531183o09apr002p06.fits 531183o10apr002p06.fits 531183o11apr002p06.fits 531288o00apr002p06.fits 531288o01apr002p06.fits 531288o02apr002p06.fits 531288o03apr002p06.fits 531288o04apr002p06.fits 531288o05apr002p06.fits 531288o06apr002p06.fits 531288o07apr002p06.fits 531288o08apr002p06.fits 531288o09apr002p06.fits 531288o10apr002p06.fits 531288o11apr002p06.fits 531289o00apr002p06.fits 531289o01apr002p06.fits 531289o02apr002p06.fits 531289o03apr002p06.fits 531289o04apr002p06.fits 531289o05apr002p06.fits 531289o06apr002p06.fits 531289o07apr002p06.fits 531289o08apr002p06.fits 531289o09apr002p06.fits 531289o10apr002p06.fits 531289o11apr002p06.fits 531290o00apr002p06.fits 531290o01apr002p06.fits 531290o02apr002p06.fits 531290o03apr002p06.fits 531290o04apr002p06.fits 531290o05apr002p06.fits 531290o06apr002p06.fits 531290o07apr002p06.fits 531290o08apr002p06.fits 531290o09apr002p06.fits 531290o10apr002p06.fits 531290o11apr002p06.fits 531630o00apr002p11.fits 531630o01apr002p11.fits 531630o02apr002p11.fits 531630o03apr002p11.fits 531630o04apr002p11.fits 531630o05apr002p11.fits 531630o06apr002p11.fits 531630o07apr002p11.fits 531630o08apr002p11.fits 531630o09apr002p11.fits 531630o10apr002p11.fits 531630o11apr002p11.fits 531631o00apr002p11.fits 531631o01apr002p11.fits 531631o02apr002p11.fits 531631o03apr002p11.fits 531631o04apr002p11.fits 531631o05apr002p11.fits 531631o06apr002p11.fits 531631o07apr002p11.fits 531631o08apr002p11.fits 531631o09apr002p11.fits 531631o10apr002p11.fits 531631o11apr002p11.fits 531632o00apr002p11.fits 531632o01apr002p11.fits 531632o02apr002p11.fits 531632o03apr002p11.fits 531632o04apr002p11.fits 531632o05apr002p11.fits 531632o06apr002p11.fits 531632o07apr002p11.fits 531632o08apr002p11.fits 531632o09apr002p11.fits 531632o10apr002p11.fits 531632o11apr002p11.fits 550361o00sep000p03.fits 550361o01sep000p03.fits 550361o02sep000p03.fits 550361o03sep000p03.fits 550361o04sep000p03.fits 550361o05sep000p03.fits 550361o06sep000p03.fits 550361o07sep000p03.fits 550361o08sep000p03.fits 550361o09sep000p03.fits 550361o10sep000p03.fits 550361o11sep000p03.fits 550362o00sep000p03.fits 550362o01sep000p03.fits 550362o02sep000p03.fits 550362o03sep000p03.fits 550362o04sep000p03.fits 550362o05sep000p03.fits 550362o06sep000p03.fits 550362o07sep000p03.fits 550362o08sep000p03.fits 550362o09sep000p03.fits 550362o10sep000p03.fits 550362o11sep000p03.fits 550363o00sep000p03.fits 550363o01sep000p03.fits 550363o02sep000p03.fits 550363o03sep000p03.fits 550363o04sep000p03.fits 550363o05sep000p03.fits 550363o06sep000p03.fits 550363o07sep000p03.fits 550363o08sep000p03.fits 550363o09sep000p03.fits 550363o10sep000p03.fits 550363o11sep000p03.fits 550366o00sep000p03.fits 550366o01sep000p03.fits 550366o02sep000p03.fits 550366o03sep000p03.fits 550366o04sep000p03.fits 550366o05sep000p03.fits 550366o06sep000p03.fits 550366o07sep000p03.fits 550366o08sep000p03.fits 550366o09sep000p03.fits 550366o10sep000p03.fits 550366o11sep000p03.fits 550367o00sep000p04.fits 550367o01sep000p04.fits 550367o02sep000p04.fits 550367o03sep000p04.fits 550367o04sep000p04.fits 550367o05sep000p04.fits 550367o06sep000p04.fits 550367o07sep000p04.fits 550367o08sep000p04.fits 550367o09sep000p04.fits 550367o10sep000p04.fits 550367o11sep000p04.fits 550368o00sep000p04.fits 550368o01sep000p04.fits 550368o02sep000p04.fits 550368o03sep000p04.fits 550368o04sep000p04.fits 550368o05sep000p04.fits 550368o06sep000p04.fits 550368o07sep000p04.fits 550368o08sep000p04.fits 550368o09sep000p04.fits 550368o10sep000p04.fits 550368o11sep000p04.fits 550369o00sep000p04.fits 550369o01sep000p04.fits 550369o02sep000p04.fits 550369o03sep000p04.fits 550369o04sep000p04.fits 550369o05sep000p04.fits 550369o06sep000p04.fits 550369o07sep000p04.fits 550369o08sep000p04.fits 550369o09sep000p04.fits 550369o10sep000p04.fits 550369o11sep000p04.fits 550373o00sep000p04.fits 550373o01sep000p04.fits 550373o02sep000p04.fits 550373o03sep000p04.fits 550373o04sep000p04.fits 550373o05sep000p04.fits 550373o06sep000p04.fits 550373o07sep000p04.fits 550373o08sep000p04.fits 550373o09sep000p04.fits 550373o10sep000p04.fits 550373o11sep000p04.fits 550374o00sep000p04.fits 550374o01sep000p04.fits 550374o02sep000p04.fits 550374o03sep000p04.fits 550374o04sep000p04.fits 550374o05sep000p04.fits 550374o06sep000p04.fits 550374o07sep000p04.fits 550374o08sep000p04.fits 550374o09sep000p04.fits 550374o10sep000p04.fits 550374o11sep000p04.fits 550375o00sep000p05.fits 550375o01sep000p05.fits 550375o02sep000p05.fits 550375o03sep000p05.fits 550375o04sep000p05.fits 550375o05sep000p05.fits 550375o06sep000p05.fits 550375o07sep000p05.fits 550375o08sep000p05.fits 550375o09sep000p05.fits 550375o10sep000p05.fits 550375o11sep000p05.fits 550376o00sep000p05.fits 550376o01sep000p05.fits 550376o02sep000p05.fits 550376o03sep000p05.fits 550376o04sep000p05.fits 550376o05sep000p05.fits 550376o06sep000p05.fits 550376o07sep000p05.fits 550376o08sep000p05.fits 550376o09sep000p05.fits 550376o10sep000p05.fits 550376o11sep000p05.fits 550377o00sep000p05.fits 550377o01sep000p05.fits 550377o02sep000p05.fits 550377o03sep000p05.fits 550377o04sep000p05.fits 550377o05sep000p05.fits 550377o06sep000p05.fits 550377o07sep000p05.fits 550377o08sep000p05.fits 550377o09sep000p05.fits 550377o10sep000p05.fits 550377o11sep000p05.fits 550378o00sep000p05.fits 550378o01sep000p05.fits 550378o02sep000p05.fits 550378o03sep000p05.fits 550378o04sep000p05.fits 550378o05sep000p05.fits 550378o06sep000p05.fits 550378o07sep000p05.fits 550378o08sep000p05.fits 550378o09sep000p05.fits 550378o10sep000p05.fits 550378o11sep000p05.fits 550379o00sep000p05.fits 550379o01sep000p05.fits 550379o02sep000p05.fits 550379o03sep000p05.fits 550379o04sep000p05.fits 550379o05sep000p05.fits 550379o06sep000p05.fits 550379o07sep000p05.fits 550379o08sep000p05.fits 550379o09sep000p05.fits 550379o10sep000p05.fits 550379o11sep000p05.fits 551371o00oct001p06.fits 551371o01oct001p06.fits 551371o02oct001p06.fits 551371o03oct001p06.fits 551371o04oct001p06.fits 551371o05oct001p06.fits 551371o06oct001p06.fits 551371o07oct001p06.fits 551371o08oct001p06.fits 551371o09oct001p06.fits 551371o10oct001p06.fits 551371o11oct001p06.fits 551372o00oct001p06.fits 551372o01oct001p06.fits 551372o02oct001p06.fits 551372o03oct001p06.fits 551372o04oct001p06.fits 551372o05oct001p06.fits 551372o06oct001p06.fits 551372o07oct001p06.fits 551372o08oct001p06.fits 551372o09oct001p06.fits 551372o10oct001p06.fits 551372o11oct001p06.fits 551373o00oct001p06.fits 551373o01oct001p06.fits 551373o02oct001p06.fits 551373o03oct001p06.fits 551373o04oct001p06.fits 551373o05oct001p06.fits 551373o06oct001p06.fits 551373o07oct001p06.fits 551373o08oct001p06.fits 551373o09oct001p06.fits 551373o10oct001p06.fits 551373o11oct001p06.fits 551374o00oct001p06.fits 551374o01oct001p06.fits 551374o02oct001p06.fits 551374o03oct001p06.fits 551374o04oct001p06.fits 551374o05oct001p06.fits 551374o06oct001p06.fits 551374o07oct001p06.fits 551374o08oct001p06.fits 551374o09oct001p06.fits 551374o10oct001p06.fits 551374o11oct001p06.fits 551375o00oct001p06.fits 551375o01oct001p06.fits 551375o02oct001p06.fits 551375o03oct001p06.fits 551375o04oct001p06.fits 551375o05oct001p06.fits 551375o06oct001p06.fits 551375o07oct001p06.fits 551375o08oct001p06.fits 551375o09oct001p06.fits 551375o10oct001p06.fits 551375o11oct001p06.fits 551376o00oct001p06.fits 551376o01oct001p06.fits 551376o02oct001p06.fits 551376o03oct001p06.fits 551376o04oct001p06.fits 551376o05oct001p06.fits 551376o06oct001p06.fits 551376o07oct001p06.fits 551376o08oct001p06.fits 551376o09oct001p06.fits 551376o10oct001p06.fits 551376o11oct001p06.fits 551377o00oct001p06.fits 551377o01oct001p06.fits 551377o02oct001p06.fits 551377o03oct001p06.fits 551377o04oct001p06.fits 551377o05oct001p06.fits 551377o06oct001p06.fits 551377o07oct001p06.fits 551377o08oct001p06.fits 551377o09oct001p06.fits 551377o10oct001p06.fits 551377o11oct001p06.fits 551379o00oct001p07.fits 551379o01oct001p07.fits 551379o02oct001p07.fits 551379o03oct001p07.fits 551379o04oct001p07.fits 551379o05oct001p07.fits 551379o06oct001p07.fits 551379o07oct001p07.fits 551379o08oct001p07.fits 551379o09oct001p07.fits 551379o10oct001p07.fits 551379o11oct001p07.fits 551380o00oct001p07.fits 551380o01oct001p07.fits 551380o02oct001p07.fits 551380o03oct001p07.fits 551380o04oct001p07.fits 551380o05oct001p07.fits 551380o06oct001p07.fits 551380o07oct001p07.fits 551380o08oct001p07.fits 551380o09oct001p07.fits 551380o10oct001p07.fits 551380o11oct001p07.fits 551381o00oct001p07.fits 551381o01oct001p07.fits 551381o02oct001p07.fits 551381o03oct001p07.fits 551381o04oct001p07.fits 551381o05oct001p07.fits 551381o06oct001p07.fits 551381o07oct001p07.fits 551381o08oct001p07.fits 551381o09oct001p07.fits 551381o10oct001p07.fits 551381o11oct001p07.fits 551382o00oct001p07.fits 551382o01oct001p07.fits 551382o02oct001p07.fits 551382o03oct001p07.fits 551382o04oct001p07.fits 551382o05oct001p07.fits 551382o06oct001p07.fits 551382o07oct001p07.fits 551382o08oct001p07.fits 551382o09oct001p07.fits 551382o10oct001p07.fits 551382o11oct001p07.fits 551383o00oct001p07.fits 551383o01oct001p07.fits 551383o02oct001p07.fits 551383o03oct001p07.fits 551383o04oct001p07.fits 551383o05oct001p07.fits 551383o06oct001p07.fits 551383o07oct001p07.fits 551383o08oct001p07.fits 551383o09oct001p07.fits 551383o10oct001p07.fits 551383o11oct001p07.fits 551384o00oct001p07.fits 551384o01oct001p07.fits 551384o02oct001p07.fits 551384o03oct001p07.fits 551384o04oct001p07.fits 551384o05oct001p07.fits 551384o06oct001p07.fits 551384o07oct001p07.fits 551384o08oct001p07.fits 551384o09oct001p07.fits 551384o10oct001p07.fits 551384o11oct001p07.fits 552573o01oct002p13.fits 552573o02oct002p13.fits 552573o03oct002p13.fits 552573o04oct002p13.fits 552573o05oct002p13.fits 552573o06oct002p13.fits 552573o07oct002p13.fits 552573o08oct002p13.fits 552573o09oct002p13.fits 552573o10oct002p13.fits 552573o11oct002p13.fits 552574o01oct002p13.fits 552574o02oct002p13.fits 552574o03oct002p13.fits 552574o04oct002p13.fits 552574o05oct002p13.fits 552574o06oct002p13.fits 552574o07oct002p13.fits 552574o08oct002p13.fits 552574o09oct002p13.fits 552574o10oct002p13.fits 552574o11oct002p13.fits 552576o01oct002p13.fits 552576o02oct002p13.fits 552576o03oct002p13.fits 552576o04oct002p13.fits 552576o05oct002p13.fits 552576o06oct002p13.fits 552576o07oct002p13.fits 552576o08oct002p13.fits 552576o09oct002p13.fits 552576o10oct002p13.fits 552576o11oct002p13.fits 552577o01oct002p13.fits 552577o02oct002p13.fits 552577o03oct002p13.fits 552577o04oct002p13.fits 552577o05oct002p13.fits 552577o06oct002p13.fits 552577o07oct002p13.fits 552577o08oct002p13.fits 552577o09oct002p13.fits 552577o10oct002p13.fits 552577o11oct002p13.fits 552578o01oct002p13.fits 552578o02oct002p13.fits 552578o03oct002p13.fits 552578o04oct002p13.fits 552578o05oct002p13.fits 552578o06oct002p13.fits 552578o07oct002p13.fits 552578o08oct002p13.fits 552578o09oct002p13.fits 552578o10oct002p13.fits 552578o11oct002p13.fits 552579o01oct002p13.fits 552579o02oct002p13.fits 552579o03oct002p13.fits 552579o04oct002p13.fits 552579o05oct002p13.fits 552579o06oct002p13.fits 552579o07oct002p13.fits 552579o08oct002p13.fits 552579o09oct002p13.fits 552579o10oct002p13.fits 552579o11oct002p13.fits 552797o00oct003p08.fits 552797o01oct003p08.fits 552797o02oct003p08.fits 552797o03oct003p08.fits 552797o04oct003p08.fits 552797o05oct003p08.fits 552797o06oct003p08.fits 552797o07oct003p08.fits 552797o08oct003p08.fits 552797o09oct003p08.fits 552797o10oct003p08.fits 552797o11oct003p08.fits 552798o00oct003p08.fits 552798o01oct003p08.fits 552798o02oct003p08.fits 552798o03oct003p08.fits 552798o04oct003p08.fits 552798o05oct003p08.fits 552798o06oct003p08.fits 552798o07oct003p08.fits 552798o08oct003p08.fits 552798o09oct003p08.fits 552798o10oct003p08.fits 552798o11oct003p08.fits 552799o00oct003p08.fits 552799o01oct003p08.fits 552799o02oct003p08.fits 552799o03oct003p08.fits 552799o04oct003p08.fits 552799o05oct003p08.fits 552799o06oct003p08.fits 552799o07oct003p08.fits 552799o08oct003p08.fits 552799o09oct003p08.fits 552799o10oct003p08.fits 552799o11oct003p08.fits 552800o00oct003p08.fits 552800o01oct003p08.fits 552800o02oct003p08.fits 552800o03oct003p08.fits 552800o04oct003p08.fits 552800o05oct003p08.fits 552800o06oct003p08.fits 552800o07oct003p08.fits 552800o08oct003p08.fits 552800o09oct003p08.fits 552800o10oct003p08.fits 552800o11oct003p08.fits 552801o00oct003p08.fits 552801o01oct003p08.fits 552801o02oct003p08.fits 552801o03oct003p08.fits 552801o04oct003p08.fits 552801o05oct003p08.fits 552801o06oct003p08.fits 552801o07oct003p08.fits 552801o08oct003p08.fits 552801o09oct003p08.fits 552801o10oct003p08.fits 552801o11oct003p08.fits 552802o00oct003p09.fits 552802o01oct003p09.fits 552802o02oct003p09.fits 552802o03oct003p09.fits 552802o04oct003p09.fits 552802o05oct003p09.fits 552802o06oct003p09.fits 552802o07oct003p09.fits 552802o08oct003p09.fits 552802o09oct003p09.fits 552802o10oct003p09.fits 552802o11oct003p09.fits 552803o00oct003p09.fits 552803o01oct003p09.fits 552803o02oct003p09.fits 552803o03oct003p09.fits 552803o04oct003p09.fits 552803o05oct003p09.fits 552803o06oct003p09.fits 552803o07oct003p09.fits 552803o08oct003p09.fits 552803o09oct003p09.fits 552803o10oct003p09.fits 552803o11oct003p09.fits 552804o00oct003p09.fits 552804o01oct003p09.fits 552804o02oct003p09.fits 552804o03oct003p09.fits 552804o04oct003p09.fits 552804o05oct003p09.fits 552804o06oct003p09.fits 552804o07oct003p09.fits 552804o08oct003p09.fits 552804o09oct003p09.fits 552804o10oct003p09.fits 552804o11oct003p09.fits 552805o00oct003p09.fits 552805o01oct003p09.fits 552805o02oct003p09.fits 552805o03oct003p09.fits 552805o04oct003p09.fits 552805o05oct003p09.fits 552805o06oct003p09.fits 552805o07oct003p09.fits 552805o08oct003p09.fits 552805o09oct003p09.fits 552805o10oct003p09.fits 552805o11oct003p09.fits 552806o00oct003p09.fits 552806o01oct003p09.fits 552806o02oct003p09.fits 552806o03oct003p09.fits 552806o04oct003p09.fits 552806o05oct003p09.fits 552806o06oct003p09.fits 552806o07oct003p09.fits 552806o08oct003p09.fits 552806o09oct003p09.fits 552806o10oct003p09.fits 552806o11oct003p09.fits 552850o00oct003p10.fits 552850o01oct003p10.fits 552850o02oct003p10.fits 552850o03oct003p10.fits 552850o04oct003p10.fits 552850o05oct003p10.fits 552850o06oct003p10.fits 552850o07oct003p10.fits 552850o08oct003p10.fits 552850o09oct003p10.fits 552850o10oct003p10.fits 552850o11oct003p10.fits 552851o00oct003p10.fits 552851o01oct003p10.fits 552851o02oct003p10.fits 552851o03oct003p10.fits 552851o04oct003p10.fits 552851o05oct003p10.fits 552851o06oct003p10.fits 552851o07oct003p10.fits 552851o08oct003p10.fits 552851o09oct003p10.fits 552851o10oct003p10.fits 552851o11oct003p10.fits 552852o00oct003p10.fits 552852o01oct003p10.fits 552852o02oct003p10.fits 552852o03oct003p10.fits 552852o04oct003p10.fits 552852o05oct003p10.fits 552852o06oct003p10.fits 552852o07oct003p10.fits 552852o08oct003p10.fits 552852o09oct003p10.fits 552852o10oct003p10.fits 552852o11oct003p10.fits 552853o00oct003p10.fits 552853o01oct003p10.fits 552853o02oct003p10.fits 552853o03oct003p10.fits 552853o04oct003p10.fits 552853o05oct003p10.fits 552853o06oct003p10.fits 552853o07oct003p10.fits 552853o08oct003p10.fits 552853o09oct003p10.fits 552853o10oct003p10.fits 552853o11oct003p10.fits 552854o00oct003p10.fits 552854o01oct003p10.fits 552854o02oct003p10.fits 552854o03oct003p10.fits 552854o04oct003p10.fits 552854o05oct003p10.fits 552854o06oct003p10.fits 552854o07oct003p10.fits 552854o08oct003p10.fits 552854o09oct003p10.fits 552854o10oct003p10.fits 552854o11oct003p10.fits 552855o00oct003p08.fits 552855o01oct003p08.fits 552855o02oct003p08.fits 552855o03oct003p08.fits 552855o04oct003p08.fits 552855o05oct003p08.fits 552855o06oct003p08.fits 552855o07oct003p08.fits 552855o08oct003p08.fits 552855o09oct003p08.fits 552855o10oct003p08.fits 552855o11oct003p08.fits 552856o00oct003p09.fits 552856o01oct003p09.fits 552856o02oct003p09.fits 552856o03oct003p09.fits 552856o04oct003p09.fits 552856o05oct003p09.fits 552856o06oct003p09.fits 552856o07oct003p09.fits 552856o08oct003p09.fits 552856o09oct003p09.fits 552856o10oct003p09.fits 552856o11oct003p09.fits 552858o00oct003p11.fits 552858o01oct003p11.fits 552858o02oct003p11.fits 552858o03oct003p11.fits 552858o04oct003p11.fits 552858o05oct003p11.fits 552858o06oct003p11.fits 552858o07oct003p11.fits 552858o08oct003p11.fits 552858o09oct003p11.fits 552858o10oct003p11.fits 552858o11oct003p11.fits 552859o00oct003p11.fits 552859o01oct003p11.fits 552859o02oct003p11.fits 552859o03oct003p11.fits 552859o04oct003p11.fits 552859o05oct003p11.fits 552859o06oct003p11.fits 552859o07oct003p11.fits 552859o08oct003p11.fits 552859o09oct003p11.fits 552859o10oct003p11.fits 552859o11oct003p11.fits 552860o00oct003p11.fits 552860o01oct003p11.fits 552860o02oct003p11.fits 552860o03oct003p11.fits 552860o04oct003p11.fits 552860o05oct003p11.fits 552860o06oct003p11.fits 552860o07oct003p11.fits 552860o08oct003p11.fits 552860o09oct003p11.fits 552860o10oct003p11.fits 552860o11oct003p11.fits 552953o00oct003p11.fits 552953o01oct003p11.fits 552953o02oct003p11.fits 552953o03oct003p11.fits 552953o04oct003p11.fits 552953o05oct003p11.fits 552953o06oct003p11.fits 552953o07oct003p11.fits 552953o08oct003p11.fits 552953o09oct003p11.fits 552953o10oct003p11.fits 552953o11oct003p11.fits 552954o00oct003p11.fits 552954o01oct003p11.fits 552954o02oct003p11.fits 552954o03oct003p11.fits 552954o04oct003p11.fits 552954o05oct003p11.fits 552954o06oct003p11.fits 552954o07oct003p11.fits 552954o08oct003p11.fits 552954o09oct003p11.fits 552954o10oct003p11.fits 552954o11oct003p11.fits 552955o00oct003p12.fits 552955o01oct003p12.fits 552955o02oct003p12.fits 552955o03oct003p12.fits 552955o04oct003p12.fits 552955o05oct003p12.fits 552955o06oct003p12.fits 552955o07oct003p12.fits 552955o08oct003p12.fits 552955o09oct003p12.fits 552955o10oct003p12.fits 552955o11oct003p12.fits 552956o00oct003p12.fits 552956o01oct003p12.fits 552956o02oct003p12.fits 552956o03oct003p12.fits 552956o04oct003p12.fits 552956o05oct003p12.fits 552956o06oct003p12.fits 552956o07oct003p12.fits 552956o08oct003p12.fits 552956o09oct003p12.fits 552956o10oct003p12.fits 552956o11oct003p12.fits 552957o00oct003p12.fits 552957o01oct003p12.fits 552957o02oct003p12.fits 552957o03oct003p12.fits 552957o04oct003p12.fits 552957o05oct003p12.fits 552957o06oct003p12.fits 552957o07oct003p12.fits 552957o08oct003p12.fits 552957o09oct003p12.fits 552957o10oct003p12.fits 552957o11oct003p12.fits 552958o00oct003p12.fits 552958o01oct003p12.fits 552958o02oct003p12.fits 552958o03oct003p12.fits 552958o04oct003p12.fits 552958o05oct003p12.fits 552958o06oct003p12.fits 552958o07oct003p12.fits 552958o08oct003p12.fits 552958o09oct003p12.fits 552958o10oct003p12.fits 552958o11oct003p12.fits 552959o00oct003p12.fits 552959o01oct003p12.fits 552959o02oct003p12.fits 552959o03oct003p12.fits 552959o04oct003p12.fits 552959o05oct003p12.fits 552959o06oct003p12.fits 552959o07oct003p12.fits 552959o08oct003p12.fits 552959o09oct003p12.fits 552959o10oct003p12.fits 552959o11oct003p12.fits -c /samplix1/gwyn/F22/swarpF22I/F22I.swarp -VMEM_DIR /samplix1/gwyn/F22/swarpF22I/vmem/ -RESAMPLE_DIR /samplix1/gwyn/F22/swarpF22I/resamp/ -VERBOSE_TYPE NORMAL -IMAGEOUT_NAME /samplix1/gwyn/F22/swarpF22I/stacks/F22Ihole.fits -WEIGHTOUT_NAME /samplix1/gwyn/F22/swarpF22I/stacks/F22Ihole.weight.fits -WEIGHT_TYPE MAP_WEIGHT
Coadding 925 CCDs... here goes!
shall I proceed (y/n)
** 7 October 2002
The Swarp failed due to lack of disk space.
Start again on mapix
cd /mapxi2/gwyn
mkdir F22
ncftp -u gwyn samplix.iap.fr
cd /samplix1/gwyn/F22
get -R swarpF22I
bye
cd /mapix4/gwyn/
mkdir F22
cd F22
mkdir swarpF22I
cd swarpF22I
mkdir stacks
mkdir vmem
mkdir resamp
cd
ln -s /mapix4/gwyn/F22/swarpF22I/vmem/ .
ln -s /mapix4/gwyn/F22/swarpF22I/stacks
ln -s /mapix4/gwyn/F22/swarpF22I/resamp
doswarp.pl F22I.swarplist F22Ihole.fits >swarp.log
"no ident" unweighted no ext. header 37573x34187 32 bits (floats)
Center: 22:17:50.69 +00:23:04.7 129'x117' Scale: 0.206 ''/pixel
Gain: 0 e-/ADU Flux scaling (astrom/photom): 1 X / 1 X
> Framing and scaling images...
** 9 October 2002
-------------- File 552959o11oct003p12.fits:
"no ident" WEIGHTED EXT. HEADER 2048x4096 16 bits (integers)
Center: 22:21:40.66 -00:25:36.8 6.92'x13.8' Scale: 0.2027 ''/pixel
Gain: 1.45 e-/ADU Flux scaling (astrom/photom): 1.032 X / 0.03112 X
Background: 23245.43 RMS: 168.4411 Weight scale: 168.4411
> Resampling line: 4076 / 4078
> Not enough memory for multibuf (nlines*multiwidth elements) !
FUCK!
Henry claims the virtual memory limits need to be set larger...
Gilles reboots, and we submit the job again...
** 10 October 2002
Swarp fails again:
"no ident" WEIGHTED EXT. HEADER 2048x4096 16 bits (integers)
Center: 22:21:40.66 -00:25:36.8 6.92'x13.8' Scale: 0.2027 ''/pixel
Gain: 1.5 e-/ADU Flux scaling (astrom/photom): 1 X / 0.03112 X
Background: 0 RMS: 0 Weight scale: 1
> Co-adding line: 944 / 34187
/mapix2/gwyn/F22/swarpF22I #13 >
/mapix2/gwyn/F22/swarpF22I #13 >
> Co-adding line: 4096 / 34187
> *Error* while reading 486370o00may990p03.fits
/u3/hjmcc/bin/toFITS2d -x 19000 -X 1000 -y 19000 -Y 1000 F22Ihole1.fits >! fred.fits
** 14 October 2002
The last swarp didn't work
once again....
doswarp.pl F22I.swarplist F22Ihole.fits > & swarp.err
** 15 October 2002
resampled images and stacks should be on different disks.
** 17 October 2002
/u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 1 -Y 17500 F22Ihole.fits >! F22Iholesw.fits
/u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 1 -Y 17500 F22Ihole.fits >! F22Iholese.fits
/u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 17500 -Y 16686 F22Ihole.fits >! F22Iholenw.fits
/u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 17500 -Y 16686 F22Ihole.fits >! F22Iholene.fits
/u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 1 -Y 17500 F22Ihole.weight.fits >! F22Iholesw.weight.fits
/u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 1 -Y 17500 F22Ihole.weight.fits >! F22Iholese.weight.fits
/u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 17500 -Y 16686 F22Ihole.weight.fits >! F22Iholenw.weight.fits
/u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 17500 -Y 16686 F22Ihole.weight.fits >! F22Iholene.weight.fits
sex -c bigimage.sex F22Iholese.fits -CATALOG_NAME F22Iholese.cat -WEIGHT_IMAGE F22Iholese.weight.fits
sex -c bigimage.sex F22Iholene.fits -CATALOG_NAME F22Iholene.cat -WEIGHT_IMAGE F22Iholene.weight.fits
sex -c bigimage.sex F22Iholenw.fits -CATALOG_NAME F22Iholenw.cat -WEIGHT_IMAGE F22Iholenw.weight.fits
sex -c bigimage.sex F22Iholesw.fits -CATALOG_NAME F22Iholesw.cat -WEIGHT_IMAGE F22Iholesw.weight.fits
** 18 October 2002
on lasa28:
working in
/home/gwyn/red/F22/F22I/hole
get the catalogs...
exsex F22Iholene.cat alpha_j2000 delta_j2000 mag_auto magerr_auto > F22Iholene.rdm
exsex F22Iholenw.cat alpha_j2000 delta_j2000 mag_auto magerr_auto > F22Iholenw.rdm
exsex F22Iholese.cat alpha_j2000 delta_j2000 mag_auto magerr_auto > F22Iholese.rdm
exsex F22Iholesw.cat alpha_j2000 delta_j2000 mag_auto magerr_auto > F22Iholesw.rdm
ln -s /home/gwyn/red/F22/F22I/swarpcheck/F22I??.cat
exsex F22Ine.cat alpha_j2000 delta_j2000 mag_auto magerr_auto > F22Ine.rdm
exsex F22Inw.cat alpha_j2000 delta_j2000 mag_auto magerr_auto > F22Inw.rdm
exsex F22Ise.cat alpha_j2000 delta_j2000 mag_auto magerr_auto > F22Ise.rdm
exsex F22Isw.cat alpha_j2000 delta_j2000 mag_auto magerr_auto > F22Isw.rdm
142430 9254635 105928561 F22Ine.cat
181433 11789830 134946785 F22Inw.cat
123094 7997795 91542551 F22Ise.cat
157748 10250305 117325141 F22Isw.cat
154528 10041005 114929450 F22Iholene.cat
181551 11797500 135034562 F22Iholenw.cat
132558 8612955 98583752 F22Iholese.cat
158213 10280530 117671077 F22Iholesw.cat
matchex F22Iholene.rdm F22Ine.rdm F22Ine.mat 1 4
matchex F22Iholenw.rdm F22Inw.rdm F22Inw.mat 1 4
matchex F22Iholese.rdm F22Ise.rdm F22Ise.mat 1 4
matchex F22Iholesw.rdm F22Isw.rdm F22Isw.mat 1 4
about .07" and .01 magntitudes... marginal
check out F22Idswarp.gif
generated with /home/gwyn/red/F22/F22I/hole/sm.2
exsex F22Iholene.cat alpha_j2000 delta_j2000 mag_auto magerr_auto | awk '$3<22.5' > F22Iholene.rdm
exsex F22Iholenw.cat alpha_j2000 delta_j2000 mag_auto magerr_auto | awk '$3<22.5' > F22Iholenw.rdm
exsex F22Iholese.cat alpha_j2000 delta_j2000 mag_auto magerr_auto | awk '$3<22.5' > F22Iholese.rdm
exsex F22Iholesw.cat alpha_j2000 delta_j2000 mag_auto magerr_auto | awk '$3<22.5' > F22Iholesw.rdm
exsex F22Ine.cat alpha_j2000 delta_j2000 mag_auto magerr_auto | awk '$3<22.5' > F22Ine.rdm
exsex F22Inw.cat alpha_j2000 delta_j2000 mag_auto magerr_auto | awk '$3<22.5' > F22Inw.rdm
exsex F22Ise.cat alpha_j2000 delta_j2000 mag_auto magerr_auto | awk '$3<22.5' > F22Ise.rdm
exsex F22Isw.cat alpha_j2000 delta_j2000 mag_auto magerr_auto | awk '$3<22.5' > F22Isw.rdm
30708 122832 1197612 F22Iholene.rdm
33344 133376 1300416 F22Iholenw.rdm
28899 115596 1127061 F22Iholese.rdm
30364 121456 1184196 F22Iholesw.rdm
28009 112036 1092351 F22Ine.rdm
33044 132176 1288716 F22Inw.rdm
26250 105000 1023750 F22Ise.rdm
29657 118628 1156623 F22Isw.rdm
matchex F22Iholene.rdm F22Ine.rdm F22Ine.mat 1 4
matchex F22Iholenw.rdm F22Inw.rdm F22Inw.mat 1 4
matchex F22Iholese.rdm F22Ise.rdm F22Ise.mat 1 4
matchex F22Iholesw.rdm F22Isw.rdm F22Isw.mat 1 4
ave_sep RMS 68%
F22Iholene.rdm F22Ine.rdm 23800 0.044 0.075 0.047
F22Iholenw.rdm F22Inw.rdm 28202 0.066 0.100 0.110
F22Iholese.rdm F22Ise.rdm 22681 0.046 0.081 0.096
F22Iholesw.rdm F22Isw.rdm 25630 0.070 0.109 0.110
doswarp.pl F22I.swarplist F22Iholem.fits > & swarpmatch.err
using a a configuration file
F22I.matchswarp
which looks like this:
# Default configuration file for SWarp 1.33# EB 30/04/2002
#
#----------------------------------- Output -----------------------------------
IMAGEOUT_NAME F22Isw.fits # Output filename
WEIGHTOUT_NAME F22Isw.weight.fits # Output weight-map filename
HEADER_ONLY N # Only a header as an output file (Y/N)?
HEADER_SUFFIX .head # Filename extension for additional headers
#------------------------------- Input Weights --------------------------------
WEIGHT_TYPE MAP_WEIGHT # "BACKGROUND","MAP_RMS","MAP_VARIANCE"
# or "MAP_WEIGHT"
WEIGHT_SUFFIX .weight.fits # Suffix to use for weight-maps
WEIGHT_IMAGE # Weightmap filename if suffix not used
# (all or for each weight-map)
#------------------------------- Co-addition ----------------------------------
COMBINE Y # Combine resampled images (Y/N)?
COMBINE_TYPE MEDIAN # "MEDIAN","AVERAGE","MIN","MAX"
# or "WEIGHTED"
#-------------------------------- Astrometry ----------------------------------
CELESTIAL_TYPE NATIVE # "NATIVE", "EQUATORIAL", "GALACTIC",
# or "ECLIPTIC"
PROJECTION_TYPE TAN # Any WCS projection code or "NONE"
CENTER_TYPE MANUAL # "MANUAL", "ALL" or "MOST"
CENTER 22:17:50.6894,0:23:04.67271 # Coordinates of the image center
#CENTER_TYPE MANUAL
#CENTER 22:20:00.0 +00:00:00.0
PIXELSCALE_TYPE MANUAL # "MANUAL", "FIT", "MIN", "MAX"
# or "MEDIAN"
PIXEL_SCALE 0.206 # Pixel scale
IMAGE_SIZE 37573,34187 # Image size (0 = AUTOMATIC)
#IMAGE_SIZE 1000,1000
#-------------------------------- Resampling ----------------------------------
RESAMPLE Y # Resample input images (Y/N)?
RESAMPLE_DIR /raid/gwyn/F22/F22I/swarp/resamp/ # Directory path for resampled images
RESAMPLE_SUFFIX .resamp.fits # filename extension for resampled images
RESAMPLING_TYPE LANCZOS3 # "NEAREST","BILINEAR","LANCZOS2",
# "LANCZOS3" or "LANCZOS4" (1 per axis)
OVERSAMPLING 0 # Oversampling in each dimension
# (0 = automatic)
INTERPOLATE N # Interpolate bad input pixels (Y/N)?
# (all or for each image)
FSCALASTRO_TYPE FIXED # "NONE" or "FIXED"
FSCALE_KEYWORD FSCALE # FITS keyword for the multiplicative
# factor to apply to each input image
FSCALE_DEFAULT 1.0 # Default FSCALE value if not in header
GAIN_KEYWORD GAIN # FITS keyword for eff. gain (e-/ADU)
GAIN_DEFAULT 0.0 # Default gain if no FITS keyword found
# 0 = infty (all or for each image)
#--------------------------- Background subtraction ---------------------------
SUBTRACT_BACK Y # Subtraction sky background (Y/N)?
# (all or for each image)
BACK_TYPE AUTO # "AUTO" or "MANUAL"
# (all or for each image)
BACK_DEFAULT 0.0 # Default background value in MANUAL
# (all or for each image)
BACK_SIZE 128 # Background mesh size (pixels)
# (all or for each image)
BACK_FILTERSIZE 3 # Background map filter range (meshes)
# (all or for each image)
#------------------------- Virtual memory management --------------------------
VMEM_DIR /raid/gwyn/F22/F22I/swarp/vmem/ # Directory path for swap files
VMEM_MAX 2047 # Maximum amount of virtual memory (MB)
MEM_MAX 512 # Maximum amount of usable RAM (MB)
#------------------------------ Miscellaneous ---------------------------------
DELETE_TMPFILES N # Delete temporary resampled FITS files
# (Y/N)?
WRITE_FILEINFO N # Write information about each input
# file in the output image header?
VERBOSE_TYPE FULL # "QUIET","NORMAL" or "FULL"
NTHREADS 1 # Number of simultaneous threads for
# the SMP version of SWarp
** 21 October 2002
The result of the last swarp was:
-rw-r--r-- 1 gwyn system 5138040960 Oct 21 08:10 F22Iholem.weight.fits
-rw-r--r-- 1 gwyn system 5138040960 Oct 21 08:10 F22Iholem.fits
gethead CRVAL1 CRVAL2 F22Iholem.fits
3.344612058e+02 3.846313083e-01
CENTER 22:17:50.6894,0:23:04.67271 # Coordinates of the image center
Try yet again with:
CENTER 3.344612061e+02,3.846313087e-01 # Coordinates of the image center
doswarp.pl F22I.swarplist F22Iholem2.fits >& swarpmatch2.err
meanwhile...
/mapix2/gwyn/F22/swarpF22I/stacks
mkdir try2
mv F22Iholem.weight.fits try2
mv F22Iholem.fits try2
cd try2
/u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 1 -Y 17500 F22Iholem.fits >! F22Iholemse.fits
/u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 1 -Y 17500 F22Iholem.fits >! F22Iholemsw.fits
/u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 17500 -Y 16686 F22Iholem.fits >! F22Iholemnw.fits
/u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 17500 -Y 16686 F22Iholem.fits >! F22Iholemne.fits
/u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 1 -Y 17500 F22Iholem.weight.fits >! F22Iholemse.weight.fits
/u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 1 -Y 17500 F22Iholem.weight.fits >! F22Iholemsw.weight.fits
/u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 17500 -Y 16686 F22Iholem.weight.fits >! F22Iholemnw.weight.fits
/u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 17500 -Y 16686 F22Iholem.weight.fits >! F22Iholemne.weight.fits
sex -c bigimage.sex F22Iholemse.fits -CATALOG_NAME F22Iholemse.cat -WEIGHT_IMAGE F22Iholemse.weight.fits
sex -c bigimage.sex F22Iholemne.fits -CATALOG_NAME F22Iholemne.cat -WEIGHT_IMAGE F22Iholemne.weight.fits
sex -c bigimage.sex F22Iholemnw.fits -CATALOG_NAME F22Iholemnw.cat -WEIGHT_IMAGE F22Iholemnw.weight.fits
sex -c bigimage.sex F22Iholemsw.fits -CATALOG_NAME F22Iholemsw.cat -WEIGHT_IMAGE F22Iholemsw.weight.fits
exsex F22Iholemne.cat alpha_j2000 delta_j2000 mag_auto magerr_auto | awk '$3<22.5' > F22Iholemne.rdm
exsex F22Iholemnw.cat alpha_j2000 delta_j2000 mag_auto magerr_auto | awk '$3<22.5' > F22Iholemnw.rdm
exsex F22Iholemse.cat alpha_j2000 delta_j2000 mag_auto magerr_auto | awk '$3<22.5' > F22Iholemse.rdm
exsex F22Iholemsw.cat alpha_j2000 delta_j2000 mag_auto magerr_auto | awk '$3<22.5' > F22Iholemsw.rdm
matchex F22Iholemne.rdm F22Ine.rdm F22Imne.mat 1 4
matchex F22Iholemnw.rdm F22Inw.rdm F22Imnw.mat 1 4
matchex F22Iholemse.rdm F22Ise.rdm F22Imse.mat 1 4
matchex F22Iholemsw.rdm F22Isw.rdm F22Imsw.mat 1 4
** 22 October 2002
/mapix2/gwyn/F22/swarpF22I/stacks #3 > gethead CRVAL1 CRVAL2 F22Iholem2.fits
3.344612061e+02 3.846313087e-01
/raid/gwyn/F22/F22Imapix #6 > gethead CRVAL1 CRVAL2 F22I.fits
3.344612061e+02 3.846313087e-01
this means that:
CENTER 3.344612061e+02,3.846313087e-01 # Coordinates of the image center
works....
/u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 1 -Y 17500 F22Iholem2.fits >! F22Iholem2se.fits
/u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 1 -Y 17500 F22Iholem2.fits >! F22Iholem2sw.fits
/u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 17500 -Y 16686 F22Iholem2.fits >! F22Iholem2nw.fits
/u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 17500 -Y 16686 F22Iholem2.fits >! F22Iholem2ne.fits
/u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 1 -Y 17500 F22Iholem2.weight.fits >! F22Iholem2se.weight.fits
/u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 1 -Y 17500 F22Iholem2.weight.fits >! F22Iholem2sw.weight.fits
/u3/hjmcc/bin/toFITS2d -x 1 -X 20000 -y 17500 -Y 16686 F22Iholem2.weight.fits >! F22Iholem2nw.weight.fits
/u3/hjmcc/bin/toFITS2d -x 20000 -X 17573 -y 17500 -Y 16686 F22Iholem2.weight.fits >! F22Iholem2ne.weight.fits
sex -c bigimage.sex F22Iholem2.fits -CATALOG_NAME F22Iholem2.cat -WEIGHT_IMAGE F22Iholem2.weight.fits
sex -c bigimage.sex F22Iholem2se.fits -CATALOG_NAME F22Iholem2se.cat -WEIGHT_IMAGE F22Iholem2se.weight.fits
sex -c bigimage.sex F22Iholem2ne.fits -CATALOG_NAME F22Iholem2ne.cat -WEIGHT_IMAGE F22Iholem2ne.weight.fits
sex -c bigimage.sex F22Iholem2nw.fits -CATALOG_NAME F22Iholem2nw.cat -WEIGHT_IMAGE F22Iholem2nw.weight.fits
sex -c bigimage.sex F22Iholem2sw.fits -CATALOG_NAME F22Iholem2sw.cat -WEIGHT_IMAGE F22Iholem2sw.weight.fits
sex -c bigimage.sex F22Iholem2ce.fits -CATALOG_NAME F22Iholem2ce.cat -WEIGHT_IMAGE F22Iholem2ce.weight.fits
exsex F22Iholem2ne.cat alpha_j2000 delta_j2000 mag_auto magerr_auto | awk '$3<22.5' > F22Iholem2ne.rdm
exsex F22Iholem2nw.cat alpha_j2000 delta_j2000 mag_auto magerr_auto | awk '$3<22.5' > F22Iholem2nw.rdm
exsex F22Iholem2se.cat alpha_j2000 delta_j2000 mag_auto magerr_auto | awk '$3<22.5' > F22Iholem2se.rdm
exsex F22Iholem2sw.cat alpha_j2000 delta_j2000 mag_auto magerr_auto | awk '$3<22.5' > F22Iholem2sw.rdm
area near the hole: 334.506 0.156675
334.794 0.665195
22:18:01.44
22.19:10.56
00:09:24.03
00:39:54.702
This includes all the borders
image section --> [19569:24603,13110:21998]
or
/u3/hjmcc/bin/toFITS2d -x 19569 -X 5034 -y 13110 -Y F22Iholem2.fits >! F22Iholem2ce.fits
** 29 October 2002
sex -c bigimage.sex F22Iholem2ce.fits -CATALOG_NAME F22Iholem2ce.cat -WEIGHT_IMAGE F22Iholem2ce.weight.fits
move to Morgan:/home/gwyn/red/F22/F22I/hole
exsex F22Iholem2ce.cat X_IMAGE Y_IMAGE ALPHA_j2000 DELTA_j2000 MAG_AUTO MAGERR_AUTO 'flux_radius(2)' flags | grep -v 99.0000 > F22Iholem2ce.xyrdm
** 30 October 2002
Made boxes for F22Iholem2ce.fits
Make a big box:
cat > F22Iceout.reg
# Region file format: DS9 version 3.0
global color=green font="helvetica 10 normal" select=1 edit=1 move=1 delete=1 include=1 fixed=0 source
fk5;box(334.65,0.410935,0.288,0.50852,0) # color=red
(Actually, this is off by one pixel)
Everything inside this box should be deleted.
The objects to be deleted are:
Morgan:/home/gwyn/red/F22/F22I/hole/vdb.rdxym
Morgan:/home/gwyn/mask/seehole.f
Calculates seeing for the hole alone.
The fairly ad hoc perl program
Morgan:/home/gwyn/mask/patchzone.pl
takes F22I.zones.old (the old seeing list)
and replaces the hole parts with the new values
from the new F22I.zones.
patchzone.pl >F22Ihole.zones
Hmm. Apparently /home/gwyn/mask/flag.f wasn't documented.
lasa28:/home/gwyn/mask/flag.f merges all the catalogs and makes nice flags
started on Morgan,
edited on lasa28: this version is the definitive with-hole catalog generator
edited again on Morgan: becomes:
flaghole.f: This version will be the definitive hole-filling versin
** 31 October 2002
Morgan:/home/gwyn/mask/dedub.f now also outputs x1,y1,x2,y2 of the masks.
It remains to convert them to RA,Dec
awk '{print $1,$2}' F22Ice.boxy >fred
xy2skypv F22Iholem2ce.head fred | awk '{print $1,$2}' >fred1
awk '{print $3,$4}' F22Ice.boxy >fred
xy2skypv F22Iholem2ce.head fred | awk '{print $1,$2}' >fred2
paste -d ' ' fred1 fred2 >F22Ice.box
awk '$1<334.506 || $2<0.156675 || $3>334.794 || $4>0.665195' F22I.box >F22Ib.box
mv F22I.box F22Iold.box
cat F22Ib.box F22Ice.box >F22I.box
** 1 November 2002
exsex F22Ice.cat ALPHA_J2000 DELTA_J2000 MAG_AUTO MAGERR_AUTO 'MAG_APER(1)' 'MAGERR_APER(2)' FLAGS > F22Ice.rdmmf
exsex F22Iholem2.cat ALPHA_J2000 DELTA_J2000 MAG_AUTO MAGERR_AUTO 'MAG_APER(1)' 'MAGERR_APER(2)' FLAGS > F22I.rdmmf
** 4 November 2002
working in virmos:/vcom/gwyn/F22Iastro
ln -s ../F22I/F22I.cat
ln -s ../F22Ipatch/F22Ice.cat
awk '$4<334.506 || $4>334.794 || $5<0.156675 || $5>0.665195' F22I.cat > F22Iex.cat
awk '($4<334.506 || $4>334.794 || $5<0.156675 || $5>0.665195) && $56<22.5 ' F22I.cat >F22Iexm.cat
awk '$56<22.5' F22Ice.cat > F22Icem.cat
grep -v # F22Icem.cat >fred
cat F22Iexm.cat fred >F22Im.cat
grep -v # F22Ice.cat > fred2
cat F22Iex.cat fred2 > F22Ipa.cat
awk '$56<22.5 && $48>$56-1 && $64==0' F22Ipa.cat > F22Ia.cat
(that is to say,
MAG_APER(1) <22.5
MAG_AUTO > MAG_APER(1) -1
FLAGS=0
)
alternatively, one could say that it has reasonably compact
fairly bright uncrowded, unsaturated sources.
awk '$0!~/#/ {print $4,$5}' /vcom/gwyn/F22Iastro/F22Ia.cat > F22Ia.rd
awk '$0!~/#/ {print $4,$5,$56}' /vcom/gwyn/F22Iastro/F22Ia.cat > F22Ia.rdm