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!