#Instrucciones para la construccion de un mosaico con imagenes de
#OSIRIS (va sin tildes)
#
#Dadas las extensiones de una imagen de OSIRIS, CCD1, CCD2, binning 2x2,
#proceder de la siguiente forma (en tareas de IRAF):
#
rotate input=CCD1.fits output=CCD1_rot.fits rotation=0.002858 xin=1049.
yin=996. xout=1049. yout=996. ncols=1049. nlines=2054. interpo="linear"
boundar="constant" constan=0. nxblock=30 nyblock=30 verbose=yes
#
rotate input=CCD2.fits output=CCD2_rot.fits rotation=-0.034044 xin=1.
yin=994. xout=13.2508 yout=995.6719 ncols=1062. nlines=2054. interpo="linear"
boundar="constant" constan=0. nxblock=30 nyblock=30 verbose=yes
#
imcreate image=CCD_mosaic.fits naxis=2 naxis1=2112 naxis2=2054
#
imcopy CCD1_rot.fits CCD_mosaic.fits[1:1049,1:2054]
#
imcopy CCD2_rot.fits CCD_mosaic.fits[1050:2112,1:2054]
#
Incluyase en CCD_mosaic.fits la extension 0 de la imagen original.
#
#Estas instrucciones estan codificadas en las herramientas de OSIRIS
#elaboradas por Alessandro Ederoclite y pueden ser obtenidas en
#http://www.iac.es/galeria/ale/OSIRISPipeline.html (OSIToS)
#
#A titulo de ejemplo, para componer el mosaico de la pre-imagen
#usada en el Caso_1 del Cookbook del MD, procedimos de la siguiente forma
#(todos los ficheros de entrada y salida indicados a continuacion pueden
#ser encontrados en el directorio que contiene este readme):
#
#Partimos de las extensiones 1 y 2 de la pre-imagen: 
#b0000039008-20100425-OSIRIS-OsirisBroadBandImage_1.fits
#b0000039008-20100425-OSIRIS-OsirisBroadBandImage_2.fits
#
#Componemos el mosaico con la herramienta make_mosaic de A. Ederoclite,
#a partir de las extensiones 1 y 2. El resultado es la imagen 
#V8412Oph_2x2_IPA240_preimaging_mosaic.fits
#
#Para incluir la astrometria en la cabecera del mosaico procedase de la
#siguiente forma:
#
#Con noao.digiphot.daophot.daofind encontramos el fichero
#V8412Oph_2x2_IPA240_preimaging_mosaic.xyxy
#(coordenadas fisicas de fuentes resueltas a 3-sigma)
#
#Hacemos un apareamiento de las fuentes con un catalogo en coordenadas
#ecuatoriales suficientemente preciso (e.g. USNO-B1.0) mediante la
#siguiente linea de comando de IRAF:
#
ccxymatch V8412Oph_2x2_IPA240_preimaging_mosaic.xyxy catalog_astrometry_only.tsv V8412Oph_2x2_IPA240_preimaging_mosaic.raxy xin=674.98 yin=944.23 xmag=0.254 ymag=0.254 xrotation=90. yrotation=270. projection="tan" lngref=254.871648 latref=-12.947478 lngcolumn=1 latcolumn=2 xcolumn=1 ycolumn=2 lngunits="degrees" latunits="degrees" nmatch=40
#
#donde catalog_astrometry_only.tsv es el fichero que contiene el catalogo de
#la region de interes. De esto resulta el fichero
#V8412Oph_2x2_IPA240_preimaging_mosaic.raxy
#
#A continuacion, copiamos el mosaico haciendo
#
imcopy V8412Oph_2x2_IPA240_preimaging_mosaic.fits V8412Oph_2x2_IPA240_preimaging_mosaic_order3.fits
#
#y procedemos a buscar e incorporar las soluciones astrometricas (orden 3) a
#la imagen resultante mediante el siguiente comando:
#
ccmap V8412Oph_2x2_IPA240_preimaging_mosaic.raxy V8412Oph_2x2_IPA240_preimaging_mosaic_order3.db solutio="V8412Oph_2x2_IPA240_preimaging_mosaic_order3.sol" images="V8412Oph_2x2_IPA240_preimaging_mosaic_order3.fits" results=V8412Oph_2x2_IPA240_preimaging_mosaic_order3.res xcol=3 ycol=4 lngcol=1 latcol=2 lngunits="degrees" latunits="degrees" insystem="j2000" refsystem="j2000" projection="tnx" fitgeometry="general" function="polynomial" xxorder=3 yyorder=3 yxorder=3 xyorder=3 xxterms="full" yxterms="full" maxiter=100 reject=3.0 update=yes pixsyst="physical"
#
#Con esto, la cabecera de V8412Oph_2x2_IPA240_preimaging_mosaic_order3.fits
#contendra las componentes lineal y de orden superior de la solucion astrometrica 
#(RMS~0.2 arscsec). Deformamos la imagen atendiendo a las componentes de orden 
#superior, usando mscred.mscimage:
#
mscimage V8412Oph_2x2_IPA240_preimaging_mosaic_order3.fits V8412Oph_2x2_IPA240_preimaging_mosaic_order3_warped.fits
#
#Notese que en las imagenes del mosaico
#DEC     = '-12:57:47.228'      / d:m:s                                          
#RA      = '+16:59:29.905'      / h:m:s      
#
#Estas coordenadas corresponden al apuntado del telescopio y estan proyectadas
#alrededor del centro del CCD1 de OSIRIS. Por lo pronto, el "apuntado" del 
#telescopio para el MD corresponde a la proyeccin del origen (0,0) de
#la mascara en el mosaico, es decir, el pixel fisico (1061.5595,995.8728) a
#binning 2x2. Por tanto, se editaron RA y DEC en todas las imagenes anteriores
#de forma que el JSky en el MD proporcione coordenadas a los targets consistentes 
#con esta definicion del centro del mosaico (verifiquese esto haciendo
#
hsel V8412Oph_2x2_IPA240_preimaging_mosaic*fits $I,RA,DEC
#
#y comparense con las coordenadas indicadas arriba).
#

