Obs. Auto : algo reconnaissance du champ

[Marc SERRAU]

Salut,
Serait-il possible de remplacer l'algorithme de reconnaissance du ciel actuellement dans la fonction ObsAuto par celui de la fonction de reconnaissance sur tout le ciel qui semble bien plus performant ?

Merci

[Patrice Le Guen]

Bonjour,

Je suis d'accord avec cette requête. J'ai bien noté à l'usage cette différence de performance entre la fonction "Analyse/Astrométrie/Etalonnage automatique" et la fonction de rencentrage automatique dans "Télescope-Observatoire/Observations automatique" Je suis tout de même surpris que ces deux fonctions dont l'objectif est le même aient un comportement différent.

Par contre la performance de la fonction "Analyse/Astrométrie/Etalonnage automatique" de PRISM me semble perfectible. Je souhaiterais qu'elle soit améliorée par l'implémentation d'un Plugin de type "PINPOINT LE".

En effet, depuis de nombreuse années j'utilise MAXIMDL. J'ai acheté PRISM8 la semaine dernière car le module "Observation automatique" de MAXIMDL me semblait dépassé comparé à celui de PRISM. Celui de MAXIMDL ne permet que de lancer une séquence d'acquisition de poses en passant les différents filtres, les bias, dark et Flats. Ce qui m'intéressait dans PRISM était de pouvoir demander au télescope de pointer un objet automatiquement. La faible performance du module de recentrage automatique du module "Observation automatique" de PRISM place ce dernier au même niveau que celui de MAXIMDL et ne représente plus autant d'intérêt pour moi.

Voici la preuve de la supériorité de MAXIMDL sur PRISM pour le calcul d'étalonnage astrométrique automatique.
Avec MAXIMDL j'ai fait une pose en binning 2 de 20s sur un champs proche de l'étoile SAIPH dans ORION sans que l'étoile soit dans le champs de la CCD (1.5° environ). Le log de MAXIMDL donne ceci :
Analyze, PinPoint Astrometry: test astrometrie MAXIMDL.fits
Trying center at RA 05 47 48.0, Dec -09 39.0 in test astrometrie MAXIMDL.fits
Trying center at RA 05 55 40.7, Dec -09 39.0
Trying center at RA 05 47 48.0, Dec -08 11 14.8
Trying center at RA 05 39 55.3, Dec -09 39.0
Matched 342 of 402 image and 1841 catalog stars;
Average residual 0.6 arcsec; order 4
RA 05h 39m 26s, Dec -09° 39' 46"
Pos Angle +355° 51', FL 530.7 mm, 4.20"/Pixel
J'ai volontairement demandé une recherche difficile et forcé MAXIMDL à faire une recherche en spirale autour de la position théorique de SAIPH. Il trouve cependant le résultat au bout d'environ 2 minutes de calculs.
Avec PRISM j'ai chargé l'image faite par MAXIMDL et j'ai lancé la fonction "Analyse/Astrométrie/Etalonnage automatique". Malgré les changements de paramètres des plus serrés aux plus larges PRISM n'a jamais réussi à étalonner la position réelle de cette image. En pièce jointe j'ai placé le LOG de PRISM. Le format FITS est interdit en PJ, je peux transmettre l'image par mail ou via un autre moyen. J'ai mis l'entête de l'image FITS dans un fichier en PJ.

Ma requête a t-elle l'espoir d'aboutir ? Sinon comment faire pour s'assurer d'un étalonnage astrométrique automatique réussi via le module d'observations automatiques ?

PS : il semble que je ne soit pas autorisé à transférer des pièces jointes sur ce forum. Je recopie donc ci-dessous les LOG dont je parle dans mon texte :

Entête de l'image lue par PRISM et acquise par MAXIMDL
01 : SIMPLE = T
02 : BITPIX = 16 /8 unsigned int, 16 & 32 int, -32 & -64 real
03 : NAXIS = 2 /number of axes
04 : NAXIS1 = 1663 /fastest changing axis
05 : NAXIS2 = 1252 /next to fastest changing axis
06 : BSCALE = 1.0000000000000000 /physical = BZERO + BSCALE*array_value
07 : BZERO = 32768.000000000000 /physical = BZERO + BSCALE*array_value
08 : DATE-OBS= '2014-03-13T20:30:40' /YYYY-MM-DDThh:mm:ss observation start, UT
09 : EXPTIME = 20.000000000000000 /Exposure time in seconds
10 : EXPOSURE= 20.000000000000000 /Exposure time in seconds
11 : CCD-TEMP= -25.000000000000000 /CCD temperature at start of exposure in C
12 : XPIXSZ = 10.800000000000001 /Pixel Width in microns (after binning)
13 : YPIXSZ = 10.800000000000001 /Pixel Height in microns (after binning)
14 : XBINNING= 2 /Binning factor in width
15 : YBINNING= 2 /Binning factor in height
16 : XORGSUBF= 0 /Subframe X position in binned pixels
17 : YORGSUBF= 0 /Subframe Y position in binned pixels
18 : READOUTM= 'Image Quality' / Readout mode of image
19 : ISOSPEED= 'High Gain'
20 : FILTER = 'Luminance' / Filter used when taking image
21 : IMAGETYP= 'Light Frame' / Type of image
22 : OBJCTRA = '05 47 48' / Nominal Right Ascension of center of image
23 : OBJCTDEC= '-09 38 57' / Nominal Declination of center of image
24 : OBJCTALT= ' 26.9188' / Nominal altitude of center of image
25 : OBJCTAZ = '222.8091' / Nominal azimuth of center of image
26 : OBJCTHA = ' 2.5283' / Nominal hour angle of center of image
27 : PIERSIDE= 'EAST ' / Side of pier telescope is on
28 : SITELAT = '43 09 17' / Latitude of the imaging location
29 : SITELONG= '05 53 41' / Longitude of the imaging location
30 : JD = 2456730.3546296298 /Julian Date at start of exposure
31 : JD-HELIO= 2456730.3550325693 /Heliocentric Julian Date at exposure midpoint
32 : AIRMASS = 2.1939324766282007 /Relative optical path length through atmosphere
33 : FOCALLEN= 530.00000000000000 /Focal length of telescope in mm
34 : APTDIA = 106.00000000000000 /Aperture diameter of telescope in mm
35 : APTAREA = 8824.7340095043182 /Aperture area of telescope in mm^2
36 : EGAIN = 0.50000000000000000 /Electronic gain in e-/ADU
37 : SWCREATE= 'MaxIm DL Version 5.15' /Name of software that created the image
38 : SBSTDVER= 'SBFITSEXT Version 1.0' /Version of SBFITSEXT standard in effect
39 : OBJECT = ' '
40 : TELESCOP= ' ' / telescope used to acquire this image
41 : INSTRUME= 'QSI 583ws S/N 00504218 HW 06.00.00 FW 05.02.06 PI 6.4.962.4'
42 : OBSERVER= ' '
43 : NOTES = ' '
44 : FLIPSTAT= ' '
45 : CSTRETCH= 'Medium ' / Initial display stretch mode
46 : CBLACK = 11260 /Initial display black level in ADUs
47 : CWHITE = 12416 /Initial display white level in ADUs
48 : PEDESTAL= 0 /Correction to add for zero-based ADU
49 : SWOWNER = 'pat ' / Licensed owner of software
50 : END

Log de l'étalonnage de PRISM tronqué pour ne pas dépasser les 60000 caractères :
Débute le processus dans le "thread" : $00000094
Extrait du catalogue les étoiles référence...
Extrait du catalogue : 606 étoiles de référence
Extrait les étoiles de l'image..
Extraction rapide des étoiles terminée, 726 trouvées, débute l'extraction avancée...
Nombre d'étoiles extraites de l'image : 543 -> 0.392 sec
Extraction des étoiles de l'image CCD
# X Y Mg Alpha = Delta = -> Carte du ciel
1 836.7 609.4 2.07 05h47m45.379s -09°40'10.81''
2 1567 311.1 7.09 05h44m17.562s -10°00'58.97''
.....
199 747.2 848.1 13.01 05h48m25.020s -09°23'27.62''
200 544.4 698.2 13.01 05h49m08.378s -09°33'56.93''

# Xccd Yccd Flux -> Image CCD
1 200 610.8 3.93E5
2 1262 1227 3.354E5
...
199 1614 667.2 7881
200 797 1177 7877
Nbre d'étoiles (les plus brillantes) limité à 200 pour la comparaison champ-catalogue.
Compare les étoiles de l'image et de la carte du ciel...
XCCD YCCD Xcat YCat CCD_Index Cat_index nbVotes RA_Cat DEC_Cat Flux_Image Magnitude
982.5 695 617.7 767.6 26 185 510 05h48m48s -09°29'06'' 78996.93 12.88
865 552 848.8 704.8 164 142 505 05h47m42s -09°33'30'' 9631.949 12.43
...
1386 933.3 115.6 910 192 177 153 05h51m10s -09°19'04'' 8237.954 12.82
1274 124.3 639.6 95 62 6 148 05h48m41s -10°16'12'' 33546.95 7.7
Nombre d'étoiles pour le calcul du premier polynome de passage : 200
TransFormation "Trigmatch" de passage pour validation des centroides, degre : 1, étoiles utilisées : 5, basé sur résidu max en pixel : 0.5
Verification avec centroides de : 9 pixels de demi-largeur
Nombre d'étoiles à tester au centroide : 200
XCentro Ycentro Centro OK ?
1062 880 Echec du centroide
826 719 Echec du centroide
...
79 209 OK xc=78.9 yc=209.0 (19x19) -> psf= 1.6x1.6 pixel
...
1581 1234 coordonnées identiques
1043 869 Echec du centroide
13 étoiles après test, 200 avant test de centroide...
Nombre d'étoiles validées par centroide : 13
Degré : 3, nombre d'étoiles prises en compte : 12
OK [5] X=78.928 Y=209.001 res.=4381.261" dist=83%
OK [7] X=165.677 Y=267.689 res.=4484.744" dist=73%
...
NO [8] X=165.676 Y=267.689 res.=1322358.109" dist=73%
Référence -> distance max = 92%
Résidu : 477236.97"

[Cyril Cavadore]

A mon avis il faudrait que j'ai une image pour voir.
Merci

[Patrice Le Guen]

Bonjour Cyril,

Voici l'image au format FITS dont je parle dans mon texte
http://www.astrosurf.com/leguenp/photos/tests_logiciel/

Je peux prendre beaucoup d'autres photos que MAXIMDL arrive à étalonner en une fraction de secondes et que PRISM n'arrive pas à étalonner du tout.

Je suis disponible pour t'apporter toute mon expérience sur l'utilisation de MAXIMDL afin d'améliorer PRISM.

Cordialement.

[Cyril Cavadore]

Il faudra apprendre a quitter MAXIMDL, car PRISM est largement plus puissant

Ton image fait 3.6M et semble incomplete.... impossible de la charger.

Cdt.

[Cyril Cavadore]

Serait-il possible de remplacer l'algorithme de reconnaissance du ciel actuellement dans la fonction ObsAuto par celui de la fonction de reconnaissance sur tout le ciel qui semble bien plus performant ?

C'est exactement le meme.

Sf que les gens tendent a augmenter ce nombre , et l'algorithme perd la performance en puissance 4 de ce chiffre... Donc plus il est petit, plus c'est rapide.

[Patrice Le Guen]

Il faudra apprendre a quitter MAXIMDL, car PRISM est largement plus puissant

Ton image fait 3.6M et semble incomplete.... impossible de la charger.

Cdt.

Effectivement le téléchargement de mon image est incomplet. Je viens de remettre l'image sur mon site (http://www.astrosurf.com/leguenp/photos ... XIMDL.fits). J'ai fais un test de téléchargement de l'image. Pourrais-tu réessayer de faire un test d'étalonnage avec cette image ?
L'orientation de la caméra est NORD/SUD (en fait à -4°) le champs X= 1°56 et le champs Y=1°27

MAXIMDL fait l'étalonnage astrométrique en moins de 3 minutes et trouve RA 05h 39m 26s, Dec -09° 39' 46" (cf. mon message du 13/03/2014 23:40) alors que le télescope indiquait RA 05 47 48.0, Dec -09 39.0 soit un écart de distance d'environ 32 minutes d'arc entre la position réelle sur le ciel et les codeurs de la monture.

PRISM n'arrive pas à étalonner l'image il affiche "Une erreur est intervenue pendant la calibration astrométrique..". Même en remettant à 60 les paramètres du nombre d'étoiles du catalogue de référence et de l'image PRISM n'arrive pas à étalonner l'image. Si j'indique dans PRISM les coordonnées Alpha et Delta trouvée par MAXIMDL alors effectivement PRISM réalise l'étalonnage rapidement et retrouve exactement les mêmes valeurs que MAXIMDL à moins d'une seconde d'arc. J'ai diminué la valeur en Alpha, de minute en minute de 05h47 à 05h44, pour aider PRISM à se rapprocher de la valeur à trouver. Il a réussi à étalonner pour Alpha = 05h44m47s. Soit à 5 minutes d'heures ou 20 minutes d'arcs entre la position réelle sur le ciel et les codeurs de la monture..

Dois-je en conclure que PRISM n'est pas capable d'étalonner s'il y a trop d'écarts entre la position réelle et la position des codeurs de la monture ? Pour mon cas cela représenterai 20minutes d'arc sur un champ de 1,5°. C'est dans ce sens que je trouve que l'algorithme manque de performance, j'aurais du écrire "manque de précision". Dans ce cas particulier, MAXIMDL est largement plus puissant que PRISM.

Que faut-il indiquer dans les paramètres "Réglage de la reconnaissance auto" pour que PRISM soit plus puissant que MAXIMDL en étalonnage? Serait-il possible d'implémenter une recherche en spirale autour de la position théorique si PRISM n'arrive pas à étalonner aux coordonnées indiquées ?

Dans l'état actuel je ne peux pas utiliser le module "Observations automatiques" de PRISM si la fonction d'étalonnage astrométrique ne tolère pas plus de 30minutes d'arc d'erreur de position du télescope.

Cordialement.

[Cyril Cavadore]

J'ai pris ton image et si j'utilise la fonction "telescope/Retrouver la position du telescope dans tout le le ciel"

Ca marche en 1.5sec

Curieux que tu aies 30 arcmin d'erreur c'est bcp qd meme sur ton pointage.
C'est genant...
Il y a le modele de pointage pour reduire ces erreurs de pointage.

[Patrice Le Guen]

Bonjour,

La fonction "telescope/Retrouver la position du telescope dans tout le le ciel" trouve effectivement les coordonnées de l'image à 1,5° des coordonnées du télescope. La fonction est puissante et de loin plus rapide que celle de MAXIMDL. Cependant comment faire pour l'utiliser dans le mode "Observations automatiques" ?

Mon installation n'est pas fixe et je sors la monture sur ma terrasse uniquement quand il fait beau et qu'il n'y a pas de pollen dans l'air et rarement plus de 8 jours d'affiler. Je fais une mise en station sur la polaire à 1' d'arc environ. Lorsque je "park" la monture (AH=6h, Dec=90°) et que je "dépark" pour aller vers le sud j'ai facilement une erreur de pointage importante. J'ai bien noté l'existence de la fonction "modèle de pointage" je crois que je ne vais pas tarder à la tester.

Cordialement.

[Patrice Le Guen]

Serait-il possible de remplacer l'algorithme de reconnaissance du ciel actuellement dans la fonction ObsAuto par celui de la fonction de reconnaissance sur tout le ciel qui semble bien plus performant ?

C'est exactement le meme.

Sf que les gens tendent a augmenter ce nombre , et l'algorithme perd la performance en puissance 4 de ce chiffre... Donc plus il est petit, plus c'est rapide.

AstromA.JPG

Après la démonstration qui vient d'être faite sur les différences de performances entre les fonctions "telescope/Retrouver la position du telescope dans tout le le ciel" et la reconnaissance automatique du champs du module "observations automatiques" je note que ma démarche converge tout à fait vers la demande de Marc SERREAU. Je propose donc que l'algorithme de la fonction "telescope/Retrouver la position du telescope dans tout le le ciel"" soit utilisée dans de module "observations automatiques" pour la reconnaissance automatique du champs avec possibilité de renseigner le paramètre "Distance d'exploration maximum". Même si les algorithmes sont les mêmes, les paramètres sont différents (orientation de la caméra supposé au nord), en tout cas le paramètre de "Distance d'exploration maximum" semble fixé à une valeur très faible et en tout cas à moins de 1,5°.