!
! astrometry
!
! jwf jul 98
!
! computes astrometric transformation from standard stars
! applies it to objects
!



defi/par p1 ? ima standardstar_coord_table=?

if p1(1:4) .eq. "help" then
        write/out Uses the astrometric data for the standards in table ref_in to
        write/out calculate astrometrical transformation using the x/y
        write/out positions for these same standards in table meas_in and applies it to
        write/out objects and ref.stars (also given in table meas_in).
	write/out
        write/out use as @@ astro  ref_in meas_in astro_out
        write/out "    where ref_in      = table with astrometric data of standard stars "
	write/out "                               (output from @@ std)" 
        write/out "          meas_in = table with x/y/positions measured on the frame to be calibrated. "
	write/out "                column :type contains 2 for standard stars"
	write/out "                                      1 for the ref.stars "
	write/out "                                  and 0 for objects"
	write/out "          astro_out  = table with astrometric positions of ref.stars and objects"
	goto exit
endif

defi/par p2 ? ima standardstar_xy_table=?
defi/par p3 ? ima output_table=?


! select standards
sel/tab {p2} :type.eq.2
!cop/tab {p2} z
astro/transfo {p1} {p2} mean 2000,2000 n,n          !3 111000000,111000000 std
write/out "pausing... write down plate center (scroll upwards)... enter cont when done"
pause



sel/tab {p2} all
astro/com {p2} 2 {p3}
sel/tab {p3} all
cop/tt {p2} :type {p3} :type
comp/tab {p3} :ad = :R_A
comp/tab {p3} :dd = :DEC
cop/tt {p2} :xcen {p3} :xcen
cop/tt {p2} :ycen {p3} :ycen
cop/tt {p2} :ident {p3} :ident
cop/tt {p2} :type {p3} :type
!sho/tab {p3}
write/out The outputtable >> {p3} << may be used as input for compmask
write/out That's it!
 
goto exit



pause

write/out copying columns :xcen,ycen,ident,mag and type...
cop/tt {p2} :xcen z :xcen
cop/tt {p2} :ycen z :ycen
cop/tt {p2} :ident z :ident
cop/tt {p2} :type z :type
cop/tt {p2} :mag z :mag 
crea/col z :ad r*8 f15.6
crea/col z :dd r*8 f15.6

comp/tab z :ad = :R_A
comp/tab z :dd = :DEC

sel/tab z :type.ne.2
cop/tab z {p3}


write/out Info: the outputtable >> {p3} << may be used as input for compmask
write/out that's it!

exit: