-----------------------------------------------------------
	Cálculo de polinomios tras cambio de Detector
-----------------------------------------------------------
No es necesario construir el mosaico ya que viene construido
por definición.
==================================
Detector --> Máscara (Polinomio A)
==================================
Se usan todos los pinholes salvo el inicio de cada fila
[debajo del viñeteo de las ruedas], es decir los pinholes: 1,
23, 45, 67, 89.


***********************************************
stdev(RES) X,Y (mm): 
     0.021220425     0.023752805
***********************************************


==================================
Máscara --> Detector (Polinomio B)
==================================

De igual modo se hace el cambio inverso. Asimismo se encuentra
el punto cero que está calculado incluyendo el prescan.

***********************************************
stdev(RES) X,Y (pixeles): 
      0.20396642      0.22564609
***********************************************

***********************************************
Punto cero del mosaico en pixeles (binning 1x1): 
       2100.8192       2010.3491
***********************************************

Para el cálculo de polinomios de cielo a detector necesitamos hacer la
astrometría precisa de un campo poblado de estrellas.

Se toman las imágenes de diciembre de 2022 [NGC1027 campo poblado]

Los resultados de la astrometría fueron

Para el campo 1:
# Coords File: /tmp/._output.ccm  Image: astrofield1.fits
#     Database: astro1_m.db  Solution: astro1_m.sol
# Coordinate mapping status
#     Ra/Dec or Long/Lat fit rms: 0.0575  0.0491   (arcsec  arcsec)
# Coordinate mapping parameters
#     Sky projection geometry: tnx
#     Reference point: 40:41:29.98  61:36:05.82  (degrees  degrees)
#     Reference point: 2039.465  1726.282  (pixels  pixels)
#     X and Y scale: 0.128  0.125  (arcsec/pixel  arcsec/pixel)
#     X and Y axis rotation: 359.046  359.183  (degrees  degrees)
# Wcs mapping status
#     Ra/Dec or Long/Lat wcs rms: 0.0575  0.0491   (arcsec  arcsec)

Con un total de 119 estrellas.



Para el campo 2:
# Coords File: /tmp/._output.ccm  Image: astrofield2.fits
#     Database: astro2_m.db  Solution: astro2_m.sol
# Coordinate mapping status
#     Ra/Dec or Long/Lat fit rms: 0.0813  0.0641   (arcsec  arcsec)
# Coordinate mapping parameters
#     Sky projection geometry: tnx
#     Reference point: 40:41:44.80  61:36:53.42  (degrees  degrees)
#     Reference point: 2251.334  1627.235  (pixels  pixels)
#     X and Y scale: 0.126  0.125  (arcsec/pixel  arcsec/pixel)
#     X and Y axis rotation: 358.840  358.726  (degrees  degrees)
# Wcs mapping status
#     Ra/Dec or Long/Lat wcs rms: 0.0813  0.0641   (arcsec  arcsec)

Con un total de 175 estrellas



Para el campo 3:
# Coords File: /tmp/._output.ccm  Image: astrofield3.fits
#     Database: astro3_m.db  Solution: astro3_m.sol
# Coordinate mapping status
#     Ra/Dec or Long/Lat fit rms: 0.0641  0.0403   (arcsec  arcsec)
# Coordinate mapping parameters
#     Sky projection geometry: tnx
#     Reference point: 40:42:43.64  61:35:34.93  (degrees  degrees)
#     Reference point: 2166.456  1953.428  (pixels  pixels)
#     X and Y scale: 0.127  0.127  (arcsec/pixel  arcsec/pixel)
#     X and Y axis rotation: 358.609  358.805  (degrees  degrees)
# Wcs mapping status
#     Ra/Dec or Long/Lat wcs rms: 0.0641  0.0403   (arcsec  arcsec)


Aplicando la congregación de los resultados (RES) de las distintas
astrometría y reconvirtiéndolas a un único plano astrométrico

Se encuentra un RMS en X de 0.22 y de 0.09 en Y

Con esta combinación entre estrellas más la conversión ccd2sky se puede
calcular el polinomio C


==================================
Cielo --->  Máscara  (Polinomio C)
==================================

Que presenta un RMS de

***********************************************
stdev(RES) X,Y (mm): 
     0.024353957     0.023653422
***********************************************

==================================
Máscara --->  Cielo  (Polinomio D)
==================================

Que presenta un RMS de


Aplicando la rutina mask2sky_pinholes_mosaic.pro se encontrará el polinomio
de conversión de Máscara a Cielo.

***********************************************
stdev(RES) X,Y (arcsec): 
     0.027503892     0.028279144
***********************************************

