I have pulled the EFIGI g images from a website and here I discuss how I cross-match these (PGC) galaxies to the RC3 and create a table file that can be used to view these 4458 galaxy images with point_selector.
I do this work in: $scohtm/scocodes/finding_images/three_sets/efigi/work1
I start with the file listing the galaxy images:
% fitsfind /media/sco/DataDisk1/sco/AD/EFIGI > list.efigi
**** This make a list of the fullpath image names.
I make a list of the names using a custom script named cat PG_LIST.bash:
% cat PG_LIST.bash
#!/bin/bash
printf "Making a list of EFIGI galaxy names \n"
cp list.efigi A
sed 's/\/media\/sco\/DataDisk1\/sco\/AD\/EFIGI\/efigi-1.6\/ima_g\///g' A > B
sed 's/_g.fits//g' B > C
mv C PGC.Names
\rm -f A B
This creates the list of PGC names in the file = PGC.Names
I cross-match these names to get the RC3 index values
% rc3_name_match.sh PGC.Names N
*** This takes just a few seconds to run on sco2019 and the 4458 indices are in the file: rc3_name_match.indexrc3
I assemble the file of RC3 parameeters based on the RC3 index values from above
The quanitites I want are in the file = PARAMS
% cat PARAMS
B(T)_RC3
(B-V)t
(U-B)t
logD25
logR25
Type_RC3
% rc3_table_make.sh rc3_name_match.indexrc3 PARAMS N
*** This also takes just a few seconds to run on sco2019 and creates the table file: RC3TAB
% ls
list.efigi PGC.Names rc3_name_match.indexrc3 RC3TAB.parlab RC3.variables tabmake.out ttt
PARAMS PG_LIST.bash* RC3TAB.images RC3TAB.table S/ tabmake.params
**** At thsi stage I can make a plot and view images with point_selector:
% point_selector RC3TAB Type_RC3 B-Vt N
=============================
I make a table file of images based on general FITS header information:
% survey_headers.sh list.efigi EX1 N
*** This takes about 36 seconds for the 4458 EFIGI galaxies on sco2019
Here are the results:
% cat EX1.parlab
im image table number
uthrs UT time (floating point hours)
fil filter name
source image source
isoobs observation time in ISO8601 format
object object name
imbase FITS image basename
% head -5 EX1.table
# data
1 -9.000000 g efigi 2006-04-22T08:14:15.00 PGC0035314 PGC0035314_g
2 -9.000000 g efigi 2006-04-22T14:43:05.00 PGC0046372 PGC0046372_g
3 -9.000000 g efigi 2006-04-22T12:14:00.00 PGC0040964 PGC0040964_g
4 -9.000000 g efigi 2006-04-22T14:55:54.00 PGC0047060 PGC0047060_g
**** I could use point_selector on EX1, but this is less intersting.
Notice that the RC3 names are different in the rc3_table_make.sh code and
the names in my parlab file (recognized by point_selector). Unfortunately
I can use characters like "(" and ")" as arguments on the bash script command
lines. So, I chnage the names on the fly in rc3_table_make.sh.
Notice in the last table above that the values for "uthrs" are all -9.0 (i.e. empty). This
is because the efigi images do noot have the UT card in their headers. I have a special
subroutine in survey_headers.f that uses the DATE card to build a time in the ISO8601 format.
This last point above is what I'm struggling with. How many special subroutines
should I place in survey_headers.f? This can get away from me quickly!